Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Heat Index”

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 37 records · Page 2

CROCUS Low Cost All-in-One Weather Station AMB-004 Data Argonne National Laboratory Prairie Site

The Ambient Weather WS-2902D (AMB) is a low cost weather station that has become very useful for filling data gaps in harder to deploy locations. These low cost weather stations collect 13 second data, which is averaged to a five minute data output available to users through an Application Programming Interface (API) key. The data files contain measurements for precipitation, temperature, wind chill/heat index, relative humidity, dew point, UV index, solar radiation, wind speed, wind direction, wind gust, and with an external particulate matter 2.5 (PM 2.5) sensor. Having all of these measurements in one condense system allows for fast deploying and dense network capabilities. Three of the AMB weather stations were deployed at the Argonne Testbed for Multiscale Observational Science (ATMOS), a 20-acre prairie site at Argonne National Laboratory in Lemont, Illinois. The instruments are denoted by their three digit identifier (CMS-AMB-xxx) format. The data is presented as daily NetCDF (.nc) files, each containing approximately 24 hours of observations. Files follow the naming convention of: the project (CROCUS), location (atmos), instrument name (CMS-AMB-004), data level (raw, a1), and date (year, month, day). The NetCDF format can be accessed using common scientific software such as Python using xarray, netCDF4 or ACT-DOE.

EARTH SCIENCE > ATMOSPHERE > AEROSOLS > PARTICULAT↗

Nexus of thermal resilience and energy efficiency in buildings: A case study of a nursing home

Extreme weather events become more frequent and severe due to climate change. Although energy efficiency technologies can influence thermal resilience of buildings, they are traditionally studied separately, and their interconnections are rarely quantified. This study developed a methodology of modeling and analysis to provide insights into the nexus of thermal resilience and energy efficiency of buildings. In this work, we conducted a case study of a real nursing home in Florida, where 12 patients died during Hurricane Irma in 2017 due to HVAC system power loss, to understand and quantify how passive and active energy efficiency measures (EEMs) can improve thermal resilience to reduce heat-exposure risk of patients. Results show that passive measures of opening windows and doors for natural ventilation, as well as miscellaneous load reduction, are very effective in eliminating the extreme dangerous occasions. However, to maintain safe conditions, active measures such as on-site power generators and thermal storage are also needed. The nursing home was further studied by changing its location to two other cities: San Francisco (mild climate) and Chicago (cold winter and hot summer). Results revealed that the EEMs' impacts on thermal resilience vary significantly by climate and building characteristics. The study also estimated the costs of EEMs to help stakeholders prioritize the measures. Passive measures that may not save energy may greatly improve thermal resilience, and thus should be considered in building design or retrofit. Findings from this study indicate energy efficiency technologies should be evaluated not only by their energy savings performance but also by their influence on a building's resilience to extreme weather events.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

A Preliminary Thermal Model of the LHe-based SCAPE Cryostat

The SCAPE (SuperConducting Arbitrarily Polarizing Emitter) undulator is under development at the Advanced Photon Source (APS). This new undulator requires a cryostat that will be designed based on expected heat loads. For instance, the expected heating of the beam chamber by electron beam is estimated to be at a level of 182 W – much higher than in planar SCUs. This and other challenges require careful thermal analysis of the LHe-based SCAPE cryostat. A detailed thermal model of the LHe-based SCAPE cryostat has been created in ANSYS. This paper presents calculated cooling capacity and temperatures of the SCAPE cryostat for the static and dynamic heat loads. Index Terms—SCAPE (SuperConducting Arbitrarily Polarizing Emitter), cryocooler, LHe, thermal isolation, thermal conductance.

Shiroyanagi, Y.↗

Differential credibility assessment for statistical downscaling

Climate science is increasingly using (i) ensembles of climate projections from multiple models derived using different assumptions and/or scenarios and (ii) process-oriented diagnostics of model fidelity. Efforts to assign differential credibility to projections and/or models are also rapidly advancing. A framework to quantify and depict the credibility of statistically downscaled model output is presented and demonstrated. Here, the approach employs transfer functions in the form of robust and resilient generalized linear models applied to downscale daily minimum and maximum temperature anomalies at 10 locations using predictors drawn from ERA-Interim reanalysis and two global climate models (GCM; GFDL-ESM2M and MPI-ESM-LR). The downscaled time series are used to derive several impact relevant CLIMDEX temperature indices that are assigned credibility based on (1) the reproduction of relevant large-scale predictors by the GCMs (i.e. fraction of regression beta-weights derived from predictors that are well-reproduced) and (2) the degree of variance in the observations reproduced in the downscaled series following application of a new variance inflation technique. Credibility of the downscaled predictands varies across locations, between the two GCM and is generally higher for minimum temperature than maximum temperature. The differential credibility assessment framework demonstrated here is easy to use and flexible. It can be applied as is to inform decision makers regarding projection confidence, and/or extended to include other components of the transfer functions, and/or used to weight members of a statistically downscaled ensemble.

54 ENVIRONMENTAL SCIENCES↗

Evaluation of Global Climate Models for Use in Energy Analysis

The interplay between energy, climate, and weather is becoming more complex due to increasing contributions of renewable energy generation, energy storage, electrified end uses, and the increasing frequency of extreme weather events. Energy system analyses commonly rely on meteorological inputs to estimate renewable energy generation and energy demand; however, these inputs rarely represent the estimated impacts of future climate change. Climate models and publicly available climate change datasets can be used for this purpose, but the selection of inputs from the myriad of available models and datasets is a nuanced and subjective process. In this work, we assess datasets from various global climate models (GCMs) from the Coupled Model Intercomparison Project Phase 6 (CMIP6). We present evaluations of their skills with respect to the historical climate and comparisons of their future projections of climate change for two climate change scenarios. We present the results for different climatic and energy system regions and include interactive figures in the accompanying software repository. Previous work has presented similar GCM evaluations, but none have presented variables and metrics specifically intended for comprehensive energy systems analysis including impacts on energy demand, thermal cooling, hydropower, water availability, solar energy generation, and wind energy generation. We focus on GCM output meteorological variables that directly affect these energy system components including the representation of extreme values that can drive grid resilience events. The objective of this work is not to recommend the best climate model and dataset for a given analysis, but instead to provide a reference to facilitate the selection of climate models and scenarios in subsequent work.

14 SOLAR ENERGY↗

Selection of Global Climate Model Data for Downscaling With Generative Machine Learning and Use in the Power Planning for Alignment of Climate and Energy Systems Project

The range of results from climate models and scenarios is important to the understanding of uncertainty in power planning analysis. A U.S. Department of Energy-funded analytic project called Power Planning for Alignment of Climate and Energy Systems is developing data and analytic methods to reflect the effects of climate change on key variables for power system planning, as part of the Grid Modernization Lab Consortium. This project will select and prepare global climate model results for use in power system planning models. A related report (Evaluation of Global Climate Models for Use in Energy Analysis) assesses the performance of various global climate models from the Coupled Model Intercomparison Project Phase 6 data archive for their historical skill with respect to energy system performance and for their future projections under multiple climate change scenarios. Building from that report, we describe the selection of a climate scenario (Shared Socioeconomic Pathway [SSP] 2-4.5) and five climate models: TaiESM1, EC-Earth3-CC, GFDL-CM4, EC-Earth3-Veg, and MPI-ESM1-2-HR. We describe the model selection criteria, which were based on the quality of the match between model results under historical conditions and on the representation of the range of future values for several variables. These results will be downscaled via an open-source generative machine learning method called Super-Resolution for Renewable Energy Resource Data with Climate Change Impacts.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Feedstock to Function (F2F) v1

The Feedstock to Function (F2F) tool was designed to help scientists and companies explore viable biofuels and bioproducts early in the R&D cycle to support more productive experimentation, while reducing early-stage exploration from months/years to days/weeks (feedstock-to-function.lbl.gov). The tool focuses on using machine learning to predict biomass-derived molecule properties, while evaluating the cost, benefits, and risks of promising molecules for sustainable aviation fuels. The tool successfully predicts (within 15% of experimental values) high-throughput aviation properties for over 10,000 molecules while enabling users to explore new possibilities and opportunities rapidly and effortlessly. It also links to lightweight life-cycle analysis and techno-economic tools for cost and emissions analyses. Predicted molecule properties include melting point, boiling point, flash point, yield sooting index, and heat of combustion. To date, F2F is more expansive and outperforms several other molecule property prediction models while enabling users (scientists, companies, and policy makers) to explore new possibilities and opportunities rapidly and effortlessly. F2F provides the foundation for developing an adaptive computational tool that predicts properties, cost, benefits, and risk of promising new and uncertified alternative jet fuel pathways and their blending effects.

Rapp, Vi↗

Global impacts of vegetation clumping on regulating land surface heat fluxes

The clumping index (CI) quantifies the non-random distribution of vegetation across space, which regulates the canopy radiative transfer processes and land surface carbon, water, and energy cycles. However, its impact on global surface energy budget, particularly sensible heat fluxes and surface temperature, is not well understood. Additionally, while there have been studies showing significant seasonal variations in CI, the impacts of these variations on surface energy fluxes remain unclear. In this study, we incorporated satellite-derived spatially and temporally explicit CI data into the Community Land Model version 5 (CLM5) to evaluate the effects of CI on global land energy fluxes. Our results showed that including CI increased the global mean sensible heat flux dissipated from ground by 3.9 W m -2 (~18%), while decreasing the global mean vegetation sensible heat flux by 4.9 W m -2 (~65%), resulting in a total sensible heat decrease of 1.0 W m -2 (~3%). In contrast, CI increased the global mean latent heat flux by 0.8 W m -2 (~2%), primarily due to increased evapotranspiration (up to 11 W m -2 ) in tropical regions. We also found considerable impacts of seasonal variations in CI, particularly on sensible heat fluxes from ground and vegetation in evergreen needleleaf forests and deciduous needleleaf forests. Using constant CI rather than considering seasonal variations resulted in significant overestimation and underestimation of the sensible heat fluxes from vegetation and ground, respectively, in boreal summer. In conclusion, these changes in surface energy fluxes caused by CI and its seasonal variations led to up to 1.7 and 0.5 K differences in simulated mean ground temperature. These findings highlight the importance of including CI and considering its seasonal variations in modeling land surface energy fluxes.

54 ENVIRONMENTAL SCIENCES↗

Photothermal heterodyne imaging of micron-sized objects

Micron-sized dye-doped polymer beads were imaged using transmitted/reflected light microscopy and photothermal heterodyne imaging (PHI) measurements. The transmitted/reflected light images show distinct ring patterns that are attributed to diffraction effects and/or internal reflections within the beads. In the PHI experiments pump laser induced heating changes the refractive index and size of the bead, which causes changes in the diffraction pattern and internal reflections. This creates an analogous ring pattern in the PHI images. The ring pattern disappears in both the reflected light and PHI experiments when an incoherent light source is used as a probe. When the beads are imaged in an organic medium heat transfer changes the refractive index of the environment, and gives rise to a ring pattern external to the beads in the PHI images. This causes the beads to appear larger than their physical dimensions in PHI experiments. This external signal does not appear when the beads are imaged in air because the refractive index changes in air are very small.

Bhandari, Janak (ORCID:0000000188423259)↗

A systematic method for selecting molecular descriptors as features when training models for predicting physiochemical properties

Machine learning has proven to be a powerful tool for accelerating biofuel development. Although numerous models are available to predict a range of properties using chemical descriptors, there is a trade-off between interpretability and performance. Neural networks provide predictive models with high accuracy at the expense of some interpretability, while simpler models such as linear regression often lack in accuracy. In addition to model architecture, feature selection is also critical for developing interpretable and accurate predictive models. We present a method for systematically selecting molecular descriptor features and developing interpretable machine learning models without sacrificing accuracy. Our method simplifies the process of selecting features by reducing feature multicollinearity and enables discoveries of new relationships between global properties and molecular descriptors. To demonstrate our approach, we developed models for predicting melting point, boiling point, flash point, yield sooting index, and net heat of combustion with the help of the Tree-based Pipeline Optimization Tool (TPOT). For training, we used publicly available experimental data for up to 8351 molecules. Our models accurately predict various molecular properties for organic molecules (mean absolute percent error (MAPE) ranges from 3.3% to 10.5%) and provide a set of features that are well-correlated to the property. This method enables researchers to explore sets of features that significantly contribute to the prediction of the property, offering new scientific insights. To help accelerate early stage biofuel research and development, we also integrated the data and models into a open-source, interactive web tool.

09 BIOMASS FUELS↗

Fluid dynamic and thermal performance of a slotted cylinder at low Reynolds number

The fluid dynamic and thermal performance of a circular cylinder with a slot parallel to the flow is numerically investigated. The study utilized the semi-implicit finite volume multi-material algorithm MPM-ICE, a component of the Uintah framework. The normalized slot width s/D ranges from 0.1 - 0.3, introducing an additional heat transfer surface area between ~ 10 and ~ 50%, and a mass reduction between ~ 13 and ~ 38% in the cylinder. We assumed two-dimensional incompressible flow and simulated a Reynolds number Re D between 100 and 1000. The slotted cylinders are found to have a total drag force reduction up to ~ 45%, compared to a solid cylinder despite the additional viscous drag force in the slot. Convection heat transfer is enhanced up to ~ 70%. Further, the slotted cylinder performance index, defined as the ratio of the heat rate to the drag force, increases up to maximum of ~ 3, indicating better overall thermal fluid performance. An entropy analysis showed the best performance index occurs at the highest Re D . Correlations for drag coefficient and Nusselt number are proposed along with an entropy optimization method.

42 ENGINEERING↗

Evaluation of the DC bus link capacitors and power transistor modules in the qualification testing of PV inverters

To achieve useful, fair, and cost-effective design qualification and type approval testing for PV inverters, two procedures were developed and demonstrated to evaluate the bus link capacitors and power transistor—and where applicable, its module—for long-term reliability. We evaluate the inverter under simulated use conditions with heat run tests to find the highest operating temperatures reached by these components. The procedures are demonstrated in this study with a 4-kVA class string inverter. The maximum temperature of the DC bus link capacitors was determined to be 76.4°C when the ambient temperature was at the inverter's maximum rated temperature of 60°C with derating occurring. On the other hand, the power transistor module's maximum temperature using the heat sink temperature as the index, 92.6°C, occurred in an ambient of 40.8°C when derating was not occurring. The conditions found for upper temperatures exhibited by these components are, respectively, proposed for implementation in the dry heat and thermal cycling tests in the IEC 62093 “Power conversion equipment for photovoltaic systems – Design qualification testing” international standard draft so that the levels applied in the tests are in proportion to those experienced in the field. Finally as a result, fairer evaluation of inverter reliability is realized.

14 SOLAR ENERGY↗

Attribution of the Australian bushfire risk to anthropogenic climate change

Abstract. Disastrous bushfires during the last months of 2019 and January 2020 affected Australia, raising the question to what extent the risk of these fires was exacerbated by anthropogenic climate change. To answer the question for southeastern Australia, where fires were particularly severe, affecting people and ecosystems, we use a physically based index of fire weather, the Fire Weather Index; long-term observations of heat and drought; and 11 large ensembles of state-of-the-art climate models. We find large trends in the Fire Weather Index in the fifth-generation European Centre for Medium-Range Weather Forecasts (ECMWF) Atmospheric Reanalysis (ERA5) since 1979 and a smaller but significant increase by at least 30 % in the models. Therefore, we find that climate change has induced a higher weather-induced risk of such an extreme fire season. This trend is mainly driven by the increase of temperature extremes. In agreement with previous analyses we find that heat extremes have become more likely by at least a factor of 2 due to the long-term warming trend. However, current climate models overestimate variability and tend to underestimate the long-term trend in these extremes, so the true change in the likelihood of extreme heat could be larger, suggesting that the attribution of the increased fire weather risk is a conservative estimate. We do not find an attributable trend in either extreme annual drought or the driest month of the fire season, September–February. The observations, however, show a weak drying trend in the annual mean. For the 2019/20 season more than half of the July–December drought was driven by record excursions of the Indian Ocean Dipole and Southern Annular Mode, factors which are included in the analysis here. The study reveals the complexity of the 2019/20 bushfire event, with some but not all drivers showing an imprint of anthropogenic climate change. Finally, the study concludes with a qualitative review of various vulnerability and exposure factors that each play a role, along with the hazard in increasing or decreasing the overall impact of the bushfires.

van Oldenborgh, Geert Jan (ORCID:0000000268989535)↗

Enhanced surface urban heat islands due to divergent urban-rural greening trends

Abstract Satellite observations show that the surface urban heat island intensity (SUHII) has been increasing over the last two decades. This is often accompanied by an increased urban-rural contrast of vegetation greenness. However, the contribution of uneven vegetation trends in urban and rural areas to the trend of SUHII is unclear, due to the confounding effects of climate change and changes in man-made infrastructures and anthropogenic heat sources. Here we use a data-model fusion approach to quantify such contributions during the peak growing season. We show that the LAI dif (the urban-rural difference of leaf area index) is increasing ( P < 0.05) in 189 of the selected 228 global megacities. The increasing trend of LAI dif from 2000 to 2019 accounts for about one quarter of the trend in satellite-derived SUHII, and the impact is particularly evident in places with rapid urbanization and rural cropland intensification. The marginal sensitivity of SUHII to LAI dif is the strongest in hot-humid megacities surrounded by croplands and in hot-dry megacities surrounded by mixed woody and herbaceous vegetation. Our study highlights the role of long-term vegetation trends in modulating the trends of urban-rural temperature differences.

54 ENVIRONMENTAL SCIENCES↗

Evaluating Diesel/Biofuel Blends Using Artificial Neural Networks and Linear/Nonlinear Equations

Abstract The use of biomass-derived additives in diesel fuel mixtures has the potential to increase the fuel’s efficiency, decrease the formation of particulate matter during its combustion, and retain the fuel’s behavior in cold weather. To this end, identifying compounds that enable these behaviors is paramount. The present work utilizes a series of linear and non-linear equations in series with artificial neural networks to predict the cetane number, yield sooting index, kinematic viscosity, cloud point, and lower heating value of multi-component blends. Property values of pure components are predicted using artificial neural networks trained with existing experimental data, and these predictions and their expected errors are propagated through linear and non-linear equations to obtain property predictions for multi-component blends. Individual component property prediction errors, defined by blind prediction median absolute error, are 4.91 units, 7.84 units, 0.06 cSt, 4.00 °C, and 0.55 MJ/kg for cetane number, yield sooting index, kinematic viscosity, cloud point, and lower heating value respectively. On average, property predictions for blends are shown to be accurate to within 6% of the blends’ experimental values. Further, a multitude of compounds expected to be produced from catalytically upgrading products of fast pyrolysis are evaluated with respect to their behavior in diesel fuel blends.

09 BIOMASS FUELS↗

DEEPEN 3D PFA Favorability Models and 2D Favorability Maps at Newberry Volcano

DEEPEN stands for DE-risking Exploration of geothermal Plays in magmatic ENvironments. Part of the DEEPEN project involved developing and testing a methodology for a 3D play fairway analysis (PFA) for multiple play types (conventional hydrothermal, superhot EGS, and supercritical). This was tested using new and existing geoscientific exploration datasets at Newberry Volcano. This GDR submission includes images, data, and models related to the 3D favorability and uncertainty models and the 2D favorability and uncertainty maps. The DEEPEN PFA Methodology is based on the method proposed by Poux et al. (2020), which uses the Leapfrog Geothermal software with the Edge extension to conduct PFA in 3D. This method uses all available data to build a 3D geodata model which can be broken down into smaller blocks and analyzed with advanced geostatistical methods. Each data set is imported into a 3D model in Leapfrog and divided into smaller blocks. Conditional queries can then be used to assign each block an index value which conditionally ranks each block's favorability, from 0-5 with 5 being most favorable, for each model (e.g., lithologic, seismic, magnetic, structural). The values between 0-5 assigned to each block are referred to as index values. The final step of the process is to combine all the index models to create a favorability index. This involves multiplying each index model by a given weight and then summing the resulting values. The DEEPEN PFA Methodology follows this approach, but split up by the specific geologic components of each play type. These components are defined as follows for each magmatic play type: 1. Conventional hydrothermal plays in magmatic environments: Heat, fluid, and permeability 2. Superhot EGS plays: Heat, thermal insulation, and producibility (the ability to create and sustain fractures suitable for and EGS reservoir) 3. Supercritical plays: Heat, supercritical fluid, pressure seal, and producibility (the proper permeability and pressure conditions to allow production of supercritical fluid) More information on these components and their development can be found in Kolker et al., 2022. For the purposes of subsurface imaging, it is easier to detect a permeable fluid-filled reservoir than it is to detect separate fluid and permeability components. Therefore, in this analysis, we combine fluid and permeability for conventional hydrothermal plays, and supercritical fluid and producibility for supercritical plays. More information on this process is described in the following sections. We also project the 3D favorability volumes onto 2D surfaces for simplified joint interpretation, and we incorporate an uncertainty component. Uncertainty was modeled using the best approach for the dataset in question, for the datasets where we had enough information to do so. Identifying which subsurface parameters are the least resolved can help qualify current PFA results and focus future efforts in data collection. Where possible, the resulting uncertainty models/indices were weighted using the same weights applied to the respective datasets, and summed, following the PFA methodology above, but for uncertainty. There are two different versions of the Leapfrog model and associated favorability models: - v1.0: The first release in June 2023 - v2.1: The second release, with improvements made to the earthquake catalog (included additional identified events, removed duplicate events), to the temperature model (fixed a deep BHT), and to the index models (updated the seismicity-heat source index models for supercritical and EGS, and the resistivity-insulation index models for all three play types). Also uses the jet color map rather than the magma color map for improved interpretability. - v2.1.1: Updated to include v2.0 uncertainty results (see below for uncertainty model versions) There are two different versions of the associated uncertainty models: - v1.0: The first release in June 2023 - v2.0: The se...

15 GEOTHERMAL ENERGY↗

UR 2 : Ultra-rapid reactivity test for real-time, low-cost quality control of calcined clays

To reduce cement's carbon footprint, there is growing interest in commercial adoption of sustainable SCMs such as calcined clays. However, the existing ASTM standard (R 3 test, C1897) to test the reactivity of such clays takes up to 7 days and cannot be used for real-time quality control in an industrial setting. We address this issue by introducing a 5-min Ultra-Rapid Reactivity (UR 2 ) test. By dissolving 47 clay specimens in 4 M NaOH solutions at 90°C, we report that a dissolution index of 1.54Al + Si correlates strongly to the 7-day R 3 heat (R 2 = 0.92, RMSE = 94.1 J/g). This dissolution index also correlates to the 28-day compressive strength for 14 clay mixtures (R 2 = 0.94, RMSE = 1.7 MPa). This UR 2 test relies on colorimetry and can be conducted via off-the-shelf, low-cost cameras. Overall, our new UR 2 test opens a pathway for real-time, low-cost quality control of calcined clays.

36 MATERIALS SCIENCE↗

DEEPEN 3D PFA Index Models for Exploration Datasets at Newberry Volcano

DEEPEN stands for DE-risking Exploration of geothermal Plays in magmatic ENvironments. As part of the development of the DEEPEN 3D play fairway analysis (PFA) methodology for magmatic plays (conventional hydrothermal, superhot EGS, and supercritical), index models needed to be developed to map values in geoscientific exploration datasets to favorability index values. This GDR submission includes those index models. Index models were created by binning values in exploration datasets into chunks based on their favorability, and then applying a number between 0 and 5 to each chunk, where 0 represents very unfavorable data values and 5 represents very favorable data values. To account for differences in how exploration methods are used to detect each play component, separate index models are produced for each exploration method for each component of each play type. Index models were created using histograms of the distributions of each exploration dataset in combination with literature and input from experts about what combinations of geophysical, geological, and geochemical signatures are considered favorable at Newberry. This is in attempt to create similar sized bins based on the current understanding of how different anomalies map to favorable areas for the different types of geothermal plays (i.e., conventional hydrothermal, superhot EGS, and supercritical). For example, an area of partial melt would likely appear as an area of low density, high conductivity, low vp, and high vp/vs. This means that these target anomalies would be given high (4 or 5) index values for the purpose of imaging the heat source. To account for differences in how exploration methods are used to detect each play component, separate index models are produced for each exploration method for each component of each play type. Index models were produced for the following datasets: - Geologic model - Alteration model - vp/vs - vp - vs - Temperature model - Seismicity (density*magnitude) - Density - Resistivity - Fault distance - Earthquake cutoff depth model

15 GEOTHERMAL ENERGY↗