SEARCH · Engineering Papers

Results for “two level system loss”

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.

SourceResearch connectionFrom yearThrough year

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

At least 109 records · Page 6

All sources NASA NTRS DOE OSTI

DOE OSTI2020

Experimental evolution reveals nitrate tolerance mechanisms in Desulfovibrio vulgaris

Elevated nitrate in the environment inhibits sulfate reduction by important microorganisms of sulfate-reducing bacteria (SRB). Several SRB may respire nitrate to survive under elevated nitrate, but how SRB that lack nitrate reductase survive to elevated nitrate remains elusive. To understand nitrate adaptation mechanisms, we evolved 12 populations of a model SRB (i.e., Desulfovibrio vulgaris Hildenborough, DvH) under elevated NaNO 3 for 1000 generations, analyzed growth and acquired mutations, and linked their genotypes with phenotypes. Nitrate-evolved (EN) populations significantly (p < 0.05) increased nitrate tolerance, and whole-genome resequencing identified 119 new mutations in 44 genes of 12 EN populations, among which six functional gene groups were discovered with high mutation frequencies at the population level. In this work, we observed a high frequency of nonsense or frameshift mutations in nitrosative stress response genes (NSR: DVU2543, DVU2547, and DVU2548), nitrogen regulatory protein C family genes (NRC: DVU2394-2396, DVU2402, and DVU2405), and nitrate cluster (DVU0246-0249 and DVU0251). Mutagenesis analysis confirmed that loss-of-functions of NRC and NSR increased nitrate tolerance. Also, functional gene groups involved in fatty acid synthesis, iron regulation, and two-component system (LytR/LytS) known to be responsive to multiple stresses, had a high frequency of missense mutations. Additionally, mutations in those gene groups could increase nitrate tolerance through regulating energy metabolism, barring entry of nitrate into cells, altering cell membrane characteristics, or conferring growth advantages at the stationary phase. This study advances our understanding of nitrate tolerance mechanisms and has important implications for linking genotypes with phenotypes in DvH.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2022

Construction of bias-preserving operations for pair-cat codes

Fault-tolerant quantum computation with depolarization error often requires demanding error threshold and resource overhead. If the operations can maintain high noise bias—dominated by dephasing error with small bit-flip error—we can achieve hardware-efficient fault-tolerant quantum computation with a more favorable error threshold. Distinct from two-level physical systems, multilevel systems (such as harmonic oscillators) can achieve a desirable set of bias-preserving quantum operations while using continuous engineered dissipation or Hamiltonian protection to stabilize to the encoding subspace. For example, cat codes stabilized with driven-dissipation or Kerr nonlinearity can possess a set of bias-preserving gates while continuously correcting bosonic dephasing error. However, cat codes are not compatible with continuous quantum error correction against excitation loss error, because it is challenging to continuously monitor the parity to correct photon loss errors. In this work, we generalize the bias-preserving operations to pair-cat codes, which can be regarded as a multimode generalization of cat codes, to be compatible with continuous quantum error correction against both bosonic loss and dephasing errors. In conclusion, our results open the door towards hardware-efficient robust quantum information processing with both bias-preserving operations and continuous quantum error correction simultaneously correcting bosonic loss and dephasing errors.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2022

Semi-Annual Report for Modular Integrated Gas High Temperature Reactor Development during Performance Period October 2021 - March 2022

Modular Integrated Gas High Temperature Reactor (MIGHTR) is being designed by a multidisciplinary team of nuclear, mechanical, and structural engineers under the support of a DOE-NE Advanced Reactor Demonstration Program’s Advanced Reactor Concepts-20 (ARC-20) award. The objective of this ARC-20 project is to deliver a conceptual design for the proposed MIGHTR in 3 years and support its commercialization as a safe and low-cost HTGR. Argonne National Laboratory is responsible for the design and analysis of the reactor cavity cooling system (RCCS) as a safety system for passive decay heat removal of the reactor concept. Work performed by Argonne also includes the thermal hydraulics analysis of the primary heat transport system including reactor pressure vessel (RPV) internals for normal operation, shutdown, and accident conditions. This first semi-annual report summarizes the progress made at Argonne on the two tasks since the start of the project. For the RCCS design task, high-level design work was performed including identification of design requirements for the conceptual MIGHTR RCCS and initial design calculations. The baseline dimensions of the MIGHTR RCCS were derived from scoping analysis results and preliminary steady state performance was estimated. At the initial stage of the primary system analysis task, the primary system thermal fluids model for the MIGHTR core was outlined to use a simplified two-dimensional model to predict the temperature distribution and coolant pressure losses in the fuel assembly blocks.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

DOE OSTI2023

A Planetary Defense System For Responsive Space Missions

As King et al succinctly write, “Collision between the Earth and small objects in the solar system are known to occur regularly, and can result in extinction-level events at their most energetic or more commonly widespread loss of life or economic destruction even for relatively low-energy impacts.” Currently, the JPL Sentry System and the NASA NEO Surveyor mission track, identify, and monitor potential asteroid threats. The JPL Sentry System4 is an automated collision monitoring system which scans the current asteroid catalog for potential impact trajectories with Earth. A 2017 study led by P. Tricarico titled “The near-Earth asteroid population from two decades of observations” analyzes the combined observations of nine leading asteroid surveys. The NASA NEO Surveyor mission seeks to identify 90% of near-Earth objects (NEO) with diameters greater than 140 meters (460 feet) within 10 years. The NEO Surveyor combines data from two infrared imaging detectors to measure NEO size, composition, shape, rotational state, and orbits. Recently, the NASA Double Asteroid Redirection Test (DART) mission successfully altered the orbit of asteroid moon Dimorphos in 2022. This demonstrated that orbit alteration is possible. Thus, collisions between the Earth and asteroids are both likely and demonstrate potential for great harm. There is ongoing work to identify, catalog, and track potential threats. However, when a threat is identified there is not yet a system in place to target the threatening asteroid and defend the planet from potentially existential threats. Building on the DART Mission, we explore a planetary defense system to mitigate asteroid threats.

42 ENGINEERING↗

DOE OSTI2021

Tri-Level Scheduling Model Considering Residential Demand Flexibility of Aggregated HVACs and EVs Under Distribution LMP

Residential loads, especially heating, ventilation and air conditioners (HVACs) and electric vehicles (EVs), have great potentials to provide demand flexibility which is an attribute of grid-interactive efficient buildings (GEB). Under this new paradigm, first, EV and HVAC aggregator models are developed in this paper to represent the fleet of GEBs, in which the aggregated parameters are obtained based on a new approach of data generation and least squares parameter estimation (DG-LSPE), which can deal with heterogeneous HVACs. Then, a tri-level bidding and dispatching framework is established based on competitive distribution operation with distribution locational marginal price (DLMP). Furthermore, the first two levels form a bilevel model to optimize the aggregators’ payment and to represent the interdependency between load aggregators and the distribution system operator (DSO) using DLMP, and the third level is to dispatch the optimal load aggregation to all residents by the proposed priority list-based demand dispatching algorithm. Finally, case studies on a modified IEEE 33-Bus system illustrate three main technical reasons of payment reduction due to demand flexibility: load shift, DLMP step changes, and power losses. They can be used as general guidelines for better decision-making for future planning and operation of demand response programs.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2023

Performance Analyses of Supercritical Carbon Dioxide-Based Parabolic Trough Collectors with Double-Glazed Receivers

Supercritical carbon dioxide is becoming a hot research topic as a potential heat transfer fluid in parabolic trough concentrators since it enables operating the solar system at high temperatures for a higher quality of energy. However, the corresponding inflated thermal losses necessitate alternative receiver designs. This work examines four double-glazed receivers, with each annular space being evacuated or non-evacuated, in terms of the absorber tube's diameter (53-80 mm) and the diameter ratios of the two glass shells (1.2-2.0). An analytical model is developed and validated for this purpose, and the four designs are further examined using ground-level solar and meteorological measurements. The results demonstrate higher performance in the case of fully evacuating the receiver and using the smallest possible diameters of the three concentric cylinders, where the energy and exergy efficiencies reach 65.3 and 40.3%, respectively. Yet, evacuating only the inner annular space is sufficient to achieve virtually the same performance level. This energy efficiency decreases to 62% in case of increasing the tube diameter to 80 mm. As the operating temperature increases from 423 to 850 K, the specific thermal losses increase by 3.98-4.34 folds, depending on the receiver design. Double glazing the receivers is favorable only at high operating temperatures of sCO2, where the reduction in thermal losses overcomes the drop in optical efficiency. For an inlet sCO2 temperature of 850 K, thermal losses are reduced by 33.64 and 53.92%, compared to evacuated and non-evacuated single-glazed receivers, respectively. Throughout the year, the fully evacuated and fully non-evacuated double-glazed receivers have energy efficiencies of 54.96 and 52.39%, exergy efficiencies of 33.64 and 32.06%, and thermal losses of 348.8 and 402.5 W/m, respectively.

analytical model↗

DOE OSTI2026

Evolution of the Antarctic Ice Sheet from 2000–2300 and beyond: model sensitivity and uncertainty analysis using MPAS-Albany Land Ice

We present a description of the Antarctic Ice Sheet model configuration submitted to the ISMIP6-Antarctica-2300 experiment using the MPAS-Albany Land Ice model, along with three new sets of simulations: (1) a set of extended simulations to 2500 for three forced experiments and to 2775 for the control experiment; (2) a sensitivity analysis of our model configuration to parameters controlling basal sliding and sub-shelf melt, and to model structural choices including the choice of the energy and stress balances; and (3) a 72-member ensemble run on graphics processing units (GPUs) and analysis of variance to determine the primary sources of uncertainty in our ice-sheet model projections. Our extended simulations predict rapid retreat beginning after 2300 for SSP1-2.6 forcing and after 2500 for present-day (control) forcing, primarily in the Amundsen Sea Embayment. We find that varying the sub-shelf melt parameter between the 5th to 95th percentile values for a mean-Antarctic calibration target results in an up to ∼ ± 40 % change in sea-level contribution relative to our baseline simulations that used the median value. Using a linear basal sliding law reduces sea-level contribution by 51 %–73 % relative to our baseline nonlinear sliding law with an exponent of 1/5. When using basal sliding law exponents of 1/3 and 1/10, the overall difference from our baseline simulations at 2300 is on the order of 10 %. The Amundsen Sea Embayment region displays a strongly non-linear dependence of mass loss on the sliding law exponent, with no discernible relationship between the sliding law exponent and the mass loss by 2300, while the sectors feeding the Ross and Filchner-Ronne ice shelves exhibit more mass loss with a more-plastic sliding law. Our model fidelity sensitivity experiments reveal a 9 %–31 % increase in sea-level contribution when using a depth-integrated stress balance approximation relative to our three-dimensional solver, while using a fixed-in-time temperature field increases sea-level contribution by 14 %–88 % relative to two thermomechanically coupled configurations. Our 72-member ensemble and analysis of variance show that the uncertainty in long-term projections is dominated by the choice of Earth system model forcing and the presence or absence of hydrofracture forcing, rather than uncertainty in sliding and sub-shelf melt parameters.

58 GEOSCIENCES↗

DOE OSTI2020

Reconstruction and uncertainty quantification of lattice Hamiltonian model parameters from observations of microscopic degrees of freedom

The emergence of scanning probe and electron beam imaging techniques has allowed quantitative studies of atomic structure and minute details of electronic and vibrational structure on the level of individual atomic units. These microscopic descriptors, in turn, can be associated with local symmetry breaking phenomena, representing the stochastic manifestation of the underpinning generative physical model. In this work, we explore the reconstruction of exchange integrals in the Hamiltonian for a lattice model with two competing interactions from observations of microscopic degrees of freedom and establish the uncertainties and reliability of such analysis in a broad parameter-temperature space. In contrast to other approaches, we specifically specify a loss function inherent to thermodynamic systems and utilize it to estimate uncertainty in simulated realizations of different models. As an ancillary task, we develop a machine learning approach based on histogram clustering to predict phase diagrams efficiently using a reduced descriptor space. We further demonstrate that reconstruction is possible well above the phase transition and in the regions of parameter space when the macroscopic ground state of the system is poorly defined due to frustrated interactions. This suggests that this approach can be applied to the traditionally complex problems of condensed matter physics such as ferroelectric relaxors and morphotropic phase boundary systems, spin and cluster glasses, and quantum systems once the local descriptors linked to the relevant physical behaviors are known.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2025

Accelerating Embedding Potential Optimization by Reconstructing the Pseudo-Valence Electron Density

Density functional embedding theory (DFET) enables use of electronic structure methods with higher accuracy than density functional theory in a local region, with applications thus far ranging from (photo/electro)catalysis to reactions in solution. DFET partitions a large collection of atoms into smaller groups that interact via a shared embedding (interaction) potential V emb , determined via functional optimization. The optimized effective potential (OEP) process used to optimize V emb is time-consuming and becomes a computational bottleneck due to sharp, oscillating features of V emb near nuclei. Here, similar to pseudopotential theory, by reconstructing electron densities used in the OEP process from smoother pseudo-valence-only (PVO) electron densities as proxies for total densities of the full system and subsystems, we can retain accuracy in the embedded electronic structure calculations while potentially reducing the overhead of V emb construction, within the projector augmented-wave (PAW) formalism. We explore three different chemical reactions as exemplars to test PVO–DFET, namely, H 2 dissociative adsorption on a Cu(111) surface, H 2 O adsorption on a Pt(111) surface, and aqueous [Ca 2+ –SO 4 2– ] ion-pair formation. The PVO approximation works well for all three systems with minimal loss of accuracy (∼10–70 meV error relative to the original exact-derivative (ED) approach) while accelerating V emb generation for the Cu and Pt systems respectively by 20× and 5×. Given proper numerical convergence parameters, the spatial distributions of differences between PVO- and ED-based V emb outside the core regions are small, explaining the exceptional agreement between the two approaches. Finally, we anticipate that this more efficient PVO–DFET approximation will be useful whenever computation of V emb is much more expensive than subsequent embedded high-level electron correlation calculations.

approximation↗

DOE OSTI2023

Iterated Gauss-Seidel GMRES

The GMRES algorithm of Saad and Schultz [SIAM J. Sci. Stat. Comput., 7 (1986), pp. 856-869] is an iterative method for approximately solving linear systems Ax = b, with initial guess x0 and residual r0 = b Ax0. The algorithm employs the Arnoldi process to generate the Krylov basis vectors (the columns of Vk ). It is well known that this process can be viewed as a QR factorization of the matrix Bk = [r0, AVk] at each iteration. Despite an O (..epsilon..)..kappa.. (Bk ) loss of orthogonality, for unit roundoff ..epsilon..and condition number ..kappa.. , the modified Gram-Schmidt formulation was shown to be backward stable in the seminal paper by Paige et al. [SIAM J. Matrix Anal.Appl., 28 (2006), pp. 264-284]. We present an iterated Gauss-Seidel formulation of the GMRES algorithm (IGS-GMRES) based on the ideas of Ruhe [Linear Algebra Appl., 52 (1983), pp. 591-601] and Swirydowicz et al. [Numer. Linear Algebra Appl., 28 (2020), pp. 1-20]. IGS-GMRES maintains orthogonality to the level O (..epsilon..)..kappa.. (Bk ) or O (..epsilon..), depending on the choice of one or two iterations; for two Gauss-Seidel iterations, the computed Krylov basis vectors remain orthogonal to working accuracy and the smallest singular value of Vk remains close to one. The resulting GMRES method is thus backward stable. We show that IGS-GMRES can be implemented with only a single synchronization point per iteration, making it relevant to large-scale parallel computing environments. We also demonstrate that, unlike MGS-GMRES, in IGS-GMRES the relative Arnoldi residual corresponding to the computed approximate solution no longer stagnates above machine precision even for highly nonnormal systems.

Arnoldi-QR↗

DOE OSTI2026

Quasiparticle spectroscopy in tantalum films with different Ta/sapphire interfaces

One of the crucial aspects of current research in quantum information science is the identification and control of loss mechanisms in superconducting circuits. Although microwave measurements directly quantify device performance, additional techniques that probe quasiparticle excitations in superconducting films are needed to understand the microscopic mechanisms underlying dissipation and decoherence. Here, we present results from quasiparticle spectroscopy of Ta/sapphire films by measuring the Meissner-state magnetic susceptibility using a precision frequency-domain resonator specifically designed for thin films. We find direct evidence for additional low-energy excitations in samples with lower internal quality factors. These excitations are consistent with deep subgap states due to two-level systems, Yu-Shiba-Rusinov states near the gap edge, and perhaps other pair-breaking mechanisms. The developed non-destructive frequency-domain quasiparticle spectroscopy is a valuable addition to the quantum materials toolbox.

Moirangthem, Bicky Singh [Ames Lab; Iowa State U.]↗

DOE OSTI2026

Quasiparticle spectroscopy in tantalum films with different Ta/sapphire interfaces

Moirangthem, Bicky Singh [Ames Lab; Iowa State U.]↗

DOE OSTI2026

Machine learning based prediction of airflow maldistribution in air-to-refrigerant heat exchangers

Flow maldistribution is a common challenge in heat exchanger (HX) design and particularly important for air-to-refrigerant geometries where capacity losses can approach 65%. This has a major impact on central air conditioning systems, as compact duct design motivates the use of A-type HXs which are known to be affected by airflow maldistribution. Because velocity profiles are difficult to predict, components are often oversized leading to increased material cost, system footprint, and refrigerant charge. Several studies detail airflow maldistribution for individual HXs and packages, but findings cannot always be extrapolated to new designs. In this work, a machine learning (ML) based flow profile prediction framework is developed and applied to two common package configurations: (i) A-type and (ii) U-type HXs, across a broad range of HX geometries and flow rates. Porous media CFD simulations are validated against independent data for both package types as well as comprehensive in house measurements for a finless geometry with shape optimized non-round tubes, which validates the framework for new heat transfer surfaces. The ML models are trained on the porous media CFD simulations, predicting volumetric flow rate (VFR) within 1.1% and 1.9% with maximum relative L 2 norm errors of 0.48 and 0.65, respectively, while also delivering 10 5 speed up factor compared to full porous media CFD. HX level simulations show an up to 9% reduction in heat transfer from flow maldistribution, with greater losses occurring at smaller half apex angles. This framework enables rapid and highly accurate prediction of airflow maldistribution induced capacity degradation.

42 ENGINEERING↗

DOE OSTI2026

Quasiparticle spectroscopy in tantalum films with different Ta/sapphire interfaces

One of the crucial aspects of current research in quantum information science is the identification and control of loss mechanisms in superconducting (SC) circuits. Although microwave measurements directly quantify device performance, additional techniques that probe quasiparticle excitations in SC films are needed to understand the microscopic mechanisms underlying dissipation and decoherence. Here, we present results from quasiparticle spectroscopy of Ta/sapphire films by measuring the Meissner-state magnetic susceptibility using a precision frequency-domain resonator specifically designed for thin films. We find direct evidence for additional low-energy excitations in samples with lower internal quality factors. These excitations are consistent with deep subgap states due to two-level systems, Yu–Shiba–Rusinov states near the gap edge, and perhaps other pair-breaking mechanisms. The developed non-destructive frequency-domain quasiparticle spectroscopy is a valuable addition to the quantum materials toolbox.

metal/substrate interlayer↗

DOE OSTI2022

Highly Efficient Layer-by-Layer Processed Quaternary Organic Solar Cells with Improved Charge Transport and Reduced Energy Loss

We report how to improve charge transport properties and reduce energy losses are the two main challenges of organic solar cells (OSCs). Herein, a new device fabrication strategy that combines the layer-by-layer (LBL) method and a quaternary active layer is demonstrated to realize highly efficient photovoltaic performance of the OSCs based on PM6:PM7/Y6:O1-2F. The introduction of the second donor, PM7,forms a cascade energy-level alignment with PM6 that facilitates effective charge transfer. The second acceptor, O1-2F, which is poorly miscible with the donors, helps to reduce the miscibility between the donors and acceptors, and thus delivers a more desirable vertical phase separation of the active layer. Furthermore, the alloy acceptor of O1-2F and Y6 can optimize the horizontal and vertical morphology of the active layer, forming an effective charge separation and continuous charge transport channel. As a result, the PM6:PM7/Y6:O1-2F system shows a distinct increase in charge mobility and a reduced nonradiative lossΔVnrof 0.231 V compared to the PM6/Y6 system, achieving a high-power conversion efficiency of18.23%. The results indicate that the quaternary LBL OSCs are promising for future large-scale and industrial manufacture of the OSCs.

14 SOLAR ENERGY↗

DOE OSTI2026

Machine Protection Studies of Fast Instability-Driven Beam Losses in the EIC Electron Storage Ring

Fast transverse beam instabilities represent one of the most challenging protection scenarios for the Electron-Ion Collider (EIC) electron storage ring (ESR). Coherent betatron oscillations driven by collective effects can grow rapidly and lead to substantial beam losses within only a few tens of turns, posing a threat to collimators, superconducting (SC) magnets, cryogenic systems, and detector components. A dedicated simulation framework based on distributed transverse excitations was developed to evaluate instability-driven beam losses in the ESR. A conservative instability rise time on the order of 10 turns was used to define a bounding protection scenario. Simulations show that the original two-collimator configuration is insufficient to localize losses under all operating conditions, whereas an optimized four-collimator configuration reduces loss leakage into protected regions by more than two orders of magnitude. Nevertheless, the most severe instability scenarios produce catastrophic beam impacts, with nearly half of the stored beam being intercepted by the collimators within a single turn, exceeding the preliminary estimated survivability limits of the baseline collimator design. The results demonstrate that collimation alone is insufficient to guarantee safe ESR operation and must be complemented by a fast machine-protection system. Dedicated turn-by-turn beam-position monitors and fast beam-loss monitors provide sufficient advance warning to detect instability growth and initiate a beam abort before equipment-damaging loss levels are reached. These findings establish quantitative requirements for the ESR machine-protection architecture and support the implementation of a multilayer protection strategy for both accelerator and detector systems.

43 PARTICLE ACCELERATORS↗

DOE OSTI2023

Label-free quantum super-resolution imaging using entangled multi-mode squeezed light

Abstract In this study, we explore the theoretical application of entangled multi-mode squeezed light for label-free optical super-resolution imaging. By generating massively entangled multi-mode squeezed light through an array of balanced beam splitters, using a single-mode squeezed light input, we create a multi-mode quantum light state with exceptional entanglement and noise suppression below the shot noise level. This significantly reduces imaging measurement errors compared to classical coherent state light imaging when the same number of photons are used on the imaging sample. We demonstrate how to optimize the imaging system’s parameters to achieve the Heisenberg imaging error limit, taking into account the number of entangled modes and photons used. We also examine the effects of optical losses in the imaging system, necessitating adjustments to the optimized parameters based on the degree of optical loss. In practical applications, this new quantum imaging approach reduces the number of photons needed to achieve the same image quality by two orders of magnitude compared to classical imaging methods that use non-entangled, non-squeezed coherent state light.

Soh, Daniel (ORCID:0000000325625675)↗

DOE OSTI2021

Interpretation of hydrothermal conditions, production-injection induced effects, and evidence for enhanced geothermal system-type heat exchange in response to >30 years of production at Roosevelt Hot Springs, Utah, USA

The Roosevelt Hot Springs hydrothermal system is located at the base of the Mineral Mountains in southwestern Utah on the eastern side of the Basin and Range. Hydrothermal activity is related to relatively recent bimodal magmatism, and the system is hosted in coarsely crystalline rock made of Oligocene–Miocene granitoids and Precambrian gneiss. The hydrothermal plume covers ~5 km 2 , with a maximum temperature of 268 °C at ~750 m depth, and a vertically extensive fault-fracture mesh east of the Opal Mound fault controls the upflow of hydrothermal fluids. Power generation (currently 38 MWe gross) began in 1984, and up through 2016, four wells were used for fluid production, and three wells were used for edge-field injection. Chemical analyses of produced fluids show that modern reservoir fluid compositions are similar to but more concentrated than those at the start of production, having near-neutral pH, total dissolved solids of 7000–10,000 mg/kg, and ionic ratios of Cl/HCO 3 ~50–100, Cl/SO4 ~50–100, and Na/K ~4–5. Chemical geothermometers indicate equilibration temperatures that mainly range between 240° and 300 °C. Early production induced a steep drop in pressure (~3.0–3.5 MPa), which was accompanied by a 250–300 m lowering of piezometric levels in wells and development of a shallow steam zone across the system. Hydrothermal fluid compositions evolved continuously in response to production-related steam-loss and injection breakthrough, which is reflected by gradual increases in chloride of up to 35% and stable isotope ratios of up to ~2‰δ 18 O and ~10‰δD. Simple mixing model calculations suggest that there has been a significant amount, ~10–20 MW th , of sustained multi-decadal heat mining and enhanced geothermal system (EGS)–type heat transfer by the injectate as it returns to the production zone. Overall, the two factors that have sustained long-term power production (currently 38 MWe gross) are the increased upflow of deep chloride water and, to a lesser extent, the mining of heat at <1 km depth.

58 GEOSCIENCES↗