SEARCH · USGS Science

Results for “Solutions Journal”

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 451 records · Page 25Linked to original sources

USGS1974

Stabilities of calcite and aragonite

A revaluation of the 25° C activity-product constants of calcite ( K C ) and aragonite ( K A ) was made on the basis of the known solubilities of these phases for which the activity of total dissolved calcium was corrected for the presence of the ion pair CaHCO 3 + in the aqueous phase. The value of the dissociation constant of CaHCO 3 + was taken to be 10 -1.225±0.02 . This value, combined with values of the analytical concentrations in solutions with partial pressure P CO2 =0.97 atmosphere, leads to K C =l0 -8.52±0.04 and K A = 10 -8.36±0.04 . Based on these K values, standard free energies of formation of calcite and aragonite were calculated to be -270,144±375 and -269,926±375 calories mole -1 , (-1,130,282±1,569 and -1,129,370±1,569 joules mole -1 ), respectively. From the 25°C K values, using appropriate entropy and heat capacity data, values of K C and K A were calculated over the temperature range 0° to 200°C. Possible errors in interpretation of measured pH values and inferred P CO2 values and the bearing of these errors on calculations of K values are discussed.

Journal of Research of the U.S. Geological Survey↗

USGS1978

Models for calculating density and vapor pressure of geothermal brines

In a model for estimating density of a brine, the density of a natural brine at a known temperature, pressure, and composition can be calculated from the densities of the component salt solutions in the complex brine. A model for estimating vapor pressure requires two direct or indirect estimates of the vapor pressure and uses a published "reference-substance principle" to extrapolate,; these data 200° to 300° above the temperature of the highest observation. The methods can be used to estimate the partial molal volume and partial molal enthalpy of H 2 O in the natural brine and the partial molal volume change and partial molal heat change for the H2O component in the vaporization process. The equations for density and vapor pressure and the derived thermochemical properties will help in designing turbines for operation with geothermal and geopressured brines, in optimizing production conditions, and in physical and chemical modeling of a geothermal reservoir.

Journal of Research of the U.S. Geological Survey↗

USGS1970

Comparison of chemical hydrogeology of the carbonate peninsulas of Florida and Yucatan

Aquifers of the peninsulas of Florida and northern Yucatan are Tertiary marine carbonate formations showing many lithologic and faunal similarities. In addition, the tropical to subtropical climates of the two areas are similar, each having annual rainfall of about 1000 to 1500 mm. Despite similarities in these fundamental controls, contrasts in the hydrologic and geochemical systems are numerous and striking. For example, Florida has many rivers; Yucatan has none. Maximum thickness of fresh ground water in Florida is about 700 meters; in the Yucatan it is less than 70 meters. In Florida the gradient of the potentiometric surface averages about 1 meter per kilometer; in the Yucatan it is exceedingly low, averaging about 0.02 meter per kilometer. In Florida the chemical character of water changes systematically downgradient, owing to solution of minerals of the aquifer and corresponding increases in total dissolved solids, sulfate, calcium, and Mg-Ca ratio; in the Yucatan no downgradient change exists, and dominant processes controlling the chemical character of the water are solution of minerals and simple mixing of the fresh water and the body of salt water that underlies the peninsula at shallow depth. Hydrologic and chemical differences are caused in part by the lower altitude of the Yucatan plain. More important, however, these differences are due to the lack of an upper confining bed in Yucatan that is hydrologically equivalent to the Hawthorn Formation of Florida. The Hawthorn cover prevents recharge and confines the artesian water except where it is punctured by sinkholes, but sands and other unconsolidated sediments fill sinkholes and cavities and impede circulation. In the Yucatan the permeability of the entire section is so enormous that rainfall immediately infiltrates to the water table and then moves laterally to discharge areas along the coasts.

Florida Peninsula, Yucatan Peninsula↗

USGS2008

GSTARS computer models and their applications, part I: theoretical development

GSTARS is a series of computer models developed by the U.S. Bureau of Reclamation for alluvial river and reservoir sedimentation studies while the authors were employed by that agency. The first version of GSTARS was released in 1986 using Fortran IV for mainframe computers. GSTARS 2.0 was released in 1998 for personal computer application with most of the code in the original GSTARS revised, improved, and expanded using Fortran IV/77. GSTARS 2.1 is an improved and revised GSTARS 2.0 with graphical user interface. The unique features of all GSTARS models are the conjunctive use of the stream tube concept and of the minimum stream power theory. The application of minimum stream power theory allows the determination of optimum channel geometry with variable channel width and cross-sectional shape. The use of the stream tube concept enables the simulation of river hydraulics using one-dimensional numerical solutions to obtain a semi-two- dimensional presentation of the hydraulic conditions along and across an alluvial channel. According to the stream tube concept, no water or sediment particles can cross the walls of stream tubes, which is valid for many natural rivers. At and near sharp bends, however, sediment particles may cross the boundaries of stream tubes. GSTARS3, based on FORTRAN 90/95, addresses this phenomenon and further expands the capabilities of GSTARS 2.1 for cohesive and non-cohesive sediment transport in rivers and reservoirs. This paper presents the concepts, methods, and techniques used to develop the GSTARS series of computer models, especially GSTARS3. ?? 2008 International Research and Training Centre on Erosion and Sedimentation and the World Association for Sedimentation and Erosion Research.

International Journal of Sediment Research↗

USGS1998

Closing the gap between regional and global travel time tomography

Recent global travel time tomography studies by Zhou [1996] and van der Hilst et al . [1997] have been performed with cell parameterizations of the order of those frequently used in regional tomography studies (i.e., with cell sizes of 1°–2°). These new global models constitute a considerable improvement over previous results that were obtained with rather coarse parameterizations (5° cells). The inferred structures are, however, of larger scale than is usually obtained in regional models, and it is not clear where and if individual cells are actually resolved. This study aims at resolving lateral heterogeneity on scales as small as 0.6° in the upper mantle and 1.2°–3° in the lower mantle. This allows for the adequate mapping of expected small-scale structures induced by, for example, lithosphere subduction, deep mantle upwellings, and mid-ocean ridges. There are three major contributions that allow for this advancement. First, we employ an irregular grid of nonoverlapping cells adapted to the heterogeneous sampling of the Earth's mantle by seismic waves [ Spakman and Bijwaard , 1998]. Second, we exploit the global data set of Engdahl et al . [1998], which is a reprocessed version of the global data set of the International Seismological Centre. Their reprocessing included hypocenter redetermination and phase reidentification. Finally, we combine all data used ( P , pP , and pwP phases) into nearly 5 million ray bundles with a limited spatial extent such that averaging over large mantle volumes is prevented while the signal-to-noise ratio is improved. In the approximate solution of the huge inverse problem we obtain a variance reduction of 57.1%. Synthetic sensitivity tests indicate horizontal resolution on the scale of the smallest cells (0.6° or 1.2°) in the shallow parts of subduction zones decreasing to approximately 2°–3° resolution in well-sampled regions in the lower mantle. Vertical resolution can be worse (up to several hundreds of kilometers) in subduction zones with rays predominantly pointing along dip. Important features of the solution are as follows: 100–200 km thick high-velocity slabs beneath all major subduction zones, sometimes flattening in the transition zone and sometimes directly penetrating into the lower mantle; large high-velocity anomalies in the lower mantle that have been attributed to subduction of the Tethys ocean and the Farallon plate; and low-velocity anomalies continuing across the 660 km discontinuity to hotspots at the surface under Iceland, east Africa, the Canary Islands, Yellowstone, and the Society Islands. Our findings corroborate that the 660 km boundary may resist but not prevent (present day) large-scale mass transfer from upper to lower mantle or vice versa. This observation confirms the results of previous, global mantle studies that employed coarser parameterizations.

Journal of Geophysical Research B: Solid Earth↗

USGS1987

Origins of seawater intrusion in a coastal aquifer - A case study of the Pajaro Valley, California

Seawater may enter and contaminate stratified coastal aquifers through a number of different pathways. These pathways and their relative contribution are examined in the Pajaro Valley, California, a coastal area with extensive groundwater development. This study considers three pathways of possible intrusion of the primary confined aquifer: (1) onshore leakage from brackish sources, the estuary and sloughs, through the confining layer; (2) near-shore leakage from the ocean through the confining layer; and (3) offshore flow from the ocean through the submarine canyon outcrop of the aquifer. Groundwater flow and seawater intrusion are simulated using an areal, two-dimensional solute-transport computer model. This analysis indicates that leakage through confining layers is the principal mechanism of recharge to the aquifer. Although lateral flow through the offshore outcrop contaminates the aquifer, as a whole, at a higher rate, vertical leakage through the sea floor initially is the main pathway of seawater intrusion to the onshore portion of the aquifer. It is likely that leakage generally is the dominant mechanism of recharge and initial cause of seawater intrusion for poorly-confined, stratified coastal aquifers. This analysis suggests that a significant time interval follows the initial observation of seawater intrusion, during which remedial action can be taken to control lateral flow through the offshore outcrop, which ultimately will be the largest component of future intrusion in these aquifers.

California↗

USGS2000

Mineralogic and compositional properties of Martian soil and dust: results from Mars Pathfinder

Mars Pathfinder obtained multispectral, elemental, magnetic, and physical measurements of soil and dust at the Sagan Memorial Station during the course of its 83 sol mission. We describe initial results from these measurements, concentrating on multispectral and elemental data, and use these data, along with previous Viking, SNC meteorite, and telescopic results, to help constrain the origin and evolution of Martian soil and dust. We find that soils and dust can be divided into at least eight distinct spectral units, based on parameterization of Imager for Mars Pathfinder (IMP) 400 to 1000 nm multispectral images. The most distinctive spectral parameters for soils and dust are the reflectivity in the red, the red/blue reflectivity ratio, the near-IR spectral slope, and the strength of the 800 to 1000 nm absorption feature. Most of the Pathfinder spectra are consistent with the presence of poorly crystalline or nanophase ferric oxide(s), sometimes mixed with small but varying degrees of well-crystalline ferric and ferrous phases. Darker soil units appear to be coarser-grained, compacted, and/or mixed with a larger amount of dark ferrous materials relative to bright soils. Nanophase goethite, akaganeite, schwertmannite, and maghemite are leading candidates for the origin of the absorption centered near 900 nm in IMP spectra. The ferrous component in the soil cannot be well-constrained based on IMP data. Alpha proton X-ray spectrometer (APXS) measurements of six soil units show little variability within the landing site and show remarkable overall similarity to the average Viking-derived soil elemental composition. Differences exist between Viking and Pathfinder soils, however, including significantly higher S and Cl abundances and lower Si abundances in Viking soils and the lack of a correlation between Ti and Fe in Pathfinder soils. No significant linear correlations were observed between IMP spectral properties and APXS elemental chemistry. Attempts at constraining the mineralogy of soils and dust using normative calculations involving mixtures of smectites and silicate and oxide minerals did not yield physically acceptable solutions. We attempted to use the Pathfinder results to constrain a number of putative soil and dust formation scenarios, including palagonitization and acid-fog weathering. While the Pathfinder soils cannot be chemically linked to the Pathfinder rocks by palagonitization, this study and McSween et al. [1999] suggest that palagonitic alteration of a Martian basaltic rock, plus mixture with a minor component of locally derived andesitic rock fragments, could be consistent with the observed soil APXS and IMP properties.

Journal of Geophysical Research E: Planets↗

USGS2022

The potential of wave energy conversion to mitigate coastal erosion from hurricanes

Wave energy conversion technologies have recently attracted more attention as part of global efforts to replace fossil fuels with renewable energy resources. While ocean waves can provide renewable energy, they can also be destructive to coastal areas that are often densely populated and vulnerable to coastal erosion. There have been a variety of efforts to mitigate the impacts of wave- and storm-induced erosion; however, they are either temporary solutions or approaches that are not able to adapt to a changing climate. This study explores a green and sustainable approach to mitigating coastal erosion from hurricanes through wave energy conversion. A barrier island, Dauphin Island, off the coast of Alabama, is used as a test case. The potential use of wave energy converter farms to mitigate erosion due to hurricane storm surges while simultaneously generating renewable energy is explored through simulations that are forced with storm data using the XBeach model. It is shown that wave farms can impact coastal morphodynamics and have the potential to reduce dune and beach erosion, predominantly in the western portion of the island. The capacity of wave farms to influence coastal morphodynamics varies with the storm intensity.

Alabama↗

USGS2003

Post-seismic relaxation theory on laterally heterogeneous viscoelastic model

Investigation was carried out into the problem of relaxation of a laterally heterogeneous viscoelastic Earth following an impulsive moment release event. The formal solution utilizes a semi-analytic solution for post-seismic deformation on a laterally homogeneous Earth constructed from viscoelastic normal modes, followed by application of mode coupling theory to derive the response on the aspherical Earth. The solution is constructed in the Laplace transform domain using the correspondence principle and is valid for any linear constitutive relationship between stress and strain. The specific implementation described in this paper is a semi-analytic discretization method which assumes isotropic elastic structure and a Maxwell constitutive relation. It accounts for viscoelastic-gravitational coupling under lateral variations in elastic parameters and viscosity. For a given viscoelastic structure and minimum wavelength scale, the computational effort involved with the numerical algorithm is proportional to the volume of the laterally heterogeneous region. Examples are presented of the calculation of post-seismic relaxation with a shallow, laterally heterogeneous volume following synthetic impulsive seismic events, and they illustrate the potentially large effect of regional 3-D heterogeneities on regional deformation patterns.

Geophysical Journal International↗

USGS1973

Concepts of karst development in relation to interpretation of surface runoff

Some unusual characteristics of streamflow occur in regions underlain by carbonate rocks. The streamflow characteristics are related to processes of karstification, these processes being dependent on circulation of subsurface water and solution of the rock to form characteristic topography and underground cavern systems. Very highly cavernous and permeable unsaturated zones tend to keep the water table depressed below land surface in many karst regions, a condition that leads to a low density of perennial streams. The uneven distribution of permeability beneath surface karst streams causes them to lose or gain water, depending on the position of the water table with reference to stream level. The conventional techniques of interpolation and extrapolation that have been reasonably successful in approximating streamflow of ungaged sites in nonkarstic regions have only limited use in karst regions. An understanding of principles of karstification and an understanding of the hydrogeologic framework of a carbonate terrane provide a useful basis for evaluating the streamflow characteristics.

Journal of Research of the U.S. Geological Survey↗

USGS1999

Climatic impact of glacial cycle polar motion: Coupled oscillations of ice sheet mass and rotation pole position

Precessional motion of Earth's rotation axis relative to its orbit is a well-known source of long-period climatic variation. It is less well appreciated that growth and decay of polar ice sheets perturb the symmetry of the global mass distribution enough that the geographic location of the rotation axis will change by at least 15 km and possibly as much as 100 km during a single glacial cycle. This motion of the pole will change the seasonal and latitudinal pattern of temperatures. We present calculations, based on a diurnal average energy balance, which compare the summer and winter temperature anomalies due to a 1° decrease in obliquity with those due to a 1° motion of the rotation pole toward Hudson Bay. Both effects result in peak temperature perturbations of about 1° Celsius. The obliquity change primarily influences the amplitude of the seasonal cycle, while the polar motion primarily changes the annual mean temperatures. The polar motion induced temperature anomaly is such that it will act as a powerful negative feedback on ice sheet growth. We also explore the evolution of the coupled system composed of ice sheet mass and pole position. Oscillatory solutions result from the conflicting constraints of rotational and thermal stability. A positive mass anomaly on an otherwise featureless Earth is in rotational equilibrium only at the poles or the equator. The two polar equilibria are rotationally unstable, and the equatorial equilibrium, though rotationally stable, is thermally unstable. We find that with a plausible choice for the strength of coupling between the thermal and rotational systems, relatively modest external forcing can produce significant response at periods of 10 4 –10 6 years, but it strongly attenuates polar motion at longer periods. We suggest that these coupled oscillations may contribute to the observed dominance of 100 kyr glacial cycles since the mid-Pleistocene and will tend to stabilize geographic patterns that are suitable to glaciations.

Journal of Geophysical Research B: Solid Earth↗

USGS1987

Toxicity of TFM lampricide to early life stages of walleye

We studied the effects of the lampricide 3-trifluoromethyl-4-nitrophenol (TFM) on gametes, newly fertilized eggs, eyed eggs, larvae, and swim-up fry of the walleye Stizostedion vitreum . When gametes from sexually mature walleyes were stripped into solutions of TFM, no effects were observed during the fertilization process at concentrations up to 3.0 mg/L—three times the concentration lethal to 99.9% of larval sea lampreys Petromyzon marinus held 12 h (LC99.9) under the same test conditions. Newly fertilized eggs likewise were unaffected during water hardening by concentrations of TFM that were lethal to sea lamprey ammocoetes. Eyed eggs, sac fry, and swim-up fry yielded LC25 values that were 2.5 to 5 times greater than the 12-h LC99.9 for sea lamprey ammocoetes. The data thus indicated that all of the early life stages of walleyes tested were considerably more resistant than sea lamprey ammocoetes to TFM, and that it is unlikely they would be adversely affected by standard stream treatments to kill sea lamprey ammocoetes.

North American Journal of Fisheries Management↗

USGS1990

Solid-solution aqueous-solution equilibria: Thermodynamic theory and representation

Thorstenson and Plummer's (1977) "stoichiometric saturation' model is reviewed, and a general relation between stoichiometric saturation Kss constants and excess free energies of mixing is derived for a binary solid-solution B1-xCxA: GE = RT[ln Kss - xln(xKCA) - (l-x)ln((l-x)KBA)]. This equation allows a suitable excess free energy function, such as Guggenheim's (1937) sub-regular function, to be fitted from experimentally determined Kss constants. Solid-phase free energies and component activity-coefficients can then be determined from one or two fitted parameters and from the endmember solubility products KBA and KCA. A general form of Lippmann's (1977,1980) "solutus equation is derived from an examination of Lippmann's (1977,1980) "total solubility product' model. Lippmann's II or "total solubility product' variable is used to represent graphically not only thermodynamic equilibrium states and primary saturation states but also stoichiometric saturation and pure phase saturation states.

American Journal of Science↗

USGS2009

Offshore double-planed shallow seismic zone in the NE Japan forearc region revealed by sP depth phases recorded by regional networks

We detected the sP depth phase at small epicentral distances of about 150 km or more in the seismograms of shallow earthquakes in the NE Japan forearc region. The focal depths of 1078 M > 3 earthquakes that occurred from 2000 to 2006 were precisely determined using the time delay of the sP phase from the initial P-wave arrival. The distribution of relocated hypocentres clearly shows the configuration of a double-planed shallow seismic zone beneath the Pacific Ocean. The upper plane has a low dip angle near the Japan Trench, increasing gradually to ???30?? at approximately 100 km landward of the Japan Trench. The lower plane is approximately parallel to the upper plane, and appears to be the near-trench counterpart of the lower plane of the double-planed deep seismic zone beneath the land area. The distance between the upper and lower planes is 28-32 km, which is approximately the same as or slightly smaller than that of the double-planed deep seismic zone beneath the land area. Focal mechanism solutions of the relocated earthquakes are determined from P-wave initial motion data. Although P-wave initial motion data for these offshore events are not ideally distributed on the focal sphere, we found that the upper-plane events that occur near the Japan Trench are characterized by normal faulting, whereas lower-plane events are characterized by thrust faulting. This focal mechanism distribution is the opposite to that of the double-planed deep seismic zone beneath the land area. The characteristics of these focal mechanisms for the shallow and deep doubled-planed seismic zones can be explained by a bending-unbending model of the subducting Pacific plate. Some of relocated earthquakes took place in the source area of the 1933 Mw8.4 Sanriku earthquake at depths of 10-23 km. The available focal mechanisms for these events are characterized by normal faulting. Given that the 1933 event was a large normal-fault event that occurred along a fault plane dipping landward, the earthquakes that currently occur just beneath or oceanwards of the Japan Trench are probably its aftershocks, suggesting that aftershock activity continues to the present day, 70 years after the main shock. ?? 2009 The Authors, Journal compilation ?? 2009 RAS.

Geophysical Journal International↗

USGS2011

Seismic swarm associated with the 2008 eruption of Kasatochi Volcano, Alaska: earthquake locations and source parameters

An energetic seismic swarm accompanied an eruption of Kasatochi Volcano in the central Aleutian volcanic arc in August of 2008. In retrospect, the first earthquakes in the swarm were detected about 1 month prior to the eruption onset. Activity in the swarm quickly intensified less than 48 h prior to the first large explosion and subsequently subsided with decline of eruptive activity. The largest earthquake measured as moment magnitude 5.8, and a dozen additional earthquakes were larger than magnitude 4. The swarm exhibited both tectonic and volcanic characteristics. Its shear failure earthquake features were b value = 0.9, most earthquakes with impulsive P and S arrivals and higher-frequency content, and earthquake faulting parameters consistent with regional tectonic stresses. Its volcanic or fluid-influenced seismicity features were volcanic tremor, large CLVD components in moment tensor solutions, and increasing magnitudes with time. Earthquake location tests suggest that the earthquakes occurred in a distributed volume elongated in the NS direction either directly under the volcano or within 5-10 km south of it. Following the M W 5.8 event, earthquakes occurred in a new crustal volume slightly east and north of the previous earthquakes. The central Aleutian Arc is a tectonically active region with seismicity occurring in the crusts of the Pacific and North American plates in addition to interplate events. We postulate that the Kasatochi seismic swarm was a manifestation of the complex interaction of tectonic and magmatic processes in the Earth's crust. Although magmatic intrusion triggered the earthquakes in the swarm, the earthquakes failed in context of the regional stress field.

Alaska↗

USGS2012

W phase source inversion for moderate to large earthquakes (1990-2010)

Rapid characterization of the earthquake source and of its effects is a growing field of interest. Until recently, it still took several hours to determine the first-order attributes of a great earthquake (e.g. M w ≥ 7.5), even in a well-instrumented region. The main limiting factors were data saturation, the interference of different phases and the time duration and spatial extent of the source rupture. To accelerate centroid moment tensor (CMT) determinations, we have developed a source inversion algorithm based on modelling of the W phase, a very long period phase (100–1000 s) arriving at the same time as the P wave. The purpose of this work is to finely tune and validate the algorithm for large-to-moderate-sized earthquakes using three components of W phase ground motion at teleseismic distances. To that end, the point source parameters of all M w ≥ 6.5 earthquakes that occurred between 1990 and 2010 (815 events) are determined using Federation of Digital Seismograph Networks, Global Seismographic Network broad-band stations and STS1 global virtual networks of the Incorporated Research Institutions for Seismology Data Management Center. For each event, a preliminary magnitude obtained from W phase amplitudes is used to estimate the initial moment rate function half duration and to define the corner frequencies of the passband filter that will be applied to the waveforms. Starting from these initial parameters, the seismic moment tensor is calculated using a preliminary location as a first approximation of the centroid. A full CMT inversion is then conducted for centroid timing and location determination. Comparisons with Harvard and Global CMT solutions highlight the robustness of W phase CMT solutions at teleseismic distances. The differences in M w rarely exceed 0.2 and the source mechanisms are very similar to one another. Difficulties arise when a target earthquake is shortly (e.g. within 10 hr) preceded by another large earthquake, which disturbs the waveforms of the target event. To deal with such difficult situations, we remove the perturbation caused by earlier disturbing events by subtracting the corresponding synthetics from the data. The CMT parameters for the disturbed event can then be retrieved using the residual seismograms. We also explore the feasibility of obtaining source parameters of smaller earthquakes in the range 6.0 ≤M w < 6.5. Results suggest that the W phase inversion can be implemented reliably for the majority of earthquakes of M w = 6 or larger.

Geophysical Journal International↗

USGS2009

Three-dimensional benchmark for variable-density flow and transport simulation: matching semi-analytic stability modes for steady unstable convection in an inclined porous box

This benchmark for three-dimensional (3D) numerical simulators of variable-density groundwater flow and solute or energy transport consists of matching simulation results with the semi-analytical solution for the transition from one steady-state convective mode to another in a porous box. Previous experimental and analytical studies of natural convective flow in an inclined porous layer have shown that there are a variety of convective modes possible depending on system parameters, geometry and inclination. In particular, there is a well-defined transition from the helicoidal mode consisting of downslope longitudinal rolls superimposed upon an upslope unicellular roll to a mode consisting of purely an upslope unicellular roll. Three-dimensional benchmarks for variable-density simulators are currently (2009) lacking and comparison of simulation results with this transition locus provides an unambiguous means to test the ability of such simulators to represent steady-state unstable 3D variable-density physics.

Hydrogeology Journal↗

USGS1996

Analytical solutions to non-Fickian subsurface dispersion in uniform groundwater flow

Analytical solutions are obtained by the Fourier transform technique for the one-, two-, and three-dimensional transport of a conservative solute injected instantaneously in a uniform groundwater flow. These solutions account for dispersive non-linearity caused by the heterogeneity of the hydraulic properties of aquifer systems and can be used as building blocks to construct solutions by convolution (principle of superposition) for source conditions other than slug injection. The dispersivity is assumed to vary parabolically with time and is thus constant for the entire system at any given time. Two approaches for estimating time-dependent dispersion parameters are developed for two-dimensional plumes. They both require minimal field tracer test data and, therefore, represent useful tools for assessing real-world aquifer contamination sites. The first approach requires mapped plume-area measurements at two specific times after the tracer injection. The second approach requires concentration-versus-time data from two sampling wells through which the plume passes. Detailed examples and comparisons with other procedures show that the methods presented herein are sufficiently accurate and easier to use than other available methods.

Journal of Hydrology↗