Search USGS⌕ Search

SEARCH · Search USGS

Results for “The Leading Edge”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 73 records · Page 4Linked to original sources

Channel response to a dam‐removal sediment pulse captured at high‐temporal resolution using routine gage data

In this study, we captured how a river channel responds to a sediment pulse originating from a dam removal using multiple lines of evidence derived from streamflow gages along the Patapsco River, Maryland, USA. Gages captured characteristics of the sediment pulse, including travel times of its leading edge (~7.8 km yr −1 ) and peak (~2.6 km yr −1 ) and suggest both translation and increasing dispersion. The pulse also changed local hydraulics and energy conditions, increasing flow velocities and Froude number, due to bed fining, homogenization and/or slope adjustment. Immediately downstream of the dam, recovery to pre‐pulse conditions occurred within the year, but farther downstream recovery was slower, with the tail of the sediment pulse working through the lower river by the end of the study 7 years later. The patterns and timing of channel change associated with the sediment pulse were not driven by large flow or suspended sediment‐transporting events, with change mostly occurring during lower flows. This suggests pulse mobility was controlled by process‐factors largely independent of high flow. In contrast, persistent changes occurred to out‐of‐channel flooding dynamics. Stage associated with flooding increased during the arrival of the sediment pulse, 1 to 2 years after dam removal, suggesting persistent sediment deposition at the channel margins and nearby floodplain. This resulted in National Weather Service‐indicated flood stages being attained by 3–43% smaller discharges compared to earlier in the study period. This study captured a two‐signal response from the sediment pulse: (1) short‐ to medium‐term (weeks to months) translation and dispersion within the channel, resulting in aggradation and recovery of bed elevations and changing local hydraulics; and (2) dispersion and persistent longer‐term (years) effects of sediment deposition on overbank surfaces. This study further demonstrated the utility of US Geological Survey gage data to quantify geomorphic change, increase temporal resolution, and provide insights into trajectories of change over varying spatial and temporal scales.

Maryland↗

Geology, distribution, and classification of gold deposits in the western Qinling belt, central China

Gold deposits of the western Qinling belt occur within the western part of the Qinling-Dabie-Sulu orogen, which is located between the Precambrian North China and Yangtze cratons and east of the Songpan-Ganzi basin. The early Paleozoic to early Mesozoic orogen can be divided into northern, central, and southern zones, separated by the Shangdan and Lixian-Shanyang thrust fault systems. The northern zone consists of an early Paleozoic arc accreted to the North China craton by ca. 450 Ma. The central zone, which contains numerous orogenic gold deposits, is dominated by clastic rocks formed in a late Paleozoic basin between the converging cratonic blocks. The southern zone is characterized by the easternmost exposure of Triassic sedimentary rocks of the Songpan-Ganzi basin. These Early to Late Triassic turbidities, in part calcareous, of the immense Songpan-Ganzi basin also border the western Qinling belt to the west. Carlinlike gold deposits are abundant (1) along a westward extension of the southern zone defined by a window of early Paleozoic clastic rocks extending into the basin, and (2) within the easternmost margin of the basinal rocks to the south of the extension, and in adjacent cover rocks of the Yangtze craton. Triassic and Early Jurassic synkinematic granitoids are widespread across the western Qinling belt, as well as in the Songpan-Ganzi basin. Orogenic lode gold deposits along brittle-ductile shear zones occur within greenschist-facies, highly deformed, Devonian and younger clastic rocks of the central zone. Mainly coarse-grained gold, along with pyrite, pyrrhotite, arsenopyrite, and minor base metal sulfides, occur in networks of quartz veinlets, brecciated wall rock, and are dissminated in altered wall rock. Isotopic dates suggest that the deposits formed during the Late Triassic to Middle Jurassic as the leading edge of the Yangtze craton was thrust beneath rocks of the western Qinling belt. Many gold-bearing placers are distributed along the river systems that flow south from the lode-bearing central zone. Carlin-like gold deposits have only been identified during the last decade in the southern zone of the western Qinling and in the northeastern corner of the Songpan-Ganzi basin. The deposits mainly contain micron-diameter gold in arsenical pyrite; are characterized by the common occurrence of cinnabar, stibnite, realgar, and orpiment; exhibit strong silicification, carbonatization, pyritization, and decalcification dissolution textures; and are structurally controlled. The lack of reactive host lithologies may have prevented development of large (> 100 tones of gold), stratigraphically-controlled orebodies, which are typical of the Carlin deposits in the western USA. These deposits are hosted by Triassic turbidities and shallow-water carbonates, and an early Paleozoic inlier in the Songpan-Ganzi basin that extends in an east-west belt for about 300 km. Rather than true "Carlin" deposits, these Carlin-like deposits may be some type of shallow-crustal (i.e., epithermal) hybrid with features intermediate to Nevada-style Carlin deposits and the orogenic gold deposits to the immediate north. These Carlin-like deposits also overlap in age with the early Mesozoic orogenic gold deposits and, therefore, also formed during the final stages of collision between the cratons and intermediate basin closure.

Mineralium Deposita↗

Pathogen transport in groundwater systems: Contrasts with traditional solute transport

Water quality affects many aspects of water availability, from precluding use to societal perceptions of fit-for-purpose. Pathogen source and transport processes are drivers of water quality because they have been responsible for numerous outbreaks resulting in large economic losses due to illness and, in some cases, loss of life. Outbreaks result from very small exposure (e.g., less than 20 viruses) from very strong sources (e.g., trillions of viruses shed by a single infected individual). Thus, unlike solute contaminants, an acute exposure to a very small amount of contaminated water can cause immediate adverse health effects. Similarly, pathogens are larger than solutes. Thus, interactions with surfaces and settling become important even as processes important for solutes such as diffusion become less important. These differences are articulated in “Colloid Filtration Theory”, a separate branch of pore-scale transport. Consequently, understanding pathogen processes requires changes in how groundwater systems are typically characterized, where the focus is on the leading edges of plumes and preferential flow paths, even if such features move only a very small fraction of the aquifer flow. Moreover, the relatively short survival times of pathogens in the subsurface require greater attention to very fast (<10 year) flow paths. By better understanding the differences between pathogen and solute transport mechanisms discussed here, a more encompassing view of water quality and source water protection is attained. With this more holistic view and theoretical understanding, better evaluations can be made regarding drinking water vulnerability and the relation between groundwater and human health.

Hydrogeology Journal↗

Heating of a fully saturated Darcian half-space: Pressure generation, fluid expulsion, and phase change

Analytical solutions are developed for the pressurization, expansion, and flow of one- and two-phase liquids during heating of fully saturated and hydraulically open Darcian half-spaces subjected to a step rise in temperature at its surface. For silicate materials, advective transfer is commonly unimportant in the liquid region; this is not always the case in the vapor region. Volume change is commonly more important than heat of vaporization in determining the position of the liquid-vapor interface, assuring that the temperatures cannot be determined independently of pressures. Pressure increases reach a maximum near the leading edge of the thermal front and penetrate well into the isothermal region of the body. Mass flux is insensitive to the hydraulic properties of the half-space.

International Journal of Heat and Mass Transfer↗

The geomorphology of the Mississippi River chenier plain

The chenier plain of the Mississippi River is a shore-parallel zone of alternating transgressive clastic ridges separated by progradational mudflats. The term chenier is derived from the cajun term chene for oak, the tree species that colonizes the crests of the higher ridges. The Mississippi River chenier plain stretches 200 km from Sabine Pass, Texas, to Southwest Point, Louisiana and ranges between 20 and 30 km wide, with elevations of 2–6 m. The timing and the process of formation could be re-evaluated in the light of new chronostratigraphic findings in the Mississippi River delta plain. The stratigraphic relationship between the Teche and Lafourche delta complexes and Ship Shoal offshore indicates that these delta complexes belong to different delta plains that developed at different sealevels. It appears that the Teche delta complex is associated with the late Holocene delta plain which developed 7000 to 3000 yrs B.P. when sealevel stood 5–6 m lower than present. A regional transgression occurred between approximately 3000 BP and 2500 yrs B.P., leading to the transgressive submergence of the late Holocene delta plain, producing the regional Teche shoreline. The timing of this transgression conforms to the age of the most landward ridge in the chenier plain, the Little Chenier-Little Pecan Island trend, which dates at about 2500 yrs B.P. This ridge trend was originally interpreted as representing the Teche delta complex switching event with the landward Holocene/Pleistocene contact representing the high stand shoreline. The implication of this new interpretation is that the Little Chenier-Little Pecan Island trend represents the high stand shoreline, a continuation of the Teche shoreline separating the late Holocene and Recent delta plains, and that the Holocene/Pleistocene contact represents the leading edge of the marshes transgressing onto the Prairie Terrace. Significant mudflat progradation seems to require a westerly position of the Mississippi River, but the numerous different forms and ages of cheniers do not correspond well to the timing of major delta complex switching. Progradation of the chenier plain appears to be associated with building of the Recent delta plain and not the Teche complex of the late Holocene delta plain. The occurrence of individual ridges appears to be primarily tied to delta lobe switching within the Lafourche complex and variations in sediment supply from local rivers. The recent development of the Atchafalaya delta complex to the west is the closest position of an active distributary to the chenier plain since sealevel stabilization; a new episode of rapid mudflat progradation is thus taking place.

Marine Geology↗

Tectonic erosion at the front of the Japan trench convergent margin

The imaging of a multichannel seismic record was improved by reprocessing using pre-stack techniques. The reprocessed record shows structures that indicate tectonic erosion and gravity collapse at the front of the Japan trench margin. Much of the lower slope appears to be underlain by a detached, coherent block of continental crust. The lower slope has failed by mass wasting and the resulting apron of slump debris at the base of the slope has become involved in thrust faulting at the front of the subduction zone. Slumping continues as long as debris is removed from the front of the margin by subduction, and the apron cannot build up sufficiently to stabilize the failing lower slope. Truncated beds at the base of the upper plate indicate subcrustal erosion as well, this probably being the main cause of massive subsidence of the margin. Subsidence was the cause of oversteepening, destabilization and subsequent gravity collapse of the leading edge of the upper plate.

Tectonophysics↗

Composite transform-convergent plate boundaries: description and discussion

The leading edge of the overriding plate at an obliquely convergent boundary is commonly sliced by a system of strike-slip faults. This fault system is often structurally complex, and may show correspondingly uneven strain effects, with great vertical and translational shifts of the component blocks of the fault system. The stress pattern and strain effects vary along the length of the system and change through time. These margins are considered to be composite transform-convergent (CTC) plate boundaries. Examples are given of structures formed along three CTC boundaries: the Aleutian Ridge, the Solomon Islands, and the Philippines. The dynamism of the fault system along a CTC boundary can enhance vertical tectonism and basin formation. This concept provides a framework for the evaluation of petroleum resources related to basin formation, and mineral exploration related to igneous activity associated with transtensional processes. ?? 1992.

Marine and Petroleum Geology↗

Crude oil in a shallow sand and gravel aquifer-II. Organic geochemistry

Crude oil spilled from a pipeline break in a remote area of north-central Minnesota has contaminated a shallow glacial outwash aquifer. Part of the oil was sprayed over a large area to the west of the pipeline and part of it accumulated in an oil body that floats at the water table to the east of the point of discharge. Total dissolved organic carbon (TDOC) concentrations in shallow groundwater collected in the oil spray area reach 16 mg/l. This is nearly an order of magnitude higher than the TDOC concentrations of native groundwater (∼2–3mg/l). The additional TDOC derives from the partial degradation of petroleum residues deposited at the land surface and transported to the aquifer by vertical recharge. In the vicinity of the oil body, TDOC concentrations in groundwater are 48 mg/l, 58% of the TDOC being composed of non-volatile organic C. The majority of the volatile DOC (63%) is a mixture of low-molecular-weight saturated, aromatic and alicyclic hydrocarbons derived from the oil. Downgradient from the oil body along the direction of groundwater flow, concentrations of all measured constituents of the TDOC pool decrease. Concentrations begin to decline most rapidly, however, in the zone where dissolved O 2 concentrations begin to increase, ∼50m downgradient from the leading edge of the oil. Within the anoxic zone near the oil body, removal rates of isometric monoaromatic hydrocarbons vary widely. This indicates that the removal processes are mediated mainly by microbiological activity. Molecular and spectroscopic characterization of the TDOC and its spatial and temporal variation provide evidence of the importance of biogeochemical processes in attenuating petroleum contaminants in this perturbed subsurface environment.

Applied Geochemistry↗

Stratigraphy and structural setting of Upper Cretaceous Frontier Formation, western Centennial Mountains, southwestern Montana and southeastern Idaho

Stratigraphic, sedimentologic, and palynologic data were used to correlate the Frontier Formation of the western Centennial Mountains with time-equivalent rocks in the Lima Peaks area and other nearby areas in southwestern Montana. The stratigraphic interval studied is in the middle and upper parts (but not uppermost) of the formation based on a comparison of sandstone petrography, palynologic age data, and our interpretation of the structure using a seismic line along the frontal zone of the Centennial Mountains and the adjacent Centennial Valley. The Frontier Formation is comprised of sandstone, siltstone, mudstone, limestone, and silty shale in fluvial and coastal depositional settings. A distinctive characteristic of these strata in the western Centennial Mountains is the absence of conglomerate and conglomeratic sandstone beds. Absence of conglomerate beds may be due to lateral facies changes associated with fluvial systems, a distal fining of grain size, and the absence of both uppermost and lower Frontier rocks in the study area. Palynostratigraphic data indicate a Coniacian age for the Frontier Formation in the western Centennial Mountains. These data are supported by a geochronologic age from the middle part of the Frontier at Lima Peaks indicating a possible late Coniacian-early Santonian age (86.25 ?? 0.38 Ma) for the middle Frontier there. The Frontier Formation in the western Centennial Mountains is comparable in age and thickness to part of the Frontier at Lima Peaks. These rocks represent one of the thickest known sequences of Frontier strata in the Rocky Mountain region. Deposition was from about 95 to 86 Ma (middle Cenomanian to at least early Santonian), during which time, shoreface sandstone of the Telegraph Creek Formation and marine shale of the Cody Shale were deposited to the east in the area now occupied by the Madison Range in southwestern Montana. Frontier strata in the western Centennial Mountains are structurally isolated from other Cretaceous rocks in the region and are part of the Lima thrust sheet that lies at the leading edge of the Sevier-style overthrusting in this part of southwestern Montana and adjacent southeastern Idaho.

Cretaceous Research↗

Foundational uncertainties in terminal Ediacaran chronostratigraphy revealed by high-precision zircon U-Pb geochronology of the Nama Group, Namibia

The Nama Group of southern Namibia and northwestern South Africa hosts the best-dated mixed carbonate-siliciclastic foreland basin succession of the terminal Ediacaran [ca. 551 million years (Ma) ago to <538 Ma] and is key for resolving the chronology of early metazoan evolution. Numerous silicified volcanic tuff interbeds are present, but differing interpretations regarding the fidelity of their ages lead to different regional stratigraphic correlations, especially for the Urusis Formation of the Schwarzrand Subgroup. An expanded record of the Urusis Formation is found in the Swartpunt area of southern Namibia, which has yielded an important metazoan biota. But the succession in this area is preserved as a series of thrusts at the leading edge of the Gariep orogenic belt and zircon U-Pb data show systematic age repetition. We use regional stratigraphic and structural mapping, integrated with carbonate carbon isotope (δ 13 C carb ) chemostratigraphy and high-precision radioisotope U-Pb zircon geochronology from outcrop and recently acquired drill core to develop a temporally calibrated basin-wide depositional model. This integrated dataset either reflects complex zircon reworking, inheritance, or potential analytical biases (Scenario 1) or the presence of a Gariep-related cryptic décollement within the Spitskop Member that has resulted in stratigraphic repetition (Scenario 2). We investigate the evidence for and against both scenarios and consider their implications for stratigraphic and δ 13 C carb correlations between the Swartpunt area and coeval exposures along the Orange River border with South Africa. Given that these issues are in an area that hosts numerous silicified ash beds and extensive exposure, an inability to confidently discount either scenario highlights a level of compounding uncertainty in zircon U-Pb geochronology that must be considered when attempting to build global chronostratigraphic frameworks. Scenario 1 implies that some of the weighted mean ages and Bayesian eruption ages from the Swartpunt area may be >1 Myr older than the depositional age of their respective ash beds when assuming existing stratigraphic correlations. If this scenario is preferred, then a cautious approach would be to consider all weighted mean zircon U-Pb ages from ash beds to reflect maximum depositional ages. Both scenarios support deposition of the Huns Member >540 Ma in the Swartpunt area if the oldest weighted mean age reported here represents a near-depositional age, which has significant implications for the temporal calibration of important terminal Ediacaran ichnofossil assemblages and future cyclostratigraphic studies. Stratigraphic correlations common to both scenarios allow us to temporally calibrate a basin evolution model for the Nama Group. Temporal trends in initial hafnium isotope (εHf) compositions of zircon grains from ash beds throughout the succession may support progressive crustal thickening associated with underplating of the Damara orogenic belt along the northern periphery of the Kalahari craton from ca. 547 Ma to ca. 538 Ma. The compilation of new and published zircon U-Pb ages may also imply that the locus of carbonate platform development migrated from north to south (present co-ordinates), tracking the migration of foredeep subsidence.

Earth-Science Reviews↗

A mass balance approach to investigate arsenic cycling in a petroleum plume

Natural attenuation of organic contaminants in groundwater can give rise to a series of complex biogeochemical reactions that release secondary contaminants to groundwater. In a crude oil contaminated aquifer, biodegradation of petroleum hydrocarbons is coupled with the reduction of ferric iron (Fe(III)) hydroxides in aquifer sediments. As a result, naturally occurring arsenic (As) adsorbed to Fe(III) hydroxides in the aquifer sediment is mobilized from sediment into groundwater. However, Fe(III) in sediment of other zones of the aquifer has the capacity to attenuate dissolved As via resorption. In order to better evaluate how long-term biodegradation coupled with Fe-reduction and As mobilization can redistribute As mass in contaminated aquifer, we quantified mass partitioning of Fe and As in the aquifer based on field observation data. Results show that Fe and As are spatially correlated in both groundwater and aquifer sediments. Mass partitioning calculations demonstrate that 99.9% of Fe and 99.5% of As are associated with aquifer sediment. The sediments act as both sources and sinks for As, depending on the redox conditions in the aquifer. Calculations reveal that at least 78% of the original As in sediment near the oil has been mobilized into groundwater over the 35-year lifespan of the plume. However, the calculations also show that only a small percentage of As (∼0.5%) remains in groundwater, due to resorption onto sediment. At the leading edge of the plume, where groundwater is suboxic, sediments sequester Fe and As, causing As to accumulate to concentrations 5.6 times greater than background concentrations. Current As sinks can serve as future sources of As as the plume evolves over time. The mass balance approach used in this study can be applied to As cycling in other aquifers where groundwater As results from biodegradation of an organic carbon point source coupled with Fe reduction.

Environmental Pollution↗

Microbial characterization of nitrification in a shallow, nitrogen-contaminated aquifer, Cape Cod, Massachusetts and detection of a novel cluster associated with nitrifying Betaproteobacteria

Groundwater nitrification is a poorly characterized process affecting the speciation and transport of nitrogen. Cores from two sites in a plume of contamination were examined using culture-based and molecular techniques targeting nitrification processes. The first site, located beneath a sewage effluent infiltration bed, received treated effluent containing O 2 (> 300 µM) and NH 4 + (51–800 µM). The second site was 2.5 km down-gradient near the leading edge of the ammonium zone within the contaminant plume and featured vertical gradients of O 2 , NH 4 + , and NO 3 − (0–300, 0–500, and 100–200 µM with depth, respectively). Ammonia- and nitrite-oxidizers enumerated by the culture-based MPN method were low in abundance at both sites (1.8 to 350 g − 1 and 33 to 35,000 g − 1 , respectively). Potential nitrifying activity measured in core material in the laboratory was also very low, requiring several weeks for products to accumulate. Molecular analysis of aquifer DNA (nested PCR followed by cloning and 16S rDNA sequencing) detected primarily sequences associated with the Nitrosospira genus throughout the cores at the down-gradient site and a smaller proportion from the Nitrosomonas genus in the deeper anoxic, NH 4 + zone at the down-gradient site. Only a single Nitrosospira sequence was detected beneath the infiltration bed. Furthermore, the majority of Nitrosospira -associated sequences represent an unrecognized cluster. We conclude that an uncharacterized group associated with Nitrosospira dominate at the geochemically stable, down-gradient site, but found little evidence for Betaproteobacteria nitrifiers beneath the infiltration beds where geochemical conditions were more variable.

Journal of Contaminant Hydrology↗

Submarine geology of Hana Ridge and Haleakala Volcano's northeast flank, Maui

We present a morphostructural analysis of the submarine portions of Haleakala Volcano and environs, based upon a 4-year program of geophysical surveys and submersible explorations of the underwater flanks of Hawaiian volcanoes that was conducted by numerous academic and governmental research organizations in Japan and the U.S. and funded primarily by the Japan Agency for Marine–Earth Science and Technology. A resulting reconnaissance geologic map features the 135-km-long Hana Ridge, the 3000 km 2 Hana slump on the volcano's northeast flank, and island-surrounding terraces that are the submerged parts of volcanic shields. Hana Ridge below 2000 m water depth exhibits the lobate morphology typical of the subaqueously erupted parts of Hawaiian rift zones, with some important distinctions: namely, subparallel crestlines, which we propose result from the down-rift migration of offsets in the dike intrusion zone, and an amphitheater at its distal toe, where a submarine landslide has embayed the ridge tip. Deformation of Haleakala's northeast flank is limited to that part identified as the Hana slump, which lies downslope from the volcano's submerged shield, indicating that flank mobility is also limited in plan, inconsistent with hypothesized volcanic spreading driven by rift-zone dilation. The leading edge of the slump has transverse basins and ridges that resemble the thrust ramps of accretionary prisms, and we present a model to describe the slump's development that emphasizes the role of coastally generated fragmental basalt on gravitational instability of Haleakala's northeast flank and that may be broadly applicable to other ocean-island slumps.

Hawai'i↗

Attenuation of barium, strontium, cobalt, and nickel plumes formed during microbial iron-reduction in a crude-oil-contaminated aquifer

We assessed the spatial distribution of 35 elements in aquifer sediments and groundwater of a crude-oil-contaminated aquifer and show evidence of the dissolution of barium (Ba), strontium (Sr), cobalt (Co), and nickel (Ni) during hydrocarbon oxidation coupled to historic microbial Fe(III)-reduction near the oil. Trace element plumes occur in the crude-oil-contaminated aquifer, where 50% Co, 47% Ni, 24% Ba, and 15% Sr have been mobilized from the sediment near the oil into groundwater, resulting in dissolved masses >33, 18, three, and two times greater than estimated dissolved masses prior to contamination, respectively. Ba 2+ and Ni 2+ concentrations exceeded the World Health Organization’s drinking-water guidelines of 700 and 20 μg/L, respectively. Sediments attenuate trace element plumes in two geochemically distinct zones, resulting in <0.01% total trace element masses dissolved in groundwater, despite the substantial mobilization near the oil body. Geochemical modeling of the modern Fe(III)-reducing zone suggests trace elements are likely attenuated via coprecipitation with/without sorption on iron carbonate precipitates. In the suboxic transition zone at the leading edge of the plume, Fe(III)-hydroxides sorb Ba 2+ , Sr 2+ , Co 2+ , and Ni 2+ . This study emphasizes that slow but persistent biogeochemical activity can substantially alter aquifer chemistry over decadal timeframes, a phenomenon we term biogeochemical gradualism.

Minnesota↗

Transport and degradation of semivolatile hydrocarbons in a petroleum-contaminated aquifer, Bemidji, Minnesota

Polycyclic aromatic hydrocarbons (PAH) were used as probes to identify the processes controlling the transport and fate of aqueous semivolatile hydrocarbons (SVHCs) in a petroleum-contaminated aquifer near Bemidji, Minnesota. PAH and other SVHCs were isolated from ground water by field solid-phase extraction and analyzed using gas chromatography/mass spectrometry. Close to the oil body, aqueous aliphatic hydrocarbon compositions are substantially different from the parent oil, suggesting microbial alteration prior to or during dissolution. Aqueous PAH concentrations are elevated above oil-water equilibrium concentrations directly beneath the oil and decrease dramatically at distances ranging from the 25 to 65 m downgradient from the leading edge of the oil body. Variations in downgradient distributions of naphthalene, fluorene and phenanthrene, coupled with their biodegradation, partitioning and volatility characteristics, suggest that the PAH are useful probes for distinguishing between the biogeochemical processes affecting SVHC transport and persistence in ground water.

Minnesota↗

Modeling the influence of variable pH on the transport of zinc in a contaminated aquifer using semiempirical surface complexation models

Land disposal of sewage effluent resulted in contamination of a sand and gravel aquifer (Cape Cod, Massachusetts) with zinc (Zn). The distribution of Zn was controlled by pH‐dependent adsorption; the Zn extended 15 m into the 30‐m‐thick sewage plume within approximately 100 m of the source but only 2–4 m into the plume between 100 and 400 m downgradient. A two‐dimensional vertical cross section model coupling groundwater flow with solute transport and equilibrium adsorption is used to simulate the influence of pH on Zn transport. Adsorption is described using semiempirical surface complexation models (SCM) by writing chemical reactions between dissolved Zn and mineral surface sites. SCM parameters were determined in independent laboratory experiments. A 59‐year simulation with a one‐site SCM describes the influence of pH on Zn transport well, with greater mobility at the low pH values near the upper sewage plume boundary than at the higher pH values deeper in the sewage‐contaminated zone. Simulation with a two‐site SCM describes both the sharpness and approximate location of the leading edge of the Zn‐contaminated region. Temporal variations in pH of incoming groundwater can result in large increases in Zn concentration and mobility. The influence of spatial and temporal variability in pH on adsorption and transport of Zn was accomplished much more easily with the semiempirical SCM approach than could be achieved with distribution coefficients or adsorption isotherms.

Water Resources Research↗

Source materials for inception stage Hawaiian magmas: Pb‐He isotope variations for early Kilauea

New noble gas and radiogenic isotopic compositions are presented for tholeiitic, transitional, and alkalic rocks from the submarine Hilina region on the south flank of Kilauea, Hawaii. The 3 He/ 4 He ratios for undegassed glass and olivine separates (11–26 Ra) contrast with those of postshield and rejuvenated alkalic lavas, consistent with the alkalic and transitional basalts at Hilina corresponding to early Kilauea magmas. Most early Kilauea samples contain highly radiogenic Pb isotopes compared with other Hawaiian rocks and therefore derive from a Hawaiian plume end‐member source (here referred to as the Hilina component) distinctive in that respect. Besides radiogenic Pb isotopes, the Hilina component has relatively low 3 He/ 4 He (<12 Ra) among the Hawaiian magmas. Hawaiian inception stage magmas, including Hilina, Loihi, and deep Hana Ridge (east Maui), define a linear array in 206 Pb/ 204 Pb‐ 3 He/ 4 He isotope space, indicating that mixing between the Hilina and Loihi components (or their melts) dominates magmatism at the leading edge of the Hawaiian plume. The Hilina component's isotopic characteristics can be derived from young subduction‐recycled crust or metasomatised mantle. The isotopic differences between the geographically discriminated Kea and Loa trend volcanic chains, observed in shield stage lavas, are also seen in the inception stage magmas, suggesting that proportions of melts derived from the Hilina and Loihi components were different between the Kea and Loa trend volcanoes.

Hawaii↗

Unsaturated flow processes and the onset of seasonal deformation in slow-moving landslides

Predicting rainfall-induced landslide motion is challenging because shallow groundwater flow is extremely sensitive to the preexisting moisture content in the ground. Here, we use groundwater hydrology theory and numerical modeling combined with five years of field monitoring to illustrate how unsaturated groundwater flow processes modulate the seasonal pore water pressure rise and therefore the onset of motion for slow-moving landslides. The onset of landslide motion at Oak Ridge earthflow in California’s Diablo Range occurs after an abrupt water table rise to near the landslide surface 52–129 days after seasonal rainfall commences. Model results and theory suggest that this abrupt rise occurs from the advection of a nearly saturated wetting front, which marks the leading edge of the integrated downward flux of seasonal rainfall, to the water table. Prior to this abrupt rise, we observe little measured pore water pressure response within the landslide due to rainfall. However, once the wetting front reaches the water table, we observe nearly instantaneous pore water pressure transmission within the landslide body that is accompanied by landslide acceleration. We cast the timescale to reach a critical pore water pressure threshold using a simple mass balance model that considers variable moisture storage with depth and explains the onset of seasonal landslide motion with a rainfall intensity-duration threshold. Our model shows that the seasonal response time of slow-moving landslides is controlled by the dry season vadose zone depth rather than the total landslide thickness.

California↗