Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Interfaces”

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 415 records · Page 23Linked to original sources

Seismic gaps and source zones of recent large earthquakes in coastal Peru

The earthquakes of central coastal Peru occur principally in two distinct zones of shallow earthquake activity that are inland of and parallel to the axis of the Peru Trench. The interface-thrust (IT) zone includes the great thrust-fault earthquakes of 17 October 1966 and 3 October 1974. The coastal-plate interior (CPI) zone includes the great earthquake of 31 May 1970, and is located about 50 km inland of and 30 km deeper than the interface thrust zone. The occurrence of a large earthquake in one zone may not relieve elastic strain in the adjoining zone, thus complicating the application of the seismic gap concept to central coastal Peru. However, recognition of two seismic zones may facilitate detection of seismicity precursory to a large earthquake in a given zone; removal of probable CPI-zone earthquakes from plots of seismicity prior to the 1974 main shock dramatically emphasizes the high seismic activity near the rupture zone of that earthquake in the five years preceding the main shock. Other conclusions on the seismicity of coastal Peru that affect the application of the seismic gap concept to this region are: (1) Aftershocks of the great earthquakes of 1966, 1970, and 1974 occurred in spatially separated clusters. Some clusters may represent distinct small source regions triggered by the main shock rather than delimiting the total extent of main-shock rupture. The uncertainty in the interpretation of aftershock clusters results in corresponding uncertainties in estimates of stress drop and estimates of the dimensions of the seismic gap that has been filled by a major earthquake. (2) Aftershocks of the great thrust-fault earthquakes of 1966 and 1974 generally did not extend seaward as far as the Peru Trench. (3) None of the three great earthquakes produced significant teleseismic activity in the following month in the source regions of the other two earthquakes. The earthquake hypocenters that form the basis of this study were relocated using station adjustments computed by the method of joint hypocenter determination. ?? 1979 Birkha??user Verlag.

Pure and Applied Geophysics PAGEOPH↗

Magmatic inclusions in rhyolites, contaminated basalts, and compositional zonation beneath the Coso volcanic field, California

Basaltic lava flows and high-silica rhyolite domes form the Pleistocene part of the Coso volcanic field in southeastern California. The distribution of vents maps the areal zonation inferred for the upper parts of the Coso magmatic system. Subalkalic basalts (<50% SiO 2 ) were erupted well away from the rhyolite field at any given time. Compositional variation among these basalts can be ascribed to crystal fractionation. Erupted volumes of these basalts decrease with increasing differentiation. Mafic lavas containing up to 58% SiO 2 , erupted adjacent to the rhyolite field, formed by mixing of basaltic and silicic magma. Basaltic magma interacted with crustal rocks to form other SiO 2 -rich mafic lavas erupted near the Sierra Nevada fault zone. Several rhyolite domes in the Coso volcanic field contain sparse andesitic inclusions (55–61% SiO 2 ). Pillow-like forms, intricate commingling and local diffusive mixing of andesite and rhyolite at contacts, concentric vesicle distribution, and crystal morphologies indicative of undercooling show that inclusions were incorporated in their rhyolitic hosts as blobs of magma. Inclusions were probably dispersed throughout small volumes of rhyolitic magma by convective (mechanical) mixing. Inclusion magma was formed by mixing (hybridization) at the interface between basaltic and rhyolitic magmas that coexisted in vertically zoned igneous systems. Relict phenocrysts and the bulk compositions of inclusions suggest that silicic endmembers were less differentiated than erupted high-silica rhyolite. Changes in inferred endmembers of magma mixtures with time suggest that the steepness of chemical gradients near the silicic/mafic interface in the zoned reservoir may have decreased as the system matured, although a high-silica rhyolitic cap persisted. The Coso example is an extreme case of large thermal and compositional contrast between inclusion and host magmas; lesser differences between intermediate composition magmas and inclusions lead to undercooling phenomena that suggest smaller ΔT. Vertical compositional zonation in magma chambers has been documented through study of products of voluminous pyroclastic eruptions. Magmatic inclusions in volcanic rocks provide evidence for compositional zonation and mixing processes in igneous systems when only lava is erupted.

California↗

Lead and strontium isotopic evidence for crustal interaction and compositional zonation in the source regions of Pleistocene basaltic and rhyolitic magmas of the Coso volcanic field, California

The isotopic compositions of Pb and Sr in Pleistocene basalt, high-silica rhyolite, and andesitic inclusions in rhyolite of the Coso volcanic field indicate that these rocks were derived from different levels of compositionally zoned magmatic systems. The 2 earliest rhyolites probably were tapped from short-lived silicic reservoirs, in contrast to the other 36 rhyolite domes and lava flows which the isotopic data suggest may have been leaked from the top of a single, long-lived magmatic system. Most Coso basalts show isotopic, geochemical, and mineralogic evidence of interaction with crustal rocks, but one analyzed flow has isotopic ratios that may represent mantle values (87Sr/86Sr=0.7036,206Pb/204Pb=19.05,207Pb/204Pb=15.62,208Pb/204Pb= 38.63). The (initial) isotopic composition of typical rhyolite (87Sr/86Sr=0.7053,206Pb/204Pb=19.29,207Pb/204Pb= 15.68,208Pb/204Pb=39.00) is representative of the middle or upper crust. Andesitic inclusions in the rhyolites are evidently samples of hybrid magmas from the silicic/mafic interface in vertically zoned magma reservoirs. Silicic end-member compositions inferred for these mixed magmas, however, are not those of erupted rhyolite but reflect the zonation within the silicic part of the magma reservoir. The compositional contrast at the interface between mafic and silicic parts of these systems apparently was greater for the earlier, smaller reservoirs.

Contributions to Mineralogy and Petrology↗

Surface electromagnetic geophysical exploration of the ground-water resources of Isla de Mona, Puerto Rico, a caribbean carbonate island

Two electromagnetic surface geophysical techniques were used to explore the ground-water resources of Isla de Mona, Puerto Rico&mdash;a 55-square-kilometer island located between Puerto Rico and Hispa&ntilde;iola, Isla de Mona is a tectonically uplifted carbonate plateau of Neogene age that has an average elevation of about 50 meters above mean sea level. This plateau is bounded by vertical cliffs except on the southwest where there is a narrow, 3-square-kilometer coastal plain. The coastal plain is composed of Quaternary carbonate deposits, and has a maximum elevation of less than 10 meters above mean sea level. No large-scale surface-water features are found on the plateau or on the coastal plain. To better understand the aquifer characteristics of the island, terrain conductivity and transient electromagnetic data were collected on the coastal plain and on the plateau. Computer programs were used to analyze quantitatively the electromagnetic data. Geoelectric models were produced to approximate the depth below land surface of the saline-freshwater interface underlying both the coastal plain and the plateau. Because the geophysical methods could not discern the water-table, it was assumed that it was near sea level. The thickness of the freshwater lens was estimated by subtracting the elevation of the land surface above mean sea level from the depth to the saline-freshwater interface as determined from the geophysical data. Results from the geophysical methods and water-level observations indicate that a freshwater lens with a maximum thickness of about 10 meters exists under the coastal plain. This lens thins towards the ocean, and it also thins away from the ocean toward the plateau. The model produced from the transient electromagnetic data indicates that the freshwater lens under the plateau has a maximum thickness of about 14 meters, which is a much thinner lens than previously estimated. A freshwater lens thickness of 14 meters is similar to direct measurements by divers in Cueva de Agua at Punta los Ingleses (located on the southeast coast), and observations by divers of freshwater seeps into the ocean at depths of 8 to 10 meters along the north coast cliffs after a major rainstorm. Ground-water flow in the coastal plain appears to be radial from the center of the freshwater mound. At the intersection between the coastal plain and the plateau, the flow is parallel to the coastline. The direction of flow on the rest of the plateau could not be determined accurately with the available data.

Carbonates and Evaporites↗

Microbial communities in standing dead trees in ghost forests are largely aerobic, saprophytic, and methanotrophic

Standing dead trees (snags) are recognized for their influence on methane (CH 4 ) cycling in coastal wetlands, yet the biogeochemical processes that control the magnitude and direction of fluxes across the snag-atmosphere interface are not fully elucidated. Herein, we analyzed microbial communities and fluxes at one height from ten snags in a ghost forest wetland. Snag-atmosphere CH 4 fluxes were highly variable (− 0.11–0.51 mg CH 4 m −2 h −1 ). CH 4 production was measured in three out of ten snags; whereas, CH 4 consumption was measured in two out of ten snags. Potential CH 4 production and oxidation in one core from each snag was assayed in vitro. A single core produced CH 4 under anoxic and oxic conditions, at measured rates of 0.7 and 0.6 ng CH 4 g −1 h −1 , respectively. Four cores oxidized CH 4 under oxic conditions, with an average rate of − 1.13 ± 0.31 ng CH 4 g −1 h −1 . Illumina sequencing of the V3/V4 region of the 16S rRNA gene sequence revealed diverse microbial communities and indicated oxidative decomposition of deadwood. Methanogens were present in 20% of the snags, with a mean relative abundance of < 0.0001%. Methanotrophs were identified in all snags, with a mean relative abundance of 2% and represented the sole CH 4 -cycling communities in 80% of the snags. These data indicate potential for microbial attenuation of CH 4 emissions across the snag-atmosphere interface in ghost forests. A better understanding of the environmental drivers of snag-associated microbial communities is necessary to forecast the response of CH 4 cycling in coastal ghost forest wetlands to a shifting coastal landscape.

Current Microbiology↗

Two mechanisms of aquatic and terrestrial habitat change along an Alaskan Arctic coastline

Arctic habitats at the interface between land and sea are particularly vulnerable to climate change. The northern Teshekpuk Lake Special Area (N-TLSA), a coastal plain ecosystem along the Beaufort Sea in northern Alaska, provides habitat for migratory waterbirds, caribou, and potentially, denning polar bears. The 60-km coastline of N-TLSA is experiencing increasing rates of coastline erosion and storm surge flooding far inland resulting in lake drainage and conversion of freshwater lakes to estuaries. These physical mechanisms are affecting upland tundra as well. To better understand how these processes are affecting habitat, we analyzed long-term observational records coupled with recent short-term monitoring. Nearly the entire coastline has accelerating rates of erosion ranging from 6 m/year from 1955 to 1979 and most recently peaking at 17 m/year from 2007 to 2009, yet an intensive monitoring site along a higher bluff (3–6 masl) suggested high interannual variability. The frequency and magnitude of storm events appears to be increasing along this coastline and these patterns correspond to a greater number of lake tapping and flooding events since 2000. For the entire N-TLSA, we estimate that 6% of the landscape consists of salt-burned tundra, while 41% is prone to storm surge flooding. This offset may indicate the relative frequency of low-magnitude flood events along the coastal fringe. Monitoring of coastline lakes confirms that moderate westerly storms create extensive flooding, while easterly storms have negligible effects on lakes and low-lying tundra. This study of two interacting physical mechanisms, coastal erosion and storm surge flooding, provides an important example of the complexities and data needs for predicting habitat change and biological responses along Arctic land–ocean interfaces.

Alaska↗

Seismo-acoustic signals associated with degassing explosions recorded at Shishaldin Volcano, Alaska, 2003-2004

In summer 2003, a Chaparral Model 2 microphone was deployed at Shishaldin Volcano, Aleutian Islands, Alaska. The pressure sensor was co-located with a short-period seismometer on the volcano&rsquo;s north flank at a distance of 6.62 km from the active summit vent. The seismo-acoustic data exhibit a correlation between impulsive acoustic signals (1&ndash;2 Pa) and long-period (LP, 1&ndash;2 Hz) earthquakes. Since it last erupted in 1999, Shishaldin has been characterized by sustained seismicity consisting of many hundreds to two thousand LP events per day. The activity is accompanied by up to &sim;200 m high discrete gas puffs exiting the small summit vent, but no significant eruptive activity has been confirmed. The acoustic waveforms possess similarity throughout the data set (July 2003&ndash;November 2004) indicating a repetitive source mechanism. The simplicity of the acoustic waveforms, the impulsive onsets with relatively short (&sim;10&ndash;20 s) gradually decaying codas and the waveform similarities suggest that the acoustic pulses are generated at the fluid&ndash;air interface within an open-vent system. SO 2 measurements have revealed a low SO 2 flux, suggesting a hydrothermal system with magmatic gases leaking through. This hypothesis is supported by the steady-state nature of Shishaldin&rsquo;s volcanic system since 1999. Time delays between the seismic LP and infrasound onsets were acquired from a representative day of seismo-acoustic data. A simple model was used to estimate source depths. The short seismo-acoustic delay times have revealed that the seismic and acoustic sources are co-located at a depth of 240&plusmn;200 m below the crater rim. This shallow depth is confirmed by resonance of the upper portion of the open conduit, which produces standing waves with f =0.3 Hz in the acoustic waveform codas. The infrasound data has allowed us to relate Shishaldin&rsquo;s LP earthquakes to degassing explosions, created by gas volume ruptures from a fluid&ndash;air interface.

Alaska↗

Pyroclast/snow interactions and thermally driven slurry formation. Part 1: Theory for monodisperse grain beds

Lahars are often produced as pyroclastic flows move over snow. This phenomenon involves a complicated interplay of mechanical and thermal processes that need to be separated to get at the fundamental physics. The thermal physics of pyroclast/snow interactions form the focus of this paper. A theoretical model is developed of heat- and mass transfer at the interface between a layer of uniformly sized pyroclasts and an underlying bed of snow, for the case in which there is no relative shear motion between pyroclasts and snow. A microscale view of the interface is required to properly specify boundary conditions. The physical model leads to the prediction that the upward flux of water vapor - which depends upon emplacement temperature, pyroclast grain size, pyroclast-layer thickness, and snow permeability - is sometimes sufficient to fluidize the pyroclasts. Uniform fluidization is usually unstable to bubble formation, which leads to vigorous convection of the pyroclasts themselves. Thus, predicted threshold conditions for fluidization are tantamount to predicted thresholds for particle convection. Such predictions are quantitatively in good agreement with results of experiments described in part 2 of this paper. Because particle convection commonly causes scour of the snow bed and transformation of the pyroclast layer to a slurry, there exists a 'thermal scour' process for generating lahars from pyroclastic flows moving over snow regardless of the possible role of mechanical scour.

Bulletin of Volcanology↗

Three-dimensional variable-density flow simulation of a coastal aquifer in southern Oahu, Hawaii, USA

Three-dimensional modeling of groundwater flow and solute transport in the Pearl Harbor aquifer, southern Oahu, Hawaii, shows that the readjustment of the freshwater–saltwater transition zone takes a long time following changes in pumping, irrigation, or recharge in the aquifer system. It takes about 50 years for the transition zone to move 90% of the distance to its new steady position. Further, the Ghyben–Herzberg estimate of the freshwater/saltwater interface depth occurred between the 10 and 50% simulated seawater concentration contours in a complex manner during 100 years of the pumping history of the aquifer. Thus, it is not a good predictor of the depth of potable water. Pre-development recharge was used to simulate the 1880 freshwater-lens configuration. Historical pumpage and recharge distributions were used and the resulting freshwater-lens size and position were simulated through 1980. Simulations show that the transition zone moved upward and landward during the period simulated. Previous groundwater flow models for Oahu have been limited to areal models that simulate a sharp interface between freshwater and saltwater or solute-transport models that simulate a vertical aquifer section. The present model is based on the US Geological Survey’s three-dimensional solute transport (3D SUTRA) computer code. Using several new tools for pre- and post-processing of model input and results have allowed easy model construction and unprecedented visualization of the freshwater lens and underlying transition zone in Hawaii’s most developed aquifer.

Hawaii↗

In situ allelopathic expression by the invasive amphibious plant, Ludwigia hexapetala (water primrose) across habitat types, seasons, and salinities

Broad infestations of invasive, non-native vegetation have transformed wetlands around the world. Ludwigia hexapetala is a widespread, amphibious invasive plant with a creeping growth habit in open water and an erect growth habit in terrestrial habitats. In the upper San Francisco Estuary of California, L. hexapetala is increasingly terrestrializing into marshes and this expansion may be facilitated by allelopathy. We conducted the first field-based study on L. hexapetala allelopathy to determine whether (1) three allelochemicals known to be exuded by L. hexapetala are expressed in situ, (2) the allelochemicals are detectable in leaves, soil, and water, and (3) allelopathic expression varies by season, salinity, and growth habit (open water “patch” vs. terrestrial marsh “interface” locations). Water, soil, and L. hexapetala leaves were collected in two freshwater sites and two oligohaline sites in the upper San Francisco Estuary in summer 2021, fall 2021, and spring 2022. Myricitrin and quercitrin, known allelochemicals, and salipurposid, a newly identified polyphenol, were detected in water, soil, and leaves. There were significant differences in allelochemical concentrations under fresh versus oligohaline conditions in water and soil, but not leaves. All three allelochemicals generally had higher concentrations in patch versus interface locations, suggesting that L. hexapetala allelopathy plays a greater competitive role in open water than terrestrial habitats. Leaf concentrations of each allelochemical varied seasonally; however, both myricitrin and salipurposid had heightened concentrations in spring. These results suggest that herbicide application in early spring may be most effective in controlling L. hexapetala terrestrialization from open water to marshes.

California↗

Influence of groundwater recharge and well characteristics on dissolved arsenic concentrations in southeastern Michigan groundwater

Arsenic concentrations exceeding 10 ??g/l, the United States maximum contaminant level and the World Health Organization guideline value, are frequently reported in groundwater from bedrock and unconsolidated aquifers of southeastern Michigan. Although arsenic-bearing minerals (including arsenian pyrite and oxide/hydroxide phases) have been identified in Marshall Sandstone bedrock of the Mississippian aquifer system and in tills of the unconsolidated aquifer system, mechanisms responsible for arsenic mobilization and subsequent transport in groundwater are equivocal. Recent evidence has begun to suggest that groundwater recharge and characteristics of well construction may affect arsenic mobilization and transport. Therefore, we investigated the relationship between dissolved arsenic concentrations, reported groundwater recharge rates, well construction characteristics, and geology in unconsolidated and bedrock aquifers. Results of multiple linear regression analyses indicate that arsenic contamination is more prevalent in bedrock wells that are cased in proximity to the bedrock-unconsolidated interface; no other factors were associated with arsenic contamination in water drawn from bedrock or unconsolidated aquifers. Conditions appropriate for arsenic mobilization may be found along the bedrock-unconsolidated interface, including changes in reduction/oxidation potential and enhanced biogeochemical activity because of differences between geologic strata. These results are valuable for understanding arsenic mobilization and guiding well construction practices in southeastern Michigan, and may also provide insights for other regions faced with groundwater arsenic contamination. ?? Springer-Verlag 2008.

Environmental Geochemistry and Health↗

Riders on the storm: selective tidal movements facilitate the spawning migration of threatened delta smelt in the San Francisco Estuary

Migration strategies in estuarine fishes typically include behavioral adaptations for reducing energetic costs and mortality during travel to optimize reproductive success. The influence of tidal currents and water turbidity on individual movement behavior were investigated during the spawning migration of the threatened delta smelt, Hypomesus transpacificus , in the northern San Francisco Estuary, California, USA. Water current velocities and turbidity levels were measured concurrently with delta smelt occurrence at sites in the lower Sacramento River and San Joaquin River as turbidity increased due to first-flush winter rainstorms in January and December 2010. The presence/absence of fish at the shoal-channel interface and near the shoreline was quantified hourly over complete tidal cycles. Delta smelt were caught consistently at the shoal-channel interface during flood tides and near the shoreline during ebb tides in the turbid Sacramento River, but were rare in the clearer San Joaquin River. The apparent selective tidal movements by delta smelt would facilitate either maintaining position or moving upriver on flood tides, and minimizing advection down-estuary on ebb tides. These movements also may reflect responses to lateral gradients in water turbidity created by temporal lags in tidal velocities between the near-shore and mid-channel habitats. This migration strategy can minimize the energy spent swimming against strong river and tidal currents, as well as predation risks by remaining in turbid water. Selection pressure on individuals to remain in turbid water may underlie population-level observations suggesting that turbidity is a key habitat feature and cue initiating the delta smelt spawning migration.

California↗

Rates of dissolution of aluminosilicates in seawater

Dissolution of eight clay minerals, four zeolites, and quartz in seawater has been monitored for 8 1 2 years. For most of the minerals, dissolution can be described as a first-order reaction in which dissolved silica approaches from undersaturation steady concentration values with time. Characteristic reaction rate constants (k 1 ) are of the order of 10 -7 sec -1 . One of the zeolites, clinoptilolite, shows a different dissolution behavior: SiO 2 concentration in solution reaches a high value within one year, followed by a decline to a lower value, suggestive of precipitation of another silicate phase (possibly sepiolite). A mathematical solution is given for a kinetic equation combining the parabolic-rate and first-order rate processes. It is shown that in a wide range of silicate dissolution reactions taking place over long periods of time, the presence of the parabolic-rate dissolution processes cannot be detected, thereby making its inclusion in the kinetic equations unnecessary. The experimental rates of dissolution are comparable to the SiO 2 - dissolution rates in oceanic sediments near the sediment/water interface. But deeper in the sediment, the calculated dissolution rates are significantly lower than the near-interface and experimental values. © 1975.

Earth and Planetary Science Letters↗

Biogeochemical cycling in an organic-rich coastal marine basin. 7. Sulfur mass balance, oxygen uptake and sulfide retention

Sulfur and oxygen fluxes were quantified in the seasonally varying anoxic marine sedimentary system of Cape Lookout Bight, N.C., U.S.A. Over the three year study period, 1981–1983, the mean annual sulfate reduction rate was determined to be 18.2 ± 1.6 moles · m −2 · y −1 . This value, added to the estimate of the detrital sulfur input of 1.2 ± 4.4 gave a total sulfur input of 19.4 ± 4.7 moles · m −2 · y −1 . The sulfide flux to the sediment-water interface, measured in anaerobic benthic chambers was 4.6 ± 0.5 moles · m −2 · y −1 , and represented 37% of the annual oxygen uptake rate of 25.2 ± 2.8 moles · m −2 · y −1 . The sulfide burial rate, determined to be 15.5 ± 3.1 moles · m −2 · y −1 , was within 5% of the value predicted by summing the fluxes at the sediment-water interface. The C S "> CS weight ratio of sediment below the depth of diagenetic reaction was determined to be 2.75. The sulfide retention rate in these rapidly accumulating sediments (10 cm/yr) was 77 ± 19%. Comparison of this result with previous studies shows that rapid sediment accumulation and the lack of bioturbation control this unusually high degree of sulfide retention.

Geochimica et Cosmochimica Acta↗

Organic geochemistry and brine composition in Great Salt, Mono, and Walker Lakes

Samples of Recent sediments, representing up to 1000 years of accumulation, were collected from three closed basin lakes (Mono Lake, CA, Walker Lake, NV, and Great Salt Lake, UT) to assess the effects of brine composition on the accumulation of total organic carbon, the concentration of dissolved organic carbon, humic acid structure and diagenesis, and trace metal complexation. The Great Salt Lake water column is a stratified Na-Mg-Cl-SO 4 brine with low alkalinity. Algal debris is entrained in the high density (1.132–1.190 g/cc) bottom brines, and in this region maximum organic matter decomposition occurs by anaerobic processes, with sulfate ion as the terminal electron acceptor. Organic matter, below 5 cm of the sediment-water interface, degrades at a very slow rate in spite of very high pore-fluid sulfate levels. The organic carbon concentration stabilizes at 1.1 wt%. Mono Lake is an alkaline (Na-CO 3 -Cl-SO 4 ) system. The water column is stratified, but the bottom brines are of lower density relative to the Great Salt Lake, and sedimentation of algal debris is rapid. Depletion of pore-fluid sulfate, near l m of core, results in a much higher accumulation of organic carbon, approximately 6 wt%. Walker Lake is also an alkaline system. The water column is not stratified, and decomposition of organic matter occurs by aerobic processes at the sediment-water interface and by anaerobic processes below. Total organic carbon and dissolved organic carbon concentrations in Walker Lake sediments vary with location and depth due to changes in input and pore-fluid sulfate concentrations. Nuclear magnetic resonance studies ( 13 C) of humic substances and dissolved organic carbon provide information on the source of the Recent sedimentary organic carbon (aquatic vs. terrestrial), its relative state of decomposition, and its chemical structure. The spectra suggest an algal origin with little terrestrial signature at all three lakes. This is indicated by the ratio of aliphatic to aromatic carbon and the absence of chemical structures indicative of the lignin of vascular plants. The dissolved organic carbon of the Mono Lake pore fluids is structurally related to humic acid and is also related to carbohydrate metabolism. The alkaline pore fluids, due to high pH, solubilize high molecular weight organic matter from the sediments. This hydrophilic material is a metal complexing agent. Despite very high algal productivities, organic carbon accumulation can be low in stratified lakes if the anoxic bottom waters are hypersaline with high concentrations of sulfate ion. Labile organic matter is recycled to the water column and the sedimentary organic matter is relatively nonsusceptible to bacterial metabolism. As a result, pore-fluid dissolved organic carbon and metal-organic complexation are low.

Geochimica et Cosmochimica Acta↗

Crystallization of accessory phases in magmas by local saturation adjacent to phenocrysts

Accessory minerals commonly occur attached to or included in the major crystalline phases of felsic and some intermediate igneous rocks. Apatite is particularly common as inclusions, but Fe-Ti oxides, pyrrhotite, zircon, monazite, chevkinite and xenotime are also known from silicic rocks. Accessories may nucleate near the host crystal/ liquid interface as a result of local saturation owing to formation of a differentiated chemical boundary layer in which accessory mineral solubility would be lower than in the surrounding liquid. Differentiation of this boundary layer would be greatest adjacent to ferromagnesian phenocrysts, especially Fe-Ti oxides; it is with oxides that accessories are most commonly associated in rocks. A boundary layer may develop if the crystal grows more rapidly than diffusion can transport incorporated and rejected elements to and from the phenocryst. Diffusion must dominate over convection as a mode of mass transfer near the advancing crystal/liquid interface in order for a boundary layer to exist. Accumulation of essential structural constituent elements of accessory minerals owing to their slow diffusion in evolved silicate melt also may force local saturation, but this is not a process that applies to all cases. Local saturation is an attractive mechanism for enhancing fractionation during crystallization differentiation. If accessory minerals attached to or included in phenocrysts formed because of local saturation, their host phenocrysts must have grown rapidly when accessories nucleated in comparison to lifetimes of magma reservoirs. Some inconsistencies remain in a local saturation origin for accessory phases that cannot be evaluated without additional information.

Geochimica et Cosmochimica Acta↗

Salt-water-freshwater transient upconing - An implicit boundary-element solution

The boundary-element method is used to solve the set of partial differential equations describing the flow of salt water and fresh water separated by a sharp interface in the vertical plane. In order to improve the accuracy and stability of the numerical solution, a new implicit scheme was developed for calculating the motion of the interface. The performance of this scheme was tested by means of numerical simulation. The numerical results are compared to experimental results for a salt-water upconing under a drain problem.

Journal of Hydrology↗

Theory and application of an approximate model of saltwater upconing in aquifers

Motion and mixing of salt water and fresh water are vitally important for water-resource development throughout the world. An approximate model of saltwater upconing in aquifers is developed, which results in three non-linear coupled equations for the freshwater zone, the saltwater zone, and the transition zone. The description of the transition zone uses the concept of a boundary layer. This model invokes some assumptions to give a reasonably tractable model, considerably better than the sharp interface approximation but considerably simpler than a fully three-dimensional model with variable density. We assume the validity of the Dupuit-Forchheimer approximation of horizontal flow in each layer. Vertical hydrodynamic dispersion into the base of the transition zone is assumed and concentration of the saltwater zone is assumed constant. Solute in the transition zone is assumed to be moved by advection only. Velocity and concentration are allowed to vary vertically in the transition zone by using shape functions. Several numerical techniques can be used to solve the model equations, and simple analytical solutions can be useful in validating the numerical solution procedures. We find that the model equations can be solved with adequate accuracy using the procedures presented. The approximate model is applied to the Smoky Hill River valley in central Kansas. This model can reproduce earlier sharp interface results as well as evaluate the importance of hydrodynamic dispersion for feeding salt water to the river. We use a wide range of dispersivity values and find that unstable upconing always occurs. Therefore, in this case, hydrodynamic dispersion is not the only mechanism feeding salt water to the river. Calculations imply that unstable upconing and hydrodynamic dispersion could be equally important in transporting salt water. For example, if groundwater flux to the Smoky Hill River were only about 40% of its expected value, stable upconing could exist where hydrodynamic dispersion into a transition zone is the primary mechanism for moving salt water to the river. The current model could be useful in situations involving dense saltwater layers.

Journal of Hydrology↗