Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “magnitude”

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 937 records · Page 52
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Fast response of global monsoon area and precipitation to regional carbonaceous aerosols

Monsoon exerts a significant social and economic impact on billions of people around the world. Changes in global monsoon (GM) rainfall under global warming have been well studied; however, the understanding of GM connection with regional aerosol emissions is limited. Using an atmospheric general circulation model (AGCM), we attempt to quantify the response of global monsoon area and precipitation to the removal of carbonaceous aerosols over North America (AM), North Africa (AF), Europe (EU), and Asia (AS). Response to the emission removal, the global average aerosol optical depths (AOD) are reduced by 2.21%, 1.8%, 1.16%, and 3.8% respectively. Additionally, a southward shift in the GM domain is observed, responding to changes in magnitude and position of the Inter-Tropical Convergence Zone, when carbonaceous aerosols are removed over these regions in the Northern Hemisphere. The percentage departure in the area-averaged GM precipitation is found to be 4-fold in magnitude for the removal of emissions over EU when compared to that of AM and AF. A shrinkage in the GM area is also observed, which is maximum for EU (~-30%) compared to AM and AF (~-25%). The change in GM precipitation and area are found to be the least and opposite in nature for AS emissions. This study provides insight into the regional nature of carbonaceous aerosol effects on global precipitation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Applicability of the Delplot method for the determination of catalytic reaction sequences: Hydrodeoxygenation of $\gamma$-valerolactone on Ni 2 P/MCM-41

Contact-time measurements based on sites are used to investigate the reaction sequence of the hydrodeoxygenation of γ-valerolactone (GVL) on Ni 2 P/MCM-41. Here the study is carried out at 300 °C and 0.5 MPa with site times of 0.29–6.1 s. A reaction network is proposed in which pentanoic acid and 2-methyltetrahydrofuran (2-MTHF) are formed in parallel initially and pentanal, butane and pentane are formed subsequently. A first-rank Delplot analysis (selectivity) gives a good fit, but higher-rank Delplots do not. The lack of applicability is explained by the magnitude of the rate constants, in which the rate constant for one of the later steps (the consumption of pentanal) is substantially larger than those for the preceding steps. For this reason, pentanal demonstrated the characteristics of a primary product and consequently affected the discrimination of higher rank products. A simulation of a model sequence with assumed rate constants confirms the reason why the Delplot method does not work for this reaction. Since the rate constants cannot be determined a priori to be of close magnitude, the Delplot method cannot be generally used to ascertain reaction sequences.

42 ENGINEERING↗
DOE OSTI2022

A comparative study of mechanical and chemical durability of non-wetting superhydrophobic and lubricant-infused surfaces

Fabrication of bioinspired non-wetting superhydrophobic surfaces (SHS) and lubricant-infused surfaces (LIS) has been studied extensively on a variety of materials. In contrast, durability of the surfaces exposed to harsh mechanical and chemical environments has been the subject of little attention. This study considers the mechanical and chemical durability of SHS and LIS copper surfaces fabricated via facile electrodeposition and chemical etching methods. The as-fabricated surfaces demonstrate excellent non-wetting characteristics with water contact angle of 160° and sliding angle below 5°. The surfaces are subject to mechanical wear through scratch test and water jet impingement at 15 psi and 20 psi as well as accelerated corrosion following the ASTM-E407 standard. The performance of the electrodeposited and etched non-wetting surfaces is systematically assessed in terms of contact and sliding angles and corrosion rate in a simulated marine environment. All surfaces are shown to be robust to mechanical wear after scratch test, with excellent stability of contact and sliding angles, and up to two orders of magnitude reduced corrosion rate compared to bare copper surface. SHS retained steadfast non-wetting characteristics under high-pressure water jet tests compared to the other surfaces while LIS, regardless of texturing method, showed one to two orders of magnitude reduced corrosion rate compared to bare copper surface throughout water jet impingement and chemical durability tests. The study presents for the first time a systematic comparison of durability of SHS and LIS through a common set of fabrication and testing protocol and helps identify appropriate non-wetting surfaces and fabrication methods based on the use environment.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2021

A molecular dynamics study of domain switching in BiFeO 3 nanofilm under DC electric field

The performance of ferroelectric materials is deeply rooted in the domain switching process triggered by an applied electric field. To understand the effect of controlling factors on the domain switching characteristics, nanoscale studies are essential. In this work, an atomistic-scale computational framework is developed to study the domain switching process of bismuth ferrite (BiFeO 3 ) nanofilms. Through a series of molecular dynamics simulations, the effects of duration and magnitude of an electric field, domain size and nanofilm thickness on the domain switching behavior are investigated. The evolution of domains under different magnitudes and durations of the applied electric field is studied, and the results of unswitched and switched cases are summarized in a switching map. Analysis of various nanofilms reveals that by increasing the film thickness, the switching time is increased. Additionally, the nanofilms with different domain sizes and configurations are examined, and the minimum domain size for complete reversal switching is obtained. Our results are helpful to predict and tune the switching process in ferroelectric materials.

36 MATERIALS SCIENCE↗
DOE OSTI2022

The effect of stress on the migration of He gas bubbles under a thermal gradient in Fe by phase-field modeling

Here a phase-field model is parameterized to study the effect of elastic stresses on the migration of He gas bubbles in Fe under a temperature gradient. Stresses caused by the gas bubble pressure and residual stress in the Fe matrix are considered. The dependence of He bubble migration velocity on the magnitude of the residual stress, average temperature, temperature gradient, and bubble size is measured. In agreement with a theoretical model based on surface diffusion, simulation results demonstrate that He bubbles move towards the high temperature region with velocities in Fe that are orders of magnitude faster than previously reported in UO 2 . It is found that local stresses in the matrix caused by the He bubble have negligible effect on the bubble migration process; however, residual stresses in the Fe matrix, potentially caused by processing or irradiation, can modestly modify bubble kinetics through pressure dependence of the He diffusion coefficients. Compressive residual stress decreases diffusion coefficients for bulk and surface diffusion mechanisms, thus reducing the migration velocity of the gas bubble. In contrast, tensile residual stress increases the diffusion coefficients, resulting in an increase in the gas bubble migration velocity. This pressure dependence is also consistent with a theoretical model. This phase field model lays the foundation for analysis of bubble coalescence-induced fracture in He bubble-containing steels.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Surrogate modeling and optimization of the leaching process in a rare earth elements recovery plant

Critical minerals (CMs) and Rare Earth Elements (REEs) play a vital role in crucial infrastructure technologies such as renewable energy generation and batteries. Recovering them from waste materials has recently been found to significantly reduce environmental impact and supply chain costs related to these materials. In this work, we investigate surrogate modeling techniques aimed to simplify the modeling, simulation, and optimization of the leaching processes involved in CM and REE recovery flowsheets. As there is currently a lack of systematic studies on this topic, we perform extensive computational testing to ascertain which surrogate models are easier to construct and offer high predictive accuracy. Further, our results suggest that sparse quadratic models balance predictive accuracy and computational efficiency. Training and using these surrogates for global optimization of the leaching process requires two orders of magnitude fewer measurements and is up to four orders of magnitude faster than optimizing the original simulation using equation-oriented optimization or derivative-free optimization.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

A review of black carbon in snow and ice and its impact on the cryosphere

Black carbon (BC) has emerged as an important short-lived climate forcer. Due to its light absorption properties, BC can darken the snow/ice surface, affect the energy balance, and further lead to acceleration of the cryospheric melting. By reviewing the recent published literatures, we present an overview of the historical changes, spatial distribution of BC in snow/ice, and how these changes are related to the melting of the cryosphere (e.g., glaciers, snow cover, and sea ice). Ice core records show an abrupt increase of BC concentrations that began in the 1850s and continued throughout the 20th century, which is consistent with an increase in BC emissions owing to industrialization. A decrease in the amount of BC since 1970s in Arctic and European ice cores has been partially attributed to the Clean Air Act. However, in the Himalayas, BC records show a continuous increase during this period. BC concentrations in snow and ice in the mid-latitude regions are one to two orders of magnitude higher than those in the polar regions. In particular, BC concentrations in aged snow and granular ice in ablation areas of mountain glaciers are one to three orders of magnitude higher than those in fresh snow or snowpits in the glacier accumulation areas due to BC accumulation during melting season. BC in the surface snow/ice is responsible for about 20% of the albedo reduction in the Tibetan Plateau during glacier melt season. Globally, observations and modeling results indicate that radiative forcing (RF) induced by BC in snow and ice is highest in the mid-latitudes, ranging from several W m–2 in fresh snow to hundreds of W m–2 in aged snow and granular ice in glacier ablation areas. The large BC-in-snow RF and associated snow albedo feedback lead to an acceleration in the total glacier melt (approximately 20%) and/or a reduction in the duration of the snow cover by several days, resulting in an increase of glacier discharge. Given our limited understanding of the role of BC in cryospheric melting, it is important to synthesize the existing research on the multi-scale processes related to BC in snow and ice to identify the gaps in our understanding of these processes and to propose a path forward to improve the quality of our observations of the aforementioned phenomena to fill these gaps.

Kang, Shichang↗
DOE OSTI2021

Miscellaneous Electric Loads: Characterization and Energy Savings Potential

Over time, miscellaneous electric loads (MELs) are expected to increase both in magnitude and share of residential and commercial building energy consumption. This trend is most apparent in North America, but it is also occurring in Japan and Europe. However, the contribution of MELs to building energy use is not currently well understood, both because the products in this category are transforming rapidly and the definition and classification of MELs is ambiguous. This study estimated the national energy consumption of 36 MELs using best-available data and found them to comprise 12% of delivered electricity to the U.S. residential and commercial building sectors. If 26 of these MELs were replaced with the most energy-efficient product models available on the market, their energy consumption could be halved to 6% of delivered electricity. National energy models will better account for building energy consumption by incorporating the MELs data collected and analyzed for this study, leading to improved policy decisions.

Miscellaneous electric loads, Plug loads, Taxonomy↗
DOE OSTI2021

What are the Impacts on the HVAC System when it Provides Frequency Regulation? – A Comprehensive Case Study with a Multi-Zone Variable Air Volume (VAV) System

Heating, ventilation, and air-conditioning (HVAC) systems in commercial buildings have proven to be a great asset for providing secondary frequency regulation (FR) service to the power grid. However, how the provision of FR service affects the HVAC system Quality of Service (QOS) is still largely unknown. In particular, there exists a controversy regarding the energy efficiency impact. This paper investigates how a HVAC system performance and its associated controls are impacted during and after the FR service. For the assessment, we developed a dynamic Modelica-based medium office building model as a virtual testbed with the Air Handling Unit (AHU) fans providing the FR service. Different regulation capacities (demand levels) and standard test signals under various load profiles in the U.S. are considered. The results show that FR has little or no impact on the HVAC operation and occupant thermal comfort while the building is providing the service if an appropriate regulation capacity is pre-determined. However, if the regulation capacity is overestimated, the quality of the FR service provided by the HVAC system and the HVAC energy efficiency will be affected due to the equipment operating constraints. In addition, the AHU supply air temperature fluctuates, and the controller cannot maintain the desired setpoint as the FR magnitude is increased. After the FR service time, the system takes a certain period to return to normal operation. The recovery period extends with the increase in the FR magnitude. For this case study, there exists a control conflict between the FR control and the existing HVAC control (from ASHRAE Guideline 36-2018) during high load conditions, causing the fan speed to saturate. The fan power consumption could increase significantly; up to 77% increase compared to the baseline was observed. This case study demonstrates that the HVAC system controls need to be enhanced and adapted to minimize the interaction between the FR controller and the indoor environmental controls in order to preserve the HVAC system QOS while (during and after) providing satisfactory FR service.

Lu, Xing↗
DOE OSTI2020

Plate motion in sheared granular fault system

Plate motion near the fault gouge layer, and the elastic interplay between the gouge layer and the plate under stick-slip conditions, is key to understanding the dynamics of sheared granular fault systems. Here, a two-dimensional implementation of the combined finite-discrete element method (FDEM), which merges the finite element method (FEM) and the discrete element method (DEM), is used to explicitly simulate a sheared granular gouge fault system. In this work we focus on investigating the influence of normal load, driving shear velocity and plate stiffness on the velocities and displacements in the direction parallel to the shear direction (x-direction) measured at locations on the upper and lower plates just adjacent to the gouge. The simulations show that during slip phases the magnitudes of the measured velocities on the upper and lower plates are proportional to the normal load and may be inversely proportional to the square root of the plate's shear modulus. Whereas, the driving shear velocity does not show distinct influence on the measured velocities. Additionally, large slip velocities are generally associated with large macroscopic friction coefficient drops. For the models subjected to smaller normal loads, larger shear velocities and with stiffer shear plates, the same magnitude of slip velocity could cause a larger drop of macroscopic friction coefficient. During stick phases, the velocities of the upper and lower plates are respectively slightly greater and slightly smaller than half of the driving shear velocity and are both in the same direction of shear. The shear strain rate of the gouge is calculated from this velocity difference between the upper and lower plate during stick phases and thus the gouge effective shear modulus can be calculated. The results show that the gouge effective shear modulus increases proportionally with normal load, while the influence of shear velocity and plate stiffness on gouge effective shear modulus is minor. The simulations address the dynamics of a laboratory-scale fault gouge system and may aid in revealing the complexities of earthquake frictional dynamics.

58 GEOSCIENCES↗
DOE OSTI2022

Thermal energy storage composites with preformed expanded graphite matrix and paraffin wax for long-term cycling stability and tailored thermal properties

Harvesting solar energy, preventing hot spots in electronics, transport of temperature-sensitive materials, and capture and repurposing of thermal energy require a latent heat thermal energy storage (TES) system to store/discharge heat repeatedly. For the practical application of phase change material (PCM) composites within TES systems, reliable thermal performance throughout its operational lifetime is essential. Nevertheless, the reliability of thermal conductivity in multi-phase composites over relevant numbers (>10 3 ) of melt/freeze cycles has barely been studied, particularly for composites containing fillers for thermal conductivity enhancement. Here, we introduce a preform-type expanded graphite (EG)/paraffin wax composite possessing highly robust heat transfer and storage properties even after 10,000 melt/freeze cycles. To achieve such excellent reliability, comparative studies on the combined influence of fabrication process, particle size, EG vol%, binder amount, and compaction on both magnitude and robustness of thermal conductivity were undertaken. Our parametric study has yielded a trade-off between thermal conductivity and latent heat. Based on our modeling, 20 vol% EG approaches the case where all EG particles are well-connected thermally while 10 vol% EG is close to loosely connected fillers in the matrix. Thermal conductivity of our paraffin composites containing 20 vol% EG (25.1 W·m -1 ·K -1 ) is highest among other EG/paraffin composites without aligned EG in the literature. After 10,000 thermal cycling, negligible conductivity fading was observed for the 10 vol% EG composite, while reduction in latent heat remained within 10% for all 10, 14, 17 and 20 vol% EG samples. Here we anticipate this work provides insight on suitable recipe for desirable magnitude and robustness of thermal conductivity of EG/paraffin composites.

25 ENERGY STORAGE↗
DOE OSTI2025

Adaptive Interface-PINNs (AdaI-PINNs) for transient diffusion: Applications to forward and inverse problems in heterogeneous media

We model transient diffusion in heterogeneous materials using a novel physics-informed neural networks framework (PINNs) termed Adaptive interface physics-informed neural networks or AdaI-PINNs (Roy et al. arXiv preprint arXiv:2406.04626, 2024). AdaI-PINNs utilize different activation functions with trainable slopes tailored to each material region within the computational domain, allowing for a fully automated and adaptive PINNs approach to model interface problems with strongly and weakly discontinuous solutions. To enhance its performance in highly heterogeneous transient diffusion systems, we prescribe a suite of robust practices, including appropriate non-dimensionalization of equations, a biased sampling method, Glorot initialization, and the hard enforcement of boundary and initial conditions. Here we evaluate the efficacy of the proposed method on several benchmark forward and inverse problems. Comparative studies on one-dimensional and two-dimensional benchmark problems reveal that the modified AdaI-PINNs outperform its unmodified counterpart, achieving root-mean-square errors that are at least two orders of magnitude better in forward problems. For inverse problems, the maximum errors in the approximated diffusion coefficients by modified AdaI-PINNs are four orders of magnitude better than those of the unmodified version. Additionally, modified AdaI-PINNs demonstrate improved stability in problems with large material mismatches.

42 ENGINEERING↗
DOE OSTI2024

Concurrent measurement of strain and chemical reaction rates in a calcite grain pack undergoing pressure solution: Evidence for surface-reaction controlled dissolution

Pressure solution is inferred to be a significant contributor to sediment compaction and lithification, especially in carbonate sediments. For a sediment deforming primarily by pressure solution, the compaction rate should be directly related to the rate of calcite dissolution, transport along grain contacts, and calcite reprecipitation. Previous experimental work has shown that there is evidence that deformation in wet calcite grain packs is consistent with control by pressure solution, but considerable ambiguity remains regarding the rate limiting mechanism. We present the results of laboratory compaction experiments designed to directly measure calcite dissolution and precipitation rates (recrystallization rates) concurrently with strain rate to test whether measured rates are consistent with predicted rates both in absolute magnitude and time evolution. Recrystallization rates are measured using trace element chemistry (Sr/Ca, Mg/Ca) and isotopes (87Sr/86Sr) of fluids flowing slowly through a compacting grain pack as it is being triaxially compressed. Imaging techniques are used to characterize the grain contacts and strain effects in the post-experiment grain pack. Our data show that calcite recrystallization rates calculated from all three geochemical parameters are in approximate agreement and that the rates closely track strain rate. The geochemically inferred rates are close to predicted rates in absolute magnitude. Uncertainty in grain contact dimensions makes distinguishing between surface reaction control and diffusion control difficult. Measured reaction rates decrease faster than predicted from standard pressure solution creep flow laws. This inconsistency may indicate that calcite dissolution rates at grain contacts are more complex, and more time-dependent, than suggested by geometric models designed to predict grain contact stresses.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

The solubility of La hydroxide and stability of La 3+ and La hydroxyl complexes at acidic to mildly acidic pH from 25 to 250 °C

The mobility of rare earth elements (REE) in natural hydrothermal systems can be assessed using geochemical modeling, which requires reliable thermodynamic data of relevant aqueous species. In this study, we evaluate the controls of pH and temperature on La speciation and the role of hydroxyl complexes in REE transport at hydrothermal conditions. Batch-type hydrothermal solubility experiments were conducted using synthetic La hydroxide powders equilibrated in perchloric acid-based aqueous solutions at temperatures between 150 and 250 °C and starting pH of 2 to 5. The La hydroxide solubility is retrograde with temperature and displays a strong pH dependence with a decrease in La concentrations from acidic to mildly acidic pH spanning between 3 and 5 orders of magnitude (e.g. log La molality of –2.5 to –7.2 at 250 °C). Thermodynamic optimizations using GEMSFITS allow to retrieve the standard partial molal Gibbs energies for the La 3+ aqua ion and the formation constants for the La hydroxyl species (i.e., LaOH 2+ , La(OH) 2 + , La(OH) 3 0 ) between 25 and 250 °C. A comparison between the experimentally derived thermodynamic properties with the calculated values from the Helgeson-Kirkham-Flowers equation of state parameters indicates an increased divergence with temperature. Discrepancies in standard partial molal Gibbs energies range between ~ 1 – 12 kJ/mol and result in a predicted La hydroxide solubility differing by up to 3 orders of magnitude at 250 °C. Speciation calculations indicate a higher stability of La 3+ and LaOH 2+ over the other La hydroxyl species in the studied pH range of 3.4 to 6. Here, the optimized thermodynamic properties for La aqueous species have important implications for modeling the solubility of REE minerals such as monazite and the mobility of REE in hydrothermal systems.

58 GEOSCIENCES↗
DOE OSTI2025

NdPO 4 solubility and aqueous Neodymium speciation in supercritical fluids: An experimental study at 500–700 °C and 1.7 kbar

A key aspect in the formation of rare earth elements (REE) deposits is the role of REE transport as aqueous REE complexes in supercritical hydrothermal solutions, where the nature of the aqueous complex is controlled by solution composition, temperature and pressure. Despite chloride being considered as one of the most abundant transporting ligands in magmatic-hydrothermal fluids, experimental investigations on the stability of aqueous REE chloride complexes are scarce above 300 °C. In this study, synthetic NdPO 4 crystals were reacted with non-saline and saline (0, 0.05 and 0.5 mNaCl), acidic (0.01 mHCl) aqueous solutions in a series of solubility experiments conducted at 500–700 °C and 1.7 kbar, where the solubilities were determined using a stable Nd isotope ( 145 Nd isotope spike) dilution technique. NdPO 4 solubility ranges between 28 ppm and 10,858 ppm, where solubility increases with both temperature and salinity. At 500 °C, log mNdPO 4 increases from –3.93 to –1.60 and there is a strong correlation between NdPO 4 solubility and NaCl concentrations (slope of 1.2 ± 0.3), indicating stabilization of the Nd chloride aqueous complexes with a stoichiometry corresponding to NdCl 2+ . At 600 °C, this correlation is weaker (slope of 0.4, log mNdPO 4 increases from –2.63 to –1.88) indicating the stabilization of both Nd chloride and hydroxyl species controlling solubility. At 700 °C, NdPO 4 solubility is largely independent of NaCl concentration indicating that solubility is controlled by Nd hydroxyl complexes, where stoichiometry suggests the neutral Nd(OH) 3 0 species is dominant. The solubility product (Ksp) of NdPO4 is derived from experimental data with the relation: log K sp = -41.81 – 0.057T – 20987/T, with T temperature in Kelvin. Comparison of the measured Nd phosphate solubility to thermodynamic predictions using the available Helgeson-Kirkham-Flowers equation of state parameters for aqueous Nd complexes indicate that predictions are up to three orders of magnitude lower compared to experimental observations. This discrepancy is most pronounced in saline solutions, suggesting that thermodynamic properties of the REE chloride species in supercritical fluids require revision. Numerical simulations of fluid-rock interaction between acidic, saline fluids and a Strange Lake felsic mineral assemblage demonstrates that NdPO 4 solubility predictions from models are four to six orders of magnitude lower than those calculated based on empirical fits from experiments, which suggests that acidic, saline fluids may play an important role in mobilizing large amounts of light REE from 450 to 700 °C.

58 GEOSCIENCES↗
DOE OSTI2020

Analysis of curtailment at The Geysers geothermal Field, California

Geothermal energy has traditionally been viewed as a baseload energy source, but the rapid growth of intermittent renewable energy has led to a need for more flexibility in power generation to avoid mandatory curtailment imposed by grid operators. This study of curtailment at The Geysers provides insights into the magnitude, duration, frequency, temporal and spatial distribution, and potential causes of curtailment events between 2013 and 2018. Annual levels of curtailment range during this period from 9 to 47 GW h, representing 0.15 to 0.81 % of the net generation. Most curtailments occurred at the power plants connected to a lower capacity transmission line and may result from transmission constriction. There is a clear link between negative pricing and economic curtailment, especially when solar production is higher. Economic curtailment events tend to be only a few hours and vary in magnitude up to almost 300 MW, whereas transmission-related curtailment events can be up to several weeks in duration. It is likely that curtailment of geothermal power will be an increasing concern, and could be mitigated by flexible generation strategies and increases in energy storage. It is critical to know the nature of curtailment events so that flexible generation options can be assessed properly.

15 GEOTHERMAL ENERGY↗
DOE OSTI2022

Induced seismicity red-light thresholds for enhanced geothermal prospects in the Netherlands

Prospects for geothermal energy in the Netherlands have renewed concerns around induced earthquakes. Risks from induced earthquakes are managed by traffic light protocols (TLPs), where the red-light is chosen as the stop-point before exceeding a tolerance to risk. Here, we simulate post-shut-in earthquake scenarios based on realistic information for the Netherlands. We focus on three risk metrics: aggregates like nuisance and damage impacts and also local personal risk (LPR) – a likelihood of building collapse fatality for an individual. Our results show that the severity of these risks varies spatially by orders of magnitude. Prior induced seismicity (e.g., the 2012 Huizinge event) provides a reference baseline to calibrate the Dutch earthquake risk tolerances. We find that these calibrated risk tolerances are similar to those observed in North America, suggesting an underlying sociological ‘license to operate.’ Furthermore, the use of calibrated risk tolerances results in nuisance concerns completely eclipsing the other two metrics. We compare our results to a hypothetical Groningen geothermal operation and find that our approach sets red-light thresholds approximately one magnitude unit below the ML 3.6 Huizinge event. Overall, our results provide a first-order recommendation for red-light thresholds and proactive management of Dutch enhanced geothermal induced seismicity.

15 GEOTHERMAL ENERGY↗