Zircon U–Pb ages and Hf isotope compositions of migmatites from the North Qinling terrane and their geological implications

Migmatites are predominant in the North Qinling (NQ) orogen, but their formation ages are poorly constrained. This paper presents a combined study of cathodoluminescence imaging, U–Pb age, trace element and Hf isotopes of zircon in migmatites from the NQ unit. In the migmatites, most zircon grains occur as new, homogeneous crystals, while some are present as overgrowth rims around inherited cores. Morphological and trace element features suggest that the zircon crystals are metamorphic and formed during partial melting. The inherited cores have oscillatory zoning and yield U–Pb ages of c. 900 Ma, representing their protolith ages. The early Neoproterozoic protoliths probably formed in an active continental margin, being a response to the assembly of the supercontinent Rodinia. The migmatite zircon yields Hf model ages of 1911 ± 20 to 990 ± 22 Ma, indicating that the protoliths were derived from reworking of Palaeoproterozoic to Neoproterozoic crustal materials. The anatexis zircon yields formation ages ranging from 455 ± 5 to 420 ± 4 Ma, with a peak at c. 435 Ma. Combined with previous results, we suggest that the migmatization of the NQ terrane occurred at c. 455–400 Ma. The migmatization was c. 50 Ma later than the c. 490 Ma ultra‐high‐P (UHP) metamorphism, indicating that they occurred in two independent tectonic events. By contrast, the migmatization was coeval with the granulite facies metamorphism and the granitic magmatism in the NQ unit, which collectively argue for their formation due to the northward subduction of the Shangdan Ocean. UHP rocks were distributed mainly along the northern margin and occasionally in the inner part of the NQ unit, indicating that they were exhumed along the northern edge and detached from the basement by the subsequent migmatization process.     

Hydrologic response of engineered media in living roofs and bioretention to large rainfalls: experiments and modeling

The hydrologic response of engineered media plays an important role in determining a stormwater control measure's ability to reduce runoff volume, flow rate, timing, and pollutant loads. Five engineered media, typical of living roof and bioretention stormwater control measures, were investigated in laboratory column experiments for their hydrologic responses to steady, large inflow rates. The inflow, medium water content response, and outflow were all measured. The water flow mechanism (uniform flow vs. preferential flow) was investigated by analyzing medium water content response in terms of timing, magnitude, and sequence with depth. Modeling the hydrologic process was conducted in the HYDRUS‐1D software, applying the Richards equation for uniform flow modeling, and a mobile–immobile model for preferential flow modeling. Uniform flow existed in most cases, including all initially dry living roof media with bimodal pore size distributions and one bioretention medium with unimodal pore size distribution. The Richards equation can predict the outflow hydrograph reasonably well for uniform flow conditions when medium hydraulic properties are adequately represented by appropriate functions. Preferential flow was found in two media with bimodal pore size distributions. The occurrence of preferential flow is more likely due to the interaction between the bimodal pore structure and the initial water content rather than the large inflow rate.     

Integration of aquifer geology,groundwater flow and arsenic distribution in deltaic aquifers – A unifying concept

Groundwater arsenic (As) presents a public health risk of great magnitude in densely populated Asian delta regions, most acutely in the Bengal Basin (West Bengal, India and Bangladesh). Research has focused on the sources, mobilisation, and heterogeneity of groundwater As, but a consistent explanation of As distribution from local to basin scale remains elusive. We show for the Bengal Aquifer System that the numerous, discontinuous silt‐clay layers together with surface topography impose a hierarchical pattern of groundwater flow, which constrains As penetration into the aquifer and controls its redistribution towards discharge zones, where it is re‐sequestered to solid phases. This is particularly so for the discrete periods of As release to groundwater in the shallow subsurface associated with sea level high‐stand conditions of Quaternary inter‐glacial periods. We propose a hypothesis concerning groundwater flow ( S ilt‐clay layers I mpose H ierarchical groundwater flow patterns constraining A rsenic progression [SIHA]), which links consensus views on the As source and history of sedimentation in the basin to the variety of spatial and depth distributions of groundwater As reported in the literature. SIHA reconciles apparent inconsistencies between independent, in some cases contrasting, field observations. We infer that lithological and topographic controls on groundwater flow, inherent to SIHA, apply more generally to deltaic aquifers elsewhere. The analysis suggests that groundwater As may persist in the aquifers of Asian deltas over thousands of years, but in certain regions, particularly at deeper levels, As will not exceed low background concentrations unless groundwater flow systems are short‐circuited by excessive pumping.     

Understory and small trees contribute importantly to stemflow of a lower montane cloud forest

Stemflow (Sf) measurements in tropical rain and montane forests dominated by large trees rarely include the understory and small trees. In this study, contributions of lower (1‐ to 2‐m height) and upper (>2‐m height and <5‐cm diameter at breast height [DBH]) woody understory, small trees (5 < DBH < 10 cm), and canopy trees (>10‐cm DBH) to Sf per unit ground area (Sf_a) of a Mexican lower montane cloud forest were quantified for 32 days with rainfall (P) during the 2014 wet season. Rainfall, stemflow yield (Sf_y), vegetation height, density, and basal area were measured. Subsequently, stemflow funneling ratios (SFRs) were calculated, and three common methods to scale up Sf_y from individual trees to the stand level (tree‐Sf_y correlation, P‐Sf_y correlation, and mean‐Sf_y extrapolation) were used to calculate Sf_a. Understory woody plants, small trees, and upper canopy trees represented 96%, 2%, and 2%, respectively, of the total density. Upper canopy trees had the lowest SFRs (1.6 ± 0.5 Standard Error (SE) on average), although the lower understory had the highest (36.1 ± 6.4). Small trees and upper understory presented similar SFRs (22.9 ± 5.4 and 20.2 ± 3.9, respectively). Different Sf scaling methods generally yielded similar results. Overall Sf_a during the study period was 22.7 mm (4.5% of rainfall), to which the understory contributed 70.1% (15.9 mm), small trees 10.6% (2.4 mm), and upper canopy trees 19.3% (4.4 mm). Our results strongly suggest that for humid tropical forests with dense understory of woody plants and small trees, Sf of these groups should be measured to avoid an underestimation of overall Sf at the stand level.     

Using hysteretic behaviour and hydrograph classification to identify hydrological function across the “hillslope–depression–stream” continuum in a karst catchment

In cockpit karst landscapes, fluxes from upland areas contribute large volumes of water to low-lying depressions and stream flow. Hydrograph hysteresis and similarity between monitoring sites is important for understanding the space–time variability of hydrologic responses across the “hillslope–depression–stream” continuum. In this study, the hysteretic feature of hydrographs was assessed by characterizing the loop-like relationships between responses at upstream sites relative to subsurface discharge at the outlet of a small karst catchment. A classification of hydrograph responses based on the multi-scale smoothing Kernel -derived distance classifies the hydrograph responses on the basis of similarities between hillslope and depression sites, and those at the catchment outlet. Results demonstrate that the temporal and spatial variability of hydrograph hysteresis and similarity between hillslope flow and outlet stream flow can be explained by the local heterogeneity of depression aquifer. Large depression storage deficits emerging in the highly heterogeneous aquifer produce strong hysteresis and multiple relationships of upstream hydrographs relative to the outlet subsurface discharge. In contrast, when depression storage deficits are filled during consecutive rainfall events, depression hydrographs at the high permeability sites are almost synchronous or exhibit a monotonous function with the hydrographs at the outlet. This reduced hydrograph hysteresis enhances preferential flow paths in fractured rocks and conduits that can accelerate the hillslope flow to the outlet. Therefore, classification of hydrograph similarities between any upstream sites and the catchment outlet can help to identify the dominant hydrological functions in the heterogeneous karst catchment.     

The June 2013 Alberta catastrophic flooding event – part 2: fine‐scale precipitation and associated features

Data obtained from a variety of sources including the Canadian Lightning Detection Network, weather radars, weather stations and operational numerical weather model analyses were used to address the evolution of precipitation during the June 2013 southern Alberta flood. The event was linked to a mid‐level closed low pressure system to the west of the region and a surface low pressure region initially to its south. This configuration brought warm, moist unstable air into the region that led to dramatic, organized convection with an abundance of lightning and some hail. Such conditions occurred in the southern parts of the region whereas the northern parts were devoid of lightning. Initially, precipitation rates were high (extreme 15‐min rainfall rates up to 102 mm h^?1 were measured) but decreased to lower values as the precipitation shifted to long‐lived stratiform conditions. Both the convective and stratiform precipitation components were affected by the topography. Similar flooding events, such as June 2002, have occurred over this region although the 2002 event was colder and precipitation was not associated with substantial convection over southwest Alberta. Copyright © 2016 Her Majesty the Queen in Right of Canada. Hydrological Processes. © John Wiley & Sons, Ltd.     

P–T–t–d evolution of orogenic middle crust of the Roc de Frausa Massif (Eastern Pyrenees): a result of horizontal crustal flow and Carboniferous doming?

Structural, petrological and textural studies are combined with phase equilibria modelling of metapelites from different structural levels of the Roc de Frausa Massif in the Eastern Pyrenees. The pre‐Variscan lithological succession is divided into the Upper, Intermediate and Lower series by two orthogneiss sheets and intruded by Variscan igneous rocks. Structural analysis reveals two phases of Variscan deformation. D1 is marked by tight to isoclinal small‐scale folds and an associated flat‐lying foliation (S1) that affects the whole crustal section. D2 structures are characterized by tight upright folds facing to the NW with steep NE–SW axial planes. D2 heterogeneously reworks the D1 fabrics, leading to an almost complete transposition into a sub‐vertical foliation (S2) in the high‐grade metamorphic domain. All structures are affected by late open to tight, steeply inclined south‐verging NW–SE folds (F3) compatible with steep greenschist facies dextral shear zones of probable Alpine age. In the micaschists of the Upper series, andalusite and sillimanite grew during the formation of the S1 foliation indicating heating from 580 to 640 °C associated with an increase in pressure. Subsequent static growth of cordierite points to post‐D1 decompression. In the Intermediate series, a sillimanite–biotite–muscovite‐bearing assemblage that is parallel to the S1 fabric is statically overgrown by cordierite and K‐feldspar. This sequence points to ~1 kbar of post‐D1 decompression at 630–650 °C. The Intermediate series is intruded by a gabbro–diorite stock that has an aureole marked by widespread migmatization. In the aureole, the migmatitic S1 foliation is defined by the assemblage biotite–sillimanite–K‐feldspar–garnet. The microstructural relationships and garnet zoning are compatible with the D1 pressure peak at ~7.5 kbar and ~750 °C. Late‐ to post‐S2 cordierite growth implies that F2 folds and the associated S2 axial planar leucosomes developed during nearly isothermal decompression to <5 kbar. The Lower series migmatites form a composite S1–S2 fabric; the garnet‐bearing assemblage suggests peak P–T conditions of >5 kbar at suprasolidus conditions. Almost complete consumption of garnet and late cordierite growth points to post‐D2 equilibration at <4 kbar and <750 °C. The early metamorphic history associated with the S1 fabric is interpreted as a result of horizontal middle crustal flow associated with progressive heating and possible burial. The upright F2 folding and S2 foliation are associated with a pressure decrease coeval with the intrusion of mafic magma at mid‐crustal levels. The D2 tectono‐metamorphic evolution may be explained by a crustal‐scale doming associated with emplacement of mafic magmas into the core of the dome.     

A network theory approach for a better understanding of overland flow connectivity

Hydrological connectivity describes the physical coupling (linkages) of different elements within a landscape regarding (sub‐) surface flows. A firm understanding of hydrological connectivity is important for catchment management applications, for example, habitat and species protection, and for flood resistance and resilience improvement. Thinking about (geomorphological) systems as networks can lead to new insights, which has also been recognized within the scientific community, seeing the recent increase in the use of network (graph) theory within the geosciences. Network theory supports the analysis and understanding of complex systems by providing data structures for modelling objects and their linkages, and a versatile toolbox to quantitatively appraise network structure and properties. The objective of this study was to characterize and quantify overland flow connectivity dynamics on hillslopes in a humid sub‐Mediterranean environment by using a combination of high‐resolution digital‐terrain models, overland flow sensors and a network approach. Results showed that there are significant differences between overland flow connectivity on agricultural areas and semi‐natural shrubs areas. Significant positive correlations between connectivity and precipitation characteristics were found. Significant negative correlations between connectivity and soil moisture were found, most likely because of soil water repellency and/or soil surface crusting. The combination of structural networks and dynamic networks for determining potential connectivity and actual connectivity proved a powerful tool for analysing overland flow connectivity. Copyright © 2016 John Wiley & Sons, Ltd.     

A nonstationary extreme value distribution for analysing the cessation of karst spring discharge

The effects of climate change and population growth in recent decades are leading us to consider their combined and potentially extreme consequences, particularly regarding hydrological processes, which can be modeled using a generalized extreme value (GEV) distribution. Most of the GEV models were based on a stationary assumption for hydrological processes, in contrast to the nonstationary reality due to climate change and human activities. In this paper, we present the nonstationary generalized extreme value (NSGEV) distribution and use it to investigate the risk of Niangziguan Springs discharge decreasing to zero. Rather than assuming the location, scale, and shape parameters to be constant as one might do for a stationary GEV distribution analysis, the NSGEV approach can reflect the dynamic processes by defining the GEV parameters as functions of time. Because most of the GEV model is designed to evaluate maxima (e.g. flooding, represented by positive numbers), and spring discharge cessation is a ?minima’, we deduced an NSGEV model for minima by applying opposite numbers, i.e. negative instead of positive numbers. The results of the model application to Niangziguan Springs showed that the probability of zero discharge at Niangziguan Springs will be 1/80 in 2025, and 1/10 in 2030. After 2025, the rate of decrease in spring discharge will accelerate, and the probability that Niangziguan Springs will cease flowing will dramatically increase. The NSGEV model is a robust method for analysing karst spring discharge. Copyright © 2013 John Wiley & Sons, Ltd.     

GEOPHYSICAL EXPLORATION

正20141074Bao Xijie(Research Institute of Exploration and Development,Daqing Oilfield Company,PetroChina,Daqing 163712,China)Gather Optimal Processing and Application Effect of Prestack AVA Instantaneous Inversion