Case study of a moisture intrusion over the Arctic with the ICOsahedral Non-hydrostatic (ICON) model: resolution dependence of its representation

The Arctic is warming faster than the global average and any other region of a similar size. One important factor in this is the poleward atmospheric transport of heat and moisture, which contributes directly to the surface and air warming. In this case study, the atmospheric circulation and spatio-temporal structure of a moisture intrusion event is assessed, which occurred from 5 to 7 June 2017 over the Nordic seas during an intensive measurement campaign over Svalbard. This analysis focuses on high-spatial-resolution simulations with the ICON (ICOsahedral Non-hydrostatic) model which is put in context with coarser-resolution runs as well the ERA5 reanalysis. A variety of observations including passive microwave satellite measurements is used for evaluation. The global operational ICON forecasts from the Deutscher Wetterdienst (DWD) at 13 km horizontal resolution are used to drive high-resolution Limited-Area Mode (LAM) ICON simulations over the Arctic with 6 and 3 km horizontal resolutions. The results show the skilful capacity of the ICON-LAM model to represent the observed spatio-temporal structure of the selected moisture intrusion event and its signature in the temperature, humidity and wind profiles, and surface radiation. In several aspects, the high-resolution simulations offer a higher accuracy than the global simulations and the ERA5 reanalysis when evaluated against observations. One feature where the high-resolution simulations demonstrated an advanced skill is the representation of the changing vertical structure of specific humidity and wind associated with the moisture intrusion passing Ny-Ålesund (western Svalbard); the humidity increase at 1–2 km height topped by a dry layer and the development of a low-level wind jet are best represented by the 3 km simulation. The study also demonstrates that such moisture intrusions can have a strong impact on the radiative and turbulent heat fluxes at the surface. A drastic decrease in downward shortwave radiation by ca. 500 W m−2 as well as an increase in downward longwave radiation by ca. 100 W m−2 within 3 h have been determined. These results highlight the importance of both moisture and clouds associated with this event for the surface energy budget.

Comparison of elevation-dependent warming and its drivers in the tropical and subtropical Andes

Mountain regions have been recognized to be more sensitive to climate and environmental changes, and in particular to global warming. Several studies report on elevation-dependent warming (EDW), i.e., when warming rates are different in different altitude ranges, particularly focusing on the enhancement of warming rates with elevation. The Andean chain proved to be a relevant climate change hot-spot with positive temperature trends and a widespread glacier retreat over the recent decades. To assess and to better understand elevation dependent warming in this mountain region and to identify its possible dependence on latitude, the Andean Cordillera was split into five domains, three pertaining to the tropical zone and two pertaining to the Subtropics. Further, for each area the eastern and western faces of the mountain range were separately analyzed. An ensemble of regional climate model (RCM) simulations participating in the Coordinated Regional Climate Downscaling Experiment (CORDEX), consisting of one RCM nested into eight different global climate models from the CMIP5 ensemble was considered in this study. EDW was assessed by calculating the temperature difference between the end of the century (2071–2100) and the period 1976–2005 and relating it to the elevation. Future projections refer to the RCP 8.5 high-emission scenario. Possible differences in EDW mechanisms were identified using correlation analyses between temperature changes and all the variables identified as possible EDW drivers. For the maximum temperatures, a positive EDW signal (i.e. enhancement of warming rates with elevation) was identified in each side of both the tropical and subtropical Andes and in all seasons. For the minimum temperatures, on the contrary, while a positive EDW was identified in the Subtropics (particularly evident in the western side of the chain), the Tropics are characterized by a negative EDW throughout the year. Therefore, the tropical boundary marks a transition between discordant EDW behaviours in the minimum temperature. In the Tropics and particularly in the inner Tropics, different EDW drivers were identified for the minimum temperature, whose changes are mostly associated with changes in downward longwave radiation, and for the maximum temperature, whose changes are mainly driven by changes in downward shortwave radiation. This might explain the opposite EDW signal found in the tropical Andes during daytime and nighttime. Changes in albedo are an ubiquitous driver for positive EDW in the Subtropics, for both the minimum and the maximum temperature. Changes in longwave radiation and humidity are also EDW drivers in the Subtropics but with different relevance throughout the seasons and during daytime and nighttime. Also, the western and eastern sides of the Cordillera might be influenced by different EDW drivers.

Toward the Estimation of All-Weather Daytime Downward Longwave Radiation over the Tibetan Plateau

Downward longwave radiation (DLR) is a critical parameter for radiation balance, energy budget, and water cycle studies at regional and global scales. Accurate estimation of the all-weather DLR with a high temporal resolution is important for the estimation of the surface net radiation and evapotranspiration. However, most DLR products involve instantaneous DLR estimates based on polar orbiting satellite data under clear-sky conditions. To obtain an in-depth understanding of the performances of different models in the estimation of DLR over the Tibetan Plateau, which is a focus area of climate change study, this study tests eight methods for clear-sky conditions and six methods for cloudy conditions based on ground-measured data. It is found that the Dilley and O’Brien model and the Lhomme model are most suitable for clear-sky conditions and cloudy conditions, respectively. For the Dilley and O’Brien model, the average root mean square error (RMSE) of DLR under clear-sky conditions is approximately 22.5 W/m2 for nine ground sites; for the Lhomme model, the average RMSE is approximately 23.2 W/m2. Based on the estimated cloud fraction and meteorological data provided by the China Land Surface Data Assimilation System (CLDAS), hourly all-weather daytime DLR with a 0.0625° resolution over the Tibetan Plateau is estimated. Results demonstrate that the average RMSE of the estimated hourly all-weather DLR is approximately 26.4 W/m2. With the combined all-weather DLR model, the hourly all-weather daytime DLR dataset with a 0.0625° resolution from 2008 to 2016 over the Tibetan Plateau is generated. This dataset can contribute to studies associated with the radiation balance and energy budget, water cycle, and climate change over the Tibetan Plateau.

An Intraseasonal Mode Linking Wintertime Surface Air Temperature over Arctic and Eurasian Continent

Key processes associated with the leading intraseasonal variability mode of wintertime surface air temperature (SAT) over Eurasia and the Arctic region are investigated in this study. Characterized by a dipole distribution in SAT anomalies centered over north Eurasia and the Arctic, respectively, and coherent temperature anomalies vertically extending from the surface to 300hPa, this leading intraseasonal SAT mode and associated circulation have pronounced influences on global surface temperature anomalies including the East Asian winter monsoon region. By taking advantage of realistic simulations of the intraseasonal SAT mode in a global climate model, it is illustrated that temperature anomalies in the troposphere associated with the leading SAT mode are mainly due to dynamic processes, especially via the horizontal advection of winter mean temperature by intraseasonal circulation. While the cloud-radiative feedback is not critical in sustaining the temperature variability in the troposphere, it is found to play a crucial role in coupling temperature anomalies at the surface and in the free-atmosphere through anomalous surface downward longwave radiation. The variability in clouds associated with the intraseasonal SAT mode is closely linked to moisture anomalies generated by similar advective processes as for temperature anomalies. Model experiments suggest that this leading intraseasonal SAT mode can be sustained by internal atmospheric processes in the troposphere over the mid-to-high latitudes by excluding forcings from Arctic sea ice variability, tropical convective variability, and the stratospheric processes.

Characteristics of Turbulence and Aerosol Optical and Radiative Properties during Haze–Fog Episodes in Shenyang, Northeast China

The characteristics of turbulence in the planetary boundary layer (PBL) and the aerosol optical and radiative properties during haze and haze–fog mixed episodes on 22–27 January 2021, in Shenyang, a provincial city in Northeast China, were analyzed using meteorological and aerosol observations. During the haze episode, the hourly mean PM2.5 concentration reached a maximum of 337 µg m−3 and visibility decreased to 1.6 km. The PM2.5 concentration decreased gradually during the haze–fog mixed episode as a result of the scavenging effects of fog, but visibility mostly remained below 1 km owing to high ambient relative humidity (>90%). During the haze–fog mixed episode, an increasing proportion of PM2.5 led to a higher ratio of the backward to the total scattering coefficient. As fog occurred, downward shortwave radiation arriving at the surface was significantly reduced, and upward longwave radiation increased and almost equaled the downward longwave radiation, which can be used as a good indicator for distinguishing haze and fog. Mechanical turbulence was weak during both episodes, and latent heat flux varied within a wider range during the haze–fog mixed episode. The PBL dynamic structure affected the vertical distribution of aerosols/fog droplets. Aerosol-rich layers appeared at altitudes below 0.5 km and above 0.6 km during the haze episode. The elevated aerosol layer was related to the aerosol transport from upstream polluted areas caused by strong upper-level turbulence, and it began to mix vertically after sunrise because of convective turbulence. Aerosols and fog droplets were mostly trapped in a shallower PBL with a height of 0.2–0.4 km during the haze–fog mixed episode because of weaker turbulence.

Analysis of the Radiation Fluxes over Complex Surfaces on the Tibetan Plateau

Analysis of long-term, ground-based observation data on the Tibetan Plateau help to enhance our understanding of land-atmosphere interactions and their influence on weather and climate in this region. In this paper, the daily, monthly, and annual averages of radiative fluxes, surface albedo, surface temperature, and air temperature were calculated for the period of 2006 to 2019 at six research stations on the Tibetan Plateau. The surface energy balance characteristics of these six stations, which include alpine meadow, alpine desert, and alpine steppe, were then compared. The downward shortwave radiation at stations BJ, QOMS, and NAMORS was found to decrease during the study period, due to increasing cloudiness. Meanwhile, the upward shortwave radiation and surface albedo at all stations were found to have decreased overall. Downward longwave radiation, upward longwave radiation, net radiation, surface temperature, and air temperature showed increasing trends on inter-annual time scales at most stations. Downward shortwave radiation was maximum in spring at BJ, QOMS, NADORS, and NAMORS, due to the influence of the summer monsoon. Upward shortwave radiation peaked in October and November due to the greater snow cover. BJ, QOMS, NADORS, and NAMORS showed strong sensible heat fluxes in the spring while MAWORS showed strong sensible heat fluxes in the summer. The monthly and diurnal variations of surface albedo at each station were “U” shaped. The diurnal variability of downward longwave radiation at each station was small, ranging from 220 to 295 W·m−2.The diurnal variation in surface temperature at each station slightly lagged behind changes in downward shortwave radiation, and the air temperature, in turn, slightly lagged behind the surface temperature.

Influence of water vapor and aerosols on downward longwave radiation in the high mountain region of Musala peak, Bulgaria

This study reveals the effect of aerosols and water vapor on downward longwave radiation in the high mountain region of Musala peak, Bulgaria. The investigated period is 01/01/2017 (Jan. 1st 2017) – 30/09/2019 (Sep. 30th 2019). Statistical methods are the main tool for discovering the relationships between the different elements. The results indicate that air temperature is the leading factor for downward longwave radiation, specific humidity, and amount of aerosols in the air. That is why, in order to reveal the pure relationship between downward longwave radiation, specific humidity and amount of aerosols in the atmosphere, the air temperature was cleared from the data series. After this procedure, the results show that specific humidity has a significant influence on the downward longwave radiation flux, and an increase of 1% of the specific humidity results in an increase of about 12-15% in the values   of the downward longwave radiation. At air temperatures around 0ºC the influence of water vapor on the downward longwave flux is highest, which is due to the phase transitions of the water – a process during which release/absorption of radiation in the longwave spectrum occurs. The amount of aerosols in the atmosphere also has a significant effect on this type of radiation, and an increase of 1% of the amount of aerosols in the air at air temperatures above –5.5°C results in an increase of the downward longwave radiation of about 2-4%. The findings of this study show that coarser and opaque aerosol particles have a stronger effect on downward longwave radiation. In the area of Musala peak, as the air temperature rises, there is an increase in the amount of aerosols in the air, a decrease in their size, and a transition from transparent to opaque aerosols. The combination of these different tendencies causes the influence of aerosols on downward longwave radiation to be strongest in the middle temperature interval – air temperatures between –5.5°C and +5.5°C. Due to the increased total amount of aerosols and increased amount of opaque aerosols, their influence on downward longwave radiation is significant also at air temperatures above 5.5°C.

Warm and moist atmospheric flow caused a record minimum July sea ice extent of the Arctic in 2020

The satellite observations unveiled that the July sea ice extent of the Arctic shrank to the lowest value in 2020 since 1979, with a major ice retreat in the Eurasian shelf seas including Kara, Laptev, and East Siberian Seas. Based on the ERA-5 reanalysis products, we explored the impacts of warm and moist air-mass transport on this extreme event. The results reveal that anomalously high energy and moisture converged into these regions in the spring months (April to June) of 2020, leading to a burst of high moisture content and warming within the atmospheric column. The convergence is accompanied by local enhanced downward longwave radiation and turbulent fluxes, which is favorable for initiating an early melt onset in the areas with severe ice loss. Once the melt begins, solar radiation played a decisive role in leading to further sea ice depletion due to ice-albedo positive feedback. The typical trajectories of the synoptic cyclones that occurred on the Eurasian side in spring 2020 agree well with the path of atmospheric flow. Assessments suggest that variations in characteristics of the spring cyclones are conducive to the severe melt of sea ice. We argue that large-scale atmospheric circulation and synoptic cyclones act in concert to trigger the exceptional poleward transport of total energy and moisture from April to June to cause this new record minimum of sea ice extent in the following July.

Convection-Aerosol Interactions in the United Arab Emirates: A Sensitivity Study

The Weather Research and Forecasting (WRF) model is used to investigate convection-aerosol interactions in the United Arab Emirates for a summertime convective event. Both an idealised and scaled versions of a 7-year climatological aerosol distribution are considered. The convection on 14 August 2013 was triggered by the low-level convergence of the circulation associated with the Arabian Heat Low (AHL) and the daytime sea-breeze circulation. The cold pools associated with the convective events, as well as the low-level wind convergence along the Intertropical Discontinuity (ITD) earlier in the day, explain the dustier environment, with Aerosol Optical Depths (AODs) in excess of two. Due to a colder surface and air temperature, the AHL is incorrectly represented in WRF, which leads to a mismatch between the observed and modelled clouds and precipitation. Employing interior nudging in the outermost grids of the three-nested simulation has a small but positive impact on the model predictions of the innermost nest. This is because the higher temperatures from more accurate boundary conditions are offset by colder temperatures from locally enhanced precipitation, the latter arising from a shift in the position of the AHL. Numerical experiments revealed a high sensitivity to the aerosol properties. In particular, replacing 20 % of the rural aerosols by carbonaceous particles has an impact on the surface radiative fluxes comparable to increasing the aerosol loading by a factor of 10, with a daily-averaged reduction in the UAE-averaged net shortwave radiation flux of ~90 W m−2 and an increase in the net longwave radiation flux of ~51 W m−2. However, in the former, WRF generates 20 % more precipitation than in the latter, due to a broader and weaker AHL. The surface downward and upward shortwave and upward longwave radiation fluxes are found to scale linearly with the aerosol loading, while the downward longwave radiation flux varies by less than ±12 W m−2 when the aerosol amount and/or properties are changed. An increase in the aerosol loading also leads to drier conditions due to a shift in the position of the AHL and rainfall occurring in a drier region, with a domain-wise decrease in the daily accumulated rainfall of 16 % when the aerosol loading is increased by a factor of 10. In addition, the onset of convection is also delayed.