Articles | Volume 14, issue 5
https://doi.org/10.5194/gmd-14-2289-2021
© Author(s) 2021. This work is distributed under
the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
https://doi.org/10.5194/gmd-14-2289-2021
© Author(s) 2021. This work is distributed under
the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
Extending legacy climate models by adaptive mesh refinement for single-component tracer transport: a case study with ECHAM6-HAMMOZ (ECHAM6.3-HAM2.3-MOZ1.0)
Department of Mathematics, Universität Hamburg, Bundesstrasse 55
20146 Hamburg, Germany
Center for Earth System Research and Sustainability (CEN), Universität Hamburg, Grindelberg 5, 20144 Hamburg, Germany
Department of Meteorology and National Centre for Earth Observation, University of Reading, RG6 6ET, Reading, UK
Konrad Simon
Department of Mathematics, Universität Hamburg, Bundesstrasse 55
20146 Hamburg, Germany
Center for Earth System Research and Sustainability (CEN), Universität Hamburg, Grindelberg 5, 20144 Hamburg, Germany
Jörn Behrens
Department of Mathematics, Universität Hamburg, Bundesstrasse 55
20146 Hamburg, Germany
Center for Earth System Research and Sustainability (CEN), Universität Hamburg, Grindelberg 5, 20144 Hamburg, Germany
Related authors
Jozef Skákala, David Ford, Keith Haines, Amos Lawless, Matthew J. Martin, Philip Browne, Marcin Chrust, Stefano Ciavatta, Alison Fowler, Daniel Lea, Matthew Palmer, Andrea Rochner, Jennifer Waters, Hao Zuo, Deep S. Banerjee, Mike Bell, Davi M. Carneiro, Yumeng Chen, Susan Kay, Dale Partridge, Martin Price, Richard Renshaw, Georgy Shapiro, and James While
Ocean Sci., 21, 1709–1734, https://doi.org/10.5194/os-21-1709-2025, https://doi.org/10.5194/os-21-1709-2025, 2025
Short summary
Short summary
UK marine data assimilation (MDA) involves a closely collaborating research community. In this paper, we offer both an overview of the state of the art and a vision for the future across all of the main areas of UK MDA, ranging from physics to biogeochemistry to coupled DA. We discuss the current UK MDA stakeholder applications, highlight theoretical developments needed to advance our systems, and reflect upon upcoming opportunities with respect to hardware and observational missions.
Yumeng Chen, Lars Nerger, and Amos S. Lawless
EGUsphere, https://doi.org/10.5194/egusphere-2024-1078, https://doi.org/10.5194/egusphere-2024-1078, 2024
Short summary
Short summary
In this paper, we present pyPDAF, a Python interface to the parallel data assimilation framework (PDAF) allowing for coupling with Python-based models. We demonstrate the capability and efficiency of pyPDAF under a coupled data assimilation setup.
Yumeng Chen, Polly Smith, Alberto Carrassi, Ivo Pasmans, Laurent Bertino, Marc Bocquet, Tobias Sebastian Finn, Pierre Rampal, and Véronique Dansereau
The Cryosphere, 18, 2381–2406, https://doi.org/10.5194/tc-18-2381-2024, https://doi.org/10.5194/tc-18-2381-2024, 2024
Short summary
Short summary
We explore multivariate state and parameter estimation using a data assimilation approach through idealised simulations in a dynamics-only sea-ice model based on novel rheology. We identify various potential issues that can arise in complex operational sea-ice models when model parameters are estimated. Even though further investigation will be needed for such complex sea-ice models, we show possibilities of improving the observed and the unobserved model state forecast and parameter accuracy.
Bo Dong, Ross Bannister, Yumeng Chen, Alison Fowler, and Keith Haines
Geosci. Model Dev., 16, 4233–4247, https://doi.org/10.5194/gmd-16-4233-2023, https://doi.org/10.5194/gmd-16-4233-2023, 2023
Short summary
Short summary
Traditional Kalman smoothers are expensive to apply in large global ocean operational forecast and reanalysis systems. We develop a cost-efficient method to overcome the technical constraints and to improve the performance of existing reanalysis products.
Tobias Sebastian Finn, Charlotte Durand, Alban Farchi, Marc Bocquet, Yumeng Chen, Alberto Carrassi, and Véronique Dansereau
The Cryosphere, 17, 2965–2991, https://doi.org/10.5194/tc-17-2965-2023, https://doi.org/10.5194/tc-17-2965-2023, 2023
Short summary
Short summary
We combine deep learning with a regional sea-ice model to correct model errors in the sea-ice dynamics of low-resolution forecasts towards high-resolution simulations. The combined model improves the forecast by up to 75 % and thereby surpasses the performance of persistence. As the error connection can additionally be used to analyse the shortcomings of the forecasts, this study highlights the potential of combined modelling for short-term sea-ice forecasting.
Sukun Cheng, Yumeng Chen, Ali Aydoğdu, Laurent Bertino, Alberto Carrassi, Pierre Rampal, and Christopher K. R. T. Jones
The Cryosphere, 17, 1735–1754, https://doi.org/10.5194/tc-17-1735-2023, https://doi.org/10.5194/tc-17-1735-2023, 2023
Short summary
Short summary
This work studies a novel application of combining a Lagrangian sea ice model, neXtSIM, and data assimilation. It uses a deterministic ensemble Kalman filter to incorporate satellite-observed ice concentration and thickness in simulations. The neXtSIM Lagrangian nature is handled using a remapping strategy on a common homogeneous mesh. The ensemble is formed by perturbing air–ocean boundary conditions and ice cohesion. Thanks to data assimilation, winter Arctic sea ice forecasting is enhanced.
Yumeng Chen, Alberto Carrassi, and Valerio Lucarini
Nonlin. Processes Geophys., 28, 633–649, https://doi.org/10.5194/npg-28-633-2021, https://doi.org/10.5194/npg-28-633-2021, 2021
Short summary
Short summary
Chaotic dynamical systems are sensitive to the initial conditions, which are crucial for climate forecast. These properties are often used to inform the design of data assimilation (DA), a method used to estimate the exact initial conditions. However, obtaining the instability properties is burdensome for complex problems, both numerically and analytically. Here, we suggest a different viewpoint. We show that the skill of DA can be used to infer the instability properties of a dynamical system.
Maša Ann, Jörn Behrens, and Jana Sillmann
EGUsphere, https://doi.org/10.5194/egusphere-2025-3505, https://doi.org/10.5194/egusphere-2025-3505, 2025
This preprint is open for discussion and under review for Nonlinear Processes in Geophysics (NPG).
Short summary
Short summary
We present a new framework based on Dynamic Mode Decomposition (DMD) to better detect outliers and model extremes. Unlike standard DMD, which focuses on average system behavior, our approach targets rare, exceptional dynamics. Applied to climate data, it improves extreme event approximation and reveals meaningful spatio-temporal patterns. The method may generalize to other types of extremes.
Jozef Skákala, David Ford, Keith Haines, Amos Lawless, Matthew J. Martin, Philip Browne, Marcin Chrust, Stefano Ciavatta, Alison Fowler, Daniel Lea, Matthew Palmer, Andrea Rochner, Jennifer Waters, Hao Zuo, Deep S. Banerjee, Mike Bell, Davi M. Carneiro, Yumeng Chen, Susan Kay, Dale Partridge, Martin Price, Richard Renshaw, Georgy Shapiro, and James While
Ocean Sci., 21, 1709–1734, https://doi.org/10.5194/os-21-1709-2025, https://doi.org/10.5194/os-21-1709-2025, 2025
Short summary
Short summary
UK marine data assimilation (MDA) involves a closely collaborating research community. In this paper, we offer both an overview of the state of the art and a vision for the future across all of the main areas of UK MDA, ranging from physics to biogeochemistry to coupled DA. We discuss the current UK MDA stakeholder applications, highlight theoretical developments needed to advance our systems, and reflect upon upcoming opportunities with respect to hardware and observational missions.
Yumeng Chen, Lars Nerger, and Amos S. Lawless
EGUsphere, https://doi.org/10.5194/egusphere-2024-1078, https://doi.org/10.5194/egusphere-2024-1078, 2024
Short summary
Short summary
In this paper, we present pyPDAF, a Python interface to the parallel data assimilation framework (PDAF) allowing for coupling with Python-based models. We demonstrate the capability and efficiency of pyPDAF under a coupled data assimilation setup.
Yumeng Chen, Polly Smith, Alberto Carrassi, Ivo Pasmans, Laurent Bertino, Marc Bocquet, Tobias Sebastian Finn, Pierre Rampal, and Véronique Dansereau
The Cryosphere, 18, 2381–2406, https://doi.org/10.5194/tc-18-2381-2024, https://doi.org/10.5194/tc-18-2381-2024, 2024
Short summary
Short summary
We explore multivariate state and parameter estimation using a data assimilation approach through idealised simulations in a dynamics-only sea-ice model based on novel rheology. We identify various potential issues that can arise in complex operational sea-ice models when model parameters are estimated. Even though further investigation will be needed for such complex sea-ice models, we show possibilities of improving the observed and the unobserved model state forecast and parameter accuracy.
Bo Dong, Ross Bannister, Yumeng Chen, Alison Fowler, and Keith Haines
Geosci. Model Dev., 16, 4233–4247, https://doi.org/10.5194/gmd-16-4233-2023, https://doi.org/10.5194/gmd-16-4233-2023, 2023
Short summary
Short summary
Traditional Kalman smoothers are expensive to apply in large global ocean operational forecast and reanalysis systems. We develop a cost-efficient method to overcome the technical constraints and to improve the performance of existing reanalysis products.
Tobias Sebastian Finn, Charlotte Durand, Alban Farchi, Marc Bocquet, Yumeng Chen, Alberto Carrassi, and Véronique Dansereau
The Cryosphere, 17, 2965–2991, https://doi.org/10.5194/tc-17-2965-2023, https://doi.org/10.5194/tc-17-2965-2023, 2023
Short summary
Short summary
We combine deep learning with a regional sea-ice model to correct model errors in the sea-ice dynamics of low-resolution forecasts towards high-resolution simulations. The combined model improves the forecast by up to 75 % and thereby surpasses the performance of persistence. As the error connection can additionally be used to analyse the shortcomings of the forecasts, this study highlights the potential of combined modelling for short-term sea-ice forecasting.
Sukun Cheng, Yumeng Chen, Ali Aydoğdu, Laurent Bertino, Alberto Carrassi, Pierre Rampal, and Christopher K. R. T. Jones
The Cryosphere, 17, 1735–1754, https://doi.org/10.5194/tc-17-1735-2023, https://doi.org/10.5194/tc-17-1735-2023, 2023
Short summary
Short summary
This work studies a novel application of combining a Lagrangian sea ice model, neXtSIM, and data assimilation. It uses a deterministic ensemble Kalman filter to incorporate satellite-observed ice concentration and thickness in simulations. The neXtSIM Lagrangian nature is handled using a remapping strategy on a common homogeneous mesh. The ensemble is formed by perturbing air–ocean boundary conditions and ice cohesion. Thanks to data assimilation, winter Arctic sea ice forecasting is enhanced.
Yumeng Chen, Alberto Carrassi, and Valerio Lucarini
Nonlin. Processes Geophys., 28, 633–649, https://doi.org/10.5194/npg-28-633-2021, https://doi.org/10.5194/npg-28-633-2021, 2021
Short summary
Short summary
Chaotic dynamical systems are sensitive to the initial conditions, which are crucial for climate forecast. These properties are often used to inform the design of data assimilation (DA), a method used to estimate the exact initial conditions. However, obtaining the instability properties is burdensome for complex problems, both numerically and analytically. Here, we suggest a different viewpoint. We show that the skill of DA can be used to infer the instability properties of a dynamical system.
Cited articles
Becker, R. and Rannacher, R.: An optimal control approach to a posteriori error
estimation in finite element methods, Acta Numer., 10, 1–102, 2001. a
Behrens, J.: An adaptive semi-Lagrangian advection scheme and its
parallelization, Mon. Weather Rev., 124, 2386–2395, 1996. a
Behrens, J.: Data Structures for Computational Efficiency, Springer
Berlin Heidelberg, Berlin, Heidelberg, 49–69, https://doi.org/10.1007/3-540-33383-5_4,
2006a. a
Behrens, J.: Adaptive atmospheric modeling: key techniques in grid generation,
data structures, and numerical operations with applications, vol. 207,
Lecture Notes in Computational Science and Engineering, Springer-Verlag Berlin Heidelberg, 2006b. a
Carpenter Jr., R. L., Droegemeier, K. K., Woodward, P. R., and Hane, C. E.:
Application of the piecewise parabolic method (PPM) to meteorological
modeling, Mon. Weather Rev., 118, 586–612, 1990. a
Carslaw, K. S., Boucher, O., Spracklen, D. V., Mann, G. W., Rae, J. G. L., Woodward, S., and Kulmala, M.: A review of natural aerosol interactions and feedbacks within the Earth system, Atmos. Chem. Phys., 10, 1701–1737, https://doi.org/10.5194/acp-10-1701-2010, 2010. a
Chen, Y., Simon, K., and Behrens, J.: yumengch/AMRTransport: Extending Legacy Climate Models by Adaptive Mesh Refinement for Single Component Tracer Transport – GMD (Version 0.01), Zenodo, https://doi.org/10.5281/zenodo.4013277, 2020. a
Colella, P. and Woodward, P. R.: The piecewise parabolic method (PPM) for
gas-dynamical simulations, J. Comput. Phys., 54, 174–201,
https://doi.org/10.1016/0021-9991(84)90143-8, 1984. a, b, c
Eliasen, E., Machenhauer, B., and Rasmussen, E.: On a numerical method for
integration of the hydrodynamical equations with a spectral representation of
the horizontal fields, Report No. 2, Institut for Teoretisk Meteorologi, University of Copenhagen, 1970. a
Ferguson, J. O., Jablonowski, C., Johansen, H., McCorquodale, P., Colella, P.,
and Ullrich, P. A.: Analyzing the adaptive mesh refinement (AMR)
characteristics of a high-order 2D cubed-sphere shallow-water model, Mon. Weather Rev., 144, 4641–4666, 2016. a
Gläser, G., Kerkweg, A., and Wernli, H.: The Mineral Dust Cycle in EMAC 2.40: sensitivity to the spectral resolution and the dust emission scheme, Atmos. Chem. Phys., 12, 1611–1627, https://doi.org/10.5194/acp-12-1611-2012, 2012. a
Herrington, A. R., Lauritzen, P. H., Reed, K. A., Goldhaber, S., and Eaton,
B. E.: Exploring a Lower-Resolution Physics Grid in CAM-SE-CSLAM, J.
Adv. Model. Earth Syst., 11, 1894–1916, 2019. a
Iske, A. and Käser, M.: Conservative semi-Lagrangian advection on adaptive
unstructured meshes, Numer. Meth. Part. D. E., 20, 388–411, 2004. a
Jablonowski, C., Herzog, M., Penner, J. E., Oehmke, R. C., Stout, Q. F.,
Van Leer, B., and Powell, K. G.: Block-structured adaptive grids on the
sphere: Advection experiments, Mon. Weather Rev., 134, 3691–3713, 2006. a
Jöckel, P., von Kuhlmann, R., Lawrence, M. G., Steil, B., Brenninkmeijer,
C. A., Crutzen, P. J., Rasch, P. J., and Eaton, B.: On a fundamental problem
in implementing flux-form advection schemes for tracer transport in
3-dimensional general circulation and chemistry transport models, Q.
J. Roy. Meteor. Soc., 127, 1035–1052, 2001. a, b
Kessler, M.: Development and analysis of an adaptive transport scheme,
Atmos. Environ., 33, 2347–2360, 1999. a
Kopera, M. A. and Giraldo, F. X.: Mass conservation of the unified continuous
and discontinuous element-based Galerkin methods on dynamically adaptive
grids with application to atmospheric simulations, J. Comput. Phys., 297, 90–103, 2015. a
Lauritzen, P. H.: A Stability Analysis of Finite-Volume Advection Schemes
Permitting Long Time Steps, Mon. Weather Rev., 135, 2658–2673,
https://doi.org/10.1175/MWR3425.1, 2007. a
Lauritzen, P. H., Nair, R. D., and Ullrich, P. A.: A conservative
semi-Lagrangian multi-tracer transport scheme (CSLAM) on the cubed-sphere
grid, J. Comput. Phys., 229, 1401–1424, 2010. a
Leonard, B., Lock, A., and Macvean, M.: The nirvana scheme applied to
one-dimensional advection, Int. J. Numer. Methods Heat Fluid Flow, 5,
341–377, https://doi.org/10.1108/EUM0000000004120, 1995. a
Leonard, B., Lock, A., and MacVean, M.: Conservative Explicit
Unrestricted-Time-Step Multidimensional Constancy-Preserving Advection
Schemes, Mon. Weather Rev., 124, 2588–2606,
https://doi.org/10.1175/1520-0493(1996)124<2588:CEUTSM>2.0.CO;2, 1996. a
Lin, S.-J.: A “Vertically Lagrangian” Finite-Volume Dynamical Core for
Global Models, Mon. Weather Rev., 132, 2293–2307,
https://doi.org/10.1175/1520-0493(2004)132<2293:AVLFDC>2.0.CO;2, 2004. a, b, c
Lin, S.-J. and Rood, R. B.: Multidimensional Flux-Form Semi-Lagrangian
Transport Schemes, Mon. Weather Rev., 124, 2046–2070,
https://doi.org/10.1175/1520-0493(1996)124<2046:MFFSLT>2.0.CO;2, 1996. a, b, c, d
Liu, M. and Westphal, D. L.: A study of the sensitivity of simulated mineral
dust production to model resolution, J. Geophys. Res.-Atmos., 106, 18099–18112, 2001. a
Nair, R. D. and Lauritzen, P. H.: A class of deformational flow test cases for
linear transport problems on the sphere, J. Comput. Phys.,
229, 8868–8887, https://doi.org/10.1016/j.jcp.2010.08.014, 2010. a
Rodríguez, S., Alastuey, A., Alonso-Pérez, S., Querol, X., Cuevas, E., Abreu-Afonso, J., Viana, M., Pérez, N., Pandolfi, M., and de la Rosa, J.: Transport of desert dust mixed with North African industrial pollutants in the subtropical Saharan Air Layer, Atmos. Chem. Phys., 11, 6663–6685, https://doi.org/10.5194/acp-11-6663-2011, 2011. a
Schultz, M. G., Stadtler, S., Schröder, S., Taraborrelli, D., Franco, B., Krefting, J., Henrot, A., Ferrachat, S., Lohmann, U., Neubauer, D., Siegenthaler-Le Drian, C., Wahl, S., Kokkola, H., Kühn, T., Rast, S., Schmidt, H., Stier, P., Kinnison, D., Tyndall, G. S., Orlando, J. J., and Wespes, C.: The chemistry–climate model ECHAM6.3-HAM2.3-MOZ1.0, Geosci. Model Dev., 11, 1695–1723, https://doi.org/10.5194/gmd-11-1695-2018, 2018. a
Simmons, A. J. and Burridge, D. M.: An energy and angular-momentum conserving
vertical finite-difference scheme and hybrid vertical coordinates, Mon. Weather Rev., 109, 758–766, 1981. a
Simon, K. and Behrens, J.: Multiscale finite elements through advection-induced
coordinates for transient advection-diffusion equations, arXiv preprint
arXiv:1802.07684, 2018. a
Skamarock, W. C. and Klemp, J. B.: Adaptive grid refinement for two-dimensional
and three-dimensional nonhydrostatic atmospheric flow, Mon. Weather Rev., 121, 788–804, 1993. a
St-Cyr, A., Jablonowski, C., Dennis, J. M., Tufo, H. M., and Thomas, S. J.: A
comparison of two shallow-water models with nonconforming adaptive grids,
Mon. Weather Rev., 136, 1898–1922, 2008. a
Stevens, B., Giorgetta, M., Esch, M., Mauritsen, T., Crueger, T., Rast, S., Salzmann, M., Schmidt, H., Bader, J., Block, K., Brokopf, R., Fast, I., Kinne, S., Kornblueh, L., Lohmann, U., Pincus, R., Reichler, T., and Roeckner, E.: Atmospheric
component of the MPI-M Earth System Model: ECHAM6, J. Adv.
Model. Earth Syst., 5, 146–172, 2013. a, b
Stier, P., Feichter, J., Kinne, S., Kloster, S., Vignati, E., Wilson, J., Ganzeveld, L., Tegen, I., Werner, M., Balkanski, Y., Schulz, M., Boucher, O., Minikin, A., and Petzold, A.: The aerosol-climate model ECHAM5-HAM, Atmos. Chem. Phys., 5, 1125–1156, https://doi.org/10.5194/acp-5-1125-2005, 2005. a
Vignati, E., Wilson, J., and Stier, P.: M7: An efficient size-resolved aerosol
microphysics module for large-scale aerosol transport models, J.
Geophys. Res.-Atmos., 109, D22202, 2004. a
Weller, H., Ringler, T., Piggott, M., and Wood, N.: Challenges facing adaptive
mesh modeling of the atmosphere and ocean, B. Am.
Meteorol. Soc., 91, 105–108, 2010. a
Williamson, D. L., Drake, J. B., Hack, J. J., Jakob, R., and Swarztrauber,
P. N.: A standard test set for numerical approximations to the shallow water
equations in spherical geometry, J. Comput. Phys., 102,
211–224, 1992. a
Short summary
Mesh adaptivity can reduce overall model error by only refining meshes in specific areas where it us necessary in the runtime. Here we suggest a way to integrate mesh adaptivity into an existing Earth system model, ECHAM6, without having to redesign the implementation from scratch. We show that while the additional computational effort is manageable, the error can be reduced compared to a low-resolution standard model using an idealized test and relatively realistic dust transport tests.
Mesh adaptivity can reduce overall model error by only refining meshes in specific areas where...