Hydraulic Tomography: 3D Hydraulic Conductivity,Fracture Network,and Connectivity in Mudstone

We present the first demonstration of hydraulic tomography (HT) to estimate the three-dimensional (3D) hydraulic conductivity (K) distribution of a fractured aquifer at high-resolution field scale (HRFS), including the fracture network and connectivity through it. We invert drawdown data collected from packer-isolated borehole intervals during 42 pumping tests in a wellfield at the former Naval Air Warfare Center, West Trenton, New Jersey, in the Newark Basin. Five additional tests were reserved for a quality check of HT results. We used an equivalent porous medium forward model and geostatistical inversion to estimate 3D K at high resolution (K blocks <1 m³), using no strict assumptions about K variability or fracture statistics. The resulting 3D K estimate ranges from approximately 0.1 (highest-K fractures) to approximately 10⁻¹³ m/s (unfractured mudstone). Important estimated features include: (1) a highly fractured zone (HFZ) consisting of a sequence of high-K bedding-plane fractures; (2) a low-K zone that disrupts the HFZ; (3) several secondary fractures of limited extent; and (4) regions of very low-K rock matrix. The 3D K estimate explains complex drawdown behavior observed in the field. Drawdown tracing and particle tracking simulations reveal a 3D fracture network within the estimated K distribution, and connectivity routes through the network. Model fit is best in the shallower part of the wellfield, with high density of observations and tests. The capabilities of HT demonstrated for 3D fractured aquifer characterization at HRFS may support improved in situ remediation for contaminant source zones, and applications in mining, repository assessment, or geotechnical engineering.     

Assessing urban rainfall-runoff response to stormwater management extent

We report an empirical analysis of the hydrologic response of three small, highly impervious urban watersheds to pulse rainfall events, to assess how traditional stormwater management (SWM) alters urban hydrographs. The watersheds vary in SWM coverage from 3% to 61% and in impervious cover from 45% to 67%. By selecting a set of storm events that involved a single rainfall pulse with >96% of total precipitation delivered in 60 min, we reduced the effect of differences between storms on hydrograph response to isolate characteristic responses attributable to watershed properties. Watershed-average radar rainfall data were used to generate local storm hyetographs for each event in each watershed, thus compensating for the extreme spatial and temporal heterogeneity of short-duration, intense rainfall events. By normalizing discharge values to the discharge peak and centring each hydrograph on the time of peak we were able to visualize the envelope of hydrographs for each group and to generate representative composite hydrographs for comparison across the three watersheds. Despite dramatic differences in the fraction of watershed area draining to SWM features across these three headwater tributaries, we did not find strong evidence that SWM causes significant attenuation of the hydrograph peak. Hydrograph response for the three watersheds is remarkably uniform despite contrasts in SWM, impervious cover and spatial patterns of land cover type. The primary difference in hydrograph response is observed on the recession limb of the hydrograph, and that change appears to be associated with higher storm-total runoff in the watersheds with more area draining to SWM. Our findings contribute more evidence to the work of previous authors suggesting that SWM is less effective at attenuating urban hydrographs than is commonly assumed. Our findings also are consistent with previous work concluding that percent impervious cover may have greater influence on runoff volume than percent SWM coverage.     

Using high‐resolution isotope data and alternative calibration strategies for a tracer‐aided runoff model in a nested catchment

Testing hydrological models over different spatio‐temporal scales is important for both evaluating diagnostics and aiding process understanding. High‐frequency (6‐hr) stable isotope sampling of rainfall and runoff was undertaken during 3‐week periods in summer and winter within 12 months of daily sampling in a 3.2‐km² catchment in the Scottish Highlands. This was used to calibrate and test a tracer‐aided model to assess the (a) information content of high‐resolution data, (b) effect of different calibration strategies on simulations and inferred processes, and (c) model transferability to <1‐km² subcatchment. The 6‐hourly data were successfully incorporated without loss of model performance, improving the temporal resolution of the modelling, and making it more relevant to the time dynamics of the isotope and hydrometric response. However, this added little new information due to old‐water dominance and riparian mixing in this peatland catchment. Time variant results, from differential split sample testing, highlighted the importance of calibrating to a wide range of hydrological conditions. This also provided insights into the nonstationarity of catchment mixing processes, in relation to storage and water ages, which varied markedly depending on the calibration period. Application to the nested subcatchment produced equivalent parameterization and performance, highlighting similarity in dominant processes. The study highlighted the utility of high‐resolution data in combination with tracer‐aided models, applied at multiple spatial scales, as learning tools to enhance process understanding and evaluation of model behaviour across nonstationary conditions. This helps reveal more fully the catchment response in terms of the different mechanistic controls on both wave celerites and particle velocities.     

Alkaline phosphatase activity in the western English Channel: Elevations induced by high summertime rainfall

Alkaline phosphatase activity (APA) was determined in bulk particulate material and in a single-cell (ELF) assay at station L4 in the western English Channel during the summer of 2007. Throughout this period, the UK experienced its heaviest summertime rainfall since records began in 1914; with the result that riverine run-off into coastal waters was also elevated relative to long-term averages. Between May and August 2007, three distinct periods of elevated river run-off were observed which resulted in salinity minima at L4 on days 141, 190 and 232. An extended period of high river run-off between days 170 and 210 was responsible for decreases in near-surface salinity at L4 from 35.2068 to a minimum on day 190 of 34.7422. This contributed to the development of haline stratification which supported the development of an intense bloom of the centric diatom Chaetoceros debelis, with maximum observed chlorophyll a concentration of 8.69 μg l⁻¹. Minima in salinity, and maxima in chlorophyll concentration on day 190 were coincident with a peak in river-derived dissolved inorganic nitrogen (DIN) of 1.9 μmol l⁻¹ which was >5 times greater than the summertime mean and 24 times the concentrations experienced at L4 on weeks immediately before and after. There was no accompanying increase in dissolved inorganic phosphorus (DIP), and the DIN:DIP ratio increased to 49. With the inherent phosphorus stress that this caused, rates of APA increased from <4 to 42.4 nmolP l⁻¹ h⁻¹. ELF analysis on day 197 identified two taxa actively expressing alkaline phosphatase: the dinoflagellate Prorocentrum micans and ciliate Tiarana sp.     

Nitrogen dynamics within a wind-driven eddy

Wind-driven cyclonic eddies are hypothesized to relieve nutrient stress and enhance primary production by the upward displacement of nutrient-rich deep waters into the euphotic zone. In this study, we measured nitrate (NO₃⁻), particulate carbon (PC), particulate nitrogen (PN), their stable isotope compositions (δ¹⁵N-NO₃⁻, δ¹³C-PC and δ¹⁵N-PN, respectively), and dissolved organic nitrogen (DON) within Cyclone Opal, a mature wind-driven eddy generated in the lee of the Hawaiian Islands. Sampling occurred in March 2005 as part of the multi-disciplinary E-Flux study, approximately 4–6 weeks after eddy formation. Integrated NO₃⁻ concentrations above 110 m were 4.8 times greater inside the eddy (85.8±6.4 mmol N m⁻²) compared to the surrounding water column (17.8±7.8 mmol N m⁻²). Using N-isotope derived estimates of NO₃⁻ assimilation, we estimated that 213±59 mmol m⁻² of NO₃⁻ was initially injected into the upper 110 m Cyclone Opal formation, implying that NO₃⁻ was assimilated at a rate of 3.75±0.5 mmol N m⁻² d⁻¹. This injected NO₃⁻ supported 68±19% and 66±9% of the phytoplankton N demand and export production, respectively. N isotope data suggest that 32±6% of the initial NO₃⁻ remained unassimilated. Self-shading, inefficiency in the transfer of N from dissolved to particulate export, or depletion of a specific nutrient other than N may have led to a lack of complete NO₃⁻ assimilation. Using a salt budget approach, we estimate that dissolved organic nitrogen (DON) concentrations increased from eddy formation (3.8±0.4 mmol N m⁻²) to the time of sampling (4.0±0.09 mmol N m⁻²), implying that DON accumulated at rate of 0.83±1.3 mmol N m⁻² d⁻¹, and accounted for 22±15% of the injected NO₃⁻. Interestingly, no significant increase in suspended PN and PC, or export production was observed inside Cyclone Opal relative to the surrounding water column. A simple N budget shows that if 22±15% of the injected NO₃⁻ was shunted into the DON pool, and 32±6% is unassimilated, then 46±16% of the injected NO₃⁻ remains undocumented. Alternative loss processes within the eddy include lateral exchange of injected NO₃⁻ along isopycnal surfaces, remineralization of PN at depth, as well as microzooplankton grazing. A 9-day time series within Cyclone Opal revealed a temporal depletion in δ¹⁵N-PN, implying a rapid change in the N source. A change in NO₃⁻ assimilation, or a shift from NO₃⁻ fueled growth to assimilation of a ¹⁵N-deplete N source, may be responsible for such observations.     

Release and transformations of inorganic iodine by marine macroalgae

A number of field and laboratory studies on the impact of marine macroalgae on dissolved inorganic iodine speciation are presented. Within tidally isolated rock pools, the brown macroalga Fucus serratus was found to both release stored iodide and to facilitate the reduction of iodate to iodide. In contrast, no discernible changes in iodine speciation were observed in rock pools containing green macroalgae of the genus Ulva. Incubation experiments confirmed that the macroalgae Laminaria digitata, F. serratus and Kallymenia antarctica release iodide, though the rate of release varied between species and between specimens of the same species. Application of oxidative stress by treatment with cell wall derived oligoguluronate elicitors increased the efflux of iodide by L. digitata approximately 20-fold. The release of iodide by macroalgae may impact upon the formation of volatile iodine species (molecular iodine and iodocarbons) that are of importance in the coastal atmosphere.     

Short-term (≤2 yrs) estuarine mudflat and saltmarsh sedimentation: High-resolution data from ultrasonic altimetery, rod surface-elevation table, and filter traps

Rates of short-term (up to 2 years) bed elevation change and sedimentation from mudflats to salt marshes were measured in a rapidly infilling macrotidal estuary using an original combination of three high-resolution techniques: an ultrasonic altimeter, the Rod Surface-Elevation Table (RSET) method, and filter traps. The Authie estuary is located on a straight, sand-rich coast and is undergoing rapid infill under the influence of flood-dominant tides reinforced by wave action. The estuarine sediment suite consists of both mud and sand derived from the sea, of sand derived from storm wave erosion of dunes lining the north bank of the estuary, and, to a much smaller extent, of mud from the river catchment. Bed elevation change and sedimentation rates show an expected increase with the duration of tidal flooding (hydroperiod) in both space and time. The estuarine bed sediment suite changes from sandy at the mouth to muddy within the low-energy inner estuary, where mudflats are rapidly accreting, paving the way for the formation of increasingly denser and mature salt marshes from the high-sedimentation pioneer zone to the upper marsh where annual sedimentation is very low. Recorded variability in rates of bed elevation change and sedimentation reflect the influence of estuarine macro-scale and local sediment transport and depositional processes in a macrotidal context dominated by high inputs of allochthonous sediments.     

Long-term evolution (1988–2008) of Zostera spp. meadows in Arcachon Bay (Bay of Biscay)

The spatial variability of seagrass meadows in Arcachon Bay, was studied between 1988 and 2008 using a combination of mapping techniques based on aerial photographs for intertidal dwarf-grass (Zostera noltii) beds and acoustic sonar for permanently submerged eelgrass (Zostera marina) populations. The results show a severe decline over the period for both species, as well as an acceleration of the decline since 2005 for Z. noltii. The total surface regression over the studied period is estimated to be 22.8 km² for Z. noltii and 2.7 km² for Z. marina, which represent declines of 33 and 74% respectively.