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.     

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.     

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

Formulation of stability-dependent empirical relations for turbulent intensities from surface layer turbulence measurements for dispersion parameterization in a lagrangian particle dispersion model

Season- and stability-dependent turbulence intensity (σ _u /u _*, σ _v /u _*, σ _w /u _*) relationships are derived from experimental turbulence measurements following surface layer scaling and local stability at the tropical coastal site Kalpakkam, India for atmospheric dispersion parameterization. Turbulence wind components (u′, v′, w′) measured with fast response UltraSonic Anemometers during an intense observation campaign for wind field modeling called Round Robin Exercise are used to formulate the flux–profile relationships using surface layer similarity theory and Fast Fourier Transform technique. The new relationships (modified Hanna scheme) are incorporated in a Lagrangian Particle Dispersion model FLEXPART-WRF and tested by conducting simulations for a field tracer dispersion experiment at Kalpakkam. Plume dispersion analysis of a ground level hypothetical release indicated that the new turbulent intensity formulations provide slightly higher diffusivity across the plume relative to the original Hanna scheme. The new formulations for σ _u , σ _v , σ _w are found to give better agreement with observed turbulent intensities during both stable and unstable conditions under various seasonal meteorological conditions. The simulated concentrations using the two methods are compared with those obtained from a classical Gaussian model and the observed SF₆ concentration. It has been found that the new relationships provide comparatively higher diffusion across the plume relative to the model default Hanna scheme and provide downwind concentration results in better agreement with observations.     

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.     

Wind-Speed Undulations Over Swell: Field Experiment and Interpretation

Results of field measurements of the swell-induced undulation of the wind speed taken from a Black Sea platform are presented. The wind speed and its fluctuations were measured at several heights between 1.3 and 21 m above the mean sea level under various wind and swell conditions. Parameters of the swell-induced undulations were derived from cross spectra of the wind-speed fluctuations and the sea-surface displacement. As found, the phase and the amplitude of the wind speed undulation in the layer from k _p z = 0.1 to k _p z = 3 (k _p is the swell wavenumber) are in good agreement with the theory of inviscid shear flow over a wavy surface. The main feature of the vertical profile of the swell-induced undulation is the exponential attenuation of its amplitude with height typical for the potential flow over the fast running waves. At the lowest levels the potential undulations are significantly distorted by the wind-speed variations caused by the vertical displacements of the shear airflow relative to a fixed sensor. No direct impact of swell on the mean properties of the turbulent boundary layer at k _p z > 0.1 is revealed. In particular, the mean wind-speed profile and spectra of the horizontal velocity in the inertial subrange obey Monin-Obukhov similarity theory.     

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.     

Field Experiments and Numerical Modeling of Wind Speed and Surface Waves in Medium-size Inland Reservoirs

An attempt is made to apply the modern methods of surface wave simulation developed for oceanic conditions to the modeling of waves in medium-size inland reservoirs (10–100 km). The results of field measurements of wind speed and waves are described, and on their basis the parameterization C _D (U ₁₀) is proposed. WAVEWATCH III spectral wave model was adapted to the conditions of a medium-size inl and reservoir. The simulated data are compared with the field data. The use of the new parameterization C _D (U ₁₀) allowed reducing the values of the wind wave growth rate that improved consistency in data from the field experiment and numerical modeling concerning the height of significant waves. Further steps towards improving the quality of prediction of the adapted WAVEWATCH III model are discussed.     

A scattering-based over-land rainfall retrieval algorithm for South Korea using GCOM-W1/AMSR-2 data

Heavy summer rainfall is a primary natural disaster affecting lives and properties in the Korean Peninsula. This study presents a satellite-based rainfall rate retrieval algorithm for the South Korea combining polarization-corrected temperature (PCT) and scattering index (SI) data from the 36.5 and 89.0 GHz channels of the Advanced microwave Scanning Radiometer 2 (AMSR-2) onboard the Global Change Observation Mission (GCOM)-W1 satellite. The coefficients for the algorithm were obtained from spatial and temporal collocation data from the AMSR-2 and groundbased automatic weather station rain gauges from 1 July - 30 August during the years, 2012-2015. There were time delays of about 25 minutes between the AMSR-2 observations and the ground raingauge measurements. A new linearly-combined rainfall retrieval algorithm focused on heavy rain for the PCT and SI was validated using ground-based rainfall observations for the South Korea from 1 July - 30 August, 2016. The validation presented PCT and SI methods showed slightly improved results for rainfall > 5 mm h^-1 compared to the current ASMR-2 level 2 data. The best bias and root mean square error (RMSE) for the PCT method at AMSR-2 36.5 GHz were 2.09 mm h^-1 and 7.29 mm h^-1, respectively, while the current official AMSR-2 rainfall rates show a larger bias and RMSE (4.80 mm h^-1 and 9.35 mm h^-1, respectively). This study provides a scatteringbased over-land rainfall retrieval algorithm for South Korea affected by stationary front rain and typhoons with the advantages of the previous PCT and SI methods to be applied to a variety of spaceborne passive microwave radiometers.     

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.     

Correlation between temperature and salinity variations as a characteristic of the North Atlantic waters

Based on the data of deep-ocean ship observations of temperature T and salinity S, analysis is carried out of the fields of pair correlation coefficients between T and S at different depths as an additional characteristic of water masses in the layer 0–1000 m in the North Atlantic. As a result of analysis, surface, subsurface, and the upper part of North Atlantic intermediate waters are classified according to a degree of correlation between temperature and salinity. The emphasis was given to regions with low correlations, because they indicate the prevalence of the interaction processes that differ most from the typical mixing of two water masses with entirely different characteristics.     

Aerodynamic Roughness Length of Fresh Snow

This study presents the results from a series of wind-tunnel experiments designed to investigate the aerodynamic roughness length z ₀ of fresh snow under no-drift conditions. A two-component hot-film anemometer was employed to obtain vertical profiles of velocity statistics in a zero pressure gradient turbulent boundary layer for flow over naturally deposited snow surfaces. The roughness of these snow surfaces was measured by means of digital photography to capture characteristic length scales that can be related to z ₀. Our results show that, under aerodynamically rough conditions, the mean value of the roughness length for fresh snow is \({\langle{z}_{0}\rangle= 0.24}\) mm with a standard deviation σ(z ₀) = 0.05 mm. In this study, we show that variations in z ₀ are associated with variations in the roughness geometry. The roughness measurements suggest that the estimated values of z ₀ are consistent with the presence of irregular roughness structures that develop during snowfalls that mimic ballistic deposition processes.     

Statistics of the <Emphasis Type="Italic">Z</Emphasis>–<Emphasis Type="Italic">R</Emphasis> Relationship for Strong Convective Weather over the Yangtze–Huaihe River Basin and Its Application to Radar Reflectivity Data Assimilation for a Heavy Rain Event

The relationship between the radar reflectivity factor (Z) and the rainfall rate (R) is recalculated based on radar observations from 10 Doppler radars and hourly rainfall measurements at 6529 automatic weather stations over the Yangtze–Huaihe River basin. The data were collected by the National 973 Project from June to July 2013 for severe convective weather events. The Z–R relationship is combined with an empirical q_r–R relationship to obtain a new Z–q_r relationship, which is then used to correct the observational operator for radar reflectivity in the three-dimensional variational (3DVar) data assimilation system of the Weather Research and Forecasting (WRF) model to improve the analysis and prediction of severe convective weather over the Yangtze–Huaihe River basin. The performance of the corrected reflectivity operator used in the WRF 3DVar data assimilation system is tested with a heavy rain event that occurred over Jiangsu and Anhui provinces and the surrounding regions on 23 June 2013. It is noted that the observations for this event are not included in the calculation of the Z–R relationship. Three experiments are conducted with the WRF model and its 3DVar system, including a control run without the assimilation of reflectivity data and two assimilation experiments with the original and corrected reflectivity operators. The experimental results show that the assimilation of radar reflectivity data has a positive impact on the rainfall forecast within a few hours with either the original or corrected reflectivity operators, but the corrected reflectivity operator achieves a better performance on the rainfall forecast than the original operator. The corrected reflectivity operator extends the effective time of radar data assimilation for the prediction of strong reflectivity. The physical variables analyzed with the corrected reflectivity operator present more reasonable mesoscale structures than those obtained with the original reflectivity operator. This suggests that the new statistical Z–R relationship is more suitable for predicting severe convective weather over the Yangtze–Huaihe River basin than the Z–R relationships currently in use.     

Semi-idealized modeling of lightning initiation related to vertical air motion and cloud microphysics

A three-dimensional charge–discharge numerical model is used, in a semi-idealized mode, to simulate a thunder-storm cell. Characteristics of the graupel microphysics and vertical air motion associated with the lightning initiation are revealed, which could be useful in retrieving charge strength during lightning when no charge–discharge model is available. The results show that the vertical air motion at the lightning initiation sites (W_ini) has a cubic polynomial correlation with the maximum updraft of the storm cell (W_cell-max), with the adjusted regression coefficient R² of approximately 0.97. Meanwhile, the graupel mixing ratio at the lightning initiation sites (q_g-ini) has a linear correlation with the maximum graupel mixing ratio of the storm cell (q_g-cell-max) and the initiation height (z_ini), with the coefficients being 0.86 and 0.85, respectively. These linear correlations are more significant during the middle and late stages of lightning activity. A zero-charge zone, namely, the area with very low net charge density between the main positive and negative charge layers, appears above the area of q_g-cell-max and below the upper edge of the graupel region, and is found to be an important area for lightning initiation. Inside the zero-charge zone, large electric intensity forms, and the ratio of q_ice (ice crystal mixing ratio) to q_g (graupel mixing ratio) illustrates an exponential relationship to q_g-ini. These relationships provide valuable clues to more accurately locating the high-risk area of lightning initiation in thunderstorms when only dual-polarization radar data or outputs from numerical models without charging/discharging schemes are available. The results can also help understand the environmental conditions at lightning initiation sites.     

Raindrop Size Distribution Parameters Retrieved from Guangzhou S-band Polarimetric Radar Observations

According to the statistical shape–slope (μ–Λ) relationship observed for the first time by several 2D-Video-Distrometers (2DVD) in southern China, a constrained gamma (C-G) model was proposed for the retrieval of rain drop size distributions (DSDs) from Guangzhou S-band polarimetric radar observations. Two typical precipitation processes were selected to verify the accuracy of the retrieval scheme. The μ–Λ relationship: Λ = 0.0241μ² + 0.867μ + 2.453 was obtained based on the 2DVD observation results from at Huizhou Longmen station, which is a very representative location in the area. Relying on the Guangzhou polarimetric radar measurements of radar reflectivity (Z_HH) and differential reflectivity (Z_DR), the gamma (Γ) size distribution parameters (N₀, μ, and Λ) can be retrieved by the C-G model retrieval scheme. The results show that the Guangzhou polarimetric radar retrievals of DSDs were close to the 2DVD observations at Guangzhou Maofengshan station. The rain rate, mass mean diameter, and normalized intercept parameter of radar retrievals were in good agreement with the 2DVD observations, and the relative errors were less than 10%. The overall accuracy of the retrieval scheme was high. The retrieval scheme has established the relationship between the polarimetric radar measurements and gamma size distribution parameters. It will be helpful to in-depth research and application of the dual-polarization radar data in microphysical precipitation processes analysis, as well as convection-resolved numerical model data assimilation and prediction effect evaluation.     

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.     

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.     

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

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.     

On the extension of the wind profile over homogeneous terrain beyond the surface boundary layer

Analysis of profiles of meteorological measurements from a 160 m high mast at the National Test Site for wind turbines at Høvsøre (Denmark) and at a 250 m high TV tower at Hamburg (Germany) shows that the wind profile based on surface-layer theory and Monin-Obukhov scaling is valid up to a height of 50–80 m. At higher levels deviations from the measurements progressively occur. For applied use an extension to the wind profile in the surface layer is formulated for the entire boundary layer, with emphasis on the lowest 200–300 m and considering only wind speeds above 3 m s^?1 at 10 m height. The friction velocity is taken to decrease linearly through the boundary layer. The wind profile length scale is composed of three component length scales. In the surface layer the first length scale is taken to increase linearly with height with a stability correction following Monin-Obukhov similarity. Above the surface layer the second length scale (L _MBL ) becomes independent of height but not of stability, and at the top of the boundary layer the third length scale is assumed to be negligible. A simple model for the combined length scale that controls the wind profile and its stability dependence is formulated by inverse summation. Based on these assumptions the wind profile for the entire boundary layer is derived. A parameterization of L _MBL is formulated using the geostrophic drag law, which relates friction velocity and geostrophic wind. The empirical parameterization of the resistance law functions A and B in the geostrophic drag law is uncertain, making it impractical. Therefore an expression for the length scale, L _MBL , for applied use is suggested, based on measurements from the two sites.