Search NASA⌕ Search

SEARCH · Search NASA

Results for “deep generative models”

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 811 records · Page 45

Detailed Statistical Models of Host-Based Data for Detection of Malicious Activity

The cybersecurity research community has focused primarily on the analysis and automation of intrusion detection systems by examining network traffic behaviors. Expanding on this expertise, advanced cyber defense analysis is turning to host-based data to use in research and development to produce the next generation network defense tools. The ability to perform deep packet inspection of network traffic is increasingly harder with most boundary network traffic moving to HTTPS. Additionally, network data alone does not provide a full picture of end-to-end activity. These are some of the reasons that necessitate looking at other data sources such as host data. We outline our investigation into the processing, formatting, and storing of the data along with the preliminary results from our exploratory data analysis. In writing this report, it is our goal to aid in guiding future research by providing foundational understanding for an area of cybersecurity that is rich with a variety of complex, categorical, and sparse data, with a strong human influence component. Including suggestions for guiding potential directions for future research.

97 MATHEMATICS AND COMPUTING↗

OmniXAS: A universal deep-learning framework for materials x-ray absorption spectra

X-ray absorption spectroscopy (XAS) is a powerful characterization technique for probing the local chemical environment of absorbing atoms. However, analyzing XAS data presents significant challenges, often requiring extensive, computationally intensive simulations, as well as significant domain expertise. These limitations hinder the development of fast, robust XAS analysis pipelines that are essential in high-throughput studies and for autonomous experimentation. Here, we address these challenges with OmniXAS, a framework that contains a suite of transfer learning approaches for XAS prediction, each uniquely contributing to improved accuracy and efficiency, as demonstrated on the K-edge spectra database covering eight 3⁢d transition metals (Ti–Cu). The OmniXAS framework is built upon three distinct strategies. First, we use M3GNet [Nat. Comput. Sci. 2, 718 (2022)] to derive latent representations of the local chemical environment of absorption sites as input for XAS prediction, achieving significant improvements over conventional featurization techniques. Second, we employ a hierarchical transfer learning strategy, training a universal multitask model across elements before fine-tuning for element-specific predictions. Models based on this cascaded approach after elementwise fine-tuning outperform element-specific models by up to 69%. Third, we implement cross-fidelity transfer learning, adapting a universal model to predict spectra generated by simulation of a different fidelity with a much higher computational cost. This approach improves prediction accuracy by up to 11% over models trained on the target fidelity alone. Our approach significantly boosts the throughput of XAS modeling by orders of magnitude as compared to first-principles simulations and is extendable to XAS prediction for a broader range of elements. The proposed transfer learning framework is generalizable to enhance deep-learning models that target other properties in materials research.

36 MATERIALS SCIENCE↗

WRF Modelling of Deep Convection and Hail for Wind Power Applications

Deep convection and the related occurrence of hail, intense precipitation and wind gusts represent a hazard to a range of energy infrastructure including wind turbine blades. Wind turbine blade leading edge erosion (LEE) is caused by the impact of falling hydrometeors onto rotating wind turbine blades. It is a major source of wind turbine maintenance costs and energy losses from wind farms. In the United States Southern Great Plains (SGP), where there is widespread wind energy development, deep convection and hail events are common, increasing the potential for precipitation-driven LEE. A 25-day Weather Research and Forecasting (WRF) model simulation conducted at convection-permitting resolution and using a detailed microphysics scheme is carried out for the SGP to evaluate the effectiveness in modelling the wind and precipitation conditions relevant to LEE potential. WRF output for these properties is evaluated using radar observations of precipitation (including hail) and reflectivity, in situ wind speed measurements and wind power generation. This research demonstrates some skill for the primary drivers of LEE. Wind speeds, rainfall rates and precipitation totals show good agreement with observations. Here, the occurrence of precipitation during power producing wind speeds is also shown to exhibit fidelity. Hail events frequently occur during periods when wind turbines are rotating and are especially important to LEE in the SGP. The presence of hail is modelled with a mean proportion correct of 0.77 and an odds ratio 4.55. Further research is needed to demonstrate sufficient model performance to be actionable for the wind energy industry and there is evidence for positive model bias in cloud reflectivity.

17 WIND ENERGY↗

A Deep Exon Cryptic Splice Site Promotes Aberrant Intron Retention in a Von Willebrand Disease Patient

A translationally silent single nucleotide mutation in exon 44 (E44) of the von Willebrand factor (VWF) gene is associated with inefficient removal of intron 44 in a von Willebrand disease (VWD) patient. This intron retention (IR) event was previously attributed to reordered E44 secondary structure that sequesters the normal splice donor site. We propose an alternative mechanism: the mutation introduces a cryptic splice donor site that interferes with the function of the annotated site to favor IR. We evaluated both models using minigene splicing reporters engineered to vary in secondary structure and/or cryptic splice site content. Analysis of splicing efficiency in transfected K562 cells suggested that the mutation-generated cryptic splice site in E44 was sufficient to induce substantial IR. Mutations predicted to vary secondary structure at the annotated site also had modest effects on IR and shifted the balance of residual splicing between the cryptic site and annotated site, supporting competition among the sites. Further studies demonstrated that introduction of cryptic splice donor motifs at other positions in E44 did not promote IR, indicating that interference with the annotated site is context dependent. We conclude that mutant deep exon splice sites can interfere with proper splicing by inducing IR.

60 APPLIED LIFE SCIENCES↗

Electrical flicker-noise generated by filling and emptying of impurity states in injectors of quantum-cascade lasers

Free running line-widths (>100 kHz), much broader than intrinsic line-widths ∼100 Hz, of existing quantum-cascade lasers are governed by strong flicker frequency-noise originating from electrical flicker noise. Understanding of microscopic origins of the electrical flicker noises in quantum-cascade lasers is crucially important for the reduction of strength of flicker frequency-noise without assistances of any type of feedback schemes. In this article, an ad hoc model that is based on fluctuating charge-dipoles induced by electron trappings and de-trappings at indispensable impurity states in injector super-lattices of a quantum-cascade laser is proposed, developing theoretical framework based on the model. The validity of the present model is evaluated by comparing theoretical voltage-noise power spectral densities based on the model with experimental ones obtained by using mid-infrared quantum-cascade lasers with designed impurity-positioning. The obtained experimental results on flicker noises, in comparison with the theoretical ones, shed light on physical mechanisms, such as the inherent one due to impurity states in their injectors and extrinsic ones due to surface states on the ridge-walls and due to residual deep traps, for electrical flicker-noise generation in existing mid-infrared quantum-cascade lasers. It is shown theoretically that quasi-delta doping of impurities in their injectors leads to strong suppression of electrical flicker noise by minimization of the dipole length at a certain temperature, for instance ∼300 K and, in turn, is expected to result in substantial narrowing of the free running line-width down below 10 kHz.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Dynamic coevolution of baseflow and multiscale groundwater flow system during prolonged droughts

Field and numerical studies suggest that baseflow is composed of waters from a spectrum of groundwater flow paths termed the Groundwater Flow System (GWFS) – from shallow hillslope contributions to watershed-scale deep circulation originating in headwaters and discharging into lowland rivers. Here, we explore the evolution of the GWFS under prolonged droughts to understand its dynamics and multiscale nature, and to elucidate its role in baseflow generation and recession at the watershed scale. In this work, we consider three drought scenarios of varying severity and simulate groundwater flow in a 2-D cross-section of an idealized watershed with deep permeable bedrock, tracking the evolution of flow paths, baseflow, and residence times during the recession process. We find that baseflow generation at different drainage stages, and within different subwatersheds, is influenced distinctly by flow paths of different scales, depending on the relative strength of the flow paths and the position of the subwatersheds relative to the recharge/discharge zones of the deeper watershed-scale groundwater circulation. Despite having the same local relief, geology, and climate, baseflow from each subwatershed has a distinct recession behavior and time-dependent residence time distribution. Also, the hydraulic and transport characteristics of baseflow generation co-evolve and are strongly affected by the connection state of the water table to subwatersheds. These findings suggest that asynchrony and dissimilarity of baseflow generation from hillslopes under the impact of the watershed-scale groundwater flow, and interactions with local-scale and intermediate-scale groundwater flow, must be taken into account when interpreting baseflow recession data and building conceptual baseflow models at the watershed scale.

54 ENVIRONMENTAL SCIENCES↗

Inverse mapping of properties to composition through generative modeling for designing molten salts

Generative modeling (GM) has been increasingly used for the inverse design and optimization of materials, yet its application to molten salt mixtures remains unexplored despite how a successful approach to the inverse design of molten salts would contribute to efficiently exploiting their customizability and unlocking their advantages in applications, such as energy production and energy storage. This work presents a workflow for the inverse design of molten salts with targeted density values, addressing the challenge of representing these complex mixtures in GM. A dataset of critically evaluated molten salt densities is used to train a variational autoencoder coupled with a predictive deep neural network, which then can be used to generate new molten salt compositions with desired density values. The effectiveness of the approach is demonstrated by designing mixtures with distinct densities and validating the predicted values using ab initio molecular dynamics simulations.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Evaluation of Nuclear Spent Fuel Disposal in Clay-Bearing Rock - Process Model Development and Experimental Studies

This report represents the milestone deliverable M2SF-23SN010301072 “Evaluation of Nuclear Spent Fuel Disposal in Clay-Bearing Rock - Process Model Development and Experimental Studies” The report provides a status update of FY23 activities for the work package Argillite Disposal work packages for the DOE-NE Spent Fuel Waste Form Science and Technology (SFWST) Program. Clay-rich geological media (often referred as shale or argillite) are among the most abundant type of sedimentary rock near the Earth’s surface. Argillaceous rock formations have the following advantageous attributes for deep geological nuclear waste disposal: widespread geologic occurrence, found in stable geologic settings, low permeability, self-sealing properties, low effective diffusion coefficient, high sorption capacity, and have the appropriate depth and thickness to host nuclear waste repository concepts. The DOE R&D program under the Spent Fuel Waste Science Technology (SFWST) campaign has made key progress (through experiment, modeling, and testing) in the study of chemical and physical phenomena that could impact the long-term safety assessment of heat-generating nuclear waste disposition in clay/shale/argillaceous rock. International collaboration activities comprising field-scale heater tests, field data monitoring, and laboratory-scale experiments provide key information on changes to the engineered barrier system (EBS) material exposed high thermal loads. Moreover, consideration of direct disposal of large capacity dual-purpose canisters (DPCs) as part of the back-end SNF waste disposition strategy has generated interest in improving our understanding of the effects of elevated temperatures on the engineered barrier system (EBS) design concepts. Chemical and structural analyses of sampled bentonite material from laboratory tests at elevated temperatures are key to the characterization of thermal effects affecting bentonite clay barrier performance. The knowledge provided by these experiments is crucial to constrain the extent of sacrificial zones in the EBS design during the thermal period. Thermal, hydrologic, mechanical, and chemical (THMC) data collected from heater tests and laboratory experiments have been used in the development, validation, and calibration of THMC simulators to model near-field coupled processes. This information leads to the development of simulation approaches to assess issues on coupled processes involving porous media flow, transport, geomechanical phenomena, chemical interactions with barrier/geologic materials, and the development of EBS concepts. These lines of knowledge are central to the design of deep geological backfilled repository concepts where temperature plays a key role in the EBS behavior, potential interactions with host rock, and long-term performance in the safety assessment.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Topographic controls on a phreatomagmatic maar-diatreme eruption: field and numerical results from the Holocene Dotsero volcano (Colorado, USA)

We analyzed the eruptive products of Holocene Dotsero volcano (Colorado, USA; ~4.15 ka), which recorded evidence of a progression from effusive magmatic activity through explosive phreatomagmatic maar-forming activity to a final explosive magmatic phase. The nature of the deposits suggests that the irregular and mountainous pre-eruptive topography strongly influenced the formation of directed pyroclastic density currents (PDCs) during the phreatomagmatic phase, where the topographically high northern crater rim acted as a barrier and promoted transport to the south. Furthermore, the crater shape strongly controlled the grain size of the final deposits, providing further evidence that deposit grain size can be a misleading proxy for fragmentation or explosion intensity. We test these hypotheses, as gleaned from fieldwork, by implementing a multiphase numerical model to simulate discrete explosions and their resulting PDCs at maar-diatreme volcanoes. We confirm that a crater shape where one rim is higher than the other promotes multiple pulses of particle transport in the direction of the lower rim. The initial PDCs contain both coarse and fine particles, whereas PDCs generated by later pulses contain, predominately, finer particles. This mechanism provides an alternative explanation to the formation of lapilli tuff/tuff couplets often found at maar-diatreme tephra rings (including Dotsero volcano). Additional modeling shows how the crater depth can control grain size of the resulting PDC deposits. In the case of a deep crater and a relatively weak explosion, it is possible that only the fine particles escape following collapse of the eruptive column, resulting in a fine-grained deposit. This work allows us to improve models developed to describe emplacement of PDCs for different crater geometries and informs on the dynamics of maar-diatreme eruptions, but also highlights the limitations of our knowledge of key parameters, especially those related to the processes of magma-water interaction and syn-eruptive hydrogeology, during such events.

58 GEOSCIENCES↗

A Stochastic Study of Flow Anisotropy and Channelling in Open Rough Fractures

The quantification of fluid flow in rough fractures is of high interest for reservoir engineering, especially for deep geothermal applications. Herein, rough self-affine fractures are stochastically generated with incremental shear displacement and geometrically described by two aperture definitions, the vertical aperture a vert and the effective aperture a eff . In order to compare their effect on fracture flow, such as anisotropy and channelling, Local Cubic Law (LCL) model-based 2D fluid flow is simulated. The particularity of this approach is the combination of a stochastic generation of self-affine fractures with a statistical analysis (560 individual realizations) of the impact of the LCL’s aperture constraint on fracture flow. The results show that aperture definition affects the quantitative interpretation of flow anisotropy and channeling as well as the aperture distribution of the fractures with shearing. Higher values of mean aperture for a given fracture are found using a vert , whereas the aperture standard deviation is larger with a eff . In addition, flow anisotropy is significantly sensitive to aperture definition for small shear displacements and shows a relative higher dispersion with a eff . Thus, LCL prediction models based on a vert are expected to lead to higher dispersion of anisotropy results with a higher uncertainty (factor ~ 2). Realizations based on avert lead to an enhanced clustering of high flow rates for higher shearing displacements. This channeling development results in higher total flow rates for these simulations. Finally, these findings support the direct calibration of pre-existing LCL anisotropy simulations based on avert towards more representative results using a eff .

58 GEOSCIENCES↗

SO(3)-invariance of informed-graph-based deep neural network for anisotropic elastoplastic materials

This work examines the frame-invariance (and the lack thereof) exhibited in simulated anisotropic elasto-plastic responses generated from supervised machine learning of classical multi-layer and informed-graph-based neural networks, and proposes different remedies to fix this drawback. The inherent hierarchical relations among physical quantities and state variables in an elasto-plasticity model are first represented as informed, directed graphs, where three variations of the graph are tested. While feed-forward neural networks are used to train path-independent constitutive relations (e.g., elasticity), recurrent neural networks are used to replicate responses that depends on the deformation history, i.e. or path dependent. In dealing with the objectivity deficiency, we use the spectral form to represent tensors and, subsequently, three metrics, the Euclidean distance between the Euler Angles, the distance from the identity matrix, and geodesic on the unit sphere in Lie algebra, can be employed to constitute objective functions for the supervised machine learning. In this, the aim is to minimize the measured distance between the true and the predicted 3D rotation entities. Following this, we conduct numerical experiments on how these metrics, which are theoretically equivalent, may lead to differences in the efficiency of the supervised machine learning as well as the accuracy and robustness of the resultant models. Neural network models trained with tensors represented in component form for a given Cartesian coordinate system are used as a benchmark. Our numerical tests show that, even given the same amount of information and data, the quality of the anisotropic elasto-plasticity model is highly sensitive to the way tensors are represented and measured. The results reveal that using a loss function based on geodesic on the unit sphere in Lie algebra together with an informed, directed graph yield significantly more accurate rotation prediction than the other tested approaches.

42 ENGINEERING↗

Development of a plug-and-play anti-noise module for fault diagnosis of rotating machines in nuclear power plants

The health of rotating machines is crucial to the stable operation and safety of nuclear power plants. However, research on machine learning-based fault diagnosis of rotating machines in the nuclear industry is still in its infancy. The signal noise generated in the plant may negatively affect the effectiveness of the analysis. In this paper, a plug-and-play anti-noise machine learning module is proposed to fill the knowledge and capability gap. The modules are loaded into a convolutional neural network called deep residual network (ResNet) to obtain a new model with noise reduction capability. The basic idea is that the module is able to identify noise features and include them in the subsequent analysis, effectively filtering noise at the feature level of the network. Nine variants of the new model are compared with the original ResNet as well as four classical machine learning models to test the effectiveness of the module and to examine the impact of the module's loading modes on the performance of the new model. Finally, this research helps facilitate the application of machine learning in the plant noise environment.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Comparative evaluation of deep learning workloads for leadership-class systems

Deep learning (DL) workloads and their performance at scale are becoming important factors to consider as we design, develop and deploy next-generation high-performance computing systems. Since DL applications rely heavily on DL frameworks and underlying compute (CPU/GPU) stacks, it is essential to gain a holistic understanding from compute kernels, models, and frameworks of popular DL stacks, and to assess their impact on science-driven, mission-critical applications. At Oak Ridge Leadership Computing Facility (OLCF), we employ a set of micro and macro DL benchmarks established through the Collaboration of Oak Ridge, Argonne, and Livermore (CORAL) to evaluate the AI readiness of our next-generation supercomputers. In this paper, we present our early observations and performance benchmark comparisons between the Nvidia V100 based Summit system with its CUDA stack and an AMD MI100 based testbed system with its ROCm stack. We take a layered perspective on DL benchmarking and point to opportunities for future optimizations in the technologies that we consider.

Yin, Junqi↗

SURVIVAL DEPTH OF ORGANICS IN ICES UNDER LOW-ENERGY ELECTRON RADIATION ({<=}2 keV)

Icy surfaces in our solar system are continually modified and sputtered with electrons, ions, and photons from solar wind, cosmic rays, and local magnetospheres in the cases of Jovian and Saturnian satellites. In addition to their prevalence, electrons specifically are expected to be a principal radiolytic agent on these satellites. Among energetic particles (electrons and ions), electrons penetrate by far the deepest into the ice and could cause damage to organic material of possible prebiotic and even biological importance. To determine if organic matter could survive and be detected through remote sensing or in situ explorations on these surfaces, such as water ice-rich Europa, it is important to obtain accurate data quantifying electron-induced chemistry and damage depths of organics at varying incident electron energies. Experiments reported here address the quantification issue at lower electron energies (100 eV-2 keV) through rigorous laboratory data analysis obtained using a novel methodology. A polycyclic aromatic hydrocarbon molecule, pyrene, embedded in amorphous water ice films of controlled thicknesses served as an organic probe. UV-VIS spectroscopic measurements enabled quantitative assessment of organic matter survival depths in water ice. Eight ices of various thicknesses were studied to determine damage depths more accurately. The electron damage depths were found to be linear, approximately 110 nm keV{sup -1}, in the tested range which is noticeably higher than predictions by Monte Carlo simulations by up to 100%. We conclude that computational simulations underestimate electron damage depths in the energy region {<=}2 keV. If this trend holds at higher electron energies as well, present models utilizing radiation-induced organic chemistry in icy solar system bodies need to be revisited. For interstellar ices of a few micron thicknesses, we conclude that low-energy electrons generated through photoionization processes in the interstellar medium could penetrate through ice grains significantly and trigger organic reactions several hundred nanometers deep-bulk chemistry thus competing with surface chemistry of astrophysical ice grains.

79 ASTRONOMY AND ASTROPHYSICS↗

The SciDAC QuantOm Framework: A composable Workflow

As part of the Scientific Discovery through Advanced Computing (SciDAC) program, the Quantum Chromodynamics Nuclear Tomography (QuantOM) project aims to analyze data from Deep Inelastic Scattering (DIS) experiments conducted at Jefferson Lab and the upcoming Electron Ion Collider. The DIS data analysis is performed on an event-level by leveraging nuclear theory models and accounting for experimental conditions. In order to efficiently run multiple analyses under varying conditions, a composable workflow was designed where each section (theory, experiment, objective minimization, etc.) has its own dedicated module. The optimization, i.e. the fit of theory to experimental data is carried out by deep learning techniques, such as Generative Adversarial Networks (GANs) or Reinforcement Learning (RL). This presentation gives an overview of the current status of the workflow, highlights present and future challenges, and highlights possible extensions to other projects with similar requirements.

Lersch, Daniel↗

Virtual Reality for Shoot/No-Shoot Decision Training in Law Enforcement: A Literature Review and Research Agenda

Virtual reality (VR) can materially improve “shoot / no-shoot” (SNS) training by giving officers realistic, repeatable practice making high-stakes decisions under pressure. Traditional tools—live-fire ranges and video simulators—build basics, but they cannot adapt to each officer in real time or fully mirror the complexity of the field. VR closes that gap by creating immersive scenarios that are safer, more flexible, easier to scale across units, and able to capture objective performance data. SNS decisions are not just about marksmanship; they rely on perception, judgment, memory, and the ability to hold fire when a threat is uncertain. Effective training therefore needs realism, decision complexity, and branching outcomes that reflect the true consequences of choices. These elements strengthen recognition of hostile intent while reducing false positives and building the self-control required in ambiguous situations. VR brings specific advantages: dynamic environments, full-body interaction, and the ability to measure performance with precision—enabling targeted feedback and better transfer of learning to the street. At the same time, responsible deployment must address scenario quality (credible environments and behaviors), lawful decision models, and user wellbeing (appropriate stress levels, comfort, and safety). Sandia’s VIPER Lab is positioned to lead this work. The team combines human-performance science, AI/ML, and VR/AR development with a deep equipment bench (e.g., omnidirectional treadmill, eye-tracking, haptics, multiple HMDs). This ecosystem supports building and validating next-generation SNS training that is immersive, measurable, and trustworthy. Bottom line: Investment in VR-enabled SNS training that blends evidence-based design with careful validation and legal safeguards is expected to pay off in safer, more consistent decision-making and improved community trust, delivered through training that is practical to deploy at scale.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

Methane Hydrate Formation from Enhanced Organic Carbon Burial During Glacial Lowstands: Examples from the Gulf of Mexico

Methane hydrates in fine-grained marine sediments are often found within veins and fractures occupying discrete depth intervals that are surrounded by hydrate-free sediments. As they are not connected with gas sources beneath the base of the methane hydrate stability zone (MHSZ), these isolated hydrate-bearing intervals have been interpreted as formed by in situ microbial methane. We investigate here the hypothesis that these hydrate deposits form in sediments that were deposited during glacial lowstands and contain higher amounts of labile particulate organic carbon (POC), leading to enhanced microbial methanogenesis. During Pleistocene lowstands, river loads are deposited near the steep top of the continental slope and turbidity currents transport organic-rich, fine-grained sediments to deep waters. Faster sedimentation rates during glacial periods result in better preservation of POC because of decreased exposure times to oxic conditions. The net result is that more labile POC enters the methanogenic zone and more methane is generated in these sediments. To test this hypothesis, we apply an advection-diffusion-reaction model with a time-dependent deposition of labile POC at the seafloor controlled by glacioeustatic sea level variations in the last 250 kyr. The model is run for parameters estimated at three sites drilled by the 2009 Gulf of Mexico Joint Industry Project: Walker Ridge in the Terrebonne Basin (WR313-G and WR313-H) and Green Canyon near the canyon embayment into the Sigsbee Escarpment (GC955-H). In the model, gas hydrate forms in sediments with higher labile POC content deposited during the glacial cycle between 230 and 130 kyr (marine isotope stages 6 and 7). The corresponding depth intervals in the three sites contain hydrates, as shown by high bulk electrical resistivities and resistive subvertical fracture fills. This match supports the hypothesis that enhanced POC burial during glacial lowstands can result in hydrate formation from in situ microbial methanogenesis. Our results have implications for carbon cycling during glacial/interglacial cycles and for hydrate accumulation in the MHSZ. In particular, once hydrate-bearing intervals formed during glacial periods are buried beneath the MHSZ and dissociate, gas bubbles can rise and recycle microbial methane into the MHSZ.

Malinverno, Alberto↗