The role of heterogeneous strain in the development and preservation of a polymetamorphic record in high-P granulites, western Canadian Shield

Mafic rocks in the Chipman domain of the Athabasca granulite terrane, western Canadian Shield, provide the first well‐documented record of two distinct high‐P granulite facies events in the same domain in this region. Textural relations and the results of petrological modelling (NCFMASHT system) of mafic granulites are interpreted in terms of a three‐stage tectonometamorphic history. Stage 1 involved development of the assemblage Grt + Cpx + Qtz ± Pl (M₁) from a primary Opx‐bearing igneous precursor at conditions of 1.3 GPa, 850–900 °C. Field and microstructural observations suggest that M₁ developed synchronously with an early S₁ gneissic fabric. Stage 2 is characterized by heterogeneous deformation (D₂) and synkinematic partial retrogression of the peak assemblage to an amphibole‐bearing assemblage (M₂). Stage 3 involved a third phase of deformation and a return to granulite facies conditions marked by the prograde breakdown of amphibole (Amph₂) to produce matrix garnet (Grt_3a) and the coronitic assemblage Cpx_3b + Opx_3b + Ilm_3b + Pl_3b (M_3b) at 1.0 GPa, 800–900 °C. M₁ and M_3b are correlated with 2.55 and 1.9 Ga metamorphic generations of zircon, respectively, which were dated in a separate study. Heterogeneous strain played a crucial role in both the development and preservation of these rare examples of multiple granulite facies events within single samples. Without this fortuitous set of circumstances, the apparent reaction history could have incorrectly led to an interpretation involving a single‐cycle high‐grade event. The detailed P–T–t–D history constructed for these rocks provides the best evidence to date that much of the east Lake Athabasca region experienced long‐term lower crustal residence from 2.55 to 1.9 Ga, and thus the region represents a rare window into the reactivation and ultimate stabilization processes of cratonic lithosphere.     

地下非均匀非饱和带中地下洞室的渗流问题数值模拟--介质参数的灵敏度分析

影响任何一个地下废物处置场所长期行为的一个重要因素是进入地下废物贮置洞室内的渗流量的大小，预测地下洞室中的渗流量是困难的，特别是当地下废物贮置洞室位于非均质非饱和带中，三维数值模型用于研究地下非均质中非饱和流动及流入地下洞室中的渗流量，讨论了非均质模型与均质模型的比较以及数值剖分尺度对计算结果的影响，进入洞室中的渗漏率随着引入系统中入渗量的增加而增大，选用4个参数来衡量渗流场的非均质程度：(1)介质的平均渗透率K0；(2)Van Genuchten参数α；(3)渗透率空间分布相关尺度；(4)渗透率空间分布变化的标准差σ，根据一个随机实现的渗透率分布，通过改变平均渗透率来研究其对流入洞室中渗透量的影响，对一个固定的入渗率而言，流入洞室中的渗漏率将随着Van Genuchten参数α的减小而减小，模拟结果表明流入洞室的渗漏率与介质的平均渗透率相关，即随平均渗透率的增大，流入洞室的渗漏率亦增大，流入洞室的渗漏率还高度依赖于非均质渗透率场的空间分布相关长度与标准差，一个大的相关长度或高的标准差均能导致流入洞室的渗漏率增大。     

Generalized Coarse Graining Procedures for Flow in Porous Media

This paper focuses on heterogeneous soil conductivities and on the impact their resolution has on a solution of the piezometric head equation: owing to spatial variations of the conductivity, the flow properties at larger scales differ from those found for experiments performed at smaller scales. The method of coarse graining is proposed in order to upscale the piezometric head equation on arbitrary intermediate scales. At intermediate scales large scale fluctuations of the conductivities are resolved, whereas small scale fluctuations are smoothed by a partialy spatial filtering procedure. The filtering procedure is performed in Fourier space with the aid of a low-frequency cut-off function. We derive the partially upscaled head equations. In these equations, the impact of the small scale variability is modeled by scale dependent effective conductivities which are determined by additional differential equations. Explicit results for the scale dependent conductivity values are presented in lowest order perturbation theory. The perturbation theory contributions are summed up with using a renormalisation group analysis yielding explicit results for the effective conductivity in isotropic media. Therefore, the results are also valid for highly heterogeneous media. The results are compared with numerical simulations performed by Dykaar and Kitanidis (1992). The method of coarse graining combined by a renormalisation group analysis offers a tool to derive exact and explicit expressions for resolution dependent conductivity values. It is, e.g., relevant for the interpretation of measurement data on different scales and for reduction of grid-block resolution in numerical modeling. This revised version was published online in July 2006 with corrections to the Cover Date.     

非均质土坝稳定性的渗流场和应力场耦合分析

綦于多孔介质的渗流特性和土的非线性本构关系，研究了渗流场与应力场的耦合效应，并针对一实际十坝建立了非线性耦合分析模型，探讨了位于软基之上的非均质土坝的渗流与应力耦合规律。分析结果表明，在渗流作用下，坝休和坝基中渗流与应力的相互作用是不应忽视的，其作用结果将会使土坝、特别足防渗心端附近应力水平和渗透坡降均有不同程度的提高，表明考虑渗流场和应力场耦合作用的计算结果对土坝稳定安全性不利。     

Scaling procedures for heterogeneous unconfined aquifers

Complex flows in heterogeneous confined and unconfined aquifers is a phenomenon that continues to present difficulties in flow mapping and modelling in the field, laboratory, and through numerical simulations. It is often the case with complicated phenomena that transformative scaling and reduction of the problem through symmetry is of great efficacy in the formation of predictive models in both the laboratory and computational settings. A detailed a study of the application of a broad class of Lie scaling transformations on a set of equations representing the groundwater flows in heterogeneous confined and unconfined aquifers has produced a set of scaling relationships between the spatial variables, hydrologic variables, and parameters. The set of scaling transformations preserve the structure of the equations in the sense that the scaling transformations leave the initial‐boundary value system representing the invariant groundwater flows. This theoretical approach elucidates not only the scaling relationships but also the properties that hydrologic variables and parameters must satisfy in order for calling to be possible. Validation of the theory developed is carried out through a series of four numerical simulations using the USGS modflow ‐2005 software package. The results of these experiments demonstrate that the derived scaling transformations can effectively form predictive models of large‐scale phenomena at small scales with negligible error in many cases. Comments on the limitations of the approach and directions for future research are made in the closing sections. Copyright © 2016 John Wiley & Sons, Ltd.     

Ensemble simulation of land evapotranspiration in China based on a multi-forcing and multi-model approach

In order to reduce the uncertainty of offline land surface model(LSM) simulations of land evapotranspiration(ET), we used ensemble simulations based on three meteorological forcing datasets [Princeton, ITPCAS(Institute of Tibetan Plateau Research, Chinese Academy of Sciences), Qian] and four LSMs(BATS, VIC, CLM3.0 and CLM3.5), to explore the trends and spatiotemporal characteristics of ET, as well as the spatiotemporal pattern of ET in response to climate factors over mainland China during 1982–2007. The results showed that various simulations of each member and their arithmetic mean(Ens Mean) could capture the spatial distribution and seasonal pattern of ET sufficiently well, where they exhibited more significant spatial and seasonal variation in the ET compared with observation-based ET estimates(Obs MTE). For the mean annual ET, we found that the BATS forced by Princeton forcing overestimated the annual mean ET compared with Obs MTE for most of the basins in China, whereas the VIC forced by Princeton forcing showed underestimations. By contrast, the Ens Mean was closer to Obs MTE, although the results were underestimated over Southeast China. Furthermore, both the Obs MTE and Ens Mean exhibited a significant increasing trend during 1982–98; whereas after 1998, when the last big EI Ni ?no event occurred, the Ens Mean tended to decrease significantly between 1999 and 2007, although the change was not significant for Obs MTE. Changes in air temperature and shortwave radiation played key roles in the long-term variation in ET over the humid area of China, but precipitation mainly controlled the long-term variation in ET in arid and semi-arid areas of China.     

基于多模式集成技术的地面气温精细化预报

基于ECMWF、JMA、T639、WRF四个数值模式2012年6月1日—9月30日地面气温3—60 h预报资料和郑州加密自动站资料,利用多模式集合平均(EMN)、消除偏差集合平均(BREM)、加权消除偏差集合(WBREM)及多模式超级集合(SUP)4种方法,对2012年8月29日—9月27日郑州城区11个站点地面逐3 h气温进行多模式集成预报试验,采用绝对误差对预报结果进行检验评估,结果表明:在30天的预报期内,BREM、WBREM及SUP对于大多数站气温预报效果有明显改善,而EMN方案对11个站预报效果改善则不太明显;4种方案中,BREM和WBREM预报效果相对较好且稳定,各个站上3—60 h预报的绝对误差均在2℃附近或以下;SUP方案虽然对个别站预报误差较低,但是其预报效果并不稳定,一些站点的个别预报时效误差大于2℃。对于郑州观测站的气温预报而言,4种集成方案20时起报的气温误差明显小于08时起报的误差,并且20时起报的SUP集成方案绝对误差明显小于其他方案的绝对误差。总体而言,BREM、WBREM及SUP三种集成方案能够给郑州精细化预报业务提供较好的参考。