Multivariate drought frequency estimation using copula method in Southwest China

This study reveals the impacts of climatic variable trends on drought severity in Xinjiang, China. Four drought indices, including the self-calibrating Palmer drought severity index (sc-PDSI), Erinç’s index (I _m), Sahin’s index (I _sh), and UNEP aridity index (AI), were used to compare drought severity. The ensemble empirical mode decomposition and the modified Mann-Kendall trend test were applied to analyze the nonlinear components and trends of the climatic variable and drought indices. Four and six climatic scenarios were generated in sc-PDSI, I _m, I _sh, and AI with different combinations of the observed and detrended climatic variables, respectively. In Xinjiang, generally increasing trends in minimal, average, and maximal air temperature (T _min, T _ave, T _max) and precipitation (P) were found, whereas a decreasing trend in wind speed at 2 m height (U ₂) was observed. There were significantly increasing trends in all of the four studied drought indices. Drought relief was more obvious in northern Xinjiang than in southern Xinjiang. The strong influences of increased P on drought relief and the weak influences of increased T _min, T _ave, and T _max on drought aggravation were shown by comparing four drought indices under different climate scenarios. Decreased U ₂ had a weak influence on drought, as shown by the AI in different climate scenarios. The weak influences of T and U ₂ were considered to be masked by the strong influences of P on droughts. Droughts were expected to be more severe if P did not increase, but were likely milder without an increase in air temperature and with a decrease in U ₂.     

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.     

Drivers of the Severity of the Extreme Hot Summer of 2015 in Western China

Western China experienced an extreme hot summer in 2015, breaking a number of temperature records. The summer mean surface air temperature (SAT) anomaly was twice the interannual variability. The hottest daytime temperature (T_Xx) and warmest night-time temperature (T_Nx) were the highest in China since 1964. This extreme hot summer occurred in the context of steadily increasing temperatures in recent decades. We carried out a set of experiments to evaluate the extent to which the changes in sea surface temperature (SST)/sea ice extent (SIE) and anthropogenic forcing drove the severity of the extreme summer of 2015 in western China. Our results indicate that about 65%–72% of the observed changes in the seasonal mean SAT and the daily maximum (T_max) and daily minimum (T_min) temperatures over western China resulted from changes in boundary forcings, including the SST/SIE and anthropogenic forcing. For the relative role of individual forcing, the direct impact of changes in anthropogenic forcing explain about 42% of the SAT warming and 60% (40%) of the increase in T_Nx and Tmin (T_Xx and T_max) in the model response. The changes in SST/SIE contributed to the remaining surface warming and the increase in hot extremes, which are mainly the result of changes in the SST over the Pacific Ocean, where a super El Niño event occurred. Our study indicates a prominent role for the direct impact of anthropogenic forcing in the severity of the extreme hot summer in western China in 2015, although the changes in SST/SIE, as well as the internal variability of the atmosphere, also made a contribution.     

Climate Change driven evolution of hazards to Europe’s transport infrastructure throughout the twenty-first century

Air circulation due to the urban heat island (UHI) effect can influence the dispersion of air pollutants in a metropolis. This study focusses on the influence of the UHI effect on particulate matter (PM; including PM_2.5 and PM_2.5–10) between May and September 2010–2012 in the Taipei basin. Meteorological and PM data were obtained from the sites, owned by the governmental authorities. The analysis was carried out using t test, relative indices (RIs), Pearson product–moment correlation and stepwise regression. The results show that the RI values for PM were the highest at moderate UHI intensity (MUI; 2 °C ≤ UHI < 4 °C) rather than at strong UHI intensity (SUI; 4 °C ≤ UHI) during the peak time for anthropogenic emissions (20:00 LST). Neither the accumulation of PM nor the surface convergence occurred in the hot centre, as shown by the case study. At MUI, more than 89 % of the synoptic weather patterns showed that the weather was clear and hot or that the atmosphere was stable. The variation in PM was associated with horizontal and vertical air dispersion. Poor horizontal air dispersion, with subsidence, caused an increase in PM at MUI. However, the updraft motion diluted the PM at SUI. The stepwise regression models show that the cloud index and surface air pressure determined the variation in PM_2.5–10, while cloud index, wind speed and mixing height influenced the variation in PM_2.5. In conclusion, a direct relationship between UHI effect and PM was not obvious.     

Effect of urbanization on the urban meteorology and air pollution in Hangzhou

Urbanization has a substantial effect on urban meteorology. It can alter the atmospheric diffusion capability in urban areas and therefore affect pollutant concentrations. To study the effects of Hangzhou’s urban development in most recent decade on its urban meteorological characteristics and pollutant diffusion, 90 weather cases were simulated, covering 9 weather types, with the Nanjing University City Air Quality Prediction System and high-resolution surface-type data and urban construction data for 2000 and 2010. The results show that the most recent decade of urban development in Hangzhou substantially affected its urban meteorology. Specifically, the average urban wind speed decreased by 1.1 m s ^?1; the average intensity of the heat island increased by 0.5°C; and the average urban relative humidity decreased by 9.7%. Based on one case for each of the nine weather types, the impact of urbanization on air pollution diffusion was investigated, revealing that the changes in the meteorological environment decreased the urban atmosphere’s diffusion capability, and therefore increased urban pollutant concentrations. For instance, the urban nitrogen oxides concentration increased by 2.1 μg m ^?3 on average; the fine particulate matter (diameter of 2.5 μm or less; PM_2.5) pollution concentration increased by 2.3 μg m ^?3 on average; in highly urbanized areas, the PM_2.5 concentration increased by 30 μg m ^?3 and average visibility decreased by 0.2 km, with a maximum decrease of 1 km; the average number of daily hours of haze increased by 0.46 h; and the haze height lifted by 100–300 m. The “self-cleaning time” of pollutants increased by an average of 1.5 h.     

Simulation of leaf blast infection in tropical rice agro-ecology under climate change scenario

Assessing disease risk has become an important component in the development of climate change adaptation strategies. Here, the infection ability of leaf blast (Magnaporthe oryzae) was modeled based on the epidemiological parameters of minimum (T _min), optimum (T _opt), and maximum (T _max) temperatures for sporulation and lesion development. An infection ability response curve was used to assess the impact of rising temperature on the disease. The simulated spatial pattern of the infection ability index (IAI) corresponded with observed leaf blast occurrence in Indo-Gangetic plains (IGP). The IAI for leaf blast is projected to increase during the winter season (December–March) in 2020 (2010–2039) and 2050 (2040–2069) climate scenarios due to temperature rise, particularly in lower latitudes. However, during monsoon season (July–October), the IAI is projected to remain unchanged or even reduce across the IGP. The results show that the response curve may be successfully used to assess the impact of climate change on leaf blast in rice. The model could be further extended with a crop model to assess yield loss.     

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.”     

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.     

Seasonal and annual trends of carbonaceous species of PM<Subscript>10</Subscript> over a megacity Delhi,India during 2010–2017

PM₁₀ samples were collected to characterize the seasonal and annual trends of carbonaceous content in PM₁₀ at an urban site of megacity Delhi, India from January 2010 to December 2017. Organic carbon (OC) and elemental carbon (EC) concentrations were quantified by thermal-optical transmission (TOT) method of PM₁₀ samples collected at Delhi. The average concentrations of PM₁₀, OC, EC and TCA (total carbonaceous aerosol) were 222?±?87 (range: 48.2–583.8 μg m^?3), 25.6?±?14.0 (range: 4.2–82.5 μg m^?3), 8.7?±?5.8 (range: 0.8–35.6 μg m^?3) and 54.7?±?30.6 μg m^?3 (range: 8.4–175.2 μg m^?3), respectively during entire sampling period. The average secondary organic carbon (SOC) concentration ranged from 2.5–9.1 μg m^?3 in PM₁₀, accounting from 14 to 28% of total OC mass concentration of PM₁₀. Significant seasonal variations were recorded in concentrations of PM₁₀, OC, EC and TCA with maxima during winter and minima during monsoon seasons. In the present study, the positive linear trend between OC and EC were recorded during winter (R²?=?0.53), summer (R²?=?0.59) and monsoon (R²?=?0.78) seasons. This behaviour suggests the contribution of similar sources and common atmospheric processes in both the fractions. OC/EC weight ratio suggested that vehicular emissions, fossil fuel combustion and biomass burning could be the major sources of carbonaceous aerosols of PM₁₀ at the megacity Delhi, India. Trajectory analysis indicates that the air mass approches to the sampling site is mainly from Indo Gangetic plain (IGP) region (Uttar Pradesh, Haryana and Punjab etc.), Thar desert, Afghanistan, Pakistan and surrounding areas.     

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.     

Satellite bulk tropospheric temperatures as a metric for climate sensitivity

We identify and remove the main natural perturbations (e.g. volcanic activity, ENSOs) from the global mean lower tropospheric temperatures (T _LT ) over January 1979 - June 2017 to estimate the underlying, potentially human-forced trend. The unaltered value is +0.155 K dec^?1 while the adjusted trend is +0.096 K dec^?1, related primarily to the removal of volcanic cooling in the early part of the record. This is essentially the same value we determined in 1994 (+0.09 K dec^?1, Christy and McNider, 1994) using only 15 years of data. If the warming rate of +0.096 K dec^?1 represents the net T _LT response to increasing greenhouse radiative forcings, this implies that the T _LT tropospheric transient climate response (ΔT _LT at the time CO₂ doubles) is +1.10 ± 0.26 K which is about half of the average of the IPCC AR5 climate models of 2.31 ± 0.20 K. Assuming that the net remaining unknown internal and external natural forcing over this period is near zero, the mismatch since 1979 between observations and CMIP-5 model values suggests that excessive sensitivity to enhanced radiative forcing in the models can be appreciable. The tropical region is mainly responsible for this discrepancy suggesting processes that are the likely sources of the extra sensitivity are (a) the parameterized hydrology of the deep atmosphere, (b) the parameterized heat-partitioning at the oceanatmosphere interface and/or (c) unknown natural variations.     

High-stability algorithm for the three-pattern decomposition of global atmospheric circulation

The study investigates the reliable computation time (RCT, termed as T _c) by applying a double-precision computation of a variable parameters logistic map (VPLM). Firstly, by using the proposed method, we obtain the reliable solutions for the logistic map. Secondly, we construct 10,000 samples of reliable experiments from a time-dependent non-stationary parameters VPLM and then calculate the mean T _c. The results indicate that, for each different initial value, the T _cs of the VPLM are generally different. However, the mean T _c trends to a constant value when the sample number is large enough. The maximum, minimum, and probable distribution functions of T _c are also obtained, which can help us to identify the robustness of applying a nonlinear time series theory to forecasting by using the VPLM output. In addition, the T _c of the fixed parameter experiments of the logistic map is obtained, and the results suggest that this T _c matches the theoretical formula-predicted value.     

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.     

Seasonal and interannual variability of oceanographic conditions in the southwestern part of the Sea of Okhotsk

The seasonal variability of oceanographic conditions in the southern part of the Sea of Okhotsk is described based on long-term mean temperature T and salinity S from observations along a standard oceanographic section Cape Aniva-Cape Dokuchaev (May–November). It is shown that the Soya Current is relatively weak in spring, with low temperature and salinity gradients along the section. The Sea of Okhotsk low-salinity water mass is observed in the upper layer. It was formed as a result of melting of a large amount of ice brought here with the East Sakhalin Current from the northwestern part of the Sea of Okhotsk. A cold intermediate layer (CIL) at depths of 50–150 m extends along the entire section. The cold intermediate layer core with a temperature at the edge of the Sakhalin shelf of about ?1.3°C is retained during a period of maximum warming in August; however, in October–November the intensified flow of the East Sakhalin Current (up to 50 cm/s) results in a situation when relatively warm low-salinity waters, connected with the Amur River runoff, dissipate CIL. The results of 12 surveys conducted by the Sakhalin Research Institute for Fisheries and Oceanography in 1998–2004 show significant deviations of T and S [10] in different years from the calculated values. Generally, maximum anomalies (ΔT > 4°C and ΔS > 0.55‰) are observed in the surface layer. Their values and statistical significance decrease with depth. However, the situation is opposite in some cases. The maximum deviation from normal was observed in June 1999, when warm and salt waters were located much further seaward from the Kunashir shelf, which is most likely connected with the Soya Current meandering.     

Characteristics of Turbulent Transfer during Episodes of Heavy Haze Pollution in Beijing in Winter 2016/17

We analyzed the structure and evolution of turbulent transfer and the wind profile in the atmospheric boundary layer in relation to aerosol concentrations during an episode of heavy haze pollution from 6 December 2016 to 9 January 2017. The turbulence data were recorded at Peking University’s atmospheric science and environment observation station. The results showed a negative correlation between the wind speed and the PM_2.5 concentration. The turbulence kinetic energy was large and showed obvious diurnal variations during unpolluted (clean) weather, but was small during episodes of heavy haze pollution. Under both clean and heavy haze conditions, the relation between the non-dimensional wind components and the stability parameter z/L followed a 1/3 power law, but the normalized standard deviations of the wind speed were smaller during heavy pollution events than during clean periods under near-neutral conditions. Under unstable conditions, the normalized standard deviation of the potential temperature σ _θ /|θ_*| was related to z/L, roughly following a –1/3 power law, and the ratio during pollution days was greater than that during clean days. The three-dimensional turbulence energy spectra satisfied a –2/3 power exponent rate in the high-frequency band. In the low-frequency band, the wind velocity spectrum curve was related to the stability parameters under clear conditions, but was not related to atmospheric stratification under polluted conditions. In the dissipation stage of the heavy pollution episode, the horizontal wind speed first started to increase at high altitudes and then gradually decreased at lower altitudes. The strong upward motion during this stage was an important dynamic factor in the dissipation of the heavy haze.     

Velocity Adjustment and Passive Scalar Diffusion in and Above an Urban Canopy in Response to Various Approach Flows

We used wind-tunnel experiments to investigate velocity-field adjustment and scalar diffusion behaviour in and above urban canopies located downwind of various roughness elements. Staggered arrays of rectangular blocks of various heights H and plan area ratios λ_p were used to model the urban canopies. The velocity field in the roughness sublayer (height \({z \lesssim 2H}\)) reached equilibrium at distances proportional to \({\sqrt{L_{\rm c}H}}\) where L _c is the canopy-drag length scale determined as a function of λ_p and the block side length L. A distance of about \({20\sqrt{L_{\rm c}H}}\) was required for adjustment at z = H/2 (in the canopy), and a distance of about \({10\sqrt{L_{\rm c}H}}\) was required at z = 2H (near the top of the roughness sublayer). Diffusion experiments from a ground emission source revealed that differences in upwind roughness conditions had negligible effects on the plume growth near the source (up to a few multiples of L from the source) if the source was located at a fetch F larger than about \({10\sqrt{L_{\rm c}H}}\) from the upwind edge of the canopy. However, at locations farther downwind (more than several multiples of L from the source), upwind conditions had considerable effects on the plume growth. For a representative urban canopy, it was shown that a much larger fetch than required for velocity-field adjustment in the roughness sublayer was necessary to eliminate the effects of upwind conditions on plume widths at 24L downwind from the source.     

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.     

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.     

Long-term variation of rainfall erosivity in Calabria (Southern Italy)

The changes in rainfall erosivity have been investigated using the rainfall erosivity factor (R) proposed for USLE by Wischmeier and Smith (R _W-S ) and some simplified indexes (the Fournier index modified by Arnoldus, F, a regional index spatial independent, R _Fr , and a regional index spatial dependent, R _Fs ) estimated by indirect approaches. The analysis has been carried out over 48 rainfall stations located in Calabria (Southern Italy) using data collected in the period 1936–2012 and divided in three sub-periods. The series of the erosivity indexes and of some precipitation variables have been analyzed for evidence of trends using standard methods. The simplified indexes suggested a general underestimation of the rainfall erosivity with respect to R _W-S . The mean underestimation ranged between 23 and 54 % for R _Fr and from 10 to 15 % for R _Fs . Both the sign and the magnitude of the trends were different for the different stations depending on the variable and sub-period considered. In general, the erosivity increased during the period 1936–1955 (1^st sub-period) and during the more recent sub-period (1992–2012, 3^rd sub-period), whereas it decreased during 1958–1977 (2^nd sub-period). The evidence of trends was generally higher for R _W-S than for R _Fr and R _Fs . Focusing on the most recent sub-period (3^rd sub-period), all the variables analyzed showed mainly increasing trends but with different magnitude. More particularly, R _W-S showed a mean increment of 29 %; F, R _Fr and R _Fs increased by 11, 15 and 18 %, respectively; the maximum intensity of 0.5-h precipitation increased by 5 %; and the annual precipitation increased by 22 %. Consequently, it remains difficult to define which precipitation variable plays the dominant role in the temporal variation of rainfall erosivity in the region. However, the overall results suggest that the indexes estimated by indirect procedures (F, R _Fr , and R _Fs ) should be used with caution for climate change analysis, despite they are used for practical purposes considering they are based on easily available information.     

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

We evaluated two methods to estimate evapotranspiration (ETo) from minimal weather records (daily maximum and minimum temperatures) in Mexico: a modified reduced set FAO-Penman-Monteith method (Allen et al. 1998, Rome, Italy) and the Hargreaves and Samani (Appl Eng Agric 1(2): 96–99, 1985) method. In the reduced set method, the FAO-Penman-Monteith equation was applied with vapor pressure and radiation estimated from temperature data using two new models (see first and second articles in this series): mean temperature as the average of maximum and minimum temperature corrected for a constant bias and constant wind speed. The Hargreaves-Samani method combines two empirical relationships: one between diurnal temperature range ΔT and shortwave radiation Rs, and another one between average temperature and the ratio ETo/Rs: both relationships were evaluated and calibrated for Mexico. After performing a sensitivity analysis to evaluate the impact of different approximations on the estimation of Rs and ETo, several model combinations were tested to predict ETo from daily maximum and minimum temperature alone. The quality of fit of these models was evaluated on 786 weather stations covering most of the territory of Mexico. The best method was found to be a combination of the FAO-Penman-Monteith reduced set equation with the new radiation estimation and vapor pressure model. As an alternative, a recalibration of the Hargreaves-Samani equation is proposed.