SEARCH · Geology Reports

Results for “Southeastern Geology”

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.

1,594 records · Page 89Linked to original sources

USGS2007

Organic-carbon sequestration in soil/sediment of the Mississippi River deltaic plain — Data; landscape distribution, storage, and inventory; accumulation rates; and recent loss, including a post-Katrina preliminary analysis

Soil/sediment of the Mississippi River deltaic plain (MRDP) in southeastern Louisiana is rich in organic carbon (OC). The MRDP contains about 2 percent of all OC in the surface meter of soil/sediment in the Mississippi River Basin (MRB). Environments within the MRDP differ in soil/sediment organic carbon (SOC) accumulation rate, storage, and inventory. The focus of this study was twofold: (1) develop a database for OC and bulk density for MRDP soil/sediment; and (2) estimate SOC storage, inventory, and accumulation rates for the dominant environments (brackish, intermediate, and fresh marsh; natural levee; distributary; backswamp; and swamp) in the MRDP. Comparative studies were conducted to determine which field and laboratory methods result in the most accurate and reproducible bulk-density values for each marsh environment. Sampling methods included push-core, vibracore, peat borer, and Hargis1 sampler. Bulk-density data for cores taken by the "short push-core method" proved to be more internally consistent than data for samples collected by other methods. Laboratory methods to estimate OC concentration and inorganic-constituent concentration included mass spectrometry, coulometry, and loss-on-ignition. For the sampled MRDP environments, these methods were comparable. SOC storage was calculated for each core with adequate OC and bulk-density data. SOC inventory was calculated using core-specific data from this study and available published and unpublished pedon data linked to SSURGO2 map units. Sample age was estimated using isotopic cesium ( 37 Cs), lead ( 210 Pb), and carbon ( 14 C), elemental Pb, palynomorphs, other stratigraphic markers, and written history. SOC accumulation rates were estimated for each core with adequate age data. Cesium-137 profiles for marsh soil/sediment are the least ambiguous. Levee and distributary 137 Cs profiles show the effects of intermittent allochthonous input and/or sediment resuspension. Cesium-137 and 210 Pb data gave the most consistent and interpretable information for age estimations of soil/sediment deposited during the 1900s. For several cores, isotopic 14 C and 137 Cs data allowed the 1963-64 nuclear weapons testing (NWT) peak-activity datum to be placed within a few-centimeter depth interval. In some cores, a too old 14 C age (when compared to 137 Cs and microstratigraphic-marker data) is the probable result of old carbon bound to clay minerals incorporated into the organic soil/sediment. Elemental Pb coupled with Pb source-function data allowed age estimation for soil/sediment that accumulated during the late 1920s through the 1980s. Exotic pollen (for example, Vigna unguiculata and Alternanthera philoxeroides ) and other microstratigraphic indicators (for example, carbon spherules) allowed age estimations for marsh soil/sediment deposited during the settlement of New Orleans (1717-20) through the early 1900s. For this study, MRDP distributary and swamp environments were each represented by only one core, backswamp environment by two cores, all other environments by three or more cores. MRDP core data for the surface meter soil/sediment indicate that (1) coastal marshes, abandoned distributaries, and swamps have regional SOC-storage values >16 kg m -2 ; (2) swamps and abandoned distributaries have the highest SOC storage values (swamp, 44.8 kg m -2 ; abandoned distributary, 50.9 kg m -2 ); (3) fresh-to-brackish marsh environments have the second highest site-specific SOC-storage values; and (4) site-specific marsh SOC storage values decrease as the salinity of the environment increases (fresh-marsh, 36.2 kg m -2 ; intermediate marsh, 26.2 kg m -2 ; brackish marsh, 21.5 kg m -2 ). This inverse relation between salinity and SOC storage is opposite the regional systematic increase in SOC storage with increasing salinity that is evident when SOC storage is mapped by linking pedon data to SSURGO map units (fresh marsh, 47 kg m -2 ; intermediate marsh, 67 kg m -2 ; brackish marsh, 75 kg m -2 ; and salt marsh, 80 kg m -2 ). MRDP core data for this study also indicate that levees and backswamp have regional SOC-storage values <16 kg m -2 . Group-mean SOC storage for cores from these environments are natural levee (17.0 kg m -2 ) and backswamp (14.1 kg m -2 ). An estimate for the SOC inventory in the surface meter of soil/sediment in the MRDP can be made using the SSURGO mapped portion of the coastal-marsh vegetative-type map (13,236 km 2 , land-only area) published by the Louisiana Department of Wildlife and Fisheries and U.S. Geological Survey (1997). This area has a SOC inventory (surface meter) of 677 Tg (slightly more than 2 percent of the 30,289 Tg SOC inventory for the MRB). The MRDP (6,180 km2, land-only area) has an estimated SOC inventory of 397 Tg. Most of the MRDP is located within the SSURGO mapped coastal marshlands. The entire MRDP, including water, has an area of about 10,800 km 2 . Using the ratio of total MRDP area to SSURGO mapped MRDP area as an adjustment, the MRDP SOC inventory is estimated at 694 Tg. This larger estimate of 694 Tg for the SOC inventory is probably more realistic, because it is reasonable to assume that the marsh sediments overlain by shallow water have comparable SOC storage to that of the adjacent land areas. MRDP core data for this study indicate that there is some variability in long-term SOC mass-accumulation rates for centuries and millennia and that this variability may indicate important geologic changes or changes in land use. However, the consistency of the range in rates of SOC accumulation through time suggests a remarkable degree of marsh sustainability throughout the Holocene, including the recent period of significant marsh modification/channelization for human use. One example of marsh sustainability is its present ability to function as a SOC sink even with Louisiana's large-scale coastal land loss during the last several decades. With coastal-marsh restoration efforts, this sink potential will increase. Looking to the future, a total of 1,101 g m -2 yr -1 SOC is projected to be lost from all of coastal Louisiana (U.S. Army Corps of Engineers, Louisiana Coastal Area (LCA) subprovinces 1-4; not just the MRDP) through coastal erosion from year 2000 to 2050. This translates to a projected SOC-loss rate of about 0.20 percent per year. The recent Hurricanes Katrina and Rita, which devastated the Louisiana coast during late August and late September 2005, transformed about 259 km 2 (100 mi 2 ) of marsh to open water (U.S. Geological Survey, 2005). To the extent that some or all of this land loss is permanent, this result equates to a SOC loss of about 15 Tg. This estimate is based on the year-2000 15,153-km 2 land area for the LCA study area that includes LCA subprovince 4. Using the year-2000 land area, the LCA study area had an estimated SOC inventory of 858 Tg. The estimated 15 Tg SOC loss attributable to Hurricanes Katrina and Rita is 1.7 percent of the year-2000 LCA inventory and 2.3 percent of the year-2000 MRDP inventory. If this SOC loss is included in the projection for the year 2050, then the MRDP would either remain a source with a net SOC loss of 3 Tg or become a weak sink with a net SOC gain of 4 Tg. These estimates are lower bounds for potential SOC flux because they are only for the surface meter of landmass.

Louisiana↗

USGS1941

Geologic structure and occurrence of gas in part of southwestern New York. Part 1, Structure and gas possibilities of the Oriskany sandstone in Steuben, Yates, and parts of the adjacent counties

The area covered by this report is in southwestern New York and includes a little more than 3,000 square miles in Steuben and Yates counties and parts of the six adjacent counties. This area has been mapped to determine the structural attitude of the exposed rocks, so as to aid those interested in prospecting for natural gas in the Oriskany sandstone of Lower Devonian age. Because of the gentle regional dip toward the southwest, the youngest beds are exposed in the southwest corner of the area, and progressively older beds crop out northeastward in successive bands that strike generally northwest. All the exposed rocks are of Upper Devonian age except those in a narrow belt at the extreme north edge of the area, where a small thickness of Middle Devonian rocks crops out. The maximum thickness of beds so exposed is nearly 4,000 feet, of which the lower part is predominantly soft dark shale and the upper part predominantly fine-grained sandstone and gray shale. All the beds are marine except a few tongues of continental deposits red shale and sandstone and gray mudstone in the youngest beds. All the beds thicken southeastward, so that there is a northwestward convergence between any two lithologic units in the series. More than 30 key horizons that are persistent and distinctive were mapped, and altitudes on these key horizons served as a basis for constructing the structure contour map. Many of the key horizons are formation or member boundaries, but others are the tops or bottoms of limestone or sandstone beds within formations. All the stratigraphic units mapped are purely lithologic. (See pi. 2.) The Tully limestone, which crops out along the northern border of the area, is an easily recognizable and therefore valuable key bed for subsurface correlations in this part of the State. Below the Tully limestone is a thick body of Middle Devonian shales of the Hamilton group which rests on another valuable key bed, the hard, cherty Onondaga limestone, also of Middle Devonian age. Below the Onondaga limestone is the Lower Devonian Oriskany sandstone, which is the gas-producing bed. Unlike the Onondaga, the Oriskany is locally thin or absent. The structure of most of the area is shown by contour lines at 25-foot intervals, but where key horizons are lacking the structure is indicated by dip symbols. Upon the regional south and southwest dip are superposed numerous gentle folds whose axes trend approximately northeastward in the greater part of the area but more nearly eastward in the eastern part. The folds generally tend to become narrower and steeper, and therefore more closely spaced, southwestward. Many of the anticlines fork southwestward, whereas the synclines. tend to fork northeastward. All the folds have a westward or southwestward plunge. Throughout the area the rocks are jointed in two dominant sets one that trends northwest and the other east or northeast. No evident relation .between these joints, which were measured only in the hard, relatively brittle beds, and the individual folds or domes was discernible. The faults are concentrated in the northeastern and southwestern parts of the area and trend either northeastward or northwestward. Some are nearly vertical normal faults; others are steep reverse faults. Subsurface data show that most of the faults increase in throw downward and also that many subsurface faults do not reach the surface. A group of faults in the northwestern part of the Greenwood quadrangle and the southwestern part of the Hornell quadrangle were active during Upper Devonian time, while the Gowanda shale and overlying beds were being deposited. At this stratigraphic horizon the beds in a zone a few hundred feet thick are highly deformed in a wide belt on both sides of the faults. Sandstone layers are thinned out into long stringers or swollen into thick masses and in places are bent acutely without fracture. Thin layers of shale, coquina, and sand have flowed together into intricately plicated zones that lack cleavage and joints. These features show that the sediments were deformed while wet and plastic and buried only a little way below the sea floor. The beds that were laid down over these disturbed zones were not involved in this deformation. Many of the sharper flexures and most of the faults are not evident in the beds several hundred feet stratigraphically higher. Accordingly, broad, gentle folds in these higher beds in parts of the area south and west of the northwest corner of the Greenwood quadrangle may conceal, at considerable depths below them, narrow folds separated by abrupt flexures or faults. Several of the larger streams and rivers occupy strike valleys, and their j courses swing to follow the changing strike of the rocks where they cross ( successive folds. But, with few exceptions, the small streams are not adjusted to the bedrock structure. Domes likely to serve as traps for natural gas are concentrated in the northeastern and southwestern parts of the area. The Wayne-Dundee gas field is in the northeastern part. All the other potentially valuable domes in this part of the area have been drilled and found valueless except one small structural feature in the southern part of the Ovid quadrangle, which, if the Oriskany is present, may trap a small quantity of gas. In the Greenwood quadrangle in the southwestern part of the area there is one gas field and four well-defined domes, all of which may be productive if the Oriskany sandstone is present. In the northwest corner of the quadrangle the dips indicate at least two domes that can be adequately defined and evaluated only by geophysical prospecting. The State Line gas field is in tbe Wellsville quadrangle. In the southeast corner of this quadrangle there are three other domes of comparable size that may also be productive if underlain by the Oriskany sandstone. At other places in the Wellsville quadrangle the dips suggest several anticlinal axes on which analogous productive domes maybe found. The structural features in this quadrangle, however, are defined by contours only in the southeastern part. In the Woodhull quadrangle a large dome east of Jasper may be productive, and the western top of the large Woodhull dome in the southwestern part of the quadrangle seems to warrant drilling, despite the absence of the Oriskany in a well on the eastern top. Two wells drilled in 1936 and 1937 a little northeast of a broad, nearly flat-topped dome in the Hornell quadrangle, a few miles east of Hornell,, struck small flows of gas, suggesting that wells drilled higher on this dome may be productive. In much of the southwestern part of the area seismograph surveys should be of great value in determining the structure at the Tully and Onondaga horizons. Without abundant subsurface control of this sort, the danger of drilling into subsurface faults can hardly be overemphasized. Three closed or nearly closed synclines in the Greenwood and Wellsville quadrangles appear to be favorable places to drill for oil in the shallow sands presumably parts of the Dunkirk sandstone.

New York↗

USGS2016

The Chahnaly low sulfidation epithermal gold deposit, western Makran volcanic arc, southeastern Iran

The Chahnaly low-sulfidation epithermal Au deposit and nearby Au prospects are located northwest of the intermittently active Bazman stratovolcano on the western end of the Makran volcanic arc, which formed as the result of subduction of the remnant Neo-Tethyan oceanic crust beneath the Lut block. The arc hosts the Siah Jangal epithermal and Kharestan porphyry prospects, near Taftan volcano, as well as the Saindak Cu-Au porphyry deposit and world-class Reko Diq Cu-Au porphyry deposit, near Koh-i-Sultan volcano to the east-northeast in Pakistan. The host rocks for the Chahnaly deposit include early Miocene andesite and andesitic volcaniclastic rocks that are intruded by younger dacitic domes. Unaltered late Miocene dacitic ignimbrites overlie these rocks. Laser ablation-inductively coupled plasma-mass spectrometry (LA-ICP-MS) U-Pb zircon geochronology data yield ages between 21.8 and 9.9 Ma for the acidic-intermediate regional volcanism. The most recent volcanic activity of the Bazman stratovolcano involved extrusion of an olivine basalt during Pliocene to Quaternary times. Interpretation of geochemical data indicate that the volcanic rocks are synsubduction and calc-alkaline to subalkaline. The lack of a significant negative Eu anomaly, a listric-shaped rare earth element pattern, and moderate La/Yb ratios of host suites indicate a high water content of the source magma. Gold and electrum are temporally and spatially related to a series of structurally controlled, 030°-trending, subvertical hydrothermal breccias with chalcedony-adularia that cut porphyritic andesite and andesitic volcaniclastic rocks. Gold is associated with pyrite, a siliceous matrix of hydrothermal breccia, and previously formed vein clasts, as well as with iron oxides and hydroxides in oxidized zones. Rare silver minerals include Ag-bearing electrum and naumannite, iodargyrite, an unnamed silver diiodide, and hessite. Hydrothermal alteration is generally well developed surrounding the ore-bearing hydrothermal breccia. The main types of alteration in the area include an inner ~0.5- to 20-m-thick gold-bearing hydrothermal breccia composed of quartz-chalcedonyadularia-illite-pyrite, a ~5- to 50-m-thick zone of quartz, chalcedony, pyrite, illitic phengite, phengite, illitic muscovite, illite, illitic paragonite, paragonite, muscovite, montmorillonite and, rarely, siderite, and a 30- to 70-m outer propylitic zone of Fe-Mg chlorite, calcite, ankerite, dolomite, epidote, palygorskite, and pyrite. The Chahnaly Au deposit formed during the early stages of magmatism. LA-ICP-MS zircon U-Pb geochronology of host andesite and 40Ar/39Ar dating of two samples of gold-associated adularia show that the ore-stage adularia (19.83 ± 0.10 and 19.2 ± 0.5 Ma) is younger, by as much as 1.5 million years, than the volcanic host rock (20.32 ± 0.4 Ma). Therefore, either hydrothermal activity continued well after volcanism or a second magmatic event rejuvenated hydrothermal activity. This second magmatic event may be related to eruption of porphyritic andesite at ~20.32 ± 0.40 Ma, which is within error of ~19.83 ± 0.10 Ma adularia. The new LA-ICP-MS zircon U-Pb host rock and vein adularia 40Ar/39Ar ages suggest that early Miocene magmatism and mineralization in the Bazman area is of a similar age to that of the Saindak porphyry and Tanjeel porphyry center of the giant Reko Diq deposit. This confirms the existence of early Miocene arc magmatism and mineralization along the Iranian part of the Makran volcanic arc. Ore, alteration mineralogy, and alteration patterns indicate that the Chahnaly deposit is a typical low-sulfidation epithermal Au deposit, located in a poorly explored part of the Makran volcanic arc in Iran.

Economic Geology↗

USGS2006

Hydrogeologic framework refinement, ground-water flow and storage, water-chemistry analyses, and water-budget components of the Yuma area, southwestern Arizona and southeastern California

The ground-water and surface-water system in the Yuma area in southwestern Arizona and southeastern California is managed intensely to meet water-delivery requirements of customers in the United States, to manage high ground-water levels in the valleys, and to maintain treaty-mandated water-quality and quantity requirements of Mexico. The following components in this report, which were identified to be useful in the development of a ground-water management model, are: (1) refinement of the hydrogeologic framework; (2) updated water-level maps, general ground-water flow patterns, and an estimate of the amount of ground water stored in the mound under Yuma Mesa; (3) review and documentation of the ground-water budget calculated by the Bureau of Reclamation, U.S. Department of the Interior (Reclamation); and (4) water-chemistry characterization to identify the spatial distribution of water quality, information on sources and ages of ground water, and information about the productive-interval depths of the aquifer. A refined three-dimensional digital hydrogeologic framework model includes the following hydrogeologic units from bottom to top: (1) the effective hydrologic basement of the basin aquifer, which includes the Pliocene Bouse Formation, Tertiary volcanic and sedimentary rocks, and pre-Tertiary metamorphic and plutonic rocks; (2) undifferentiated lower units to represent the Pliocene transition zone and wedge zone; (3) coarse-gravel unit; (4) lower, middle, and upper basin fill to represent the upper, fine-grained zone between the top of the coarse-gravel unit and the land surface; and (5) clay A and clay B. Data for the refined model includes digital elevation models, borehole lithology data, geophysical data, and structural data to represent the geometry of the hydrogeologic units. The top surface of the coarse-gravel unit, defined by using borehole and geophysical data, varies similarly to terraces resulting from the down cutting of the Colorado River. Clay A is nearly the same as the previous conceptual hydrogeologic model definition (Olmsted and others, 1973), except for a minor westward extension from the city of Yuma. Clay B is extended to the southerly international boundary and increased in areal extent by about two-thirds of the original extent (Olmsted and others, 1973). The other hydrogeologic units generally are the same as in the previous conceptual hydrogeologic model. Before development, the Colorado and Gila Rivers were the sources of nearly all the ground water in the Yuma area through direct infiltration of water from river channels and annual overbank flooding. After construction of upstream reservoirs and clearing and irrigation of the floodplains, the rivers now act as drains for the ground water. Ground-water levels in most of the Yuma area are higher now than they were in predevelopment time. A general gradient of ground-water flow toward the natural discharge area south of the Yuma area still exists, but many other changes in flow are evident. Ground water in Yuma Valley once flowed away from the Colorado River, but now has a component of flow towards the river and Mexicali Valley. A ground-water mound has formed under Yuma Mesa from long-term surface-water irrigation; about 600,000 to 800,000 acre-ft of water are stored in the mound. Ground-water withdrawals adjacent to the southerly international boundary have resulted in water-level declines in that area. The reviewed and documented water budget includes the following components: (1) recharge in irrigated areas, (2) evapotranspiration by irrigated crops and phreatophytes, (3) ground-water return flow to the Colorado River, and (4) ground-water withdrawals (including those in Mexicali Valley). Recharge components were calculated by subtracting the amount of water used by crops from the amount of water delivered. Evapotranspiration rates were calculated on the basis of established methods, thus were appropriate for input to the ground-water flow model developed by the Bureau of Reclamation (William Greer, hydrologist, Bureau of Reclamation, written commun., 2005). Evapotranspiration by crops and phreatophytes were calculated by using crop coefficient methods and meteorological data. Other methods of calculating evapotranspiration rates by using combinations of satellite imagery and ground-based data could be used for higher spatial and temporal resolution. Ground-water return flow during years of low flow on the Colorado River (1972–82, 1987–92, and 1994–96) averaged 79,000 acre-ft per year. Ground-water withdrawal data for 1970–99 were similar to other estimates made by the U.S. Geological Survey for the Yuma area. New water-chemistry data were collected in 12 wells and 8 canals/drains to characterize spatial patterns in chemical constituents, determine isotopic ages of water, infer possible sources of ground water, and locate the vertical intervals of the aquifer that contribute most water to wells. Depth-dependent samples were collected at one of the wells (YM-10). A large quantity of water-quality data were compiled from Bureau of Reclamation and U.S. Geological Survey records and merged into the U.S. Geological Survey National Water Information System database. New samples were analyzed for major ions, nutrients, stable isotopes of oxygen and hydrogen, tritium ( 3 H), and carbon-14 ( 14 C) (along with C 13 /C 12 ratios). Light values of oxygen-18 ( 18 O) and deuterium ( 2 H, D) in well 242-2 indicate recharge from the Colorado River. Heavy water samples from wells 242-22, CADC, and Mesa del Sol indicate local recharge sources. Tritium data indicate there is young water in wells in the valleys and near the edge of Yuma Mesa, while older water is found far from the Colorado River. 14 C data indicate that water from wells near the southerly international boundary is at least several thousand years old.

Arizona, California↗

USGS2016

U-Pb, Re-Os, and Ar/Ar geochronology of rare earth element (REE)-rich breccia pipes and associated host rocks from the Mesoproterozoic Pea Ridge Fe-REE-Au deposit, St. Francois Mountains, Missouri

Rare earth element (REE)-rich breccia pipes (600,000 t @ 12% rare earth oxides) are preserved along the margins of the 136-million metric ton (Mt) Pea Ridge magnetite-apatite deposit, within Mesoproterozoic (~1.47 Ga) volcanic-plutonic rocks of the St. Francois Mountains terrane in southeastern Missouri, United States. The breccia pipes cut the rhyolite-hosted magnetite deposit and contain clasts of nearly all local bedrock and mineralized lithologies. Grains of monazite and xenotime were extracted from breccia pipe samples for SHRIMP U-Pb geochronology; both minerals were also dated in one polished thin section. Monazite forms two morphologies: (1) matrix granular grains composed of numerous small (<50 μ m) crystallites intergrown with rare xenotime, thorite, apatite, and magnetite; and (2) coarse euhedral, glassy, bright-yellow grains similar to typical igneous or metamorphic monazite. Trace element abundances (including REE patterns) were determined on selected grains of monazite (both morphologies) and xenotime. Zircon grains from two samples of host rhyolite and two late felsic dikes collected underground at Pea Ridge were also dated. Additional geochronology done on breccia pipe minerals includes Re-Os on fine-grained molybdenite and 40 Ar/ 39 Ar on muscovite, biotite, and K-feldspar. Ages (±2 σ errors) obtained by SHRIMP U-Pb analysis are as follows: (1) zircon from the two host rhyolite samples have ages of 1473.6 ± 8.0 and 1472.7 ± 5.6 Ma; most zircon in late felsic dikes is interpreted as xenocrystic (age range ca. 1522–1455 Ma); a population of rare spongy zircon is likely of igneous origin and yields an age of 1441 ± 9 Ma; (2) pale-yellow granular monazite—1464.9 ± 3.3 Ma (no dated xenotime); (3) reddish matrix granular monazite—1462.0 ± 3.5 Ma and associated xenotime—1453 ± 11 Ma; (4) coarse glassy-yellow monazite—1464.8 ± 2.1, 1461.7 ± 3.7 Ma, with rims at 1447.2 ± 4.7 Ma; and (5) matrix monazite (in situ)—1464.1 ± 3.6 and 1454.6 ± 9.6 Ma, and matrix xenotime (in situ)—1468.0 ± 8.0 Ma. Two slightly older ages of cores are about 1478 Ma. The young age of rims on the coarse glassy monazite coincides with an Re-Os age of 1440.6 ± 9.2 Ma determined in this study for molybdenite intergrown with quartz and allanite, and with the age of monazite inclusions in apatite from the magnetite ore ( Neymark et al., 2016 ). A 40 Ar/ 39 Ar age of 1473 ± 1 Ma was obtained for muscovite from a breccia pipe sample. Geochronology and trace element geochemical data suggest that the granular matrix monazite and xenotime (in polygonal texture), and cores of coarse glassy monazite precipitated from hydrothermal fluids during breccia pipes formation at about 1465 Ma. The second episode of mineral growth at ca. 1443 Ma may be related to faulting and fluid flow that rebrecciated the pipes. The ca. 10-m.y. gap between the ages of host volcanic rocks and breccia pipe monazite and xenotime suggests that breccia pipe mineral formation cannot be related to the felsic magmatism represented by the rhyolitic volcanic rocks, and hence is linked to a different magmatic-hydrothermal system.

Missouri↗

USGS2000

SHRIMP U-Pb zircon ages for Big Creek gneiss, Wyoming and Boulder Creek batholith, Colorado: Implications for timing of Paleoproterozoic accretion of the northern Colorado province

Sensitive, high-resolution, ion microprobe (SHRIMP) U-Pb zircon ages from a sample of the high-grade, hornblende-feldspathic Big Creek gneiss of the southeastern Sierra Madre, along with samples of a quartz monzonitic phase of the Boulder Creek batholith, help define timing of three major Paleoproterozoic thermo-tectonic events within the northern Colorado province at approximately 1810, 1710, and 1610 Ma. Previous ages determined for these key rock units were problematic; they hindered regional tectonic interpretations of the Paleoproterozoic crustal accretion history of the Colorado province that extends from the Cheyenne belt of southern Wyoming to north-central New Mexico. The Colorado province has been popularly modelled as a series of accreted oceanic volcano-plutonic arc systems and associated sediments, although alternative interpretations suggest that the series represents continental-margin arc rocks. The Big Creek gneiss has been interpreted as a high-grade basement equivalent of the oldest arc volcanic rocks exposed within the Green Mountain arc terrane, but it also has been suspected of being either an older block of pre-arc basement or perhaps an allochthonous piece of crust from slightly older orogens to the east and north. Previous ID-TIMS work on mg-size zircon fractions yielded U-Pb concordia upper-intercept ages of 1618 ± 22 and 1684 ± 5 Ma as well as negative lower-intercept ages, indicating complex U-Pb isotopic systematics involving at least two ages of zircon growth overprinted by at least one episode of Pb-loss. Zircons from this gneiss were analyzed using the SHRIMP, and a total of 32 spot analyses on both centers and rims produced a range of different 207 Pb/ 206 Pb ages between ∼1840 and ∼1560 Ma. The weighted mean of the oldest 207 Pb/ 206 Pb ages is 1812 ± 12 Ma and is interpreted to estimate the age of the protolith that appears to be slightly older than lower-grade metabasalts and associated plutons at ∼1790–1775 Ma. This protolith age of 1812 Ma further implies that significantly older crust (> 1820 Ma; e.g., Penokean orogeny) is not found in the Green Mountain magmatic arc. The youngest 207 Pb/ 206 Pb ages of ∼1610 Ma are interpreted to represent a time of new zircon growth during highly localized high-grade metamorphism—an event that also produced local granitic magmatism at ∼1625 Ma. The Boulder Creek batholith had been dated previously using the ID-TIMS, U-Pb zircon technique that yielded ages at ∼1670 and ∼1714 Ma, a 45-m.y. discrepancy that left the true age of the batholith in doubt. Zircons from two samples, previously dated using the ID-TIMS method, were analyzed using SHRIMP, and yielded concordia upper-intercept ages of 1713 ± 10 and 1721 ± 15 Ma. These results, combined with two earlier U-Pb zircon determinations, help to establish the age of the Boulder Creek batholith at 1714.4 ± 4.6 Ma (weighted mean), an age more compatible with those for the other large, tonalitic to quartz monzonitic, syntectonic plutons within the northern Colorado province. The new Boulder Creek age helps to establish a discrete period of plutonism (∼1735–1705 Ma) that is syn- to post-tectonic with respect to major regional structures of deformation and metamorphism in the northern Colorado province. Assuming the multiple oceanic arc accretion model, the new age for the mid-crustal emplacement of this batholith into a deforming composite back-arc basin may date the closure of that basin during crustal shortening.

Colorado, Wyoming↗

USGS1963

Two pollen diagrams from southeastern Minnesota: Problems in the regional late-glacial and postglacial vegetational history

Kirchner Marsh and Lake Carlson are located 3 miles apart in Dakota County about 15 miles south of Minneapolis in the St. Croix moraine, which was formed by the Superior lobe during the Gary phase of the Wisconsin glaciation. During the Mankato phase that followed, the Des Moines lobe advanced to within a few miles of the sites. The region today is in a mixed-oak forest, with a maplebasswood forest 15 miles to the west and a re-entrant of the prairie on the sand plain south of the moraine. The general limit of coniferous trees is about 50 miles northeast of the sites, although outliers, especially of Pinus strobus, may be found along the Mississippi Valley a few miles to the east. One sediment core 12-13 m long from each site was analyzed for pollen content at 5-25-cm intervals. Diagrams based on percentage of total pollen (trees, shrubs, wind-pollinated herbs) show essentially identical sequences at the two sites, starting with the late-glacial phase of ice retreat. The diagrams have been subdivided into pollen zones according to the A-B-C sequence introduced by Deevey for New England. The late-glacial pollen record starts at Kirchner Marsh with a short Picea-Cyperaceae-Gramineae phase (Zone K), believed to represent a spruce parkland. Its C-14 date of 13,270 BP and the stratigraphy indicate a pre- Two Creeks and post- Gary correlation. Apparently the Kirchner site did not become established as a lake until this time owing to persistence of dead ice in the moraine. The absence of pollen of specific tundra indicators and the presence of pollen of such thermophilous plants as Fraxinus, Quercus, Corylus, Ambrosia, Humulus, and Typha latifolia imply that the climate was cool rather than cold. Zone A-a, which follows, correlates with the Two Creeks interstade. It is marked by the dominance of Picea, with appreciable percentages of Fraxinus and Ambrosia and with minor amounts of other thermophilous plants and the normal boreal associates of spruce like Betula, Larix, and Salix. Zone A-b, starting 12,050 C-14 years ago, correlates with the Valders ice advance. It is represented at both Kirchner and Carlson and shows the withdrawal of Fraxinus and Ambrosia and the slight rise of Artemisia. Except for the absence of pine in the late-glacial assemblage the vegetation implied by these three zones seems to have its closest modern counterpart in the southern fringe of the Boreal Forest of the Riding Mountain region of southwest Manitoba. It is concluded that pine did not migrate southward with the spruce during the Wisconsin glaciation, at least in the western Great Lakes region, and was thus eliminated from this region. During the lateglacial phases of ice retreat, herbs and spruce pioneered on the deglaciated terrain; pine did not follow until the destruction of the spruce forest at the end of the late-glacial phase. Zone B introduces postglacial time. It represents the time of rapid Vegetational succession following the deterioration of the spruce forest. Simultaneous maxima of Betula, Alnus, Fraxinus, and Abies occurred 10,230 years ago at Kirchner Marsh. These were followed rapidly by a Pinus maximum and then a rise of Ulmus, Quercus, and other deciduous types, dated as 9300 years ago at the correlative site of Madeha. This succession may represent differential rates of migration from refuges south and east of Minnesota . Deciduous trees dominate the C Zones. Zone C-a shows Ulmus and Ostrya /Carpinus followed by Quercus; it probably represents principally a mesic maple-basswood forest changing to oak. Zone C-b represents the advance of prairie into the region at the expense of the oak woodland or savanna. The large and abrupt fluctuations in the curves for Ambrosia-type and Chenopodiineae, especially at the Carlson site, may record encroachment of annual weeds onto intermittently dried lake bottoms. C-14 dates place Zone C-b between 7100 and 5100 years ago. In Zone C-c the Quercus again dominates until the abrupt increase in Ambrosiatype and Chenopodiineae that marks the time of forest clearance and land settlement 50-75 years ago.

Minnesota↗

USGS2023

Integrated science strategy for assessing and monitoring water availability and migratory birds for terminal lakes across the Great Basin, United States

Executive Summary In 2022, the U.S. Geological Survey (USGS) established the Saline Lake Ecosystems Integrated Water Availability Assessment (IWAAs) to monitor and assess the hydrology of terminal lakes in the Great Basin and the migratory birds and other wildlife dependent on those habitats. Scientists from across the USGS (with specialties in water quantity, water quality, limnology, avian biology, data science, landscape ecology, and science communication) formed the Saline Lake Ecosystems IWAAs Team. The team has developed this regional strategic science plan to guide data collection and assessment activities at terminal lakes in the Great Basin. The U.S. Congress requested the USGS to establish the Saline Lake Ecosystems IWAAs in response to historically low water levels at terminal lakes and associated wetlands across the Great Basin. Not all Great Basin terminal lakes have high salinity; however, all terminal lakes occur in endorheic, closed, basins with no surface-water outflow. Low lake levels across the Great Basin are the result of increased water use for agriculture and municipalities, drought conditions, and a warming climate. Great Basin terminal lake water extents have decreased by as much as 90 percent over the last 150 years, and terminal lake wetlands have decreased in area by as much as 47 percent since 1984. Lake elevations and wetland areas are primarily supported by freshwater inputs from snowmelt feeding upgradient rivers, streams, and springs. These freshwater inputs have been severely reduced because of continued and increased surface-water diversions and surface-water capture through groundwater pumping for agriculture, mining, and public supply as well as unprecedented drought conditions and warming temperatures related to climate change. Water quality, specifically salinity, is highly variable for terminal lakes of the Great Basin, and this variability is a result of the balance between freshwater inflow and evaporation. Variability of salinity at each of the terminal lakes can be affected by lake morphology, hydrogeologic features of the basin, annual variability in weather patterns, and changes in upgradient water use. Hypersaline terminal lakes provide abundant food resources such as brine shrimp and brine flies that support nesting and migrating birds. The density and composition of invertebrates are closely tied to lake salinity. Increased salinity can exceed the tolerance of invertebrates, severely limiting their biomass. In contrast, decreased salinity can lead to altered invertebrate community composition, reducing the abundance of optimal avian prey resources. Great Basin terminal lake ecosystems, including open-water and adjacent aquatic and terrestrial environments, provide resources necessary to sustain many animal populations throughout the year. Although a variety of taxa use terminal lakes, these ecosystems are of acute importance for the millions of migratory waterbirds (for example, shorebirds, wading birds, and waterfowl) dependent on the network of terminal lakes and their associated wetlands. Migratory birds transiting the Pacific and Central Flyways use Great Basin terminal lake ecosystems throughout the year to feed, nest, and transit between wintering and breeding ranges. As such, successful conservation of birds and their habitats requires coordinated management of water and habitats across the Great Basin network of terminal lakes and wetlands. The linkages between water availability and ecosystem vulnerability of terminal lakes in the Great Basin are not well understood. The vulnerability of terminal lakes is related to the factors driving change and adaptive capacity of the lake ecosystem. Saline lake ecosystems are vulnerable when changes in water quantity affect ecosystem function. Water quantity affects salinity, which affects food webs and habitat; these linkages can be investigated with water-quality and food web monitoring. Water quantity also affects inundated habitat, which can be quantified through remote sensing. It is necessary to quantify hydroclimatic and water use controls on water availability to terminal lakes to assess the response of the ecosystems. Remotely sensed data can provide a broad-scale and long-term synoptic view of terminal lake hydrologic characteristics, but ground observations are required to interpret changes in water quality and ecological functions. Some terminal lake basins have ongoing monitoring and modeling efforts within the Great Basin (for example, Great Salt Lake, Carson River Basin), yet most monitoring locations are hydrologically upgradient and too far away from lake inflows to provide an accurate assessment of hydrological trends for the lake ecosystems. Other terminal lakes have no long-term hydrological monitoring in their respective watersheds (for example, Lake Abert). Ecological data collection in the Great Basin is also insufficient to understand how many birds exist on the landscape, how birds use the mosaic of terminal-lake habitats as an interconnected system, and how Great Basin terminal lakes are linked to the larger continental system of the Pacific and Central Flyways. Across agencies and organizations, tracking bird movement, abundance, and diversity is inconsistent, with some lakes having once- or twice-a-year bird survey efforts and a few locations having more intensive ecological data-gathering efforts (for example, Great Salt Lake, Lake Abert). Bridging hydrological and ecological information gaps will improve understanding of the trends in water supply and water quality, habitat availability and usage, and impacts on vulnerable waterbird species, all of which would be used by managers in coordinated conservation of this unique network of terminal-lake habitats. The terminal lakes of the Great Basin are part of the Basin and Range physiographic province that extends from the Colorado Plateau on the east to the Sierra Nevada on the west, and from the Snake River Plain on the north to the Garlock fault and the Mojave block on the south. The Great Basin is larger than 650,000 square kilometers and encompasses most of the State of Nevada but also extends to western Utah, eastern California, southeastern Idaho, southwestern Wyoming, and southeastern Oregon. The climate is arid to semiarid with a hydrologic regime that is snowmelt dominated, providing as much as 75 percent of total annual runoff for the region. Terminal lakes of the Great Basin occupy the lowest areas of closed (endorheic) drainage basins, such that lake levels and water quality respond rapidly to surface-water inflow. Terminal lakes provide local and regional economic value to the States in the Great Basin, including mineral extraction, aquaculture, public works, and recreational uses. As an example, assessments of Great Salt Lake’s ecological health and economic impact find hemispheric importance for the former and regional importance for the latter. Great Salt Lake creates about 7,000 jobs and $2 billion of economic output per year, most of which would be lost with further declines in lake level. The objectives of this Science Strategy are threefold: (1) to identify how changing water availability affects the quality, diversity, and abundance of habitats supporting continental waterbird populations; (2) to highlight the scientific monitoring and assessment needs of Great Basin terminal lakes; and (3) to support coordinated management and conservation actions to benefit those ecosystems, migratory birds, and other wildlife. There are long-term hydrological, ecological, and societal challenges associated with terminal lakes ecosystems in the Great Basin. This Science Strategy benefits partners by providing a conceptual model, nested at different spatial extents, that identifies key scientific information needs to inform coordinated implementation of management and conservation plans within and among hydrologic basins to address these complex challenges.

California, Idaho, Nevada, Oregon, Utah, Wyoming↗

USGS2013

The use of process models to inform and improve statistical models of nitrate occurrence, Great Miami River Basin, southwestern Ohio

Statistical models of nitrate occurrence in the glacial aquifer system of the northern United States, developed by the U.S. Geological Survey, use observed relations between nitrate concentrations and sets of explanatory variables—representing well-construction, environmental, and source characteristics— to predict the probability that nitrate, as nitrogen, will exceed a threshold concentration. However, the models do not explicitly account for the processes that control the transport of nitrogen from surface sources to a pumped well and use area-weighted mean spatial variables computed from within a circular buffer around the well as a simplified source-area conceptualization. The use of models that explicitly represent physical-transport processes can inform and, potentially, improve these statistical models. Specifically, groundwater-flow models simulate advective transport—predominant in many surficial aquifers— and can contribute to the refinement of the statistical models by (1) providing for improved, physically based representations of a source area to a well, and (2) allowing for more detailed estimates of environmental variables. A source area to a well, known as a contributing recharge area, represents the area at the water table that contributes recharge to a pumped well; a well pumped at a volumetric rate equal to the amount of recharge through a circular buffer will result in a contributing recharge area that is the same size as the buffer but has a shape that is a function of the hydrologic setting. These volume-equivalent contributing recharge areas will approximate circular buffers in areas of relatively flat hydraulic gradients, such as near groundwater divides, but in areas with steep hydraulic gradients will be elongated in the upgradient direction and agree less with the corresponding circular buffers. The degree to which process-model-estimated contributing recharge areas, which simulate advective transport and therefore account for local hydrologic settings, would inform and improve the development of statistical models can be implicitly estimated by evaluating the differences between explanatory variables estimated from the contributing recharge areas and the circular buffers used to develop existing statistical models. The larger the difference in estimated variables, the more likely that statistical models would be changed, and presumably improved, if explanatory variables estimated from contributing recharge areas were used in model development. Comparing model predictions from the two sets of estimated variables would further quantify—albeit implicitly—how an improved, physically based estimate of explanatory variables would be reflected in model predictions. Differences between the two sets of estimated explanatory variables and resultant model predictions vary spatially; greater differences are associated with areas of steep hydraulic gradients. A direct comparison, however, would require the development of a separate set of statistical models using explanatory variables from contributing recharge areas. Area-weighted means of three environmental variables—silt content, alfisol content, and depth to water from the U.S. Department of Agriculture State Soil Geographic (STATSGO) data—and one nitrogen-source variable (fertilizer-application rate from county data mapped to Enhanced National Land Cover Data 1992 (NLCDe 92) agricultural land use) can vary substantially between circular buffers and volume-equivalent contributing recharge areas and among contributing recharge areas for different sets of well variables. The differences in estimated explanatory variables are a function of the same factors affecting the contributing recharge areas as well as the spatial resolution and local distribution of the underlying spatial data. As a result, differences in estimated variables between circular buffers and contributing recharge areas are complex and site specific as evidenced by differences in estimated variables for circular buffers and contributing recharge areas of existing public-supply and network wells in the Great Miami River Basin. Large differences in areaweighted mean environmental variables are observed at the basin scale, determined by using the network of uniformly spaced hypothetical wells; the differences have a spatial pattern that generally is similar to spatial patterns in the underlying STATSGO data. Generally, the largest differences were observed for area-weighted nitrogen-application rate from county and national land-use data; the basin-scale differences ranged from -1,600 (indicating a larger value from within the volume-equivalent contributing recharge area) to 1,900 kilograms per year (kg/yr); the range in the underlying spatial data was from 0 to 2,200 kg/yr. Silt content, alfisol content, and nitrogen-application rate are defined by the underlying spatial data and are external to the groundwater system; however, depth to water is an environmental variable that can be estimated in more detail and, presumably, in a more physically based manner using a groundwater-flow model than using the spatial data. Model-calculated depths to water within circular buffers in the Great Miami River Basin differed substantially from values derived from the spatial data and had a much larger range. Differences in estimates of area-weighted spatial variables result in corresponding differences in predictions of nitrate occurrence in the aquifer. In addition to the factors affecting contributing recharge areas and estimated explanatory variables, differences in predictions also are a function of the specific set of explanatory variables used and the fitted slope coefficients in a given model. For models that predicted the probability of exceeding 1 and 4 milligrams per liter as nitrogen (mg/L as N), predicted probabilities using variables estimated from circular buffers and contributing recharge areas generally were correlated but differed significantly at the local and basin scale. The scale and distribution of prediction differences can be explained by the underlying differences in the estimated variables and the relative weight of the variables in the statistical models. Differences in predictions of exceeding 1 mg/L as N, which only includes environmental variables, generally correlated with the underlying differences in STATSGO data, whereas differences in exceeding 4 mg/L as N were more spatially extensive because that model included environmental and nitrogen-source variables. Using depths to water from within circular buffers derived from the spatial data and depths to water within the circular buffers calculated from the groundwater-flow model, restricted to the same range, resulted in large differences in predicted probabilities. The differences in estimated explanatory variables between contributing recharge areas and circular buffers indicate incorporation of physically based contributing recharge area likely would result in a different set of explanatory variables and an improved set of statistical models. The use of a groundwater-flow model to improve representations of source areas or to provide more-detailed estimates of specific explanatory variables includes a number of limitations and technical considerations. An assumption in these analyses is that (1) there is a state of mass balance between recharge and pumping, and (2) transport to a pumped well is under a steady state flow field. Comparison of volumeequivalent contributing recharge areas under steady-state and transient transport conditions at a location in the southeastern part of the basin shows the steady-state contributing recharge area is a reasonable approximation of the transient contributing recharge area after between 10 and 20 years of pumping. The first assumption is a more important consideration for this analysis. A gradient effect refers to a condition where simulated pumping from a well is less than recharge through the corresponding contributing recharge area. This generally takes place in areas with steep hydraulic gradients, such as near discharge locations, and can be mitigated using a finer model discretization. A boundary effect refers to a condition where recharge through the contributing recharge area is less than pumping. This indicates other sources of water to the simulated well and could reflect a real hydrologic process. In the Great Miami River Basin, large gradient and boundary effects—defined as the balance between pumping and recharge being less than half—occurred in 5 and 14 percent of the basin, respectively. The agreement between circular buffers and volume-equivalent contributing recharge areas, differences in estimated variables, and the effect on statisticalmodel predictions between the population of wells with a balance between pumping and recharge within 10 percent and the population of all wells were similar. This indicated process-model limitations did not affect the overall findings in the Great Miami River Basin; however, this would be model specific, and prudent use of a process model needs to entail a limitations analysis and, if necessary, alterations to the model.

Ohio↗

USGS2017

Geophysics- and geochemistry-based assessment of the geochemical characteristics and groundwater-flow system of the U.S. part of the Mesilla Basin/Conejos-Médanos aquifer system in Doña Ana County, New Mexico, and El Paso County, Texas, 2010–12

One of the largest rechargeable groundwater systems by total available volume in the Rio Grande/Río Bravo Basin (hereinafter referred to as the “Rio Grande”) region of the United States and Mexico, the Mesilla Basin/Conejos-Médanos aquifer system, supplies water for irrigation as well as for cities of El Paso, Texas; Las Cruces, New Mexico; and Ciudad Juárez, Chihuahua, Mexico. The U.S. Geological Survey in cooperation with the Bureau of Reclamation assessed the groundwater resources in the Mesilla Basin and surrounding areas in Doña Ana County, N. Mex., and El Paso County, Tex., by using a combination of geophysical and geochemical methods. The study area consists of approximately 1,400 square miles in Doña Ana County, N. Mex., and 100 square miles in El Paso County, Tex. The Mesilla Basin composes most of the study area and can be divided into three parts: the Mesilla Valley, the West Mesa, and the East Bench. The Mesilla Valley is the part of the Mesilla Basin that was incised by the Rio Grande between Selden Canyon to the north and by a narrow valley (about 4 miles wide) to the southeast near El Paso, Tex., named the Paso del Norte, which is sometimes referred to in the literature as the “El Paso Narrows.” Previously published geophysical data for the study area were compiled and these data were augmented by collecting additional geophysical and geochemical data. Geophysical resistivity measurements from previously published helicopter frequency domain electromagnetic data, previously published direct-current resistivity soundings, and newly collected (2012) time-domain electromagnetic soundings were used in the study to detect spatial changes in the electrical properties of the subsurface, which reflect changes that occur within the hydrogeology. The geochemistry of the groundwater system was evaluated by analyzing groundwater samples collected in November 2010 for physicochemical properties, major ions, trace elements, nutrients, pesticides (reported but not used in the assessment), and environmental tracers. The data obtained from these samples (with the exception of the pesticide data) were used to gain insights into processes controlling the groundwater movement through the groundwater system in the study area. Results from the geophysical and geochemical assessments facilitated the interpretation of the geochemical characteristics of the groundwater sources and geochemical groups within the groundwater system. The groundwater-flow system in the study area consists primarily of the Mesilla Basin aquifer system, which can be divided into four hydrogeologic units by using an informal classification scheme based on basin-fill stratigraphy and sedimentology with an emphasis on aquifer characteristics. The four hydrogeologic units are (1) the Rio Grande alluvium, which is the shallow aquifer of the Mesilla Basin within the confines of the Mesilla Valley, and the three hydrogeologic units that compose the Santa Fe Group: (2) the lower part of the Santa Fe Group, which is the least productive zone, (3) the middle part of the Santa Fe Group, which is the primary water-bearing hydrogeologic unit in the basin and is generally saturated, and (4) the upper part of the Santa Fe Group, which is the most productive water-bearing unit within the Santa Fe Group but is only partially saturated in the north and largely unsaturated in the south and western parts of the Mesilla Basin. The helicopter frequency domain electromagnetic survey results indicated that approximately half of the resistivity values were less than 10 ohm-meters at depths of 50 and 100 feet with a transition where the resistivity values changed from relatively high values (greater than 20 ohm-meters) to relatively low resistivity values (less than 10 ohm-meters) near Vado, New Mexico. Slightly more than 25 percent of the gridded resistivity values from the three-dimensional grid of the combined inverse modeling results of the direct-current resistivity and time-domain electromagnetic soundings were equal to or less than 10 ohm-meters with large regions of low resistivity becoming apparent in the southernmost part of the study area near the Paso Del Norte where these low resistivity features are spatially the widest at or below the top of the bedrock. These low resistivity values might represent clayey deposits, sediments composed largely of sand and gravel saturated with saline water, or both. Historical dissolved-solids-concentration data within the surface geophysical subset area of the study area were compiled and compared to the inverse modeling results of the combined direct-current resistivity and time-domain soundings; this comparison was done to strengthen the interpretation made from the combined inverse modeling results that the low resistivity features were representative of sand and gravel deposits saturated with saline water and not clayey deposits. Water-level altitudes within the Rio Grande alluvium generally decreased from north to south, with a west to east decrease in water-level altitudes near Las Cruces, New Mexico, as a result of groundwater pumping. Groundwater flow within the Santa Fe Group is more complex than the groundwater flow within the Rio Grande alluvium because of the larger lateral and vertical extent of the Santa Fe Group compared to the Rio Grande alluvium. Groundwater from the Organ Mountains flows directly south towards the Paso del Norte. Groundwater from the Robledo Mountains, the Rough and Ready Hills, and the Sleeping Lady Hills generally flows to the southeast. Groundwater flowing near the north end of the midbasin uplift generally continues east towards the Rio Grande and then flows south on the east side of the midbasin uplift. Groundwater flowing near the west side of the midbasin uplift generally continues south parallel to the faults that make up the midbasin uplift and then flows east towards the Paso del Norte when it reaches the south end of the midbasin uplift. Groundwater from the Aden Hills and the East and West Potrillo Mountains flows to the south end of the midbasin uplift and then continues east towards the Paso del Norte. Throughout most of the Mesilla Valley, the vertical hydraulic gradient was downward because the water-level altitude in the Rio Grande alluvium was higher than it was in the Santa Fe Group, but in some areas (typically in the middle and southern parts of the Mesilla Valley), the vertical hydraulic gradient was substantially reduced or even reversed to an upward hydraulic gradient. The geochemistry data indicate that there was a complex system of multiple geochemical endmembers and mixing between these endmembers with recharge to the Rio Grande alluvium and Santa Fe Group composed mostly of seepage from the Rio Grande, inflows from deeper or neighboring water systems, and mountain-front recharge. Five distinct geochemical groups were identified in the Mesilla Basin study area: (1) ancestral Rio Grande (pre-Pleistocene) geochemical group, (2) modern Rio Grande (Pleistocene to present) geochemical group, (3) mountain-front geochemical group, (4) deep groundwater upwelling geochemical group, and (5) unknown freshwater geochemical group. The ancestral Rio Grande groundwater was water that recharged into the system as seepage losses from the ancestral Rio Grande; this groundwater generally flows from north to south-southeast towards the Paso del Norte. Groundwater on the west side of the midbasin uplift generally flows south until it reaches the southern part of the study area; from the southern part of the study area, the groundwater flows east towards the Paso del Norte. Groundwater on the east side of the midbasin uplift flows south-southeast towards the Paso del Norte where it mixes with groundwater from the modern Rio Grande, uplifted areas in the west, and the deep saline source. The water type of the modern Rio Grande geochemical group ranged from calcium-sulfate water type in the northern part of the study area to sodium-chloride-sulfate water type in the southern part of the study area; from north to south there was a substantial increase in specific conductance, strontium-87/strontium-86 ratio, potassium, and the trace metals of iron and lithium, changing the water chemistry such that it became similar to the water chemistry of the deep groundwater upwelling geochemical group. From age-dating results, water in the modern Rio Grande geochemical group was recharged to the Rio Grande alluvium within the past 10 years. The mountain-front geochemical group was generally old water (apparent age was greater than 10,000 carbon-14 years before present) that was somewhat mineralized and has relatively high concentrations of fluoride and silica, which might indicate longer exposure to volcanic and siliciclastic rocks or aluminosilicate minerals. There were five different locations of recharge determined from the groundwater geochemistry within the mountain-front geochemical group, all having a slightly different geochemical signature: (1) the Rough and Ready Hills, Robledo Mountains, and the Sleeping Lady Hills, (2) the Doña Ana Mountains, (3) the Aden Hills and West Potrillo Mountains, (4) the East Potrillo Mountains, and (5) the Sierra Juárez in Mexico. The groundwater from the Rough and Ready Hills, Robledo Mountains, the Sleeping Lady Hills, and the Doña Ana Mountains generally flows toward the Rio Grande and eventually mixes together and with the modern Rio Grande groundwater. The groundwater originating from the Aden Hills and East and West Potrillo Mountains generally flows east to southeast at a slow rate and eventually mixes and continues east, where it mixes with groundwater from the ancestral Rio Grande geochemical group and with the groundwater from the Sierra Juárez. The groundwater from the Sierra Juárez flows north and then east towards the Paso del Norte where it mixes with groundwater from the uplifted areas in the west, ancestral and modern Rio Grande groundwater, and the upwelling groundwater from a deep saline source. The deep groundwater upwelling geochemical group had the highest concentrations of bicarbonate, potassium, silica, aluminum, iron, and lithium within the study area, indicating that it had been in contact with carbonate and siliciclastic rocks for a much longer period of time and at higher temperatures compared to the other geochemical groups, and was most likely ancient marine groundwater originating from the Paleozoic and Cretaceous carbonate rocks which was upwelling into the Mesilla Basin aquifer system in the southeastern part of the study area through the extensive fault systems. Direct-current resistivity and time-domain electromagnetic soundings support the interpretation of ancient marine groundwater upwelling into the Mesilla Basin aquifer system, as do the analytical results from wells, and the helicopter frequency domain electromagnetic data collected along the Rio Grande. The hydrogen-2/hydrogen-1 ratio and oxygen-18/oxygen-16 ratio isotopic results for samples in the unknown freshwater geochemical group did not plot on the Rio Grande evaporation line, indicating this group did not have a Rio Grande signature (that is, there was no isotopic evidence of a component of Rio Grande water) and it also had the lowest mineralized content of any geochemical group in the study area.

New Mexico, Texas↗