Search USGS⌕ Search

SEARCH · Search USGS

Results for “Journal of Applied Mechanics”

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 91 records · Page 5Linked to original sources
USGS2020

Semiautomated estimates of directivity and related source properties of small to moderate southern California earthquakes using second seismic moments

We develop a semiautomated method for estimating with second seismic moments the directivity, rupture area, duration, and centroid velocity of earthquakes. The method is applied to 41 southern California earthquakes with magnitude in the range 3.5–5.2 and provides stable results for 28 events. Apparent source time functions (ASTFs) of P and S phases are derived using deconvolution with three stacked empirical Green's functions (seGf). The use of seGf suppresses nongeneric source effects, improves the focal mechanism correspondence to the analyzed earthquakes, and typically allows inclusion of 5 to 15 more ASTFs compared with analysis using a single eGf. Most analyzed earthquakes in the Trifurcation area of the San Jacinto Fault have directivities toward the northwest, while events around Cajon Pass and San Gabriel Mountain tend to propagate toward the southeast. These results are generally consistent with predictions for dynamic rupture on bimaterial interfaces associated with the imaged velocity contrasts in the area. The second moment inversions also provide constraints on the upper and lower bounds of rupture areas in our data set. Stress drops and uncertainties are estimated for elliptical ruptures using the derived characteristic rupture length and width. The semiautomated second moment method with seGfs can be used for routine application to moderate earthquakes in locations with good station coverage.

California↗
USGS2012

Implications of regional gravity for state of stress in the earth's crust and upper mantle

Topography is maintained by stress differences within the earth. Depending on the distribution of the stress we classify the support as either local or regional compensation. In general, the stresses implied in a regional compensation scheme are an order of magnitude larger than those corresponding to local isostasy. Gravity anomalies, a measure of the earth's departure from hydrostatic equilibrium, can be used to distinguish between the two compensation mechanisms and thus to estimate the magnitude of deviatoric stress in the crust and upper mantle. Topography created at an oceanic ridge crest or in a major continental orogenic zone appears to be locally compensated. Such features were formed on weak crust incapable of maintaining stress differences much greater than the stress from the applied load. Oceanic volcanoes formed on an already cooled, thickened lithosphere are regionally supported with elastic stresses. Similarly, the broad topographic rise seaward of subduction zones is elastically supported as the lithosphere is bent near the plate margin. Although the implied stress is to some degree dependent on the rheological model assumed, the gravity anomalies and surface deformation produced by these features demonstrate that the upper 30–40 km of the oceanic lithosphere is capable of regionally supporting stress differences in the 100-MPa range. Given certain conditions of load emplacement, continental crust can also support loads regionally over 100-m.y. time scales, but the effects of erosion only allow an estimate of a lower bound on stress. Data from space probes indicate that the upper layers of other terrestrial planets also support topographic-induced stress differences in excess of 100 MPa.

Journal of Geophysical Research Solid Earth↗
USGS2000

Numerical model of frazil ice and suspended sediment concentrations and formation of sediment laden ice in the Kara Sea

A one-dimensional (vertical) numerical model of currents, mixing, frazil ice concentration, and suspended sediment concentration has been developed and applied in the shallow southeastern Kara Sea. The objective of the calculations is to determine whether conditions suitable for turbid ice formation can occur during times of rapid cooling and wind- and wave-induced sediment resuspension. Although the model uses a simplistic approach to ice particles and neglects ice-sediment interactions, the results for low-stratification, shallow (∼20-m) freeze-up conditions indicate that the coconcentrations of frazil ice and suspended sediment in the water column are similar to observed concentrations of sediment in turbid ice. This suggests that wave-induced sediment resuspension is a viable mechanism for turbid ice formation, and enrichment mechanisms proposed to explain the high concentrations of sediment in turbid ice relative to sediment concentrations in underlying water may not be necessary in energetic conditions. However, salinity stratification found near the Ob' and Yenisey Rivers damps mixing between ice-laden surface water and sediment-laden bottom water and probably limits incorporation of resuspended sediment into turbid ice until prolonged or repeated wind events mix away the stratification. Sensitivity analyses indicate that shallow (≤20 m), unstratified waters with fine bottom sediment (settling speeds of ∼1 mm s −1 or less) and long open water fetches (>25 km) are ideal conditions for resuspension.

Kara Sea↗
USGS2017

An accurate method for measuring triploidy of larval fish spawns

A standard flow cytometric protocol was developed for estimating triploid induction in batches of larval fish. Polyploid induction treatments are not guaranteed to be 100% efficient, thus the ability to quantify the proportion of triploid larvae generated by a particular treatment helps managers to stock high-percentage spawns and researchers to select treatments for efficient triploid induction. At 3 d posthatch, individual Grass Carp Ctenopharyngodon idella were mechanically dissociated into single-cell suspensions; nuclear DNA was stained with propidium iodide then analyzed by flow cytometry. Following ploidy identification of individuals, aliquots of diploid and triploid cell suspensions were mixed to generate 15 levels (0–100%) of known triploidy ( n = 10). Using either 20 or 50 larvae per level, the observed triploid percentages were lower than the known, actual values. Using nonlinear regression analyses, quadratic equations solved for triploid proportions in mixed samples and corresponding estimation reference plots allowed for predicting triploidy. Thus, an accurate prediction of the proportion of triploids in a spawn can be made by following a standard larval processing and analysis protocol with either 20 or 50 larvae from a single spawn, coupled with applying the quadratic equations or reference plots to observed flow cytometry results. Due to the universality of triploid DNA content being 1.5 times the diploid level and because triploid fish consist of fewer cells than diploids, this method should be applicable to other produced triploid fish species, and it may be adapted for use with bivalves or other species where batch analysis is appropriate.

North American Journal of Aquaculture↗
USGS2002

Estimating recharge at Yucca Mountain, Nevada, USA: Comparison of methods

Obtaining values of net infiltration, groundwater travel time, and recharge is necessary at the Yucca Mountain site, Nevada, USA, in order to evaluate the expected performance of a potential repository as a containment system for high-level radioactive waste. However, the geologic complexities of this site, its low precipitation and net infiltration, with numerous mechanisms operating simultaneously to move water through the system, provide many challenges for the estimation of the spatial distribution of recharge. A variety of methods appropriate for arid environments has been applied, including water-balance techniques, calculations using Darcy's law in the unsaturated zone, a soil-physics method applied to neutron-hole water-content data, inverse modeling of thermal profiles in boreholes extending through the thick unsaturated zone, chloride mass balance, atmospheric radionuclides, and empirical approaches. These methods indicate that near-surface infiltration rates at Yucca Mountain are highly variable in time and space, with local (point) values ranging from zero to several hundred millimeters per year. Spatially distributed net-infiltration values average 5 mm/year, with the highest values approaching 20 mm/year near Yucca Crest. Site-scale recharge estimates range from less than 1 to about 12 mm/year. These results have been incorporated into a site-scale model that has been calibrated using these data sets that reflect infiltration processes acting on highly variable temporal and spatial scales. The modeling study predicts highly non-uniform recharge at the water table, distributed significantly differently from the non-uniform infiltration pattern at the surface.

Hydrogeology Journal↗
USGS2018

Mechanisms of wave‐driven water level variability on reef‐fringed coastlines

Wave‐driven water level variability (and runup at the shoreline) is a significant cause of coastal flooding induced by storms. Wave runup is challenging to predict, particularly along tropical coral reef‐fringed coastlines due to the steep bathymetric profiles and large bottom roughness generated by reef organisms, which can violate assumptions in conventional models applied to open sandy coastlines. To investigate the mechanisms of wave‐driven water level variability on a reef‐fringed coastline, we performed a set of laboratory flume experiments on an along‐shore uniform bathymetric profile with and without bottom roughness. Wave setup and waves at frequencies lower than the incident sea‐swell forcing (infragravity waves) were found to be the dominant components of runup. These infragravity waves were positively correlated with offshore wave groups, signifying they were generated in the surf zone by the oscillation of the breakpoint. On the reef flat and at the shoreline, the low‐frequency waves formed a standing wave pattern with energy concentrated at the natural frequencies of the reef flat, indicating resonant amplification. Roughness elements used in the flume to mimic large reef bottom roughness reduced low frequency motions on the reef flat and reduced wave run up by 30% on average, compared to the runs over a smooth bed. These results provide insight into sea‐swell and infragravity wave transformation and wave setup dynamics on steep‐sloped coastlines, and the effect that future losses of reef bottom roughness may have on coastal flooding along reef‐fringed coasts.

Journal of Geophysical Research C: Oceans↗
USGS2002

Waveform inversion of oscillatory signatures in long-period events beneath volcanoes

The source mechanism of long-period (LP) events is examined using synthetic waveforms generated by the acoustic resonance of a fluid-filled crack. We perform a series of numerical tests in which the oscillatory signatures of synthetic LP waveforms are used to determine the source time functions of the six moment tensor components from waveform inversions assuming a point source. The results indicate that the moment tensor representation is valid for the odd modes of crack resonance with wavelengths 2 L / n , 2 W / n , n = 3, 5, 7, …, where L and W are the crack length and width, respectively. For the even modes with wavelengths 2 L / n , 2 W / n , n = 2, 4, 6, …, a generalized source representation using higher-order tensors is required, although the efficiency of seismic waves radiated by the even modes is expected to be small. We apply the moment tensor inversion to the oscillatory signatures of an LP event observed at Kusatsu-Shirane Volcano, central Japan. Our results point to the resonance of a subhorizontal crack located a few hundred meters beneath the summit crater lakes. The present approach may be useful to quantify the source location, geometry, and force system of LP events, and opens the way for moment tensor inversions of tremor.

Journal of Geophysical Research B: Solid Earth↗
USGS2016

Three-dimensional surface deformation derived from airborne interferometric UAVSAR: Application to the Slumgullion Landslide

In order to provide surface geodetic measurements with “landslide-wide” spatial coverage, we develop and validate a method for the characterization of 3-D surface deformation using the unique capabilities of the Uninhabited Aerial Vehicle Synthetic Aperture Radar (UAVSAR) airborne repeat-pass radar interferometry system. We apply our method at the well-studied Slumgullion Landslide, which is 3.9 km long and moves persistently at rates up to ∼2 cm/day. A comparison with concurrent GPS measurements validates this method and shows that it provides reliable and accurate 3-D surface deformation measurements. The UAVSAR-derived vector velocity field measurements accurately capture the sharp boundaries defining previously identified kinematic units and geomorphic domains within the landslide. We acquired data across the landslide during spring and summer and identify that the landslide moves more slowly during summer except at its head, presumably in response to spatiotemporal variations in snowmelt infiltration. In order to constrain the mechanics controlling landslide motion from surface velocity measurements, we present an inversion framework for the extraction of slide thickness and basal geometry from dense 3-D surface velocity fields. We find that the average depth of the Slumgullion Landslide is 7.5 m, several meters less than previous depth estimates. We show that by considering a viscoplastic rheology, we can derive tighter theoretical bounds on the rheological parameter relating mean horizontal flow rate to surface velocity. Using inclinometer data for slow-moving, clay-rich landslides across the globe, we find a consistent value for the rheological parameter of 0.85 ± 0.08.

Colorado↗
USGS2019

Effects of infiltration characteristics on the spatial-temporal evolution of stability of an interstate highway embankment

Infiltration-induced landslides are among the most common natural disasters threatening modern civilization, but conventional methods for studying the triggering mechanisms and predicting the occurrence of these slides are limited by incomplete consideration of underlying physical processes and the lack of precision inherent in limit-equilibrium analyses. To address this problem the spatial-temporal evolution of failure is investigated in a seasonally unstable section of interstate highway embankment, known as the Straight Creek landslide, Colorado. The study includes multi-year site investigation, monitoring, and numerical simulation using a rigorous hydromechanical framework along with a field of local factor of safety method. The sensitivity of episodic landslide reactivation to infiltration characteristics is evaluated. Results indicate that annual cumulative snowmelt infiltration, which typically accounts for approximately 75% of total annual cumulative infiltration and occurs over a short period in the spring, has the most substantial impact on slide activation. The rate of snowmelt infiltration varies independently of annual cumulative snowmelt infiltration and cumulative infiltration in the previous year, but still affects antecedent soil moisture conditions at the onset of snowmelt infiltration and therefore also the level of slide activation. These findings are used to establish specific thresholds for exacerbated slide movement using annual snowpack accumulation, forecasted snowmelt rate, and the previous year’s snowmelt, an approach which may be applied for predicting movement at this and other recurring or potential slide sites.

Colorado↗
USGS2014

Mechanisms driving recruitment variability in fish: comparisons between the Laurentian Great Lakes and marine systems

In his seminal work, Hjort (in Fluctuations in the great fisheries of Northern Europe. Conseil Parmanent International Pour L'Exploration De La Mar. Rapports et Proces-Verbaux, 20: 1–228, 1914) observed that fish population levels fluctuated widely, year-class strength was set early in life, and egg production by adults could not alone explain variability in year-class strength. These observations laid the foundation for hypotheses on mechanisms driving recruitment variability in marine systems. More recently, researchers have sought to explain year-class strength of important fish in the Laurentian Great Lakes and some of the hypotheses developed for marine fisheries have been transferred to Great Lakes fish. We conducted a literature review to determine the applicability of marine recruitment hypotheses to Great Lakes fish. We found that temperature, interspecific interactions, and spawner effects (abundance, age, and condition of adults) were the most important factors in explaining recruitment variability in Great Lakes fish, whereas relatively fewer studies identified bottom-up trophodynamic factors or hydrodynamic factors as important. Next, we compared recruitment between Great Lakes and Baltic Sea fish populations and found no statistical difference in factors driving recruitment between the two systems, indicating that recruitment hypotheses may often be transferable between Great Lakes and marine systems. Many recruitment hypotheses developed for marine fish have yet to be applied to Great Lakes fish. We suggest that future research on recruitment in the Great Lakes should focus on forecasting the effects of climate change and invasive species. Further, because the Great Lakes are smaller and more enclosed than marine systems, and have abundant fishery-independent data, they are excellent candidates for future hypothesis testing on recruitment in fish.

Laurentian Great Lakes↗
USGS2002

Topographic stress perturbations in southern Davis Mountains, west Texas 1. Polarity reversal of principal stresses

We have applied a previously developed analytical stress model to interpret subsurface stress conditions inferred from acoustic televiewer logs obtained in two municipal water wells located in a valley in the southern Davis Mountains near Alpine, Texas. The appearance of stress-induced breakouts with orientations that shift by 90° at two different depths in one of the wells is explained by results from exact solutions for the effects of valleys on gravity and tectonically induced subsurface stresses. The theoretical results demonstrate that above a reference depth termed the hinge point, a location that is dependent on Poisson's ratio, valley shape, and magnitude of the maximum horizontal tectonic stress normal to the long axis of the valley, horizontal stresses parallel to the valley axis are greater than those normal to it. At depths below this hinge point the situation reverses and horizontal stresses normal to the valley axis are greater than those parallel to it. Application of the theoretical model at Alpine is accommodated by the fact that nearby earthquake focal mechanisms establish an extensional stress regime with the regional maximum horizontal principal stress aligned perpendicular to the valley axis. We conclude that the localized stress field associated with a valley setting can be highly variable and that breakouts need to be examined in this context when estimating the orientations and magnitudes of regional principal stresses.

Texas↗
USGS2015

A new method for discovering behavior patterns among animal movements

Advanced satellite tracking technologies enable biologists to track animal movements at fine spatial and temporal scales. The resultant data present opportunities and challenges for understanding animal behavioral mechanisms. In this paper, we develop a new method to elucidate animal movement patterns from tracking data. Here, we propose the notion of continuous behavior patterns as a concise representation of popular migration routes and underlying sequential behaviors during migration. Each stage in the pattern is characterized in terms of space (i.e., the places traversed during movements) and time (i.e. the time spent in those places); that is, the behavioral state corresponding to a stage is inferred according to the spatiotemporal and sequential context. Hence, the pattern may be interpreted predictably. We develop a candidate generation and refinement framework to derive all continuous behavior patterns from raw trajectories. In the framework, we first define the representative spots to denote the underlying potential behavioral states that are extracted from individual trajectories according to the similarity of relaxed continuous locations in certain distinct time intervals. We determine the common behaviors of multiple individuals according to the spatiotemporal proximity of representative spots and apply a projection-based extension approach to generate candidate sequential behavior sequences as candidate patterns. Finally, the candidate generation procedure is combined with a refinement procedure to derive continuous behavior patterns. We apply an ordered processing strategy to accelerate candidate refinement. The proposed patterns and discovery framework are evaluated through conceptual experiments on both real GPS-tracking and large synthetic datasets.

International Journal of Geographical Information ↗
USGS2000

Stress resistance strategy in an arid land shrub: interactions between developmental instability and fractal dimention

This paper investigates allocation of energy to mechanisms that generate and preserve architectural forms (i.e. developmental stability, complexity of branching patterns) and productivity (growth and reproduction) in response to environmental disturbances (i.e. grazing and resource availability). The statistical error in translational symmetry was used to detect random intra-individual variability during development. This can be thought of as a measure of developmental instability caused by stress. Additionally, we use changes in fractal complexity and shoot distribution of branch structures as an alternate indicator of stress. These methods were applied to Anthyllis cytisoides L., a semi-arid environment shrub, to ascertain the effect of grazing and slope exposure on developmental traits in a 2×2 factorial design. The results show that A. cytisoides maintains developmental stability at the expense of productivity. Anthyllis cytisoides was developmentally more stable when grazed and when on south-facing, as opposed to north-facing slopes. On the contrary, shoot length, leaf area, fractal dimension and reproductive-to-vegetative allocation ratio were larger in north- than in south-facing slopes. As a consequence, under extreme xeric conditions, shrub mortality increased in north-facing slopes, especially when not grazed. The removal of transpiring area and the reduction of plant competition favoured developmental stability and survival in grazed plants. Differences between grazed and ungrazed plants were most evident in more mesic (north-facing) areas.

Journal of Fish Diseases↗
USGS2020

Conceptualizing alternate regimes in a large floodplain-river ecosystem

Regime shifts –persistent changes in the structure and function of an ecosystem - are well-documented in many ecosystems but remain poorly understood in floodplain-river ecosystems. We apply a resilience perspective to large floodplain-river ecosystems by presenting three examples of plausible sets of alternate regimes that are relevant to natural resource management interests within the Upper Mississippi River and Illinois River. These alternate regimes include: 1) a clear water and abundant vegetation regime vs. a turbid water and sparse vegetation regime in lentic, off-channel areas, 2) a diverse native fish community regime vs. an invasive-dominated fish community regime, and 3) a regime characterized by a diverse and dynamic mosaic of floodplain vegetation types vs. one characterized as a persistent invasive wet meadow monoculture. For each set of potential alternate regimes, we synthesize known or hypothesized feedback mechanisms that reinforce regimes, controlling variables that drive regime transitions, and restoration pathways. The conceptual models presented here provide a framework for synthesizing our understanding of the dynamics of this ecosystem and are relevant to other large floodplain-river ecosystems that face similar human pressures across the world. The models are currently being used to prioritize future research, test hypotheses, and inform restoration and management on the Upper Mississippi River and Illinois River. Through sharing our approach, we provide a case study in which we document an important step in operationalizing resilience concepts for the management of natural resources.

Minnesota, Wisconsin, Iowa, Illinois, Missouri↗
USGS1985

Three-dimensional trend mapping from wire-line logs

Mapping of lithofacies and porosities of stratigraphic units is complicated because these properties vary in three dimensions. The method of moments was proposed by Krumbein and Libby (1957) as a technique to aid in resolving this problem. Moments are easily computed from wireline logs and are simple statistics which summarize vertical variation in a log trace. Combinations of moment maps have proved useful in understanding vertical and lateral changes in lithology of sedimentary rock units. Although moments have meaning both as statistical descriptors and as mechanical properties, they also define polynomial curves which approximate lithologic changes as a function of depth. These polynomials can be fitted by least-squares methods, partitioning major trends in rock properties from finescale fluctuations. Analysis of variance yields the degree of fit of any polynomial and measures the proportion of vertical variability expressed by any moment or combination of moments. In addition, polynomial curves can be differentiated to determine depths at which pronounced expressions of facies occur and to determine the locations of boundaries between major lithologic subdivisions. Moments can be estimated at any location in an area by interpolating from log moments at control wells. A matrix algebra operation then converts moment estimates to coefficients of a polynomial function which describes a continuous curve of lithologic variation with depth. If this procedure is applied to a grid of geographic locations, the result is a model of variability in three dimensions. Resolution of the model is determined largely by number of moments used in its generation. The method is illustrated with an analysis of lithofacies in the Simpson Group of south-central Kansas; the three-dimensional model is shown as cross sections and slice maps. In this study, the gamma-ray log is used as a measure of shaliness of the unit. However, the method is general and can be applied, for example, to suites of neutron, density, or sonic logs to produce three-dimensional models of porosity in reservoir rocks. ?? 1985 Plenum Publishing Corporation.

Journal of the International Association for Mathe↗
USGS2025

Concentration-discharge relations and transient metal loads reveal spatiotemporal variability in solute-generation mechanisms in a mine-affected watershed

Concentration-discharge (CQ) relations are commonly used to understand geochemical and hydrologic controls on the generation of solutes in watersheds. Despite the widespread application of CQ relations, this technique has been infrequently applied to acid mine drainage (AMD) sites, but the CQ framework may allow mechanistic understanding of remedial outcomes such as impoundment of water within underground mines. Results of CQ analyses and changes in metal loads in an AMD affected watershed in Colorado, USA indicate that dissolved loads increased at many individual locations following water impoundment within mine workings. Although increased loads were observed at most individual locations, these increases were offset by a large decrease in loading from the largest mine. A loading analysis that included data from an instream monitoring location showed a statistically significant decrease in Fe and Zn after bulkhead emplacement, indicating a net positive effect of bulkheads. Streams generally displayed dilution CQ patterns whereas mines and springs showed either flushing or chemostatic patterns prior to bulkheading, which transitioned to chemostatic patterns following bulkheading, indicating a transition from dynamic to equilibrium geochemical processes. Saturation indices for sulfide and secondary minerals indicated that mines and springs were near equilibrium for phases including schwertmannite, fluorite, and gypsum. Saturation indices vary through time for mines suggesting progressive leaching of sulfide minerals as the mass of available minerals in the mine workings decreases. Together, these diverse analyses provide an integrated understanding of the variability in solute generating processes in this watershed and may inform remediation plans for similarly affected sites by indicating the nature of mineralogic controls on water quality.

Colorado↗
USGS2003

Application of two hydrologic models with different runoff mechanisms to a hillslope dominated watershed in the northeastern US: A comparison of HSPF and SMR

Differences in the simulation of hydrologic processes by watershed models directly affect the accuracy of results. Surface runoff generation can be simulated as either: (1) infiltration-excess (or Hortonian) overland flow, or (2) saturation-excess overland flow. This study compared the Hydrological Simulation Program - FORTRAN (HSPF) and the Soil Moisture Routing (SMR) models, each representing one of these mechanisms. These two models were applied to a 102 km2 watershed in the upper part of the Irondequoit Creek basin in central New York State over a seven-year simulation period. The models differed in both the complexity of simulating snowmelt and baseflow processes as well as the detail in which the geographic information was preserved by each model. Despite their differences in structure and representation of hydrologic processes, the two models simulated streamflow with almost equal accuracy. Since streamflow is an integral response and depends mainly on the watershed water balance, this was not unexpected. Model efficiency values for the seven-year simulation period were 0.67 and 0.65 for SMR and HSPF, respectively. HSPF simulated winter streamflow slightly better than SMR as a result of its complex snowmelt routine, whereas SMR simulated summer flows better than HSPF as a result of its runoff and baseflow processes. An important difference between model results was the ability to predict the spatial distribution of soil moisture content. HSPF aggregates soil moisture content, which is generally related to a specific pervious land unit across the entire watershed, whereas SMR predictions of moisture content distribution are geographically specific and matched field observations reasonably well. Important is that the saturated area was predicted well by SMR and confirmed the validity of using saturation-excess mechanisms for this hillslope dominated watershed. ?? 2003 Elsevier B.V. All rights reserved.

Journal of Hydrology↗
USGS2003

Automated calibration of a stream solute transport model: Implications for interpretation of biogeochemical parameters

The hydrologic processes of advection, dispersion, and transient storage are the primary physical mechanisms affecting solute transport in streams. The estimation of parameters for a conservative solute transport model is an essential step to characterize transient storage and other physical features that cannot be directly measured, and often is a preliminary step in the study of reactive solutes. Our study used inverse modeling to estimate parameters of the transient storage model OTIS (One dimensional Transport with Inflow and Storage). Observations from a tracer injection experiment performed on Uvas Creek, California, USA, are used to illustrate the application of automated solute transport model calibration to conservative and nonconservative stream solute transport. A computer code for universal inverse modeling (UCODE) is used for the calibrations. Results of this procedure are compared with a previous study that used a trial-and-error parameter estimation approach. The results demonstrated 1) importance of the proper estimation of discharge and lateral inflow within the stream system; 2) that although the fit of the observations is not much better when transient storage is invoked, a more randomly distributed set of residuals resulted (suggesting non-systematic error), indicating that transient storage is occurring; 3) that inclusion of transient storage for a reactive solute (Sr2+) provided a better fit to the observations, highlighting the importance of robust model parameterization; and 4) that applying an automated calibration inverse modeling estimation approach resulted in a comprehensive understanding of the model results and the limitation of input data.

Journal of the North American Benthological Societ↗