A. N. Souza, J. He, T. Bischoff, M. Waruszewski, L. Novak, V. Barra, T. Gibson, A. Sridhar, S. Kandala, S. Byrne, L. C. Wilcox, J. Kozdon, F. X. Giraldo, O. Knoth, J. Marshall, R. Ferrari, T. Schneider
{"title":"The Flux-Differencing Discontinuous Galerkin Method Applied to an Idealized Fully Compressible Nonhydrostatic Dry Atmosphere","authors":"A. N. Souza, J. He, T. Bischoff, M. Waruszewski, L. Novak, V. Barra, T. Gibson, A. Sridhar, S. Kandala, S. Byrne, L. C. Wilcox, J. Kozdon, F. X. Giraldo, O. Knoth, J. Marshall, R. Ferrari, T. Schneider","doi":"10.1029/2022MS003527","DOIUrl":"https://doi.org/10.1029/2022MS003527","url":null,"abstract":"<p>Dynamical cores used to study the circulation of the atmosphere employ various numerical methods ranging from finite-volume, spectral element, global spectral, and hybrid methods. In this work, we explore the use of Flux-Differencing Discontinuous Galerkin (FDDG) methods to simulate a fully compressible dry atmosphere at various resolutions. We show that the method offers a judicious compromise between high-order accuracy and stability for large-eddy simulations and simulations of the atmospheric general circulation. In particular, filters, divergence damping, diffusion, hyperdiffusion, or sponge-layers are not required to ensure stability; only the numerical dissipation naturally afforded by FDDG is necessary. We apply the method to the simulation of dry convection in an atmospheric boundary layer and in a global atmospheric dynamical core in the standard benchmark of Held and Suarez (1994, https://doi.org/10.1175/1520-0477(1994)075〈1825:apftio〉2.0.co;2).</p>","PeriodicalId":14881,"journal":{"name":"Journal of Advances in Modeling Earth Systems","volume":"15 4","pages":""},"PeriodicalIF":6.8,"publicationDate":"2023-04-23","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1029/2022MS003527","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"5758318","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Vincent Verjans, Alexander Robel, Andrew F. Thompson, Helene Seroussi
{"title":"Bias Correction and Statistical Modeling of Variable Oceanic Forcing of Greenland Outlet Glaciers","authors":"Vincent Verjans, Alexander Robel, Andrew F. Thompson, Helene Seroussi","doi":"10.1029/2023MS003610","DOIUrl":"https://doi.org/10.1029/2023MS003610","url":null,"abstract":"<p>Variability in oceanic conditions directly impacts ice loss from marine outlet glaciers in Greenland, influencing the ice sheet mass balance. Oceanic conditions are available from Atmosphere-Ocean Global Climate Model (AOGCM) output, but these models require extensive computational resources and lack the fine resolution needed to simulate ocean dynamics on the Greenland continental shelf and close to glacier marine termini. Here, we develop a statistical approach to generate ocean forcing for ice sheet model simulations, which incorporates natural spatiotemporal variability and anthropogenic changes. Starting from raw AOGCM ocean heat content, we apply: (a) a bias-correction using ocean reanalysis, (b) an extrapolation accounting for on-shelf ocean dynamics, and (c) stochastic time series models to generate realizations of natural variability. The bias-correction reduces model errors by ∼25% when compared to independent in-situ measurements. The bias-corrected time series are subsequently extrapolated to fjord mouth locations using relations constrained from available high-resolution regional ocean model results. The stochastic time series models reproduce the spatial correlation, characteristic timescales, and the amplitude of natural variability of bias-corrected AOGCMs, but at negligible computational expense. We demonstrate the efficiency of this method by generating >6,000 time series of ocean forcing for >200 Greenland marine-terminating glacier locations until 2100. As our method is computationally efficient and adaptable to any ocean model output and reanalysis product, it provides flexibility in exploring sensitivity to ocean conditions in Greenland ice sheet model simulations. We provide the output and workflow in an open-source repository, and discuss advantages and future developments for our method.</p>","PeriodicalId":14881,"journal":{"name":"Journal of Advances in Modeling Earth Systems","volume":"15 4","pages":""},"PeriodicalIF":6.8,"publicationDate":"2023-04-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1029/2023MS003610","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"5961821","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Mélissa Cholette, Jason A. Milbrandt, Hugh Morrison, Danahé Paquin-Ricard, Dominik Jacques
{"title":"Combining Triple-Moment Ice With Prognostic Liquid Fraction in the P3 Microphysics Scheme: Impacts on a Simulated Squall Line","authors":"Mélissa Cholette, Jason A. Milbrandt, Hugh Morrison, Danahé Paquin-Ricard, Dominik Jacques","doi":"10.1029/2022MS003328","DOIUrl":"https://doi.org/10.1029/2022MS003328","url":null,"abstract":"<p>The Predicted Particle Properties (P3) bulk microphysics scheme has been recently modified to combine the two major innovations. The triple-moment approach to represent ice, allowing for a freely evolving spectral dispersion of the size distribution, is combined with the predicted liquid fraction, which enables an explicit representation of mixed-phase particles. The impacts of this combination are examined in the context of high-resolution (1-km horizontal grid spacing) simulations of an observed mid-latitude squall line using the Global Environmental Multiscale atmospheric model. The simulation of mixed-phase particles results in a faster squall line propagation speed and stronger cold pool due to greater cooling from the microphysical processes of sublimation, melting and evaporation. There is a reduction in the mass of ice reaching the surface resulting from a decrease in the mean size of melting ice particles aloft with the predicted liquid fraction. Compared to the original double-moment configuration, triple-moment P3 configuration results in larger mean ice sizes at the surface. The reflectivity structure is improved with the new version, now with a more pronounced bright band in the melting zone with the predicted liquid fraction.</p>","PeriodicalId":14881,"journal":{"name":"Journal of Advances in Modeling Earth Systems","volume":"15 4","pages":""},"PeriodicalIF":6.8,"publicationDate":"2023-04-20","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1029/2022MS003328","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"5944382","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"The Dynamics of Megafire Smoke Plumes in Climate Models: Why a Converged Solution Matters for Physical Interpretations","authors":"S. R. Guimond, J. Reisner, M. Dubey","doi":"10.1029/2022MS003432","DOIUrl":"https://doi.org/10.1029/2022MS003432","url":null,"abstract":"<p>As the climate system warms, megafires have become more frequent with devastating effects. A byproduct of these events is the creation of smoke plumes that can rise into the stratosphere and spread across the globe where they reside for many months. To gain a deeper understanding of the plume dynamics, global climate simulations of a megafire were performed at a wide range of grid spacings from 2.0° down to 7 km, including a 7 km nonhydrostatic experiment. The analysis focuses on how the resolved dynamics affects the specification of the plume characteristics such as injection height and black carbon (BC) mass. Prior studies initialize the smoke plume at one or a few grid points and this is shown here to produce severely dissipative dynamics. In order to validate such simulations with observations, enhancements of the plume characteristics to offset the dissipation is necessary. Using a numerically converged simulation, sensitivity tests show that to approximate the observed stratospheric lifetime, a reduction in BC fraction by 50% is necessary for external mixtures. The vorticity dynamics of the plume is also analyzed with a Lagrangian budget to understand the mechanisms responsible for the evolution of a collocated anticyclonic vortex. The results can be distilled down into a simple conceptual model. As the plume rises, the air diverges at the top of the updraft where the largest concentrations of smoke are found. This divergence induces a dilution of the background cyclonic absolute vorticity producing an anticyclonic vortex. Vortex decay occurs from opposite arguments.</p>","PeriodicalId":14881,"journal":{"name":"Journal of Advances in Modeling Earth Systems","volume":"15 4","pages":""},"PeriodicalIF":6.8,"publicationDate":"2023-04-19","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1029/2022MS003432","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"5923104","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
C. Spencer Jones, Qiyu Xiao, Ryan P. Abernathey, K. Shafer Smith
{"title":"Using Lagrangian Filtering to Remove Waves From the Ocean Surface Velocity Field","authors":"C. Spencer Jones, Qiyu Xiao, Ryan P. Abernathey, K. Shafer Smith","doi":"10.1029/2022MS003220","DOIUrl":"https://doi.org/10.1029/2022MS003220","url":null,"abstract":"<p>The Surface Water and Ocean Topography satellite will measure altimetry on scales down to about 15 km: at these scales, the sea-surface-height signature of inertia-gravity waves, including barotropic tides and internal tides, will be visible. However, tides and inertia-gravity waves have little impact on tracer transport. Recent work has shown that Lagrangian filtering can be used to isolate the inertia-gravity wave part of the flow. This manuscript presents a recipe for removing barotropic motions and inertia-gravity waves from the surface velocities and from the sea-surface height (SSH), to estimate the non-wave part of the flow in the Agulhas region of a high-resolution ocean model (LLC4320). First, two methods for removing the barotropic component of SSH variability are presented. Then Lagrangian filtering, a method that accounts for Doppler shifting of high-frequency motions by the low-frequency velocity field, is applied to both the SSH and the ocean surface velocity field. The results of Lagrangian filtering are presented in spectral space. Lagrangian filtering preserves motions that appear super-inertial in the reference frame of the Earth, while other methods do not preserve these motions as effectively. In some locations most of the energy at high frequencies comes from these Doppler shifted balanced motions. We show that the non-wave part of the velocity field that is preserved more effectively by Lagrangian filtering includes convergent motions near regions of frontogenesis.</p>","PeriodicalId":14881,"journal":{"name":"Journal of Advances in Modeling Earth Systems","volume":"15 4","pages":""},"PeriodicalIF":6.8,"publicationDate":"2023-04-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1029/2022MS003220","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"5694845","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Dylan Schlichting, Lixin Qu, Daijiro Kobashi, Robert Hetland
{"title":"Quantification of Physical and Numerical Mixing in a Coastal Ocean Model Using Salinity Variance Budgets","authors":"Dylan Schlichting, Lixin Qu, Daijiro Kobashi, Robert Hetland","doi":"10.1029/2022MS003380","DOIUrl":"https://doi.org/10.1029/2022MS003380","url":null,"abstract":"<p>Numerical mixing, the spurious mixing primarily generated by the discretization of advection, is often significant in estuarine and coastal models due to sharp, energetic fronts. We compare on- and offline estimates of numerical mixing in a submesoscale-resolving realistic simulation of the ocean state over the Texas-Louisiana continental shelf. While offline estimates of numerical mixing differ from online estimates, offline methods may be the only analysis available. We use two methods to estimate numerical mixing offline based on the residuals of the salinity squared <i>s</i><sup>2</sup> and volume-mean salinity variance budgets. The budget overestimates the time-averaged online numerical mixing by 60% at hourly output. The <i>s</i><sup>2</sup> budget compares poorly due to large truncation errors associated with the tendency and advection terms. The residual of the <i>s</i><sup>2</sup> budget starts to converge to the budget as output frequency increases to 10 min—an unrealistic frequency for long-term coastal ocean simulations—but neither method unconditionally converges to the online method and therefore cannot be recommended for generic analysis of numerical mixing. We also investigate the effects of horizontal resolution on numerical mixing using a two-way nested grid with the online method. The volume-integrated numerical mixing constitutes 57% of the bulk physical mixing—the mixing prescribed by the turbulence closure scheme—in the coarse model and may exceed the physical mixing by half an order of magnitude. We find numerical mixing is reduced by 35% on average in the nested model, likely due to new dynamical processes that emerge in the nested simulation.</p>","PeriodicalId":14881,"journal":{"name":"Journal of Advances in Modeling Earth Systems","volume":"15 4","pages":""},"PeriodicalIF":6.8,"publicationDate":"2023-04-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1029/2022MS003380","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"5888804","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Ocean Barotropic Vorticity Balances: Theory and Application to Numerical Models","authors":"Robin Waldman, Hervé Giordani","doi":"10.1029/2022MS003276","DOIUrl":"https://doi.org/10.1029/2022MS003276","url":null,"abstract":"<p>The barotropic vorticity (BV) balance is fundamental when interpreting the ocean gyre circulation. Here we propose an intercomparison of vorticity equations for the depth-integrated flow applied to ocean models. We review four distinct variants of the BV balances, each giving access to diagnostic equations for the depth-integrated ocean circulation, either meridional, across geostrophic contours or its divergence. We then formulate those balances in the Vorticity Balances in NEMO (VoBiN) diagnostic package aimed at the NEMO ocean platform and more generally C-grid ocean models. We show that spatial discretization of the equations of motion have profound implications for those vorticity balances. Finally, we diagnose the main balances of a global ocean climate simulation. In all vorticity balances, topographic torques arise from interactions of the flow with slanting topography. We identify significant spurious topographic torques related to the model's C-grid discretizations, and we suggest ways to address them. In the depth-integrated and BV balances, bottom vortex stretching and bottom pressure torque drive the flow interaction with topography, respectively. Contrary to Sverdrup theory, the wind stress curl, although dominant in the interior Subtropics, becomes a minor player anywhere significant bottom velocities prevail. The geostrophic contour vorticity balance highlights the limits of barotropic models of the ocean circulation through the so-called JEBAR term. Finally, the transport divergence vorticity balance stresses the limitations of Ekman plus geostrophic dynamics for the mass balance closure in ocean models. This framework should encourage ocean modellers to diagnose more routinely momentum and vorticity equations.</p>","PeriodicalId":14881,"journal":{"name":"Journal of Advances in Modeling Earth Systems","volume":"15 4","pages":""},"PeriodicalIF":6.8,"publicationDate":"2023-04-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1029/2022MS003276","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"5693163","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
{"title":"Modifications to Three-Dimensional Turbulence Parameterization for Tropical Cyclone Simulation at Convection-Permitting Resolution","authors":"Gengjiao Ye, Xu Zhang, Hui Yu","doi":"10.1029/2022MS003530","DOIUrl":"https://doi.org/10.1029/2022MS003530","url":null,"abstract":"<p>Adequate representation of the subgrid-scale (SGS) turbulent fluxes associated with convective clouds in the eyewall and rainbands above the boundary layer is important for simulating the formation of tropical cyclone (TC) dynamic and thermal structure, as well as the evolution and intensification of the TC. Two sets of benchmark large-eddy simulations (LESs) for an idealized TC during the rapid intensification and mature stages were conducted. The turbulent transport above the boundary layer in the TC eyewall and rainbands exhibits a remarkable countergradient characteristic, which is poorly represented by the traditional eddy-diffusivity closure. In contrast, the <i>H</i>-gradient closure based on the horizontal gradients of the resolved variables is capable of accurately capturing the countergradient features and exhibiting a spatial distribution of SGS fluxes that mimics much better the coarse-grained fluxes from the LES benchmarks. Moreover, the <i>H</i>-gradient closure allows for the backscatter transfer of energy. By implementing the <i>H</i>-gradient closure into a three-dimensional turbulence parameterization, the TC simulated using the modified parameterization bears closer resemblance to the LES benchmarks in terms of the spatial distribution of SGS fluxes, TC intensity, primary and secondary circulations, and cloud morphology.</p>","PeriodicalId":14881,"journal":{"name":"Journal of Advances in Modeling Earth Systems","volume":"15 4","pages":""},"PeriodicalIF":6.8,"publicationDate":"2023-04-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1029/2022MS003530","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"5695201","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Timothy B. Higgins, Aneesh C. Subramanian, Andre Graubner, Lukas Kapp-Schwoerer, Peter A. G. Watson, Sarah Sparrow, Karthik Kashinath, Sol Kim, Luca Delle Monache, Will Chapman
{"title":"Using Deep Learning for an Analysis of Atmospheric Rivers in a High-Resolution Large Ensemble Climate Data Set","authors":"Timothy B. Higgins, Aneesh C. Subramanian, Andre Graubner, Lukas Kapp-Schwoerer, Peter A. G. Watson, Sarah Sparrow, Karthik Kashinath, Sol Kim, Luca Delle Monache, Will Chapman","doi":"10.1029/2022MS003495","DOIUrl":"https://doi.org/10.1029/2022MS003495","url":null,"abstract":"<p>There is currently large uncertainty over the impacts of climate change on precipitation trends over the US west coast. Atmospheric rivers (ARs) are a significant source of US west coast precipitation and trends in ARs can provide insight into future precipitation trends. There are already a variety of different methods used to identify ARs, but many are used in contexts that are often difficult to apply to large climate datasets due to their computational cost and requirement of integrated vapor transport as an input variable, which can be expensive to compute in climate models at high temporal frequencies. Using deep learning (DL) to track ARs is a unique approach that can alleviate some of the computational challenges that exist in more traditional methods. However, some questions still remain regarding its flexibility and robustness. This research investigates the consistency of a DL methodology of tracking ARs with more established algorithms to demonstrate its high-level performance for future studies.</p>","PeriodicalId":14881,"journal":{"name":"Journal of Advances in Modeling Earth Systems","volume":"15 4","pages":""},"PeriodicalIF":6.8,"publicationDate":"2023-04-15","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1029/2022MS003495","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"5645676","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pierre Boudier, Anthony Fillion, Serge Gratton, Selime Gürol, Sixin Zhang
{"title":"Data Assimilation Networks","authors":"Pierre Boudier, Anthony Fillion, Serge Gratton, Selime Gürol, Sixin Zhang","doi":"10.1029/2022MS003353","DOIUrl":"https://doi.org/10.1029/2022MS003353","url":null,"abstract":"<p>Data Assimilation aims at estimating the posterior conditional probability density functions based on error statistics of the noisy observations and the dynamical system. State of the art methods are sub-optimal due to the common use of Gaussian error statistics and the linearization of the non-linear dynamics. To achieve a good performance, these methods often require case-by-case fine-tuning by using explicit regularization techniques such as inflation and localization. In this paper, we propose a <i>fully data driven deep learning framework</i> generalizing recurrent Elman networks and data assimilation algorithms. Our approach approximates a sequence of prior and posterior densities conditioned on noisy observations using a <i>log-likelihood cost function</i>. By construction our approach can then be used for general nonlinear dynamics and non-Gaussian densities. As a first step, we evaluate the performance of the proposed approach by using fully and partially observed Lorenz-95 system in which the outputs of the recurrent network are fitted to Gaussian densities. We numerically show that our approach, <i>without using any explicit regularization technique</i>, achieves comparable performance to the state-of-the-art methods, IEnKF-Q and LETKF, across various ensemble size.</p>","PeriodicalId":14881,"journal":{"name":"Journal of Advances in Modeling Earth Systems","volume":"15 4","pages":""},"PeriodicalIF":6.8,"publicationDate":"2023-04-13","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"https://onlinelibrary.wiley.com/doi/epdf/10.1029/2022MS003353","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"6111994","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"OA","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}