Metamorphic inheritance of Rheic passive margin evolution and its early‐Variscan overprint in the Teplá‐Barrandian Unit,Bohemian Massif

Rift‐related regional metamorphism of passive margins is usually difficult to observe on the surface, mainly due to its strong metamorphic overprint during the subsequent orogenic processes that cause its exposure. However, recognition of such a pre‐orogenic evolution is achievable by careful characterization of the polyphase tectono‐metamorphic record of the orogenic upper plate. A multidisciplinary approach, involving metamorphic petrology, P–T modelling, structural geology and in situ U‐Pb monazite geochronology using laser‐ablation split‐stream inductively coupled plasma mass spectrometry, was applied to unravel the polyphase tectono‐metamorphic record of metapelites at the western margin of the Teplá‐Barrandian domain in the Bohemian Massif. The study resulted in discovery of three tectono‐metamorphic events. The oldest event M1 is LP–HT regional metamorphism with a geothermal gradient between 30 and 50 °C km^?1, peak temperatures up to 650 °C and of Cambro‐Ordovician age (c. 485 Ma). The M1 event was followed by M2‐D2, which is characterized by a Barrovian sequence of minerals from biotite to kyanite and a geothermal gradient of 20–25 °C km^?1. D2‐M2 is associated with a vertical fabric S2 and was dated as Devonian (c. 375 Ma). Finally, the vertical fabric S2 was overprinted by a D3‐M3 event that formed sillimanite to chlorite bearing gently inclined fabric S3 also of Devonian age. The high geothermal gradient of the M1 event can be explained as the result of an extensional, rift‐related tectonic setting. In addition, restoration of the deep architecture and polarity of the extended domain before the Devonian history – together with the supracrustal sedimentary and magmatic record – lead us to propose a model for formation of an Ordovician passive continental margin. The subsequent Devonian evolution is interpreted as horizontal shortening of the passive margin at the beginning of Variscan convergence, followed by detachment‐accommodated exhumation of lower‐crustal rocks. Both Devonian shortening and detachment occurred in the upper plate of a Devonian subduction zone. The tectonic evolution presented in this article modifies previous models of the tectonic history of the western margin of the Teplá‐Barrandian domain, and also put constraints on the evolution of the southern margin of the Rheic ocean from the passive margin formation to the early phases of Variscan orogeny.     

Petrology,phase equilibria modelling and zircon U–Pb geochronology of Paleoproterozoic mafic granulites from the Fuping Complex,North China Craton

The Fuping Complex is one of the important basement terranes within the central segment of the Trans‐North China Orogen (TNCO) where mafic granulites are exposed as boudins within tonalite–trondhjemite–granodiorite (TTG) gneisses. Garnet in these granulites shows compositional zoning with homogeneous cores formed in the peak metamorphic stage, surrounded by thin rims with an increase in almandine and decrease in grossular contents suggesting retrograde decompression and cooling. Petrological and phase equilibria studies including pseudosection calculation using thermocalc define a clockwise P–T path. The peak mineral assemblages comprise garnet+clinopyroxene+amphibole+quartz+plagioclase+K‐feldspar+ilmenite±orthopyroxene±magnetite, with metamorphic P–T conditions estimated at 8.2–9.2 kbar, 870–882 °C (15FP‐02), 9.6–11.3 kbar, 855–870 °C (15FP‐03) and 9.7–10.5 kbar, 880–900 °C (15FP‐06) respectively. The pseudosections for the subsequent retrograde stages based on relatively higher H₂O contents from P/T–M(H₂O) diagrams define the retrograde P–T conditions of <6.1 kbar, <795 °C (15FP‐02), 5.6–5.8 kbar, <795 °C (15FP‐03), and <9 kbar, <865 °C (15FP‐06) respectively. Data from LA‐ICP‐MS zircon U–Pb dating show that the mafic dyke protoliths of the granulite were emplaced at c. 2327 Ma. The metamorphic zircon shows two groups of ages at 1.96–1.90 Ga (peak at 1.93–1.92 Ga) and 1.89–1.80 Ga (peak at 1.86–1.83 Ga), consistent with the two metamorphic events widely reported from different segments of the TNCO. The 1.93–1.92 Ga ages are considered to date the peak granulite facies metamorphism, whereas the 1.86–1.83 Ga ages are correlated with the retrograde event. Thus, the collisional assembly of the major crustal blocks in the North China Craton (NCC) might have occurred during 1.93–1.90 Ga, marking the final cratonization of the NCC.     

The P–T–t evolution of the exhumed Himalayan metamorphic core in the Likhu Khola region,East Central Nepal

The recent identification of multiple strike‐parallel discontinuities within the exhumed Himalayan metamorphic core has helped revise the understanding of convergence accommodation processes within the former mid‐crust exposed in the Himalaya. Whilst the significance of these discontinuities to the overall development of the mountain belt is still being investigated, their identification and characterization has become important for potential correlations across regions, and for constraining the kinematic framework of the mid‐crust. The result of new phase equilibria modelling, trace element analysis and high‐precision Lu–Hf garnet dating of the metapelites from the Likhu Khola region in east central Nepal, combined with the previously published monazite petrochronology data confirms the presence of one of such cryptic thrust‐sense tectonometamorphic discontinuities within the lower portion of the exhumed metamorphic core and provides new constraints on the P–T estimates for that region. The location of the discontinuity is marked by an abrupt change in the nature of P–T–t paths of the rocks across it. The rocks in the footwall are characterized by a prograde burial P–T path with peak metamorphic conditions of ~660°C and ~9.5 kbar likely in the mid‐to‐late Miocene, which are overlain by the hanging wall rocks, that preserve retrograde P–T paths with P–T conditions of >700°C and ~7 kbar in the early Miocene. The occurrence of this thrust‐sense structure that separates rock units with unique metamorphic histories is compatible with orogenic models that identify a spatial and temporal transition from early midcrustal deformation and metamorphism in the deeper hinterland to later deformation and metamorphism towards the shallower foreland of the orogen. Moreover, these observations are comparable with those made across other discontinuities at similar structural levels along the Himalaya, confirming their importance as important orogen‐scale structures.     

Pseudosection modelling and garnet Lu–Hf geochronology of HP amphibole schists constrain the closure of an ocean basin between the northern and southern Lhasa blocks,central Tibet

The P–T–t path of high‐P metamorphic rocks in subduction zones may reveal valuable information regarding the tectonic processes along convergent plate boundaries. Herein, we present a detailed petrological, pseudosection modelling and radiometric dating study of several amphibole schists of oceanic affinity from the Lhasa Block, Tibet. The amphibole schists experienced an overall clockwise P–T path that was marked by post‐P_max heating–decompression and subsequent isothermal decompression following the attainment of peak high‐P and low‐T conditions (~490°C and 1.6 GPa). Pseudosection modelling shows that the amphibole schists underwent water‐unsaturated conditions during prograde metamorphism, and the stability field of the assemblage extends to lower temperatures and higher pressures within the water‐unsaturated condition relative to water‐saturated model along the prograde path. The high‐P amphibole schists were highly reduced during retrograde metamorphism. Precise evaluation of the ferric iron conditions determined from the different compositions of epidote inclusions in garnet and matrix epidote is crucial for a true P–T estimate by garnet isopleth thermobarometry. Lu–Hf isotope analyses on garnet size separates from a garnet‐bearing amphibole schist yield four two‐point garnet–whole‐rock isochron ages from 228.2 ± 1.2 Ma to 224.3 ± 1.2 Ma. These Lu–Hf dates are interpreted to constrain the period of garnet growth and approximate the timing of prograde metamorphism because of the low peak metamorphic temperature of the rock and the well‐preserved Mn/Lu growth zoning in garnet. The majority of zircon U–Pb dates provide no constraints on the timing of metamorphism; however, two concordant U–Pb dates of 222.4 ± 3.9 Ma and 223.3 ± 4.2 Ma were obtained from narrow and uncommon metamorphic rims. Coexistence of zircon and sphene in the samples implies that the metamorphic zircon growth was likely assisted by retrogression of rutile to sphene during exhumation. The near coincident radiometric dates of zircon U–Pb and garnet Lu–Hf indicate rapid burial and exhumation of the amphibole schists, suggesting a closure time of c. 224–223 Ma for the fossil ocean basin between the northern and southern Lhasa blocks.     

A curious case of agreement between conventional thermobarometry and phase equilibria modelling in granulites: New constraints on P–T estimates in the Antarctica segment of the Musgrave–Albany–Fraser–Wilkes Orogen

The Windmill Islands region in Wilkes Land, east Antarctica, preserves granulite facies metamorphic mineral assemblages that yield seemingly comparable P–T estimates from conventional thermobarometry and mineral equilibria modelling. This is uncommon in granulite facies terranes, where conventional thermobarometry and phase equilibria modelling generally produce conflicting P–T estimates because peak mineral compositions tend to be modified by retrograde diffusion processes. In situ U–Pb monazite geochronology and calculated metamorphic phase diagrams show that the Windmill Islands experienced two phases of high thermal gradient metamorphism during the Mesoproterozoic. The first phase of metamorphism is recorded by monazite ages in two widely separated samples and occurred at c. 1,305 Ma. This event was regional in extent, involved crustally derived magmatism and reached conditions of ~3.2–5 kbar and 690–770°C corresponding to very high thermal gradients of >150°C/kbar. The elevated thermal regime is interpreted to reflect a period of extension or increased extension in a back‐arc setting that existed prior to c. 1,330 Ma. The first metamorphic event was overprinted by granulite facies metamorphism at c. 1,180 Ma that was coeval with the intrusion of charnockite. This event involved peak temperatures of ~840–850°C and pressures of ~4–5 kbar. A phase of granitic magmatism at c. 1,250–1,210 Ma, prior to the intrusion of the charnockite, is interpreted to reflect a phase of compression within an overall back‐arc setting. Existing conventional thermobarometry suggests conditions of ~4 kbar and 750°C for M₁ and 4–7 kbar and 750–900°C for M₂. The apparent similarities between the phase equilibria modelling and existing conventional thermobarometry may suggest either that the terrane cooled relatively quickly, or that the P–T ranges obtained from conventional thermobarometry are sufficiently imprecise that they cover the range of P–T conditions obtained in this study. However, without phase equilibria modelling, the veracity of existing conventional P–T estimates cannot be evaluated. The calculated phase diagrams from this study allow the direct comparison of P–T conditions in the Windmill Islands with phase equilibria models from other regions in the Musgrave–Albany–Fraser–Wilkes Orogen. This shows that the metamorphic evolution of the Wilkes Land region is very similar to that of the eastern Albany–Fraser Orogen and Musgrave Province in Australia, and further demonstrates the remarkable consistency in the timing of metamorphism and the thermal gradients along the ~5,000 km strike length of this system.     

Analysis of incorporating groundwater exchanges in hydrological models

In many catchments, the geographical demarcation does not coincide with the limits of the aquifers, so groundwater may be exchanged beyond their topographic boundaries. By studying groundwater exchanges, the natural resources of a catchment can be better assessed, and the divergences between hydrological models and measurements can be explained. The aim of this work is to reveal the importance of including groundwater exchanges in the hydrological modelling of some catchments, using a water balance model. For this purpose, a simple example is conducted. The so‐called parent model scheme is modified to only allow groundwater exchanges, and it is applied to the headwater of the Segura River Basin District, located in the southeast of Spain. This area is selected because groundwater plays an important role in surface hydrology. The results reveal that groundwater exchanges cannot be neglected in some catchments when assessing water resources because their integration in the hydrological model corrects errors in the water balance. Moreover, this paper proves that water balance models are a useful tool for estimating groundwater exchanges between catchments, which can be contrasted with more complex distributed models or isotopic tracers if there is enough information available. Copyright © 2015 John Wiley & Sons, Ltd.     

Assessment of rainfall spatial variability and its influence on runoff modelling: A case study in the Brue catchment,UK

This study explores rainfall spatial variability and its influence on runoff modelling. A novel assessment scheme integrated with coefficient of variance and Moran's I is introduced to describe effective rainfall spatial variability. Coefficient of variance is widely accepted to identify rainfall variability through rainfall intensity, whereas Moran's I reflects rainfall spatial autocorrelation. This new assessment framework combines these two indicators to assess the spatial variability derived from both rainfall intensity and distribution, which are crucial in determining the time and magnitude of runoff generation. Four model structures embedded in the Variable Infiltration Capacity model are adopted for hydrological modelling in the Brue catchment of England. The models are assigned with 1, 3, 8, and 27 hydrological response units, respectively, and diverse rainfall spatial information for 236 events are extracted from 1995. This study investigates the model performance of different partitioning based on rainfall spatial variability through peak volume (Q_p) and time to peak (T_p), along with the rainfall event process. The results show that models associated with dense spatial partitioning are broadly capable of capturing more spatial information with better performance. It is unnecessary to utilize models with high spatial density for simple rainfall events, though they show distinct advantages on complex events. With additional spatial information, Q_p experiences a notable improvement over T_p. Moreover, seasonal patterns signified by the assessment scheme imply the feasibility of seasonal models.     

Evaluation of 2D models for the prediction of surface depression storage using realistic reference values

Depression storage (DS) is the maximum storage of precipitation and runoff in the soil surface at a given slope. The DS is determined by soil roughness that in agricultural soils is largely affected by tillage. The direct measurement of DS is not straightforward because of the natural permeability of the soil. Therefore, DS has generally been estimated from 2D/3D empirical relationships and numerical algorithms based on roughness indexes and height measurements of the soil surface, respectively. The objective of this work was to evaluate the performance of some 2D models for DS, using direct and reliable measurements of DS in an agricultural soil as reference values. The study was carried out in experimental microplots where DS was measured in six situations resulting from the combination of three types of tillage carried out parallel and perpendicular to the main slope. Those data were used as reference to evaluate four empirical models and a numerical method. Longitudinal altitudinal profiles of the relief were obtained by a laser profilometer. Infiltration measurements were carried out before and after tillage. The DS was largely affected by tillage and its direction. Highest values of DS are found on rougher surfaces mainly when macroforms cut off the dominant slope. The empirical models had a limited performance while the numerical method was the most effective, even so, with an important variability. In addition, a correct hydrological management should take into account that each type of soil tillage affects infiltration rate differently. Copyright © 2016 John Wiley & Sons, Ltd.     

Understanding fen hydrology across multiple scales

Fens, which are among the most biodiverse of wetland types in the USA, typically occur in glacial landscapes characterized by geo‐morphologic variability at multiple spatial scales. As a result, the hydrologic systems that sustain fens are complex and not well understood. Traditional approaches for characterizing such systems use simplifying assumptions that cannot adequately capture the impact of variability in geology and topography. In this study, a hierarchical, multi‐scale groundwater modelling approach coupled with a geologic model is used to understand the hydrology of a fen in Michigan. This approach uses high‐resolution data to simulate the multi‐scale topographic and hydrologic framework and lithologic data from more than 8500 boreholes in a statewide water well database to capture the complex geology. A hierarchy of dynamically linked models is developed that simulates groundwater flow at all scales of interest and to delineate the areas that contribute groundwater to the fen. The results show the fen receiving groundwater from multiple sources: an adjacent wetland, local recharge, a nearby lake and a regional groundwater mound. Water from the regional mound flows to an intermediate source before reaching the fen, forming a ‘cascading’ connection, while other sources provide water through ‘direct’ connections. The regional mound is also the source of water to other fens, streams and lakes in this area, thus creating a large, interconnected hydrologic system that sustains the entire ecosystem. In order to sustainably manage such systems, conservation efforts must include both site‐based protection and management, as well as regional protection and management of groundwater source areas. Copyright © 2016 John Wiley & Sons, Ltd.