Nature of sediment dispersal off the Sepik River, Papua New Guinea: preliminary sediment budget and implications for margin processes

A sediment budget is constructed for the slope and narrow continental shelf off the Sepik River in order to estimate the relative importance of turbid plumes versus bottom gravity transport through a near-shore submarine canyon in the dispersal of sediment across this collision margin. ²¹⁰Pb geochronology and inventories of Kasten cores are consistent with the northwestward dispersal of sediment from the river mouth via hypopycnal and possible isopycnal plumes. Sediment accumulation rates are 5 cm yr⁻¹ on the upper slope just off of the Sepik mouth, decreasing gradually to 1 cm yr⁻¹ toward the northwest, and decreasing abruptly offshore (<0.2 cm yr⁻¹ at 1200 m water depth). A sediment budget indicates that only about 7–15% of the Sepik River sediment discharge accumulates on the adjacent open shelf and slope. The remainder presumably escapes offshore via gravity flows through a submarine canyon, the head of which extends into the river mouth. The divergent sediment pathways observed off the Sepik River (i.e., surface and subsurface plumes versus sediment gravity flows through a canyon) may be common along high-yield collision margins of the Indo–Pacific archipelago, and perhaps are analogous to most margins during Late Quaternary low sea-level conditions.     

Organic biomarkers for tracing carbon cycling in the Gulf of Papua (Papua New Guinea)

Sediment traps were deployed in the Gulf of Papua in June–July 1997, to determine fluxes of organic matter and inorganic elements from the photic zone to deeper waters at the base of the continental slope and in the northern Coral Sea. Three stations, ranging from 900 to 1500 m depth, had “shallow” traps at 300 m below the water surface and “deep” traps set 100 m above the bottom. Infiltrex II water samplers collected particulate and dissolved organic matter from the Fly, Purari and Kikori rivers, and near-surface water from the shelf of the Gulf of Papua. Samples were analysed for molecular organic biomarkers to estimate the sources of organic carbon and its cycling processes.Dry weight fluxes from the shallow traps ranged from 115 to 181 mg m⁻² day⁻¹ and particulate organic carbon (POC) fluxes ranged from 1.2 to 1.9 mM OC m⁻² d⁻¹ with molar organic carbon to particulate nitrogen ratios (C/N) ranging from 6.0 to 6.5. Fluxes in deep traps were likely influenced by both early diagenesis and entrapment of resuspended shelf sediments. Dry weight fluxes in deep traps ranged from 106 to 574 mg m⁻² day⁻¹ and POC fluxes ranged from 0.6 to 1.5 mM OC m⁻² d⁻¹, with C/N ratios ranging from 8.5 to 10.8. ¹³C/¹²C ratios were −20.2‰ to −21.7‰ in all trap samples, indicating that most of the settling POC was “marine-derived”. Shallow traps had δ¹⁵N values of 6.3‰ to 7.2‰ while the values in deep traps were 4.9–5.0‰, indicating the N-rich near-surface OC was less degraded than that in the deep traps. The biogenic lipids consisted of hydrocarbon, sterol and fatty acid biomarkers indicative of marine zooplankton, phytoplankton and bacteria. Sterol markers for diatoms and dinoflagellates were abundant in the water samples. Highly branched isoprenoid alkenes, usually attributable to diatoms, were also detected in both water and shallow traps. Traces of C₂₆–C₃₄ n-alcohols indicative of land–plant biomarkers, were found in river water samples and in the shallow sediment traps. A large unresolved complex mixture (UCM) of hydrocarbons, and a uniform distribution of n-alkanes, indicative of petroleum hydrocarbons, were also detected in the traps. Hopane and sterane biomarkers detected in the trap oil were characteristic of a marine carbonate source, and the aromatic hydrocarbon composition distinguished at least two different oil signatures.We concluded that mass and POC fluxes were similar to those reported for other continental shelves and marginal oceans in tropical and subtropical regions. There was a dramatic decrease in POC as particles sank, due to zooplankton repackaging and photochemical and bacterial decomposition. Carbon isotopic and biomarker patterns showed most of the POC in the sediment traps was marine-sourced with only traces of terrestrial input. There was a significant flux of petroleum, which may signal the existence of natural petroleum seeps in this region.     

Surface measurement system for the atmospheric electrical vertical conduction current density, with displacement current density correction

Global thunderstorm and shower cloud activity generate the global electric potential difference between the Earth's surface and the lower ionosphere. The finite conductivity of atmospheric air, which arises from cosmic ray and natural radioactive ionisation, permits a vertical conduction current density (1 pA m⁻²) between the lower ionosphere and the surface during fair-weather conditions; this current provides a physical link between the upper and lower atmospheres. A new instrument system is described to measure the conduction current density at the surface (the “air–Earth current”), which operates on a novel principle using two collecting electrodes of different geometry. Simultaneous measurements from two independent co-located systems using the geometrical principle show close agreement (correlation of 0.96 during 2.5 h of 5 min measurements). The sensor design described is durable and successful measurements in fair and disturbed weather have been obtained in air temperatures between −6 and 35 °C, relative humidity between 44% and 100%, fog, rain and snowfall. The uncertainty in conduction current density determinations is 0.20 pA m⁻².     

Efficient trapping of organic carbon in sediments on the continental margin with high fluvial sediment input off southwestern Taiwan

Accumulation rates of marine and terrigenous organic carbon in the continental margin sediments off southwestern Taiwan were estimated from the measured concentrations and isotopic compositions of total organic carbon (TOC) and previously reported sedimentation rates. Surficial sediments were collected from the study area spanning from the narrow shelf near the Kaoping River mouth to the deep slope with depths reaching almost 3000 m. The average sediment loading of Kaoping River is 17 Mt/yr, which yields high sediment accumulation rates ranging from 0.08 to 1.44 g cm⁻² yr⁻¹ in the continental margin. About half of the discharged sediments were deposited on the margin within 120 km of the river mouth. Carbon isotopic compositions of terrestrial and marine end-members of organic matter were determined, respectively, based on suspended particulate matter (SPM) collected from three major rivers in the southwestern Taiwan and from an offshore station. All samples were analyzed for the TOC content and its isotopic composition (δ¹³C_org). The SPM samples were also analyzed for the total nitrogen (TN) content. TOC content in marine sediments ranges from 0.45% to 1.35% with the highest values on the upper slope near the Kaoping River mouth. The TOC/TN ratio of the SPM samples from the offshore station is 6.8±0.6, almost identical to the Redfield ratio, indicating their predominantly marine origin; their δ¹³C_org values are also typically marine with a mean of −21.5±0.3‰. The riverine SPM samples exhibit typical terrestrial δ¹³C_org values around −25‰. The δ¹³C_org values of surficial sediments range from −24.8‰ to −21.2‰, showing a distribution pattern influenced by inputs from the Kaoping River. The relative contributions from marine and terrestrial sources to sedimentary organic carbon were determined by the isotope mixing model with end-member compositions derived from the riverine and marine SPM. High fluvial sediment inputs lead to efficient trapping of organic carbon over a wide range of water depth in this continental margin. The marine organic accumulation rate ranges from 1.6 to 70 g C m⁻² yr⁻¹ with an area weighted mean of 4.2 g C m⁻² yr⁻¹, which is on a par with the mean terrestrial contribution and accounts for 2.3% of mean primary production. The depth-dependent accumulation rate of marine organic carbon can be simulated with a function involving primary productivity and mineral accumulation rate, which may be applicable to other continental margins with high sedimentation rates. Away from the nearshore area, the content of terrigenous organic carbon in surficial sediments decreases with distance from the river mouth, indicating its degradation in marine environments.     

Hemolytic activity and fatty acids composition in the ichthyotoxic dinoflagellate Cochlodinium polykrikoides isolated from Bahía de La Paz, Gulf of California

The hemolytic activity of the dinoflagellate Cochlodinium polykrikoides from Bahía de La Paz, Gulf of California was investigated as part of the ichthyotoxic mechanism of this microalga. Two different kinds of erythrocytes, fish and human, were tested for the hemolytic assay. Since fatty acids have been associated with hemolytic activity in C. polykrikoides, the composition of fatty acids of this dinoflagellate was also analyzed. The concentration of C. polykrikoides causing 50% hemolysis (HE₅₀) was 4.88 and 5.27 × 10⁶ cells L⁻¹, for fish and human erythrocytes, respectively. According to the standard curve of saponin, an equivalence between the hemolytic activity of saponin and the dinoflagellate concentration was found with 1 μg saponin mL⁻¹ equivalent to 1 × 10⁶ cells L⁻¹ of C. polykrikoides. The polyunsaturated fatty acids: hexadecaenoic (16:0), docosahexaenoic (22:6 n3) and octadecapentaenoic (18:5 n3) were found in an abundance of 62% of total fatty acids.     

Seasonal variability in mesopause region temperatures over Fort Collins, CO (41°N, 105°W) based on lidar observations, 1990 through 2007

Nearly 900 nocturnal temperature profiles (85–105 km) from the Colorado State University Na lidar at Fort Collins, CO (40.59N, 105.14W) from 1990 to 2007. After the removal of an episodic warming attributable to Mt. Pinatubo eruption, the time series is analyzed as the sum of the climatological mean, annual and semiannual oscillation, solar cycle effect and trends along with possible annual/semiannual modulation of the latter two. The direct seasonal variation is consistent with the concept of the two-level mesopause. The trends in summer and winter are comparable 90–96 km at −0.15±0.1 K/year. The summer trend turns positive above 96 km. The winter trend is negative with minimum of −0.3 K/year at 100 km but positive at 104 km. The negative trend values are a factor of five smaller than an earlier analysis of the early part of this data due to removal of an episodic event.     

Groundwater recharge in an arid grassland as indicated by soil chloride profile and multiple tracers

Previous studies have shown that shallow groundwater in arid regions is often not in equilibrium with near‐surface boundary conditions due to human activities and climate change. This is especially the case where the unsaturated zone is thick and recharge rate is limited. Under this nonequilibrium condition, the unsaturated zone solute profile plays an important role in estimating recent diffuse recharge in arid environments. This paper combines evaluation of the thick unsaturated zone with the saturated zone to investigate the groundwater recharge of a grassland in the arid western Ordos Basin, NW China, using the soil chloride profiles and multiple tracers (²H, ¹⁸O, ¹³C, ¹⁴C, and water chemistry) of groundwater. Whereas conventional water balance and Darcy flux measurements usually involve large errors in recharge estimations for arid areas, chloride mass balance has been widely and generally successfully used. The results show that the present diffuse recharge beneath the grassland is 0.11–0.32 mm/year, based on the chloride mass balance of seven soil profiles. The chloride accumulation age is approximately 2,500 years at a depth of 13 m in the unsaturated zone. The average Cl content in soil moisture in the upper 13 m of the unsaturated zone ranges from 2,842 to 7,856 mg/L, whereas the shallow groundwater Cl content ranges from 95 to 351 mg/L. The corrected ¹⁴C age of shallow groundwater ranges from 4,327 to 29,708 years. Stable isotopes show that the shallow groundwater is unrelated to modern precipitation. The shallow groundwater was recharged during the cold and wet phases of the Late Pleistocene and Holocene humid phase based on palaeoclimate, and consequently, the groundwater resources are nonrenewable. Due to the limited recharge rate and thick unsaturated zone, the present shallow groundwater has not been in hydraulic equilibrium with near‐surface boundary conditions in the past 2,500 years.     

Postglacial isostatic adjustment in a self-gravitating spherical earth with power-law rheology

Since microphysics cannot say definitively whether the rheology of the mantle is linear or non-linear, the aim of this paper is to constrain mantle rheology from observations related to the glacial isostatic adjustment (GIA) process—namely relative sea-levels (RSLs), land uplift rate from GPS and gravity-rate-of-change from GRACE. We consider three earth model types that can have power-law rheology (n = 3 or 4) in the upper mantle, the lower mantle or throughout the mantle. For each model type, a range of A parameter in the creep law will be explored and the predicted GIA responses will be compared to the observations to see which value of A has the potential to explain all the data simultaneously. The coupled Laplace finite-element (CLFE) method is used to calculate the response of a 3D spherical self-gravitating viscoelastic Earth to forcing by the ICE-4G ice history model with ocean loads in self-gravitating oceans. Results show that ice thickness in Laurentide needs to increase significantly or delayed by 2 ka, otherwise the predicted uplift rate, gravity rate-of-change and the amplitude of the RSL for sites inside the ice margin of Laurentide are too low to be able to explain the observations. However, the ice thickness elsewhere outside Laurentide needs to be slightly modified in order to explain the global RSL data outside Laurentide. If the ice model is modified in this way, then the results of this paper indicate that models with power-law rheology in the lower mantle (with A  10⁻³⁵ Pa⁻³ s⁻¹ for n = 3) have the highest potential to simultaneously explain all the observed RSL, uplift rate and gravity rate-of-change data than the other model types.     

Characterising deep vadose zone water movement and solute transport under typical irrigated cropland in the North China Plain

Understanding the dynamics and mechanisms of soil water movement and solute transport is essential for accurately estimating recharge rates and evaluating the impacts of agricultural activities on groundwater resources. In a thick vadose zone (0–15 m) under irrigated cropland in the piedmont region of the North China Plain, soil water content, matric potential, and solute concentrations were measured. Based on these data, the dynamics of soil water and solutes were analysed to investigate the mechanisms of soil water and solute transport. The study showed that the 0–15‐m vadose zone can be divided into three layers: an infiltration and evaporation layer (0–2 m), an unsteady infiltration layer (2–6 m), and a quasi‐steady infiltration layer (6–15 m). The chloride, nitrate, and sulphate concentrations all showed greater variations in the upper soil layer (0–1 m) compared to values in the deep vadose zone (below 2 m). The average concentrations of these three anions in the deep vadose zone varied insignificantly with depth and approached values of 125, 242, and 116 mg/L. The accumulated chloride, sulphate, and nitrate were 2,179 ± 113, 1,760 ± 383, and 4,074 ± 421 kg/ha, respectively. The soil water potential and solute concentrations indicated that uniform flow and preferential flow both occurred in the deep vadose zone, and uniform flow was the dominant mechanism of soil water movement in this study. The piston‐like flow velocity of solute transport was 1.14 m per year, and the average value of calculated leached nitrate nitrogen was 107 kg/ha?year below the root zone. The results can be used to better understand recharge processes and improve groundwater resources management.     

Advantages and challenges of using soil water isotopes to assess groundwater recharge dominated by snowmelt at a field study located in Canada

Seasonal signals of stable isotopes in precipitation, combined with measurements of isotope ratios in soil water, can be used for quantitative estimation of groundwater recharge rates. This study investigates the applicability of using the piston flow principle and the peak shift displacement method to estimate actual groundwater recharge rates in a humid Nordic region located in the province of Quebec, Canada. Two different sites with and without vegetation (C1 and C2) in an unconfined aquifer were tested by measuring soil water isotope ratios (¹⁸O/¹⁶O and ²H/¹H) and volumetric pore water content. Core samples were obtained along the vadose zone down to the groundwater table at the two sites (2.45 m for Site C1 and 4.15 m for Site C2). The peak shift method to estimate groundwater recharge rates was shown to be accurate only in certain specific conditions inherent to the soil properties and the topographical situation of the investigated sites. Indeed, at Site C2, recharge from the snowmelt could not be estimated because of heterogeneity in the lower part of the vadose zone. At this same site the later recharge after the snowmelt (in the period from late spring to early autumn) could be estimated accurately because the upper part of the vadose zone was homogeneous. Furthermore, at site C1, runoff/runon phenomena hampered calculations of actual infiltration and thus produced inaccurate results for recharge. These two different site effects (heterogeneity in the first site and runoff/runon in the other site) were identified as being limiting factors in the accurate assessment of actual recharge. This study therefore recommends the use of the peak shift method for (1) humid Nordic regions, (2) homogeneous and thick vadose zones, and (3) areas with few or limited site effects (runoff/runon).     

Multivariate statistical analysis of geochemical data as indicative of the hydrogeochemical evolution of groundwater in a sedimentary rock aquifer system

The study of groundwater hydrogeochemistry of the Paleozoic Basses-Laurentides sedimentary rock aquifer system in Québec produced a large geochemical dataset. Groundwater samples were collected at 153 sites over a 1500 km² study area and analyzed for major and minor ions. The large number of data can lead to difficulties in the integration, interpretation and representation of the results. Two multivariate statistical methods, hierarchical cluster analysis (HCA) and principal components analysis (PCA), were applied to a subgroup of the dataset to evaluate their usefulness to classify the groundwater samples, and to identify geochemical processes controlling groundwater geochemistry. This subgroup consisted of 144 samples and 14 parameters (Ca²⁺, Mg²⁺, Na⁺, K⁺, , Cl⁻, , Fe²⁺, Mn²⁺, Br⁻, Sr²⁺, F⁻, Ba²⁺, HS⁻). Seven geochemically distinct clusters, C1–C7, resulted from the HCA. Samples from clusters C3, C4, C6 and C7 are mostly located in preferential recharge areas. The majority of these samples have Ca–Mg–HCO₃ recharge groundwater (C3, C6, C7) and Na–HCO₃ evolved groundwater (C4). Samples from the other three clusters (C1, C2, C5) are characteristic of an aquifer system under confined conditions. The majority of these samples have Na–HCO₃ evolved groundwater (C1, C5) and Na–Cl ancient groundwater that exhibits elevated concentrations in Br⁻ (C2). In addition to recognizing the importance of hydrogeological conditions on groundwater geochemistry, the distribution of clusters also showed the importance of the geological formations on minor and trace elements, such as Fe²⁺, Mn²⁺, Sr²⁺, F⁻ and Ba²⁺. The first five components of the PCA account for 78.3% of the total variance in the dataset. Component 1 is defined by highly positive loadings in Na⁺, Cl⁻ and Br⁻ and is related to groundwater mixing with Champlain Sea water and solute diffusion from the marine clay aquitard. The high positive loadings in Ca²⁺ and Mg²⁺ of component 2 suggest the importance of dissolution of carbonate rocks in this aquifer system. From their characteristic loadings, the first two components are defined as the “salinity” and “hardness” components, respectively. Components 3–5 are related to more local and geological effects. The integration of the HCA and the PCA, with conventional classification of groundwater types, as well as with the hydrogeological and geological contexts, allowed the division of the region into four main geochemical areas, providing an improved regional picture of the aquifer system dynamics and hydrogeochemical evolution of groundwater. The following factors were recognized as influencing the evolution of groundwater identified in every geochemical area: (1) geological characteristics including sedimentary rock type and till mineralogy; (2) hydrogeological characteristics represented by the level of confinement and the hydraulic gradient; and (3) the geological history including the latest glaciation and the Champlain Sea invasion. With its integrated approach, this hydrogeochemical study contributes to the characterization and understanding of complex groundwater flow systems, and provides an example of the long-term geochemical evolution of hydrogeological systems after a major perturbation, in this case seawater invasion.     

Transport of oceanic nitrate from the continental shelf to the coastal basin in relation to the path of the Kuroshio

Hydrographic and biogeochemical observations were conducted along the longitudinal section from Ise Bay to the continental margin (southern coast of Japan) to investigate changes according to the Kuroshio path variations during the summer. The strength of the uplift of the cold deep water was influenced by the surface intrusion of the Kuroshio water to the shelf region. When the intrusion of the Kuroshio surface water to the shelf region was weak in 2006, the cold and NO₃⁻-rich shelf water intruded into the bottom layer in the bay from the shelf. This bottom intrusion was intensified by the large river discharge. The nitrogen isotope ratio (δ¹⁵N) of NO₃⁻ (4–5‰) in the bottom bay water was same as that in the deeper NO₃⁻ over the shelf, indicating the supply of new nitrogen to the bay. The warm and NO₃⁻-poor shelf water intruded into the middle layer via the mixing region at the bay mouth when the Kuroshio water distributed in the coastal areas off Ise Bay in 2005. The regenerated NO₃⁻ with isotopically light nitrogen (δ¹⁵N=−1‰) was supplied from the shelf to the bay. This NO₃⁻ is regenerated by the nitrification in the upper layer over the shelf. The contribution rate of regenerated NO₃⁻ over the shelf to the total NO₃⁻ in the subsurface chlorophyll maximum layer in the bay was estimated at 56% by a two-source mixing model coupled with the Rayleigh equation.     

Hydrothermal activity on the southern Mid-Atlantic Ridge: Tectonically- and volcanically-controlled venting at 4–5°S

We report results from an investigation of the geologic processes controlling hydrothermal activity along the previously-unstudied southern Mid-Atlantic Ridge (3–7°S). Our study employed the NOC (UK) deep-tow sidescan sonar instrument, TOBI, in concert with the WHOI (USA) autonomous underwater vehicle, ABE, to collect information concerning hydrothermal plume distributions in the water column co-registered with geologic investigations of the underlying seafloor. Two areas of high-temperature hydrothermal venting were identified. The first was situated in a non-transform discontinuity (NTD) between two adjacent second-order ridge-segments near 4°02′S, distant from any neovolcanic activity. This geologic setting is very similar to that of the ultramafic-hosted and tectonically-controlled Rainbow vent-site on the northern Mid-Atlantic Ridge. The second site was located at 4°48′S at the axial-summit centre of a second-order ridge-segment. There, high-temperature venting is hosted in an  18 km² area of young lava flows which in some cases are observed to have flowed over and engulfed pre-existing chemosynthetic vent-fauna. In both appearance and extent, these lava flows are directly reminiscent of those emplaced in Winter 2005−06 at the East Pacific Rise, 9°50′N and reference to global seismic catalogues reveals that a swarm of large (M 4.6−5.6) seismic events was centred on the 5°S segment over a  24 h period in late June 2002, perhaps indicating the precise timing of this volcanic eruptive episode. Temperature measurements at one of the vents found directly adjacent to the fresh lava flows at 5°S MAR (Turtle Pits) have subsequently revealed vent-fluids that are actively phase separating under conditions very close to the Critical Point for seawater, at  3000 m depth and 407 °C: the hottest vent-fluids yet reported from anywhere along the global ridge crest.     

Analysis of radionuclide migration through a 200-m Vadose zone following a 16-year infiltration event

Atomic weapons testing at the Nevada Test Site has introduced many tracers for quantifying subsurface hydrologic transport processes in arid climates. In 1975, groundwater adjacent to the Cambric test, conducted beneath Frenchman Flat 10 years earlier, was pumped steadily for 16 years to elicit information on the migration of residual radioactivity through the saturated zone. Radionuclides in the pumping well effluent, including tritium, ¹⁴C, ³⁶Cl, and ⁸⁵Kr, were extensively monitored prior to its discharge in an unlined ditch, where approximately a third of the flow infiltrated over a distance of 1 km. Radionuclide infiltration through a 220-m thick vadose zone created a second, and rather unique long-term field experiment. Effluent data have been utilized in conjunction with geologic data, new radionuclide measurements, isotopic age-dating estimates, and vadose zone flow and transport models to better understand the movement of radionuclides between the ditch, the water table and a nearby groundwater monitoring well. Detection of tritium in the monitoring well occurred approximately 16 years after its initial discharge into the ditch. Modeling and tritium age dating have suggested 3–5 years of this 16-year transit time occurred solely in the vadose zone. They also suggest considerable recirculation of the pumping well discharge back into the original pumping well. Notably, there have been no observations of ¹⁴C or ⁸⁵Kr in the monitoring well, suggesting their preferential retention or volatilization during transit to the water table.     

Rapid Cleanup of a Perchlorate Plume from Fireworks

Perchlorate was detected in a municipal wellfield in Evart, Michigan in April 2015. Perchlorate concentrations were detected initially in six of the City's wells at concentrations ranging up to 20 μg/L. An investigation to identify the source determined that the perchlorate was from fireworks launched during the annual 4th of July show held at the fairgrounds located upgradient from the wellfield. The use of approximately 600 kg of fireworks during the annual display resulted in an annual loading of approximately 4 kg of perchlorate to groundwater. An aggressive groundwater extraction system began operation in June 2016 to restore water quality in the affected aquifer, and the 2016 fireworks display was relocated to a location outside the capture zone of the water supply wells. Within 18 months average perchlorate concentrations in the water supply wells had been reduced to about 0.6 μg/L. The extraction system continued to operate through the end of 2019, by which time the average perchlorate concentrations in the water supply wells were reduced to 0.2 μg/L. In 2019, approximately 0.4 kg of perchlorate were removed from the aquifer, about one-half of the amount removed in 2018, reflecting the slow leaching of perchlorate of fireworks residuals from vadose zone soils.     

Mixing at the interface between fresh and salt waters in 3D steady flow with application to a pumping well in a coastal aquifer

We consider 3D steady flow of fresh water over a salt water body in a confined aquifer of constant thickness D, with application to a pumping well in a coastal aquifer. With neglect of mixing, a sharp interface separates the two fluid bodies and an existing analytical solution, based on the Dupuit assumption, is adopted. The aim is to solve for the mixing between the fresh and salt waters for α_T/D  1 (α_T transverse dispersivity), as field studies indicate that α_T = O(10⁻³ − 10⁻² m). The mixing zone around the interface is narrow and solutions by existing codes experience numerical difficulties. The problem is solved by the boundary layer (BL) approximation, extending a method, applied previously to two-dimensional flows. The BL equations of variable-density flow are solved by using the Von Karman integral method, to determine the BL thickness and the rate of entrainment of salt water along the interface. Application to the pumping well problem yields the salinity of the pumped water, as function of the parameters of the problem (well discharge, seaward discharge, well distance from the coast and density difference).     

The emplacement of an obsidian dyke through thin ice: Hrafntinnuhryggur, Krafla Iceland

An eruption along a 2.5 km-long rhyolitic dyke at Krafla volcano, northern Iceland during the last glacial period formed a ridge of obsidian (Hrafntinnuhryggur). The ridge rises up to 80 m above the surrounding land and is composed of a number of small-volume lava bodies with minor fragmental material. The total volume is < 0.05 km³. The lava bodies are flow- or dome-like in morphology and many display columnar-jointed sides typical of magma–ice interaction, quench-fragmented lower margins indicative of interaction with meltwater and pumiceous upper surfaces typical of subaerial obsidian flows. The fragmental material compromises poorly-sorted perlitic quench hyaloclastites and poorly-exposed pumiceous tuffs. Lava bodies on the western ridge flanks are columnar jointed and extensively hydrothermally altered. At the southern end of the ridge the feeder dyke is exposed at an elevation  95 m beneath the ridge crest and flares upwards into a lava body.Using the distribution of lithofacies, we interpret that the eruption melted through ice only 35–55 m thick, which is likely to have been dominated by firn. Hrafntinnuhryggur is therefore the first documented example of a rhyolitic fissure eruption beneath thin ice/firn. The eruption breached the ice, leading to subaerial but ice/firn-contact lava effusion, and only minor explosive activity occurred. The ridge appears to have been well-drained during the eruption, aided by the high permeability of the thin ice/firn, which appears not to have greatly affected the eruption mechanisms. We estimate that the eruption lasted between 2 and 20 months and would not have generated a significant jökulhlaup (< 70 m³ s^− 1).     

Progressive cooling of the hyaloclastite ridge at Gjálp, Iceland, 1996–2005

In the subglacial eruption at Gjálp in October 1996 a 6 km long and 500 m high subglacial hyaloclastite ridge was formed while large volumes of ice were melted by extremely fast heat transfer from magma to ice. Repeated surveying of ice surface geometry, measurement of inflow of ice, and a full Stokes 2-D ice flow model have been combined to estimate the heat output from Gjálp for the period 1996–2005. The very high heat output of order 10⁶ MW during the eruption was followed by rapid decline, dropping to  2500 MW by mid 1997. It remained similar until mid 1999 but declined to 700 MW in 1999–2001. Since 2001 heat output has been insignificant, probably of order 10 MW. The total heat carried with the 1.2 × 10¹² kg of basaltic andesite erupted (0.45 km³ DRE) is estimated to have been 1.5 × 10¹⁸ J. About two thirds of the thermal energy released from the 0.7 km³ edifice in Gjálp occurred during the 13-day long eruption, 20% was released from end of eruption until mid 1997, a further 10% in 1997–2001, and from mid 2001 to present, only a small fraction remained. The post-eruption heat output history can be reconciled with the gradual release of 5 × 10¹⁷ J thermal energy remaining in the Gjálp ridge after the eruption, assuming single phase liquid convection in the cooling edifice. The average temperature of the edifice is found to have been approximately 240 °C at the end of the eruption, dropping to  110 °C after 9 months and reaching  40 °C in 2001. Although an initial period of several months of very high permeability is possible, the most probable value of the permeability from 1997 onwards is of order 10^− 12 m². This is consistent with consolidated/palagonitized hyaloclastite but incompatible with unconsolidated tephra. This may indicate that palagonitization had advanced sufficiently in the first 1–2 years to form a consolidated hyaloclastite ridge, resistant to erosion. No ice flow traversing the Gjálp ridge has been observed, suggesting that it has effectively been shielded from glacial erosion in its first 10 years of existence.     

Changes in groundwater bacterial community during cyclic groundwater-table variations

Column experiments containing an aquifer sand were subjected to static and oscillating water tables to investigate the impact of natural fluctuations and rainfall infiltration on the groundwater bacterial community just below the phreatic surface, and its association with the geochemistry. Once the columns were established, the continuously saturated zone was anoxic in all three columns. The rate of soil organic matter (SOM) mineralization was higher when the water table varied cyclically than when it was static due to the greater availability of NO₃⁻ and SO₄²⁻. Natural fluctuations in the water table resulted in a similar NO₃⁻ concentration to that observed with a static water table but the cyclic wetting of the intermittently saturated zone resulted in a higher SO₄²⁻ concentration. Rainfall infiltration induced cyclic water-table variations resulted in a higher NO₃⁻ concentration than those in the other two columns, and a SO₄²⁻ concentration intermediate between those columns. As rainwater infiltration resulted in slow downward displacement of the groundwater, it is inferred that NO₃⁻ and SO₄²⁻ were being mobilized from the vadose zone. NO₃⁻ was mainly released by SOM mineralization (which was enhanced by the infiltration of oxygenated rainwater), but the larger amount of SO₄²⁻ release required a second mechanism (possibly desorption). Different groundwater bacterial communities evolved from initially similar populations due to the different groundwater histories.     

Geological controls on groundwater chemistry and arsenic mobilization: Hydrogeochemical study along an E–W transect in the Meghna basin, Bangladesh

Hydrogeochemical investigations along an E–W transect in the middle Meghna basin show groundwater chemistry and redox condition vary considerably with the change in geology. Groundwater in the Holocene shallow (<150 m bgl) alluvial aquifer in western part of the transect is affected by high arsenic concentration (As > 10 μg/l) and salinity. On the other hand, groundwater from the Pliocene Dupi Tila sandy aquifer in the eastern part is fresh and low in As (<10 μg/l). The Holocene shallow aquifers are high in dissolved As, , Fe and dissolved organic carbon (DOC), but generally low in and . High concentrations (250–716 mg/l) together with high DOC concentrations (1.4–21.7 mg/l) in these aquifers reflect active sources of degradable natural organic matter that drives the biogeochemical process. There is generally de-coupling of As from other redox-sensitive elements. In contrast, the Pliocene aquifers are low in As, and DOC. Molar ratio of /H₄SiO₄ suggests that silicate weathering is dominant in the deeper Holocene aquifers and in the Pliocene aquifers. Molar ratios of Cl⁻/ and Na⁺/Cl⁻ suggest mixing of relict seawater with the fresh water as the origin of groundwater salinity. Speciation calculations show that saturation indices for siderite and rhodochrosite vary significantly between the Holocene and Pliocene aquifers. Stable isotopes (δ²H and δ¹⁸O) in groundwater indicate rapid infiltration without significant effects of evaporation. The isotopic data also indicates groundwater recharge from monsoonal precipitation with some impact of altitude effect at the base of the Tripura Hills in the east. The results of the study clearly indicate geological control (i.e. change in lithofacies) on groundwater chemistry and distribution of redox-sensitive elements such as As along the transect.