Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “functional modeling”

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 379 records · Page 21

Analyzing and Exploring Training Recipes for Large-Scale Transformer-Based Weather Prediction

Abstract The rapid rise of deep learning (DL) in numerical weather prediction (NWP) has led to a proliferation of models which forecast atmospheric variables with comparable or superior skill than traditional physics-based NWP. However, among these leading DL models, there is a wide variance in both the training settings and architecture used. Further, the lack of thorough ablation studies makes it hard to discern which components are most critical to success. In this work, we show that it is possible to attain high forecast skill even with relatively off-the-shelf architectures, simple training procedures, and moderate compute budgets. Specifically, we train a minimally modified Swin Transformer V2 (SwinV2) on ERA5 data and find that it attains superior skill in terms of mean-square errors of deterministic forecasts when compared against the European Centre for Medium-Range Weather Forecasts’ Integrated Forecasting System (IFS). Almost all DL–NWP systems share a core set of hyperparameters and design decisions. To aid and expedite future DL–NWP research, we present an in-depth, systematic exploration of different loss functions, model sizes and depths, patch sizes, and multistep training objectives. We also examine the model performance with metrics beyond the typical accuracy (ACC) and RMSE and investigate how the performance scales with model size. Through our open-source code, scoring pipelines, and models, we share our findings on key aspects of the training pipeline. These ablations reduce the necessity for expensive hyperparameter tuning and lower the barrier to entry for future DL–NWP research. Significance Statement This study investigates the potential of using large-scale transformer-based models for weather prediction, showing that it is possible to achieve high forecast accuracy with simpler, off-the-shelf architectures. By training a minimally modified SwinV2 transformer on ERA5 data, we show that the model achieves competitive forecast skill in terms of mean-square error for key variables, outperforming the European Centre for Medium-Range Weather Forecasts’ Integrated Forecasting System (IFS) at all lead times. Our findings suggest that effective training strategies, such as multistep fine-tuning and channel-weighted losses, significantly enhance the model’s performance. However, we also highlight that these improvements come with trade-offs in other areas, such as ensemble spread and high-frequency spatial detail. This work highlights the promise of deep learning in improving weather forecasts, which could lead to better preparedness and response to weather events, ultimately benefiting society by providing more reliable weather predictions.

Willard, Jared D. [Lawrence Berkeley National Labo↗

US East Coast synthetic aperture radar wind atlas for offshore wind energy

We present the first synthetic aperture radar (SAR) offshore wind atlas of the US East Coast from Georgia to the Canadian border. Images from RADARSAT-1, Envisat, and Sentinel-1A/B are processed to wind maps using the geophysical model function (GMF) CMOD5.N. Extensive comparisons with 6008 collocated buoy observations of the wind speed reveal that biases of the individual systems range from -0.8 to 0.6 m s -1 . Unbiased wind retrievals are crucial for producing an accurate wind atlas, and intercalibration of the SAR observations is therefore applied. Wind retrievals from the intercalibrated SAR observations show biases in the range of to -0.2 to 0.0 m s -1 , while at the same time improving the root-mean-squared error from 1.67 to 1.46 m s -1 . The intercalibrated SAR observations are, for the first time, aggregated to create a wind atlas at the height 10 m a.s.l. (above sea level). The SAR wind atlas is used as a reference to study wind resources derived from the Wind Integration National Dataset Toolkit (WTK), which is based on 7 years of modelling output from the Weather Research and Forecasting (WRF) model. Comparisons focus on the spatial variation in wind resources and show that model outputs lead to lower coastal wind speed gradients than those derived from SAR. Areas designated for offshore wind development by the Bureau of Ocean Energy Management are investigated in more detail; the wind resources in terms of the mean wind speed show spatial variations within each designated area between 0.3 and 0.5 m s -1 for SAR and less than 0.2m s -1 for the WTK. Our findings indicate that wind speed gradients and variations might be underestimated in mesoscale model outputs along the US East Coast.

17 WIND ENERGY↗

Vertical canopy gradients of respiration drive plant carbon budgets and leaf area index

Despite its importance for determining global carbon fluxes, leaf respiration remains poorly constrained in land surface models (LSMs). We tested the sensitivity of the Energy Exascale Earth System Model Land Model – Functionally Assembled Terrestrial Ecosystem Simulator (ELM-FATES) to variation in the canopy gradients of leaf maintenance respiration (R dark ). We ran global and point simulations varying the canopy gradient of R dark to explore the impacts on forest structure, composition, and carbon cycling. In global simulations, steeper canopy gradients of R dark lead to increased understory survival and leaf biomass. Leaf area index (LAI) increased up to 77% in tropical regions compared with the default parameterization, improving alignment with remotely sensed benchmarks. Global vegetation carbon varied from 308 Pg C to 449 Pg C across the ensemble. In tropical forest simulations, steeper gradients of R dark had a large impact on successional dynamics. Results show the importance of canopy gradients in leaf traits and fluxes for determining plant carbon budgets and emergent ecosystem properties such as competitive dynamics, LAI, and vegetation carbon. The high-model sensitivity to canopy gradients in R dark highlights the need for more observations of how leaf traits and fluxes vary along light micro-environments to inform critical dynamics in LSMs.

59 BASIC BIOLOGICAL SCIENCES↗

A Novel Active Optimization Approach for Rapid and Efficient Design Space Exploration Using Ensemble Machine Learning

In this work, a novel design optimization technique based on active learning, which involves dynamic exploration and exploitation of the design space of interest using an ensemble of machine learning algorithms, is presented. In this approach, a hybrid methodology incorporating an explorative weak learner (regularized basis function model) that fits high-level information about the response surface and an exploitative strong learner (based on committee machine) that fits finer details around promising regions identified by the weak learner is employed. For each design iteration, an aristocratic approach is used to select a set of nominees, where points that meet a threshold merit value as predicted by the weak learner are selected for evaluation. In addition to these points, the global optimum as predicted by the strong learner is also evaluated to enable rapid convergence to the actual global optimum once the most promising region has been identified by the optimizer. Additionally, this methodology is first tested by applying it to the optimization of a two-dimensional multi-modal surface and, subsequently, to a complex internal combustion (IC) engine combustion optimization case with nine control parameters related to fuel injection, initial thermodynamic conditions, and in-cylinder flow. It is found that the new approach significantly lowers the number of function evaluations that are needed to reach the optimum design configuration (by up to 80%) when compared to conventional optimization techniques, such as particle swarm and genetic algorithm-based optimization techniques.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Analysis of Ku-band cross section at low incidence angles

This study is using airborne Ku-band data to address questions which have implications for both model function development and for advancing our physical understanding of the sea surface. Concurrent measurements of ocean directional spectra, significant wave height, and mean surface roughness are made using the capabilities of the radar ocean wave spectrometer (ROWS). The NASA/GSFC's ROWS is a 15 GHz pulse compressed radar which is a radar sensor designed to measure the direction of the long wave components using spectral analysis of the tilt induced reflectively modulation. The ROWS are modified to cycle at 50 Hz for the scanning spectrometer antenna and a wide beamwidth nadir altimeter mode. This change allows the sensor to simultaneously measure directional wave spectra, wave height, mean square slope parameter, and small scale surface roughness. The surface stress caused by the wind is widely believed to be the predominant quantity related to the Ku-band radar cross section for a wide range of incidence angles. The complete coverage in the quasi specular region provided by one sensor is essential to understand the uncertainties between the scattering model and what is happening on the surface. For this presentation, special attention is devoted to sort out some measurement of the anisotropy associated with the band of high frequencies. Using the other geophysical parameters, comparisons are then made with the classic spectral form currentlyused to describe the wind impact on the sea surface.

Chapron, B.↗

Towards an Optimal Noise Versus Resolution Trade-Off in Wind Scatterometry

A scatterometer is a radar that measures the normalized radar cross section sigma(sup 0) of the Earth's surface. Over the ocean this signal is related to the wind via the geophysical model function (GMF). The objective of wind scatterometry is to estimate the wind vector field from sigma(sup 0) measurements; however, there are many subtleties that complicate this problem-making it difficult to obtain a unique wind field estimate. Conventionally, wind estimation is split into two stages: a wind retrieval stage in which several ambiguous solutions are obtained, and an ambiguity removal stage in which ambiguities are chosen to produce an appropriate wind vector field estimate. The most common approach to wind field estimation is to grid the scatterometer swath into wind vector cells and estimate wind vector ambiguities independently for each cell. Then, field wise structure is imposed on the solution by an ambiguity selection routine. Although this approach is simple and practical, it neglects field wise structure in the retrieval step and does not account for the spatial correlation imposed by the sampling. This makes it difficult to develop a theoretically appropriate noise versus resolution trade-off using pointwise retrieval. Fieldwise structure may be imposed in the retrieval step using a model-based approach. However, this approach is generally only practical if a low order wind field model is applied, which may discard more information than is desired. Furthermore, model-based approaches do not account for the structure imposed by the sampling. A more general fieldwise approach is to estimate all the wind vectors for all the WVCs simultaneously from all the measurements. This approach can account for structure of the wind field as well as structure imposed by the sampling in the wind retrieval step. Williams and Long in 2010 developed a fieldwise retrieval method based on maximum a posteriori estimation (MAP). This MAP approach can be extended to perform a noise versus resolution trade-off, and deal with ambiguity selection. This paper extends the fieldwise MAP estimation approach and investigates both the noise versus resolution trade-off as well as ambiguity removal in the fieldwise wind retrieval step. The method is then applied to the Sea Winds scatterometer and the results are analyzed. This paper extends the fieldwise MAP estimation approach and investigates both the noise versus resolution trade-off as well as ambiguity removal in the fieldwise wind retrieval step. The method is then applied to the Sea Winds scatterometer and the results are analyzed.

soil moisture↗

Characterization of Humanized Mouse Model of Organophosphate Poisoning and Detection of Countermeasures via MALDI-MSI

Organophosphoate (OP) chemicals are known to inhibit the enzyme acetylcholinesterase (AChE). Studying OP poisoning is difficult because common small animal research models have serum carboxylesterase, which contributes to animals’ resistance to OP poisoning. Historically, guinea pigs have been used for this research; however, a novel genetically modified mouse strain (KIKO) was developed with nonfunctional serum carboxylase (Es1 KO) and an altered acetylcholinesterase (AChE) gene, which expresses the amino acid sequence of the human form of the same protein (AChE KI). KIKO mice were injected with 1xLD50 of an OP nerve agent or vehicle control with or without atropine. After one to three minutes, animals were injected with 35 mg/kg of the currently fielded Reactivator countermeasure for OP poisoning. Postmortem brains were imaged on a Bruker RapifleX ToF/ToF instrument. Data confirmed the presence of increased acetylcholine in OP-exposed animals, regardless of treatment or atropine status. More interestingly, we detected a small amount of Reactivator within the brain of both exposed and unexposed animals; it is currently debated if reactivators can cross the blood–brain barrier. Further, we were able to simultaneously image acetylcholine, the primary affected neurotransmitter, as well as determine the location of both Reactivator and acetylcholine in the brain. This study, which utilized sensitive MALDI-MSI methods, characterized KIKO mice as a functional model for OP countermeasure development.

2-PAM↗

Analysis of Experimental Ice Accretion Data and Assessment of a Thermodynamic Model During Ice Crystal Icing

This paper evaluates a thermodynamic ice crystal icing model that has been previously presented to describe the possible mechanisms of icing within the core of a turbofan jet engine. The model functions between two distinct ice accretions based on a surface energy balance: freeze-dominated icing and melt-dominated icing. Freeze-dominated icing occurs when liquid water (from melted ice crystals) freezes and accretes on a surface along with the existing ice of the impinging water and ice mass. This freeze-dominated icing is characterized as having strong adhesion to the surface. The amount of ice accretion is partially dictated by a freeze fraction, which is the fraction of impinging liquid water that freezes. Melt-dominated icing occurs as unmelted ice on a surface accumulates. This melt-dominated icing is characterized by weakly bonded surface adhesion. The amount of ice accumulation is partially dictated by a melt fraction, which is the fraction of impinging ice crystals that melts. Experimentally observed ice growth rates suggest that only a small fraction of the impinging ice remains on the surface, implying a mass loss mechanism such as splash, runback, bounce, or erosion. The fraction of mass loss must be determined in conjunction with the fraction of freezing liquid water or fraction of melting ice on an icing surface for a given ice growth rate. This mass loss parameter, however, along with the freeze fraction and melt fraction, are the only experimental parameters that are currently not measured directly. Using icing growth rates from ice crystal icing experiments, a methodology that has been previously proposed is used to determine these unknown parameters. This work takes ice accretion data from tests conducted by the National Aeronautics and Space Administration (NASA) at the Glenn Research Center in 2018 that examined the fundamental physics of ice crystal icing. This paper continues evaluation of the thermodynamic model from a previous effort, with additions to the model that account for sub-freezing temperatures that have been observed at the leading edge of the airfoil during icing. The predicted temperatures were generally in good agreement with measured temperatures. Other key findings include the total wet-bulb temperature being a good first order indicator of whether icing is freeze-dominated (sub-freezing values) or melt-dominated (above freezing). Maximum sticking efficiency values, the fraction of impinging mass that adheres to a surface, was calculated to be about 0.2, and retained this maximum value for a range of melt ratios (0.3 to 0.65 and possibly higher), which is defined as the ratio of liquid water content to total water content. Higher air velocities reduced the maximum sticking efficiency and shifted the icing regime to higher melt ratio values. Finally, the leading edge ice accretion angle was found to be related to ice growth (lower growth rates for smaller angles) and melt ratio (smaller melt ratios resulted in smaller angles, likely due to erosion effects).

Ice-Crystal Icing↗

Thermoelectric heating and cooling–integrated dishwasher with thermal energy storage

Household dishwashers must address several performance goals: maximize washing and drying performance, while minimizing cycle duration, energy consumption, and water consumption. This study develops and examines a novel thermoelectric heating and cooling (TEHC) system with thermal energy storage applied to a household dishwasher (DW), aiming to improve the energy and drying performance of a commercially-available dishwasher while maintaining its washing performance, water consumption, and cycle duration. Experimental testing conducted on the novel TEHC-DW system demonstrates an 8.7% reduction in total energy consumption, lowering the per-cycle usage to 0.952 kWh, and a 40% reduction in energy consumption for internal water heating. The novel TEHC-DW system also demonstrates better drying performance, shortening the drying time by 42% to reach the same remaining moisture content as a commercially-available system. In addition, a resistance-capacitance network model is developed that predicts total drying time, total energy consumption, and the highest temperature reached in the tub (53.4 °C). The model accuracy is validated with experimental data (within ±1 K), and the model functions as a design tool via a parametric study to evaluate the next-generation design and the effect of the number of thermoelectric modules and thermoelectric driving force (i.e., current). Overall, this study demonstrates the potential for TEHC technology to improve energy efficiency and drying performance in household dishwashers.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

G-Function Library for Modeling Vertical Bore Ground Heat Exchanger

This library contains g-functions (thermal response functions) for standard, regularly spaced vertical borehole ground heat exchangers. In total, it contains 34, 321 configurations. To permit interpolation, each configuration has g-functions for heights of 24, 48, 96, 192, and 384 m. All the g-functions were calculated with burial depths of 2m, and borehole diameters of 15 to 17.5 cm, depending on height. In configurations with uniform spacing, the spacing between the boreholes is set to 5m, though it can be scaled to other horizontal spacings.

15 GEOTHERMAL ENERGY↗

Methods for Color Center Preserving Hydrogen‐Termination of Diamond

Abstract Chemical functionalization of diamond surfaces by hydrogen is an important method for controlling the charge state of near‐surface fluorescent color centers, an essential process in fabricating devices such as diamond field‐effect transistors and chemical sensors, and a required first step for realizing families of more complex terminations through subsequent chemical processing. In all these cases, termination is typically achieved using hydrogen plasma sources that can etch or damage the diamond, as well as deposited materials or embedded color centers. This work explores alternative methods for lower‐damage hydrogenation of diamond surfaces, specifically the annealing of diamond samples in high‐purity, non‐explosive mixtures of nitrogen and hydrogen gas, and the exposure of samples to microwave hydrogen plasmas in the absence of intentional stage heating. The effectiveness of these methods are characterized by x‐ray photoelectron spectroscopy (XPS), and comparison of the results to density‐functional modelling of the surface hydrogenation energetics implicates surface oxygen ligands as the primary factor limiting the termination quality of annealed samples. Finally, photoluminescence (PL) spectroscopy is used to verify that both the annealing and reduced sample temperature plasma methods are non‐destructive to near‐surface ensembles of nitrogen‐vacancy (NV) centers, in stark contrast to plasma treatments that use heated sample stages.

36 MATERIALS SCIENCE↗

Pure quark and gluon observables in collinear drop

We construct a class of pure quark and gluon observables by using the collinear drop grooming technique. The construction is based on linear combinations of multiple cumulative distributions of the jet mass in collinear drop, whose specific weights are fully predicted perturbatively. This yields observables which obtain their values purely from quarks (or purely from gluons) in a wide region of phase space. We demonstrate this by showing that these observables are effective in two phase space regions, one dominated by perturbative resummation and one dominated by nonperturbative effects. The nonperturbative effects are included using shape functions which only appear as a common factor in the linear combinations constructed. We test this construction using a numerical analysis with next-to-leading logarithmic resummation and various shape function models, as well as analyzing these observables with PYTHIA and VINCIA. Choices for the collinear drop parameters are optimized for experimental use.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Inversion-based correction of Double-Torsion (DT) subcritical crack growth tests for crack profile geometry

Because of its simplicity and the ability to produce a stable, slow-propagating crack, the Double-Torsion (DT) method has been used widely for investigating the critical and subcritical propagation of a slow-propagating tensile (mode-I) crack. However, to determine the complex relationship between the crack velocity $\mathcal{v_c}$ vs. the strain energy release rate $\mathscr{G}$ (or the stress intensity factor K) from laboratory measurements, several corrections must be made to account for the impact of sample and crack geometry. Particularly, DT test typically produces a crack with a curved edge profile instead of a straight line, causing the local $\mathcal{v_c}$ and $\mathscr{G}$ vary along the crack front. The experimentally measured $\mathcal{v_c}$ and $\mathscr{G}$ data merely reflect collective, averaged behavior of the crack. This makes inversion for the intrinsic, “true” crack growth kinetics necessary, based upon the knowledge of the crack geometry. Simple and effective correction methods have been proposed and validated for the slow, chemical-reaction-controlled part (Region I) of the $\mathcal{v_c}$-$\mathscr{G}$ curve. However, reliable methods for the highly nonlinear, transport-dominated part (Region II) and its sudden transition to the dynamic propagation part (Region III) are still lacking. Here we propose a method for determining the intrinsic $\mathcal{v_c}$-$\mathscr{G}$ relationship across all three Regions based upon DT test data, using a simple model function and its numerical inversion. The performance of this approach is examined and demonstrated using both synthetic and laboratory data for subcritical crack growth in soda lime glass.

42 ENGINEERING↗

Squeezing Out Interfacial Solvation: The Role of Hydrogen-Bonding in the Structural and Orientational Freedom of Molecular Self-Assembly

We report bioinspired membrane molecules with improved physical properties and enhanced stability can serve as functional models for conventional lipid or amphiphilic species. Importantly, these molecules can also provide new insights into emergent phenomena that manifest during self-assembly at interfaces. Here, we elucidate the structural response and mechanistic steps underlying the self-assembly of the amphiphilic, charged oligodimethylsiloxane imidazolium cation (ODMS-MIM+) at the air–aqueous interface using Langmuir trough methods with coincident surface-specific vibrational sum-frequency generation (SFG) spectroscopy. We find evidence for a new compression-induced desolvation step that precedes commonly known disordered-to-ordered phase transitions to form nanoscopic assemblies. The experimental data was supported by atomistic molecular dynamics (MD) simulations to provide a detailed mechanistic picture underlying the assembly and the role of water in these phase transitions. The sensitivity of the hydrophobic ODMS tail conformations to compression-owing to distinct water–ODMS interactions and tail–tail solvation properties-offers new strategies for the design of interfaces that can be further used to develop soft-matter electronics and low-dimensional materials using physical and chemical controls.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A Post-translational Histidine–Histidine Cross-Link Enhances Enzymatic Oxygen Reduction Activity with Greater pH Adaptability

Cross-linked protein residues exist as enzyme cofactors to enable or enhance catalytic activities. Despite their importance in nature, the chemical identity of the cross-links is limited to certain amino acid combinations, whose function and the formation mechanism remain insufficiently understood due to the difficulty in isolating native enzymes without the cross-links. Herein, we report the formation and characterization of both His-Tyr and His-His cross-links under oxidative enzymatic turnover conditions in L29H/F33Y/F43H Mb, a structural and functional model of heme-copper oxidase (HCO). The connectivity of the cross-link was characterized as Ne ε2 (His29)-Cd δ2 (His43) by mass spectrometry (LC-MS/MS) and nuclear magnetic resonance (NMR). Interestingly, formation of the cross-link significantly enhances the oxygen reduction activity of the enzyme at neutral or basic pH with higher product specificity. X-ray crystallography has identified a novel Tyr-His cross-link through a Tyr-O-His linkage. Our mechanistic studies indicate the involvement of high-valent heme-iron and the neighboring tyrosine in an oxidative self-processing pathway to generate the cross-link. In conclusion, this work serves as a new example while providing insights into the enzyme cross-link formation, allowing the design of artificial biocatalysts containing these novel cross-links with higher activity and pH adaptability.

Liu, Yiwei [Univ. of Texas, Austin, TX (United Sta↗

Gaussian processes for autonomous data acquisition at large-scale synchrotron and neutron facilities

The execution and analysis of complex experiments are challenged by the vast dimensionality of the underlying parameter spaces. Although an increase in data-acquisition rates should allow broader querying of the parameter space, the complexity of experiments and the subtle dependence of the model function on input parameters remains daunting owing to the sheer number of variables. New strategies for autonomous data acquisition are being developed, with one promising direction being the use of Gaussian process regression (GPR). GPR is a quick, non-parametric and robust approximation and uncertainty quantification method that can be applied directly to autonomous data acquisition. We review GPR-driven autonomous experimentation and illustrate its functionality using real-world examples from large experimental facilities in the USA and France. We introduce the basics of a GPR-driven autonomous loop with a focus on Gaussian processes, and then shift the focus to the infrastructure that needs to be built around GPR to create a closed loop. Finally, the case studies we discuss show that Gaussian-process-based autonomous data acquisition is a widely applicable method that can facilitate the optimal use of instruments and facilities by enabling the efficient acquisition of high-value datasets.

36 MATERIALS SCIENCE↗

Uniqueness of relaxation times determined by dielectric spectroscopy

Dielectric spectroscopy is extremely powerful to study molecular dynamics, because of the very broad frequency range. Often multiple processes superimpose resulting in spectra that expand over several orders of magnitude, with some of the contributions partially hidden. For illustration, we selected two examples, (i) normal mode of high molar mass polymers partially hidden by conductivity and polarization and (ii) contour length fluctuations partially hidden by reptation using the well-studied polyisoprene melts as example. The intuitive approach to describe experimental spectra and to extract relaxation times is the addition of two or more model functions. Here, we use the empirical Havriliak-Negami function to illustrate the ambiguity of the extracted relaxation time, despite an excellent agreement of the fit with experimental data. We show that there are an infinite number of solutions for which a perfect description of experimental data can be achieved. However, a simple mathematical relationship indicates uniqueness of the pairs of the relaxation strength and relaxation time. Sacrificing the absolute value of the relaxation time enables to find the temperature dependence of the parameters with a high accuracy. For the specific cases studied here, the time temperature superposition (TTS) is very useful to confirm the principle. However, the derivation is not based on a specific temperature dependence, hence, independent from the TTS. We compare new and traditional approaches and find the same trend for the temperature dependence. The important advantage of the new technology is the knowledge of the accuracy of the relaxation times. Relaxation times determined from data for which the peak is clearly visible are the same within the experimental accuracy for traditional and new technology. However, for data where a dominant process hides the peak, substantial deviations can be observed. Finally, we conclude that the new approach is particularly helpful for cases in which relaxation times need to be determined without having access to the associated peak position.

36 MATERIALS SCIENCE↗

A Science Gateway for the Repeatable Analysis of Machine Learning Predicted Gravity Anomalies

In recent years, deep learning has become an increasingly popular alternative for modeling in geoscience applications due to its scalability and efficiency. However, the interpretability, compute, data volume, and hyperparameter tuning requirements of deep learning models make development and monitoring difficult. Furthermore, model explainability and communicating results obtained by these models to users or domain experts is a challenge, as domain experts in geoscience also need to have a deep understanding of how those models function in order to support their scientific works. Here, we describe a science gateway and machine learning pipeline for predicting gravity anomalies from geophysical data. The gateway, built on open-source technologies, provides a holistic view of the pipeline through interactive visualizations aimed at enabling efficient exploratory data analysis. The repeatability, reproducibility, and monitoring capabilities of this overall system allow us to iterate and analyze at scale. Using this pipeline and gateway, we can repeatedly produce accurate high-resolution gravity anomaly datasets. By describing the underlying technologies, implementation, and results, here we provide a foundation for the broader adoption of science gateways into cross-cutting geoscience and machine learning research projects as a means to improve the scientific discovery and collaboration in the geophysics and computational sciences community.

58 GEOSCIENCES↗