SEARCH · Geology Reports

Results for “Contribution”

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

SourceResearch connectionFrom yearThrough year

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

At least 1,693 records · Page 94Linked to original sources

USGS2025

Light absorbing particles deposited to snow cover across the Upper Colorado River Basin, Colorado Rocky Mountains, 2013-16: Interannual variations from multiple natural and anthropogenic sources

Atmospheric particulate matter (PM) as light-absorbing particles (LAPs) deposited to snow cover can result in early onset and rapid snow melting, challenging management of downstream water resources. We identified LAPs in 38 snow samples (water years 2013–2016) from the mountainous Upper Colorado River basin by comparing among laboratory-measured spectral reflectance, chemical, physical, and magnetic properties. Dust sample reflectance, averaged over the wavelength range of 0.35–2.50 μm, varied by a factor of 1.9 (range, 0.2300–0.4444) and was suppressed mainly by three components: (a) carbonaceous matter measured as total organic carbon (1.6–22.5 wt. %) including inferred black carbon, natural organic matter, and carbon-based synthetic, black road-tire-wear particles, (b) dark rock and mineral particles, indicated by amounts of magnetite (0.11–0.37 wt. %) as their proxy, and (c) ferric oxide minerals identified by reflectance spectroscopy and magnetic properties. Fundamental compositional differences were associated with different iron oxide groups defined by dominant hematite, goethite, or magnetite. These differences in iron oxide mineralogy are attributed to temporally varying source-area contributions implying strong interannual changes in regional source behavior, dust-storm frequency, and (or) transport tracks. Observations of dust-storm activity in the western U.S. and particle-size averages for all samples (median, 25 μm) indicated that regional dust from deserts dominated mineral-dust masses. Fugitive contaminants, nevertheless, contributed important amounts of LAPs from many types of anthropogenic sources.

Colorado↗

USGS2025

Organic matter availability and the production of methane-derived dissolved organic carbon at methane seeps

Methane seeps located along continental margins and slopes export sedimentary methane and dissolved organic carbon (DOC) into the ocean. The flux of these reduced carbon molecules from the seafloor into the ocean impacts ocean chemistry and supports deep-sea life. While significant effort has been made to understand how the anaerobic oxidation of methane (AOM) regulates the release of methane from the seafloor, little is known about the production of DOC in association with AOM or its flux and fate in the ocean. We hypothesize a mechanism for methane incorporation into DOC at seeps and investigate the relationship between sediment total organic carbon (TOC) availability and the incorporation of methane-derived carbon into DOC at four methane seep regions along the Cascadia margin, with a range of microbial and thermogenic methane sources. At sites with <2.0 wt.% TOC (Hydrate Ridge and Bullseye Vent), up to 60%–80% of carbon in DOC is methane-carbon, much more than sites with >2.0 wt.% TOC (Astoria Canyon and Barkley Canyon). We attribute the greater methane contribution at the more TOC-limited sites to a greater role of AOM in the carbon cycle, whereas at the organic matter-rich sites, microbial competition for sulfate as an electron acceptor for organic matter decomposition limits AOM and hence the transfer of carbon from methane to DOC. We estimate that the global diffusive flux of methane-derived DOC from the seafloor is 0.07–10.1 Tg C/yr, contributing to the stock of DOC present in the deep ocean and/or fueling the deep-sea microbial loop.

California, Oregon, Washington↗

USGS2025

Wetlands, groundwater and seasonality influence the spatial distribution of stream chemistry in a low-relief catchment

Evaluating stream water chemistry patterns provides insight into catchment ecosystem and hydrologic processes. Spatially distributed patterns and controls of stream solutes are well-established for high-relief catchments where solute flow paths align with surface topography. However, the controls on solute patterns are poorly constrained for low-relief catchments where hydrogeologic heterogeneities and river corridor features, like wetlands, may influence water and solute transport. Here, we provide a data set of solute patterns from 58 synoptic surveys across 28 sites and over 32 months in a low-relief wetland-rich catchment to determine the major surface and subsurface controls along with wetland influence across the catchment. In this low-relief catchment, the expected wetland storage, processing, and transport of solutes is only apparent in solute patterns of the smallest subcatchments. Meanwhile, downstream seasonal and wetland influence on observed chemistry can be masked by large groundwater contributions to the main stream channel. These findings highlight the importance of incorporating variable groundwater contributions into catchment-scale studies for low-relief catchments, and that understanding the overall influence of wetlands on stream chemistry requires sampling across various spatial and temporal scales. Therefore, in low-relief wetland-rich catchments, given the mosaic of above and below ground controls on stream solutes, modeling efforts may need to include both surface and subsurface hydrological data and processes.

Michigan↗

USGS2012

A model for the generation of two-dimensional surf beat

A finite difference model predicting group-forced long waves in the nearshore is constructed with two interacting parts: an incident wave model providing time-varying radiation stress gradients across the nearshore, and a long-wave model which solves the equations of motion for the forcing imposed by the incident waves. Both shallow water group-bound long waves and long waves generated by a time-varying breakpoint are simulated. Model-generated time series are used to calculate the cross correlation between wave groups and long waves through the surf zone. The cross-correlation signal first observed by Tucker (1950) is well predicted. For the first time, this signal is decomposed into the contributions from the two mechanisms of leaky mode forcing. Results show that the cross-correlation signal can be explained by bound long waves which are amplified, though strongly modified, through the surf zone before reflection from the shoreline. The breakpoint-forced long waves are added to the bound long waves at a phase of pi/2 and are a secondary contribution owing to their relatively small size.

Journal of Geophysical Research C: Oceans↗

USGS2008

Preliminary analysis of the hydrologic and geochemical controls on acid‐neutralizing capacity in two acidic seepage lakes in Florida

In late 1988, parallel studies of Lake Five‐O ( p H 5.14) in the Florida panhandle and Lake Barco ( p H 4.50) in north central Florida were initiated to develop hydrologic and major ion budgets of these lakes as part of an overall effort to improve understanding of the hydrologic, depositional, and biogeochemical factors that control acid‐neutralizing capacity (ANC) in seepage lakes. Preliminary findings from these studies indicate that earlier perceptions of lake hydrology and mechanisms of ANC regulation in Florida seepage lakes may have to be revised. The traditional perspective of seepage lakes in the Florida panhandle views these systems as dominated by precipitation inputs and that ANC regulation is due largely to in‐lake processes. Our results for Lake Five‐O show modest to steep hydraulic gradients almost entirely around the lake. In addition, the horizontal hydraulic conductivity of the surficial aquifer is high (8–74 m day −1 ), indicating that large quantities of groundwater flow into Lake Five‐O. Calculations of net groundwater flow from hydrologic budgets also indicate that groundwater may comprise at least 38 to 46% of the total inflow. For Lake Barco, net flow estimates of the minimum groundwater inflow range from 5 to 14% of total inflow. Enrichment factor and ion flux calculations for Lake Five‐O and Lake Barco indicate that terrestrial as well as in‐lake processes contribute significantly to ANC regulation. The extent that terrestrial processes contribute to ANC generation is directly related to the magnitude of groundwater inflow as well as the degree of ion enrichment or depletion that occurs in the surficial aquifer. Net ANC generation in both study lakes was dominated by anion retention (NO 3 − and SO 4 2− ). Where previous studies concluded that in‐lake reduction was the primary sink for SO 4 2− , our preliminary calculations show that adsorption of SO 4 2− within the watershed is perhaps twice as important as in‐lake reduction as a source of ANC. Net base cation enrichment in both lakes was negligible.

Florida↗

USGS1997

Variation of the 3-μm absorption feature on Mars: observations over eastern Valles Marineris by the mariner 6 infrared spectrometer

A new approach for calibration of the shortest wavelength channel (1.8 to 6.0 μm) of the Mariner 6 infrared spectrometer was derived. This calibration provides a new description of the instrument response function from 1.8 to 3.7 μm and accounts for the thermal contribution to the signal at longer wavelengths. This allows the two segments from 1.8 to 6 μm to be merged into a single spectrum. The broad water of hydration absorption spans these two segments and is examined in these merged spectra using a method of band integration. Unlike previous analyses which rely on ratios at two wavelengths, the integration method can assess the band strength independently from the albedo in the near infrared. Spectra taken over the eastern end of the Valles Marineris are examined for variations of the band-integrated value, and three distinct clusters are found. Within the estimated uncertainty, two clusters (both low and high albedo) have approximately the same integrated band depth. The third cluster (medium albedo) has an integrated band depth about 10% higher. This difference cannot be systematically attributed to either surface or atmospheric parameters and suggests variation in the amount of water either chemically or physically bound in surface materials. Approximately one-half of the high integrated band depth cluster is associated with chaotic terrain at the source of outflow channels, the other half occurs over lower inertia plains adjacent to chasmata. This suggests both surface physical properties and mineralogy as well as water in exchange with the atmosphere contribute to the 3-μm bound water absorption.

Journal of Geophysical Research E: Planets↗

USGS1996

Metastable mantle phase transformations and deep earthquakes in subducting oceanic lithosphere

Earth's deepest earthquakes occur as a population in subducting or previously subducted lithosphere at depths ranging from about 325 to 690 km. This depth interval closely brackets the mantle transition zone, characterized by rapid seismic velocity increases resulting from the transformation of upper mantle minerals to higher-pressure phases. Deep earthquakes thus provide the primary direct evidence for subduction of the lithosphere to these depths and allow us to investigate the deep thermal, thermodynamic, and mechanical ferment inside slabs. Numerical simulations of reaction rates show that the olivine → spinel transformation should be kinetically hindered in old, cold slabs descending into the transition zone. Thus wedge-shaped zones of metastable peridotite probably persist to depths of more than 600 km. Laboratory deformation experiments on some metastable minerals display a shear instability called transformational faulting. This instability involves sudden failure by localized superplasticity in thin shear zones where the metastable host mineral transforms to a denser, finer-grained phase. Hence in cold slabs, such faulting is expected for the polymorphic reactions in which olivine transforms to the spinel structure and clinoenstatite transforms to ilmenite. It is thus natural to hypothesize that deep earthquakes result from transformational faulting in metastable peridotite wedges within cold slabs. This consideration of the mineralogical states of slabs augments the traditional largely thermal view of slab processes and explains some previously enigmatic slab features. It explains why deep seismicity occurs only in the approximate depth range of the mantle transition zone, where minerals in downgoing slabs should transform to spinel and ilmenite structures. The onset of deep shocks at about 325 km is consistent with the onset of metastability near the equilibrium phase boundary in the slab. Even if a slab penetrates into the lower mantle, earthquakes should cease at depths near 700 km, because the seismogenic phase transformations in the slab are completed or can no longer occur. Substantial metastability is expected only in old, cold slabs, consistent with the observed restriction of deep earthquakes to those settings. Earthquakes should be restricted to the cold cores of slabs, as in any model in which the seismicity is temperature controlled, via the distribution of metastability. However, the geometries of recent large deep earthquakes pose a challenge for any such models. Transformational faulting may give insight into why deep shocks lack appreciable aftershocks and why their source characteristics, including focal mechanisms indicating localized shear failure rather than implosive deformation, are so similar to those of shallow earthquakes. Finally, metastable phase changes in slabs would produce an internal source of stress in addition to those due to the weight of the sinking slab. Such internal stresses may explain the occurrence of earthquakes in portions of lithosphere which have foundered to the bottom of the transition zone and/or are detached from subducting slabs. Metastability in downgoing slabs could have considerable geodynamic significance. Metastable wedges would reduce the negative buoyancy of slabs, decrease the driving force for subduction, and influence the state of stress in slabs. Heat released by metastable phase changes would raise temperatures within slabs and facilitate the transformation of spinel to the lower mantle mineral assemblage, causing slabs to equilibrate more rapidly with the ambient mantle and thus contribute to the cessation of deep seismicity. Because wedge formation should occur only for fast subducting slabs, it may act as a “parachute” and contribute to regulating plate speeds. Wedge formation would also have consequences for mantle evolution because the density of a slab stagnated near the bottom of the transition zone would increase as it heats up and the wedge transforms to denser spinel, favoring the subsequent sinking of the slab into the lower mantle.

Reviews of Geophysics↗

USGS1997

The roles of time and displacement in velocity-dependent volumetric strain of fault zones

The relationship between measured friction μ A and volumetric strain during frictional sliding was determined using a rate and state variable dependent friction constitutive equation, a common work balance relating friction and volume change, and two types of experimental faults: initially bare surfaces of Westerly granite and rock surfaces separated by a 1 mm layer of less than 90 μm Westerly granite gouge. The constitutive equation is the sum of a constant term representing the nominal resistance to sliding and two smaller terms: a rate dependent term representing the shear viscosity of the fault surface (direct effect), and a term which represents variations in the area of contact (evolution effect). The work balance relationship requires that μ A differs from the frictional resistance that leads to shear heating by the derivative of fault normal displacement with respect shear displacement, d δ n / d δ s . An implication of this relationship is that the rate dependence of d δ n / d δ s contributes to the rate dependence of μ A . Experiments show changes in sliding velocity lead to changes in both fault strength and volume. Analysis of data with the rate and state equations combined with the work balance relationship preclude the conventional interpretation of the direct effect in the rate and state variable constitutive equations. Consideration of a model bare surface fault consisting of an undeformable indentor sliding on a deformable surface reveals a serious flaw in the work balance relationship if volume change is time-dependent. For the model, at zero slip rate indentation creep under the normal load leads to time-dependent strengthening of the fault surface but, according to the work balance relationship, no work is done because compaction or dilatancy can only be induced by shearing. Additional tests on initially bare surfaces and gouges show that fault normal strain in experiments is time-dependent, consistent with the model. This time-dependent fault normal strain, which is not accounted for in the work balance relationship, explains the inconsistency between the constitutive equations and the work balance. For initially bare surface faults, all rate dependence of volume change is due to time dependence. Similar results are found for gouge. We conclude that μ A reflects the frictional resistance that results in shear heating, and no correction needs to be made for the volume changes. The result that time-dependent volume changes do not contribute to μ A is a general result and extends beyond these experiments, the simple indentor model and particular constitutive equations used to illustrate the principle.

Journal of Geophysical Research B: Solid Earth↗

USGS1997

Comparing nocturnal eddy covariance measurements to estimates of ecosystem respiration made by scaling chamber measurements at six coniferous boreal sites

During the growing season, nighttime ecosystem respiration emits 30–100% of the daytime net photosynthetic uptake of carbon, and therefore measurements of rates and understanding of its control by the environment are important for understanding net ecosystem exchange. Ecosystem respiration can be measured at night by eddy covariance methods, but the data may not be reliable because of low turbulence or other methodological problems. We used relationships between woody tissue, foliage, and soil respiration rates and temperature, with temperature records collected on site to estimate ecosystem respiration rates at six coniferous BOREAS sites at half-hour or 1-hour intervals, and then compared these estimates to nocturnal measurements of CO 2 exchange by eddy covariance. Soil surface respiration was the largest source of CO 2 at all sites (48–71%), and foliar respiration made a large contribution to ecosystem respiration at all sites (25–43%). Woody tissue respiration contributed only 5–15% to ecosystem respiration. We estimated error for the scaled chamber predictions of ecosystem respiration by using the uncertainty associated with each respiration parameter and respiring biomass value. There was substantial uncertainty in estimates of foliar and soil respiration because of the spatial variability of specific respiration rates. In addition, more attention needs to be paid to estimating foliar respiration during the early part of the growing season, when new foliage is growing, and to determining seasonal trends of soil surface respiration. Nocturnal eddy covariance measurements were poorly correlated to scaled chamber estimates of ecosystem respiration ( r 2 =0.06–0.27) and were consistently lower than scaled chamber predictions (by 27% on average for the six sites). The bias in eddy covariance estimates of ecosystem respiration will alter estimates of gross assimilation in the light and of net ecosystem exchange rates over extended periods.

Journal of Geophysical Research D: Atmospheres↗

USGS1997

Case study modeling of turbulent and mesoscale fluxes over the BOREAS region

Results from aircraft and surface observations provided evidence for the existence of mesoscale circulations over the Boreal Ecosystem-Atmosphere Study (BOREAS) domain. Using an integrated approach that included the use of analytical modeling, numerical modeling, and data analysis, we have found that there are substantial contributions to the total budgets of heat over the BOREAS domain generated by mesoscale circulations. This effect is largest when the synoptic flow is relatively weak, yet it is present under less favorable conditions, as shown by the case study presented here. While further analysis is warranted to document this effect, the existence of mesoscale flow is not surprising, since it is related to the presence of landscape patches, including lakes, which are of a size on the order of the local Rossby radius and which have spatial differences in maximum sensible heat flux of about 300 W m-2. We have also analyzed the vertical temperature profile simulated in our case study as well as high-resolution soundings and we have found vertical profiles of temperature change above the boundary layer height, which we attribute in part to mesoscale contributions. Our conclusion is that in regions with organized landscapes, such as BOREAS, even with relatively strong synoptic winds, dynamical scaling criteria should be used to assess whether mesoscale effects should be parameterized or explicitly resolved in numerical models of the atmosphere.

Journal of Geophysical Research D: Atmospheres↗

USGS2012

Heat flow and energetics of the San Andreas fault zone

Approximately 100 heat flow measurements in the San Andreas fault zone indicate (1) there is no evidence for local factional heating of the main fault trace at any latitude over a 1000-km length from Cape Mendocino to San Bernardino, (2) average heat flow is high (∼2 HFU, ∼80 mW m −2 ) throughout the 550-km segment of the Coast Ranges that encloses the San Andreas fault zone in central California; this broad anomaly falls off rapidly toward the Great Valley to the east, and over a 200-km distance toward the Mendocino Triple Junction to the northwest. As others have pointed out, a local conductive heat flow anomaly would be detectable unless the frictional resistance allocated to heat production on the main trace were ≲100 bars. Frictional work allocated to surface energy of new fractures is probably unimportant, and hydrologic convection is not likely to invalidate the conduction assumption, since the heat discharge by thermal springs near the fault is negligible. Explanations for the low dynamic friction fall into two intergradational classes: those in which the fault is weak all of the time and those in which it is weak only during earthquakes (possibly just large ones). The first class includes faults containing anomalously weak gouge materials and faults containing materials with normal frictional properties under near-lithostatic steady state fluid pressures. In the second class, weakening is caused by the event (for example, a thermally induced increase in fluid pressure, dehydration of clay minerals, or acoustic fluidization). In this class, unlike the first, the average strength and ambient tectonic shear stress may be large, ∼1 kbar, but the stress allocated to elastic radiation (the apparent stress) must be of similar magnitude, an apparent contradiction with seismic estimates. Unless seismic radiation is underestimated for large earthquakes, it is difficult to justify average tectonic stresses on the main trace of the San Andreas fault in excess of ∼200 bars. The development of the broad Coast Range heat flow anomaly southward from Cape Mendocino suggests that heat flow increases by a factor of 2 within 4 m.y. after the passage of the Mendocino Triple Junction. This passage leaves the San Andreas transform fault zone in its wake; the depth of the anomalous sources cannot be much greater than the depth of the seismogenic layer. Some of the anomalous heat may be supplied by conduction from the warmer mantle that must occur south of the Mendocino transform (where there is no subducting slab), and some might be supplied by shear heating in the fault zone. With no contribution from shear heating, extreme mantle upwelling would be required, and asthenosphere conditions should exist today at depths of only ∼20 km in the northernmost Coast Ranges. If there is an appreciable contribution from shear heating, the heat flow constraint implies that the seismogenic layer is partially decoupled at its base and that the basal traction is in the sense that resists right lateral motion on the fault(s). As a result of these basal tractions, the average shearing stress in the seismogenic layer would increase with distance from the main fault, and the seismogenic layer would offer substantial resistance to plate motion even though resistance on the main fault might be negligible. These speculative models have testable consequences.

Journal of Geophysical Research Solid Earth↗

USGS2012

Magnetic properties of the Bay of Islands ophiolite suite and implications for the magnetization of oceanic crust

Rock magnetic properties, opaque mineralogy, and degree of metamorphism were determined for 101 unoriented samples from the North Arm and Blow-Me-Down massifs of the Bay of Islands ophiolite complex, Newfoundland. The weathered and metamorphosed extrusive basalt samples have a weak, secondary magnetization arising from oxidation and exsolution of ilmenite of unknown origin. The initial magnetization of the underlying sheeted dike complex appears to have been destroyed by hydrothermal alteration soon after formation. The magnetic intensity of the gabbroic samples increases as the degree of alteration increases, with the highly altered upper metagabbros having an average intensity of 3×10 −3 emu/c 3 . Because magnetization of the metagabbro samples is related to nonpervasive, variable alteration, these crustal units are unlikely to make a significant contribution to lineated magnetic anomalies. A compilation of our results and other studies suggests a model in which oceanic crust magnetization results from an upper extrusive basalt source layer, roughly 600 m thick, with no contribution from a deeper source layer recognizable from these Bay of Islands data.

Newfoundland↗

USGS2010

Water quality.

Significant contributions in the broad area of water quality over the quadrennium 1975-78 are highlighted. This summare is concerned primarily with physical and chemical aspects of water quality. The diversity of subject areas within the topic heading and the large volume of published research results necessitated the selection of representative contributions. Over 400 references are cited which are believed to be indicative of general trends in research and of the more important developments during this period.- from Authors

Reviews of Geophysics and Space Physics↗

USGS2010

Earthquakes of the Holocene

Geologic studies of earthquakes involve mapping of coseismic surface faulting and secondary deformation from historical events, trenching and geomorphic analysis to define the timing and size of past earthquakes, and investigations of fault zone structure and geometry in both unconsolidated sediments and bedrock. This research is now being referred to as paleoseismology, seismic geology, and earthquake geology. Since the mid-1970s, it has led to some of the most exciting and important contributions to the understanding of earthquake behavior in space and time [Hanks, 1985; Allen, 1986]. The present report is the first of what is hoped will become a regularly contributed summary of the geologic aspects of the study of earthquakes.

Reviews of Geophysics↗

USGS2014

Glacier‐measurements in the United States

It is doubtless known to many members of the Section of Hydrology that prior to the World War there existed an International Glacier Commission (la Commission Internationale des Glaciers) whose self‐appointed task it was to gather and publish the results of all available measurements showing the annual variations in length—advance or recession—of glaciers in different parts of the world. On this Commission the United States was represented by two men who made frequent contributions to the International record.

Eos, Transactions, American Geophysical Union↗

USGS2014

Recovery of ground‐water supplies by pumping from watertable ponds

This paper summarizes a study made to determine whether ground‐water storage in the glacial deposits in the vicinity of Fresh Pond, a water‐table pond used as part of the public supply of the City of Cambridge, Massachusetts, is available to augment the supply from surface‐water reservoirs. Test wells were drilled; water level, well and test boring, and quality of water data, and data on Inflow and outflow from Fresh Pond were collected, and seismic profiles were run, to determine whether the glacial deposits in the vicinity of the pond contribute appreciable quantities of water to the pond. The data show that water levels in certain wells fluctuate with the level of Fresh Pond and that certain other water levels do not, indicating ground‐water contributions to the pond when it is at low stages and increases in ground‐water storage when pond stages are high. The largest changes in ground‐water storage occur in the glacial deposits underlying the banks to the west and northwest of the pond. Mineral characteristics of the well and pond waters offer no evidence as to the direction of ground‐water flow.

Massachusetts↗

USGS2014

The international hydrological decade

Work toward establishing a program of international cooperative studies in scientific hydrology began about 4 years ago, and the IHD ( International Hydrological Decade ) will begin under international auspices in January 1965. This program will be highly important for hydrology and hydrologists, and it should contribute greatly to human welfare in the future. Fortuitously, the opening year coincides with International Cooperation Year of the United Nations Organization. American participation in the program will contribute to the purpose avowed by the President [Johnson, 1964] in a commencement address at Holy Cross early in June 1964: “to dedicate this year to finding new techniques for making Man's knowledge serve Man's welfare. Let this be a year of science.”

Eos, Transactions, American Geophysical Union↗

USGS2010

Multiobjective sampling design for parameter estimation and model discrimination in groundwater solute transport

Sampling design for site characterization studies of solute transport in porous media is formulated as a multiobjective problem. Optimal design of a sampling network is a sequential process in which the next phase of sampling is designed on the basis of all available physical knowledge of the system. Three objectives are considered: model discrimination, parameter estimation, and cost minimization. For the first two objectives, physically based measures of the value of information obtained from a set of observations are specified. In model discrimination, value of information of an observation point is measured in terms of the difference in solute concentration predicted by hypothesized models of transport. Points of greatest difference in predictions can contribute the most information to the discriminatory power of a sampling design. Sensitivity of solute concentration to a change in a parameter contributes information on the relative variance of a parameter estimate. Inclusion of points in a sampling design with high sensitivities to parameters tends to reduce variance in parameter estimates. Cost minimization accounts for both the capital cost of well installation and the operating costs of collection and analysis of field samples. Sensitivities, discrimination information, and well installation and sampling costs are used to form coefficients in the multiobjective problem in which the decision variables are binary (zero/one), each corresponding to the selection of an observation point in time and space. The solution to the multiobjective problem is a noninferior set of designs. To gain insight into effective design strategies, a one-dimensional solute transport problem is hypothesized. Then, an approximation of the noninferior set is found by enumerating 120 designs and evaluating objective functions for each of the designs. Trade-offs between pairs of objectives are demonstrated among the models. The value of an objective function for a given design is shown to correspond to the ability of a design to actually meet an objective.

Water Resources Research↗