Response of Convective Boundary Layer and Shallow Cumulus to Soil Moisture Heterogeneity: A Large‐Eddy Simulation Study

In this study, the impact of varying soil moisture heterogeneity (spatial variance and structure) on the development of the convective boundary layer and shallow cumulus clouds was investigated. Applying soil moisture heterogeneity generated via spatially correlated Gaussian random fields based on a power law model and idealized atmospheric vertical profiles as initial conditions, three sets of large‐eddy simulations provide insight in the influence of soil moisture heterogeneity on the ensuing growth of the convective boundary layer and development of shallow cumulus clouds. A sensitivity on the strong, weak, and unstructured soil moisture heterogeneity is investigated. The simulation results show that domain‐averaged land surface sensible heat and latent heat flux change strongly with changing soil moisture variance because of the interactions between surface heterogeneity and induced circulations, while domain means of soil moisture are identical. Vertical profiles of boundary layer characteristics are strongly influenced by the surface energy partitioning and induced circulations, especially the profiles of liquid water and liquid water flux. The amount of liquid water and liquid water flux increases with increasing structure. In addition, the liquid water path is higher in case of strongly‐structured heterogeneity because more available energy is partitioned into latent heat and more intensive updrafts exist. Interestingly, the increase of liquid water path with increasing soil moisture variance only occurs in the strongly structured cases, which suggests that soil moisture variance and structure work conjunctively in the surface energy partitioning and the cloud formation

Sensitivity of cloud-phase distribution to cloud microphysics and thermodynamics in simulated deep convective clouds and SEVIRI retrievals

The formation of ice in clouds is an important process in mixed-phase clouds, and the radiative properties and dynamical developments of clouds strongly depend on their partitioning between the liquid and ice phases. In this study, we investigated the sensitivities of the cloud phase to the ice-nucleating particle (INP) concentration and thermodynamics. Moreover, passive satellite retrieval algorithms and cloud products were evaluated to identify whether they could detect cloud microphysical and thermodynamical perturbations. Experiments were conducted using the ICOsahedral Nonhydrostatic (ICON) model at the convection-permitting resolution of about 1.2 km on a domain covering significant parts of central Europe, and they were compared to two different retrieval products based on Spinning Enhanced Visible and InfraRed Imager (SEVIRI) measurements. We selected a day with multiple isolated deep convective clouds, reaching a homogeneous freezing temperature at the cloud top. The simulated cloud liquid pixel fractions were found to decrease with increasing INP concentration, both within clouds and at the cloud top. The decrease in the cloud liquid pixel fraction was not monotonic and was stronger in high-INP cases. Cloud-top glaciation temperatures shifted toward warmer temperatures with an increasing INP concentration by as much as 8$^◦$C

Long-term variations in actual evapotranspiration over the Tibetan Plateau

Actual terrestrial evapotranspiration (ET$_{a}$) is a key parameter controlling land–atmosphere interaction processes and water cycle. However, spatial distribution and temporal changes in ET$_{a}$ over the Tibetan Plateau (TP) remain very uncertain. Here we estimate the multiyear (2001–2018) monthly ET$_{a}$ and its spatial distribution on the TP by a combination of meteorological data and satellite products. Validation against data from six eddy-covariance monitoring sites yielded root-mean-square errors ranging from 9.3 to 14.5 mm per month and correlation coefficients exceeding 0.9. The domain mean of annual ET$_{a}$ on the TP decreased slightly (−1.45 mm yr$^{-1}$, p90° E) but decreased significantly at a rate of −5.52 mm yr$^{-1}$ (p<0.05) in the western sector of the TP (long <90° E). In addition, the decreases in annual ET$_{a}$ were pronounced in the spring and summer seasons, while almost no trends were detected in the autumn and winter seasons. The mean annual ET$_{a}$ during 2001–2018 and over the whole TP was 496±23 mm. Thus, the total evapotranspiration from the terrestrial surface of the TP was 1238.3±57.6 km3 yr$^{-1}$. The estimated ET$_{a}$ product presented in this study is useful for an improved understanding of changes in energy and water cycle on the TP. The dataset is freely available at the Science Data Bank (https://doi.org/10.11922/sciencedb.t00000.00010; Han et al., 2020b) and at the National Tibetan Plateau Data Center (https://doi.org/10.11888/Hydro.tpdc.270995, Han et al., 2020a)

Monitoring Water and Energy Cycles at Climate Scale in the Third Pole Environment (CLIMATE-TPE)

A better understanding of the water and energy cycles at climate scale in the Third Pole Environment is essential for assessing and understanding the causes of changes in the cryosphere and hydrosphere in relation to changes of plateau atmosphere in the Asian monsoon system and for predicting the possible changes in water resources in South and East Asia. This paper reports the following results: (1) A platform of in situ observation stations is briefly described for quantifying the interactions in hydrosphere-pedosphere-atmosphere-cryosphere-biosphere over the Tibetan Plateau. (2) A multiyear in situ L-Band microwave radiometry of land surface processes is used to develop a new microwave radiative transfer modeling system. This new system improves the modeling of brightness temperature in both horizontal and vertical polarization. (3) A multiyear (2001–2018) monthly terrestrial actual evapotranspiration and its spatial distribution on the Tibetan Plateau is generated using the surface energy balance system (SEBS) forced by a combination of meteorological and satellite data. (4) A comparison of four large scale soil moisture products to in situ measurements is presented. (5) The trajectory of water vapor transport in the canyon area of Southeast Tibet in different seasons is analyzed, and (6) the vertical water vapor exchange between the upper troposphere and the lower stratosphere in different seasons is presented

Large‐eddy simulation of catchment‐scale circulation

The impact of soil moisture heterogeneity on the convective boundary layer (CBL) development was studied. Based on results from large‐eddy simulation (LES) applying soil moisture patterns along a river corridor and idealized atmospheric vertical profiles as initial conditions, this study provides insight in the influence of spatial scale of soil moisture heterogeneity on catchment‐scale circulations (CCs) and the ensuing growth of the CBL. The simulation results show that the intensity of organized circulations resulting from soil moisture heterogeneity is nonlinearly dependent upon soil moisture heterogeneity scale λ (SMHS) and horizontal gradient. Because of the large SMHS and strong soil moisture contrast, none of the simulations has reached a true steady state even after 24 h of simulation time. The intensity of organized circulations shows a sigmoidal dependence on SMHS. The optimal SMHS for horizontal transport is on the order of 19.2 km, while optimal SMHS for vertical motions occurs at 2.4 km. In these cases, the CCs also exert a strong influence on the boundary‐layer structure and the entrainment layer. The potential temperature is not constant with height due to a weak mixing in the boundary layer for large SMHS cases. Differences in sensible heat flux profiles between the heterogeneous cases increase with increasing height and reach a maximum at the top of the CBL. Interestingly, boundary‐layer height changes strongly with changing horizontal soil moisture gradient and SMHS while domain means, variances, and amplitudes of land surface energy fluxes are all almost identical. The entrainment flux and subsidence at the top of the CBL are jointly responsible for the CBL height variation

Annual and Seasonal Trends of Vegetation Responses and Feedback to Temperature on the Tibetan Plateau since the 1980s

The vegetation–temperature relationship is crucial in understanding land–atmosphere interactions on the Tibetan Plateau. Although many studies have investigated the connections between vegetation and climate variables in this region using remote sensing technology, there remain notable gaps in our understanding of vegetation–temperature interactions over different timescales. Here, we combined site-level air temperature observations, information from the global inventory modeling and mapping studies (GIMMS) dataset, and moderate-resolution imaging spectroradiometer (MODIS) products to analyze the spatial and temporal patterns of air temperature, vegetation, and land surface temperature (LST) on the Tibetan Plateau at annual and seasonal scales. We achieved these spatiotemporal patterns by using Sen’s slope, sequential Mann–Kendall tests, and partial correlation analysis. The timescale differences of vegetation-induced LST were subsequently discussed. Our results demonstrate that a breakpoint of air temperature change occurred on the Tibetan Plateau during 1994–1998, dividing the study period (1982–2013) into two phases. A more significant greening response of NDVI occurred in the spring and autumn with earlier breakpoints and a more sensitive NDVI response occurred in recent warming phase. Both MODIS and GIMMS data showed a common increase in the normalized difference vegetation index (NDVI) on the Tibetan Plateau for all timescales, while the former had a larger greening area since 2000. The most prominent trends in NDVI and LST were identified in spring and autumn, respectively, and the largest areas of significant variation in NDVI and LST mostly occurred in winter and autumn, respectively. The partial correlation analysis revealed a significant negative impact of NDVI on LST during the annual scale and autumn, and it had a significant positive impact during spring. Our findings improve the general understanding of vegetation–climate relationships at annual and seasonal scales

Trends of land surface heat fluxes on the Tibetan Plateau from 2001 to 2012

A parameterization approach of effective roughness length was introduced into the Surface Energy Balance System (SEBS) model to account for subgrid-scale topographical influences. Regional distribution of land surface heat flux values (including net radiation flux, ground heat flux, sensible heat flux, and latent heat flux) was estimated on the Tibetan Plateau (TP) based on the SEBS model, and utilizing remote sensing products and reanalysis datasets. We then investigated annual trends in these fluxes for the period 2001–2012. It was found that land surface net radiation flux increased slightly, especially in high, mountainous regions and the central TP, and was influenced by glacial retreat and topsoil wetting, respectively. Sensible heat flux decreased overall, especially in the central and northern TP. In the Yarlung Zangbo River (YZR) Basin, the sensible heat flux increased because of a rise in the ground-air temperature difference. The latent heat flux increased over the majority TP, except for areas in the YZR Basin. This can be attributed to increases in precipitation and vegetation greening

Sensitivity of cloud-phase distribution to cloud microphysics and thermodynamics in simulated deep convective clouds and SEVIRI retrievals

International audienceThe formation of ice in clouds is an important process in mixed-phase clouds, and the radiative properties and dynamical developments of clouds strongly depend on their partitioning between the liquid and ice phases. In this study, we investigated the sensitivities of the cloud phase to the ice-nucleating particle (INP) concentration and thermodynamics. Moreover, passive satellite retrieval algorithms and cloud products were evaluated to identify whether they could detect cloud microphysical and thermodynamical perturbations. Experiments were conducted using the ICOsahedral Nonhydrostatic (ICON) model at the convection-permitting resolution of about 1.2 km on a domain covering significant parts of central Europe, and they were compared to two different retrieval products based on Spinning Enhanced Visible and InfraRed Imager (SEVIRI) measurements. We selected a day with multiple isolated deep convective clouds, reaching a homogeneous freezing temperature at the cloud top. The simulated cloud liquid pixel fractions were found to decrease with increasing INP concentration, both within clouds and at the cloud top. The decrease in the cloud liquid pixel fraction was not monotonic and was stronger in high-INP cases. Cloud-top glaciation temperatures shifted toward warmer temperatures with an increasing INP concentration by as much as 8˚C. Moreover, the impact of the INP concentration on cloud-phase partitioning was more pronounced at the cloud top than within the cloud. Furthermore, initial and lateral boundary temperature fields were perturbed with increasing and decreasing temperature increments from 0 to ±3 and ±5 K between 3 and 12 km, respectively. Perturbing the initial thermodynamic state was also found to systematically affect the cloud-phase distribution. However, the simulated cloud-top liquid pixel fraction, diagnosed using radiative transfer simulations as input to a satellite forward operator and two different satellite remote-sensing retrieval algorithms, deviated from one of the satellite products regardless of perturbations in the INP concentration or the initial thermodynamic state for warmer subzero temperatures while agreeing with the other retrieval scheme much better, in particular for the high-INP and high-CAPE (convective available potential energy) scenarios. Perturbing the initial thermodynamic state, which artificially increases the instability of the mid- and upper-troposphere, brought the simulated cloud-top liquid pixel fraction closer to the satellite observations, especially in the warmer mixed-phase temperature range

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

Dominant Modes of Tibetan Plateau Summer Surface Sensible Heating and Associated Atmospheric Circulation Anomalies

Based on empirical orthogonal function (EOF) analysis, the dominant modes of variations in summer surface sensible heating (SH) over the Tibetan Plateau (TP), as well as the associated atmospheric circulation anomalies, were investigated in this study. The results show that the first dominant mode of summer SH presents a feature of decadal reduction over the whole TP on an interdecadal time scale, and the second dominant mode is characterized by a zonally asymmetric pattern with positive (negative) SH anomalies in the western (eastern) TP on an interannual time scale. The variations of summer SH are dominated by anomalies in downwelling surface shortwave radiation (DSWR), which are associated with atmospheric circulation changes. The first dominant mode of variation in SH is connected to the interdecadal variation of the Silk Road Pattern (SRP). Further analysis reveals that the interdecadal phase shift of the SRP induces anticyclone circulation to the northeast of the TP, leading to enhanced water vapor supply and convergence over the TP. This can lead to an increase in the total cloud cover, and a reduction in DSWR, contributing to the decadal reduction in SH over the TP. The second dominant mode of variation in SH is related to a stationary teleconnection pattern over the Eurasian continent named the North Atlantic-East and North Asia pattern (NAENA). Corresponding to the positive phase of the NAENA, there is a cyclone anomaly to the west TP, leading to anomalous water vapor convergence (divergence) over the eastern (western) TP. This can result in enhanced (decreased) cloud cover, reduced (increased) DSWR, and therefore, an anomalous decrease (enhancement) in SH over the east (west) of the TP. Furthermore, the southwesterly wind anomaly, which is accompanied by the anomalous cyclone to the west TP, leads to positive SH in the western TP