Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Processes”

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 919 records · Page 51Linked to original sources
USGS2015

Evaluating the importance of characterizing soil structure and horizons in parameterizing a hydrologic process model

Incorporating the influence of soil structure and horizons into parameterizations of distributed surface water/groundwater models remains a challenge. Often, only a single soil unit is employed, and soil-hydraulic properties are assigned based on textural classification, without evaluating the potential impact of these simplifications. This study uses a distributed physics-based model to assess the influence of soil horizons and structure on effective parameterization. This paper tests the viability of two established and widely used hydrogeologic methods for simulating runoff and variably saturated flow through layered soils: (1) accounting for vertical heterogeneity by combining hydrostratigraphic units with contrasting hydraulic properties into homogeneous, anisotropic units and (2) use of established pedotransfer functions based on soil texture alone to estimate water retention and conductivity, without accounting for the influence of pedon structures and hysteresis. The viability of this latter method for capturing the seasonal transition from runoff-dominated to evapotranspiration-dominated regimes is also tested here. For cases tested here, event-based simulations using simplified vertical heterogeneity did not capture the state-dependent anisotropy and complex combinations of runoff generation mechanisms resulting from permeability contrasts in layered hillslopes with complex topography. Continuous simulations using pedotransfer functions that do not account for the influence of soil structure and hysteresis generally over-predicted runoff, leading to propagation of substantial water balance errors. Analysis suggests that identifying a dominant hydropedological unit provides the most acceptable simplification of subsurface layering and that modified pedotransfer functions with steeper soil-water retention curves might adequately capture the influence of soil structure and hysteresis on hydrologic response in headwater catchments.

Hydrological Processes↗
USGS2001

A process for fire-related debris flow initiation, Cerro Grande fire, New Mexico

In this study we examine factors that pertain to the generation of debris flows from a basin recently burned by wildfire. Throughout the summer 2000 thunderstorm season, we monitored rain gauges, channel cross-sections, hillslope transects, and nine sediment-runoff traps deployed in a steep, 0??15 km2 basin burned by the May 2000 Cerro Grande fire in New Mexico. Debris flows were triggered in the monitored basin during a rainstorm on July 16, 2000, in response to a maximum 30 min rainfall intensity of 31 mm h-1 (return period of approximately 2 years). Eleven other storms occurred before and after the July storm; these storms resulted in significant runoff, but did not generate debris flows. The debris flows generated by the July 16 storm initiated on a broad, open hillslope as levee-lined rills. The levees were composed of gravel- and cobble-sized material supported by an abundant fine-grained matrix. Debris-flow deposits were observed only on the hillslopes and in the first and second-order drainages of the monitored basin. No significant amounts of channel incision were measured following the passage of the debris flows, indicating that most of the material in the flows originated from the hillslopes. Sediment-runoff concentrations of between 0??23 and 0??81 kg 1-1 (with a mean of 0-42 kg 1-1) were measured from the hillslope traps following the debris-flow-producing storm. These concentrations, however, were not unique to the July 16 storm. The materials entrained by the July 16 storm contained a higher proportion of silt- plus clay-sized materials in the <2 mm fraction than the materials collected from storms that produced comparable sediment-runoff concentrations but not debris flows. The difference in materials demonstrates the critical role of the availability of fine-grained wood ash mantling the hillslopes in the runoff-dominated generation of post-wildfire debris flows. The highest sediment-runoff concentrations, again not unique to debris-flow production, were produced from maximum 30 min rainfall intensities greater than 20 mm h-1. Copyright ?? 2001 John Wiley and Sons, Ltd.

Hydrological Processes↗
USGS2013

Infiltration and runoff generation processes in fire-affected soils

Post-wildfire runoff was investigated by combining field measurements and modelling of infiltration into fire-affected soils to predict time-to-start of runoff and peak runoff rate at the plot scale (1 m 2 ). Time series of soil-water content, rainfall and runoff were measured on a hillslope burned by the 2010 Fourmile Canyon Fire west of Boulder, Colorado during cyclonic and convective rainstorms in the spring and summer of 2011. Some of the field measurements and measured soil physical properties were used to calibrate a one-dimensional post-wildfire numerical model, which was then used as a ‘virtual instrument’ to provide estimates of the saturated hydraulic conductivity and high-resolution (1 mm) estimates of the soil-water profile and water fluxes within the unsaturated zone. Field and model estimates of the wetting-front depth indicated that post-wildfire infiltration was on average confined to shallow depths less than 30 mm. Model estimates of the effective saturated hydraulic conductivity, K s , near the soil surface ranged from 0.1 to 5.2 mm h −1 . Because of the relatively small values of K s , the time-to-start of runoff (measured from the start of rainfall), t p , was found to depend only on the initial soil-water saturation deficit (predicted by the model) and a measured characteristic of the rainfall profile (referred to as the average rainfall acceleration, equal to the initial rate of change in rainfall intensity). An analytical model was developed from the combined results and explained 92–97% of the variance of t p , and the numerical infiltration model explained 74–91% of the variance of the peak runoff rates. These results are from one burned site, but they strongly suggest that t p in fire-affected soils (which often have low values of K s ) is probably controlled more by the storm profile and the initial soil-water saturation deficit than by soil hydraulic properties.

Hydrological Processes↗
USGS2014

Floodplain biogeochemical processing of floodwaters in the Atchafalaya River Basin during the Mississippi River flood of 2011

The 2011 flood in the Lower Mississippi resulted in the second highest recorded river flow diverted into the Atchafalaya River Basin (ARB). The higher water levels during the flood peak resulted in high hydrologic connectivity between the Atchafalaya River and floodplain, with up to 50% of the Atchafalaya River water moving off channel. Water quality samples were collected throughout the ARB over the course of the flood event. Significant nitrate (NO 3 -) reduction (75%) occurred within the floodplain, resulting in a total NO 3 - reduction of 16.6% over the flood. The floodplain was a small but measurable source of dissolved reactive phosphorus (SRP) and ammonium (NH 4 +). Collectively, these results from this large flood event suggest that enhancing river-floodplain connectivity through freshwater diversions will reduce NO 3 - loads to the Gulf of Mexico during large annual floods.

Louisiana↗
USGS2014

Coupled hydrological and biogeochemical processes controlling variability of nitrogen species in streamflow during autumn in an upland forest

Autumn is a season of dynamic change in forest streams of the northeastern United States due to effects of leaf fall on both hydrology and biogeochemistry. Few studies have explored how interactions of biogeochemical transformations, various nitrogen sources, and catchment flow paths affect stream nitrogen variation during autumn. To provide more information on this critical period, we studied (1) the timing, duration, and magnitude of changes to stream nitrate, dissolved organic nitrogen (DON), and ammonium concentrations; (2) changes in nitrate sources and cycling; and (3) source areas of the landscape that most influence stream nitrogen. We collected samples at higher temporal resolution for a longer duration than typical studies of stream nitrogen during autumn. This sampling scheme encompassed the patterns and extremes that occurred during base flow and stormflow events of autumn. Base flow nitrate concentrations decreased by an order of magnitude from 5.4 to 0.7 &micro;mol L &minus;1 during the week when most leaves fell from deciduous trees. Changes to rates of biogeochemical transformations during autumn base flow explained the low nitrate concentrations; in-stream transformations retained up to 72% of the nitrate that entered a stream reach. A decrease of in-stream nitrification coupled with heterotrophic nitrate cycling were primary factors in the seasonal nitrate decline. The period of low nitrate concentrations ended with a storm event in which stream nitrate concentrations increased by 25-fold. In the ensuing weeks, peak stormflow nitrate concentrations progressively decreased over closely spaced, yet similarly sized events. Most stormflow nitrate originated from nitrification in near-stream areas with occasional, large inputs of unprocessed atmospheric nitrate, which has rarely been reported for nonsnowmelt events. A maximum input of 33% unprocessed atmospheric nitrate to the stream occurred during one event. Large inputs of unprocessed atmospheric nitrate show direct and rapid effects on forest streams that may be widespread, although undocumented, throughout nitrogen-polluted temperate forests. In contrast to a week-long nitrate decline during peak autumn litterfall, base flow DON concentrations increased after leaf fall and remained high for 2 months. Dissolved organic nitrogen was hydrologically flushed to the stream from riparian soils during stormflow. In contrast to distinct seasonal changes in base flow nitrate and DON concentrations, ammonium concentrations were typically at or below the detection limit, similar to the rest of the year. Our findings reveal couplings among catchment flow paths, nutrient sources, and transformations that control seasonal extremes of stream nitrogen in forested landscapes.

Vermont↗
USGS2015

Crustal deformation in the New Madrid seismic zone and the role of postseismic processes

Global Navigation Satellite System data across the New Madrid seismic zone (NMSZ) in the central United States over the period from 2000 through 2014 are analyzed and modeled with several deformation mechanisms including the following: (1) creep on subsurface dislocations, (2) postseismic frictional afterslip and viscoelastic relaxation from the 1811&ndash;1812 and 1450 earthquakes in the NMSZ, and (3) regional strain. In agreement with previous studies, a dislocation creeping at about 4 mm/yr between 12 and 20 km depth along the downdip extension of the Reelfoot fault reproduces the observations well. We find that a dynamic model of postseismic frictional afterslip from the 1450 and February 1812 Reelfoot fault events can explain this creep. Kinematic and dynamic models involving the Cottonwood Grove fault provide minimal predictive power. This is likely due to the smaller size of the December 1811 event on the Cottonwood Grove fault and a distribution of stations better suited to constrain localized strain across the Reelfoot fault. Regional compressive strain across the NMSZ is found to be less than 3&thinsp;&times;&thinsp;10 &minus;9 /yr. If much of the present-day surface deformation results from afterslip, it is likely that many of the earthquakes we see today in the NMSZ are aftershocks from the 1811&ndash;1812 New Madrid earthquakes. Despite this conclusion, our results are consistent with observations and models of intraplate earthquake clustering. Given this and the recent paleoseismic history of the region, we suggest that seismic hazard is likely to remain significant.

Journal of Geophysical Research B: Solid Earth↗
USGS2015

Earthquake rupture process recreated from a natural fault surface

What exactly happens on the rupture surface as an earthquake nucleates, spreads, and stops? We cannot observe this directly, and models depend on assumptions about physical conditions and geometry at depth. We thus measure a natural fault surface and use its 3D coordinates to construct a replica at 0.1 m resolution to obviate geometry uncertainty. We can recreate stick-slip behavior on the resulting finite element model that depends solely on observed fault geometry. We clamp the fault together and apply steady state tectonic stress until seismic slip initiates and terminates. Our recreated M~1 earthquake initiates at contact points where there are steep surface gradients because infinitesimal lateral displacements reduce clamping stress most efficiently there. Unclamping enables accelerating slip to spread across the surface, but the fault soon jams up because its uneven, anisotropic shape begins to juxtapose new high-relief sticking points. These contacts would ultimately need to be sheared off or strongly deformed before another similar earthquake could occur. Our model shows that an important role is played by fault-wall geometry, though we do not include effects of varying fluid pressure or exotic rheologies on the fault surfaces. We extrapolate our results to large fault systems using observed self-similarity properties, and suggest that larger ruptures might begin and end in a similar way, though the scale of geometrical variation in fault shape that can arrest a rupture necessarily scales with magnitude. In other words, fault segmentation may be a magnitude dependent phenomenon and could vary with each subsequent rupture.

Journal of Geophysical Research B: Solid Earth↗
USGS2018

Response to comment by Walker et al. on “From data to decisions: Processing information, biases, and beliefs for improved management of natural resources and environments”

Our different kinds of minds and types of thinking affect the ways we decide, take action, and cooperate (or not). The comment by Walker et al. (2018, https://doi.org/10.1002/2017EF000750 ) illustrates several points made by Glynn et al. (2017, https://doi.org/10.1002/2016EF000487 ) and many other articles. Namely, biases and beliefs often drive scientific reasoning, and scientists, just like other humans, are intimately attached to their values and heuristics. Scientists, just like many other people, also tend to read and interpret text in ways that best match their individual perceptions of a problem or issue: in many cases paraphrasing and changing the meaning of what they read to better match their initial ideas. Walker et al. are doing interesting and important research on uncertainty. Nonetheless, they misinterpret the work, assumptions, and conclusions brought forth by Glynn et al. (2017, https://doi.org/10.1002/2016EF000487 ).

Earth's Future↗
USGS2018

Processes and facies relationships in a Lower(?) Devonian rocky shoreline depositional environment, East Lime Creek Conglomerate, south‐western Colorado, USA

Rocky shorelines are relatively common features along modern coastlines, but few have been recognized in the geological record. The hard substrates of rocky shorelines telescope the width of offshore marine environments, thus the diagnostic deposits observed in such settings today have a low preservation potential due to small accommodation space and high‐energy conditions. This study recognized previously overlooked, laterally extensive Lower(?) Devonian rocky shoreline deposits in the San Juan Mountains of south‐western Colorado. The newly defined lithostratigraphic unit, the East Lime Creek Conglomerate (ELCC), is 0–23 m thick, unconformably overlying Proterozoic crystalline rocks and unconformably overlain by the Upper Devonian Ignacio Formation and/or Elbert Formation. The unit mostly consists of clast‐supported cobble‐boulder conglomerate with rounded quartzite clasts up to 1.4 m in length interbedded with thin sandstone layers and lenses. Sandstones in the ELCC are distinguished from unconformably overlying Upper Devonian sedimentary rocks because they have sericite cements. Most importantly, there are buttressing relationships between the ELCC and underlying Proterozoic crystalline rocks interpreted as palaeo‐sea cliffs, palaeo‐wave‐cut platforms and palaeo‐tombolos. A proposed rocky shoreline facies model includes headlands with upper shoreface‐beachface tabular cobble‐boulder gravels sourced from rock fall talus, nearshore subaqueous debris‐flow deposits and intervening pocket beaches with imbricated, stratified pebble‐cobble gravel sheets. Palaeocurrent data ( n = 338) from clast long‐axis orientations, imbrication and cross‐bedding indicate south‐to‐north transport roughly onshore‐offshore to a coastline consisting of alternating rocky headlands and pocket beaches. This Lower(?) Devonian unit documents a previously unrecognized episode in the geological history of south‐western Colorado.

Colorado↗
USGS2016

Small-scale barriers mitigate desertification processes and enhance plant recruitment in a degraded semiarid grassland

Anthropogenic desertification is a problem that plagues drylands globally; however, the factors which maintain degraded states are often unclear. In Canyonlands National Park on the Colorado Plateau of southeastern Utah, many degraded grasslands have not recovered structure and function >40 yr after release from livestock grazing pressure, necessitating active restoration. We hypothesized that multiple factors contribute to the persistent degraded state, including lack of seed availability, surficial soil-hydrological properties, and high levels of spatial connectivity (lack of perennial vegetation and other surface structure to retain water, litter, seed, and sediment). In combination with seeding and surface raking treatments, we tested the effect of small barrier structures (&ldquo;ConMods&rdquo;) designed to disrupt the loss of litter, seed and sediment in degraded soil patches within the park. Grass establishment was highest when all treatments (structures, seed addition, and soil disturbance) were combined, but only in the second year after installation, following favorable climatic conditions. We suggest that multiple limiting factors were ameliorated by treatments, including seed limitation and microsite availability, seed removal by harvester ants, and stressful abiotic conditions. Higher densities of grass seedlings on the north and east sides of barrier structures following the summer months suggest that structures may have functioned as artificial &ldquo;nurse-plants&rdquo;, sheltering seedlings from wind and radiation as well as accumulating wind-blown resources. Barrier structures increased the establishment of both native perennial grasses and exotic annuals, although there were species-specific differences in mortality related to spatial distribution of seedlings within barrier structures. The unique success of all treatments combined, and even then only under favorable climatic conditions and in certain soil patches, highlights that restoration success (and potentially, natural regeneration) often is contingent on many interacting factors.

Utah↗
USGS2016

Forecasting tidal marsh elevation and habitat change through fusion of Earth observations and a process model

Reducing uncertainty in data inputs at relevant spatial scales can improve tidal marsh forecasting models, and their usefulness in coastal climate change adaptation decisions. The Marsh Equilibrium Model (MEM), a one-dimensional mechanistic elevation model, incorporates feedbacks of organic and inorganic inputs to project elevations under sea-level rise scenarios. We tested the feasibility of deriving two key MEM inputs—average annual suspended sediment concentration (SSC) and aboveground peak biomass—from remote sensing data in order to apply MEM across a broader geographic region. We analyzed the precision and representativeness (spatial distribution) of these remote sensing inputs to improve understanding of our study region, a brackish tidal marsh in San Francisco Bay, and to test the applicable spatial extent for coastal modeling. We compared biomass and SSC models derived from Landsat 8, DigitalGlobe WorldView-2, and hyperspectral airborne imagery. Landsat 8-derived inputs were evaluated in a MEM sensitivity analysis. Biomass models were comparable although peak biomass from Landsat 8 best matched field-measured values. The Portable Remote Imaging Spectrometer SSC model was most accurate, although a Landsat 8 time series provided annual average SSC estimates. Landsat 8-measured peak biomass values were randomly distributed, and annual average SSC (30 mg/L) was well represented in the main channels (IQR: 29–32 mg/L), illustrating the suitability of these inputs across the model domain. Trend response surface analysis identified significant diversion between field and remote sensing-based model runs at 60 yr due to model sensitivity at the marsh edge (80–140 cm NAVD88), although at 100 yr, elevation forecasts differed less than 10 cm across 97% of the marsh surface (150–200 cm NAVD88). Results demonstrate the utility of Landsat 8 for landscape-scale tidal marsh elevation projections due to its comparable performance with the other sensors, temporal frequency, and cost. Integration of remote sensing data with MEM should advance regional projections of marsh vegetation change by better parameterizing MEM inputs spatially. Improving information for coastal modeling will support planning for ecosystem services, including habitat, carbon storage, and flood protection.

California↗
USGS2019

Simulating demography, genetics, and spatially explicit processes to inform reintroduction of a threatened char

The success of species reintroductions can depend on a combination of environmental, demographic, and genetic factors. Although the importance of these factors in the success of reintroductions is well‐accepted, they are typically evaluated independently, which can miss important interactions. For species that persist in metapopulations, movement through and interaction with the landscape is predicted to be a vital component of persistence. Simulation‐based approaches are a promising technique for evaluating the independent and combined effects of these factors on the outcome of various reintroduction and associated management actions. We report results from a simulation study of bull trout ( Salvelinus confluentus ) reintroduction to three watersheds of the Pend Oreille River system in northeastern Washington State, USA. We used an individual‐based, spatially explicit simulation model to evaluate how reintroduction strategies, life history variation, and riverscape structure (e.g., network topology) interact to influence the demographic and genetic characteristics of reintroduced bull trout populations in three watersheds. Simulation scenarios included a range of initial genetic stocks (informed by empirical bull trout genetic data), variation in migratory tendency and life history, and two landscape connectivity alternatives representing a connected network (isolation‐by‐distance) and a fragmented network (isolation‐by‐barrier, using the known existing barriers). A novel feature of these simulations was the ability to consider the interaction of both demographic and genetic (i.e., demogenetic) factors in riverscapes with implicit asymmetric movement probabilities across the barriers. We found that connectivity (presence or absence of barriers) had the largest effect on demographic and genetic outcomes over 200 yr, with a greater effect than both initial genetic diversity and life history variation. We also identified regions of the study system in which bull trout populations persisted across a wide range of demographic, life history, and environmental connectivity parameters. Finally, we found no evidence that initial neutral genetic diversity influenced genetic diversity and structure after 200 yr; instead, genetic drift due to stray rate and population isolation dominated and erased any initial differences in genetic diversity. Our results highlight the utility of spatially explicit demogenetic approaches in exploring and understanding population dynamics—and their implications for management strategies—in fresh waters.

Washington↗
USGS2023

Ecotoxicological studies indicate that sublethal and lethal processes limit insect-mediated contaminant flux

Merolimnic insects can accumulate and transport considerable amounts of aquatic contaminants to terrestrial systems. The rate of contaminant biotransport, termed insect-mediated contaminant flux (IMCF), depends on emergent insect biomass and contaminant accumulation, both functions of environmental concentration. Here we develop a mathematical model of IMCF and apply it to three ecotoxicological studies obtained through the U.S. Environmental Protection Agency's ECOTOX database to determine at which concentration maximum IMCF occurs. Model results demonstrate that the maximum IMCF depends on competing rates of biomass loss and contaminant accumulation and does not necessarily occur at the highest insect or environmental contaminant concentration. Additionally, modeling results suggest that sublethal contaminant effects (e.g., decreased growth) on insect biomass can be an important and potentially underappreciated control on IMCF.

Environmental Toxicology and Chemistry↗
USGS2018

Patterns, pace and processes of water-quality variability: Examples from a long-studied estuary

Environmental time series have rich information content that is invaluable for measuring and understanding changes over time and guiding policies to manage change. I extracted information from measurements of 10 water‐quality constituents in upper San Francisco Bay from 1975 to 2016, one of the longest observational records in a U.S. estuary. Changes were detected at every time scale captured by monthly sampling. Long‐term trends included increased ammonium (+53%), nitrate + nitrate (+50%), silicate (+14%), Secchi depth (+42%), and decreased chlorophyll a (Chl a ) (−74%) and suspended particulate matter (−45%). Changes at the decadal scale included abrupt shifts (Chl a , nitrate + nitrite) and oscillations between shorter trends of increase and decrease (Secchi depth, phosphate). Long‐term trends were not expressed equally across all seasons, and seasonal patterns of change varied across constituents. These examples illustrate key features of environmental variability at the land–sea interface: (1) water‐quality components change continually at time scales from months to decades; (2) patterns of seasonal, multiyear, and multidecadal change are complex and vary across constituents; (3) primary drivers of change are freshwater inflow, the master regulator of estuarine dynamics, and human activities such as river damming, water diversions, wastewater discharge, environmental policies, and species introductions; (4) extracting the full information content of time series requires multiple analyses, each revealing a different layer of insight into how changes develop over time; (5) water‐quality variability is nonstationary, so future changes cannot be forecast reliably; (6) repeated observation is an essential method of Earth system science with applications in the design and performance measures of environmental policies.

California↗
USGS2017

Comparison and contrast of processes of soil formation in the San Timoteo Badlands with chronosequences in California

The degree of soil development associated with geomorphic surfaces in the San Timoteo Badlands area (STB), California allows correlation of the surfaces. Soil development indices, based on field descriptions and laboratory analysis, provide a basis for comparison of these soils to each other and to dated soil chronosequences at Cajon Pass, Merced, and Anza, California. The soils in this study record a complex tectonic and climatic history, include preserved surfaces intermediate in soil development to those at Cajon Pass, and do not preserve a record of the major late Pleistocene to early Holocene aggradational event observed elsewhere in the region. Pedogenesis is similar to that of regional soils in southern California. With time, the profiles develop progressively thicker argillic horizons, and they increase in redness, clay, and secondary iron oxide content, and amount and thickness of clay films. With increasing soil age the ratio of dithionite-extractable Fe to total Fe increases, while the ratios of Fe 2+ to Fe 3+ and of oxalate- to dithionite-extractable Fe decrease. These trends are more subdued in the STB than elsewhere in the region, perhaps because the slightly higher temperature and lower precipitation of the STB may slow the rate of transformation of ferrihydrite to hematite. We use weighted means of soil development indices and rubification index values to compare the soils of this study to those at Cajon Pass, Merced, and Anza. Wide-range age estimates using the maximum limits of these values result in 27,500–305,000 yr for surface Q1, 43,000–570,000 yr for surface Q2, and 300,000–700,000 yr for surfaces Q3 and Q4. Weighing most heavily the nearest dated chronosequences, Anza and Cajon Pass, yields best estimates of age. This method constrains the time of formation of the Q1 surface to between 27,500 and 67,000 yr and the Q2 surface to 43,000–67,000 yr, and does not change the age estimates of the older surfaces.

California↗
USGS2003

Geochemical processes controlling transport of arsenic in groundwater: A review of adsorption

Adsorption is the predominate mechanism controlling transport of arsenic in many ground water systems. Hydrous oxides of iron, aluminum, and manganese, and clay minerals are commonly associated with aquifer solids and have been shown to be significant adsorbents of arsenic. The extent of arsenic adsorption is influenced by the chemistry of the aqueous phase including pH, arsenic speciation, and the presence and concentration of competing ions. Under moderately reducing conditions, trivalent arsenite is stable and adsorption increases with increasing pH. In an oxidizing environment, arsenate is stable and adsorption decreases with increasing pH. The presence of phosphate, sulfate, carbonate, silica, and other anions have been shown to decrease adsorption of arsenic to varying degrees. The effects of complex aqueous and solid phase chemistry on arsenic adsorption are best simulated using surface complexation models. Coupling of such models with hydrologic solute transport codes provide a powerful method for predicting the spatial and temporal distribution of arsenic in ground water.

Book chapter↗