Climatic change and Chinese population growth dynamics over the last millennium

The climate–population relationship has long been conceived. Although the topic has been repeatedly investigated, most of the related works are Eurocentric or qualitative. Consequently, the relationship between climate and population remains ambiguous. In this study, fine-grained temperature reconstructions and historical population data sets have been employed to statistically test a hypothesized relationship between temperature change and population growth (i.e., cooling associated with below average population growth) in China over the past millennium. The important results were: (1) Long-term temperature change significantly determined the population growth dynamics of China. However, spatial variation existed, whilst population growth in Central China was shown to be responsive to both long- and short-term temperature changes; in marginal areas, population growth was only sensitive to short-term temperature fluctuations. (2) Temporally, the temperature–population relationship was obscured in some periods, which was attributable to the factors of drought and social buffers. In summary, a temperature–population relationship was mediated by geographic factors, the aridity threshold, and social factors. Given the upcoming threat posed by climate change to human societies, this study seeks to improve our knowledge and understanding of the climate–society relationship.     

A Parameterization of Dry Thermals and Shallow Cumuli for Mesoscale Numerical Weather Prediction

For numerical weather prediction models and models resolving deep convection, shallow convective ascents are subgrid processes that are not parameterized by classical local turbulent schemes. The mass flux formulation of convective mixing is now largely accepted as an efficient approach for parameterizing the contribution of larger plumes in convective dry and cloudy boundary layers. We propose a new formulation of the EDMF scheme (for Eddy Diffusivity\Mass Flux) based on a single updraft that improves the representation of dry thermals and shallow convective clouds and conserves a correct representation of stratocumulus in mesoscale models. The definition of entrainment and detrainment in the dry part of the updraft is original, and is specified as proportional to the ratio of buoyancy to vertical velocity. In the cloudy part of the updraft, the classical buoyancy sorting approach is chosen. The main closure of the scheme is based on the mass flux near the surface, which is proportional to the sub-cloud layer convective velocity scale w _*. The link with the prognostic grid-scale cloud content and cloud cover and the projection on the non- conservative variables is processed by the cloud scheme. The validation of this new formulation using large-eddy simulations focused on showing the robustness of the scheme to represent three different boundary layer regimes. For dry convective cases, this parameterization enables a correct representation of the countergradient zone where the mass flux part represents the top entrainment (IHOP case). It can also handle the diurnal cycle of boundary-layer cumulus clouds (EUROCS\ARM) and conserve a realistic evolution of stratocumulus (EUROCS\FIRE).     

Experimental Investigation of Turbulent Momentum Transfer in a Neutral Boundary Layer over a Rough Surface

The turbulent characteristics of the neutral boundary layer developing over rough surfaces are not well predicted with operational weather-forecasting models. The problem is attributed to inadequate mixing-length models, to the anisotropy of the flow and to a lack of controlled experimental data against which to validate numerical studies. Therefore, in order to address directly the modelling difficulties for the development of a neutral boundary layer over rough surfaces, and to investigate the turbulent momentum transfer of such a layer, a set of hydraulic flume experiments were carried out. In the experiments, the mean and turbulent quantities were measured by a particle image velocimetry (PIV) technique. The measured velocity variances and fluxes \({(\overline{{u_{i}^{\prime}}{u_{j}^{\prime}}})}\) in longitudinal vertical planes allowed the vertical and longitudinal gradients (?/?z and ?/?x) of the mean and turbulent quantities (fluxes, variances and third-order moments) to be evaluated and the terms of the evolution equations for ?e/?t, \({\partial \overline{u^{\prime 2}}/\partial t}\), \({\partial \overline{w^{\prime 2}}/\partial t}\) and \({\partial \overline{{u^{\prime}}{w^{\prime}}}/\partial t}\) to be quantified, where e is the turbulent kinetic energy. The results show that the pressure-correlation terms allow the turbulent energy to be transferred equitably from \({\overline{{u^{\prime}}^{2}}}\) to \({\overline{{w^{\prime}}^{2}}}\). It appears that the repartition between the constitutive terms of the budget of e, \({\overline{{u^{\prime}}^{2}}}\), \({\overline{{w^{\prime}}^{2}}}\) and \({\overline{{u^{\prime}}{w^{\prime}}}}\) is not significantly affected by the development of the rough neutral boundary layer. For the whole evolution, the transfers of energy are governed by the same terms that are also very similar to the smooth-wall case. The PIV measurements also allowed the spatial integral scales to be computed directly and to be compared with the dissipative and mixing length scales, which were also computed from the data.     

Scaling of Wall-Normal Turbulence Intensity and Vertical Eddy Structures in the Atmospheric Surface Layer

Adequate high-quality data on three-dimensional velocities in the atmospheric surface layer (height \(\delta \)) were acquired in the field at the Qingtu Lake Observation Array. The measurement range occupies nearly the entire logarithmic layer from approximately \(0.006\delta \)–\(0.2\delta \). The turbulence intensity and eddy structures of the velocity fluctuations in the logarithmic region were primarily analyzed, and their variations in the z (wall-normal) direction were revealed. The primary finding was that the turbulent intensity of wall-normal velocity fluctuations exhibits a sharp upswing in the logarithmic region, which differs from classic scaling law and laboratory results. The upswing of the wall-normal turbulence intensity in the logarithmic region is deemed to be linear based on an ensemble of 20 sets of data. In addition, the wall-normal extent of the correlated structures and wall-normal spectra were compared to low Reynolds number results in the laboratory.     

Some Insights into Topographic,Elastic and Self-gravitation Interaction in Modelling Ground Deformation and Gravity Changes in Active Volcanic Areas

Surface displacements and gravity changes due to volcanic sources are influenced by medium properties. We investigate topographic, elastic and self-gravitation interaction in order to outline the major factors that are significant in data modelling. While elastic-gravitational models can provide a suitable approximation to problems of volcanic loading in areas where topographic relief is negligible, for prominent volcanoes the rough topography could affect deformation and gravity changes to a greater extent than self-gravitation. This fact requires the selection, depending on local relief, of a suitable model for use in the interpretation of surface precursors of volcanic activity. We use the three-dimensional Indirect Boundary Element Method to examine the effects of topography on deformation and gravity changes in models of magma chamber inflation/deflation. Topography has a significant effect on predicted surface deformation and gravity changes. Both the magnitude and pattern of the geodetic signals are significantly different compared to half-space solutions. Thus, failure to account for topographic effects in areas of prominent relief can bias the estimate of volcanic source parameters, since the magnitude and pattern of deformation and gravity changes depend on such effects.     

UV attenuation in the cloudy atmosphere

Ultraviolet (UV) energy absorption plays a very important role in the Earth–atmosphere system. Based on observational data for Beijing, we suggest that some atmospheric constituents utilize or transfer UV energy in chemical and photochemical (C&P) reactions, in addition to those which absorb UV energy directly. These constituents are primarily volatile organic compounds (VOCs) emitted from both vegetative and anthropogenic sources. The total UV energy loss in the cloudy atmosphere for Beijing in 1990 was 78.9 Wm⁻². This attenuation was caused by ozone (48.3 Wm⁻²), other compounds in the atmosphere (26.6 Wm⁻²) and a scattering factor (4.0 Wm⁻²). Our results for a cloudy atmosphere in the Beijing area show that the absorption due to these other compounds occurs largely through the mediation of water vapor. This fraction of energy loss has not been fully accounted for in previous models. Observations and previous models results suggest that 1) a cloudy atmosphere absorbs 25∼30 Wm⁻² more solar shortwave radiation than models predict; and 2) aerosols can significantly decrease the downward mean UV-visible radiation and the absorbed solar radiation at the surface by up to 28 and 23 Wm⁻², respectively. Thus, quantitative study of UV and visible absorption by atmospheric constituents involved in homogeneous and heterogeneous C&P reactions is important for atmospheric models.     

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

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

Characteristics of carbonyl compounds in ambient air of Shanghai,China

The levels of carbonyl compounds in Shanghai ambient air were measured in five periods from January 2007 to October 2007 (covering winter, high-air-pollution days, spring, summer and autumn). A total of 114 samples were collected and eighteen carbonyls were identified. Formaldehyde, acetaldehyde and acetone were the most abundant carbonyls and their mean concentrations of 19.40 ± 12.00, 15.92 ± 12.07 and 11.86 ± 7.04 μg m⁻³ respectively, in the daytime for five sampling periods. Formaldehyde and acetaldehyde showed similar diurnal profiles with peak mixing ratios in the morning and early afternoon during the daytime. Their mean concentrations were highest in summer and lowest in winter. Acetone showed reversed seasonal variation. The high molecular weight (HMW, ≥C5) carbonyls also showed obvious diurnal variations with higher concentrations in the daytime in summer and autumn, while they were all not detected in winter. Formaldehyde and acetaldehyde played an important role in removing OH radicals in the atmosphere, but the contribution of acetone was below 1%. The carbonyls levels in high-air-pollution days were reported. More carbonyl species with higher concentrations were found in high-air-pollution days than in spring. These carbonyls were transported with other pollutants from north and northwest in March 27 to April 2, 2007 and then mixed with local sources. Comparing with Beijing and Guangzhou, the concentrations of formaldehyde and acetaldehyde in Shanghai were the highest, which indicated that the air pollution in Shanghai was even worse than expected.     

Adaptation of Pressure Based CFD Solvers for Mesoscale Atmospheric Problems

General purpose Computational Fluid Dynamics (CFD) solvers are frequently used in small-scale urban pollution dispersion simulations without a large extent of ver- tical flow. Vertical flow, however, plays an important role in the formation of local breezes, such as urban heat island induced breezes that have great significance in the ventilation of large cities. The effects of atmospheric stratification, anelasticity and Coriolis force must be taken into account in such simulations. We introduce a general method for adapting pressure based CFD solvers to atmospheric flow simulations in order to take advantage of their high flexibility in geometrical modelling and meshing. Compressibility and thermal stratification effects are taken into account by utilizing a novel system of transformations of the field variables and by adding consequential source terms to the model equations of incompressible flow. Phenomena involving mesoscale to microscale coupled effects can be analyzed without model nesting, applying only local grid refinement of an arbitrary level. Elements of the method are validated against an analytical solution, results of a reference calculation, and a laboratory scale urban heat island circulation experiment. The new approach can be applied with benefits to several areas of application. Inclusion of the moisture transport phenomena and the surface energy balance are important further steps towards the practical application of the method.     

Vertical and Horizontal Transport of Energy and Matter by Coherent Motions in a Tall Spruce Canopy

In the framework of the EGER (ExchanGE processes in mountainous Regions) project, the contribution of coherent structures to vertical and horizontal transports in a tall spruce canopy is investigated. The combination of measurements done in both the vertical and horizontal directions allows us to investigate coherent structures, their temporal scales, their role in flux transport, vertical coupling between the sub-canopy, canopy and air above the canopy, and horizontal coupling in the sub-canopy layer. The temporal scales of coherent structures detected with the horizontally distributed systems in the sub-canopy layer are larger than the temporal scales of coherent structures detected with the vertically distributed systems. The flux contribution of coherent structures to the momentum and sensible heat transport is found to be dominant in the canopy layer. Carbon dioxide and latent heat transport by coherent structures increase with height and reach a maximum at the canopy height. The flux contribution of the ejection decreases with increasing height and becomes dominant above the canopy level. The flux fraction transported during the sweep increases with height and becomes the dominant exchange process at the upper canopy level. The determined exchange regimes indicate consistent decoupling between the sub-canopy, canopy and air above the canopy during evening, nighttime and morning hours, whereas the coupled states and coupled by sweep states between layers are observed mostly during the daytime. Furthermore, the horizontal transport of sensible heat by coherent structures is investigated, and the heterogeneity of the contribution of coherent events to the flux transport is demonstrated. A scheme to determine the horizontal coupling by coherent structures in the sub-canopy layer is proposed, and it is shown that the sub-canopy layer is horizontally coupled mainly in the wind direction. The vertical coupling in most cases is observed together with streamwise horizontal coupling, whereas the cross-stream direction is decoupled.