农林复合带一维非静力大气边界层模式及其模拟分析

建立了一个农林复合带地区一维非静力大气边界层能量闭合模式，对1000m以下的大气边界层内的风、温、湿作了24h的预报，并对下垫面3种不同参数化方案(农作物、森林、无植被)的输出结果与实测值进行了分析和比较；同时通过敏感性试验，突出比较了农作物和森林下垫面对大气边界层垂直流场，湍流垂直交换和湍能的影响。结果表明，本模式能改善边界层风速、位温和湿度预报的模拟效果，下垫面植被对边界层气象要素大小和分布有显著的作用，对湍能垂直分布有一定影响。     

黄土高原塬区地表辐射和热量平衡观测与分析

利用2005年夏季黄土高原塬区陆面过程野外试验(LOPEX05)的观测资料,初步分析了甘肃平凉黄土高原塬区地表辐射收支和热量平衡特征。结果显示,黄土高原塬区地面长波辐射大于大气长波辐射,典型晴天、阴天和雨天情况下两者平均差值分别为65,25和8 W.m-2;对于地气能量交换各个分量而言,黄土高原塬上和塬下在相同下垫面下的差别不大,但裸地和有植被的下垫面差别很明显;在白天,潜热在净辐射中所占的比重较大,其次是感热,最后是土壤热通量。对能量平衡中的储存项如热通量板上层土壤的热存储和植被冠层存储进行了估算,结果表明,土壤的热储存项在-30~70W.m-2之间,而植被的热能储存项在-10~25 W.m-2之间。在考虑估算的存储项之后,能量平衡散布图斜率由0.68提高到0.79,相关系数R由0.90提高到0.93,两者分别提高了11.0%和3.0%,并对能量不平衡有明显的改进,说明能量储存项在地表能量闭合中必须考虑。     

Chemical composition of PM2.5 at an urban site of Chengdu in southwestern China

PM2.5 aerosols were sampled in urban Chengdu from April 2009 to January 2010, and their chemical compositions were characterized in detail for elements, water soluble inorganic ions, and carbonaceous matter. The annual average of PM2.5 was 165g m－3, which is generally higher than measurements in other Chinese cities, suggesting serious particulate pollution issues in the city. Water soluble ions contributed 43.5% to the annual total PM2.5 mass, carbonaceous aerosols including elemental carbon and organic carbon contributed 32.0%, and trace elements contributed 13.8%. Distinct daily and seasonal variations were observed in the mass concentrations of PM2.5 and its components, reflecting the seasonal variations of different anthropogenic and natural sources. Weakly acidic to neutral particles were found for PM2.5. Major sources of PM2.5 identified from source apportionment analysis included coal combustion, traffic exhaust, biomass burning, soil dust, and construction dust emissions. The low nitrate: sulfate ratio suggested that stationary emissions were more important than vehicle emissions. The reconstructed masses of ammonium sulfate, ammonium nitrate, particulate carbonaceous matter, and fine soil accounted for 79% of the total measured PM2.5 mass; they also accounted for 92% of the total measured particle scattering.     

控制网最小独立闭合环搜索算法研究

论述了控制网中最小独立闭合环搜索的两种算法：基于生成树的算法和基于迭代加深搜索的算法,在理论上进行了推导,并结合实际算例,对两种算法进行了比较,得出了在某些情况下,后一种算法搜索结果不是太稳定,前一种算法搜索结果稳定的结论。     

Closed system behaviour of argon in osumilite records protracted high‐T metamorphism within the Rogaland–Vest Agder Sector,Norway

Here, we present results of the first ⁴⁰Ar/³⁹Ar dating of osumilite, a high‐T mineral that occurs in some volcanic and high‐grade metamorphic rocks. The metamorphic osumilite studied here is from a metapelitic rock within the Rogaland–Vest Agder Sector, Norway, an area that experienced regional granulite facies metamorphism and subsequent contact metamorphism between 1,100 Ma and 850 Ma. The large grain size (~1 cm) of osumilite in the studied rock, which preserves a nominally anhydrous assemblage, increases the potential for large portions of individual grains to have remained essentially unaffected by the effects of diffusive argon loss, potentially preserving prograde ages. Step‐heating diffusion experiments yielded a maximum activation energy of ~461 kJ/mol and a pre‐exponential factor of ~8.34 × 10⁸ cm²/s for Ar diffusion in osumilite. These parameters correspond to a relatively high closure temperature of ~620°C for a cooling rate of 10°C/Ma in an osumilite crystal with a 175 µm radius. Fragments of osumilite separated from the sample preserve a range of ages between c. 1,070 and 860 Ma. The oldest ages are inferred to date the growth of coarse‐grained osumilite during prograde granulite facies regional metamorphism, which pre‐date contact metamorphism that has historically been ascribed to the growth of osumilite in the region. The majority of fragments record ages between c. 920 and 860 Ma, inferred to reflect the growth of osumilite and/or diffusive argon loss during contact metamorphism. The retention of old ⁴⁰Ar/³⁹Ar dates was facilitated by the low diffusivity of Ar in osumilite (i.e. a closed system), large grain sizes, and anhydrous metamorphic conditions. The ability to date osumilite with the ⁴⁰Ar/³⁹Ar method provides a valuable new thermochronometer that may constrain the timing and duration of high‐T magmatic and metamorphic events.     

Eddy-Covariance Flux Measurements in the Complex Terrain of an Alpine Valley in Switzerland

We measured the surface energy budget of an Alpine grassland in highly complex terrain to explore possibilities and limitations for application of the eddy-covariance technique, also for CO₂ flux measurements, at such non-ideal locations. This paper focuses on the influence of complex terrain on the turbulent energy measurements of a characteristic high Alpine grassland on Crap Alv (Alp Weissenstein) in the Swiss Alps during the growing season 2006. Measurements were carried out on a topographic terrace with a slope of 25^◦ inclination. Flux data quality is assessed via the closure of the energy budget and the quality flag method used within the CarboEurope project. During 93% of the time the wind direction was along the main valley axis (43% upvalley and 50% downvalley directions). During the transition times of the typical twice daily wind direction changes in a mountain valley the fraction of high and good quality flux data reached a minimum of ≈50%, whereas during the early afternoon ≈70% of all records yielded good to highest quality (CarboEurope flags 0 and 1). The overall energy budget closure was 74 ± 2%. An angular correction for the shortwave energy input to the slope improved the energy budget closure slightly to 82 ± 2% for afternoon conditions. In the daily total, the measured turbulent energy fluxes are only underestimated by around 8% of net radiation. In summary, our results suggest that it is possible to yield realistic energy flux measurements under such conditions. We thus argue that the Crap Alv site and similar topographically complex locations with short-statured vegetation should be well suited also for CO₂ flux measurements.     

Large-eddy Simulations of Flow Over Forested Ridges

Large-eddy simulations (LES) of flow over a series of small forested ridges are performed, and compared with numerical simulations using a one-and-a-half order mixing length closure scheme. The qualitative and quantitative similarity between these results provides some confidence in the results of recent analytical and numerical studies of flow over forested hills using first-order mixing length schemes. Time series of model velocities at various locations within the canopy allow the application of various experimental techniques to study the turbulence in the LES. The application of conditional analysis shows that the structure of the turbulence over a forested hill is broadly similar to that over flat ground, with sweeps and ejections dominating. Differences are seen across the hill, particularly associated with regions of mean flow separation and recirculation near the summit and in the lee of the hill. Detailed comparison of derived mixing lengths from the LES with the assumed values used in mixing-length closure schemes show that the mixing length varies with location across the hill and with height in the canopy. This is consistent with previous wind-tunnel measurements, and demonstrates that a constant mixing-length assumption is not strictly valid within the canopy. Despite this, the first-order mixing-length schemes do give similar results both for the mean flow and the turbulence in such situations.     

The Temperature–Humidity Covariance in the Marine Surface Layer: A One-dimensional Analytical Model

An analytical model that predicts how much of the temperature–humidity covariance within the marine atmospheric surface layer (ASL) originates just above the ASL and just near the surface is proposed and tested using observations from the Risø Air Sea Experiment (RASEX). The model is based on a simplified budget for the two-scalar covariance that retains three basic terms: production, dissipation, and vertical transport. Standard second-order closure formulations are employed for the triple moments and the dissipation terms, and the canonical mixing length for the closure model is assumed linear with height (z) from the surface. Despite the poor performance of the gradient–diffusion closure in reproducing the measured triple moment, the overall covariance model was shown to be sufficiently robust to these assumptions. One of the main findings from the analytical treatment is the origin of the asymmetry in how the top and bottom boundary conditions affect the two-scalar covariance in the ASL. The analytical model reveals that ‘bottom-up’ boundary-condition variations scale with \(z^{-\sqrt{a}}\) , while ‘top-down’ variations scale with \(z^{\sqrt{a}}\) , where a is a constant that can be derived from similarity and closure constants. The genesis of this asymmetry stems from the flux-transport term but is modulated by the dissipation, and persists even in the absence of any inhomogeneity in the local production function. It is shown that the local production function acts to adjust the relative proportions of these two boundary conditions with weights that vary with the Obukhov length. The findings here do not provide ‘finality’ to the discussions on the covariance between humidity and temperature or the role of entrainment in modulating the turbulence within the ASL. Rather, they are intended to guide new hypotheses about interpretations of existing field data and identify needs for future field and numerical experiments.     

A One-Dimensional Finite-Element Boundary-Layer Model with a Vertical Adaptive Grid

A one-dimensional atmospheric boundary-layer model is developed using the finite element method and a 1.5-order e-l turbulence closure scheme. A vertical adaptive strategy is implemented, based upon an heuristic error estimator that depends upon properties of the layer, such as the stratification. The model is used to simulate a moderately stratified stable boundary layer as described in the GABLS (GEWEX Atmospheric Boundary-Layer Study, where GEWEX is the Global Water and Energy Cycle Experiment) First Intercomparison Project, and then a more complicated diurnal cycle, as used in the GABLS Second Intercomparison Project. In the stable boundary-layer experiment, it is shown that including the adaptive strategy can improve the performance of the model such that the error in the model is significantly less (greater than an order of magnitude with an effective resolution of 8 m) than that of the model without adaptivity. The model’s turbulence closure scheme and the adaptivity strategy also successfully simulate the different stability regimes present in the diurnal cycle simulation, and represented all of the expected features.     

Effects of non‐uniform traction and specimen height in the direct shear test on stress and deformation in a rock fracture

In the direct shear test (DST), an internal moment is distributed within the rock specimen by non‐coaxial shear loads applied to the specimen, which cause non‐uniform distributions of both the traction on the loading planes and the stress and deformation in the specimen. To examine the validity of the DST for a rock fracture and to clarify the effect of specimen height, both the stress and deformation in a fracture in the DST were analyzed for specimens with three different heights using a three‐dimensional finite element method with quadratic joint elements for a fracture model. The constitutive law of the fracture considers the dependence of the non‐linear behavior of closure on shear displacement and that of shear stiffness on normal stress and was implemented in simulation code to give a conceptional fracture with uniform mechanical properties to extract only the effect of non‐uniform traction on the stress and deformation in the fracture. The results showed that both normal and shear stresses are concentrated near the end edges of the fracture, and these stress concentrations decrease with a decrease in the specimen height according to the magnitude of the moment produced by the non‐coaxial shear loads. Furthermore, although closure is greater near the end edges of the fracture, where normal stress is concentrated, this concentration of closure is not so significant within the range of this study because of the non‐linear behavior of closure, that is, closure does not significantly increase with an increase in normal stress at large normal stresses. Copyright © 2012 John Wiley & Sons, Ltd.