Evaluating Subsurface Parameterization to Simulate Hyporheic Exchange: The Steinlach River Test Site

Hyporheic exchange is the interaction of river water and groundwater, and is difficult to predict. One of the largest contributions to predictive uncertainty for hyporheic exchange has been attributed to the representation of heterogeneous subsurface properties. Our study evaluates the trade-offs between intrinsic (irreducible) and epistemic (reducible) model errors when choosing between homogeneous and highly complex subsurface parameter structures. We modeled the Steinlach River Test Site in Southwest Germany using a fully coupled surface water-groundwater model to simulate hyporheic exchange and to assess the predictive errors and uncertainties of transit time distributions. A highly parameterized model was built, treated as a “virtual reality” and used as a reference. We found that if the parameter structure is too simple, it will be limited by intrinsic model errors. By increasing subsurface complexity through the addition of zones or heterogeneity, we can begin to exchange intrinsic for epistemic errors. Thus, the appropriate level of detail to represent the subsurface depends on the acceptable range of intrinsic structural errors for the given modeling objectives and the available site data. We found that a zonated model is capable of reproducing the transit time distributions of a more detailed model, but only if the geological structures are known. An interpolated heterogeneous parameter field (cf. pilot points) showed the best trade-offs between the two errors, indicating fitness for practical applications. Parameter fields generated by multiple-point geostatistics (MPS) produce transit time distributions with the largest uncertainties, however, these are reducible by additional hydrogeological data, particularly flux measurements.     

An Open,Object-Based Framework for Generating Anisotropy in Sedimentary Subsurface Models

The spatial distribution of hydraulic properties in the subsurface controls groundwater flow and solute transport. However, many approaches to modeling these distributions do not produce geologically realistic results and/or do not model the anisotropy of hydraulic conductivity caused by bedding structures in sedimentary deposits. We have developed a flexible object-based package for simulating hydraulic properties in the subsurface—the Hydrogeological Virtual Realities (HyVR) simulation package. This implements a hierarchical modeling framework that takes into account geological rules about stratigraphic bounding surfaces and the geometry of specific sedimentary structures to generate realistic aquifer models, including full hydraulic-conductivity tensors. The HyVR simulation package can create outputs suitable for standard groundwater modeling tools (e.g., MODFLOW), is written in Python, an open-source programming language, and is openly available at an online repository. This paper presents an overview of the underlying modeling principles and computational methods, as well as an example simulation based on the Macrodispersion Experiment site in Columbus, Mississippi. Our simulation package can currently simulate porous media that mimic geological conceptual models in fluvial depositional environments, and that include fine-scale heterogeneity in distributed hydraulic parameter fields. The simulation results allow qualitative geological conceptual models to be converted into digital subsurface models that can be used in quantitative numerical flow-and-transport simulations, with the aim of improving our understanding of the influence of geological realism on groundwater flow and solute transport.     

Movement patterns of an endangered fishery species,Lithognathus lithognathus (Sparidae), and the role of no-take marine protected areas as a management tool

Understanding movement behaviour is essential for effective management of fishery species. Dart tags were used to study coastal movement patterns of white steenbras Lithognathus lithognathus, an overexploited seabream (Sparidae) in South Africa. In total, 6 962 fish (190–1 080 mm fork length) were tagged throughout the species’ distributional range, in four long-term fish-tagging programmes. The predominant behaviour recorded was residency, with relatively short-ranging movements. More than 60% of the 351 recaptured fish were recaptured within 1 km of their tagging site, some of which were at liberty in excess of three years. Most of the juveniles (93.9%), subadults (71.7%) and adults (64.0%) were recaptured within 10 km. Occasional long-distance movements of up to 800 km were recorded. Recapture distances were positively, but weakly, correlated with fish fork length (n = 257 fish measured at recapture; r² = 0.166, p < 0.001). Low levels of connectivity among coastal areas suggest that large-scale annual spawning migrations, as previously hypothesised for this species, are unlikely, which raises the possibility of multiple spawning sites. Seventy-seven percent of L. lithognathus tagged within three marine protected areas (MPAs) were recaptured within the same MPA, suggesting that area closures provide protection for L. lithognathus through the post-estuarine juvenile, subadult and adult life stages. We confirm that the country’s current network of coastal MPAs plays a vital role in sustaining this species, and suggest that additional closures, or otherwise substantial reductions in catch rates, are necessary for the species’ recovery.     

Philopatry and dispersal of juvenile leervis Lichia amia (Teleostei: Carangidae) tagged in a warm-temperate South African estuary

Understanding and characterising movement and area-use patterns of fishes within estuaries, as well as understanding the degree of connectivity between estuaries and the marine environment, can provide important insights into a species’ ecology, which is fundamental for effective management and conservation. Mark-and-recapture data obtained from the 082 TAG FISH Project, a dart-tagging programme, were used to describe movement patterns of juvenile leervis Lichia amia in the permanently open Swartkops Estuary, South Africa, as well as the degree of connectivity with the adjacent marine environment. A total of 628 juvenile leervis were tagged from 2008 to 2015, of which 77 fish (12.3%) were recaptured, with juveniles comprising 85.7% of all measured recaptures. The majority of recaptures (75.3%) were made within the Swartkops Estuary, indicating estuarine philopatry, with 23.4% having moved less than 1 km, revealing site fidelity, but most (51.9%) displaying estuarine roaming. A smaller portion of recaptures (27.7%) had dispersed to other habitats; 16 (20.8%) displayed multiple-habitat connectivity, while only 3 (3.9%) had undertaken long-distance movements (>100 km) and were recaptured up to 825 km northeast of the Swartkops Estuary. Lichia amia that had moved greater distances were characterised by a greater mean size and age at the time of recapture. This study provides valuable information on the philopatry and dispersal of juvenile L. amia.     

Estuarine and coastal connectivity of an estuarine-dependent fishery species,Pomadasys commersonnii (Haemulidae)

Understanding the level of connectivity between estuarine and coastal waters is essential for appropriate management of estuarine-associated taxa. Most studies have focused on the role of a single estuary, while limited research exists on the importance of multiple estuaries to individuals of estuarine-associated species. This study used acoustic telemetry to assess the usage of multiple estuaries and coastal waters by the estuarine-dependent spotted grunter Pomadasys commersonnii. Twenty-six adult fish were tagged with acoustic transmitters in the Kariega and Bushmans estuaries, South Africa, and their movements along a 300-km stretch of Indian Ocean coastline were monitored for up to 17 months. Tagged individuals spent most of their time in the estuary where they were tagged (55% and 85% for fish tagged in the two estuaries, respectively), followed by time in the sea (30% and 15%) and in other estuaries (15% and <1%). The mean durations of sea trips for fish tagged in the Kariega Estuary or Bushmans Estuary, respectively, were 25 days (range 3–55) and 12 days (range 2–22). Of the fish that went to sea, 93% from the Kariega Estuary and 60% from the Bushmans Estuary visited other estuaries. Most visits were undertaken to the Swartkops, Bushmans and Kowie estuaries, although the longest durations were spent in the Sundays Estuary. Individuals moved to estuaries up to 130 km away. The total distance travelled between estuaries by an individual during the study was ～529 km, with means of 201 and 184 km, respectively, for fish tagged in the Kariega and Bushmans estuaries. Despite covering large distances between estuaries, individuals often returned to their tagging estuary. Residency in their tagging estuary, combined with frequent visits to a neighbouring estuary, highlights the importance of estuarine habitats for this popular fish species, even after reaching maturity.     

The Pushtashan juvenile suprasubduction zone assemblage of Kurdistan (northeastern Iraq): A Cretaceous (Cenomanian) Neo-Tethys missing link

The Pushtashan suprasubduction zone assemblage of volcanic rocks, gabbros, norites and peridotites occurs in the Zagros suture zone, Kurdistan region, northeastern Iraq. Volcanic rocks are dominant in the assemblage and consist mainly of basalt and basaltic andesite flows with interlayered red shale and limestone horizons. Earlier lavas tend to be MORB-like, whereas later lavas display island arc tholeiite to boninitic geochemical characteristics. Tholeiitic gabbros intrude the norites and display fractionation trends typical of crystallisation under low-pressure conditions, whereas the norites display calc-alkaline traits, suggesting their source included mantle metasomatised by fluids released from subducted oceanic crust. Enrichment of Rb, Ba, Sr, Th and the presence of negative Nb anomalies indicate generation in a suprasubduction zone setting. Trondhjemite and granodiorite intrusions are present in the volcanic rocks, gabbros and norites. SHRIMP U-Pb dating of magmatic zircons from a granodiorite yields a mean~(206)Pb/~(238)U age of 96.0 ±2.0 Ma(Cenomanian). The initial ε_(Hf) value for the zircons show a narrow range from +12.8 to+15.6, with a weighted mean of + 13.90±0.96. This initial value is within error of model depleted mantle at 96 Ma or slightly below that, in the field of arc rocks with minimal contamination by older continental crust. The compositional bimodality of the Pushtashan suprasubduction sequence suggests seafloor spreading during the initiation of subduction, with a lava stratigraphy from earlyerupted MORB transitioning into calc-alkaline lavas and finally by 96 Ma intrusion of granodioritic and trondhjemitic bodies with juvenile crustal isotopic signatures. The results confirm another Cretaceous arc remnant preserved as an allochthon within the Iraqi segment of the Cenozoic Zagros suture zone. Implications for the closure of Neo-Tethys are discussed.     

Development and Testing of a Field Screening Method Based on Bubbling Extraction and Photoionization Detection for Measurement of Benzene and Total VOCs

A field screening method was developed for rapid measurement of benzene and gasoline range total petroleum hydrocarbons (TPHg) concentrations in groundwater. The method is based on collecting photoionization detector (PID) measurements from vapor samples. The vapor samples are collected by bubbling air through groundwater samples (air sparging) with a constant volume, temperature and sparging rate. The level of accuracy, sensitivity, precision, and statistical significance of the estimated concentrations, derived from the screening method, are comparable to conventional laboratory analytical results at concentrations equal to or greater than 150 µg/L for benzene and greater than 50 µg/L for TPHg. The method's concentration estimations can assist in making real‐time decisions regarding location of dissolved plumes and light nonaqueous phase liquid (LNAPL) source zones at many fuel release sites. The screening method was tested in the laboratory and in the field with 208 and 107 samples, respectively. The study concludes that the screening method can be used as a tool to aid in completing a site conceptual model as well as analyzing groundwater from monitoring wells.     

Iron(III) accumulations in inland saline waterways,Hunter Valley,Australia: Mineralogy,micromorphology and pore-water geochemistry

Discharge of Fe(II)-rich groundwaters into surface-waters results in the accumulation of Fe(III)-minerals in salinized sand-bed waterways of the Hunter Valley, Australia. The objective of this study was to characterise the mineralogy, micromorphology and pore-water geochemistry of these Fe(III) accumulations. Pore-waters had a circumneutral pH (6.2–7.2), were sub-oxic to oxic (Eh 59–453 mV), and had dissolved Fe(II) concentrations up to 81.6 mg L⁻¹. X-ray diffraction (XRD) on natural and acid-ammonium-oxalate (AAO) extracted samples indicated a dominance of 2-line ferrihydrite in most samples, with lesser amounts of goethite, lepidocrocite, quartz, and alumino-silicate clays. The majority of Fe in the samples was bound in the AAO extractable fraction (Fe^Ox) relative to the Na-dithionite extractable fraction (Fe^Di), with generally high Fe^Ox:Fe^Di ratios (0.52–0.92). The presence of nano-crystalline 2-line ferrihydrite (Fe₅HO₃·4H₂O) with lesser amounts of goethite (α-FeOOH) was confirmed by scanning electron microscopy (SEM) coupled with energy dispersive X-ray analysis (EDX), and transmission electron microscopy (TEM) coupled with selected area electron diffraction (SAED). In addition, it was found that lepidocrocite (γ-FeOOH), which occurred as nanoparticles as little as ∼5 lattice spacings thick perpendicular to the (0 2 0) lattice plane, was also present in the studied Fe(III) deposits. Overall, the results highlight the complex variability in the crystallinity and particle-size of Fe(III)-minerals which form via oxidation of Fe(II)-rich groundwaters in sand-bed streams. This variability may be attributed to: (1) divergent precipitation conditions influencing the Fe(II) oxidation rate and the associated supply and hydrolysis of the Fe(III) ion, (2) the effect of interfering compounds, and (3) the influence of bacteria, especially Leptothrix ochracea.     

Pore‐water pressure measurements: Mississippi delta submarine sediments

Abstract

A pore‐water pressure probe (piezometer) was implanted in Mississippi delta sediments at a preselected site (Block 28, South Pass area, 29°00´N, 89°15´W) 145 m from an offshore production platform (water depth approx. 19 m) in September 1975. Total pore‐water pressures (u_w ) were monitored for extended periods of time at depths of approximately 15 and 8 m below the mudline concurrently with hydrostatic pressures (u₈ ) measured at depths of 15 m and approximately 1 m below the mudline. Relatively high excess pore‐water pressures, u_e = (u_w ‐u₈ ), were recorded at the time of probe insertion measuring 99 kPa (14.4 psi) at 15 m and 50 kPa (7.3 psi) at 8 m. Six hours after the probe was implanted, excess pore pressures were still high at 81 kPa (11.8 psi, 15 m) and 37 kPa (5.4 psi, 8 m). Pore pressures appeared to become relatively constant at the 8‐m depth after 7 h had elapsed, and at the 15 m depth after 10–12 h. Excess pore‐water pressures averaged 72 kPa (10.4 psi, 15 m) and 32 kPa (4.6 psi, 8 m) prior to the initial effects of Hurricane Eloise, which passed in close proximity to the probe site. Significant variations in pressures were recorded during storm activity. As the effects of the storm subsided, excess pore‐water pressures began to decline slightly at the 15‐m depth; however, concurrently at the 8‐m depth, pore pressures began to increase gradually. During the period of 21–25 days after the probe was implanted, excess pore pressures appeared to become more constant, averaging 24 kPa (3.5 psi) at 15 m and 43 kPa (6.2 psi) at the 8‐m depth. The presence of methane, a common occurrence in these delta muds, may have influenced, or contributed to, the total pore‐water pressures measured during this experiment.     

Initial results and progress of the Mississippi delta sediment pore water pressure experiment

Abstract

This report describes the instrumentation, initial results, and progress of an experiment designed to measure and monitor submarine sediment pore water and hydrostatic pressures in a selected area of the Mississippi Delta. The experiment also is intended to monitor significant pressure perturbations during active storm periods. Initial analysis of the data revealed excess pore water pressures in the silty clay sediment at selected depths below the mudline. Continuous monitoring of the pore water and hydrostatic pressures was expected to reveal important information regarding sediment pore water pressure variations as a function of the geological processes active in the Mississippi Delta.