Iron and Manganese in Groundwater: Using Kriging and GIS to Locate High Concentrations in Buncombe County,North Carolina

For health, economic, and aesthetic reasons, allowable concentrations (as suggested by the United States Environmental Protection Agency) of the secondary contaminants iron (Fe) and manganese (Mn) found present in drinking water are 0.3 and 0.05 mg/L, respectively. Water samples taken from private drinking wells in rural communities within Buncombe County, North Carolina contain concentrations of these metals that exceed secondary water quality criteria. This study predicted the spatial distribution of Fe and Mn in the county, and evaluated the effect of site environmental factors (bedrock geology, ground elevation, saprolite thickness, and drinking water well depth) in controlling the variability of Fe and Mn in groundwater. A statistically significant correlation between Fe and Mn concentrations, attributable to bedrock geology, was identified. Prediction models were created using ordinary kriging and cokriging interpolation techniques to estimate the presence of Fe and Mn in groundwater where direct measurements are not possible. This same procedure can be used to estimate the trend of other contaminants in the groundwater in different areas with similar hydrogeological settings.     

A regime shift in sediment export from a coastal watershed during a record wet winter,California: Implications for landscape response to hydroclimatic extremes

Small, steep watersheds are prolific sediment sources from which sediment flux is highly sensitive to climatic changes. Storm intensity and frequency are widely expected to increase during the 21st century, and so assessing the response of small, steep watersheds to extreme rainfall is essential to understanding landscape response to climate change. During record winter rainfall in 2016–2017, the San Lorenzo River, coastal California, had nine flow peaks representing 2–10‐year flood magnitudes. By the third flood, fluvial suspended sediment showed a regime shift to greater and coarser sediment supply, coincident with numerous landslides in the watershed. Even with no singular catastrophic flood, these flows exported more than half as much sediment as had a 100‐year flood 35 years earlier, substantially enlarging the nearshore delta. Annual sediment load in 2017 was an order of magnitude greater than during an average‐rainfall year, and 500‐fold greater than in a recent drought. These anomalous sediment inputs are critical to the coastal littoral system, delivering enough sediment, sometimes over only a few days, to maintain beaches for several years. Future projections of megadroughts punctuated by major atmospheric‐river storm activity suggest that interannual sediment‐yield variations will become more extreme than today in the western USA, with potential consequences for coastal management, ecosystems, and water‐storage capacity. The occurrence of two years with major sediment export over the past 35 years that were not associated with extremes of the El Niño Southern Oscillation or Pacific Decadal Oscillation suggests caution in interpreting climatic signals from marine sedimentary deposits derived from small, steep, coastal watersheds, to avoid misinterpreting the frequencies of those cycles. Published 2018. This article is a U.S. Government work and is in the public domain in the USA.     

On the frontal condition for finite amplitude α2Ω-dynamo wave trains in stellar shells

Abstract

In a previous paper, Bassom et al. (Proc. R. Soc. Lond. A, 455, 1443–1481, 1999) (BKS) investigated finite amplitude αΩ-dynamo wave trains in a thin turbulent, differentially rotating convective stellar shell; nonlinearity arose from α-quenching. There asymptotic solutions were developed based upon the small aspect ratio ε of the shell. Specifically, as a consequence of a prescribed latitudinally dependent α-effect and zonal shear flow, the wave trains have smooth amplitude modulation but are terminated abruptly across a front at some high latitude θ_F. Generally, the linear WKB-solution ahead of the front is characterised by the vanishing of the complex group velocity at a nearby point θ_f; this is essentially the Dee–Langer criterion, which determines both the wave frequency and front location.

Recently, Griffiths et al. (Geophys. Astrophys. Fluid Dynam. 94, 85–133, 2001) (GBSK) obtained solutions to the α²Ω-extension of the model by application of the Dee—Langer criterion. Its justification depends on the linear solution in a narrow layer ahead of the front on the short O(θ_f—θ_F) length scale; here conventional WKB-theory, used to describe the solution elsewhere, is inadequate because of mode coalescence. This becomes a highly sensitive issue, when considering the transition from the linear solution, which occurs when the dynamo number D takes its critical value D _c corresponding to the onset of kinematic dynamo action, to the fully nonlinear solutions, for which the Dee—Langer criterion pertains.

In this paper we investigate the nature of the narrow layer for α²Ω-dynamos in the limit of relatively small but finite α-effect Reynolds numbers R _α, explicitly ε^½ ? R ² _α ? 1. Though there is a multiplicity of solutions, our results show that the space occupied by the corresponding wave train is generally maximised by a solution with θ_f—θ_F small; such solutions are preferred as evinced by numerical simulations. This feature justifies the application by GBSK of the Dee—Langer criterion for all D down to the minimum D _min that the condition admits. Significantly, the frontal solutions are subcritical in the sense that |D _min| ≤ |D _c|; equality occurs as the α-effect Reynolds number tends to zero. We demonstrate that the critical linear solution is not connected by any parameter track to the preferred nonlinear solution associated with D _min. By implication, a complicated bifurcation sequence is required to make the connection between the linear and nonlinear states. This feature is in stark contrast to the corresponding results for αΩ-dynamos obtained by BKS valid in the limit R _{2 α} ? ε^½, which, though exhibiting a weak subcriticality, showed that the connection follows a clearly identifiable nonbifurcating track.     

High‐frequency water quality monitoring in an urban catchment: hydrochemical dynamics,primary production and implications for the Water Framework Directive

This paper describes the hydrochemistry of a lowland, urbanised river‐system, The Cut in England, using in situ sub‐daily sampling. The Cut receives effluent discharges from four major sewage treatment works serving around 190 000 people. These discharges consist largely of treated water, originally abstracted from the River Thames and returned via the water supply network, substantially increasing the natural flow. The hourly water quality data were supplemented by weekly manual sampling with laboratory analysis to check the hourly data and measure further determinands. Mean phosphorus and nitrate concentrations were very high, breaching standards set by EU legislation. Although 56% of the catchment area is agricultural, the hydrochemical dynamics were significantly impacted by effluent discharges which accounted for approximately 50% of the annual P catchment input loads and, on average, 59% of river flow at the monitoring point. Diurnal dissolved oxygen data demonstrated high in‐stream productivity. From a comparison of high frequency and conventional monitoring data, it is inferred that much of the primary production was dominated by benthic algae, largely diatoms. Despite the high productivity and nutrient concentrations, the river water did not become anoxic, and major phytoplankton blooms were not observed. The strong diurnal and annual variation observed showed that assessments of water quality made under the Water Framework Directive (WFD) are sensitive to the time and season of sampling. It is recommended that specific sampling time windows be specified for each determinand, and that WFD targets should be applied in combination to help identify periods of greatest ecological risk. © 2015 The Authors. Hydrological Processes published by John Wiley & Sons Ltd.     

Particle‐size evidence of barrier estuary regime as a new proxy for ENSO climate variability

We investigate a new proxy for ENSO climate variability based on particle‐size data from long‐term, coastal sediment records preserved in a barrier estuary setting. Corresponding ~4–8 year periodicities identified from Wavelet analysis of particle‐size data from Pescadero Marsh in Central Coast California and rainfall data from San Francisco reflect established ENSO periodicity, as further evidenced in the Multivariate ENSO Index (MEI), and thus confirms an important ENSO control on both precipitation and barrier regime variability. Despite the fact that barrier estuary mean particle size is influenced by coastal erosion, precipitation and streamflow, balanced against barrier morphology and volume, it is encouraging that considerable correspondence can also be observed in the time series of MEI, regional rainfall and site‐based mean particle size over the period 1871–2008. This correspondence is, however, weakened after c.1970 by temporal variation in sedimentation rate and event‐based deposition. These confounding effects are more likely when: (i) accommodation space may be a limiting factor; and (ii) particularly strong El Niños, e.g. 1982/1983 and 1997/1998, deposit discrete >cm‐thick units during winter storms. The efficacy of the sediment record of climate variability appears not to be compromised by location within the back‐barrier setting, but it is limited to those El Niños that lead to barrier breakdown. For wider application of this particle size index of ENSO variability, it is important to establish a well‐resolved chronology and to sample the record at the appropriate interval to characterize deposition at a sub‐annual scale. Further, the sample site must be selected to limit the influence of decreasing accommodation space through time (infilling) and event‐based deposition. It is concluded that particle‐size data from back‐barrier sediment records have proven potential for preserving evidence of sub‐decadal climate variability, allowing researchers to explore temporal and spatial patterns in phenomena such as ENSO. Copyright © 2017 John Wiley & Sons, Ltd.     

Channel response to an extreme flood and sediment pulse in a mixed bedrock and gravel‐bed river

We exploit a natural experiment caused by an extreme flood (~500 year recurrence interval) and sediment pulse derived from more than 2500 concurrent landslides to explore the influence of valley‐scale geomorphic controls on sediment slug evolution and the impact of sediment pulse passage and slug deposition and dispersion on channel stability and channel form. Sediment slug movement is a crucial process that shapes gravel‐bed rivers and alluvial valleys and is an important mechanism of downstream bed material transport. Further, increased bed material transport rates during slug deposition can trigger channel responses including increases in lateral mobility, channel width, and alluvial bar dominance. Pre‐ and post‐flood LiDAR and aerial photographs bracketing the 2007 flood on the Chehalis River in south‐western Washington State, USA, document the channel response with high spatial and temporal definition. The sediment slug behaved as a Gilbert Wave, with both channel aggradation and sequestration of large volumes of material in floodplains of headwaters' reaches and reaches where confined valleys enter into broad alluvial valleys. Differences between the valley form of two separate sub‐basins impacted by the pulse highlight the important role channel and channel‐floodplain connectivity play in governing downstream movement of sediment slug material. Finally, channel response to the extreme flood and sediment pulse illustrate the connection between bed material transport and channel form. Specifically, the channel widened, lateral channel mobility increased, and the proportion of the active channel covered by bars increased in all reaches in the study area. The response scaled tightly with the relative amount of bed material sediment transport through individual reaches, indicating that the amount of morphological change caused by the flood was conditioned by the simultaneous introduction of a sediment pulse to the channel network. Copyright © 2015 John Wiley & Sons, Ltd.     

The next few years of the ICES Journal of Marine Science

This issue marks a change in the editorial team of the ICESJournal of Marine Science, in that we have bid farewell to ourformer Editor-in-Chief, Niels Daan, after six years of sterlingservice, and to our Elsevier Publishing Editor, Andrew Richford,after an even longer period of years. Both deserve credit fortheir energy and vision that has left us the legacy of a journalas widely read and cited as the ICES Journal now is, but theyleave a gap that will be a huge challenge for us to meet. However,in my new capacity as Editor-in-Chief and in Els Bosma's capacityas Elsevier's Publishing Editor, I     

Structural controls on the formation of BSR over a diapiric anticline from a dense MCS survey offshore southwestern Taiwan

A dense seismic reflection survey with up to 250-m line-spacing has been conducted in a 15 × 15 km wide area offshore southwestern Taiwan where Bottom Simulating Reflector is highly concentrated and geochemical signals for the presence of gas hydrate are strong. A complex interplay between north–south trending thrust faults and northwest–southeast oblique ramps exists in this region, leading to the formation of 3 plunging anticlines arranged in a relay pattern. Landward in the slope basin, a north–south trending diapiric fold, accompanied by bright reflections and numerous diffractions on the seismic profiles, extends across the entire survey area. This fold is bounded to the west by a minor east-verging back-thrust and assumes a symmetric shape, except at the northern and southern edges of this area, where it actively overrides the anticlines along a west-verging thrust, forming a duplex structure. A clear BSR is observed along 67% of the acquired profiles. The BSR is almost continuous in the slope basin but poorly imaged near the crest of the anticlines. Local geothermal gradient values estimated from BSR sub-bottom depths are low along the western limb and crest of the anticlines ranging from 40 to 50 °C/km, increase toward 50–60 °C/km in the slope basin and 55–65 °C/km along the diapiric fold, and reach maximum values of 70 °C/km at the southern tip of the Good Weather Ridge. Furthermore, the local dips of BSR and sedimentary strata that crosscut the BSR at intersections of any 2 seismic profiles have been computed. The stratigraphic dips indicated a dominant east–west shortening in the study area, but strata near the crest of the plunging anticlines generally strike to southwest almost perpendicular to the direction of plate convergence. The intensity of the estimated bedding-guided fluid and gas flux into the hydrate stability zone is weaker than 2 in the slope basin and the south-central half of the diapiric fold, increases to 7 in the northern half of the diapiric fold and plunging anticlines, and reaches a maximum of 16 at the western frontal thrust system. Rapid sedimentation, active tectonics and fluid migration paths with significant dissolved gas content impact on the mechanism for BSR formation and gas hydrate accumulation. As we begin to integrate the results from these studies, we are able to outline the regional variations, and discuss the importance of structural controls in the mechanisms leading to the gas hydrate emplacements.