Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “atmospheric model”

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 595 records · Page 33
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Deforestation reshapes land-surface energy-flux partitioning

Land-use and land-cover change significantly modify local land-surface characteristics and water/energy exchanges, which can lead to atmospheric circulation and regional climate changes. In particular, deforestation accounts for a large portion of global land-use changes, which transforms forests into other land cover types, such as croplands and grazing lands. Many previous efforts have focused on observing and modeling land–atmosphere–water/energy fluxes to investigate land–atmosphere coupling induced by deforestation. However, interpreting land–atmosphere–water/energy-flux responses to deforestation is often complicated by the concurrent impacts from shifts in land-surface properties versus background atmospheric forcings. In this study, we used 29 paired FLUXNET sites, to improve understanding of how deforested land surfaces drive changes in surface-energy-flux partitioning. Each paired sites included an intact forested and non-forested site that had similar background climate. We employed transfer entropy, a method based on information theory, to diagnose directional controls between coupling variables, and identify nonlinear cause–effect relationships. Transfer entropy is a powerful tool to detective causal relationships in nonlinear and asynchronous systems. The paired eddy covariance flux measurements showed consistent and strong information flows from vegetation activity (gross primary productivity (GPP)) and physical climate (e.g. shortwave radiation, air temperature) to evaporative fraction (EF) over both non-forested and forested land surfaces. More importantly, the information transfers from radiation, precipitation, and GPP to EF were significantly reduced at non-forested sites, compared to forested sites. We then applied these observationally constrained metrics as benchmarks to evaluate the Energy Exascale Earth System Model (E3SM) land model (ELM). ELM predicted vegetation controls on EF relatively well, but underpredicted climate factors on EF, indicating model deficiencies in describing the relationships between atmospheric state and surface fluxes. Moreover, changes in controls on surface energy flux partitioning due to deforestation were not detected in the model. We highlight the need for benchmarking model simulated surface-energy fluxes and the corresponding causal relationships against those of observations, to improve our understanding of model predictability on how deforestation reshapes land surface energy fluxes.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Uncertain role of clouds in shaping summertime atmosphere-sea ice connections in reanalyses and CMIP6 models

Downwelling longwave radiation (DLR) driven by the atmospheric and cloud conditions in the troposphere is suggested to be a dominant factor to determine the summertime net surface energy budget over the Arctic Ocean and thus plays a key role to shape the September sea ice. We use reanalyses and the self-organizing map (SOM) method to distinguish CMIP6 model performance in replicating the observed strong atmosphere-DLR connection. We find all models can reasonably simulate the linkage between key atmosphere variables and the clear sky DLR but behave differently in replicating the atmosphere-DLR connection due to cloud forcing. In ERA5 and strongly coupled models, tropospheric high pressure is associated with decreased clouds in the mid- and high-levels and increased clouds near the surface. This out-of-phase structure indicates that DLR cloud forcing is nearly neutral, making the clear sky DLR more important to bridge JJA circulation to late-summer sea ice. In MERRA-2 and weakly coupled models, tropospheric clouds display a vertically homogeneous reduction; the cloud DLR is thus strongly reduced due to the cooling effect, which partially cancels out the clear sky DLR and makes the total DLR less efficient to translate circulation forcing to sea ice. The differences of cloud vertical distribution in CMIP6 appear to be differentiated by circulation related relative humidity. Therefore, a better understanding of the discrepancy of different reanalyses and remote sensing products is critical to comprehensively evaluate simulated interactions among circulation, clouds, sea ice and energy budget at the surface in summer.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Utilizing HYSPLIT for Emergency Response Modeling at SRS

At SRS, emergency responders use a variety of tools to detect, track, and mitigate hazardous material releases into the atmosphere. Two models currently used at SRS are Puff-Plume and the Lagrangian Particle Dispersion Model (LPDM), a Gaussian and Lagrangian model, respectively. A decision has been made to replace LPDM with the more widely-supported Hybrid Single-Particle Lagrangian Integrated Trajectory (HYSPLIT) model for evaluating inhalation and ingestion doses following a release. HYSPLIT is designed to compute complex dispersion and deposition simulation To achieve the implementation of HYSPLIT, we have developed a preliminary UI framework that will allow HYSPLIT to be run on ATG computers without the need for active network connections, thus avoiding the loss of capabilities in the event of a network outage during an emergency.

Earley, Ian↗
DOE OSTI2022

Evaluation of Mechanistic and Empirical Models against Existing FFRD and LOCA Experimental Databases

The desire of the nuclear industry to improve the economics of existing nuclear power plants has necessitated research into the potential of a phenomenon known as fuel fragmentation, relocation, and dispersal (FFRD). This phenomenon is possible during a loss-of-coolant accident (LOCA) at relatively high burnup. The Nuclear Energy Advanced Modeling and Simulation program has been developing simulation capabilities for FFRD and LOCA in the BISON fuel performance code for multiple years. This year, the effort has been focused on evaluating new lower length scale informed pulverization thresholds as well as updating and adding new empirical models for various phenomena. Models added or updated this year include a preliminary transient fission gas release model, new high-temperature Zircaloy creep models, a Zircaloy rupture opening area model, and a temperature-dependent emissivity during radiation from the fuel rod to the surrounding atmosphere. The models are verified through implementation tests to demonstrate code correctness after addition to BISON. The models are then assessed against a subset of the existing BISON validation suite for LOCA and FFRD cases. Two new cases, Studsvik Rods 192 and 193, have been added. The results indicate that the inclusion of a bubble pressure evolution model in the bubbles in the high-burnup structure has the largest impact compared to the 3D fracture criterion on reducing the calculated amount of pulverized fuel.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Air temperature and precipitation constraining the modelled wetland methane emissions in a boreal region in northern Europe

Wetland methane responses to temperature and precipitation are studied in a boreal wetland-rich region in northern Europe using ecosystem process models. Six ecosystem models (JSBACH-HIMMELI, LPX-Bern, LPJ-GUESS, JULES, CLM4.5, and CLM5) are compared to multi-model means of ecosystem models and atmospheric inversions from the Global Carbon Project and upscaled eddy covariance flux results for their temperature and precipitation responses and seasonal cycles of the regional fluxes. Two models with contrasting response patterns, LPX-Bern and JSBACH-HIMMELI, are used as priors in atmospheric inversions with Carbon Tracker Europe–CH4 (CTE-CH4) in order to find out how the assimilation of atmospheric concentration data changes the flux estimates and how this alters the interpretation of the flux responses to temperature and precipitation. Inversion moves wetland emissions of both models towards co-limitation by temperature and precipitation. Between 2000 and 2018, periods of high temperature and/or high precipitation often resulted in increased emissions. However, the dry summer of 2018 did not result in increased emissions despite the high temperatures. The process models show strong temperature and strong precipitation responses for the region (51 %–91 % of the variance explained by both). The month with the highest emissions varies from May to September among the models. However, multi-model means, inversions, and upscaled eddy covariance flux observations agree on the month of maximum emissions and are co-limited by temperature and precipitation. The setup of different emission components (peatland emissions, mineral land fluxes) has an important role in building up the response patterns. Considering the significant differences among the models, it is essential to pay more attention to the regional representation of wet and dry mineral soils and periodic flooding which contribute to the seasonality and magnitude of methane fluxes. The realistic representation of temperature dependence of the peat soil fluxes is also important. Furthermore, it is important to use process-based descriptions for both mineral and peat soil fluxes to simulate the flux responses to climate drivers.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Subsetted model output for surface-related variables from a LASSO-CACTI WRF simulation

Surface-related model output from a Weather Research and Forecasting (WRF) simulation from the Large-Eddy Simulation (LES) Atmospheric Radiation Measurement (ARM) Symbiotic Simulation and Observation (LASSO) deep-convection scenario for the Cloud, Aerosol, and Complex Terrain Interactions (CACTI) field campaign. The LASSO-CACTI simulations span grid spacings from 7.5 km to 100 m through the use of four nested domains, labeled d1 through d4. The simulations are of convection near the Sierras de Córdoba mountain range, roughly centered on the ARM Mobile Facility. More information can be found at https://www.arm.gov/capabilities/modeling/lasso. This version of the output is a collection of variables subsetted from the raw WRF output to reduce the file overhead for users not needing the full raw dataset. The subset files are in netCDF format with the variables in these files at the surface and sub-surface levels.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Phase Modeling of the TRAPPIST-1 Planetary Atmospheres

Transiting compact multiplanet systems provide many unique opportunities to characterize the planets, including studies of size distributions, mean densities, orbital dynamics, and atmospheric compositions. The relatively short orbital periods in these systems ensure that events requiring specific orbital locations of the planets (such as primary transit and secondary eclipse points) occur with high frequency. The orbital motion and associated phase variations of the planets provide a means to constrain the atmospheric compositions through measurement of their albedos. Here we describe the expected phase variations of the TRAPPIST-1 system and times of superior conjunction when the summation of phase effects produce maximum amplitudes. We also describe the infrared flux emitted by the TRAPPIST-1 planets and the influence on the overall phase amplitudes. We further present the results from using the global circulation model ROCKE-3D to model the atmospheres of TRAPPIST-1e and TRAPPIST-1f assuming modern Earth and Archean atmospheric compositions. These simulations are used to calculate predicted phase curves for both reflected light and thermal emission components. We discuss the detectability of these signatures and the future prospects for similar studies of phase variations for relatively faint M stars.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

An Electron-scattering Time Delay in Black Hole Accretion Disks

Abstract Universal to black hole X-ray binaries, the high-frequency soft lag gets longer during the hard-to-intermediate state transition, evolving from ≲1 to ∼10 ms. The soft lag production mechanism is thermal disk reprocessing of nonthermal coronal irradiation. X-ray reverberation models account for the light-travel time delay external to the disk, but assume instantaneous reprocessing of the irradiation inside the electron-scattering-dominated disk atmosphere. We model this neglected scattering time delay as a random walk within an α -disk atmosphere, with approximate opacities. To explain soft lag trends, we consider a limiting case of the scattering time delay that we dub the thermalization time delay, t th ; this is the time for irradiation to scatter its way down to the effective photosphere, where it gets thermalized, and then scatter its way back out. We demonstrate that t th plausibly evolves from being inconsequential for low mass accretion rates m ̇ characteristic of the hard state, to rivaling or exceeding the light-travel time delay for m ̇ characteristic of the intermediate state. However, our crude model confines t th to a narrow annulus near peak accretion power dissipation, so cannot yet explain in detail the anomalously long-duration soft lags associated with larger disk radii. We call for time-dependent models with accurate opacities to assess the potential relevance of a scattering delay.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Amazon rainforest photosynthesis increases in response to atmospheric dryness

Earth system models predict that increases in atmospheric and soil dryness will reduce photosynthesis in the Amazon rainforest, with large implications for the global carbon cycle. Using in situ observations, solar-induced fluorescence, and nonlinear machine learning techniques, we show that, in reality, this is not necessarily the case: In many of the wettest parts of this region, photosynthesis and biomass tend to increase with increased atmospheric dryness, despite the associated reductions in canopy conductance to CO 2 . These results can be largely explained by changes in canopy properties, specifically, new leaves flushed during the dry season have higher photosynthetic capacity than the leaves they replace, compensating for the negative stomatal response to increased dryness. As atmospheric dryness will increase with climate change, our study highlights the importance of reframing how we represent the response of ecosystem photosynthesis to atmospheric dryness in very wet regions, to accurately quantify the land carbon sink.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

DeepPhenoMem V1.0: deep learning modelling of canopy greenness dynamics accounting for multi-variate meteorological memory effects on vegetation phenology

Abstract. Vegetation phenology plays a key role in controlling the seasonality of ecosystem processes that modulate carbon, water and energy fluxes between the biosphere and atmosphere. Accurate modelling of vegetation phenology in the interplay of Earth's surface and the atmosphere is thus crucial to understand how the coupled system will respond to and shape climatic changes. Phenology is controlled by meteorological conditions at different timescales: on the one hand, changes in key meteorological variables (temperature, water, radiation) can have immediate effects on the vegetation development; on the other hand, phenological changes can be driven by past environmental conditions, known as memory effects. However, the processes governing meteorological memory effects on phenology are not completely understood, resulting in their limited performance of vegetation phenology represented in land surface models. A deep learning model, specifically a long short-term memory network (LSTM), has the potential to capture and model the meteorological memory effects on vegetation phenology. Here, we apply the LSTM to model the vegetation phenology using meteorological drivers and high-temporal-resolution canopy greenness observations through digital repeat photography by the PhenoCam network. We compare a multiple linear regression model, a no-memory-effect LSTM model and a full-memory-effect LSTM model to predict the whole seasonal greenness trajectory and the corresponding phenological transition dates across 50 sites and 317 site years during 2009–2018, covering deciduous broadleaf forests, evergreen needleleaf forests and grasslands. Results show that the deep learning model outperforms the multiple linear regression model, and the full-memory-effect LSTM model performs better than the no-memory-effect model for all three plant function types (median R2 of 0.878, 0.957 and 0.955 for broadleaf forests, evergreen needleleaf forests and grasslands). We also find that the full-memory-effect LSTM model is capable of predicting the seasonal dynamic variations of canopy greenness and reproducing trends in shifting phenological transition dates. We also performed a sensitivity analysis of the full-memory-effect LSTM model to assess its plausibility, revealing its coherence with established knowledge of vegetation phenology sensitivity to meteorological conditions, particularly changes in temperature. Our study highlights that (1) multi-variate meteorological memory effects play a crucial role in vegetation phenology, and (2) deep learning opens up new avenues for improving the representation of vegetation phenological processes in land surface models via a hybrid modelling approach.

Geology↗
DOE OSTI2025

Changes in the Regional Water Cycle and Their Impact on Societies

ABSTRACT Changes in “blue water”, which is the total supply of fresh water available for human extraction over land, are quite closely related to changes in runoff or equivalently precipitation minus evaporation, . This article examines how climate change‐driven recent past and future changes in the regional water cycle relate to blue water availability and changes in human blue water demand. Although at the largest scales theoretical and numerical model predictions are in broad agreement with observations, at continental scales and below models predict large ranges of possible future and runoff especially at the scale of individual river catchments and for shorter timescale subseasonal floods and droughts. Nevertheless, it is expected that the occurrence and severity of floods will increase and that of droughts may increase, possibly compounded by human‐driven non‐climatic changes such as changes in land use, dam water impoundment, irrigation and extraction of groundwater. Contemporary assessments predict that increases in 21st century human water extraction in many highly‐populated regions are unlikely to be sustainable given projections of future . To reduce uncertainty in future predictions, there is an urgent need to improve modeling of atmospheric, land surface and human processes and how these components are coupled. This should be supported by maintaining the observing network and expanding it to improve measurements of land surface, oceanic and atmospheric variables. This includes the development of satellite observations stable over multiple decades and suitable for building reanalysis datasets appropriate for model evaluation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

New efforts in improving numerical aspects of atmospheric physics parameterizations

This is an invited article for GEWEX Quarterly, the newsletter of the Global Energy and Water Exchanges (GEWEX) program which is a core project of the World Climate Research Programme (https://www.gewex.org/resources/gewex-news/). The article introduces recent research results and ongoing activities being carried out in a SciDAC project funded by DOE that focuses on identifying and addressing numerical issues in physics parameterizations in global atmospheric general circulation models.

global atmospheric model↗
DOE OSTI2025

Improved Representations of Longwave Surface Emissivity to Reduce Surface and Atmospheric Heating Biases in Earth System Models

Abstract Many Earth system models (ESMs) approximate surface emissivity as a broadband constant. This approximation reduces the computational burden, yet omits the spectral structure of emissivity and atmospheric absorption. Neglecting spectral variation in surface emission introduces biases in longwave (LW) atmospheric fluxes and heating. Biases are strongest over surfaces with strongly varying emissivity and minimal atmospheric opacity. We examine these biases over water, ice, and snow surfaces. We partition spectral emissivity into the 16 spectral bands utilized by a single‐column atmospheric radiative transfer model (RRTMG_LW) commonly used in ESMs. We quantify flux and heating biases introduced by broadband assumptions relative to the spectrally resolved case for standard atmospheric profiles over each surface type. Current assumptions tend to overestimate upwelling surface fluxes; for example, the greybody assumption overestimates flux by 1.6 W/m 2 (0.52%) at the bottom of a mid‐latitude winter atmosphere over ice, and by 2.33 (1.0%) at the top of atmosphere. The blackbody assumption tends to artificially cool Earth's surface, stabilizing the lower troposphere. Interestingly, the optimal broadband emissivity can deviate from the Planck‐weighted mean by up to 3% depending on surface type and atmospheric profile. We investigate bias sensitivity to surface temperature, cloud water path, and atmospheric water vapor. Bias is most sensitive to water vapor content, and least sensitive to cloud water path. Lastly, we show that a modified greybody method with updated broadband values can reduce total surface flux bias up to 1.69 , comparable to a five‐band approach and at a fraction of the computational cost.

Manzo, L. [Department of Earth System Science Univ↗
DOE OSTI2022

The ICON-A model for direct QBO simulations on GPUs (version icon-cscs:baf28a514)

Abstract. Classical numerical models for the global atmosphere, as used for numerical weather forecasting or climate research, have been developed for conventional central processing unit (CPU) architectures. This hinders the employment of such models on current top-performing supercomputers, which achieve their computing power with hybrid architectures, mostly using graphics processing units (GPUs). Thus also scientific applications of such models are restricted to the lesser computer power of CPUs. Here we present the development of a GPU-enabled version of the ICON atmosphere model (ICON-A), motivated by a research project on the quasi-biennial oscillation (QBO), a global-scale wind oscillation in the equatorial stratosphere that depends on a broad spectrum of atmospheric waves, which originates from tropical deep convection. Resolving the relevant scales, from a few kilometers to the size of the globe, is a formidable computational problem, which can only be realized now on top-performing supercomputers. This motivated porting ICON-A, in the specific configuration needed for the research project, in a first step to the GPU architecture of the Piz Daint computer at the Swiss National Supercomputing Centre and in a second step to the JUWELS Booster computer at the Forschungszentrum Jülich. On Piz Daint, the ported code achieves a single-node GPU vs. CPU speedup factor of 6.4 and allows for global experiments at a horizontal resolution of 5 km on 1024 computing nodes with 1 GPU per node with a turnover of 48 simulated days per day. On JUWELS Booster, the more modern hardware in combination with an upgraded code base allows for simulations at the same resolution on 128 computing nodes with 4 GPUs per node and a turnover of 133 simulated days per day. Additionally, the code still remains functional on CPUs, as is demonstrated by additional experiments on the Levante compute system at the German Climate Computing Center. While the application shows good weak scaling over the tested 16-fold increase in grid size and node count, making also higher resolved global simulations possible, the strong scaling on GPUs is relatively poor, which limits the options to increase turnover with more nodes. Initial experiments demonstrate that the ICON-A model can simulate downward-propagating QBO jets, which are driven by wave–mean flow interaction.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Evaluation of dynamical downscaling in a fully coupled regional earth system model

A set of decadal simulations has been completed and evaluated for gains using the Regional Arctic System Model (RASM) to dynamically downscale data from a global Earth system model and two atmospheric reanalyses. RASM is a fully coupled atmosphere–land–ocean–sea ice regional Earth system model. Nudging to the forcing data is applied to approximately the top half of the atmospheric domain. RASM simulations were also completed with a modification to the atmospheric physics for evaluating changes to the modeling system. The results show that for the top half of the atmosphere, the RASM simulations follow closely to that of the forcing data, regardless of the forcing data. The results for the lower half of the atmosphere, as well as the surface, show a clustering of atmospheric state and surface fluxes based on the modeling system. At all levels of the atmosphere the imprint of the weather from the forcing data is present as indicated in the pattern of the annual means. Biases, in comparison to reanalyses, are evident in the Earth system model forced simulations for the top half of the atmosphere but are not present in the lower atmosphere. This suggests that bias correction is not needed for fully coupled dynamical downscaling simulations. While the RASM simulations tended to go to the same mean state for the lower atmosphere, there are a differences in the variability and changes of weather patterns across the ensemble of simulations. These differences in the weather result in variances in the sea ice and oceanic states.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Characterization and application of a high-resolution microfluidic device in atmospheric ice nucleation research and integrated science teaching

Ambient ice-nucleating particles (INPs) represent an important subset of aerosol particles that facilitate heterogeneous formation of ice crystals in the atmosphere, which can impact Earth system and energy. In this work, we developed an affordable high-resolution platform to replicate atmospheric immersion freezing through the unique capabilities of microfluidics to measure INP concentration of known ice nucleation active materials, as well as high-latitude soil samples. The new instrument produced comparable and extended freezing efficiency data for a multitude of samples over the temperature range of -5 to -35 °C as compared to an existing freezing assay within known uncertainties. Moreover, the new assay was integrated into science teaching at a primarily undergraduate and Hispanic-serving institute, contributing to fostering a deeper understanding of ice nucleation processes among college students, instilling their enthusiasm for science, and addressing the critical need for a well-informed citizenry. Student research assistants were involved in both research and pedagogical development processes. Three students from our team visited the Pacific Northwest National Laboratory to help characterize physicochemical properties of our test samples and to test modeling of atmospheric ice nucleation using a box model as well as a global climate model (E3SM) for the regional-scale conditions over the North Slope of Alaska (NSA) and Southern Great Plains (SGP). We also assessed student learning experience via hands-on training, research processes, and development of educational materials. On pre- and post-participation surveys and in exit interviews, students reported personal growth as advancements in breadth of experience, research-related skill sets, and regarding their educational and career goals. They also affirmed that the methods enacted incorporated the NASA Minority University Research and Education Program (MUREP) STEM Engagement Best Practices.

58 GEOSCIENCES↗
DOE OSTI2021

Evaluation of Tritium Behavior in Forest Vegetation for the Purpose of Determining Appropriate Non-Zero Deposition Velocities for Tritium Oxide

Modeling and field measurements have been conducted between 2019 and 2021 to assess the practicality of using a non-zero deposition velocity for safety basis estimates of tritium oxide fate and transport modeling. A model was developed which used a complex deposition algorithm designed to assess how tritium oxide would mix based on the turbulent motions and wind speed effects that the forest canopy has on the atmosphere. The model was driven based on measurements of wind and turbulence taken from the Aiken AmeriFlux Tower which measures these properties at five levels located within and just above the forest. The model was then validated against a series of field experiments which were designed to test the model predictions and estimate the deposition velocity occurring over the forest environment at the Savannah River Site. The field releases used deuterium oxide as a surrogate for tritium oxide and was released as a fine mist which rapidly evaporated, creating a gaseous tracer in the atmosphere. Using air samplers, the elevations in deuterium concentration in the air relative to background measurements was assessed and then modeled. Generally, the numerical model tended to underpredict the amount of deuterium being mixed from above the canopy to the forest floor, indicating that the predictions it provides are still conservative relative to what was measured during the field experiments. Across a suite of modeling runs, the 95 th and 99 th percentile deposition velocities were estimated to be 1.2 and 0.7 cm s -1 , respectively. Estimated deposition velocities in the 2021 field experiments, which specifically assessed a release above the forest canopy and its mixing to the surface, predicted deposition velocities ranging from 1.75 to 6.61 cm s -1 . While these field releases do not cover all possible meteorological conditions, it seems appropriate to use a non-zero deposition velocity when performing safety-basis modeling of tritium oxide. The recommendation presented in this report is to use 1.0 cm s -1 . This is between the 95 th and 99 th percentile value estimated from the modeling study, suggesting it should be appropriate for the majority of release scenarios given the model’s apparent conservatism relative to field measurements.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

Equation‐Free Surrogate Modeling of Geophysical Flows at the Intersection of Machine Learning and Data Assimilation

Abstract There is a growing interest in developing data‐driven reduced‐order models for atmospheric and oceanic flows that are trained on data obtained either from high‐resolution simulations or satellite observations. The data‐driven models are non‐intrusive in nature and offer significant computational savings compared to large‐scale numerical models. These low‐dimensional models can be utilized to reduce the computational burden of generating forecasts and estimating model uncertainty without losing the key information needed for data assimilation (DA) to produce accurate state estimates. This paper aims at exploring an equation‐free surrogate modeling approach at the intersection of machine learning and DA in Earth system modeling. With this objective, we introduce an end‐to‐end non‐intrusive reduced‐order modeling (NIROM) framework equipped with contributions in modal decomposition, time series prediction, optimal sensor placement, and sequential DA. Specifically, we use proper orthogonal decomposition (POD) to identify the dominant structures of the flow, and a long short‐term memory network to model the dynamics of the POD modes. The NIROM is integrated within the deterministic ensemble Kalman filter (DEnKF) to incorporate sparse and noisy observations at optimal sensor locations obtained through QR pivoting. The feasibility and the benefit of the proposed framework are demonstrated for the NOAA Optimum Interpolation Sea Surface Temperature (SST) V2 data set. Our results indicate that the NIROM is stable for long‐term forecasting and can model dynamics of SST with a reasonable level of accuracy. Furthermore, the prediction accuracy of the NIROM gets improved by almost one order of magnitude by the DEnKF algorithm.

Pawar, Suraj↗