Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “fault localization”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Machine Learning to Identify Geologic Factors Associated with Production in Geothermal Fields: A Case-Study Using 3D Geologic Data from Brady Geothermal Field and NMFk

In this paper, we present an analysis using unsupervised machine learning (ML) to identify the key geologic factors that contribute to the geothermal production in Brady geothermal field. Brady is a hydrothermal system in northwestern Nevada that supports both electricity production and direct use of hydrothermal fluids. Transmissive fuid-fow pathways are relatively rare in the subsurface, but are critical components of hydrothermal systems like Brady and many other types of fuid-fow systems in fractured rock. Here, we analyze geologic data with ML methods to unravel the local geologic controls on these pathways. The ML method, non-negative matrix factorization with k-means clustering (NMFk), is applied to a library of 14 3D geologic characteristics hypothesized to control hydrothermal circulation in the Brady geothermal field. Our results indicate that macro-scale faults and a local step-over in the fault system preferentially occur along production wells when compared to injection wells and non-productive wells. We infer that these are the key geologic characteristics that control the through-going hydrothermal transmission pathways at Brady. Our results demonstrate: (1) the specific geologic controls on the Brady hydrothermal system and (2) the efficacy of pairing ML techniques with 3D geologic characterization to enhance the understanding of subsurface processes. This submission includes the published journal article detailing this work, the published 3D geologic map of the Brady Geothermal Area used as a basis to develop structural and geological variables that are hypothesized to control or effect permeability or connectivity, 3D well data, along which geologic data were sampled for PCA analyses, and associated metadata file. This work was done using the GeoThermalCloud framework, which is part of SmartTensors (both are linked below).

15 GEOTHERMAL ENERGY↗
NASA NTRS1973

Application of ERTS and EREP images to geologic investigations of the basin and range: Colorado Plateau boundary in northwestern and north central Arizona

The author has identified the following significant results. Regional mapping using ERTS-1 images has led to the recognition of two parallel northeast-trending systems of normal faults, the Bright Angel and Mesa Butte fault systems, each of which can be traced at least 100 miles. Many eruptive centers appear to be localized along these fault systems or along their extensions. The faults are chiefly observed in Phanerozoic rocks and have minor displacement, but are interpreted to reflect fault zones of major displacement in the crystalline Precambrian basement. The location of epicenters of recent earthquakes and reports of earthquakes by residents in the region suggest that the Bright Angel and Mesa Butte fault systems are currently active.

Goetz, A. F. H.↗
DOE OSTI2022

Creep fronts and complexity in laboratory earthquake sequences illuminate delayed earthquake triggering

Earthquakes occur in clusters or sequences that arise from complex triggering mechanisms, but direct measurement of the slow subsurface slip responsible for delayed triggering is rarely possible. We investigate the origins of complexity and its relationship to heterogeneity using an experimental fault with two dominant seismic asperities. The fault is composed of quartz powder, a material common to natural faults, sandwiched between 760 mm long polymer blocks that deform the way 10 meters of rock would behave. We observe periodic repeating earthquakes that transition into aperiodic and complex sequences of fast and slow events. Neighboring earthquakes communicate via migrating slow slip, which resembles creep fronts observed in numerical simulations and on tectonic faults. Utilizing both local stress measurements and numerical simulations, we observe that the speed and strength of creep fronts are highly sensitive to fault stress levels left behind by previous earthquakes, and may serve as on-fault stress meters.

58 GEOSCIENCES↗
DOE OSTI2019

Deciphering the Stress State of Seismogenic Faults in Oklahoma and Southern Kansas Based on an Improved Stress Map

Fault location and geometry are prime considerations in the reactivation of preexisting faults. Here, we combine relocated earthquake catalogs and focal mechanisms to delineate seismogenic faults in Oklahoma and southern Kansas and analyze their stress state. Initially, we identify and map seismogenic faults based on earthquake clustering. We then obtain an improved stress map using 2,047 high-quality focal mechanisms. The regional stress map shows a gradual transition from oblique normal faulting in western Oklahoma to strike-slip faulting in central and eastern Oklahoma. Stress amplitude ratio shows a strong correlation with pore pressure from hydrogeologic models, suggesting that pore pressure exhibits a measurable influence on stress patterns. Finally, we assess fault stress state via 3-D Mohr circles; a parameter understress is used to quantify the level of fault criticality (with 0 meaning critically stressed faults and 1 meaning faults with no applied shear stress). Our findings suggest that most active faults have near vertical planes (planarity >0.8 and dip >70°), and there is a strong correlation between fault length and maximum magnitude on each fault. The fault trends show prominent conjugate sets that strike [55–75°] and [105–125°]. A comparison with mapped sedimentary faults and basement fractures reveals common tectonic control. Based on 3-D Mohr circles, we find that 78% of the faults are critically stressed (understress ≤0.2), while several seismogenic faults are misoriented with high understress (>0.4). Fault geometry and local stress fields may be used to evaluate potential seismic hazard, as the largest earthquakes tend to occur on long, critically stressed faults.

58 GEOSCIENCES↗
DOE OSTI2019

Probabilistic Seismic Hazard Assessment for Georgia

Probabilistic Seismic Hazard Assessments (PSHA) form the basis for most contemporary seismic provisions in building codes around the world. The current building code of Georgia (2009) is not fully compatible with seismic hazard and design concepts of modern codes. The Ministry of Economy and Sustainable Development of Georgia is currently undertaking an update to the building code to align it with EuroCode. In conjunction with EuroCode, a National Annex specific to Georgia is needed, including seismic hazard values for the territory of Georgia. A comprehensive earthquake catalogue is fundamental to a robust PSHA in the region. The Caucasus has had an uneven seismic recording and reporting history both temporally and spatially. Therefore, there has been a strong regional effort to bring together data from all the countries in the Caucasus, to digitize and standardize paper-based data from the Soviet era, and to re-examine magnitude and location calculations. The earthquake catalogue resulting from this effort forms the basis for a reliable seismic hazard assessment. Recognizing these needs, the U.S. Department of Energy through Lawrence Livermore National Laboratory (LLNL) initiated a project in 2015 to engage and train local scientists in the Caucasus to install new equipment, to enhance the quality of seismic monitoring and reporting, and to improve and harmonize the regional earthquake catalogue. In 2016, the project scope was expanded to include the modernization of seismic hazard assessments. To this end, a series of workshops were held in Georgia, the first one of which presented an overview of PSHA and its components to various participants from the universities and other research organizations in the region. This kick-off workshop was followed by numerous hands-on training workshops in 2017 and 2018, covering all aspects of PSHA. In addition to their training goals, the workshops brought together existing and new information relevant to PSHA studies in the region, providing a platform for local scientists/engineers and experts from other parts of the world to interact, cooperate and discuss various pieces of the PSHA input. This report outlines the aspects of the earthquake catalogue development and PSHA project that were specifically undertaken to form the basis for the building code update in Georgia. Georgia straddles the Greater Caucasus and Lesser Caucasus Fold and Thrust belts, where the convergent tectonic setting generates significant earthquake activity. To the east of the country, evidence indicates the possibility of a remnant subducting slab that causes earthquakes deeper than 40km. This study systematically compiles and improves all available data on local seismicity, active faults, and ground motion attenuation characteristics of the region; and builds a framework to enable a contemporary PSHA to be carried out with the engagement of local scientists. To this end, six regional workshops were organized in total over a three-year period of time, with participation from the United States, Canada and Hungary as well as countries in the Caucasus: Georgia, Armenia, Azerbaijan and Turkey. The workshops served multiple purposes such as training, data collection, interactions between local earth scientists and engineers, and brainstorming and knowledge exchange among local and international experts.

Geosciences↗
DOE OSTI2024

Central Asia Seismic Hazard Assessment (CASHA): A Probabilistic Seismic Hazard Assessment for Kazakhstan, Kyrgyzstan and Tajikistan

Probabilistic seismic hazard assessments (PSHA) underpin the calculation of earthquake loads in most building codes around the world. In Central Asia, the building codes are slowly being updated to incorporate some of the contemporary concepts of seismic hazard representation. There is also a regional desire to coordinate hazard assessments and building code modernization. However, some challenges remain. Expertise in the region related to seismic hazard assessments is still largely compartmentalised, requiring a significant amount of training and capacity building in seismic hazard assessment related topics. In addition, there are vast amounts of seismic data (bulletin and waveforms), both from analogue and digital eras, that the region’s countries stored but until recently did not use or share among themselves or with the broader seismological community around the world. Finally, after the collapse of the Soviet Union in the 1990s, many of the countries’ seismic networks suffered a major setback with the lack of attention and budget to update existing equipment and installation of new instruments. In order to address these issues, the United States Department of Energy through Lawrence Livermore National Laboratory (LLNL) initiated a project in 2016 to engage and train local scientists in Central Asia to install new equipment, to enhance the quality of seismic monitoring and reporting, to improve and harmonise the regional earthquake catalogue, and to conduct national probabilistic seismic hazard assessments using the new and improved datasets. To achieve the seismic hazard assessment related goals, a series of workshops were held in Almaty, Kazakhstan; Bishkek, Kyrgyzstan; and Dushanbe, Tajikistan from 2016 until 2020. During the time that the COVID-19 pandemic restricted travel, workshops continued online (22 online workshops were hosted in two years). Finally, in May 2022, an in-person workshop in Istanbul, Turkey brought together all project participants along with civil engineers engaged with building code activities in their respective countries, providing a platform to discuss the implementation of the hazard models into updates of building codes in each country, as well as to discuss model parameters, sensitivity analyses and model results in terms of hazard maps, uniform hazard spectra and hazard deaggregations. The workshops were a combination of lectures and hands-on exercises, and included international participation as well as local scientists and engineers. The workshops served several purposes, including training, coordination of data collection, interactions between local earth scientists and engineers, and brainstorming and knowledge exchange among local and international experts. This report outlines the new earthquake catalogue compilation effort and the PSHA project undertaken in Kyrgyzstan, Tajikistan, and Kazakhstan as part of this initiative. The southern part of this region is tectonically active with moderate to high levels of both shallow crustal seismic activity and occurrence of deeper earthquakes under the Hindu Kush and Pamir mountain ranges. Deeper earthquakes also occur near southwestern Kazakhstan, under the eastern Greater Caucasus and Caspian Sea. Large portions of central and northern Kazakhstan, on the other hand, are in stable continental regions with low levels of seismic activity. This study systematically compiles and improves all available data on local seismicity, active faults, and ground motion attenuation characteristics of the region; and builds a framework to enable a contemporary PSHA to be carried out with the engagement of local scientists. While the project was regional, the seismic hazard assessments are primarily driven by the countries’ own national preferences and understanding of data collection, interpretation, and validation of results.

58 GEOSCIENCES↗
NASA NTRS1987

Analysis of typical fault-tolerant architectures using HARP

Difficulties encountered in the modeling of fault-tolerant systems are discussed. The Hybrid Automated Reliability Predictor (HARP) approach to modeling fault-tolerant systems is described. The HARP is written in FORTRAN, consists of nearly 30,000 lines of codes and comments, and is based on behavioral decomposition. Using the behavioral decomposition, the dependability model is divided into fault-occurrence/repair and fault/error-handling models; the characteristics and combining of these two models are examined. Examples in which the HARP is applied to the modeling of some typical fault-tolerant systems, including a local-area network, two fault-tolerant computer systems, and a flight control system, are presented.

Bavuso, Salvatore J.↗
NASA NTRS1991

Models of recurrent strike-slip earthquake cycles and the state of crustal stress

Numerical models of the strike-slip earthquake cycle, assuming a viscoelastic asthenosphere coupling model, are examined. The time-dependent simulations incorporate a stress-driven fault, which leads to tectonic stress fields and earthquake recurrence histories that are mutually consistent. Single-fault simulations with constant far-field plate motion lead to a nearly periodic earthquake cycle and a distinctive spatial distribution of crustal shear stress. The predicted stress distribution includes a local minimum in stress at depths less than typical seismogenic depths. The width of this stress 'trough' depends on the magnitude of crustal stress relative to asthenospheric drag stresses. The models further predict a local near-fault stress maximum at greater depths, sustained by the cyclic transfer of strain from the elastic crust to the ductile asthenosphere. Models incorporating both low-stress and high-stress fault strength assumptions are examined, under Newtonian and non-Newtonian rheology assumptions. Model results suggest a preference for low-stress (a shear stress level of about 10 MPa) fault models, in agreement with previous estimates based on heat flow measurements and other stress indicators.

Lyzenga, Gregory A.↗
DOE OSTI2016

Protected gates for topological quantum field theories

We study restrictions on locality-preserving unitary logical gates for topological quantum codes in two spatial dimensions. A locality-preserving operation is one which maps local operators to local operators — for example, a constant-depth quantum circuit of geometrically local gates, or evolution for a constant time governed by a geometrically local bounded-strength Hamiltonian. Locality-preserving logical gates of topological codes are intrinsically fault tolerant because spatially localized errors remain localized, and hence sufficiently dilute errors remain correctable. By invoking general properties of two-dimensional topological field theories, we find that the locality-preserving logical gates are severely limited for codes which admit non-abelian anyons, in particular, there are no locality-preserving logical gates on the torus or the sphere with M punctures if the braiding of anyons is computationally universal. Furthermore, for Ising anyons on the M-punctured sphere, locality-preserving gates must be elements of the logical Pauli group. We derive these results by relating logical gates of a topological code to automorphisms of the Verlinde algebra of the corresponding anyon model, and by requiring the logical gates to be compatible with basis changes in the logical Hilbert space arising from local F-moves and the mapping class group.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

New Insights from Legacy Seismic Data regarding Basalt Elevations and Variability on the Hanford Site

Migration of groundwater contaminants in the Gable Gap area of the Hanford Site in southeastern Washington State is strongly influenced by the distribution and permeability of basalts that lie beneath an unconfined aquifer. Locally, folding and faulting of the Columbia River Basalt associated with the Yakima fold and thrust belt followed by erosion due to the Lake Missoula floods resulted in a complex basalt surface that represents either an impermeable lower boundary to the unconfined aquifer system or localized regions of increased permeability that potentially promote communication between the unconfined aquifer system and deeper, confined aquifer systems. Paleo-channels carved into the basalt by floodwaters are thought to provide preferential flow paths for groundwater contaminants. In 2011, a seismic landstreamer campaign was carried out to image the basalt surface and produced pre-stack depth migrated p-wave reflection images. The reflection images identified two large troughs that may represent paleo-channels and several areas of possible faulting. Here, the streamer data are re-analyzed using refraction travel-time and Rayleigh wave dispersion analyses to obtain images of compressional and shear wave velocities within the suprabasalt sediment sections and the upper basalt surface. The combined interpretation of reflection and seismic velocity images shows complexity in the basalt velocity and elevation, which varies by 50 m or more within the study area. These results, along with other ongoing geophysical investigations, will be used to inform the site geologic model and potentially guide placement of future boreholes needed to quantify vertical flow between the confined and unconfined aquifers.

58 GEOSCIENCES↗
DOE OSTI2016

Ab initio study of point defects near stacking faults in 3C-SiC

Interactions between point defects and stacking faults in 3C-SiC are studied using an ab initio method based on density functional theory. The results show that the discontinuity of the stacking sequence considerably affects the configurations and behavior of intrinsic defects, especially in the case of silicon interstitials. The existence of an intrinsic stacking fault (missing a C-Si bilayer) shortens the distance between the tetrahedral-center site and its second-nearest-neighboring silicon layer, making the tetrahedral silicon interstitial unstable. Instead of a tetrahedral configuration with four C neighbors, a pyramid-like interstitial structure with a defect state within the band gap becomes a stable configuration. In addition, orientation rotation occurs in the split interstitials that has diverse effects on the energy landscape of silicon and carbon split interstitials in the stacking fault region. Moreover, our analyses of ionic relaxation and electronic structure of vacancies show that the built-in strain field, owing to the existence of the stacking fault, makes the local environment around vacancies more complex than that in the bulk.

36 MATERIALS SCIENCE↗
DOE OSTI2018

Pressure transient analysis during CO 2 push-pull tests into faults for EGS characterization

Here, with the goal of detecting and characterizing faults and fractures in enhanced geothermal systems (EGS), a new technology involving CO 2 push-pull testing, active-source geophysical imaging, and well logging has recently been proposed. This technique takes advantage of (1) the contrasting properties of supercritical CO 2 and water which cause CO 2 to appear distinct from surrounding brine in seismic and other geophysical logging approaches, (2) the non-wetting nature of CO 2 which keeps it localized to the faults and fractures to create contrast potentially sufficient for active seismic and well-logging approaches to image faults and fracture zones at EGS sites. In this study, we evaluate the feasibility of using pressure transient monitoring during CO 2 push-pull tests to complement active seismic and wireline well logging for EGS characterization. For this purpose, we developed a 2D model of a prototypical geothermal site (Desert Peak, NV) that includes a single fault. The fault zone consists of a slip plane, fault gouge, and damage zone, and is bounded by the surrounding matrix of the country rock. Through numerical simulation using iTOUGH2, we found that the pressure transient at the monitoring wells in the fault gouge shows unique traits due to the multiphase flow conditions developed by CO 2 injection, and varies sensitively on the arrival of the CO 2 plume and the degree of CO 2 saturation. A sensitivity analysis shows the pressure transient is most sensitive to the fault gouge permeability, but also depends on multiphase flow parameters and the boundary conditions of the fault. Lastly, an inversion study reveals that the fault gouge permeability can be best estimated with the pressure transient data, whereas additional CO 2 saturation data do not improve the accuracy of the inversion significantly.

58 GEOSCIENCES↗
DOE OSTI2024

IoT-Enabled Traveling Wave Microgrid Protection

Traveling Wave Protection based on the Internet of Things (TWP-IoT) is developed to enable ultra-resilient microgrids. TWP-IoT adopts a directional traveling wave approach that utilizes the voltage and current modal components obtained from Discrete Hilbert Transform to identify the fault features via ultra-lightweight IoT hardware. Furthermore, the main contributions of this work include (1) a novel TWP-IoT for microgrid and distribution networks that can achieve the time-of-arrival fault detection and localization with composite wave impedance and conductance using Discrete Hilbert Transform, (2) a cost-effective smart IoT prototype with functionalities outperforming commercial relays, and (3) a real-time hardware-in-the-loop prototype based on RTDS that verifies the robustness and efficacy of TWP-IoT subject to a diverse set of working conditions. TWP-IoT is found to have excellent performance under a wide spectrum of internal and external fault conditions across a myriad of microgrid locations, which is unattainable by today’s TWP products.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2025

A Novel Protection Scheme for Unbalanced Faults in Inverter Dominated Networks: A Computationally Efficient Algorithm for Entry-Level Relays

Microgrids are now a common practice in distribution systems to increase resilience and reliability. However, microgrid protection remains a critical challenge, considering its requirement to operate in both grid connected and islanded, and the variability in fault characteristics under each mode of operation. This paper presents unbalanced power (S unb ) based fault detection algorithm, which considers local voltage and current unbalances to determine faults in the system. S unb is a computationally efficient fault detection algorithm that is suitable for implementation in the programmable logic of entry level protective relays. In addition, the difference in current and voltage unbalance (D n ) is used to determine the fault type. The proposed method demonstrates high sensitivity and selectivity for line-to-ground (LG), line-to-line (LL), and double line-to-ground (LLG) faults, representing the most common faults in distribution systems. It also allows relay coordination with upstream and downstream protection devices in both island and grid connected operation, while preserving grading margins. The same pickup and time multiplier settings of a particular relay for both modes of operation eliminates the need for adaptive settings, which rely on communication networks. Validation was performed with a hardware-in-the-loop (HIL) setup using Typhoon HIL real time simulator interfaced with three entry-level, SEL 751 relays. Results confirmed the algorithm’s ability to discriminate fault conditions, and determine the fault type under both operating modes, maintain fast detection times, and ensure proper protection coordination.

fault classification↗
DOE OSTI2020

The nature and origins of decametre-scale porosity in Ordovician carbonate rocks, Halahatang oilfield, Tarim Basin, China

At >7 km depths in the Tarim Basin, hydrocarbon reservoirs in Ordovician rocks of the Yijianfang Formation contain large cavities ( c . 10 m or more), vugs, fractures and porous fault rocks. Although some Yijianfang Formation outcrops contain shallow (formed near surface) palaeokarst features, cores from the Halahatang oilfield lack penetrative palaeokarst evidence. Outcrop palaeokarst cavities and opening-mode fractures are mostly mineral filled but some show evidence of secondary dissolution and fault rocks are locally highly ( c . 30%) porous. Cores contain textural evidence of repeated formation of dissolution cavities and subsequent filling by cement. Calcite isotopic analyses indicate depths between c . 220 and 2000 m. Correlation of core and image logs shows abundant cement-filled vugs associated with decametre-scale fractured zones with open cavities that host hydrocarbons. A Sm–Nd isochron age of 400 ± 37 Ma for fracture-filling fluorite indicates that cavities in core formed and were partially cemented prior to the Carboniferous, predating Permian oil emplacement. Repeated creation and filling of vugs, timing constraints and the association of vugs with large cavities suggest dissolution related to fractures and faults. In the current high-strain-rate regime, corroborated by velocity gradient tensor analysis of global positioning system (GPS) data, rapid horizontal extension could promote connection of porous and/or solution-enlarged fault rock, fractures and cavities. Supplementary material: Stable isotopic analyses and the velocity gradient tensor and principal direction and magnitude calculation are available at https://doi.org/10.6084/m9.figshare.c.4946046 Thematic collection: This article is part of the The Geology of Fractured Reservoirs collection available at: https://www.lyellcollection.org/cc/the-geology-of-fractured-reservoirs

Geology↗
DOE OSTI2025

Stratigraphy‐Induced Localization of Microseismicity During CO 2 Injection in Illinois Basin

Abstract Subsurface fluid injection stimulates complex hydromechanical interaction, necessitating the integration of geomechanical data across spatial and temporal scales to consider the sophisticated behavior. Induced seismic response is usually associated with the complex reservoir architecture and pre‐existing features that are three‐dimensional, such as local stratigraphy, fractures, faults, and other discontinuities. This study encompasses laboratory characterization of the coupled hydromechanical response of cores extracted from rock formations in Illinois Basin: reservoir ‐ Mt. Simon sandstone, basal seal ‐ Argenta sandstone, and crystalline basement ‐ Precambrian rhyolite. High‐resolution numerical modeling allows considering the three‐dimensional complexity of the Illinois Basin Decatur Project with spatial resolution comparable to one of the active seismic surveys. A detailed reconstruction of the evolving state of stress in formations lacking direct stress measurements is achieved by numerical modeling that integrated laboratory‐derived hydromechanical properties, a porosity‐permeability relationship, active seismic data, and an inverted three‐dimensional porosity distribution. It appears that the microseismic clusters, mainly observed in the crystalline basement during the injection, are linked to zones experiencing more critically stressed conditions prior to injection. These zones have a potential for reactivation during the injection and are attributed to the specific local stratigraphy of the injection site, as well as transfer of triggering perturbations during the injection.

Bondarenko, N. [University of Illinois Urbana‐Cham↗
DOE OSTI2024

Identifying recharge sources and their impacts on a North Central New Mexico shallow aquifer using unsupervised machine learning

In this article, shallow aquifers are important but highly variable resources in arid to semi-arid regions. Limited shallow aquifer volume results in high sensitivity to recharge fluctuations, which can impact the local fauna and flora, and transport of contaminants in the aquifer or vadose zone. Aquifer response to external forcing (e.g., precipitation) is usually solved by estimating aquifer parameters and running physics-based models to match known fluctuations of hydraulic head. However, this technique is time and computationally expensive. Furthermore, high aquifer complexity decreases precision in physics-based models. Alternatively supervised machine learning is used to predict aquifer dynamics. However, these techniques rely on input data and struggle to interpret aquifer response for missing sources (i.e., snowpack data). To counter these problems, we propose an unsupervised machine learning technique (NMFk) to estimate the impact of different sources on aquifer recharge. NMFk is used to understand the influence of external forcing on shallow aquifer recharge in the Pajarito Plateau (Los Alamos, NM, USA). The results show how NMFk can be used to reduce the data dimension in a complex field dataset to three recharge signals that cause fluctuations within the field data. Here, the source signals are interpreted as rainfall, snowmelt, and a delayed aquifer response to the previous two signals. These results evidence how heterogeneous aquifers delimited by canyons incised into the Pajarito Plateau respond in similar ways to the source signals identified by NMFk. Furthermore, results show the importance of the local geology where faults act as sinks, and anthropogenic disturbances can facilitate infiltration amplifying the interpreted signal.

54 ENVIRONMENTAL SCIENCES↗