Pre‐Alpine discordant granitic dikes in the metamorphic core of the Betic Cordillera: tectonic implications

Pegmatite dikes bearing andalusite crosscut foliation S₂ in Alpujarride gneisses and schists. Post‐S₂ andalusite is transposed by a foliation S₃, defined by fibrolite, which affects the dikes. The dikes represent highly differentiated granitic magmas with low REE and Zr contents and a positive Eu anomaly. U‐Pb SHRIMP dating of magmatic zircons provided Pan‐African ages (cores) and late Variscan ages (rims). However, U‐rich rims also provided metamorphic Alpine ages, supporting a polyorogenic tectonometamorphic history for pre‐Mesozoic Alpujarride rocks.     

Does stress transmission in forelands depend on structural style? Distinctive stress magnitudes during Sevier thin‐skinned and Laramide thick‐skinned layer‐parallel shortening in the Bighorn Basin (USA) revealed by stylolite and calcite twinning paleopiezometry

The Sheep Mountain‐Little Sheep Mountain Anticlines, Bighorn Basin (USA) formed as basement‐cored Laramide structures in the formerly undeformed foreland of the thin‐skinned Sevier orogen. We take advantage of the well‐constrained microstructural network there to reconstruct differential stress magnitudes that prevailed during both Sevier and Laramide layer‐parallel shortening (LPS), before the onset of large‐scale folding. Differential stress magnitudes determined from tectonic stylolites are compared and combined to previous stress estimates from calcite twinning paleopiezometry in the same formations. During stress loading related to LPS, differential stress magnitudes transmitted from the distant Sevier thin‐skinned orogen into the sedimentary cover of the Bighorn basin (11–43 MPa) are higher than the differential stress magnitudes accompanying the early stage of LPS related to the thick‐skinned Laramide deformation (2–19 MPa). This study illustrates that the tectonic style of an orogen affects the transmission of early orogenic stress into the stable continental interior.     

Polygonal grain‐based distinct element modeling for mechanical behavior of brittle rock

The microstructure of rock was numerically reproduced by a polygonal grain‐based model, and its mechanical behavior was examined by performing the uniaxial compression test and Brazilian tests via the Universal Distinct Element Code. The numerical results of the model demonstrated good agreement with the experimental results obtained with rock specimens in terms of the stress–strain behavior, strength characteristics, and brittle fracture phenomenon. An encouraging result is that the grain‐based model‐Universal Distinct Element Code model can reproduce a low ratio of tensile to compressive strength of 1/20 to 1/10 without the need for an additional process. This finding is ascribed to the fact that the geometrical features of polygons can effectively capture the effects of angularity, finite rotation, and interlocking of grains that exist in reality. A numerical methodology to monitor the evolution of micro‐cracks was developed, which enabled us to examine the progressive process of the failure and distinguish the contribution of tensile cracking to the process from that of shear cracking. From the observations of the micro‐cracking process in reference to the stress–strain relation, crack initiation stress, and crack damage stress, it can be concluded that the failure process of the model closely resembles the microscopic observations of rock. We also carried out a parametric study to examine the relationships between the microscopic properties and the macroscopic behavior of the model. Depending on the micro‐properties, the model exhibited a variety of responses to the external load in terms of the strength and deformation characteristics, the evolution of micro‐cracks, and the post‐peak behavior. Copyright © 2016 John Wiley & Sons, Ltd.     

The influence of coupled physical swelling and chemical reactions on deformable geomaterials

Coupled thermo‐hydro‐mechanical‐chemical modelling has attracted attention in past decades due to many contemporary geotechnical engineering applications (e.g., waste disposal, carbon capture and storage). However, molecular‐scale interactions within geomaterials (e.g., swelling and dissolution/precipitation) have a significant influence on the mechanical behaviour, yet are rarely incorporated into existing Thermal‐Hydro‐Mechanical‐Chemical (THMC) frameworks. This paper presents a new coupled hydro‐mechanical‐chemical constitutive model to bridge molecular‐scale interactions with macro‐physical deformation by combining the swelling and dissolution/precipitation through an extension of the new mixture‐coupling theory. Entropy analysis of the geomaterial system provides dissipation energy, and Helmholtz free energy gives the relationship between solids and fluids. Numerical simulation is used to compare with the selected recognized models, which demonstrates that the swelling and dissolution/precipitation processes may have a significant influence on the mechanical deformation of the geomaterials.     

Including spatial distribution in a data‐driven rainfall‐runoff model to improve reservoir inflow forecasting in Taiwan

Multi‐step ahead inflow forecasting has a critical role to play in reservoir operation and management in Taiwan during typhoons as statutory legislation requires a minimum of 3‐h warning to be issued before any reservoir releases are made. However, the complex spatial and temporal heterogeneity of typhoon rainfall, coupled with a remote and mountainous physiographic context, makes the development of real‐time rainfall‐runoff models that can accurately predict reservoir inflow several hours ahead of time challenging. Consequently, there is an urgent, operational requirement for models that can enhance reservoir inflow prediction at forecast horizons of more than 3 h. In this paper, we develop a novel semi‐distributed, data‐driven, rainfall‐runoff model for the Shihmen catchment, north Taiwan. A suite of Adaptive Network‐based Fuzzy Inference System solutions is created using various combinations of autoregressive, spatially lumped radar and point‐based rain gauge predictors. Different levels of spatially aggregated radar‐derived rainfall data are used to generate 4, 8 and 12 sub‐catchment input drivers. In general, the semi‐distributed radar rainfall models outperform their less complex counterparts in predictions of reservoir inflow at lead times greater than 3 h. Performance is found to be optimal when spatial aggregation is restricted to four sub‐catchments, with up to 30% improvements in the performance over lumped and point‐based models being evident at 5‐h lead times. The potential benefits of applying semi‐distributed, data‐driven models in reservoir inflow modelling specifically, and hydrological modelling more generally, are thus demonstrated. Copyright © 2012 John Wiley & Sons, Ltd.     

What is the source of baseflow in agriculturally fragmented catchments? Complex groundwater/surface‐water interactions in three tributary catchments of the Wabash River,Indiana, USA

Some conceptual models suggest that baseflow in agriculturally fragmented watersheds may contain little, if any, groundwater. This has critical implications for stream quality and ecosystem functioning. Here, we (a) identify the sources and flowpaths contributing to baseflow using ²²²Rn and ⁸⁷Sr/⁸⁶Sr and (b) quantify mean apparent ages of groundwater and baseflow using multiple isotopic tracers (CFC, SF₆, ³⁶Cl, and ³H) in 4 small (0.08 to 0.64 km²) tributary catchments to the Wabash River in Indiana, USA. ²²²Rn activities and ⁸⁷Sr/⁸⁶Sr ratios indicate that baseflow in 3 catchments is sourced primarily from groundwater; baseflow in the fourth is dominated by a source similar to agricultural run‐off. CFC‐12 data indicate that springs in 1 catchment are discharging significant proportions of water that recharged between 1974 (42 ± 2 years) and 1961 (55 ± 2 years). Those same springs have ³⁶Cl/Cl ratios between 1,381.08 ± 29.37 (×10^?15) and 1,530.64 ± 27.65 (×10^?15) indicating that a substantial proportion of the discharge likely recharged between 1975 (41 years) and 1950 (66 years). Groundwater samples collected from streambed mini‐piezometers in a separate catchment have CFC‐12 concentrations indicating that a large proportion of the recharge occurred between 1948 (68 ± 2 years) and 1950 (66 ± 2 years). Repeat sampling conducted in September 2015 after above‐average summer rainfall did not show significant decreases in mean apparent age. The relatively old ages observed in 3 of the catchments can be explained by geological complexities that are likely present in all 4 catchments, but overwhelmed by flow from the shallow phreatic aquifer in the fourth catchment.     

Chemical and textural equilibration of garnet during amphibolite facies metamorphism: The influence of coupled dissolution–reprecipitation

Metamorphic equilibration requires chemical communication between minerals and may be inhibited through sluggish volume diffusion and or slow rates of dissolution in a fluid phase. Relatively slow diffusion and the perceived robust nature of chemical growth zoning may preclude garnet porphyroblasts from readily participating in low‐temperature amphibolite facies metamorphic reactions. Garnet is widely assumed to be a reactant in staurolite‐isograd reactions, and the evidence for this has been assessed in the Late Proterozoic Dalradian pelitic schists of the Scottish Highlands. The 3D imaging of garnet porphyroblasts in staurolite‐bearing schists reveals a good crystal shape and little evidence of marginal dissolution; however, there is also lack of evidence for the involvement of either chlorite or chloritoid in the reaction. Staurolite forms directly adjacent to the garnet, and its nucleation is strongly associated with deformation of the muscovite‐rich fabrics around the porphyroblasts. “Cloudy” fluid inclusion‐rich garnet forms in both marginal and internal parts of the garnet porphyroblast and is linked both to the production of staurolite and to the introduction of abundant quartz inclusions within the garnet. Such cloudy garnet typically has a Mg‐rich, Mn‐poor composition and is interpreted to have formed during a coupled dissolution–reprecipitation process, triggered by a local influx of fluid. All garnet in the muscovite‐bearing schists present in this area is potentially reactive, irrespective of the garnet composition, but very few of the schists contain staurolite. The staurolite‐producing reaction appears to be substantially overstepped during the relatively high‐pressure Barrovian regional metamorphism reflecting the limited permeability of the schists in peak metamorphic conditions. Fluid influx and hence reaction progress appear to be strongly controlled by subtle differences in deformation history. The remaining garnet fails to achieve chemical equilibrium during the reaction creating distinctive patchy compositional zoning. Such zoning in metamorphic garnet created during coupled dissolution–reprecipitation reactions may be difficult to recognize in higher grade pelites due to subsequent diffusive re‐equilibration. Fundamental assumptions about metamorphic processes are questioned by the lack of chemical equilibrium during this reaction and the restricted permeability of the regional metamorphic pelitic schists. In addition, the partial loss of prograde chemical and textural information from the garnet porphyroblasts cautions against their routine use as a reliable monitor of metamorphic history. However, the partial re‐equilibration of the porphyroblasts during coupled dissolution–reprecipitation opens possibilities of mapping reaction progress in garnet as a means of assessing fluid access during peak metamorphic conditions.     

Experimental determination of REE partition coefficients between zircon,garnet and melt: a key to understanding high‐T crustal processes

The partitioning of rare earth elements (REE) between zircon, garnet and silicate melt was determined using synthetic compositions designed to represent partial melts formed in the lower crust during anatexis. The experiments, performed using internally heated gas pressure vessels at 7 kbar and 900–1000 °C, represent equilibrium partitioning of the middle to heavy REE between zircon and garnet during high‐grade metamorphism in the mid to lower crust. The D_REE (zircon/garnet) values show a clear partitioning signature close to unity from Gd to Lu. Because the light REE have low concentrations in both minerals, values are calculated from strain modelling of the middle to heavy REE experimental data; these results show that zircon is favoured over garnet by up to two orders of magnitude. The resulting general concave‐up shape to the partitioning pattern across the REE reflects the preferential incorporation of middle REE into garnet, with D_Gd (zircon/garnet) ranging from 0.7 to 1.1, D_Ho (zircon/garnet) from 0.4 to 0.7 and D_Lu (zircon/garnet) from 0.6 to 1.3. There is no significant temperature dependence in the zircon–garnet REE partitioning at 7 kbar and 900–1000 °C, suggesting that these values can be applied to the interpretation of zircon–garnet equilibrium and timing relationships in the ultrahigh‐T metamorphism of low‐Ca pelitic and aluminous granulites.