Seasonal inhomogeneity of soot particles over the central Indo-Gangetic Plains,India: Influence of meteorology

Black carbon (BC) particles play a unique and important role in earth’s climate system. BC was measured (in-situ) in the central part of the Indo-Gangetic Plains (IGP) at Varanasi, which is a highly populated and polluted region due to its topography and extensive emission sources. The annual mean BC mass concentration was 8.92 ± 7.0 µg m ^-3, with 34% of samples exceeding the average value. Seasonally, BC was highest during the post-monsoon and winter periods (approximately 18 µg m ^-3) and lower in the premonsoon/ monsoon seasons (approximately 6 µg m ^-3). The highest frequency (approximately 46%) observed for BC mass was in the interval from 5 to 10 µg m ^-3. However, during the post-monsoon season, the most common values (approximately 23%) were between 20 and 25 µg m ^-3. The nighttime concentrations of BC were approximately twice as much as the daytime values because of lower boundary layer heights at nighttime. The Ångström exponent was significantly positively correlated (0.55) with ground-level BC concentrations, indicating the impact of BC on the columnar aerosol properties. The estimated mean absorption Ångström exponent was 1.02 ± 0.08 µg m ^-3, indicating that the major source of BC was from fossil fuel combustion. Significant negative correlations between BC mass and meteorological parameters indicate a pronounced effect of atmospheric dynamics on the BC mass in this region. The highest mean BC mass concentration (18.1 ± 6.9 µg m ^-3) as a function of wind speed was under calm wind conditions (38% of the time).     

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.     

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.     

Cold air outbreaks along a non-frozen sea channel: effects of wind on snow bands

Wintertime cold air outbreaks along a non-frozen sea channel or a long lake can become destructive if the related bands of heavy snowfall hit onto land. The forcing for such bands is studied with a 2D numerical model set across an east–west sea channel at 60^oN (‘Gulf of Finland’), varying the basic geostrophic wind V _g. Without any V _g opposite coastal land breezes emerge with convergence. This results in a quasi-steady rising motion w _max ~ 7.5 cm/s at 600 m in the middle of the gulf, which can force a snow band. During weak V _g, the rising motion is reduced but least so for winds from 60^o to 80^o (~ENE), when modest alongshore bands could exist near the downstream (Estonian) coast. During V _g of 4–6 m/s from any direction, the land breezes and rising motions are reduced more effectively, so snow bands are not expected during moderate basic flow. In contrast, during a strong V _g of 20–25 m/s from 110^o to 120^o (~ESE) the land breeze perturbations are intense with w _max up to 15–18 cm/s. The induced alongshore bands of heavy snowfall are located in these cases at the sea but quite close to the downstream (Finnish) coast. They can suddenly make a landfall if the basic wind turns clockwise.     

A new approach to identify the sensitivity and importance of physical parameters combination within numerical models using the Lund–Potsdam–Jena (LPJ) model as an example

Daily global solar irradiation (R _s) is one of the main inputs in environmental modeling. Because of the lack of its measuring facilities, high-quality and long-term data are limited. In this research, R _s values were estimated based on measured sunshine duration and cloud cover of our synoptic meteorological stations in central and southern Iran during 2008, 2009, and 2011. Clear sky solar irradiation was estimated from linear regression using extraterrestrial solar irradiation as the independent variable with normalized root mean square error (NRMSE) of 4.69 %. Daily R _s was calibrated using measured sunshine duration and cloud cover data under different sky conditions during 2008 and 2009. The 2011 data were used for model validation. According to the results, in the presence of clouds, the R _s model using sunshine duration data was more accurate when compared with the model using cloud cover data (NRMSE = 11. 69 %). In both models, with increasing sky cloudiness, the accuracy decreased. In the study region, more than 92 % of sunshine durations were clear or partly cloudy, which received close to 95 % of total solar irradiation. Hence, it was possible to estimate solar irradiation with a good accuracy in most days with the measurements of sunshine duration.     

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.     

Indirect determination of broadband turbidity coefficients over Egypt

Long-term data from diffuse and global irradiances were used to calculate direct beam irradiance which was used to determine three atmospheric turbidity coefficients (Linke T _L , Ångström β and Unsworth–Monteith δ _a ) at seven sites in Egypt in the period from 1981 to 2000. Seven study sites (Barrani, Matruh, Arish, Cairo, Asyut, Aswan and Kharga) have been divided into three categories: Mediterranean climate (MC), desert Nile climate (DNC) and urban climate (UC, Cairo). The indirect method (i.e., global irradiance minus diffuse irradiance) used here allows to estimate the turbidity coefficients with an RMSE% ≤20 % (for β, δ _a and T _L ) and ~30 % (for β) if compared with those estimated by direct beam irradiance and sunphotometeric data, respectively. Monthly averages of T _L , β and δ _a show seasonal variations with mainly maxima in spring at all stations, due to Khamsin depressions coming from Sahara. Secondary maxima is observed in summer and autumn at DNC and MC (Barrani and Arish) stations in summer due to dust haze which prevails during that season and at UC (Cairo) in autumn, due to the northern extension of the Sudan monsoon trough, which is accompanied by small-scale depressions with dust particles. The mean annual values of β, δ _a , and T _L (0.216, 0.314, and 4.6, respectively) are larger in Cairo than at MC stations (0.146, 0.216, and 3.8, respectively) and DNC stations (0.153, 0.227, and 3.8, respectively). Both El-Chichon and Mt. Pinatubo eruptions were examined for all records data at MC, UC and DNC stations. The overburden caused by Mt. Pinatubo’s eruption was larger than El-Chichon’s eruption and overburden for β, and T _L at DNC stations (0.06, and 0.58 units, respectively) was more pronounced than that at MC (0.02, and 0.26, respectively) and UC (0.05 and 0.52 units, respectively) stations. The annual variations in wind speed and turbidity parameters show high values for both low and high wind speed at all stations. The wind directions have a clear effect on atmospheric turbidity, and consequently, largest turbidities occur when the wind carries aerosols from the main particle sources, such as industrial particle sources around Cairo or to some extent from the Sahara surrounding all study stations.     

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.     

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.     

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

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.     

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.     

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.     

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.     

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.     

Estimation of ground heat flux from soil temperature over a bare soil

Ground soil heat flux, G ₀, is a difficult-to-measure but important component of the surface energy budget. Over the past years, many methods were proposed to estimate G ₀; however, the application of these methods was seldom validated and assessed under different weather conditions. In this study, three popular models (force-restore, conduction-convection, and harmonic) and one widely used method (plate calorimetric), which had well performance in publications, were investigated using field data to estimate daily G ₀ on clear, cloudy, and rainy days, while the gradient calorimetric method was regarded as the reference for assessing the accuracy. The results showed that harmonic model was well reproducing the G ₀ curve for clear days, but it yielded large errors on cloudy and rainy days. The force-restore model worked well only under rainfall condition, but it was poor to estimate G ₀ under rain-free conditions. On the contrary, the conduction-convection model was acceptable to determine G ₀ under rain-free conditions, but it generated large errors on rainfall days. More importantly, the plate calorimetric method was the best to estimate G ₀ under different weather conditions compared with the three models, but the performance of this method is affected by the placement depth of the heat flux plate. As a result, the heat flux plate was recommended to be buried as close as possible to the surface under clear condition. But under cloudy and rainy conditions, the plate placed at depth of around 0.075 m yielded G ₀ well. Overall, the findings of this paper provide guidelines to acquire more accurate estimation of G ₀ under different weather conditions, which could improve the surface energy balance in field.