Bayesian chemistry-assisted hydrograph separation (BACH) and nutrient load partitioning from monthly stream phosphorus and nitrogen concentrations

A Bayesian chemistry-assisted hydrograph separation (BACH) approach was developed, based on calibration of a three-component recursive digital filter, that requires monthly water quality data only. This enables BACH to be applied to the large number of rural catchments for which continuous flow records and monthly water chemistry time series exist from ‘state of environment’ monitoring programmes, but little supplementary data required for more sophisticated analysis techniques. As well as estimating fast, medium, and slow flow components, BACH also estimates for each flow component a time-invariant concentration of the chemical tracers chosen, allowing flow path-specific loads to be calculated. The method was demonstrated using 15 years of total phosphorus (TP) and total nitrogen (TN) data from eight mesoscale catchments in the Waikato region of New Zealand’s North Island. Calibration was done separately for three 5-year data periods, and validated against data from the following 5-year period. Flow path separation and concentration predictions were consistent between data periods, indicating that the TP–TN combination contained sufficient information to reliably identify three flow paths in each catchment; an event-response near-surface flow path with high concentrations of both phosphorus and nitrogen, a seasonal shallow groundwater flow path with lower concentrations of TP but high concentrations of TN, and a deeper slower groundwater flow path characterised by generally low concentrations of both TP and TN. Based on this analysis, the catchments were able to be grouped in three hydro-types. This shows that commonly available water quality data can support robust, objective flow separation and nutrient load apportionment, even in the absence of other supporting data, provided appropriate modelling methods are used.     

Fault slip and identification of the second fault plane in the Varzeghan earthquake doublet

An intraplate earthquake doublet, with 11-min delay between the events, devastated the city of Varzeghan in northwestern Iran on August 11, 2012. The first Mw 6.5 strike-slip earthquake, which occurred after more than 200 years of low seismicity, was followed by an Mw 6.4 oblique thrust event at an epicentral separation of about 6 km. While the first event can be associated with a distinct surface rupture, the absence of a surface fault trace and no clear aftershock signature makes it challenging to identify the fault plane of the second event. We use teleseismic body wave inversion to deduce the slip distribution in the first event. Using both P and SH waves stabilize the inversion and we further constrain the result with the surface rupture extent and the aftershock distribution. The obtained slip pattern shows two distinct slip patches with dissimilar slip directions where aftershocks avoid high-slip areas. Using the estimated slip for the first event, we calculate the induced Coulomb stress change on the nodal planes of the second event and find a preference for higher Coulomb stress on the N-S nodal plane. Assuming a simple slip model for the second event, we estimate the combined Coulomb stress changes from the two events on the focal planes of the largest aftershocks. We find that 90% of the aftershocks show increased Coulomb stress on one of their nodal planes when the N-S plane of the second event is assumed to be the correct fault plane.     

Data requirements for a 5-mm quasi-geoid in the Netherlands

We assess the surface gravity data requirements for a 5-mm quasi-geoid model for the Netherlands mainland and continental shelf in terms of omission and commission errors. The omission error critically depends on the roughness of the topography and bathymetry. For the Netherlands continental shelf, Central and Northern Netherlands, the omission error is well described by the model 0.32d mm, where d is the data spacing in km. For the more hilly Southern Netherlands, the omission error model is 0.92d mm. The commission error depends on the kernel modification, the data spacing, and the data accuracy. When using the spheroidal Stokes kernel, it is well described by 0.277 d σ_Δg mm, where σ_Δg is the noise standard deviation of surface gravity data in mGal. An upper bound of the commission error of the state-of-the-art satellite-only gravity model GOCO05S over the Netherlands is e^{0.03676L–11.419} mm, where L is the maximum degree up to which this model is used. Only if this model is truncated at a sufficiently low degree, e.g., at degree 100, its contribution to the total commission error can be neglected. We determine the total error as the sum of commission and omission errors. Hence, to realize a 5-mm quasi-geoid model for the Netherlands mainland and continental shelf, a data spacing of 3.5 km is needed when assuming a noise standard deviation of 1.5 mGal for surface gravity data. The currently available land-based gravity data fulfill this requirement. This does not apply to the situation at sea, where the density of the shipboard gravity data and the accuracy of the radar altimeter-derived data do not allow the realization of a 5-mm quasi-geoid model.     

Chemical Characterisation of NIST Silicate Glass Certified Reference Material SRM 610 by ICP-MS, TIMS, LIMS, SSMS, INAA, AAS and PIXE

National Institute of Science and Technology (NIST) silicate glass SRM 610 is widely used as a certified reference material for various micro-analytical techniques such as SIMS or laser ablation ICP-MS. SRM 610 has been nominally doped with sixty one trace elements at the 500 μg g⁻¹ level, but certified concentration data exist for only a few of these elements. This study reports concentration data for fifty nine trace elements obtained by ICP-MS, SSMS, LIMS, TIMS, INAA, AAS, and PIXE analyses of two different SRM 610 wafers. Most elements fall within a 10% band around a median value of about 440 μg g⁻¹. The REE concentrations are shown to be constant to 3% (1 σ), thus emphasizing the value of SRM 610 as a reference material for REE analyses.
Comparison of our values with published data suggests that different SRM 610 wafers are, within errors, chemically identical for most elements. Exceptions to this general rule appear to be restricted to elements which were partly lost during the production of the glass, e.g. Ag and Br. On the basis of six independent determinations of Rb concentrations, which are systematically lower by a few percent than the reported NIST value, we argue that the certified Rb concentration may not be representative for all distributed SRM 610 wafers.     

Groundwater—Surface waters interactions at slope and catchment scales: implications for landsliding in clay‐rich slopes

Understanding water infiltration and transfer in soft‐clay shales slopes is an important scientific issue, especially for landsliding. Geochemical investigations are carried out at the Super‐Sauze and Draix‐Laval landslides, both developed in the Callovo‐Oxfordian black marls, with the objective to define the origin of the groundwater. In situ investigations, soil leaching experiments and geochemical modeling are combined to identify the boundaries of the hydrological systems. At Super‐Sauze, the observations indicate that an external water flow occurs in the upper part of the landslide at the contact between the weathered black marls and the overlying formations, or at the landslide basement through a fault network. Such external origin of water is not observed at the local scale of the Draix‐Laval landslide but is detected at the catchment scale with the influence of deep waters in the streamwater quality of low river flows. Hydrogeological conceptual models are proposed emphasizing the role of the interactions between local (slope) and regional (catchment) flow systems. The observations suggest that this situation is a common case in the Alpine area. Expected consequences of the regional flows on slope stability are discussed in term of rise of pore water pressures and physicochemical weathering of the clay shales.     

Middle and Late Holocene Sea-Level Changes in Eastern Maine Reconstructed from Foraminiferal Saltmarsh Stratigraphy and AMS C Dates on Basal Peat

A relative sea-level history is reconstructed for Machiasport, Maine, spanning the past 6000 calendar years and combining two different methods. The first method establishes the long-term (10³ yr) trend of sea-level rise by dating the base of the Holocene saltmarsh peat overlying a Pleistocene substrate. The second method uses detailed analyses of the foraminiferal stratigraphy of two saltmarsh peat cores to quantify fluctuations superimposed on the long-term trend. The indicative meaning of the peat (the height at which the peat was deposited relative to mean tide level) is calculated by a transfer function based on vertical distributions of modern foraminiferal assemblages. The chronology is determined from AMS ¹⁴C dates on saltmarsh plant fragments embedded in the peat. The combination of the two different approaches produces a high-resolution, replicable sea-level record, which takes into account the autocompaction of the peat sequence. Long-term mean rates of sea-level rise, corrected for changes in tidal range, are 0.75 mm/yr between 6000 and 1500 cal yr B.P. and 0.43 mm/yr during the past 1500 years. The foraminiferal stratigraphy reveals several low-amplitude fluctuations during a relatively stable period between 1100 and 400 cal yr B.P., and a sea-level rise of 0.5 m during the past 300 years.     

Simulating compound flooding events in a hurricane

Ocean Dynamics - Compound flooding is usually induced by the concurrence of coastal storm surge and heavy precipitation induced river flooding, with the former involving oceanic processes and the...     

Density function evaluation from borehole gravity meter data – regularized spectral domain deconvolution approach

We present a new method of transforming borehole gravity meter data into vertical density logs. This new method is based on the regularized spectral domain deconvolution of density functions. It is a novel alternative to the “classical” approach, which is very sensitive to noise, especially for high‐definition surveys with relatively small sampling steps. The proposed approach responds well to vertical changes of density described by linear and polynomial functions. The model used is a vertical cylinder with large outer radius (flat circular plate) crossed by a synthetic vertical borehole profile. The task is formulated as a minimization problem, and the result is a low‐pass filter (controlled by a regularization parameter) in the spectral domain. This regularized approach is tested on synthetic datasets with noise and gives much more stable solutions than the classical approach based on the infinite Bouguer slab approximation. Next, the tests on real‐world datasets are presented. The properties and presented results make our proposed approach a viable alternative to the other processing methods of borehole gravity meter data based on horizontally layered formations.     

Geochemistry of a confirmed Precambrian impact ejecta deposit: The Grænsesø spherule layer,South Greenland

Samples from a single outcrop of the Graenseso spherule layer, Midternaes, South Greenland, consist of a spherule‐bearing dolomixtite with matrix‐supported intraclasts up to 1 m in size. In addition to field observations, we performed mineralogical and whole rock geochemical analysis, including electron microprobe, neutron activation analysis, X‐ray fluorescence, and mass spectrometry of the horizon and the overlying and underlying strata. We show that the spherules are petrographically similar to those in the Zaonega spherule layer, Karelia, Russia. Our petrographic and chemical results are consistent with the previous suggestion that the Grænsesø layer correlates with the Zaonega layer, and it is possible that both layers are related to the Vredefort impact event. The samples containing spherules, as well as the overlying rocks, have elevated REEs compared to the underlying pre‐impact layer, suggestive of a new continental source of sediment that may be coincident with the impact event. Zircons separated from the lower part of the Grænsesø spherule layer display complex age patterns suggesting that they have genetically different origins based on distinctly different Th/U ratios. Crystallization ages of all groups are ≥ 2.8 Ga, with ~2.8 Ga representing a time of major crustal growth globally. Therefore, we cannot conclusively determine in this study if the zircons are locally derived or if they are transported with the ejecta. The spherule layer was deposited by a high‐energy, subaqueous debris flow, an origin that is consistent with the hypothesis that the layer was deposited by impact‐induced waves and/or currents.