Isotopes as Tracers of Water Origin in and near a Regional Carbonate Aquifer: the Southern Sacramento Mountains, New Mexico

High-elevation groundwater sampled in 2003 in the Sacramento Mountains defines a line resembling an evaporation trend in δD-δ 18 O space. The trend results from recharge of winter precipitation into fractured limestone, with evaporation prior to recharge in broad mountain valleys. The same trend occurs in basin groundwater east and west of the range, indicating the high Sacramento Mountains as the principal regional water source, either direct from the limestone aquifers or from mountain-derived surface water. Tritium and carbon-14 indicate bulk residence times of a few decades in the high Sacramento Mountains and at Alamogordo, and of thousands of years south of Alamogordo and in the artesian aquifer near Artesia. Stable O, H isotope data fail to demonstrate the presence of Sacramento Mountains water in a saline aquifer of the Hueco Bolson (Texas).


Introduction
In the Basin and Range province of the southwest USA, stable isotope studies have proved useful in distinguishing sources of recharge where altitude effects are large, e.g., Tucson Basin [1], or where OPEN ACCESS isotope effects due to latitude/altitude and evaporation generate river water that is distinctive beside native basin groundwater [2,3].
The Sacramento Mountains of south-central New Mexico (Figure 1) include a broad area of forested, well-watered terrain, the source of perennial streams flowing east toward the Roswell Basin and the Pecos River and of intermittent streams flowing west into the Tularosa Valley.Pre-development and recent water level data indicate groundwater movement from Tularosa Valley to the Hueco Bolson [4][5][6].Mayer and Sharp [7] suggested regional flow of groundwater from the Sacramento Mountains to the Texas-New Mexico border, through karst aquifers.Groundwater samples for this study were collected in the Sacramento Mountains and the flanking basins from 2003 to 2008.The aim of the study was to use environmental isotopes to determine the relationship between water from the Sacramento Mountains and the adjacent basin aquifers.The relationship between groundwater in the high mountains and that in the Tularosa Valley, the deep alluvial basin to the west, is the first topic to be addressed.The relationship between groundwater in the high mountains and that in the hard-rock Roswell artesian basin to the east is the second topic to be addressed.In both cases, we also attempt to constrain groundwater residence times, and to determine the seasonality of recharge.Finally, we discuss whether the isotope signature of Sacramento Mountains groundwater can be recognized as far south as the Hueco Bolson in Texas (Figure 1).

Background
Figure 1 shows the location of the study areas.Area 1, encompassing the Sacramento Mountains near Cloudcroft and Weed, and the freshwater lens on the western flank of the mountains, encompasses sites sampled for the first topic described above.Area 2 stretches from the eastern flank of the mountains to the Pecos River at Artesia, and encompasses sites samples for the second topic.

Topography, Climate and Vegetation
The Sacramento Mountains rise to 2500-2800 m above sea level (m.a.s.l.) in the study area.A steep western escarpment with deep canyons abuts the Tularosa Valley, a typical fault-bounded basin of the Basin-and-Range province.Tularosa Valley continues southward into Texas where the valley is named the Hueco Bolson; Neogene alluvium fills the entire extent of the combined basin to depths of 600 to 3000 m in the basin center [6].The eastern flank of the Sacramento Mountains approximates a dip-slope, and descends gradually toward the Pecos River.No deep alluvial basin is present on the east side of the range.The climate in the basins is semi-arid; average annual precipitation is 335 mm at Alamogordo and 340 mm at Artesia.In the high mountains, precipitation is higher, e.g., 715 mm at Cloudcroft near the range crest [8].There are two wet seasons, a weak summer monsoon (June to October) providing 65%-70% of the precipitation, and a winter season of rain and snow from frontal weather systems [8].The amounts of both winter and summer precipitation vary greatly from year to year (Figure 15 of [9]).Vegetation consists of coniferous forest interspersed with grassy valleys above 2300 m.a.s.l..At lower elevation, scrubby oak forest and desert scrub predominate, except along perennial streams where riparian forest is present.Much of the study area is dry ranch land on which groundwater pumping is essential to the survival of cattle herds.Large-scale irrigated agriculture, using quarter-section and larger center-pivot and side-roll equipment, is practiced on the Pecos River flood plain near Artesia.Scattered irrigated plots are present near Tularosa.

Geology
The Sacramento Mountains constitute a tilted horst with range-bounding faults on the western side (Figure 2), and consist of Paleozoic marine sedimentary rocks, mostly Permian-Mississippian limestone and evaporite, overlying concealed Precambrian basement [9,10].The surface east of the range crest approximates a dip slope, with a discontinuous veneer of lower members of the San Andres formation overlying dolomite and anhydrite of the Yeso formation, both overlain to the east by the Queen-Grayburg anhydrite and limestone and dolomite of the San Andres formation.Thin (40 m) Quaternary alluvium overlies Paleozoic strata on the Pecos River flood plain.West of the range crest, the entire Paleozoic section of the region, mainly carbonate strata, is exposed.Neogene alluvium fills the Tularosa Valley to depths of 230 to 300 m at Alamogordo.

Geohydrology
The carbonate strata constitute a regional aquifer system conveying water from the mountains to the basins, eastward from the range crest through the Yeso formation, and westward through the highly fractured Paleozoic carbonate section.A map of the potentiometric surface east of the range crest is available in ([9], Figure 18), and indicates general eastward groundwater flow.West of the range crest, groundwater levels are less precisely known within area 1, but they decline steeply towards the west and southwest elsewhere on the escarpment [9].In the high mountains, the geohydrology is complex and governed by the detailed lithology of the Yeso formation.The following geohydrologic features are present ( [9], Figure 25): a regional aquifer, locally confined beneath impermeable interbeds of the Yeso formation, and probably continuous with regional aquifers east of the range; multiple perched aquifers overlying the regional aquifer, some discharging in small springs controlled by impermeable strata; and vadose zones above and between the perched aquifers.Large summer rain events in 2006 and 2008 caused rapid water-level response in the perched aquifers, but slower response in the regional aquifer.In the Roswell artesian basin, a shallow unconfined aquifer is present in carbonate strata overlying the Queen-Grayburg anhydrite, at depths less than 100 m below the surface at Artesia.The regional aquifer in this area, at 200-300 m below the surface, is confined beneath the Queen-Grayburg anhydrite and was artesian at the time of first development; subsequent pumping has lowered static water levels by tens of meters (Table 1) [11].Groundwater is present in an unconfined basin-fill aquifer in the alluvium of the Tularosa Valley, where supply wells pump water from the upper 50 m.

Figure 2.
Cross section X-X' (see Figure 1 for location), after Roswell Geological Society (1956).SL = sea level.The east slope of the Sacramento Mountains is a dip-slope with widespread veneer, too thin to depict here, of the Glorieta Sandstone and overlying members of the San Andres formation.

Previous Isotope Studies
Stable oxygen and hydrogen isotope and tritium data were collected for rainwater, surface water and groundwater from Roswell Basin in the late 1970s [12,13], in order to identify sources of recharge and groundwater residence times.The authors concluded that more detailed sampling was required, but were able to identify loci of local, rapid recharge using tritium data in the mountain areas and near Roswell [13].Stable oxygen and hydrogen isotope data for surface water in the Pecos River [14,15], have been used to determine the relative contributions of winter snow and monsoon precipitation to the river in Texas, the authors concluding that the latter predominates [14].Sulfur isotopes in Sacramento Mountains groundwater have been utilized to determine the relative inputs of evaporite gypsum, oxidized sulfides and rain sulfate to the dissolved sulfate inventory [16].Reference [17] provided stable O and H isotope and 14 C data for the well-field supplying water to the air-force base in Tularosa Valley, and interpreted the data to indicate water residence times greater than 1000 years.The most detailed recent work is in reference [9], which presented detailed stable isotope data for precipitation and groundwater collected in 2006-2009 between the range crest and Hope, New Mexico.The authors identified predominant summer recharge in years of heavy summer rainfall, and used tritium, 14 C and CFCs to estimate groundwater residence times of decades in the high mountains, to thousands of years in the aquifer extending east of the range.The pattern of stable O and H isotope data in groundwater presented in reference [9] differs markedly from that in our dataset, allowing for an improved understanding of the hydrology of the mountain range when both datasets are taken into account.Our study also complements reference [9] in extending spatial coverage into flanking basins east and west of the Sacramento Mountains.

Analytical Methods
Samples were taken from domestic, agricultural and municipal production wells, springs, and surface water in the Peñasco and Pecos Rivers, and from rain gauges near Weed.Isotope measurements (except accelerator mass spectrometry carbon-14) were performed at the Environmental Isotope Laboratory, University of Arizona.Stable O, H and C isotopes were measured on a Finnigan Delta S ® dual-inlet mass spectrometer equipped with an automated CO 2 equilibrator (for O) and an automated Cr-reduction furnace (for H).Stable S isotopes were measured on a Thermo Electron Delta Plus XL ® continuous flow mass spectrometer equipped with a Costech ® elemental analyzer for preparation of SO 2 .Carbon-14 was measured by accelerator mass spectrometry at the NSF-Arizona Accelerator Facility, University of Arizona.Data generated for this study are listed in Table 1.Analytical precisions (1σ) are 0.08‰ (O), 0.9‰ (H), 0.15‰ (C) and 0.15‰ (S).Detection limits are 0.6 TU ( 3 H) and 0.2 pMC ( 14 C).

Correction of Raw Carbon-14 Data
The data lack sufficient detail for chemical balance modeling; therefore a simpler method based on δ 13 C values is used, following ( [18], p. 210).Values of δ 13 C for soil gas are assumed to be −23‰ (corresponding to 100% C3 plant matter input) for the forested mountains, and −19.9‰ (corresponding to 75% C 3 input) for desert areas."Dead" rock carbonate of Guadalupian age has δ 13 C values from +1 to +5‰ [19,20]; for these strata corrections were calculated for +1‰ and +4‰.Corrected ages are given in Table 1.In the basin-fill aquifer near Alamogordo, corrections were calculated for rock δ 13 C from 0‰ to +3‰, representing the entire Paleozoic section.

Area 1: High Sacramento Mountains
Samples were collected from wells near La Luz and Fresnal canyons and areas near New Mexico Route 24 east of the range crest (Figure 3).All samples are from fractured limestone except for 1-4, which is from shallow alluvium in Fresnal Canyon.  1, where the corresponding number is 1-1 for area 1, or 2-1 for area 2. Black circles: sample sites for this study; white circles: sample sites from reference [17].

O and H Isotopes
On a plot of δD vs. δ 18 O, most of the data fall on a straight line with a slope near 5.6 (Figure 4A), henceforth called the Sacramento Mountains Trend (SMT).The straight line intersects an estimate of average winter precipitation at a station at 2790 m.a.s.l.(calculated as arithmetic means (because amount data are not available) of δ 18 O and δD for three bulk collections in March 2007, 2008 and 2009, and representing the prior 3 months; data from [9]), but does not intersect mean summer precipitation [9] for that station.In Figure 4B, three groups of δ 18 O values emerge in relation to site altitude.For the wells, collar altitude is used because static water levels are not available in all cases.Values of δ 18 O of group W (western slopes) overlap those of group H1 (high elevations), despite the large altitude difference between the two groups.The difference between groups H2 (high elevations, but generally lower than H1) and H1 is too great to attribute to altitude.Group H2 sites (1-13, 1-14, 1-16, 1-17) are adjacent to broad, flat, canyon bottoms, a typical geomorphic feature of the Sacramento Mountains.In such places, deep soil (more than 1 m near site 1-16) overlies carbonate strata, while elsewhere carbonate outcrop is widespread.The data for groundwater in 2003 differ from data for groundwater in 2006-2009 [9].The latter occupy a field between the SMT and summer rain for 2006 and 2008 (Figure 4A), and reflect rapid recharge from heavy monsoon rains in 2006 and 2008.Prior to 2003, there had been no large monsoon rain totals since 1997.[21], and 1.8‰/1000 m [13].Site numbers (e.g., 3) correspond to entries in Table 1, where the corresponding number is 1-3.

Other Parameters
Groundwater in this area generally has δ 34 S values of 10‰ to 13‰, tritium concentrations of 1 to 3 TU, and 14 C in the range 72 to 93 pMC (cf.0 to 7 TU, and 83 to 93 pMC, for samples from 2006 to 2008 [9]).Corrected 14 C data indicate post-bomb water in the west-slope canyons and in two high-elevation springs, with older groundwater (300-4500 years) at sites 1-15 and 1-18 (Table 1).

Interpretation
The SMT can be explained as an evaporation trend originating in winter precipitation.Evaporation prior to infiltration varies in degree, and is greatest in groundwater near the broad canyon bottoms, (sites 1-13, 1-14, 1-16 and 1-17), where standing water and wet soil are likely to undergo partial evaporation.Well-mixed high altitude groundwater will plot between groups H1 and H2, and this isotope signature will be found in groundwater of the limestone aquifer at lower elevations unless water of different isotope composition is added downgradient.Evaporated runoff from high elevations may plot on the SMT to the right of group H2.Addition of water from local low-elevation precipitation would shift groundwater isotopes towards the GMWL.
The difference between the 2003 and 2006-2009 data sets indicates two modes of recharge.In years with unusually wet summers, (e.g., 2006 and 2008), summer recharge with little evaporation is the dominant source of recharge.The local meteoric water line (LMWL) in Figure 4A is governed by rainfall from those years, and may not apply under drier conditions.Following a succession of dry to average summers, however, winter recharge predominates, even though there is more precipitation in summer than in winter rain on average.Such was the case from summer 1998 to 2003 when sampling for this study occurred.Under these conditions, evaporation of the infiltrating water occurs in the broad canyon bottoms east of the range crest, but is not observed between the range crest and the canyons on the steep west escarpment.
Tritium and corrected 14 C contents of high-elevation groundwater indicate the presence of post-bomb recharge, but tritium levels in 2003 were predominantly lower than average tritium in post-1992 precipitation (4-7 TU, see Table 1 and [9]; compare a better-constrained average of 5.3 TU for Tucson [22]), indicating mixing with pre-bomb meteoric water.By 2006-2008, more post-bomb recharge was present, tritium-helium dates were mainly 1-15 years, and CFC ages were largely 20-30 years [9].Values of δ 34 S indicate Permian marine gypsum (+12‰ to +13‰) as the main source of sulfate; lower values most likely reflect oxidation of sulfide present in these strata [16].

Area 1-Alamogordo and Tularosa
Sampling from supply wells in basin-fill alluvium represents lower-TDS water suitable for human consumption; brackish water is also present >5 km west of the range front.Groundwater in this area flows west at Alamogordo and Tularosa [4,5], but parallel to the range front south of Alamogordo, where no major canyons contribute water to the basin.

O and H Isotopes
Most data plot on the SMT (Figure 5), to the right of group W. Samples from Tularosa include the most and least evaporated of the set.Data for wells south of site 1-33 (δ 18 O between −9.9‰ and −9.5‰ in reference [17]) differ from data collected for the present study in the same area (δ 18 O between -9.5‰ and -9.0‰).Actual variation in δ 18 O (as opposed to measurement error) is unlikely in such old groundwater (see below); the earlier data are not used here.In southeastern Tularosa Valley, sites 1-39 and 1-40 (Figure 1) have groundwater that plots below the SMT Figure 5. Plot of δ 18 O vs. δD for groundwater samples from basin sediments near Alamogordo and Tularosa, in relation to samples from La Luz and nearby canyons and the high Sacramento Mountains.Samples 1-39 and 1-40 are from basin fill more than 40 km south of Alamogordo (see Figures 1 and 3).

Other Parameters
Tritium is present (1-3 TU) north of site 1-33, and is generally absent (below detection to 0.5 TU) south of 1-33. 14C generally decreases from near 80 pMC near Tularosa to 20 pMC south of Alamogordo (Figure 6).Corrected 14 C data indicate young groundwater (post-bomb to a few hundred years) north of Alamogordo and in La Luz canyon, and much older water (500-7500 years, considering also corrected data from [17]) south of Alamogordo.Values of δ 34 S are near +12‰.At sites 1-39 and 1-40, tritium is below detection, 14 C levels are 3 pMC, and corrected ages are 12,000 to 21,000 years (Table 1).

Interpretation
O and H isotope data plotting on the SMT indicate high-elevation precipitation as the source of groundwater in basin alluvium near Alamogordo and Tularosa.Groundwater from the high Sacramento Mts.flows to La Luz Canyon sample sites without isotopic shift.The higher degree of evaporation in samples from the alluvial aquifer could be explained: (1) as mountain-block recharge combining more-evaporated and less-evaporated recharge from high elevations; or (2) as mountain-front recharge of surface water supplied from high elevations by way of the mountain canyons.The absence of an evaporation signature in groundwater from carbonate strata in La Luz Canyon, between the range crest and the basin) argues against the first possibility, while the presence of evaporated water in the alluvium argues for the second.The higher degree of evaporation of groundwater farther from the range front (Tularosa, 12-15 km from the range front), in contrast to groundwater nearer to the range front (Alamogordo, within 6 km), suggests that the sites of infiltration of surface water extend into the basin, rather than being confined to a narrow zone at the range front.This is particularly evident in the case of site 1-22 at Tularosa, (Figure 5), where the coincidence of low δD and δ 18 O with high 14 C indicates recharge of very recent runoff at a distance of up to 15 km from the range front.Both mountain-front and mountain-block recharge seem likely, but the data do not indicate the relative amounts.In the basin fill south of Alamogordo, tritium and 14 C data are consistent with slow southward flow of groundwater, with little recharge from nearby mountain canyons.

Area 2: Peñasco to Artesia
Samples were collected between Peñasco and Artesia (Figure 3).Near Peñasco, groundwater samples were taken from a spring and a windmill in limestone, and from wells in the Rio Peñasco flood plain.East of the range front, as far as site 2-10, an unconfined aquifer (the principal aquifer of [11]) is present near the boundary of the Yeso and San Andres formations.Recharge to these strata may occur near the range front.Samples are from domestic and agricultural wells up to 260 m deep, with static water levels (SWLs) near 190 m below the surface.East of site 2-10, beneath a broad plain west of the Pecos River, two major aquifers were sampled.An unconfined aquifer with SWLs from 30 to 60 m below the surface exists in flood-plain sediments near Artesia.The eastward continuation of the principal aquifer, 275 to 300 m below the surface, is confined beneath the Queen-Grayburg anhydrite (Figure 2).It was artesian at the time of first exploitation; SWLs at present range from 20 to 60m below the surface.Surface water samples were collected from the Peñasco and Pecos Rivers.
From Peñasco to Hope, groundwater flow is east-southeast (Figure 18 of [9]).Allowing for variation due to pumping, SWLs in the principal aquifer east of Hope are close to 1000 m.a.s.l.(Table 1).Southward flow is likely in this area.

O and H Isotopes
All groundwater and surface water samples plot close to the SMT (Figure 7A).Most data cluster at the intersection of the SMT with the GMWL, where the separation of the lines is less than 2‰ in δD, and therefore impossible to resolve within measurement error.Data for the principal aquifer from Mayhill to Hope [9] match the present data set in δD but include lower values of δ 18 O.Two groundwater samples from near Artesia (2-13, 2-19) plot to the right of the main data cluster.Surface water from the Pecos River in the reaches between Artesia and Red Bluff ( [15] for 1984-1987, [14] for 2005, and data from this study) largely plot as a linear trend, close to an extrapolated SMT (Figure 7A,B).
Previous data for the principal aquifer [13] pre-date automated isotope methods, and partially overlap the main data cluster from the present study.The two data sets correspond in δ 18 O, but the older δD data have a spread >20‰, apparently spurious, and appear not to be useful.Weighted precipitation averages from [12] are for δ 18 O alone, and have been plotted on the GMWL in Figure 7.

Tritium
In 1977-1978, when bomb tritium averaged about 35 TU in local precipitation, surface water and alluvial groundwater from the Peñasco River flood plain contained about 10 TU, and tritium in the principal aquifer near Artesia was below detection [13].In samples collected for this study, tritium is present at low levels (<2 TU) in groundwater from near the range front (sites 2-2, 2-3, 2-5), at site 2-19 in the shallow aquifer, and in one deep aquifer sample (site 2-21, 0.5 TU); at other sites it is below detection (Figure 8).Groundwater from the alluvial aquifer beneath the Peñasco River (site 2-3) contains 1.5 TU, distinctly lower than the average for precipitation, and consistent (cf.[13]), with a large pre-bomb groundwater contribution to the surface water of the river.Reference [9] listed tritium contents <2 TU in groundwater between Elk and Hope.  1O vs. δD for groundwater and surface water samples from study area 2. (A) All data from this study.The field of data from [13] is for the principal aquifer from Artesia to Roswell, and encompasses all but three outlying data points.The inset shows a magnified view of clustered data.Site numbers (e.g., 1) correspond to entries in Table 1, where the corresponding number is 2-1; (B) Plot of δ 18 O vs. δD for the Pecos River between Artesia and Red Bluff, from other studies [14,15], relative to the SMT.showing well depths and static water levels (SWL) in relation to the surface."Shallow" refers to the shallow aquifer at Artesia, and "deep" to the deeper artesian aquifer.Also shown are measurements of carbon-14 (pMC) and tritium (TU).Tritium below detection is indicated as "bd".

Interpretation
Groundwater in the Pecos Slope and artesian aquifers is largely uniform in isotope content over an east-west extent of about 100 km, and lies on or close to the SMT.A dominant water source in the high Sacramento Mountains is therefore likely.East of Peñasco, a few samples (2-13, 2-19, 2-24 and 2-25) have isotope data plotting on the SMT, but to the right of the main data cluster; these may reflect local recharge of evaporated surface water.Site 2-13 is close to the ephemeral lower reach of the Rio Peñasco, where recharge of evaporated surface water may occur.The other three samples are from the shallow aquifer beneath irrigated fields, where reflux of evaporated irrigation water is probable.Addition of local rainwater is likely for site 2-10 (Figure 7).
Bulk groundwater residence times (the corrected versions of the data shown in Figure 8) in the principal aquifer are 1300 to 5600 years east of Hope, greater than those suggested in [13].
Most of the δ 18 O and δD data for the Pecos River between Artesia and Red Bluff plot on a linear trend close to an extrapolated SMT, regardless of season (Figure 7B), and can be therefore be generated as a result of evaporation of water like that in the principal aquifer at Artesia.The principal source of river water in this area is therefore most likely the Sacramento Mountains, either by natural recharge from the aquifer, or by way of irrigation on the Pecos flood plain.If this is true, mountain-derived water is discharged with an isotope signature of evaporation into the river near Roswell and Artesia.This suggests a modification to the modeling, based on deuterium excess, of river water sources in reference [14].

Water Sources in Study Area
Most water sampled for this study plots on the Sacramento Mountains trend (SMT) in δD-δ 18 O space.The SMT originates in high-altitude winter precipitation.Such precipitation is therefore the principal and ultimate source of groundwater in the area between Alamogordo and the Pecos River.Most water in the Pecos River near Artesia also appears to be of that origin.There is scant evidence for recharge of local meteoric water at low altitudes.The few exceptions include groundwater from the southeastern part of Tularosa Valley (where ancient water from high elevations appears to be present), and some unconfined-aquifer samples from Artesia (where local recharge probably occurs).

Seasonality of Recharge
The heavy monsoon rains of 2006 and 2008 generated recharge of distinctive δD-δ 18 O signature in groundwater of the high Sacramento Mountains, but in drier years, 2007 and 2009, groundwater isotopes shifted towards the SMT [9].Monsoon rainfall comparable to that in 2006 had not occurred since 1997, and in the 2003 sampling for this study, winter recharge, plotting on the SMT as a result of local evaporation prior to infiltration, was predominant.Where old groundwater is present (south of Alamogordo and in the principal aquifer of Roswell Basin), δD and δ 18 O conform largely to the SMT.In the long term, therefore, recharge in dry to average years contributes the larger volume to low-elevation aquifers around the Sacramento Mountains.Years with unusually wet summers lead to a transient (a few years) response in the high-altitude aquifers, but make little contribution to the old groundwater in basins at the foot of the mountains.
The Sacramento Mountains are therefore an unusual example of a mountain block in which the dominant season of recharge can change in response to seasonal precipitation amounts, although winter precipitation, only 35% of total precipitation on average, predominates in the long-term.Winter recharge is considered predominant in a number of other mountain ranges in the arid western USA.In the Spring Mountains, Nevada, another carbonate-rock range, winter precipitation is dominant; summer rain contributes about 30% of annual precipitation, but only about 10% of recharge [23].Winter recharge also predominates in the Huachuca Mountains, Arizona, where summer precipitation contributes 54% of the annual total on average, but winter precipitation accounts for 65% ± 25% of recharge [24].In the Santa Catalina Mountains, Arizona [25], and the ranges delimiting the Verde River watershed, Arizona [26], winter recharge is considered predominant, contributing 98% of recharge in the latter case.

Sacramento Mountain Carbonate Strata as a Karst Aquifer
A continuum of aquifers exists in carbonate rock [27].At one extreme, carbonate dissolution leads to wide solution cavities that self-organize into dendritic drainage networks discharging through large springs; water flow rates are commonly 10 2 to 10 4 m/day over distances of 10 3 to 10 4 m.At the other extreme, solution cavities are narrow and of limited interconnection, generating an aquifer with lower flow rates and discharge through many small springs.Although small-scale collapse structures are recognized [9], cavern networks are not developed in the thinly bedded strata, some impermeable, of the study area.A flow rate of 10 m/day over 30 km between the range crest and the eastern range front would result in water travel times of about 8 years.The tritium and 14 C data imply residence times >60 years in the mountain aquifers several cases, while surface water in the Rio Peñasco and associated flood-plain groundwater also contain some pre-bomb precipitation.The isotope evidence indicates widespread persistence of pre-bomb precipitation in groundwater, and flow rates typically much lower than 10 m/day.The Sacramento Mountains therefore fall at the latter end of the continuum of karst aquifers as described above.
Nonetheless, the carbonate strata in and east of the Sacramento Mountains compose a regional aquifer system over a distance of 130 km.Regional carbonate aquifers of similar extent have been demonstrated elsewhere in the region on the basis of geochemical modeling, east of the Salt Basin in West Texas [28], and in the region southwest of the Cuatrocienegas Basin of Coahuila, Mexico [29].

Mode of Mountain-Front Recharge to Basin Alluvium
The location of mountain-front recharge relative to the interface between hard-rock mountain blocks and basin alluvium in the southwest USA has been addressed in several studies.In the middle Rio Grande Basin (New Mexico), infiltration is thought to occur in a narrow zone along the range front of the Sandia and Manzano Mountains [30].In Chino Valley (Arizona) [26] and Tucson Basin (Arizona) [1], evidence indicates infiltration from stream beds downstream of the mountain fronts, at distances of 6 to 10 km in the case of Tucson Basin.Groundwater isotope data also indicate recharge of ponded surface water in the center of the Hueco Bolson (Texas) [31].The present study concurs with the possibility of infiltration as far as 15 km downstream of the range front.

Source of Hueco Bolson Groundwater
The question addressed here is the source of saline water in the center of the Hueco Bolson, an alluvial basin near El Paso, Texas, 100-130 km southwest of Alamogordo (Figure 1).Subsurface movement of groundwater from the Tularosa Valley to the Hueco Bolson is physically possible according to piezometric data [4,6].An alternative source is recharge from the Organ and Franklin Mountains (Figure 1), which supply a freshwater aquifer, the Franklin Mountains freshwater lens (FMFWL) in ancient fluvial deposits at the western edge of the Hueco Bolson [6].The catchment for the FMFWL is largely at altitudes between 1300 and 2400 m.a.s.l. in the Organ Mountains, in contrast to a catchment at 2400-2800 m.a.s.l. for the four large canyons that focus fresh water from the Sacramento Mountains into basin sediment near Alamogordo.H and O isotopes might therefore discriminate between the two sources, as discussed inconclusively in [31].
Groundwater from the FMFWL plots along the global meteoric water line with δ 18 O values between -9‰ and -11‰ (Figure 9).The upper end of the data array corresponds to groundwater of short residence time, while the lower end corresponds largely to groundwater resident for thousands of years [2].The SMT and the suggested paleo-SMT (based on samples 1-39 and 1-40) intersect the FMFWL trend near -9‰ and -11‰, respectively.On the one hand, the δ 18 O and. δD values of the saline, evaporated water in the center of the Hueco Bolson plot between the SMT and the paleo-SMT, and could therefore represent mixtures of older and younger water from the Sacramento Mountains.On the other hand, δ 18 O and. δD values define an evaporation trend that could originate in older FMFWL water, so that the water could have originated in the Frankin and Organ Mountains, perhaps as surface water ponded and evaporated in the basin center at a time of cooler, wetter climate.The stable isotopes fail to distinguish the two possibilities because of the likely presence of ancient groundwater.4) and data for groundwater at Alamogordo; (b) Data for groundwater in the Hueco Bolson in Texas, from [2] and [24], distinguished according to salinity (HB-saline at the basin center, and HB-fresh from the Franklin Mountains fresh water lens on the western side of the basin); (c) A suggested paleo-SMT based on two samples of ancient water in the southeastern part of Tularosa Valley.

Implications for Groundwater Management
High-elevation winter recharge is the principal source of groundwater over the long term in the aquifers of the Sacramento Mountains and the flanking basins.If winter precipitation declines, for instance in response to global climate change, groundwater supply will decrease.The effect would be felt initially in the high mountain communities such as Weed and Cloudcroft (but might be mitigated if occasional wet summers persist) and in areas from La Luz Canyon to Alamogordo and Tularosa where groundwater is of post-bomb age (Table 1).In the Roswell basin, where artesian water has been resident for thousands of years, there would be no short-term effect of diminished winter recharge; over-pumping for irrigation would be of more immediate concern.

Figure 4 .
Figure 4. (A) Plot of δ18 O vs. δD for groundwater samples from the high Sacramento Mountains.The green line encloses groundwater isotope data from[9].Seasonal means for precipitation and the local meteoric water line (LMWL) are for years 2006-2009[9].Data plotted as individual points were collected for this study in 2003; (B) Plot of elevation of well collars vs. δ18 O for sample sites in the high Sacramento Mountains.The diagonal lines show the long-term δ 18 O lapse-rates of −1.2‰/1000 m (Tucson Basin[21], and 1.8‰/1000 m[13].Site numbers (e.g., 3) correspond to entries in Table1, where the corresponding number is 1-3.

Figure 7 .
Figure 7. (A) Plot of δ18 O vs. δD for groundwater and surface water samples from study area 2. (A) All data from this study.The field of data from[13] is for the principal aquifer from Artesia to Roswell, and encompasses all but three outlying data points.The inset shows a magnified view of clustered data.Site numbers (e.g., 1) correspond to entries in Table1, where the corresponding number is 2-1; (B) Plot of δ18 O vs. δD for the Pecos River between Artesia and Red Bluff, from other studies[14,15], relative to the SMT.

Figure 8 .
Figure 8. A. East-west profile of Area 2 (refer to Figure2for location) showing well depths and static water levels (SWL) in relation to the surface."Shallow" refers to the shallow aquifer at Artesia, and "deep" to the deeper artesian aquifer.Also shown are measurements of carbon-14 (pMC) and tritium (TU).Tritium below detection is indicated as "bd".

Figure 9 .
Figure 9. Plot of δ 18 O vs. δD showing: (a) The Sacramento Mountains Trend (SMT, as in Figure4) and data for groundwater at Alamogordo; (b) Data for groundwater in the Hueco Bolson in Texas, from[2] and[24], distinguished according to salinity (HB-saline at the basin center, and HB-fresh from the Franklin Mountains fresh water lens on the western side of the basin); (c) A suggested paleo-SMT based on two samples of ancient water in the southeastern part of Tularosa Valley.

Table 1 .
Site information and isotope data.
Notes: S = shallow aquifer, D = Deep (Principal) aquifer in Artesia area; na = not available; A = Apparent tritium; * meters below surface.