Stable staurolite–cordierite assemblages in K‐poor metapelitic schists in Aston and Hospitalet gneiss domes of the central Pyrenees (France,Andorra)

Staurolite–cordierite assemblages are common in mica schists of the Aston and Hospitalet gneiss domes of the central Axial Zone, Pyrenees (France, Andorra). Within a 200 m wide zone, staurolite, cordierite and andalusite porphyroblasts contain inclusion trails that preserve the same stage of development of a crenulation cleavage, strongly suggesting that all three phases are contemporaneous. Their syntectonic growth occurred during a short period at the beginning of the formation of the dominant schistosity (S₂) of the domes. Staurolite and cordierite touching each other further indicates an equilibrium relationship. Whole‐rock analyses show that some staurolite–cordierite schists are depleted in K₂O compared to post‐Archean shales (PAAS) and amphibolite facies pelites. Analysis of the st‐crd paragenesis in K‐poor schists without muscovite using KFMASH and MnNCKFMASH petrogentic grids, pseudosections and AFM compatibility diagrams predicts stable conditions at pressures of ~3.5 kbar at 575 °C. For metapelites with intermediate X_Mg values (0.7 >  X_Mg >0.48) a ‘muscovite‐out window’ exists from 550–650 °C at 3.5 kbar in the KFMASH system. Conventional thermobarometry (GB‐GASP, AvT‐AvP) and petrogenetic grids show an isobaric P–T path to peak temperatures of ~650 °C, supported by the presence of sillimanite‐K‐feldspar gneiss and migmatites. LP‐HT metamorphism in the Aston dome is related to early Carboniferous (c. 339 Ma) granitic intrusions into the dome core. As metamorphism is directly linked with the formation of the main S₂ schistosity, the temporal relations demonstrated in this study conflict with previous studies which constrained LP‐HT metamorphism and the development of flat‐lying schistosity to the late Carboniferous (315–305 Ma) – at least in the eastern Axial Zone.     

Dating megafaunal extinction on the Pleistocene Darling Downs,eastern Australia: the promise and pitfalls of dating as a test of extinction hypotheses

A key to understanding Late Pleistocene megafaunal extinction dynamics is knowledge of megafaunal ecological response(s) to long-term environmental perturbations. Strategically, that requires targeting fossil deposits that accumulated during glacial and interglacial intervals both before and after human arrival, with subsequent palaeoecological models underpinned by robust and reliable chronologies. Late Pleistocene vertebrate fossil localities from the Darling Downs, eastern Australia, provide stratigraphically-intact, abundant megafaunal sequences, which allows for testing of anthropogenic versus climate change megafauna extinction hypotheses. Each stratigraphic unit at site QML796, Kings Creek Catchment, was previously shown to have had similar sampling potential, and the basal units contain both small-sized taxa (e.g., land snails, frogs, bandicoots, rodents) and megafauna. Importantly, sequential faunal horizons show stepwise decrease in taxonomic diversity with the loss of some, but not all, megafauna in the geographically-small palaeocatchment. The purpose of this paper is to present the results of our intensive, multidisciplinary dating study of the deposits (>40 dates). Dating by means of accelerator mass spectrometry (AMS) ¹⁴C (targeting bone, freshwater molluscs, and charcoal) and thermal ionisation mass spectrometry U/Th (targeting teeth and freshwater molluscs) do not agree with each other and, in the case of AMS ¹⁴C dating, lack internal consistency. Scanning electron microscopy and rare earth element analyses demonstrate that the dated molluscs are diagenetically altered and contain aragonite cements that incorporated secondary young C, suggesting that such dates should be regarded as minimum ages. AMS ¹⁴C dated charcoals provide ages that occur out of stratigraphic order, and cluster in the upper chronological limits of the technique (～40–48 ka). Again, we suggest that such results should be regarded as suspicious and only minimum ages. Subsequent OSL and U/Th (teeth) dating provide complimentary results and demonstrate that the faunal sequences actually span ～120–83 ka, thus occurring beyond the AMS ¹⁴C dating window. Importantly, the dates suggest that the local decline in biological diversity was initiated ～75,000 years before the colonisation of humans on the continent. Collectively, the data are most parsimoniously consistent with a pre-human climate change model for local habitat change and megafauna extinction, but not with a nearly simultaneous extinction of megafauna as required by the human-induced blitzkrieg extinction hypothesis. This study demonstrates the problems inherent in dating deposits that lie near the chronological limits of the radiocarbon dating technique, and highlights the need to cross-check previously-dated archaeological and megafauna deposits within the timeframe of earliest human colonisation and latest megafaunal survival.     

Influence of a coarse interlayer on seawater intrusion and contaminant migration in coastal aquifers

Vertical 2D slice laboratory experiments were carried out in homogenous and layered sand tanks to elucidate the effects of a highly permeable (coarse‐grained sand) interlayer on seawater intrusion and transport of contaminants to a coastal sea. Tidal fluctuations produced oscillations in the seawater–freshwater transition zone, fluctuations of the contaminant infiltration rate and a zigzag contaminant plume outline. The seawater wedge became discontinuous at the (vertical) edges of the interlayer because of increased lateral movement of the seawater–freshwater interface within the interlayer. The contaminant plume formed a tail within the interlayer depending on the tidal stage, and similar to the wedge, its movement was accentuated. A simple analytical model that neglected vertical flow reliably predicted steady‐state seawater intrusion into the coastal aquifer. Numerical modeling was used to gain insight into the groundwater hydrodynamics and contaminant migration. The numerical results confirmed the experimental findings, i.e. that a highly permeable interlayer can provide a rapid transit path for contaminants to reach the seaward boundary and that the interlayer amplifies the effects of tidal fluctuations, resulting in wider transition zones for the seawater wedge and contaminant plume. Numerical simulations further showed that, with increasing interlayer hydraulic conductivity, the maximum seawater intrusion distance inside the interlayer increases approximately linearly. For the fixed‐head contaminant injection condition used, the model showed that contaminant infiltration increases approximately logarithmically with increasing interlayer hydraulic conductivity (other factors held fixed). Copyright © 2013 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.     

Piecewise prediction model for watershed‐scale erosion and sediment yield of individual rainfall events on the Loess Plateau,China

Establishing a universal watershed‐scale erosion and sediment yield prediction model represents a frontier field in erosion and soil/water conservation. The research presented here was conducted on the Chabagou watershed, which is located in the first sub‐region of the hill‐gully area of the Loess Plateau, China. A back‐propagation artificial neural model for watershed‐scale erosion and sediment yield was established, with the accuracy of the model, then compared with that of multiple linear regression. The sensitivity degree of various factors to erosion and sediment yield was quantitatively analysed using the default factor test. On the basis of the sensitive factors and the fractal information dimension, the piecewise prediction model for erosion and sediment yield of individual rainfall events was established and further verified. The results revealed the back‐propagation artificial neural network model to perform better than the multiple linear regression model in terms of predicting the erosion modulus, with the former able to effectively characterize dynamic changes in sediment yield under comprehensive factor conditions. The sensitivity of runoff erosion power and runoff depth to the erosion and sediment yield associated with individual rainfall events was found to be related to the complexity of surface topography. The characteristics of such a hydrological response are thus closely related to topography. When the fractal information dimension is greater than the topographic threshold, the accuracy of prediction using runoff erosion power is higher than that of using runoff depth. In contrast, when the fractal information dimension is smaller than the topographic threshold, the accuracy of prediction using runoff depth is higher than that of using runoff erosion power. The developed piecewise prediction model for watershed‐scale erosion and sediment yield of individual rainfall events, which introduces runoff erosion power and runoff depth using the fractal information dimension as a boundary, can be considered feasible and reliable and has a high prediction accuracy. Copyright © 2013 John Wiley & Sons, Ltd.     

Water level observations from unmanned aerial vehicles for improving estimates of surface water–groundwater interaction

Integrated hydrological models are usually calibrated against observations of river discharge and piezometric head in groundwater aquifers. Calibration of such models against spatially distributed observations of river water level can potentially improve their reliability and predictive skill. However, traditional river gauging stations are normally spaced too far apart to capture spatial patterns in the water surface, whereas spaceborne observations have limited spatial and temporal resolution. Unmanned aerial vehicles can retrieve river water level measurements, providing (a) high spatial resolution; (b) spatially continuous profiles along or across the water body, and (c) flexible timing of sampling. A semisynthetic study was conducted to analyse the value of the new unmanned aerial vehicle‐borne datatype for improving hydrological models, in particular estimates of groundwater–surface water (GW–SW) interaction. Mølleåen River (Denmark) and its catchment were simulated using an integrated hydrological model (MIKE 11–MIKE SHE). Calibration against distributed surface water levels using the Differential Evolution Adaptive Metropolis algorithm demonstrated a significant improvement in estimating spatial patterns and time series of GW–SW interaction. After water level calibration, the sharpness of the estimates of GW–SW time series improves by ~50% and root mean square error decreases by ~75% compared with those of a model calibrated against discharge only.     

Diatom‐based reconstruction of summer sea‐surface salinity in the South China Sea over the last 15 000 years

We present a new reconstruction of summer sea‐surface salinity (SSS) over the past 15 000 years based on a diatom record from piston core 17940, located on the northern slope of the South China Sea (SCS). The reconstructed diatom‐based summer SSS values for the modern period are in accord with instrumental observations of summer SSS in the area. Here, the modern summer SSS is primarily controlled by river runoff, in particular from the Pearl River. The reconstruction presented in this study shows that the summer SSS varied between 33.3 and 34.2 psu over the past 15 000 years. The long‐term summer SSS trend closely followed the trend of the orbitally controlled solar insolation at 20°N, suggesting that orbital forcing was the dominant driver of changes in summer SSS in this area. Comparisons to speleothem δ¹⁸O data and studies of surface hydrography in the region suggest that changes in solar insolation affected the summer SSS through changes in the East Asian Monsoon and sea‐level changes associated with the last deglaciation. Univariate spectral analyses indicate that centennial‐scale oscillatory variations in summer SSS were superimposed on the long‐term trend. During the deglacial period (c. 12 000–9000 cal. a BP), the dominant periodicity was centred around 230–250 years, whereas a ～350‐year oscillation dominated in the period 2200–4500 cal. a BP. The balance of evidence suggests that these centennial‐scale changes in summer SSS may have been driven by solar‐induced changes in the East Asian Monsoon, but further evidence is needed to firmly establish this relationship.     

Using hydrogeochemical signatures of stream water to assess pathways for rainfall events: towards a predictive model

Thanks to its simple division into agricultural and forestry land use, the Corbeira catchment (Galicia, Spain) is used as a case study to build a predictive model using hydrogeochemical signatures. Stream data acquired under recessional flow conditions over a one year period were obtained from a sampling station near the downstream end of the catchment, and using principal component analysis, it is shown that some of the analytical parameters are covariant, and some are negatively correlated. These findings support inferences about the pathways of rainfall in the catchment. Specific signatures may be associated with the dominant hydrological source, either surface runoff or subsurface waters: additionally, the dominant land use in that part of the catchment, where the flow originated, can also be predicted. The dominant runoff shows a strong covariance between suspended solids (SS) and particulate phosphorus (PP), with a clear negative correlation with pH. Dissolved organic carbon (DOC) data are associated with this covariant set when these compounds are available in the soils in question. Dissolved phosphorus, total organic nitrogen and dissolved nitrates are also associated with the same covariant set when the runoff flows through areas of extensive agricultural use. The SS ? PP covariance is less significant at lower flows. Typical base flow regimes show a significant covariance between salinity and pH, with a marked negative correlation with SS ? PP set, confirming the dominance of subsurface waters in the baseflow, as expected. Seasonally divergent DOC ? SS behaviour proves to be a useful tracer for rainfall regimes. The DOC trend shows a sinusoidal annual variation in amplitude, determined by the rainfall regime. As a result, flow from the catchment is dominated by surface water whenever there is synchronicity between the peaks of DOC and SS. Copyright © 2013 John Wiley & Sons, Ltd.     

Measurement of stream cross section using ground penetration radar with Hilbert–Huang transform

This study presents a new method to measure stream cross section without having contact with water. Compared with conventional measurement methods which apply instruments such as sounding weight, ground penetration radar (GPR), used in this study, is a non‐contact measurement method. This non‐contact measurement method can reduce the risk to hydrologists when they are conducting measurements, particularly in high flow period. However, the original signals obtained by using GPR are very complex, different from studies in the past where the measured data were mostly interpreted by experts with special skill or knowledge of GPR so that the results obtained were less objective. This study employs Hilbert–Huang transform (HHT) to process GPR signals which are difficult to interpret by hydrologists. HHT is a newly developed signal processing method that can not only process the nonlinear and non‐stationary complex signals, but also maintain the physical significance of the signal itself. Using GPR with HHT, this study establishes a non‐contact stream cross‐section measurement method with the ability to measure stream cross‐sectional areas precisely and quickly. Also, in comparison with the conventional method, no significant difference in results is found to exist between the two methods, but the new method can considerably reduce risk, measurement time, and manpower. It is proven that the non‐contact method combining GPR with HHT is applicable to quickly and accurately measure stream cross section. Copyright © 2013 John Wiley & Sons, Ltd.