Establishment of background water quality conditions in the Great Zab River catchment: influence of geogenic and anthropogenic controls on developing a baseline for water quality assessment and resource management

The Great Zab River catchment is a major left-bank tributary of the River Tigris and drains a substantial part of the Kurdistan Region, an autonomous region of Northern Iraq. Within Kurdistan, the water resources of the Great Zab River catchment are under pressure from population increase and are utilized for potable, domestic and agricultural and industrial supply. As with many parts of the world, effective management of water resources within Kurdistan is hindered by a lack of water quality data and established background concentrations. This study therefore represents the first regional survey of river water chemistry for the Great Zab River catchment and presents data on the spatial and temporal trends in concentrations of As, Ba, Ca, Cd, Co, Cr, Cu, Fe, Hg, Li, Mn, Mo, Ni, Pb, Sr, Zn, NO₃^?, SO₄^2?, F^?, Cl^? and PO₄^3?, in addition to pH, electrical conductivity, dissolved oxygen and turbidity. As a tool for underpinning the management and monitoring of water quality, background concentrations were defined for the Great Zab catchment using three methods. The influences of geogenic and anthropogenic controls upon spatial and temporal trends in water chemistry are also evaluated. The influence of geogenic loading from underlying bedrock was identifiable within the observed spatial trends, with the most notable differences found between waters sampled from the relatively more volcanic-rich Zagros zone to the north and those sampled from the lower catchment underlain by younger clay-, sand- and siltstones. The greatest anthropogenic influence, identifiable through elements such as Cl^? and NO₃^?, is present in the more highly populated lower catchment. The background concentrations identified in the Great Zab catchment would be those expected as a result of geogenic loading with some anthropogenic influence and represent a more conservative value when compared to those such as the World Health Organization Maximum Admissible Concentration. However, background concentrations represent a powerful tool for identifying potential anthropogenic impacts on water quality and informing management of such occurrences.     

40Ar/39Ar ages and zircon petrochronology for the rear arc of the Izu-Bonin-Marianas intra-oceanic subduction zone

Long-lived intra-oceanic arcs of Izu-Bonin-Marianas (IBM)-type are built on thick, granodioritic crust formed in the absence of pre-existing continental crust. International Ocean Discovery Program Expedition 350, Site U1437, explored the IBM rear arc to better understand continental crust formation in arcs. Detailed petrochronological (U–Pb geochronology combined with trace elements, oxygen and hafnium isotopes) characterizations of zircon from Site U1437 were carried out, taking care to exclude potential contaminants by (1) comparison of zircon ages with ship-board palaeomagnetic and biostratigraphic ages and ⁴⁰Ar/³⁹Ar geochronology, (2) analysing zircon from drill muds for comparison, (3) selectively carrying out in situ analysis in petrographic thin sections, and (4) minimizing potential laboratory contamination through using pristine equipment during mineral separation. The youngest zircon ages in Site U1437 are consistent with ⁴⁰Ar/³⁹Ar and shipboard ages to a depth of ~1390 m below sea floor (mbsf) where Igneous Unit Ig 1 yielded an ⁴⁰Ar/³⁹Ar age of 12.9 ± 0.3 Ma (all errors 2σ). One single zircon (age 15.4 ± 1.0 Ma) was recovered from the deepest lithostratigraphic unit drilled, Unit VII (1459.80–1806.5 mbsf). Site U1437 zircon trace element compositions are distinct from those of oceanic and continental arc environments and differ from those generated in thick oceanic crust (Iceland-type) where low-δ¹⁸O evolved melts are produced via re-melting of hydrothermally altered mafic rocks. Ti-in-zircon model temperatures are lower than for mid-ocean ridge rocks, in agreement with low zircon saturation temperatures, suggestive of low-temperature, hydrous melt sources. Zircon oxygen (δ¹⁸O = 3.3–6.0‰) and hafnium (εHf = + 10–+16) isotopic compositions indicate asthenospheric mantle sources. Trace element and isotopic differences between zircon from Site U1437 rear-arc rocks and the Hadean detrital zircon population suggest that preserved Hadean zircon crystals were probably generated in an environment different from modern oceanic convergent margins underlain by depleted mantle.     

The Origin and Evolution of the Kaapvaal Cratonic Lithospheric Mantle

A detailed petrological and geochemical study of low-temperatureperidotite xenoliths from Kimberley and northern Lesotho ispresented to constrain the processes that led to the magmaphileelement depletion of the Kaapvaal cratonic lithospheric mantleand its subsequent re-enrichment in Si and incompatible traceelements. Whole-rocks and minerals have been characterized forRe–Os isotope compositions, and major and trace elementconcentrations, and garnet and clinopyroxene for Lu–Hfand Sm–Nd isotope compositions. Most samples are characterizedby Archaean Os model ages, low Al, Fe and Ca contents, highMg/Fe, low Re/Os, very low (< 0·1 x chondrite) heavyrare earth element (HREE) concentrations and a decoupling betweenNd and Hf isotope ratios. These features are most consistentwith initial melting at 3·2 Ga followed by metasomatismby hydrous fluids, which may have also caused additional meltingto produce a harzburgitic residue. The low HREE abundances ofthe peridotites require that extensive melting occurred in thespinel stability field, possibly preceded by some melting inthe presence of garnet. Fractional melting models suggest that30% melting in the spinel field or 20% melting in the garnetfield followed by 20% spinel-facies melting are required toexplain the most melt-depleted samples. Garnet Nd–Hf isotopecharacteristics indicate metasomatic trace element enrichmentduring the Archaean. We therefore suggest a model includingshallow ridge melting, followed by metasomatism of the Kaapvaalupper mantle in subduction zones surrounding cratonic nuclei,probably during amalgamation of smaller pre-existing terranesin the Late Archaean (2·9 Ga). The fluid-metasomatizedresidua have subsequently undergone localized silicate meltinfiltration that led to clinopyroxene ± garnet enrichment.Calculated equilibrium liquids for clinopyroxene and their Hf–Ndisotope compositions suggest that most diopside in the xenolithscrystallized from an infiltrating kimberlite-like melt, eitherduring Group II kimberlite magmatism at 200–110 Ma (Kimberley),or shortly prior to eruption of the host kimberlite around 90Ma (northern Lesotho). KEY WORDS: Kaapvaal craton; lithospheric mantle; metasomatism; Nd–Hf isotopes; Re–Os isotopes     

Response spectral relationships for rock sites derived from the component attenuation model

The seismological model was developed initially from the fundamental relationship between earthquake ground motion properties and the seismic moment generated at the source of the earthquake. Following two decades of continuous seismological research in the United States, seismological models which realistically account for both the source and path effects on the seismic shear waves have been developed and their accuracy rigorously verified (particularly in the long and medium period ranges). An important finding from the seismological research by Atkinson and Boore and their co‐investigators is the similarity of the average frequency characteristics of seismic waves generated at the source between the seemingly very different seismic environments of Eastern and Western North America (ENA and WNA, respectively). A generic definition of the average source properties of earthquakes has therefore been postulated, referred to herein as the generic source model. Further, the generic ‘hard rock’ crustal model which is characteristic of ENA and the generic ‘rock’ crustal model characteristic of WNA have been developed to combine with the generic source model, hence enabling simulations to be made of the important path‐related modifications to ground motions arising from different types of crustal rock materials. It has been found that the anelastic contribution to whole path attenuation is consistent between the ENA and WNA models, for earthquake ground motions (response spectral velocities and displacements) in the near and medium fields, indicating that differences in the ENA and WNA motions arise principally from the other forms of path‐related modifications, namely the mid‐crust amplification and the combined effect of the upper‐crust amplification and attenuation, both of which are significant only for the generic WNA ‘rock’ earthquake ground motions. This paper aims to demonstrate the effective utilization of the latest seismological model, comprising the generic source and crustal models, to develop a response spectral attenuation model for direct engineering applications. The developed attenuation model also comprises a source factor and several crustal (wave‐path modification) component factors, and thus has also been termed herein the component attenuation model (CAM). Generic attenuation relationships in CAM, which embrace both ENA and WNA conditions, have been developed using stochastic simulations. The crustal classification of a region outside North America can be based upon regional seismological and geological information. CAM is particularly useful for areas where local strong motion data are lacking for satisfactory empirical modelling. In the companion paper entitled ‘response spectrum modelling for rock sites in low and moderate seismicity regions combining velocity, displacement and acceleration predictions’, the CAM procedure has been incorporated into a response spectrum model which can be used to effectively define the seismic hazard of bedrock sites in low and moderate seismicity regions. This paper and the companion paper constitute the basis of a long‐term objective of the authors, to develop and effectively utilize the seismological model for engineering applications worldwide.     

A technique for the reduction of dynamic degrees of freedom

This paper presents a variation on the component mode technique for the dynamic substructuring of large‐scale structural analysis of building and bridge frames. The principal innovation of the proposed method of dynamic reduction is that the resulting mass matrix of the reduced substructures remains diagonal. As in the component mode technique, the reduction is accomplished by transforming the degrees of freedom in the substructure using boundary shapes and internal shapes. The diagonal mass matrix is achieved by orthogonalization of the boundary shapes to the internal shapes, and a selective row‐by‐row summation of the mass matrix into the diagonal entry (where off‐diagonal terms are unavoidable). To aid in recovering the accuracy of the rigid‐body inertias that is lost in the diagonalization process, additional pseudo‐rigid‐body‐mode shapes are proposed. Copyright © 2001 John Wiley & Sons, Ltd.     

The insoluble residue test for abundance of carbonate

Among methods for determining the abundance of carbonate minerals in solid samples, dissolution in acid is appealing because it is cheap and can be quick, but it suffers from several sources of uncertainty. The weight loss due to evaporation of water can be corrected for easily by measuring the weight loss from blank vessels as well as sample vessels. However, a further correction, for excess evaporation due to heat liberated by dissolution, is desirable; we show how to calculate it. Certain puzzling results point to a need for further investigation of interfering factors. © 1997 John Wiley & Sons, Ltd.     

Uplift rates defined by U-series and 14C ages of serpulid-encrusted speleothems from submerged caves near Siracusa,Sicily (Italy)

We have established a plausible rate of uplift near Siracusa in southeastern Sicily (Italy) over the last glacial–interglacial cycle using U-series ages of submerged speleothem calcite and ¹⁴C ages of calcite serpulid layers that encrust the speleothems during cave submergence. The precisely determined ages of these sea level benchmarks were compared with expected relative sea level position based on glacio-hydro-isostatic modeling to assess the rate of uplift in this region. When combined with the age of various late Holocene archaeological sites that have been recently described and characterized in terms of their functional position relative to sea level these data collectively define a rate of uplift ≤0.4 mm a^?1 along this portion of the Sicilian coastline. These results are consistent with an age assignment of marine isotope stage (MIS) 5.3 or 5.5 for the Akradina terrace, which in turn places temporal constraints on paleoshorelines above and below this level.     

Channel response to sediment release: insights from a paired analysis of dam removal

Dam removals with unmanaged sediment releases are good opportunities to learn about channel response to abruptly increased bed material supply. Understanding these events is important because they affect aquatic habitats and human uses of floodplains. A longstanding paradigm in geomorphology holds that response rates to landscape disturbance exponentially decay through time. However, a previous study of the Merrimack Village Dam (MVD) removal on the Souhegan River in New Hampshire, USA, showed that an exponential function poorly described the early geomorphic response. Erosion of impounded sediments there was two‐phased. We had an opportunity to quantitatively test the two‐phase response model proposed for MVD by extending the record there and comparing it with data from the Simkins Dam removal on the Patapsco River in Maryland, USA. The watershed sizes are the same order of magnitude (10² km²), and at both sites low‐head dams were removed (~3–4 m) and ~65 000 m³ of sand‐sized sediments were discharged to low‐gradient reaches. Analyzing four years of repeat morphometry and sediment surveys at the Simkins site, as well as continuous discharge and turbidity data, we observed the two‐phase erosion response described for MVD. In the early phase, approximately 50% of the impounded sediment at Simkins was eroded rapidly during modest flows. After incision to base level and widening, a second phase began when further erosion depended on floods large enough to go over bank and access impounded sediments more distant from the newly‐formed channel. Fitting functional forms to the data for both sites, we found that two‐phase exponential models with changing decay constants fit the erosion data better than single‐phase models. Valley width influences the two‐phase erosion responses upstream, but downstream responses appear more closely related to local gradient, sediment re‐supply from the upstream impoundments, and base flows. Copyright © 2017 John Wiley & Sons, Ltd.     

New statistical methods for the comparison and characterization of particle shape

This paper presents novel methods for robust statistical testing of particle shape data. Shape (the relative lengths of three orthogonal axes) is a key property of sedimentary particles, providing information on provenance, transport history and depositional environment. However, the usefulness of shape data, including the ability to make robust comparisons between samples, has been constrained by the absence of a satisfactory definition of the mean shape for a sample of particles. Such a definition is proposed and used to develop confidence regions for the population mean shape using both parametric (theoretical) and computational (bootstrap) methods. These techniques are based on a transform that permits multivariate statistical methods for the analysis of compositional data to be extended to shape. These techniques are validated with reference to a dataset of 169 clast samples and found to perform well. A statistical test on the mean – using the multivariate extension of Student's t-test, Hotelling's T² – is presented. The benefits of the methods presented are demonstrated with reference to a case study. © 2019 John Wiley & Sons, Ltd.