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.     

Application of topo-edaphic factors and remotely sensed vegetation indices to enhance biomass estimation in a heterogeneous landscape in the Eastern Arc Mountains of Tanzania

Estimating tropical biomass is critical for establishment of conservation inventories and landscape monitoring. However, monitoring biomass in a complex and dynamic environment using traditional methods is challenging. Recently, biomass estimates based on remotely sensed data and ecological variables have shown great potential. The present study explored the utility of remotely sensed data and topo-edaphic factors to improve biomass estimation in the Eastern Arc Mountains of Tanzania. Twenty-nine vegetation indices were calculated from RapidEye data, while topo-edaphic factors were taken from field measurements. Results showed that using topo-edaphic variables or vegetation indices, biomass could be predicted with an R² of 0.4. A combination of topo-edaphic variables and vegetation indices improved the prediction accuracy to an R² of 0.6. Results further showed a decrease in biomass estimates from 1162 ton ha^?1 in 1980 to 285.38 ton ha^?1 in 2012. This study demonstrates the value of combining remotely sensed data with topo-edaphic variables in biomass estimation.     

基于非均质地层参数表征的颗粒流模型构建

非均质地层参数的获取与表征方法是浅层地质体改造的重要依据。基于非均质地层的地质建模,表征具有随机不确定性和模糊不确定性的土体表观参数和物理力学参数。利用 Weibull 分布统计描述局部地层信息,并应用Diamond-Square分形插值方法,结合差分盒子维计算非均质地层分形维数,进一步合理演算整体地层参数。采用PFC3D数值模拟软件,结合马尔科夫链蒙特卡洛方法采样数据,对非均质地层块石颗粒及土颗粒进行建模。通过非均质地层模拟结果,发现其破坏过程与应力应变曲线要比均质土体复杂得多,由于块石大颗粒的含量、分布特性对土样力学性质影响显著。本文建模方法有助于非均质地层参数表征研究。