Geology topics

Research about New Mexico

Source-linked reports with geographic coverage including New Mexico.

At least 19 recordsLinked to original sources

USGS2026

Humans experienced abrupt environmental changes at White Sands during the Last Glacial Maximum

Expanded excavations at White Sands National Park have unearthed 20 stratigraphic horizons containing in situ human footprints, megafaunal trackways, and linear features created by travois. Advances in radiocarbon dating of conifer pollen, combined with ages of aquatic seeds, refines the chronology of the footprint-bearing sequence to between 23.5 and 21.7 thousand years ago (ka). Biomarker, pollen, and stratigraphic evidence indicate a complete overturning of the environment from a perennial lake to a playa and capture a large aridification episode coincident with Dansgaard-Oeschger event 2 at 23.2 ka, which shrank pluvial Lake Otero and allowed people to traverse the former shoreline area for nearly two millennia. This well-dated and highly resolved paleoclimate record provides compelling evidence that humans in North America experienced marked environmental changes during the Last Glacial Maximum.

New Mexico↗

USGS2026

Fire regime shift reveals 17th century socio-ecological change in the Rio Grande del Norte National Monument, New Mexico

Humans have influenced ecosystems for millennia, and social changes can have notable ecological effects. Recent human influence on ecosystems is exemplified by the collapse of fire regimes circa 1900 in the western United States due to increased grazing, fire suppression, and declining Indigenous fire use. Climate can also drive fire regime changes, which is reflected by increasing area burned across the West in recent decades associated with warming and increased aridity. In the southwestern United States during the mid-1600s, large social changes accompanied Spanish colonization. These changes included increases in domestic livestock (e.g., sheep and horses), declines of Indigenous populations and cultural burning, and exploitation of Puebloans that led to the 1680 Pueblo Revolt, an Indigenous uprising that expelled Spanish colonizers for 12 years. The ecological effects of these large social changes are not well studied. The goal of this study was to use tree-ring fire histories to investigate fire dynamics of a sagebrush-dominated landscape in the Rio Grande del Norte National Monument (RGDN) in northern New Mexico, including the effects of historical social and climatic changes. We reconstructed fire history on four forested mountains that rise above the sagebrush/grass plateau. We asked: (1) What was the pre-1900 fire frequency, season, and fire–climate relationship in the RGDN? (2) Was there evidence for fire spread in the sagebrush/grass flats between the mountains? and (3) Did regional social or climate changes influence fire regimes? From 1400 to 1679, we found frequent fire within forests (8–12-year median fire intervals [mfi]) and potential evidence of spreading fire through the sagebrush/grass ecosystem based on greater-than-expected intermountain fire synchrony. However, beginning in the mid-1600s and continuing after the Pueblo Revolt of 1680, there was (1) a significant decline in fire frequency, including 38–166-year mfi and no fire for >350 years at two sites; and (2) only one synchronous fire between mountains. During this 17th century period of socio-ecological change there was no persistent shift in climate. We propose that reduced Indigenous ignitions and increased domestic grazing, which reduced herbaceous fuels and fire spread, are likely causes of the extraordinary early and persistent reduction in fire.

New Mexico↗

USGS2026

Recent high-severity wildfires in a dry-conifer landscape are unprecedented over five centuries and foretell future forest loss

Dry forests of the western United States evolved with frequent, spatially extensive, low-severity fires, followed by a century of fire exclusion. Recently, these same forests have experienced large and increasingly severe fires causing widespread forest loss. To what extent do modern high-severity fires depart from pre-1900 fire regimes, and how do these departures affect the future persistence of dry-conifer forests? We address these questions using a multiscalar approach that combines the largest landscape-scale tree-ring fire-scar network in North America with neutral landscape model simulations in the Jemez Mountains of northern New Mexico. Fire scars recorded a continuous 400-y history (ca. 1500–1900 CE) of frequent, low-severity fires across individual sites, landscapes, and centuries. We evaluated whether forests exhibiting this record could plausibly persist if they had historically burned under modern patterns of high-severity fire. Across the 240,000-ha Jemez landscape, high-severity fires from 1995–2024 burned 32% of fire-scar sites and 25% of dry-conifer forests. Within the Las Conchas Fire perimeter—a 61,000 ha area reburned by multiple fires—66% of fire-scar sites and 71% of forests were burned at high severity, erasing multiple centuries of forest and fire history. Simulations indicate that the probability of the observed fire-scar network persisting under modern high-severity fire rates is near-zero ( P < 0.001). Projecting these dynamics forward indicates rapid loss of both dry-conifer forests and tree-ring fire histories—50% of modern dry-conifer forest will remain by 2055, and 24% by 2115. Modern high-severity fire regimes appear incompatible with dry-conifer forest persistence.

New Mexico↗

USGS2026

Development of a hydrologic flow model of the Zuñi Salt Lake and surrounding areas, west-central New Mexico

The terminal Zuñi Salt Lake is in a maar in west-central New Mexico and contains hypersaline water that has long been used by Native Americans for religious purposes and the collection of salt. The U.S. Geological Survey (USGS), in cooperation with the Zuni Tribe of the Zuni Reservation, New Mexico, and Bureau of Reclamation, developed the Zuñi Salt Lake hydrologic flow model to simulate the steady-state conditions of the lake and surrounding groundwater-flow system. The model is a tool that can be used to analyze the potential hydrologic responses to resource development in the region under different water-use scenarios, which can support negotiations that could help protect future use of the lake. The USGS modular hydrologic model, MODFLOW 6, was used to simulate the hydrologic system of the lake and surrounding areas. To define the layering of the model, a three-dimensional hydrogeologic framework was constructed, which included seven informal hydrogeologic units and the locations of faults that truncate and offset the units. Calibration of the model was performed using the iterative ensemble smoother included in a parameter estimation software suite (PEST++). The iterative ensemble smoother approach resulted in a posterior parameter ensemble of model realizations with parameter and simulated values that show similar success in matching observations. Groundwater flux to the Zuñi Salt Lake simulated by using the posterior parameter distribution ranged from 26 to 730 acre-feet per year with a mean, median, and standard deviation of 530, 570, and 170 acre-feet per year, respectively. The relative contribution of groundwater discharge to the lake was simulated using a particle-tracking model for MODFLOW (MODPATH) to be from, in order of decreasing contribution, the Mesaverde, Cenozoic, Mancos, and Dakota hydrogeologic units. Simulations of future conditions using the ensemble of models can provide an empirical estimate of forecast uncertainty without substantial additional analysis.

New Mexico↗

USGS2026

3D Elevation Program—New Mexico

Using 3D Elevation Program data, a modern illustration of each State's landscape showcases its elevation differences in vivid detail. Each State has been consistently adjusted and colorized according to its highest and lowest elevations, accentuating differences even in States with low relief. To further emphasize landforms, shaded relief has been added on top of the elevation colors.

New Mexico↗

USGS2026

Beaver dam analogs as nature-based solutions to mitigate snowpack loss in northern New Mexico

Reductions in snow have left many streams in northern New Mexico dry or with very low flows during the summer. Base flow, defined as the contribution of groundwater to streamflow, can sustain streamflow during the summer and during drought conditions. Beaver dam analogs can be used to “slow the flow,” or increase the infiltration of rain-based runoff to replenish shallow groundwater reservoirs and increase base flow and potentially increase summertime streamflow.

New Mexico↗

USGS2026

Mass mortality of avian migrants in New Mexico, USA, that coincided with an extreme weather event

Many birds are migratory, and this life history strategy allows for maximized access to seasonally abundant resources and favorable climates. However, migration exposes birds to threats and stressors, resulting in high mortality during migration. Anthropogenic landscape alterations and climate change have intensified threats, and mass mortality events linked to extreme weather are more common in recent decades. Documenting mass mortality is critical for predicting future occurrences and implementing effective conservation. Here, we describe a mortality event that occurred throughout New Mexico, USA in fall 2020. Carcasses began appearing in the region in mid-August, during a period of extreme heat and drought. Following an extreme cold weather event on September 8–9 th , the number of carcasses increased dramatically and expanded throughout the state. In total, we collected 628 carcasses comprising 58 species within Doña Ana and Otero Counties in New Mexico. Necropsy determined emaciation was the cause of death for 74.6% of carcasses. Live birds captured during the period of peak mortality (n = 223) were in similarly poor condition. This event provides a striking example of how multiple types of extreme stressors, in this case widespread drought and unseasonal cold, coincided with a mortality event, indicating a possible synergistic relationship between these factors and the mass mortality. Mortality events are likely to increase in frequency with intensifying climate change. Establishment of networks of biologists and researchers could improve our ability to identify and communicate developing mortality events, organize data collection, and improve understanding of the causes and consequences of mortality.

New Mexico↗

USGS2026

Fossil footprints and Ice Age ecosystems of White Sands National Park

Introduction In September 2021, National Park Service staff, U.S. Geological Survey scientists, and an international team of researchers revealed evidence in the form of human footprints at White Sands National Park, New Mexico, that showed people were present in North America between 23,000 and 21,000 years ago. This time was during the Last Glacial Maximum, when large ice sheets covered much of the continent. The results stunned the scientific community and sparked a global debate. The story of how the discoveries were made, how they upended traditional thought, and how they “rewrote the book” on the earliest phases of North American prehistory is a classic example of the process of science.

New Mexico↗

USGS2026

Relationship of satellite-derived fire severity to archaeological fire severity and fire effects in the Jemez Mountains, New Mexico, USA

Background Archaeological attributes important for determining cultural affiliation, site age, or artifact origin can be damaged or destroyed by wildfires. Satellite-derived fire severity products are used in post-fire rapid response to prioritize and guide stabilization and preservation activities, but their ability to predict archaeological fire impacts has not yet been tested. We compared the severity of fire from Burned Area Reflectance Classification (BARC) maps with ground-based observations of fire severity (“archaeological fire severity”) and archaeological impacts for 904 prehistoric or historic sites within three wildfires in the Jemez Mountains of north-central New Mexico, U.S. Results from this unique study inform archaeological preservation in the context of wildfires and fire management. Results We found a statistically significant, moderately strong association between BARC and archaeological fire severity and a high likelihood that sites with a low BARC fire severity also had a low archaeological fire severity designation. Sites of moderate or high BARC severity were less likely to have a matching severity classification in post-fire assessments, and BARC maps tended to identify as unburned areas where archaeologists determined that sites were burned at low severity. Wildfire of any severity resulted in measurable impacts to archaeological architecture, features, or artifacts at some sites, with increasing proportion of sites impacted at moderate versus low BARC severity, but fire impacts were also observed at sites classified by BARC as unburned. Conclusions Generally, BARC maps provide good assessment of the severity of fire and its effects on archaeological resources, particularly for areas burned at low severity. In areas of moderate or high fire severity, BARC maps may not consistently identify sites with a high need for assessment, and archaeological fire effects can also be present at sites classified by BARC as unburned. Adoption of standardized post-fire archaeological assessments across wildfires and regions, exploration of alternative fire mapping products, and continued research on the fire environments that produce archaeological fire effects can enhance archaeological preservation and post-fire response.

New Mexico↗

USGS2026

Living with wildfire in Santa Fe, New Mexico: 2024 data report

Community wildfire readiness includes homeowner wildfire risk mitigation and wildfire evacuation preparedness. This report presents results from a household survey distributed to homeowners in the study area around Santa Fe, New Mexico in 2024 which was a follow-up to a previous survey of homeowners of the same residences, consisting of mostly identical questions, that was conducted in 2021. This follow-up survey was motivated by the nearby 2022 Hermit’s Peak/Calf Canyon Fire. Comparison of results indicates general stability in the aggregated responses to most survey questions. As found previously, 2024 results indicate that Santa Fe survey respondents are engaged in preparing for wildfire, yet more could be done to reduce risk in study communities. Despite decreases between 2021 and 2024 in the acceptability of fuels treatment methods that concern the use of fire, respondents generally support multiple types of programs intended to reduce wildfire risk to the city, including most but not all types of fuels treatments on public lands.

New Mexico↗

USGS2026

Hydrologic dynamics of ephemerally flooded playas in a dryland environment

Ephemerally flooded playas are common in the southwestern United States and globally in drylands. Often formed in closed basins, playas are depressions which inundate infrequently from local precipitation and streamflow produced near the playa or from upland areas. Few studies have quantified the hydrologic connectivity between upland catchments and playas using observations. Here, we used rain gauge-corrected precipitation from weather radar and water level measurements in 18 playas of the Chihuahuan Desert to identify precipitation thresholds leading to playa inundation over a 6.4-year period. Geospatial data sets on topography, soil properties, and vegetation cover were employed to determine the controls on inundation. Only 9.4% of all precipitation events above 1 mm led to inundation, with 69.8% of all inundations occurring during the North American monsoon (NAM, July-September). Mean and standard deviations (Std) of runoff ratios at all playas were 2.74 ± 4.08% and 3.29 ± 5.19% for annual and NAM periods. At the annual scale, playa inundation occurred when mean precipitation thresholds of 18.3 ± 7.5 mm (event total) and 12.0 ± 4.5 mm/hr (60-min intensity) were exceeded. Across all playas, inundation occurrence and volume were related most strongly to precipitation metrics and catchment area, with secondary controls of soil and terrain properties. The explanatory power of the derived regressions describing the inundation response across the playas were significantly improved when considering their geological origin. As a result, the inundation response classification system could be applied to ephemeral playas in other arid and semiarid landscapes.

New Mexico↗

USGS2026

Patterns of rift basin development and the fidelity of the subsidence record: Insights through Bayesian modeling of rapid tectonic subsidence in a Rio Grande rift basin, Socorro, NM, U.S.A

Characterizing tectonic subsidence rates within depositional sequences provides direct insight into the driving mechanism(s) of accommodation in a basin. However, the temporal resolution of this record is often stymied by a lack of high-precision and high-resolution ages, which enable a more complete description of basin subsidence drivers. We explore the effect of high-precision and high-resolution ages in modeling accommodation for the Miocene La Jencia Basin of the central Rio Grande rift (RGR) and interpret driving mechanism behavior from these models (e.g., lithospheric thinning). We present a new geochronologic dataset of both laser ablation-inductively coupled plasma mass spectrometry (LA-ICPMS) and chemical abrasion-isotope dilution-thermal ionization mass spectrometry (CA-ID-TIMS) data and use these ages in Bayesian accommodation modeling. Models constrained by high-precision and high-accuracy TIMS ages yield peak tectonic subsidence rates exceeding 220 m/Myr, and an average Miocene subsidence of ∼120 m/Myr. While timing and magnitude vary, all models suggest two pulses of rapid Miocene tectonic subsidence, which we interpret to reflect basin-bounding fault movement. Prior to peak subsidence, there was an initial period of fault linkage and organization that occurred over <1–3 Myr that produced the basin-bounding La Jencia-Cerro Colorado fault zone. A comparison of published tectonic subsidence rates to those modeled here shows that while tectonic subsidence during continental rifting is highly variable, the La Jencia Basin rates appear relatively high. However, the significant difference between peak and average La Jencia Basin rates modeled here highlights the potential for underestimation of many records of tectonic subsidence due to a lack of high-precision and high-resolution age constraints. Furthermore, age data and modeling results presented here document fault movement and consequent rates of tectonic subsidence that lower-resolution data would not, providing a high-fidelity case study of continental rift basin development.

New Mexico↗

USGS2026

Valuing socio-economic and ecological attributes of forested watershed restoration to reduce wildfire risk in the southwestern U.S.

Forest restoration in a watershed can provide numerous ecological improvements and social benefits, including reducing the risk of extreme wildfire. Understanding the values of the accrued benefits can be used to evaluate the use of funds to support restoration. The Rio Grande watershed is a vast watershed covering approximately 335,000 mile 2 (867,646 km 2 ). The Rio Grande watershed provides a host of ecosystem services and recreation opportunities and supports municipal water supplies. We estimate the non-market values of forest restoration in the Rio Grande watershed using a choice experiment (CE) approach. While the ecological benefits are established in literature, we focus on the social characteristics, capturing the human-forest system in a comprehensive manner. Our results indicate a significant willingness to pay (WTP) for improving air quality, reducing private property damage, and creating local jobs, with the highest WTP for job creation. Split-sample analysis indicates respondents residing within the watershed highly value the socio-economic attributes of the restoration, while ecological attributes are preferred more outside of the watershed. Our results provide insights into the benefits of multi-dimensional services from forest restoration activities in a watershed.

New Mexico↗

USGS2026

Carcass size and ground substrate drive detection rates of avian carcasses by human surveyors and a dog team

Accurate avian mortality estimates are essential for understanding anthropogenic impacts to bird populations and informing conservation strategies. Carcass surveys are commonly conducted by human surveyors or by detection dogs, but the factors influencing surveyor detection abilities have not been fully explored. In this study, we conducted two years of detection trials in the semi-arid high desert of southern New Mexico, USA, testing 27 human surveyors and one conservation detection dog across 1096 trials with 238 carcasses representing 50 avian species. We directly compared detection abilities between surveyor types (human and dog) and identified key factors influencing detection probabilities. The conservation detection dog exhibited a significantly higher detection probability (mean = 0.87) than human surveyors (mean = 0.49, individuals ranged 0.25–0.71), consistent with previous studies. Detection probabilities for both surveyor types were influenced by carcass size and ground substrate; detection probability was higher for larger carcasses and areas with lower vegetative complexity. We discuss our results in the context of common tradeoffs faced by managers in designing carcass surveys and how guidance may vary under different scenarios. Broadly, our study provides valuable insight that can enhance wildlife mortality monitoring, ensuring more accurate mortality estimates to inform management and conservation efforts.

New Mexico↗

USGS2026

Cultural pollen delivery is not needed to account for consistent, trace abundances of Zea mays pollen in the sediments of Santa Fe Lake

Roos (2026) proposes an intriguing cultural mechanism to account for the rare (<1.5% of terrestrial pollen), but consistent (28/34 samples spanning 450 years), presence of Zea mays (maize) pollen in the sediments of Santa Fe Lake from 800 to 350 cal yr BP (1150 – 1600 cal yr CE; Henne et al., 2026a, 2026b). Whereas the importance of maize agriculture to Ancestral Puebloans is well established, Roos argues that the near-absence of maize pollen from soils proximate to historical maize cultivation and in peat or lake-sediment records in the Southwestern USA indicates that natural pollen delivery was not sufficient to account for the trace abundances of maize pollen we counted at Santa Fe Lake. Roos concludes that ceremonial practices provide an alternative mechanism that should be considered in the interpretation of lake-sediment archives like Santa Fe Lake. We agree that the cultural practices described by Roos are a credible source of pollen delivery and understand the significance of sacred lakes to Indigenous communities (Bodine, 1988). However, we conclude that natural delivery to the Santa Fe catchment is sufficient to account for the trace abundances of maize pollen in our record. We first place our consistent finding of maize pollen in a context beyond the Southwestern USA that includes other mountain regions, discuss why it is unlikely to find consistent maize pollen in the records presented by Roos, and expand on what makes Santa Fe Lake a superior archive for maize pollen than previous records. Based on this discussion, we conclude that natural pollen deposition is the most parsimonious explanation for consistent maize pollen over a period of 450 years.

New Mexico↗

USGS2026

Regreening, restoring, and reconnecting a southwestern wetland ecosystem – the Zeedyk wetland

Alluvial wetland ecosystems are vital as biodiversity hotspots but are increasingly threatened by anthropogenic stressors and drought. These pressures are especially acute in arid and semi-arid regions, where eco-hydrologic connectivity is fragile and recovery is slow. This study quantifies the efficacy of nature-based solutions, particularly the ‘Zeedyk approach,’ which employs low-tech Natural Infrastructure in Dryland Streams (NIDS)—including rock detention structures—to slow surface water, raise groundwater tables, and restore wetland function at a spring-fed wetland in Cebolla Canyon, New Mexico, U.S.A. Our results depict a Restoration Feedback Loop that captures stages of change from a healthy wetland in 1935, altered by 20th-century agriculture and grazing, to the re-establishment of the historical flow regime by 2024 documented through an 89-year archive of aerial imagery (1935–2024). By the end of our study period, the Spring-Fed Wetland had expanded by roughly 229% of the original 1935 area, to 4.13 ha. Using 40 years of satellite data, we assess changes in vegetation and hydrology with remote sensing indices. Spatial and temporal analyses reveal significant increases in vegetation greenness and wetness, particularly in an Expanded Wetland subregion, which exhibited ∼3.5x higher wetness and ∼1.5x higher greenness trends compared to adjacent areas. Monthly metrics highlight seasonal variability, with increases in greenness linked to monsoonal rainfall and lateral water redistribution, indicating that restoration impacts extend beyond the primary wetland. This study demonstrates the utility of cloud-based platforms like Google Earth Engine and USGS EarthExplorer for long-term monitoring of wetland restoration, while quantifying the efficacy of the ‘Zeedyk approach’ and demonstrating its potential as a scalable method to restore and conserve wetland meadows in other arid and semi-arid landscapes.

New Mexico↗

USGS2025

Geochemical and hydrological investigations of historical data collected at the Lee Acres Landfill and Giant Bloomfield Refinery, New Mexico, 1985–2020

The Lee Acres Landfill and Giant Bloomfield Refinery are adjacent properties near the City of Farmington, New Mexico, each having undergone monitoring and remediation related to historical site activities. At the landfill, site cleanup has included the installation of a capillary barrier over former liquid waste lagoons and periodic monitoring of groundwater elevations and groundwater quality. At the refinery, remediation has focused on several petrochemical and crude oil release areas and included soil excavation, groundwater treatment, and regular monitoring of groundwater elevations and quality. Groundwater at both sites has higher concentrations of volatile organic compounds and trace metals than background aquifer concentrations. In 2022, the U.S. Geological Survey compiled the Lee Acres-Giant Bloomfield Refinery Database (LAGBRD), which contains publicly available groundwater-elevation data and organic and inorganic groundwater-quality data from both sites, spanning from 1985 to 2020. Data from the LAGBRD and precipitation data from other sources were used to better understand the cause of relatively high manganese concentrations observed in some groundwater wells at the site through comparison of groundwater chemistry to chemical end members, interpretation of spatial and temporal patterns in the groundwater chemistry, and interpretation of groundwater flow properties. In this study, elevated chloride concentrations in groundwater downgradient from the landfill have been attributed to landfill leachate based on the temporal and spatial variability of chloride concentrations and chloride-to-bromide ratios. Installation of a capillary barrier and surface-water runoff controls at the landfill in 2005 appears to have altered infiltration patterns at that site, resulting in a decrease in chloride at some wells but an increase in chloride and dissolved manganese at others. The timing and relation among groundwater elevation, chloride concentration, and manganese concentration suggest that leachate stored in the vadose zone provides a continued source of contamination to groundwater.

New Mexico↗

USGS2025

Depth of magma crystallization and fluid exsolution beneath the porphyry-skarn Cu deposits at Santa Rita and Hanover-Fierro, New Mexico, USA

The depth level at which porphyry Cu–forming magmas fractionated and exsolved mineralizing fluids is actively debated. In the classic model, extensive magma fractionation occurs in large, upper crustal magma chambers, and concomitant fluid exsolution leads to forceful expulsion of residual magmas in the form of porphyry dikes, stocks, and breccia pipes, which subsequently serve as pathways for the mineralizing fluids. In contrast, some recent studies highlighting the role of deep crustal magma fractionation in the production of fertile magmas essentially deny the existence of upper crustal magma chambers at the time of mineralization. To address this, we conducted a detailed thermobarometric investigation of 13 intermediate to felsic, porphyritic intrusive rocks related to porphyry-skarn Cu mineralization at Santa Rita and Hanover-Fierro, New Mexico, United States, representing two premineralization magmas (61–60 Ma), seven synmineralization magmas (60–58 Ma), and four late- to postmineralization magmas (58–57 Ma). For each sample, the pressure of last magma crystallization before final magma ascent to the current exposure level was reconstructed based on Al-in-hornblende barometry of small hornblende inclusions trapped within quartz phenocrysts and through titanium-in-quartz (TitaniQ) thermobarometry of the host quartz phenocrysts themselves. Since quartz is one of the last crystallizing magmatic minerals, and no significant phenocryst growth could have occurred in small dikes and stocks after final magma emplacement, quartz phenocrysts and their contained hornblende inclusions record the depth of last magma crystallization before final magma ascent. When present, hornblende phenocrysts and hornblende inclusions within other major phenocrysts were also analyzed. Both quartz and hornblende barometers return consistent average pressures of 3.2 ± 0.4 kbar for the entire suite of pre- to postmineralization magmas, corresponding to depths of 11 to 14 km. The synmineralization magmas return even more consistent average pressures of 3.1 ± 0.2 kbar, corresponding to a depth of 12 ± 1 km. The volume of the mineralizing porphyry dikes and stocks at the emplacement level is far too small to have provided all the fluids and metals required to form the observed ore deposits. Therefore, the majority of the ore-forming fluids must have originated from the magmas that crystallized at 12 ± 1 km depth. The ore deposits, conversely, formed at ~5-km paleodepth. This implies that most of the mineralizing fluids traveled an average vertical distance of ~7 km from their magmatic source to the eventual site of ore precipitation. The relatively unaltered nature and low veining degree of deeper parts of mineralized porphyry dikes and stocks suggest that the fluid transport through these intrusive bodies occurred mostly at near-solidus conditions by means of fluid percolation along grain boundaries. In summary, our results suggest that (1) a large, upper crustal pluton exists ~7 km beneath the Santa Rita and Hanover-Fierro deposits; (2) abundant phenocryst crystallization occurred at this depth level; and (3) this pluton was the main source for the exsolution of ore-forming fluids. However, the investigated rocks have elevated whole-rock Sr/Y ratios, indicating magma fractionation at deep crustal levels. As a result, our preferred model is a combination of the two end-member models introduced above, with most magma fractionation having occurred in the deep crust and with residual, intermediate to felsic melts having ascended and accumulated at 11 to 14 km paleodepth, where they continued to crystallize with comparatively little crystal-liquid separation, before some of these magmas ascended further to shallow levels and quenched to porphyries.

New Mexico↗