Retrieving air humidity,global solar radiation,and reference evapotranspiration from daily temperatures: development and validation of new methods for Mexico. Part II: radiation

We propose a new model to estimate daily global radiation from daily temperature range measurements. This model combines that of Majumdar et al. (Sol Energy 13(4):383–394, 1972) to estimate clear sky radiation with a Gompertz function to estimate the relation between temperature range and cloud transmittance. Model parameters are estimated from historical weather data: maximum and minimum temperatures and, if available, relative humidity; no other calibration is required. The model was parametrized and validated using 788 weather stations in Mexico. When calibrated using historical humidity data, daily global radiation was estimated with a mean root mean square error of 3.06 MJ m^?2 day^?1. The model performed well in all situations, except for a few stations around the Gulf of Mexico and in mountain areas. When using estimated humidity, the root mean square error of prediction was only slightly degraded (3.07 MJ m^?2 day^?1). Possible theoretical basis and applicability of this model to other environments are discussed.     

Retrieving air humidity,global solar radiation,and reference evapotranspiration from daily temperatures: development and validation of new methods for Mexico. Part I: humidity

The occurrence of flood and drought frequency is highly correlated with the temporal fluctuations of streamflow series; understanding of these fluctuations is essential for the improved modeling and statistical prediction of extreme changes in river basins. In this study, the complexity of daily streamflow fluctuations was investigated by using multifractal detrended fluctuation analysis (MF-DFA) in a large heterogeneous lake basin, the Poyang Lake basin in China, and the potential impacts of human activities were also explored. Major results indicate that the multifractality of streamflow fluctuations shows significant regional characteristics. In the study catchment, all the daily streamflow series present a strong long-range correlation with Hurst exponents bigger than 0.8. The q-order Hurst exponent h(q) of all the hydrostations can be characterized well by only two parameters: a (0.354 ≤ a ≤ 0.384) and b (0.627 ≤ b ≤ 0.677), with no pronounced differences. Singularity spectrum analysis pointed out that small fluctuations play a dominant role in all daily streamflow series. Our research also revealed that both the correlation properties and the broad probability density function (PDF) of hydrological series can be responsible for the multifractality of streamflow series that depends on watershed areas. In addition, we emphasized the relationship between watershed area and the estimated multifractal parameters, such as the Hurst exponent and fitted parameters a and b from the q-order Hurst exponent h(q). However, the relationship between the width of the singularity spectrum (Δα) and watershed area is not clear. Further investigation revealed that increasing forest coverage and reservoir storage can effectively enhance the persistence of daily streamflow, decrease the hydrological complexity of large fluctuations, and increase the small fluctuations.     

Estimation of Surface-Layer Scaling Parameters in the Unstable Boundary Layer: Implications for Orographic Flow Speed-Up

A method is proposed for estimating the surface-layer depth \((z_s)\) and the friction velocity \((u_*)\) as a function of stability (here quantified by the Obukhov length, L) over the complete range of unstable flow regimes. This method extends that developed previously for stable conditions by Argaín et al. (Boundary-Layer Meteorol 130:15–28, 2009), but uses a qualitatively different approach. The method is specifically used to calculate the fractional speed-up \((\varDelta S)\) in flow over a ridge, although it is suitable for more general boundary-layer applications. The behaviour of \(z_s \left( L\right) \) and \(u_*\left( L\right) \) as a function of L is indirectly assessed via calculation of \(\varDelta S\left( L\right) \) using the linear model of Hunt et al. (Q J R Meteorol Soc 29:16–26, 1988) and its comparison with the field measurements reported in Coppin et al. (Boundary-Layer Meteorol 69:173–199, 1994) and with numerical simulations carried out using a non-linear numerical model, FLEX. The behaviour of \(\varDelta S\) estimated from the linear model is clearly improved when \(u_*\) is calculated using the method proposed here, confirming the importance of accounting for the dependences of \(z_s\left( L \right) \) and \(u_*\left( L \right) \) on L to better represent processes in the unstable boundary layer.     

New alternatives for reference evapotranspiration estimation in West Africa using limited weather data and ancillary data supply strategies.

Accurate estimation of reference evapotranspiration (ET ₀ ) is essential for the computation of crop water requirements, irrigation scheduling, and water resources management. In this context, having a battery of alternative local calibrated ET ₀ estimation methods is of great interest for any irrigation advisory service. The development of irrigation advisory services will be a major breakthrough for West African agriculture. In the case of many West African countries, the high number of meteorological inputs required by the Penman-Monteith equation has been indicated as constraining. The present paper investigates for the first time in Ghana, the estimation ability of artificial intelligence-based models (Artificial Neural Networks (ANNs) and Gene Expression Programing (GEPs)), and ancillary/external approaches for modeling reference evapotranspiration (ET ₀ ) using limited weather data. According to the results of this study, GEPs have emerged as a very interesting alternative for ET ₀ estimation at all the locations of Ghana which have been evaluated in this study under different scenarios of meteorological data availability. The adoption of ancillary/external approaches has been also successful, moreover in the southern locations. The interesting results obtained in this study using GEPs and some ancillary approaches could be a reference for future studies about ET ₀ estimation in West Africa.     

Crop modelling: towards locally relevant and climate-informed adaptation

A gap between the potential and practical realisation of adaptation exists: adaptation strategies need to be both climate-informed and locally relevant to be viable. Place-based approaches study local and contemporary dynamics of the agricultural system, whereas climate impact modelling simulates climate-crop interactions across temporal and spatial scales. Crop-climate modelling and place-based research on adaptation were strategically reviewed and analysed to identify areas of commonality, differences, and potential learning opportunities to enhance the relevance of both disciplines through interdisciplinary approaches. Crop-modelling studies have projected a 7–15% mean yield change with adaptation compared to a non-adaptation baseline (Nature Climate Change 4:1–5, 2014). Of the 17 types of adaptation strategy identified in this study as place-based adaptations occurring within Central America, only five were represented in crop-climate modelling literature, and these were as follows: fertiliser, irrigation, change in planting date, change in cultivar and area cultivated. The breath and agency of real-life adaptation compared to its representation in modelling studies is a source of error in climate impact simulations. Conversely, adaptation research that omits assessment of future climate variability and impact does not enable to provide sustainable adaptation strategies to local communities so risk maladaptation. Integrated and participatory methods can identify and reduce these sources of uncertainty, for example, stakeholder’s engagement can identify locally relevant adaptation pathways. We propose a research agenda that uses methodological approaches from both the modelling and place-based approaches to work towards climate-informed locally relevant adaptation.     

Relationship Between Fine-Particle Pollution and the Urban Heat Island in Beijing,China: Observational Evidence

Urbanization has led to a significant urban heat island (UHI) effect in Beijing in recent years. At the same time, air pollution caused by a large number of fine particles significantly influences the atmospheric environment, urban climate, and human health. The distribution of fine particulate matter (PM_2.5) concentration and its relationship with the UHI effect in the Beijing area are analyzed based on station-observed hourly data from 2012 to 2016. We conclude that, (1) in the last five years, the surface concentrations of PM_2.5 averaged for urban and rural sites in and around Beijing are 63.2 and 40.7 µg m^?3, respectively, with significant differences between urban and rural sites (ΔPM_2.5) at the seasonal, monthly and daily scales observed; (2) there is a large correlation between ΔPM_2.5 and the UHI intensity defined as the differences in the mean (ΔT_ave), minimum (ΔT_min), and maximum (ΔT_max) temperatures between urban and rural sites. The correlation between ΔPM_2.5 and ΔT_min (ΔT_max) is the highest (lowest); (3) a Granger causality analysis further shows that ΔPM_2.5 and ΔT_min are most correlated for a lag of 1–2 days, while the correlation between ΔPM_2.5 and ΔT_ave is lower; there is no causal relationship between ΔPM_2.5 and ΔT_max; (4) a case analysis shows that downwards shortwave radiation at the surface decreases with an increase in PM_2.5 concentration, leading to a weaker UHI intensity during the daytime. During the night, the outgoing longwave radiation from the surface decreases due to the presence of daytime pollutants, the net effect of which is a slower cooling rate during the night in cities than in the suburbs, leading to a larger ΔT_min.     

Modeling monthly meteorological and agronomic frost days,based on minimum air temperature,in Center-Southern Brazil

Although Brazil is predominantly a tropical country, frosts are observed with relative high frequency in the Center-Southern states of the country, affecting mainly agriculture, forestry, and human activities. Therefore, information about the frost climatology is of high importance for planning of these activities. Based on that, the aims of the present study were to develop monthly meteorological (F _MET) and agronomic (F _AGR) frost day models, based on minimum shelter air temperature (T _MN), in order to characterize the temporal and spatial frost days variability in Center-Southern Brazil. Daily minimum air temperature data from 244 weather stations distributed across the study area were used, being 195 for developing the models and 49 for validating them. Multivariate regression models were obtained to estimate the monthly T _MN, once the frost day models were based on this variable. All T _MN regression models were statistically significant (p < 0.001), presenting adjusted R ² between 0.69 and 0.90. Center-Southern Brazil is mainly hit by frosts from mid-fall (April) to mid-spring (October). The period from November to March is considered as frost-free, being very rare a frost day within that period. Monthly F _MET and F _AGR presented significant sigmoidal relationships with T _MN (p < 0.0001), with adjusted R ² above of 0.82. The residuals of the frost day models were random, which means that the sigmoidal models performed quite well for interpreting the frost day variability throughout the study area. The highlands of Santa Catarina, Rio Grande do Sul, São Paulo, and Minas Gerais had in average more than 25 and 13 frosts per year, respectively, for F _MET and F _AGR. The F _MET and F _AGR maps developed in this study for Center-Southern Brazil is a useful tool for farmers, foresters, and researchers, since they contribute to reduce frost spatial and temporal uncertainty, helping in planning project for strategic purposes. Furthermore, the monthly F _MET and F _AGR maps for this Brazilian region are the first zoning of these variables for the country.     

A compact form for the analytic description of temperature dependence of saturation vapor pressure over plane surfaces of water and ice

A possibility is studied of extending the range of action of the simple three-parameter formula (ITS-90 scale) proposed in the previous work of the author [2] for the dependence of saturation vapor pressure E on temperature T within the range of 250 to 490 K. The results demonstrated that the dependence ln[E(T)/E(T _bas)] = (T - T _bas)[A - B(T - T _bas) + C(T - T _bas)²]/T with four sets of coefficients A, B, and C obtained using one base temperature Tbas equal to the temperature of triple point of water T _t = 273.16 K and two additional base values T _bas2 = 473.16 K and T _bas3 = 623.16 K makes it possible to approximate rather accurately the initial experimental and computed data in the temperature range from the point of homogeneous freezing of 235 K to the critical temperature of 647 K for liquid water and from 193 K to T _t for ice. A procedure used for obtaining the inverse function T(E) by solving the third-degree algebraic equation is validated. A hypothesis is proposed for the physical substantiation of additional base points in the form of “a noticeable appearance of dimers at the point T _bas2 and their 100% concentration at the temperature T _bas3.”     

Effects of Unstable Thermal Stratification on Vertical Fluxes of Heat and Momentum in Urban Areas

A Reynolds-averaged Navier–Stokes microscale model is used for the simulation of the effect of unstable thermal stratification on the flow within an aligned configuration of building-like cubes as used in Santiago et al. (Urban Clim 9:115–133, 2014). The spatially-averaged results show increased dispersive fluxes, turbulent length scales and sectional drag coefficient. An extension of K-theory is presented to parametrize the sum of the turbulent and dispersive fluxes, and the length scale and drag coefficient increases are parametrized as functions of the ratio of buoyant and inertial forces. This approach improves the results of urban canopy parametrization simulations inside and above the urban canyon and represents the first attempt to account for the dispersive fluxes and the effect of solar radiation on the flow.     

Characteristics of the precipitation recycling ratio and its relationship with regional precipitation in China

A dynamic recycling model (DRM) with an analytical moisture trajectory tracking method, together with Japan Meteorological Agency 25-year reanalysis data, is used to study the regional precipitation recycling process across China, by calculating the regional recycling ratio (ρ _r ) at the daily time scale during 1979–2010. The distribution of ρ _r shows that, in western China, especially the Tibetan Plateau and its surrounding areas, precipitation is strongly dependent on the recycling process associated with regional evaporation. In Southeast China, however, the contribution from the recycling processes is much smaller due to the influence of the summer monsoon. A precipitation threshold value of about 4 mm/day is obtained from detailed analysis of both extreme and all-range ρ _r years. According to this threshold, China is classified into three types of sub-regions: low-precipitation sub-regions (mainly in the northwest), high-precipitation sub-regions (mainly in the south), and medium-precipitation sub-regions (mainly in the northeast). It is found that ρ _r correlates positively with precipitation, as well as convective precipitation (P _CP) and large-scale precipitation (P _LP) in the low-precipitation sub-regions. However, negative ρ _r ?～?P _LP correlations are found in the high-precipitation sub-regions and nonsignificant correlations exist in the medium-precipitation sub-regions. As P _CP is mainly locally generated due to mid-latitude mesoscale systems and the cumulus parameterization used in producing the reanalysis, the recycling ratio positively correlates to the ratio P _CP/P _LP in almost all sub-regions, particularly in the Tibetan Plateau and its surrounding areas. The correlation between radiation flux and ρ _r suggests more net radiation supports more evaporation and higher ρ _r , especially in the high-precipitation sub-regions. The influence of clouds on shortwave radiation is crucial, since evaporation is suppressed when the amount of cloudiness increases, especially in the high-precipitation sub-regions. Together with the consideration of soil moisture, it can be inferred that limited soil moisture inhibits evaporation in the low-precipitation sub-regions, while the energy or radiation is the dominant factor controlling evaporation in the high-precipitation sub-regions.     

Historical and experimental evidence for enhanced concentration of artemesinin,a global anti-malarial treatment,with recent and projected increases in atmospheric carbon dioxide

Although the role of rising atmospheric carbon dioxide concentration [CO₂] on plant growth and fecundity is widely acknowledged as important within the scientific community; less research is available regarding the impact of [CO₂] on secondary plant compounds, even though such compounds can play a significant role in human health. At present, Artemisia annua, an annual plant species native to China, is widely recognized as the primary source of artemesinin used in artemesinin combination therapies or ACTs. ACTs, in turn, are used globally for the treatment of simple Plasmodium falciparum malaria, the predominant form of malaria in Africa. In this study, artemesinin concentration was quantified for multiple A. annua populations in China using a free-air CO₂ enrichment (FACE) system as a function of [CO₂]-induced changes both in situ and as a function of the foliar ratio of carbon to nitrogen (C:N). The high correlation between artemesinin concentration and C:N allowed an historical examination of A. annua leaves collected at 236 locations throughout China from 1905 through 2009. Both the historical and experimental data indicate that increases in artemesinin foliar concentration are likely to continue in parallel with the ongoing increase in atmospheric [CO₂]. The basis for the [CO₂]-induced increase in artemesinin is unclear, but could be related to the carbon: nutrient hypothesis of Bryant et al. (1983). Overall, these data provide the first evidence that historic and projected increases in atmospheric [CO₂] may be associated with global changes in artemesinin chemistry, potentially allowing a greater quantity of drug available for the same area of cultivation.     

Comprehensive modeling of monthly mean soil temperature using multivariate adaptive regression splines and support vector machine

Soil temperature (T _s) and its thermal regime are the most important factors in plant growth, biological activities, and water movement in soil. Due to scarcity of the T _s data, estimation of soil temperature is an important issue in different fields of sciences. The main objective of the present study is to investigate the accuracy of multivariate adaptive regression splines (MARS) and support vector machine (SVM) methods for estimating the T _s. For this aim, the monthly mean data of the T _s (at depths of 5, 10, 50, and 100 cm) and meteorological parameters of 30 synoptic stations in Iran were utilized. To develop the MARS and SVM models, various combinations of minimum, maximum, and mean air temperatures (T _min, T _max, T); actual and maximum possible sunshine duration; sunshine duration ratio (n, N, n/N); actual, net, and extraterrestrial solar radiation data (R _s, R _n, R _a); precipitation (P); relative humidity (RH); wind speed at 2 m height (u ₂); and water vapor pressure (V_p) were used as input variables. Three error statistics including root-mean-square-error (RMSE), mean absolute error (MAE), and determination coefficient (R ²) were used to check the performance of MARS and SVM models. The results indicated that the MARS was superior to the SVM at different depths. In the test and validation phases, the most accurate estimations for the MARS were obtained at the depth of 10 cm for T _max, T _min, T inputs (RMSE = 0.71 °C, MAE = 0.54 °C, and R ² = 0.995) and for RH, V _p, P, and u ₂ inputs (RMSE = 0.80 °C, MAE = 0.61 °C, and R ² = 0.996), respectively.     

Preparation of frost atlas using different interpolation methods in a semiarid region of south of Iran

In this research, suitability of different kriging and inverse distance weighted (IDW) methods in estimating occurrence date of frost was evaluated. Data included minimum daily air temperature values from 27 meteorological stations of Fars province in southern Iran from 18 to 45 years. Data ranges of 0 to ?1.5, ?1.5 to ?3 and below ?3°C were considered as mild, moderate and severe frost intensities, respectively. Starting with the first day of autumn, iso-occurrence days for the frost intensities and occurrence probabilities (25%, 50%, 75% and 90%) were estimated using ordinary kriging, cokriging, residual kriging type 1 (RK1), residual kriging type 2 (RK2), universal kriging and IDW methods. In these models, the errors of estimated frost intensities at different probabilities were lowest in the RK2 model, but lack of establishment of spatial structure due to long distance between stations caused the predictions not to be acceptable in some cases. In a proposed method (modified inverse distance weighted, MIDW), the trend between the first and last days of frost occurrence with earth elevation was removed, and the reminder values were estimated by (IDW) method. Although, the errors for estimated frost dates by MIDW and RK2 methods were the same, but the MIDW method did not have the spatial establishment shortcoming. Furthermore, the simplicity and practicality of the MIDW method makes it a reasonable selection.     

Response to commentary by J. L. Bamber,W. P. Aspinall and R. M. Cooke (2016)

In a commentary paper, Bamber et al. (Nat Clim Change 3:424–427, 2016) respond to our recent assessment (De Vries and Van de Wal Clim Change 1–14, 2015) of their expert judgment based study on projections of future sea level rise due to the melting of the large ice sheets (Bamber and Aspinall Nat Clim Chang 3:424–427, 2013). In this response we comment on their remarks.     

Geo-spatial analysis of temporal trends of temperature and its extremes over India using daily gridded (1°×1°) temperature data of 1969–2005

Daily gridded (1°×1°) temperature data (1969–2005) were used to detect spatial patterns of temporal trends of maximum and minimum temperature (monthly and seasonal), growing degree days (GDDs) over the crop-growing season (kharif, rabi, and zaid) and annual frequencies of temperature extremes over India. The direction and magnitude of trends, at each grid level, were estimated using the Mann–Kendall statistics (α = 0.05) and further assessed at the homogeneous temperature regions using a field significance test (α=0.05). General warming trends were observed over India with considerable variations in direction and magnitude over space and time. The spatial extent and the magnitude of the increasing trends of minimum temperature (0.02–0.04 °C year^?1) were found to be higher than that of maximum temperature (0.01–0.02 °C year^?1) during winter and pre-monsoon seasons. Significant negative trends of minimum temperature were found over eastern India during the monsoon months. Such trends were also observed for the maximum temperature over northern and eastern parts, particularly in the winter month of January. The general warming patterns also changed the thermal environment of the crop-growing season causing significant increase in GDDs during kharif and rabi seasons across India. The warming climate has also caused significant increase in occurrences of hot extremes such as hot days and hot nights, and significant decrease in cold extremes such as cold days and cold nights.     

Uncertainty of global summer precipitation in the CMIP5 models: a comparison between high-resolution and low-resolution models

Worldwide, the majority of rapidly growing neighborhoods are found in the Global South. They often exhibit different building construction and development patterns than the Global North, and urban climate research in many such neighborhoods has to date been sparse. This study presents local-scale observations of net radiation (Q _* ) and sensible heat flux (Q _H ) from a lightweight low-rise neighborhood in the desert climate of Andacollo, Chile, and compares observations with results from a process-based urban energy-balance model (TUF3D) and a local-scale empirical model (LUMPS) for a 14-day period in autumn 2009. This is a unique neighborhood-climate combination in the urban energy-balance literature, and results show good agreement between observations and models for Q _* and Q _H . The unmeasured latent heat flux (Q _E ) is modeled with an updated version of TUF3D and two versions of LUMPS (a forward and inverse application). Both LUMPS implementations predict slightly higher Q _E than TUF3D, which may indicate a bias in LUMPS parameters towards mid-latitude, non-desert climates. Overall, the energy balance is dominated by sensible and storage heat fluxes with mean daytime Bowen ratios of 2.57 (observed Q _H /LUMPS Q _E )–3.46 (TUF3D). Storage heat flux (ΔQ _S ) is modeled with TUF3D, the empirical objective hysteresis model (OHM), and the inverse LUMPS implementation. Agreement between models is generally good; the OHM-predicted diurnal cycle deviates somewhat relative to the other two models, likely because OHM coefficients are not specified for the roof and wall construction materials found in this neighborhood. New facet-scale and local-scale OHM coefficients are developed based on modeled ΔQ _S and observed Q _* . Coefficients in the empirical models OHM and LUMPS are derived from observations in primarily non-desert climates in European/North American neighborhoods and must be updated as measurements in lightweight low-rise (and other) neighborhoods in various climates become available.     

Analysis of photosynthetically active radiation and applied parameterization model for estimating of PAR in the North China Plain

As photosynthetically active radiation (PAR) variability and PAR estimating methods play an important role in climate change and ecological process research, PAR variation trends and broadband global solar radiation (R _s ) ratios (PAR/R _s ) in the North China Plain (NCP) are examined using in situ PAR and R _s observed data for 2005 to 2011. The annual average PAR value found in the NCP is 22.9 mol m^?2 d^?1. The highest and lowest values were recorded at Changwu and Luancheng sites, respectively. The highest PAR/R _s value was found in Jiaozhouwan due to large water vapor volumes present in this area. PAR/R _s levels have increased in the NCP due to a decrease in fine aerosols and increase in water vapor concentration. From these analysis results, a parameterization model that can be applied to all sky conditions was checked. Empirical estimation model comparisons for obtaining PAR values indicate that model was least accurate when R _s was used independently. When the model included R _s, the clearness index (K _s) and the solar zenith angle, the model estimated PAR values with acceptable accuracy. A parameterization model was constructed by considering K _s and attenuation factors of PAR under clear weather conditions (ρ _clear). The improved parameterization model more accurately predicts values for local sites and for various observation sites.     

Effect of solubility limitation on hygroscopic growth and cloud drop activation of SOA particles produced from traffic exhausts

Hygroscopicity measurements of secondary organic aerosol (SOA) particles often show inconsistent results between the supersaturated and subsaturated regimes, with higher activity as cloud condensation nucleus (CCN) than indicated by hygroscopic growth. In this study, we have investigated the discrepancy between the two regimes in the Lund University (LU) smog chamber. Various anthropogenic SOA were produced from mixtures of different precursors: anthropogenic light aromatic precursors (toluene and m-xylene), exhaust from a diesel passenger vehicle spiked with the light aromatic precursors, and exhaust from two different gasoline-powered passenger vehicles. Three types of seed particles were used: soot aggregates from a diesel vehicle, soot aggregates from a flame soot generator and ammonium sulphate (AS) particles. The hygroscopicity of seed particles with condensed, photochemically produced, anthropogenic SOA was investigated with respect to critical supersaturation (s_c) and hygroscopic growth factor (gf) at 90% relative humidity. The hygroscopicity parameter κ was calculated for the two regimes: κ_sc and κ_gf, from measurements of s_c and gf, respectively. The two κ showed significant discrepancies, with a κ_gf /κ_sc ratio closest to one for the gasoline experiments with ammonium sulphate seed and lower for the soot seed experiments. Empirical observations of s_c and gf were compared to theoretical predictions, using modified Köhler theory where water solubility limitations were taken into account. The results indicate that the inconsistency between measurements in the subsaturated and supersaturated regimes may be explained by part of the organic material in the particles produced from anthropogenic precursors having a limited solubility in water.     

Projected Changes in Extreme High Temperature and Heat Stress in China

High temperature accompanied with high humidity may result in unbearable and oppressive weather. In this study, future changes of extreme high temperature and heat stress in mainland China are examined based on daily maximum temperature (T_x) and daily maximum wet-bulb globe temperature (T_w). T_w has integrated the effects of both temperature and humidity. Future climate projections are derived from the bias-corrected climate data of five general circulation models under the Representative Concentration Pathways (RCPs) 2.6 and 8.5 scenarios. Changes of hot days and heat waves in July and August in the future (particularly for 2020–50 and 2070–99), relative to the baseline period (1981–2010), are estimated and analyzed. The results show that the future T_x and T_w of entire China will increase by 1.5–5°C on average around 2085 under different RCPs. Future increases in T_x and T_w exhibit high spatial heterogeneity, ranging from 1.2 to 6°C across different regions and RCPs. By around 2085, the mean duration of heat waves will increase by 5 days per annum under RCP8.5. According to T_x, heat waves will mostly occur in Northwest and Southeast China, whereas based on T_w estimates, heat waves will mostly occur over Southeast China and the mean heat wave duration will be much longer than those from T_x. The total extreme hot days (T_x or T_w > 35°C) will increase by 10–30 days. Southeast China will experience the severest heat stress in the near future as extreme high temperature and heat waves will occur more often in this region, which is particularly true when heat waves are assessed based on T_w. In comparison to those purely temperature-based indices, the index T_w provides a new perspective for heat stress assessment in China.     

Atmospheric circulation epochs and climate changes

The atmospheric circulation studies allow climate changes to be diagnosed and forecasted. Variations in occurrence frequencies of the atmospheric circulation forms W, E, and C (by the Vangengeim classification) and Z, M ₁, and M ₂ (by the Girs classification), which characterize climatic conditions in most of the Northern Hemisphere, are analyzed over a period of more than 100 years. It is shown that the occurrence frequency of the forms W, C, and M ₁ continually decreased, while that of the forms E and Z increased, which indicates a significant change in atmospheric circulation in the Northern Hemisphere during the last century. The occurrence frequency of the forms C and Z demonstrates specific features at inter-decade time scales. Correlations are found between accumulated sums of anomalies of occurrence frequencies of the atmospheric circulation forms C, (W + E), Z, and (M ₁ + M ₂) and inter-decade variations of the Earth’s rotation. The causes of these relationships are discussed along with possibilities of their use for diagnosis of climatic variations in the Northern Hemisphere.