An Intercomparison of Large-Eddy Simulations of the Stable Boundary Layer

Results are presented from the first intercomparison of large-eddy simulation (LES) models for the stable boundary layer (SBL), as part of the Global Energy and Water Cycle Experiment Atmospheric Boundary Layer Study initiative. A moderately stable case is used, based on Arctic observations. All models produce successful simulations, in as much as they generate resolved turbulence and reflect many of the results from local scaling theory and observations. Simulations performed at 1-m and 2-m resolution show only small changes in the mean profiles compared to coarser resolutions. Also, sensitivity to subgrid models for individual models highlights their importance in SBL simulation at moderate resolution (6.25 m). Stability functions are derived from the LES using typical mixing lengths used in numerical weather prediction (NWP) and climate models. The functions have smaller values than those used in NWP. There is also support for the use of K-profile similarity in parametrizations. Thus, the results provide improved understanding and motivate future developments of the parametrization of the SBL.     

Evaluation of the Weather Research and Forecasting Mesoscale Model for GABLS3: Impact of Boundary-Layer Schemes,Boundary Conditions and Spin-Up

We evaluated the performance of the three-dimensional Weather Research and Forecasting (WRF) mesoscale model, specifically the performance of the planetary boundary-layer (PBL) parametrizations. For this purpose, Cabauw tower observations were used, with the study extending beyond the third GEWEX Atmospheric Boundary-Layer Study (GABLS3) one-dimensional model intercomparison. The WRF model (version 3.4.1) contains 12 different PBL parametrizations, most of which have been only partially evaluated. The GABLS3 case offers a clear opportunity to evaluate model performance, focusing on time series of near-surface weather variables, radiation and surface flux budgets, vertical structure and the nighttime inertial oscillation. The model results revealed substantial differences between the PBL schemes. Generally, non-local schemes tend to produce higher temperatures and higher wind speeds than local schemes, in particular, for nighttime. The WRF model underestimates the 2-m temperature during daytime (about \(2\) K) and substantially underestimates it at night (about \(4\) K), in contrast to the previous studies where modelled 2-m temperature was overestimated. Considering the 10-m wind speed, during the night turbulent kinetic energy based schemes tend to produce lower wind speeds than other schemes. In all simulations the sensible and latent heat fluxes were well reproduced. For the net radiation and the soil heat flux we found good agreement with daytime observations but underestimations at night. Concerning the vertical profiles, the selected non-local PBL schemes underestimate the PBL depth and the low-level jet altitude at night by about 50 m, although with the correct wind speed. The latter contradicts most previous studies and can be attributed to the revised stability function in the Yonsei University PBL scheme. The local, turbulent kinetic energy based PBL schemes estimated the low-level jet altitude and strength more accurately. Compared to the observations, all model simulations show a similar structure for the potential temperature, with a consistent cold bias ( \(\approx \) 2 K) in the upper PBL. In addition to the sensitivity to the PBL schemes, we studied the sensitivity to technical features such as horizontal resolution and domain size. We found a substantial difference in the model performance for a range of 12, 18 and 24 h spin-up times, longer spin-up time decreased the modelled wind speed bias, but it strengthened the negative temperature bias. The sensitivity of the model to the vertical resolution of the input and boundary conditions on the model performance is confirmed, and its influence appeared most significant for the non-local PBL parametrizations.     

Assessing Winter Wheat Responses to Climate Change Scenarios: A Simulation Study in the U.S. Great Plains

Hard red winter wheat (Triticum aestivum L.) is a major crop in the Great Plains region of the U.S. The goal of this assessment effort was to investigate the influence of two contrasting global climate change projections (U.K. Hadley Center for Climate Prediction and Research and Canadian Centre for Climate Modelling and Analysis) on the yield and percent kernel nitrogen content of winter wheat at three locations in Nebraska. These three locations represent sub-humid and semi arid areas and the transition between these areas and are also representative of major portions of the winter wheat growing areas of the central Great Plains. Climate scenarios based on each of the projections for each location were developed using the LARS-WG weather generator along with data from automated weather stations. CERES-Wheat was used to simulate the responses for two contrasting cultivars of wheat using two sowing dates. The first sowing date represented current sowing dates appropriate for each location. The second sowing date was later and represents the approximate date when the mean air temperature from the climate scenarios is the same as the mean air temperature from the actual climate data at the current sowing dates. The yield and percent kernel nitrogen content using the two climate scenarios generally decrease going from the sub-humid eastern to the semi arid western parts of Nebraska. Results from these simulations indicate that yield and percent kernel nitrogen content using the two climate scenarios could not both be maintained at levels currently simulated. Protein content (directly related to kernel nitrogen content) and end-use quality are the primary determinants for the use of hard red winter wheat in baked goods. Nitrogen management and new cultivars, which can enhance the uptake and translocation of nitrogen, will be proactive steps to meet the challenges of global climate change as represented by these climate scenarios.     

西藏高原南部雅鲁藏布江缝合带地区地壳电性结构研究

为了探测西藏高原南部雅鲁藏布江缝合带地区地壳浅部和深部构造沿东西和南北方向的变化特征,在雅鲁藏布江缝合带地区布设了三条南北向剖面(错那—墨竹工卡、亚东—雪古拉、吉隆—措勤),采用超宽频带大地电磁测深方法进行了地壳、上地幔电性结构探测研究,发现该区主要电性结构特征为：1. 雅鲁藏布江缝合带附近表层发育大规模的高阻体,岩体延伸最深达30km以上,是冈底斯花岗岩体的反映. 2. 雅鲁藏布江缝合带的南部有小规模的良导体,在其下方和北侧发育有大规模良导体. 3. 沿剖面从南往北壳内普遍发育良导体,各良导体主体间是不连续的,规模逐渐增大,总体北倾,在缝合带附近产状较陡. 4. 在雅鲁藏布江缝合带附近良导体由西往东规模逐渐增大,导电性逐渐变好,相对雅鲁藏布江在剖面上的位置逐渐南移. 这些重要的电性特征可能是印度板块向北俯冲所形成的,深部大规模的良导体特征沿东西向的差异可能是板块碰撞引起物质沿东西向运移作用的结果.     

Crustal thickening and crustal extension as controls on the evolution of the drainage network of the central Swiss Alps between 30 Ma and the present: constraints from the stratigraphy of the North Alpine Foreland Basin and the structural evolution of the Alps

The combined information about sedimentary petrography from the North Alpine Foreland Basin and structural geology from the Alps allows a qualitative reconstruction of the drainage network of the central Swiss Alps between 30 Ma and the present. This study suggests that crustal thickening and crustal thinning significantly controlled the location of the drainage divide. It also reveals the possible controls of crustal thickening/thinning on the change of the orientation of the drainage network from across-strike between 30 and 14 Ma to along-strike thereafter. Initial crustal thickening in the rear of the wedge is considered to have formed the drainage divide between north and south at 30 Ma. Because the location of crustal thickening shifted from east to west between ≈30–20 Ma, the catchment areas of the eastern dispersal systems reached further south than those of the western Alpine palaeorivers for the same time slice. Similarly, the same crustal dynamics appear to have controlled two phases of denudation that are reflected in the Molasse Basin by petrographic trends. Uplift in the rear of the wedge caused the Alpine palaeorivers to expand further southward. This is reflected in the foreland basin by increasing admixture of detritus from structurally higher units. However, tectonic quiescence in the rear of the wedge allowed the Alpine palaeorivers to cut down into the Alpine edifice, resulting in an increase of detritus from structurally lower units. Whereas uplift in the rear of the wedge was responsible for initiation of the Alpine drainage systems, underplating of the external massifs some 50 km further north is thought to have caused along-strike deviation of the major Alpine palaeorivers. Besides crustal thickening, extension in the rear of the wedge appears to have significantly controlled the evolution of the drainage network of the western Swiss Alps. Slip along the Simplon detachment fault exposed the core of the Lepontine dome, and caused a 50-km-northward shift of the drainage divide.     

Calcium and magnesium isotope systematics in rivers draining the Himalaya-Tibetan-Plateau region: Lithological or fractionation control?

In rivers draining the Himalaya-Tibetan-Plateau region, the ²⁶Mg/²⁴Mg ratio has a range of 2‰ and the ⁴⁴Ca/⁴²Ca ratio has a range of 0.6‰. The average δ²⁶Mg values of tributaries from each of the main lithotectonic units (Tethyan Sedimentary Series (TSS), High Himalayan Crystalline Series (HHCS) and Lesser Himalayan Series (LHS)) are within 2 standard deviation analytical uncertainty (0.14‰). The consistency of average riverine δ²⁶Mg values is in contrast to the main rock types (limestone, dolostone and silicate) which range in their average δ²⁶Mg values by more than 2‰. Tributaries draining the dolostones of the LHS differ in their values compared to tributaries from the TSS and HHCS. The chemistry of these river waters is strongly influenced by dolostone (solute Mg/Ca close to unity) and both δ²⁶Mg (−1.31‰) and (0.64‰) values are within analytical uncertainty of the LHS dolostone. These are the most elevated values in rivers and rock reported so far demonstrating that both riverine and bedrock values may show greater variability than previously thought.Although rivers draining TSS limestone have the lowest values at −1.41 and 0.42‰, respectively, both are offset to higher values compared to bedrock TSS limestone. The average δ²⁶Mg value of rivers draining mainly silicate rock of the HHCS is −1.25‰, lower by 0.63‰ than the average silicate rock. These differences are consistent with a fractionation of δ²⁶Mg values during silicate weathering. Given that the proportion of Mg exported from the Himalaya as solute Mg is small, the difference in ²⁶Mg/²⁴Mg ratios between silicate rock and solute Mg reflects the ²⁶Mg/²⁴Mg isotopic fractionation factor () between silicate and dissolved Mg during incongruent silicate weathering. The value of of 0.99937 implies that in the TSS, solute Mg is primarily derived from silicate weathering, whereas the source of Ca is overwhelmingly derived from carbonate weathering. The average value in HHCS rivers is within uncertainty of silicate rock at 0.39‰. The widespread hot springs of the High Himalaya have an average δ²⁶Mg value of −0.46‰ and an average value of 0.5‰, distinct from riverine values for δ²⁶Mg but similar to riverine values. Although rivers draining each major rock type have and δ²⁶Mg values in part inherited from bedrock, there is no correlation with proxies for carbonate or silicate lithology such as Na/Ca ratios, suggesting that Ca and Mg are in part recycled. However, in spite of the vast contrast in vegetation density between the arid Tibetan Plateau and the tropical Lesser Himalaya, the isotopic fractionation factor for Ca and Mg between solute and rocks are not systematically different suggesting that vegetation may only recycle a small amount of Ca and Mg in these catchments.The discrepancy between solute and solid Ca and Mg isotope ratios in these rivers from diverse weathering environments highlight our lack of understanding concerning the origin and subsequent path of Ca and Mg, bound as minerals in rock, and released as cations in rivers. The fractionation of Ca and Mg isotope ratios may prove useful for tracing mechanisms of chemical alteration. Ca isotope ratios of solute riverine Ca show a greater variability than previously acknowledged. The variability of Ca isotope ratios in modern rivers will need to be better quantified and accounted for in future models of global Ca cycling, if past variations in oceanic Ca isotope ratios are to be of use in constraining the past carbon cycle.     

Markov chain Monte Carlo methods for conditioning a permeability field to pressure data

Generating one realization of a random permeability field that is consistent with observed pressure data and a known variogram model is not a difficult problem. If, however, one wants to investigate the uncertainty of reservior behavior, one must generate a large number of realizations and ensure that the distribution of realizations properly reflects the uncertainty in reservoir properties. The most widely used method for conditioning permeability fields to production data has been the method of simulated annealing, in which practitioners attempt to minimize the difference between the ’ ’true and simulated production data, and “true” and simulated variograms. Unfortunately, the meaning of the resulting realization is not clear and the method can be extremely slow. In this paper, we present an alternative approach to generating realizations that are conditional to pressure data, focusing on the distribution of realizations and on the efficiency of the method. Under certain conditions that can be verified easily, the Markov chain Monte Carlo method is known to produce states whose frequencies of appearance correspond to a given probability distribution, so we use this method to generate the realizations. To make the method more efficient, we perturb the states in such a way that the variogram is satisfied automatically and the pressure data are approximately matched at every step. These perturbations make use of sensitivity coefficients calculated from the reservoir simulator.     

Spatial Variation in Diameter Structures of Forests in Lassen Volcanic National Park,California*

Variation in seedling/sapling densities and stand diameter forms for six coniferous tree species is related to stand structural development and to elevation and topography in Lassen Volcanic National Park, California. Understory density patterns reflect differences in species tolerance; densities decrease with stand development for shade intolerant pines, but increase for shade tolerant firs and mountain hemlock. Pine species exhibit reverse-J diamter structures on recently disturbed sites, and decreaser and random forms elsewhere. More tolerant fir species show topographically mediated patterns of diameter structure, with reverse-J form common on northerly exposures and upland sites, but with decreaser and random forms on more xeric slopes. Interaction among species tolerance, environmental setting, and disturbance history yields a complex mosaic of stand diameter structures in the Lassen landscape.     

Performance of a mixed-waste landfill amid geologic uncertainty—learning from a case study: Altamont Hills,California, U.S.A.

The Jharia coalfield is the most important and active minig region; it experiences groundwater inflow and affects groundwater levels in overlying aquifers, and it provides the basis for a conceptual model of the hydrogeological impacts of coal mining. The several sandstone aquifers of the overburden are separated by aquitards that limit vertical hydraulic connection, but the inflow responds to seasonal events and seems to be linked to shallow groundwater behavior. The mine drainage behavior suggests a hydraulic connection between the mine and the shallower groundwater system. The greatest declines are directly above the panels, with an immediate response to coal mining. The inflow is localized by natural and induced fracture zones and is mostly into recent workings. The groundwater behavior is controlled by hydraulic property changes caused by mine-induced fracturing. The hydrological and chemical qualities of the shallow groundwater regime in 13 mining collieries in Mukunda Block have been investigated. Water samples collected from 30 shallow monitoring dug wells were chosen for the study. Rainfall, runoff, and infiltration rates have been calculated in the area. The water-quality plottings were used to interpret the distribution of individual chemical parameters and in predicting the water quality. The underground mine water has been classified as: (1) unconfined groundwater in the calcareous siltstone and sandstone—its composition is Na, Ca, SO₄ and Na-MgHCO₃ with moderate total dissolved solids (TDS) 200–1480 ppm; (2) the deep groundwater originating from the coal seams and associated sediments in the near-surface environments—this is a Na-HCO₃ water with higher TDS; and (3) spoil dump waters are essentially Na-HCO₃ with high TDS. This article presents some hydrologic results and conclusions relating to the hydrogeological and environmental impacts of the coal mining in the Jharia coalfield.