Articles | Volume 12, issue 3
https://doi.org/10.5194/gmd-12-1165-2019
© Author(s) 2019. 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-12-1165-2019
© Author(s) 2019. This work is distributed under
the Creative Commons Attribution 4.0 License.
the Creative Commons Attribution 4.0 License.
Devito (v3.1.0): an embedded domain-specific language for finite differences and geophysical exploration
Mathias Louboutin
CORRESPONDING AUTHOR
School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, USA
Michael Lange
ECMWF, Reading, UK
Fabio Luporini
Earth Science and Engineering, Imperial College London, London, UK
Navjot Kukreja
Earth Science and Engineering, Imperial College London, London, UK
Philipp A. Witte
School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, USA
Felix J. Herrmann
School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, USA
Paulius Velesko
Earth Science and Engineering, Imperial College London, London, UK
Gerard J. Gorman
Earth Science and Engineering, Imperial College London, London, UK
Related authors
Navjot Kukreja, Jan Hückelheim, Mathias Louboutin, John Washbourne, Paul H. J. Kelly, and Gerard J. Gorman
Geosci. Model Dev., 15, 3815–3829, https://doi.org/10.5194/gmd-15-3815-2022, https://doi.org/10.5194/gmd-15-3815-2022, 2022
Short summary
Short summary
Full waveform inversion (FWI) is a partial-differential equation (PDE)-constrained optimization problem that is notorious for its high computational load and memory footprint. In this paper we present a method that combines recomputation with lossy compression to accelerate the computation with minimal loss of precision in the results. We show this using experiments running FWI with a variety of compression settings on a popular academic dataset.
Navjot Kukreja, Jan Hückelheim, Mathias Louboutin, John Washbourne, Paul H. J. Kelly, and Gerard J. Gorman
Geosci. Model Dev., 15, 3815–3829, https://doi.org/10.5194/gmd-15-3815-2022, https://doi.org/10.5194/gmd-15-3815-2022, 2022
Short summary
Short summary
Full waveform inversion (FWI) is a partial-differential equation (PDE)-constrained optimization problem that is notorious for its high computational load and memory footprint. In this paper we present a method that combines recomputation with lossy compression to accelerate the computation with minimal loss of precision in the results. We show this using experiments running FWI with a variety of compression settings on a popular academic dataset.
Andreas Müller, Willem Deconinck, Christian Kühnlein, Gianmarco Mengaldo, Michael Lange, Nils Wedi, Peter Bauer, Piotr K. Smolarkiewicz, Michail Diamantakis, Sarah-Jane Lock, Mats Hamrud, Sami Saarinen, George Mozdzynski, Daniel Thiemert, Michael Glinton, Pierre Bénard, Fabrice Voitus, Charles Colavolpe, Philippe Marguinaud, Yongjun Zheng, Joris Van Bever, Daan Degrauwe, Geert Smet, Piet Termonia, Kristian P. Nielsen, Bent H. Sass, Jacob W. Poulsen, Per Berg, Carlos Osuna, Oliver Fuhrer, Valentin Clement, Michael Baldauf, Mike Gillard, Joanna Szmelter, Enda O'Brien, Alastair McKinstry, Oisín Robinson, Parijat Shukla, Michael Lysaght, Michał Kulczewski, Milosz Ciznicki, Wojciech Piątek, Sebastian Ciesielski, Marek Błażewicz, Krzysztof Kurowski, Marcin Procyk, Pawel Spychala, Bartosz Bosak, Zbigniew P. Piotrowski, Andrzej Wyszogrodzki, Erwan Raffin, Cyril Mazauric, David Guibert, Louis Douriez, Xavier Vigouroux, Alan Gray, Peter Messmer, Alexander J. Macfaden, and Nick New
Geosci. Model Dev., 12, 4425–4441, https://doi.org/10.5194/gmd-12-4425-2019, https://doi.org/10.5194/gmd-12-4425-2019, 2019
Short summary
Short summary
This paper presents an overview of the ESCAPE project. Dwarfs (key patterns in terms of computation and communication) are identified in weather prediction models. They are optimised for different hardware architectures. New algorithms are developed that are specifically designed for better energy efficiency and improved portability through domain-specific languages. Different numerical techniques are compared in terms of energy efficiency and performance for a variety of computing technologies.
Gheorghe-Teodor Bercea, Andrew T. T. McRae, David A. Ham, Lawrence Mitchell, Florian Rathgeber, Luigi Nardi, Fabio Luporini, and Paul H. J. Kelly
Geosci. Model Dev., 9, 3803–3815, https://doi.org/10.5194/gmd-9-3803-2016, https://doi.org/10.5194/gmd-9-3803-2016, 2016
Short summary
Short summary
Unstructured meshes offer flexibility but are perceived as slow. Some applications, including atmosphere or ocean simulations, admit an extruded mesh: the horizontal mesh may be unstructured, but the vertical dimension can be traversed in a structured way. By extending the Firedrake automated simulation framework to this case, we show that an extruded mesh can be traversed as fast as a structured mesh. This paves the way for highly efficient unstructured mesh models of the ocean and atmosphere.
A. S. Candy, A. Avdis, J. Hill, G. J. Gorman, and M. D. Piggott
Geosci. Model Dev. Discuss., https://doi.org/10.5194/gmdd-7-5993-2014, https://doi.org/10.5194/gmdd-7-5993-2014, 2014
Revised manuscript has not been submitted
Related subject area
Climate and Earth system modeling
INFERNO-peat v1.0.0: a representation of northern high-latitude peat fires in the JULES-INFERNO global fire model
The 4DEnVar-based weakly coupled land data assimilation system for E3SM version 2
Continental-scale bias-corrected climate and hydrological projections for Australia
G6-1.5K-SAI: a new Geoengineering Model Intercomparison Project (GeoMIP) experiment integrating recent advances in solar radiation modification studies
Modeling the effects of tropospheric ozone on the growth and yield of global staple crops with DSSAT v4.8.0
A one-dimensional urban flow model with an eddy-diffusivity mass-flux (EDMF) scheme and refined turbulent transport (MLUCM v3.0)
DCMIP2016: the tropical cyclone test case
Interactions between atmospheric composition and climate change – progress in understanding and future opportunities from AerChemMIP, PDRMIP, and RFMIP
CD-type discretization for sea ice dynamics in FESOM version 2
CSDMS Data Components: data–model integration tools for Earth surface processes modeling
A generic algorithm to automatically classify urban fabric according to the local climate zone system: implementation in GeoClimate 0.0.1 and application to French cities
Modelling water isotopologues (1H2H16O, 1H217O) in the coupled numerical climate model iLOVECLIM (version 1.1.5)
Accurate assessment of land–atmosphere coupling in climate models requires high-frequency data output
Towards variance-conserving reconstructions of climate indices with Gaussian process regression in an embedding space
A diatom extension to the cGEnIE Earth system model – EcoGEnIE 1.1
Carbon isotopes in the marine biogeochemistry model FESOM2.1-REcoM3
Flux coupling approach on an exchange grid for the IOW Earth System Model (version 1.04.00) of the Baltic Sea region
Using EUREC4A/ATOMIC field campaign data to improve trade wind regimes in the Community Atmosphere Model
New model ensemble reveals how forcing uncertainty and model structure alter climate simulated across CMIP generations of the Community Earth System Model
Quantifying wildfire drivers and predictability in boreal peatlands using a two-step error-correcting machine learning framework in TeFire v1.0
Benchmarking GOCART-2G in the Goddard Earth Observing System (GEOS)
Energy-conserving physics for nonhydrostatic dynamics in mass coordinate models
Evaluation and optimisation of the soil carbon turnover routine in the MONICA model (version 3.3.1)
Assessing the sensitivity of aerosol mass budget and effective radiative forcing to horizontal grid spacing in E3SMv1 using a regional refinement approach
Towards the definition of a solar forcing dataset for CMIP7
ibicus: a new open-source Python package and comprehensive interface for statistical bias adjustment and evaluation in climate modelling (v1.0.1)
Disentangling the hydrological and hydraulic controls on streamflow variability in Energy Exascale Earth System Model (E3SM) V2 – a case study in the Pantanal region
Constraining the carbon cycle in JULES-ES-1.0
The utility of simulated ocean chlorophyll observations: a case study with the Chlorophyll Observation Simulator Package (version 1) in CESMv2.2
GeoPDNN 1.0: a semi-supervised deep learning neural network using pseudo-labels for three-dimensional shallow strata modelling and uncertainty analysis in urban areas from borehole data
The prototype NOAA Aerosol Reanalysis version 1.0: description of the modeling system and its evaluation
Performance and process-based evaluation of the BARPA-R Australasian regional climate model version 1
Monsoon Mission Coupled Forecast System version 2.0: model description and Indian monsoon simulations
Exploring the ocean mesoscale at reduced computational cost with FESOM 2.5: efficient modeling strategies applied to the Southern Ocean
Truly conserving with conservative remapping methods
High-resolution downscaling of CMIP6 Earth system and global climate models using deep learning for Iberia
Earth system modeling on modular supercomputing architecture: coupled atmosphere–ocean simulations with ICON 2.6.6-rc
Global Downscaled Projections for Climate Impacts Research (GDPCIR): preserving quantile trends for modeling future climate impacts
Understanding changes in cloud simulations from E3SM version 1 to version 2
WRF (v4.0)–SUEWS (v2018c) coupled system: development, evaluation and application
Scenario setup and forcing data for impact model evaluation and impact attribution within the third round of the Inter-Sectoral Impact Model Intercomparison Project (ISIMIP3a)
Deep learning model based on multi-scale feature fusion for precipitation nowcasting
The Framework for Assessing Changes To Sea-level (FACTS) v1.0: a platform for characterizing parametric and structural uncertainty in future global, relative, and extreme sea-level change
Getting the leaves right matters for estimating temperature extremes
The Southern Ocean Freshwater Input from Antarctica (SOFIA) Initiative: scientific objectives and experimental design
Modeling and evaluating the effects of irrigation on land–atmosphere interaction in southwestern Europe with the regional climate model REMO2020–iMOVE using a newly developed parameterization
The Regional Climate-Chemistry-Ecology Coupling Model RegCM-Chem (v4.6)-YIBs (v1.0): Development and Application
Process-oriented models of autumn leaf phenology: ways to sound calibration and implications of uncertain projections
An evaluation of the LLC4320 global-ocean simulation based on the submesoscale structure of modeled sea surface temperature fields
An emulation-based approach for interrogating reactive transport models
Katie R. Blackford, Matthew Kasoar, Chantelle Burton, Eleanor Burke, Iain Colin Prentice, and Apostolos Voulgarakis
Geosci. Model Dev., 17, 3063–3079, https://doi.org/10.5194/gmd-17-3063-2024, https://doi.org/10.5194/gmd-17-3063-2024, 2024
Short summary
Short summary
Peatlands are globally important stores of carbon which are being increasingly threatened by wildfires with knock-on effects on the climate system. Here we introduce a novel peat fire parameterization in the northern high latitudes to the INFERNO global fire model. Representing peat fires increases annual burnt area across the high latitudes, alongside improvements in how we capture year-to-year variation in burning and emissions.
Pengfei Shi, L. Ruby Leung, Bin Wang, Kai Zhang, Samson M. Hagos, and Shixuan Zhang
Geosci. Model Dev., 17, 3025–3040, https://doi.org/10.5194/gmd-17-3025-2024, https://doi.org/10.5194/gmd-17-3025-2024, 2024
Short summary
Short summary
Improving climate predictions have profound socio-economic impacts. This study introduces a new weakly coupled land data assimilation (WCLDA) system for a coupled climate model. We demonstrate improved simulation of soil moisture and temperature in many global regions and throughout the soil layers. Furthermore, significant improvements are also found in reproducing the time evolution of the 2012 US Midwest drought. The WCLDA system provides the groundwork for future predictability studies.
Justin Peter, Elisabeth Vogel, Wendy Sharples, Ulrike Bende-Michl, Louise Wilson, Pandora Hope, Andrew Dowdy, Greg Kociuba, Sri Srikanthan, Vi Co Duong, Jake Roussis, Vjekoslav Matic, Zaved Khan, Alison Oke, Margot Turner, Stuart Baron-Hay, Fiona Johnson, Raj Mehrotra, Ashish Sharma, Marcus Thatcher, Ali Azarvinand, Steven Thomas, Ghyslaine Boschat, Chantal Donnelly, and Robert Argent
Geosci. Model Dev., 17, 2755–2781, https://doi.org/10.5194/gmd-17-2755-2024, https://doi.org/10.5194/gmd-17-2755-2024, 2024
Short summary
Short summary
We detail the production of datasets and communication to end users of high-resolution projections of rainfall, runoff, and soil moisture for the entire Australian continent. This is important as previous projections for Australia were for small regions and used differing techniques for their projections, making comparisons difficult across Australia's varied climate zones. The data will be beneficial for research purposes and to aid adaptation to climate change.
Daniele Visioni, Alan Robock, Jim Haywood, Matthew Henry, Simone Tilmes, Douglas G. MacMartin, Ben Kravitz, Sarah J. Doherty, John Moore, Chris Lennard, Shingo Watanabe, Helene Muri, Ulrike Niemeier, Olivier Boucher, Abu Syed, Temitope S. Egbebiyi, Roland Séférian, and Ilaria Quaglia
Geosci. Model Dev., 17, 2583–2596, https://doi.org/10.5194/gmd-17-2583-2024, https://doi.org/10.5194/gmd-17-2583-2024, 2024
Short summary
Short summary
This paper describes a new experimental protocol for the Geoengineering Model Intercomparison Project (GeoMIP). In it, we describe the details of a new simulation of sunlight reflection using the stratospheric aerosols that climate models are supposed to run, and we explain the reasons behind each choice we made when defining the protocol.
Jose Rafael Guarin, Jonas Jägermeyr, Elizabeth A. Ainsworth, Fabio A. A. Oliveira, Senthold Asseng, Kenneth Boote, Joshua Elliott, Lisa Emberson, Ian Foster, Gerrit Hoogenboom, David Kelly, Alex C. Ruane, and Katrina Sharps
Geosci. Model Dev., 17, 2547–2567, https://doi.org/10.5194/gmd-17-2547-2024, https://doi.org/10.5194/gmd-17-2547-2024, 2024
Short summary
Short summary
The effects of ozone (O3) stress on crop photosynthesis and leaf senescence were added to maize, rice, soybean, and wheat crop models. The modified models reproduced growth and yields under different O3 levels measured in field experiments and reported in the literature. The combined interactions between O3 and additional stresses were reproduced with the new models. These updated crop models can be used to simulate impacts of O3 stress under future climate change and air pollution scenarios.
Jiachen Lu, Negin Nazarian, Melissa Anne Hart, E. Scott Krayenhoff, and Alberto Martilli
Geosci. Model Dev., 17, 2525–2545, https://doi.org/10.5194/gmd-17-2525-2024, https://doi.org/10.5194/gmd-17-2525-2024, 2024
Short summary
Short summary
This study enhances urban canopy models by refining key assumptions. Simulations for various urban scenarios indicate discrepancies in turbulent transport efficiency for flow properties. We propose two modifications that involve characterizing diffusion coefficients for momentum and turbulent kinetic energy separately and introducing a physics-based
mass-fluxterm. These adjustments enhance the model's performance, offering more reliable temperature and surface flux estimates.
Justin L. Willson, Kevin A. Reed, Christiane Jablonowski, James Kent, Peter H. Lauritzen, Ramachandran Nair, Mark A. Taylor, Paul A. Ullrich, Colin M. Zarzycki, David M. Hall, Don Dazlich, Ross Heikes, Celal Konor, David Randall, Thomas Dubos, Yann Meurdesoif, Xi Chen, Lucas Harris, Christian Kühnlein, Vivian Lee, Abdessamad Qaddouri, Claude Girard, Marco Giorgetta, Daniel Reinert, Hiroaki Miura, Tomoki Ohno, and Ryuji Yoshida
Geosci. Model Dev., 17, 2493–2507, https://doi.org/10.5194/gmd-17-2493-2024, https://doi.org/10.5194/gmd-17-2493-2024, 2024
Short summary
Short summary
Accurate simulation of tropical cyclones (TCs) is essential to understanding their behavior in a changing climate. One way this is accomplished is through model intercomparison projects, where results from multiple climate models are analyzed to provide benchmark solutions for the wider climate modeling community. This study describes and analyzes the previously developed TC test case for nine climate models in an intercomparison project, providing solutions that aid in model development.
Stephanie Fiedler, Vaishali Naik, Fiona M. O'Connor, Christopher J. Smith, Paul Griffiths, Ryan J. Kramer, Toshihiko Takemura, Robert J. Allen, Ulas Im, Matthew Kasoar, Angshuman Modak, Steven Turnock, Apostolos Voulgarakis, Duncan Watson-Parris, Daniel M. Westervelt, Laura J. Wilcox, Alcide Zhao, William J. Collins, Michael Schulz, Gunnar Myhre, and Piers M. Forster
Geosci. Model Dev., 17, 2387–2417, https://doi.org/10.5194/gmd-17-2387-2024, https://doi.org/10.5194/gmd-17-2387-2024, 2024
Short summary
Short summary
Climate scientists want to better understand modern climate change. Thus, climate model experiments are performed and compared. The results of climate model experiments differ, as assessed in the latest Intergovernmental Panel on Climate Change (IPCC) assessment report. This article gives insights into the challenges and outlines opportunities for further improving the understanding of climate change. It is based on views of a group of experts in atmospheric composition–climate interactions.
Sergey Danilov, Carolin Mehlmann, Dmitry Sidorenko, and Qiang Wang
Geosci. Model Dev., 17, 2287–2297, https://doi.org/10.5194/gmd-17-2287-2024, https://doi.org/10.5194/gmd-17-2287-2024, 2024
Short summary
Short summary
Sea ice models are a necessary component of climate models. At very high resolution they are capable of simulating linear kinematic features, such as leads, which are important for better prediction of heat exchanges between the ocean and atmosphere. Two new discretizations are described which improve the sea ice component of the Finite volumE Sea ice–Ocean Model (FESOM version 2) by allowing simulations of finer scales.
Tian Gan, Gregory E. Tucker, Eric W. H. Hutton, Mark D. Piper, Irina Overeem, Albert J. Kettner, Benjamin Campforts, Julia M. Moriarty, Brianna Undzis, Ethan Pierce, and Lynn McCready
Geosci. Model Dev., 17, 2165–2185, https://doi.org/10.5194/gmd-17-2165-2024, https://doi.org/10.5194/gmd-17-2165-2024, 2024
Short summary
Short summary
This study presents the design, implementation, and application of the CSDMS Data Components. The case studies demonstrate that the Data Components provide a consistent way to access heterogeneous datasets from multiple sources, and to seamlessly integrate them with various models for Earth surface process modeling. The Data Components support the creation of open data–model integration workflows to improve the research transparency and reproducibility.
Jérémy Bernard, Erwan Bocher, Matthieu Gousseff, François Leconte, and Elisabeth Le Saux Wiederhold
Geosci. Model Dev., 17, 2077–2116, https://doi.org/10.5194/gmd-17-2077-2024, https://doi.org/10.5194/gmd-17-2077-2024, 2024
Short summary
Short summary
Geographical features may have a considerable effect on local climate. The local climate zone (LCZ) system proposed by Stewart and Oke (2012) is seen as a standard approach for classifying any zone according to a set of geographic indicators. While many methods already exist to map the LCZ, only a few tools are openly and freely available. We present the algorithm implemented in GeoClimate software to identify the LCZ of any place in the world using OpenStreetMap data.
Thomas Extier, Thibaut Caley, and Didier M. Roche
Geosci. Model Dev., 17, 2117–2139, https://doi.org/10.5194/gmd-17-2117-2024, https://doi.org/10.5194/gmd-17-2117-2024, 2024
Short summary
Short summary
Stable water isotopes are used to infer changes in the hydrological cycle for different time periods in climatic archive and climate models. We present the implementation of the δ2H and δ17O water isotopes in the coupled climate model iLOVECLIM and calculate the d- and 17O-excess. Results of a simulation under preindustrial conditions show that the model correctly reproduces the water isotope distribution in the atmosphere and ocean in comparison to data and other global circulation models.
Kirsten L. Findell, Zun Yin, Eunkyo Seo, Paul A. Dirmeyer, Nathan P. Arnold, Nathaniel Chaney, Megan D. Fowler, Meng Huang, David M. Lawrence, Po-Lun Ma, and Joseph A. Santanello Jr.
Geosci. Model Dev., 17, 1869–1883, https://doi.org/10.5194/gmd-17-1869-2024, https://doi.org/10.5194/gmd-17-1869-2024, 2024
Short summary
Short summary
We outline a request for sub-daily data to accurately capture the process-level connections between land states, surface fluxes, and the boundary layer response. This high-frequency model output will allow for more direct comparison with observational field campaigns on process-relevant timescales, enable demonstration of inter-model spread in land–atmosphere coupling processes, and aid in targeted identification of sources of deficiencies and opportunities for improvement of the models.
Marlene Klockmann, Udo von Toussaint, and Eduardo Zorita
Geosci. Model Dev., 17, 1765–1787, https://doi.org/10.5194/gmd-17-1765-2024, https://doi.org/10.5194/gmd-17-1765-2024, 2024
Short summary
Short summary
Reconstructions of climate variability before the observational period rely on climate proxies and sophisticated statistical models to link the proxy information and climate variability. Existing models tend to underestimate the true magnitude of variability, especially if the proxies contain non-climatic noise. We present and test a promising new framework for climate-index reconstructions, based on Gaussian processes, which reconstructs robust variability estimates from noisy and sparse data.
Aaron A. Naidoo-Bagwell, Fanny M. Monteiro, Katharine R. Hendry, Scott Burgan, Jamie D. Wilson, Ben A. Ward, Andy Ridgwell, and Daniel J. Conley
Geosci. Model Dev., 17, 1729–1748, https://doi.org/10.5194/gmd-17-1729-2024, https://doi.org/10.5194/gmd-17-1729-2024, 2024
Short summary
Short summary
As an extension to the EcoGEnIE 1.0 Earth system model that features a diverse plankton community, EcoGEnIE 1.1 includes siliceous plankton diatoms and also considers their impact on biogeochemical cycles. With updates to existing nutrient cycles and the introduction of the silicon cycle, we see improved model performance relative to observational data. Through a more functionally diverse plankton community, the new model enables more comprehensive future study of ocean ecology.
Martin Butzin, Ying Ye, Christoph Völker, Özgür Gürses, Judith Hauck, and Peter Köhler
Geosci. Model Dev., 17, 1709–1727, https://doi.org/10.5194/gmd-17-1709-2024, https://doi.org/10.5194/gmd-17-1709-2024, 2024
Short summary
Short summary
In this paper we describe the implementation of the carbon isotopes 13C and 14C into the marine biogeochemistry model FESOM2.1-REcoM3 and present results of long-term test simulations. Our model results are largely consistent with marine carbon isotope reconstructions for the pre-anthropogenic period, but also exhibit some discrepancies.
Sven Karsten, Hagen Radtke, Matthias Gröger, Ha T. M. Ho-Hagemann, Hossein Mashayekh, Thomas Neumann, and H. E. Markus Meier
Geosci. Model Dev., 17, 1689–1708, https://doi.org/10.5194/gmd-17-1689-2024, https://doi.org/10.5194/gmd-17-1689-2024, 2024
Short summary
Short summary
This paper describes the development of a regional Earth System Model for the Baltic Sea region. In contrast to conventional coupling approaches, the presented model includes a flux calculator operating on a common exchange grid. This approach automatically ensures a locally consistent treatment of fluxes and simplifies the exchange of model components. The presented model can be used for various scientific questions, such as studies of natural variability and ocean–atmosphere interactions.
Skyler Graap and Colin M. Zarzycki
Geosci. Model Dev., 17, 1627–1650, https://doi.org/10.5194/gmd-17-1627-2024, https://doi.org/10.5194/gmd-17-1627-2024, 2024
Short summary
Short summary
A key target for improving climate models is how low, bright clouds are predicted over tropical oceans, since they have important consequences for the Earth's energy budget. A climate model has been updated to improve the physical realism of the treatment of how momentum is moved up and down in the atmosphere. By comparing this updated model to real-world observations from balloon launches, it can be shown to more accurately depict atmospheric structure in trade-wind areas close to the Equator.
Marika M. Holland, Cecile Hannay, John Fasullo, Alexandra Jahn, Jennifer E. Kay, Michael Mills, Isla R. Simpson, William Wieder, Peter Lawrence, Erik Kluzek, and David Bailey
Geosci. Model Dev., 17, 1585–1602, https://doi.org/10.5194/gmd-17-1585-2024, https://doi.org/10.5194/gmd-17-1585-2024, 2024
Short summary
Short summary
Climate evolves in response to changing forcings, as prescribed in simulations. Models and forcings are updated over time to reflect new understanding. This makes it difficult to attribute simulation differences to either model or forcing changes. Here we present new simulations which enable the separation of model structure and forcing influence between two widely used simulation sets. Results indicate a strong influence of aerosol emission uncertainty on historical climate.
Rongyun Tang, Mingzhou Jin, Jiafu Mao, Daniel M. Ricciuto, Anping Chen, and Yulong Zhang
Geosci. Model Dev., 17, 1525–1542, https://doi.org/10.5194/gmd-17-1525-2024, https://doi.org/10.5194/gmd-17-1525-2024, 2024
Short summary
Short summary
Carbon-rich boreal peatlands are at risk of burning. The reproducibility and predictability of rare peatland fire events are investigated by constructing a two-step error-correcting machine learning framework to tackle such complex systems. Fire occurrence and impacts are highly predictable with our approach. Factor-controlling simulations revealed that temperature, moisture, and freeze–thaw cycles control boreal peatland fires, indicating thermal impacts on causing peat fires.
Allison B. Collow, Peter R. Colarco, Arlindo M. da Silva, Virginie Buchard, Huisheng Bian, Mian Chin, Sampa Das, Ravi Govindaraju, Dongchul Kim, and Valentina Aquila
Geosci. Model Dev., 17, 1443–1468, https://doi.org/10.5194/gmd-17-1443-2024, https://doi.org/10.5194/gmd-17-1443-2024, 2024
Short summary
Short summary
The GOCART aerosol module within the Goddard Earth Observing System recently underwent a major refactoring and update to the representation of physical processes. Code changes that were included in GOCART Second Generation (GOCART-2G) are documented, and we establish a benchmark simulation that is to be used for future development of the system. The 4-year benchmark simulation was evaluated using in situ and spaceborne measurements to develop a baseline and prioritize future development.
Oksana Guba, Mark A. Taylor, Peter A. Bosler, Christopher Eldred, and Peter H. Lauritzen
Geosci. Model Dev., 17, 1429–1442, https://doi.org/10.5194/gmd-17-1429-2024, https://doi.org/10.5194/gmd-17-1429-2024, 2024
Short summary
Short summary
We want to reduce errors in the moist energy budget in numerical atmospheric models. We study a few common assumptions and mechanisms that are used for the moist physics. Some mechanisms are more consistent with the underlying equations. Separately, we study how assumptions about models' thermodynamics affect the modeled energy of precipitation. We also explain how to conserve energy in the moist physics for nonhydrostatic models.
Konstantin Aiteew, Jarno Rouhiainen, Claas Nendel, and René Dechow
Geosci. Model Dev., 17, 1349–1385, https://doi.org/10.5194/gmd-17-1349-2024, https://doi.org/10.5194/gmd-17-1349-2024, 2024
Short summary
Short summary
This study evaluated the biogeochemical model MONICA and its performance in simulating soil organic carbon changes. MONICA can reproduce plant growth, carbon and nitrogen dynamics, soil water and temperature. The model results were compared with five established carbon turnover models. With the exception of certain sites, adequate reproduction of soil organic carbon stock change rates was achieved. The MONICA model was capable of performing similar to or even better than the other models.
Jianfeng Li, Kai Zhang, Taufiq Hassan, Shixuan Zhang, Po-Lun Ma, Balwinder Singh, Qiyang Yan, and Huilin Huang
Geosci. Model Dev., 17, 1327–1347, https://doi.org/10.5194/gmd-17-1327-2024, https://doi.org/10.5194/gmd-17-1327-2024, 2024
Short summary
Short summary
By comparing E3SM simulations with and without regional refinement, we find that model horizontal grid spacing considerably affects the simulated aerosol mass budget, aerosol–cloud interactions, and the effective radiative forcing of anthropogenic aerosols. The study identifies the critical physical processes strongly influenced by model resolution. It also highlights the benefit of applying regional refinement in future modeling studies at higher or even convection-permitting resolutions.
Bernd Funke, Thierry Dudok de Wit, Ilaria Ermolli, Margit Haberreiter, Doug Kinnison, Daniel Marsh, Hilde Nesse, Annika Seppälä, Miriam Sinnhuber, and Ilya Usoskin
Geosci. Model Dev., 17, 1217–1227, https://doi.org/10.5194/gmd-17-1217-2024, https://doi.org/10.5194/gmd-17-1217-2024, 2024
Short summary
Short summary
We outline a road map for the preparation of a solar forcing dataset for the upcoming Phase 7 of the Coupled Model Intercomparison Project (CMIP7), considering the latest scientific advances made in the reconstruction of solar forcing and in the understanding of climate response while also addressing the issues that were raised during CMIP6.
Fiona Raphaela Spuler, Jakob Benjamin Wessel, Edward Comyn-Platt, James Varndell, and Chiara Cagnazzo
Geosci. Model Dev., 17, 1249–1269, https://doi.org/10.5194/gmd-17-1249-2024, https://doi.org/10.5194/gmd-17-1249-2024, 2024
Short summary
Short summary
Before using climate models to study the impacts of climate change, bias adjustment is commonly applied to the models to ensure that they correspond with observations at a local scale. However, this can introduce undesirable distortions into the climate model. In this paper, we present an open-source python package called ibicus to enable the comparison and detailed evaluation of bias adjustment methods, facilitating their transparent and rigorous application.
Donghui Xu, Gautam Bisht, Zeli Tan, Chang Liao, Tian Zhou, Hong-Yi Li, and L. Ruby Leung
Geosci. Model Dev., 17, 1197–1215, https://doi.org/10.5194/gmd-17-1197-2024, https://doi.org/10.5194/gmd-17-1197-2024, 2024
Short summary
Short summary
We aim to disentangle the hydrological and hydraulic controls on streamflow variability in a fully coupled earth system model. We found that calibrating only one process (i.e., traditional calibration procedure) will result in unrealistic parameter values and poor performance of the water cycle, while the simulated streamflow is improved. To address this issue, we further proposed a two-step calibration procedure to reconcile the impacts from hydrological and hydraulic processes on streamflow.
Douglas McNeall, Eddy Robertson, and Andy Wiltshire
Geosci. Model Dev., 17, 1059–1089, https://doi.org/10.5194/gmd-17-1059-2024, https://doi.org/10.5194/gmd-17-1059-2024, 2024
Short summary
Short summary
We can run simulations of the land surface and carbon cycle, using computer models to help us understand and predict climate change and its impacts. These simulations are not perfect reproductions of the real land surface, and that can make them less effective tools. We use new statistical and computational techniques to help us understand how different our models are from the real land surface, how to make them more realistic, and how well we can simulate past and future climate.
Genevieve L. Clow, Nicole S. Lovenduski, Michael N. Levy, Keith Lindsay, and Jennifer E. Kay
Geosci. Model Dev., 17, 975–995, https://doi.org/10.5194/gmd-17-975-2024, https://doi.org/10.5194/gmd-17-975-2024, 2024
Short summary
Short summary
Satellite observations of chlorophyll allow us to study marine phytoplankton on a global scale; yet some of these observations are missing due to clouds and other issues. To investigate the impact of missing data, we developed a satellite simulator for chlorophyll in an Earth system model. We found that missing data can impact the global mean chlorophyll by nearly 20 %. The simulated observations provide a more direct comparison to real-world data and can be used to improve model validation.
Jiateng Guo, Xuechuang Xu, Luyuan Wang, Xulei Wang, Lixin Wu, Mark Jessell, Vitaliy Ogarko, Zhibin Liu, and Yufei Zheng
Geosci. Model Dev., 17, 957–973, https://doi.org/10.5194/gmd-17-957-2024, https://doi.org/10.5194/gmd-17-957-2024, 2024
Short summary
Short summary
This study proposes a semi-supervised learning algorithm using pseudo-labels for 3D geological modelling. We establish a 3D geological model using borehole data from a complex real urban local survey area in Shenyang and make an uncertainty analysis of this model. The method effectively expands the sample space, which is suitable for geomodelling and uncertainty analysis from boreholes. The modelling results perform well in terms of spatial morphology and geological semantics.
Shih-Wei Wei, Mariusz Pagowski, Arlindo da Silva, Cheng-Hsuan Lu, and Bo Huang
Geosci. Model Dev., 17, 795–813, https://doi.org/10.5194/gmd-17-795-2024, https://doi.org/10.5194/gmd-17-795-2024, 2024
Short summary
Short summary
This study describes the modeling system and the evaluation results for the first prototype version of a global aerosol reanalysis product at NOAA, prototype NOAA Aerosol ReAnalysis version 1.0 (pNARA v1.0). We evaluated pNARA v1.0 against independent datasets and compared it with other reanalyses. We identified deficiencies in the system (both in the forecast model and in the data assimilation system) and the uncertainties that exist in our reanalysis.
Emma Howard, Chun-Hsu Su, Christian Stassen, Rajashree Naha, Harvey Ye, Acacia Pepler, Samuel S. Bell, Andrew J. Dowdy, Simon O. Tucker, and Charmaine Franklin
Geosci. Model Dev., 17, 731–757, https://doi.org/10.5194/gmd-17-731-2024, https://doi.org/10.5194/gmd-17-731-2024, 2024
Short summary
Short summary
The BARPA-R modelling configuration has been developed to produce high-resolution climate hazard projections within the Australian region. When using boundary driving data from quasi-observed historical conditions, BARPA-R shows good performance with errors generally on par with reanalysis products. BARPA-R also captures trends, known modes of climate variability, large-scale weather processes, and multivariate relationships.
Deepeshkumar Jain, Suryachandra A. Rao, Ramu A. Dandi, Prasanth A. Pillai, Ankur Srivastava, Maheswar Pradhan, and Kiran V. Gangadharan
Geosci. Model Dev., 17, 709–729, https://doi.org/10.5194/gmd-17-709-2024, https://doi.org/10.5194/gmd-17-709-2024, 2024
Short summary
Short summary
The present paper discusses and evaluates the new Monsoon Mission Coupled Forecast System model (MMCFS) version 2.0 which upgrades the currently operational MMCFS v1.0 at the Indian Meteorological Department, India. The individual model components have been substantially upgraded independently by their respective scientific groups. MMCFS v2.0 includes these upgrades in the operational coupled model. The new model shows significant skill improvement in simulating the Indian monsoon.
Nathan Beech, Thomas Rackow, Tido Semmler, and Thomas Jung
Geosci. Model Dev., 17, 529–543, https://doi.org/10.5194/gmd-17-529-2024, https://doi.org/10.5194/gmd-17-529-2024, 2024
Short summary
Short summary
Cost-reducing modeling strategies are applied to high-resolution simulations of the Southern Ocean in a changing climate. They are evaluated with respect to observations and traditional, lower-resolution modeling methods. The simulations effectively reproduce small-scale ocean flows seen in satellite data and are largely consistent with traditional model simulations after 4 °C of warming. Small-scale flows are found to intensify near bathymetric features and to become more variable.
Karl E. Taylor
Geosci. Model Dev., 17, 415–430, https://doi.org/10.5194/gmd-17-415-2024, https://doi.org/10.5194/gmd-17-415-2024, 2024
Short summary
Short summary
Remapping gridded data in a way that preserves the conservative properties of the climate system can be essential in coupling model components and for accurate assessment of the system’s energy and mass constituents. Remapping packages capable of handling a wide variety of grids can, for some common grids, calculate remapping weights that are somewhat inaccurate. Correcting for these errors, guidelines are provided to ensure conservation when the weights are used in practice.
Pedro M. M. Soares, Frederico Johannsen, Daniela C. A. Lima, Gil Lemos, Virgílio A. Bento, and Angelina Bushenkova
Geosci. Model Dev., 17, 229–259, https://doi.org/10.5194/gmd-17-229-2024, https://doi.org/10.5194/gmd-17-229-2024, 2024
Short summary
Short summary
This study uses deep learning (DL) to downscale global climate models for the Iberian Peninsula. Four DL architectures were evaluated and trained using historical climate data and then used to downscale future projections from the global models. These show agreement with the original models and reveal a warming of 2 ºC to 6 ºC, along with decreasing precipitation in western Iberia after 2040. This approach offers key regional climate change information for adaptation strategies in the region.
Abhiraj Bishnoi, Olaf Stein, Catrin I. Meyer, René Redler, Norbert Eicker, Helmuth Haak, Lars Hoffmann, Daniel Klocke, Luis Kornblueh, and Estela Suarez
Geosci. Model Dev., 17, 261–273, https://doi.org/10.5194/gmd-17-261-2024, https://doi.org/10.5194/gmd-17-261-2024, 2024
Short summary
Short summary
We enabled the weather and climate model ICON to run in a high-resolution coupled atmosphere–ocean setup on the JUWELS supercomputer, where the ocean and the model I/O runs on the CPU Cluster, while the atmosphere is running simultaneously on GPUs. Compared to a simulation performed on CPUs only, our approach reduces energy consumption by 45 % with comparable runtimes. The experiments serve as preparation for efficient computing of kilometer-scale climate models on future supercomputing systems.
Diana R. Gergel, Steven B. Malevich, Kelly E. McCusker, Emile Tenezakis, Michael T. Delgado, Meredith A. Fish, and Robert E. Kopp
Geosci. Model Dev., 17, 191–227, https://doi.org/10.5194/gmd-17-191-2024, https://doi.org/10.5194/gmd-17-191-2024, 2024
Short summary
Short summary
The freely available Global Downscaled Projections for Climate Impacts Research (GDPCIR) dataset gives researchers a new tool for studying how future climate will evolve at a local or regional level, corresponding to the latest global climate model simulations prepared as part of the UN Intergovernmental Panel on Climate Change’s Sixth Assessment Report. Those simulations represent an enormous advance in quality, detail, and scope that GDPCIR translates to the local level.
Yuying Zhang, Shaocheng Xie, Yi Qin, Wuyin Lin, Jean-Christophe Golaz, Xue Zheng, Po-Lun Ma, Yun Qian, Qi Tang, Christopher R. Terai, and Meng Zhang
Geosci. Model Dev., 17, 169–189, https://doi.org/10.5194/gmd-17-169-2024, https://doi.org/10.5194/gmd-17-169-2024, 2024
Short summary
Short summary
We performed systematic evaluation of clouds simulated in the Energy
Exascale Earth System Model (E3SMv2) to document model performance and understand what updates in E3SMv2 have caused changes in clouds from E3SMv1 to E3SMv2. We find that stratocumulus clouds along the subtropical west coast of continents are dramatically improved, primarily due to the retuning done in CLUBB. This study offers additional insights into clouds simulated in E3SMv2 and will benefit future E3SM developments.
Exascale Earth System Model (E3SMv2) to document model performance and understand what updates in E3SMv2 have caused changes in clouds from E3SMv1 to E3SMv2. We find that stratocumulus clouds along the subtropical west coast of continents are dramatically improved, primarily due to the retuning done in CLUBB. This study offers additional insights into clouds simulated in E3SMv2 and will benefit future E3SM developments.
Ting Sun, Hamidreza Omidvar, Zhenkun Li, Ning Zhang, Wenjuan Huang, Simone Kotthaus, Helen C. Ward, Zhiwen Luo, and Sue Grimmond
Geosci. Model Dev., 17, 91–116, https://doi.org/10.5194/gmd-17-91-2024, https://doi.org/10.5194/gmd-17-91-2024, 2024
Short summary
Short summary
For the first time, we coupled a state-of-the-art urban land surface model – Surface Urban Energy and Water Scheme (SUEWS) – with the widely-used Weather Research and Forecasting (WRF) model, creating an open-source tool that may benefit multiple applications. We tested our new system at two UK sites and demonstrated its potential by examining how human activities in various areas of Greater London influence local weather conditions.
Katja Frieler, Jan Volkholz, Stefan Lange, Jacob Schewe, Matthias Mengel, María del Rocío Rivas López, Christian Otto, Christopher P. O. Reyer, Dirk Nikolaus Karger, Johanna T. Malle, Simon Treu, Christoph Menz, Julia L. Blanchard, Cheryl S. Harrison, Colleen M. Petrik, Tyler D. Eddy, Kelly Ortega-Cisneros, Camilla Novaglio, Yannick Rousseau, Reg A. Watson, Charles Stock, Xiao Liu, Ryan Heneghan, Derek Tittensor, Olivier Maury, Matthias Büchner, Thomas Vogt, Tingting Wang, Fubao Sun, Inga J. Sauer, Johannes Koch, Inne Vanderkelen, Jonas Jägermeyr, Christoph Müller, Sam Rabin, Jochen Klar, Iliusi D. Vega del Valle, Gitta Lasslop, Sarah Chadburn, Eleanor Burke, Angela Gallego-Sala, Noah Smith, Jinfeng Chang, Stijn Hantson, Chantelle Burton, Anne Gädeke, Fang Li, Simon N. Gosling, Hannes Müller Schmied, Fred Hattermann, Jida Wang, Fangfang Yao, Thomas Hickler, Rafael Marcé, Don Pierson, Wim Thiery, Daniel Mercado-Bettín, Robert Ladwig, Ana Isabel Ayala-Zamora, Matthew Forrest, and Michel Bechtold
Geosci. Model Dev., 17, 1–51, https://doi.org/10.5194/gmd-17-1-2024, https://doi.org/10.5194/gmd-17-1-2024, 2024
Short summary
Short summary
Our paper provides an overview of all observational climate-related and socioeconomic forcing data used as input for the impact model evaluation and impact attribution experiments within the third round of the Inter-Sectoral Impact Model Intercomparison Project. The experiments are designed to test our understanding of observed changes in natural and human systems and to quantify to what degree these changes have already been induced by climate change.
Jinkai Tan, Qiqiao Huang, and Sheng Chen
Geosci. Model Dev., 17, 53–69, https://doi.org/10.5194/gmd-17-53-2024, https://doi.org/10.5194/gmd-17-53-2024, 2024
Short summary
Short summary
This study presents a deep learning architecture, multi-scale feature fusion (MFF), to improve the forecast skills of precipitations especially for heavy precipitations. MFF uses multi-scale receptive fields so that the movement features of precipitation systems are well captured. MFF uses the mechanism of discrete probability to reduce uncertainties and forecast errors so that heavy precipitations are produced.
Robert E. Kopp, Gregory G. Garner, Tim H. J. Hermans, Shantenu Jha, Praveen Kumar, Alexander Reedy, Aimée B. A. Slangen, Matteo Turilli, Tamsin L. Edwards, Jonathan M. Gregory, George Koubbe, Anders Levermann, Andre Merzky, Sophie Nowicki, Matthew D. Palmer, and Chris Smith
Geosci. Model Dev., 16, 7461–7489, https://doi.org/10.5194/gmd-16-7461-2023, https://doi.org/10.5194/gmd-16-7461-2023, 2023
Short summary
Short summary
Future sea-level rise projections exhibit multiple forms of uncertainty, all of which must be considered by scientific assessments intended to inform decision-making. The Framework for Assessing Changes To Sea-level (FACTS) is a new software package intended to support assessments of global mean, regional, and extreme sea-level rise. An early version of FACTS supported the development of the IPCC Sixth Assessment Report sea-level projections.
Gregory Duveiller, Mark Pickering, Joaquin Muñoz-Sabater, Luca Caporaso, Souhail Boussetta, Gianpaolo Balsamo, and Alessandro Cescatti
Geosci. Model Dev., 16, 7357–7373, https://doi.org/10.5194/gmd-16-7357-2023, https://doi.org/10.5194/gmd-16-7357-2023, 2023
Short summary
Short summary
Some of our best tools to describe the state of the land system, including the intensity of heat waves, have a problem. The model currently assumes that the number of leaves in ecosystems always follows the same cycle. By using satellite observations of when leaves are present, we show that capturing the yearly changes in this cycle is important to avoid errors in estimating surface temperature. We show that this has strong implications for our capacity to describe heat waves across Europe.
Neil C. Swart, Torge Martin, Rebecca Beadling, Jia-Jia Chen, Christopher Danek, Matthew H. England, Riccardo Farneti, Stephen M. Griffies, Tore Hattermann, Judith Hauck, F. Alexander Haumann, André Jüling, Qian Li, John Marshall, Morven Muilwijk, Andrew G. Pauling, Ariaan Purich, Inga J. Smith, and Max Thomas
Geosci. Model Dev., 16, 7289–7309, https://doi.org/10.5194/gmd-16-7289-2023, https://doi.org/10.5194/gmd-16-7289-2023, 2023
Short summary
Short summary
Current climate models typically do not include full representation of ice sheets. As the climate warms and the ice sheets melt, they add freshwater to the ocean. This freshwater can influence climate change, for example by causing more sea ice to form. In this paper we propose a set of experiments to test the influence of this missing meltwater from Antarctica using multiple different climate models.
Christina Asmus, Peter Hoffmann, Joni-Pekka Pietikäinen, Jürgen Böhner, and Diana Rechid
Geosci. Model Dev., 16, 7311–7337, https://doi.org/10.5194/gmd-16-7311-2023, https://doi.org/10.5194/gmd-16-7311-2023, 2023
Short summary
Short summary
Irrigation modifies the land surface and soil conditions. The effects can be quantified using numerical climate models. Our study introduces a new irrigation parameterization, which simulates the effects of irrigation on land, atmosphere, and vegetation. We applied the parameterization and evaluated the results in terms of their physical consistency. We found an improvement in the model results in the 2 m temperature representation in comparison with observational data for our study.
Nanhong Xie, Tijian Wang, Xiaodong Xie, Xu Yue, Filippo Giorgi, Qian Zhang, Danyang Ma, Rong Song, Baiyao Xu, Shu Li, Bingliang Zhuang, Mengmeng Li, Min Xie, Natalya Andreeva Kilifarska, Georgi Gadzhev, and Reneta Dimitrova
EGUsphere, https://doi.org/10.5194/egusphere-2023-1733, https://doi.org/10.5194/egusphere-2023-1733, 2023
Short summary
Short summary
For the first time, we coupled a regional climate chemistry model RegCM-Chem with a dynamic vegetation model YIBs to create a regional climate-chemistry-ecology model RegCM-Chem-YIBs. We applied it to simulate climatic, chemical and ecological parameters in East Asia and fully validated it on a variety of observational data. The research results show that RegCM-Chem-YIBs model is a valuable tool for studying terrestrial carbon cycle, atmospheric chemistry, and climate change in regional scale.
Michael Meier and Christof Bigler
Geosci. Model Dev., 16, 7171–7201, https://doi.org/10.5194/gmd-16-7171-2023, https://doi.org/10.5194/gmd-16-7171-2023, 2023
Short summary
Short summary
We analyzed >2.3 million calibrations and 39 million projections of leaf coloration models, considering 21 models, 5 optimization algorithms, ≥7 sampling procedures, and 26 climate scenarios. Models based on temperature, day length, and leaf unfolding performed best, especially when calibrated with generalized simulated annealing and systematically balanced or stratified samples. Projected leaf coloration shifts between −13 and +20 days by 2080–2099.
Katharina Gallmeier, J. Xavier Prochaska, Peter Cornillon, Dimitris Menemenlis, and Madolyn Kelm
Geosci. Model Dev., 16, 7143–7170, https://doi.org/10.5194/gmd-16-7143-2023, https://doi.org/10.5194/gmd-16-7143-2023, 2023
Short summary
Short summary
This paper introduces an approach to evaluate numerical models of ocean circulation. We compare the structure of satellite-derived sea surface temperature anomaly (SSTa) instances determined by a machine learning algorithm at 10–80 km scales to those output by a high-resolution MITgcm run. The simulation over much of the ocean reproduces the observed distribution of SSTa patterns well. This general agreement, alongside a few notable exceptions, highlights the potential of this approach.
Angus Fotherby, Harold J. Bradbury, Jennifer L. Druhan, and Alexandra V. Turchyn
Geosci. Model Dev., 16, 7059–7074, https://doi.org/10.5194/gmd-16-7059-2023, https://doi.org/10.5194/gmd-16-7059-2023, 2023
Short summary
Short summary
We demonstrate how, given a simulation of fluid and rock interacting, we can emulate the system using machine learning. This means that, for a given initial condition, we can predict the final state, avoiding the simulation step once the model has been trained. We present a workflow for applying this approach to any fluid–rock simulation and showcase two applications to different fluid–rock simulations. This approach has applications for improving model development and sensitivity analyses.
Cited articles
Alnæs, M. S., Logg, A., Ølgaard, K. B., Rognes, M. E., and Wells, G. N.:
Unified Form Language: a domain-specific language for weak formulations
of partial differential equations, ACM T. Math. Software, 40, 9, https://doi.org/10.1145/2566630, 2014. a
Andreolli, C., Thierry, P., Borges, L., Skinner, G., and Yount, C.: Chapter 23
– Characterization and Optimization Methodology Applied to Stencil
Computations, in: High Performance Parallelism Pearls, edited by: Reinders, J.
and Jeffers, J., 377–396, Morgan Kaufmann, Boston, https://doi.org/10.1016/B978-0-12-802118-7.00023-6, 2015. a
Arbona, A., Miñano, B., Rigo, A., Bona, C., Palenzuela, C., Artigues, A.,
Bona-Casas, C., and Massó, J.: Simflowny 2: An upgraded platform for
scientific modeling and simulation, arXiv preprint arXiv:1702.04715, Computer Physics Communications, 229, 170–181, 2017. a
Asanović, K., Bodik, R., Catanzaro, B. C., Gebis, J. J., Husbands, P.,
Keutzer,
K., Patterson, D. A., Plishker, W. L., Shalf, J., Williams, S. W., and Yelick, K. A.:
The landscape of parallel computing research: A view from berkeley, Tech.
rep., Technical Report UCB/EECS-2006-183, EECS Department, University of
California, Berkeley, 2006. a
Backus, J.: The history of Fortran I, II, and III, in: History of programming
languages I, ACM SIGPLAN Notices, 13, 165–180, 1978. a
Barba, L. A. and Forsyth, G. F.: CFD Python: the 12 steps to Navier-Stokes
equations, Journal of Open Source Education, 9, 21, https://doi.org/10.21105/jose.00021, 2018. a, b, c
Bondhugula, U., Hartono, A., Ramanujam, J., and Sadayappan, P.: A Practical
Automatic Polyhedral Parallelizer and Locality Optimizer, in: Proceedings of
the 2008 ACM SIGPLAN Conference on Programming Language Design and
Implementation, PLDI 2008, 101–113, ACM, New York, NY, USA, https://doi.org/10.1145/1375581.1375595, 2008. a, b
Cárdenas, A. F. and Karplus, W. J.: PDEL – a language for partial
differential equations, Commun. ACM, 13, 184–191, 1970. a
Cauchy, A.-L.: Méthode générale pour la résolution des
systèmes d'équations simultanées, Compte Rendu des Séances
de L'Académie des Sciences XXV, Série A, 25, 536–538, 1847. a
Christen, M., Schenk, O., and Burkhart, H.: PATUS: A Code Generation and
Autotuning Framework for Parallel Iterative Stencil Computations on Modern
Microarchitectures, in: Proceedings of the 2011 IEEE International Parallel
& Distributed Processing Symposium, IPDPS 2011, 676–687, IEEE Computer
Society, Washington, DC, USA, https://doi.org/10.1109/IPDPS.2011.70, 2011. a
Clayton, R. and Engquist, B.: Absorbing boundary conditions for acoustic and
elastic wave equations, B. Seismol. Soc. Am., 67,
1529–1540, 1977. a
Colella, P.: Defining Software Requirements for Scientific Computing, DARPA
HPCS, 2004. a
Cook Jr., G. O.: ALPAL: A tool for the development of large-scale simulation
codes, Tech. rep., Lawrence Livermore National Lab., CA, USA, 1988. a
Dongarra, J.: The LINPACK Benchmark: An Explanation, in: Proceedings of the
1st International Conference on Supercomputing, 456–474,
Springer Verlag, London, UK, available at: http://dl.acm.org/citation.cfm?id=647970.742568
(last access: 24 March 2019), 1988. a
Farrell, P. E., Ham, D. A., Funke, S. W., and Rognes, M. E.: Automated
Derivation of the Adjoint of High-Level Transient Finite Element Programs,
SIAM J. Sci. Comput., 35, C369–C393,
https://doi.org/10.1137/120873558, 2013. a
Fomel, S., Sava, P., Vlad, I., Liu, Y., and Bashkardin, V.: Madagascar:
open-source software project for
multidimensional data analysis and reproducible computational experiments,
Journal of Open Research Software, 1, p.e8, https://doi.org/10.5334/jors.ag, 2013. a
Griewank, A. and Walther, A.: Algorithm 799: Revolve: An Implementation of
Checkpointing for the Reverse or Adjoint Mode of Computational
Differentiation, ACM Trans. Math. Softw., 26, 19–45,
https://doi.org/10.1145/347837.347846, 2000. a
Haber, E., Chung, M., and Herrmann, F. J.: An effective method for parameter
estimation with PDE constraints with multiple right hand sides, SIAM J. Optimiz., 22, 739–757, https://doi.org/10.1137/11081126X,
2012. a, b
Hawick, K. A. and Playne, D. P.: Simulation Software Generation using a
Domain-Specific Language for Partial Differential Field Equations, in: 11th
International Conference on Software Engineering Research and Practice
(SERP '13), CSTN-187, p. SER3829, WorldComp, Las Vegas, USA, 2013. a
Henretty, T., Veras, R., Franchetti, F., Pouchet, L.-N., Ramanujam, J., and
Sadayappan, P.: A Stencil Compiler for Short-vector SIMD Architectures, in:
Proceedings of the 27th International ACM Conference on International
Conference on Supercomputing, ICS '13, 13–24, ACM, New York, NY, USA,
https://doi.org/10.1145/2464996.2467268, 2013. a, b, c
Hopper, G. M.: The education of a computer, in: Proceedings of the 1952 ACM
national meeting (Pittsburgh), 243–249, ACM, 1952. a
Igel, H.: Computational Seismology: A Practical Introduction, Oxford
University Press, 1. edn., available at: https://global.oup.com/academic/product/computational-seismology-9780198717409?cc=de&lang=en&
(last access: 24 March 2019), 2016. a
Intel Corporation: Intel VTune Performance Analyzer,
available at: https://software.intel.com/en-us/intel-vtune-amplifier-xe (last access: 24 March 2019), 2016. a
Iverson, K. E.: A Programming Language, John Wiley & Sons, Inc., New York, NY,
USA, 1962. a
Jacobs, C. T., Jammy, S. P., and Sandham, N. D.: OpenSBLI: A framework for
the automated derivation and parallel execution of finite difference solvers
on a range of computer architectures, CoRR, abs/1609.01277, available at: http://arxiv.org/abs/1609.01277
(last access: 24 March 2019), 2016. a
Jones, J. L.: A survey of automatic coding techniques for digital computers,
Ph.D. thesis, Massachusetts Institute of Technology, Boston, MA, USA, 1954. a
Köster, M., Leißa, R., Hack, S., Membarth, R., and Slusallek, P.:
Platform-Specific Optimization and Mapping of Stencil Codes through
Refinement, in: Proceedings of the 1st International Workshop on
High-Performance Stencil Computations, 21 January 2014, Vienna, Austria, 1–6, 2014. a
Kukreja, N., Hückelheim, J., Lange, M., Louboutin, M., Walther, A., Funke,
S. W., and Gorman, G.: High-level python abstractions for optimal
checkpointing in inversion problems, arXiv preprint arXiv:1802.02474, 2018. a
Kumar, R., Wason, H., and Herrmann, F. J.: Source separation for simultaneous
towed-streamer marine acquisition – a compressed sensing
approach, Geophysics, 80, WD73–WD88, https://doi.org/10.1190/geo2015-0108.1, 2015. a
Lange, M., Kukreja, N., Luporini, F., Louboutin, M., Yount, C.,
Hückelheim, J., and Gorman, G. J.: Optimised finite difference
computation from symbolic equations, in: Proceedings of the 16th Python
in Science Conference (SciPy 2017), 10–16 July, Austin, Texas, edited by: Huff, K., Lippa, D., Niederhut, D., and Pacer, M., 89–96,
2017. a
Lions, J. L.: Optimal control of systems governed by partial differential
equations, 1st edn., Springer Verlag, Berlin, Heidelberg, 1971. a
Liu, Y. and Fomel, S.: Seismic data interpolation beyond aliasing using
regularized nonstationary autoregression, Geophysics, 76, V69–V77,
https://doi.org/10.1190/geo2010-0231.1, 2011. a
Logg, A., Mardal, K.-A., Wells, and Wells, G.: Automated Solution of
Differential Equations by the Finite Element Method, The FEniCS Book, Springer, https://doi.org/10.1007/978-3-642-23099-8, 2012. a, b, c
Louboutin, M., Lange, M., Herrmann, F. J., Kukreja, N., and Gorman, G.:
Performance prediction of finite-difference solvers for different computer
architectures, Comput. Geosci., 105, 148–157,
https://doi.org/10.1016/j.cageo.2017.04.014, 2017a. a, b
Louboutin, M., Witte, P., Lange, M., Kukreja, N., Luporini, F., Gorman, G., and
Herrmann, F. J.: Full-waveform inversion, Part 1: Forward modeling, The
Leading Edge, 36, 1033–1036, https://doi.org/10.1190/tle36121033.1, 2017b. a
Louboutin, M., Luporini, F., Lange, M., Kukreja, N., Pandolfo, V., Kazakas,
P., Zhang, S., Gorman, G., Hueckelheim, J., Peng, P., Velesko, P., and
McCormick, D.: opesci/devito: Devito-3.1, https://doi.org/10.5281/zenodo.1038305, 2017. a
Luporini, F., Varbanescu, A. L., Rathgeber, F., Bercea, G.-T., Ramanujam, J.,
Ham, D. A., and Kelly, P. H. J.: Cross-Loop Optimization of Arithmetic
Intensity for Finite Element Local Assembly, ACM Trans. Archit. Code Optim.,
11, 57:1–57:25, https://doi.org/10.1145/2687415, 2015. a
Luporini, F., Lange, M., Louboutin, M., Kukreja, N., Hückelheim,
J., Yount, C., Witte, P., Kelly, P. H. J., Gorman, G. J., and
Herrmann, F. J.: Architecture and performance of Devito, a system for
automated stencil computation, CoRR, abs/1807.03032, available at: http://arxiv.org/abs/1807.03032
(last access: 24 March 2019), 2018. a, b, c, d
McCalpin, J. D.: STREAM: Sustainable Memory Bandwidth in High Performance
Computers, Tech. rep., University of Virginia, Charlottesville, Virginia, available at: http://www.cs.virginia.edu/stream/
(last access: 24 March 2019), a continually updated
technical report, 1991–2007. a
McMechan, G. A.: Migration by extrapolation of time-dependent boundary values,
Geophys. Prospect., 31, 413–420, https://doi.org/10.1111/j.1365-2478.1983.tb01060.x, 1983. a
Membarth, R., Hannig, F., Teich, J., and Köstler, H.: Towards
domain-specific computing for stencil codes in HPC, in: High Performance
Computing, Networking, Storage and Analysis (SCC), Salt Lake City, UT, 2012 SC Companion, IEEE, 1133–1138, https://doi.org/10.1109/SC.Companion.2012.136, 2012. a
Meurer, A., Smith, C. P., Paprocki, M., Čertík, O., Kirpichev, S. B.,
Rocklin, M., Kumar, A., Ivanov, S., Moore, J. K., Singh, S., Rathnayake, T.,
Vig, S., Granger, B. E., Muller, R. P., Bonazzi, F., Gupta, H., Vats, S.,
Johansson, F., Pedregosa, F., Curry, M. J., Terrel, A. R., Roučka, V.,
Saboo, A., Fernando, I., Kulal, S., Cimrman, R., and Scopatz, A.: SymPy:
symbolic computing in Python, PeerJ Computer Science, 3, e103,
https://doi.org/10.7717/peerj-cs.103, 2017. a
Mittet, R.: Implementation of the Kirchhoff integral for elastic waves in
staggered-grid modeling schemes, Geophysics, 59, 1894–1901,
https://doi.org/10.1190/1.1443576, 1994. a
Naghizadeh, M. and Sacchi, M. D.: f-x adaptive seismic-trace interpolation,
Geophysics, 74, V9–V16, https://doi.org/10.1190/1.3008547, 2009. a
Osuna, C., Fuhrer, O., Gysi, T., and Bianco, M.: STELLA: A domain-specific
language for stencil methods on structured grids, in: Poster Presentation at
the Platform for Advanced Scientific Computing (PASC) Conference, Zurich,
Switzerland, 2014. a
Patterson, D. A. and Hennessy, J. L.: Computer Organization and Design: The
Hardware/Software Interface, 3rd edn., Morgan Kaufmann Publishers Inc., San Francisco,
CA, USA, 2007. a
Peters, B. and Herrmann, F. J.: Constraints versus penalties for
edge-preserving full-waveform inversion, The Leading Edge, 36, 94–100,
https://doi.org/10.1190/tle36010094.1, 2017. a, b
Plessix, R.-E.: A review of the adjoint-state method for computing the gradient
of a functional with geophysical applications, Geophys. J. Int., 167, 495–503, https://doi.org/10.1111/j.1365-246X.2006.02978.x, 2006. a
Raknes, E. B. and Weibull, W.: Efficient 3D elastic full-waveform inversion
using wavefield reconstruction methods, Geophysics, 81, R45–R55,
https://doi.org/10.1190/geo2015-0185.1, 2016. a
Rathgeber, F., Ham, D. A., Mitchell, L., Lange, M., Luporini, F., McRae, A.
T. T., Bercea, G., Markall, G. R., and Kelly, P. H. J.: Firedrake: automating
the finite element method by composing abstractions, CoRR, abs/1501.01809, available at: http://arxiv.org/abs/1501.01809
(last access: 24 March 2019), 2015. a, b, c
Schmidt, M., van den Berg, E., Friedlander, M. P., and Murphy, K.: Optimizing
Costly Functions with Simple Constraints: A Limited-Memory Projected
Quasi-Newton Algorithm, in: Proceedings of The Twelfth International
Conference on Artificial Intelligence and Statistics (AISTATS) 2009, edited
by: van Dyk, D. and Welling, M., vol. 5, 456–463, Clearwater Beach,
Florida, 2009. a
Seongjai Kim, H. L.: High-order schemes for acoustic waveform simulation,
Appl. Numer. Math., 57, 402–414, 2007. a
Sun, D. and Symes, W. W.: IWAVE implementation of adjoint state method, Tech.
Rep. 10-06, Department of Computational and Applied Mathematics,
Rice University, Houston, Texas, USA, available at: https://pdfs.semanticscholar.org/6c17/cfe41b76f6b745c435891ea\%6ba6f4e2c2dbf.pdf
(last access: 24 March 2019), 2010. a
Symes, W. W.: IWAVE structure and basic use cases, The Rice Inversion Project,
p. 85, available at: http://www.trip.caam.rice.edu/reports/2014/book.pdf#page=89 (last access: 24 March 2019),
2015b. a
Symes, W. W., Sun, D., and Enriquez, M.: From modelling to inversion: designing
a well-adapted simulator, Geophys. Prospect., 59, 814–833,
https://doi.org/10.1111/j.1365-2478.2011.00977.x, 2011. a
Tarantola, A.: Inversion of seismic reflection data in the acoustic
approximation, Geophysics, 49, 1259, 1140–1395, https://doi.org/10.1190/1.1441754, 1984. a
Umetani, Y.: DEQSOL A numerical Simulation Language for Vector/Parallel
Processors, Proc. IFIP TC2/WG22, in: Parallel Processing for Scientific Computing, 5, 147–164, 1985. a
Unat, D., Cai, X., and Baden, S. B.: Mint: realizing CUDA performance in 3-D
stencil methods with annotated C, in: ICS '11 International Conference on Supercomputing, Tucson, AZ, USA,
214–224, ACM, New York, NY, USA, 2011. a
Van Engelen, R., Wolters, L., and Cats, G.: Ctadel: A generator of
multi-platform high performance codes for pde-based scientific applications,
in: ICS '96 Proceedings of the 10th international conference on Supercomputing, 86–93, ACM, Philadelphia, Pennsylvania, USA, 1996. a
van Leeuwen, T. and Herrmann, F. J.: A penalty method for PDE-constrained
optimization in inverse problems, Inverse Probl., 32, 015007,
https://doi.org/10.1088/0266-5611/32/1/015007, 2015. a
Versteeg, R.: The Marmousi experience; velocity model determination on a
synthetic complex data set, The Leading Edge, 13, 927–936, 1994. a
Virieux, J.: P-SV wave propagation in heterogeneous media: Velocity-stress
finite-difference method, Geophysics, 51, 889–901, https://doi.org/10.1190/1.1442147, 1986. a
Virieux, J. and Operto, S.: An overview of full-waveform inversion in
exploration geophysics, Geophysics, 74, WCC1–WCC26, https://doi.org/10.1190/1.3238367, 2009. a
Wang, S., Nissen, A., and Kreiss, G.: Convergence of finite difference methods
for the wave equation in two space dimensions, Computing Research Repository, available at:
https://arxiv.org/abs/1702.01383 (last access: 24 March 2019), 2017. a
Warner, M. and Guasch, L.: Adaptive waveform inversion: Theory, 1089–1093,
https://doi.org/10.1190/segam2014-0371.1, SEG 2014, Denver, Colorado, 2014. a
Wason, H., Oghenekohwo, F., and Herrmann, F. J.: Low-cost time-lapse seismic
with distributed compressive sensing – Part 2: impact on
repeatability, Geophysics, 82, P15–P30, https://doi.org/10.1190/geo2016-0252.1, 2017. a
Watanabe, K.: Green's Functions for Laplace and Wave Equations, pp. 33–76,
Springer International Publishing, Cham, https://doi.org/10.1007/978-3-319-17455-6_2, 2015. a
Weiss, R. M. and Shragge, J.: Solving 3D anisotropic elastic wave equations
on parallel GPU devices., Geophysics, 78, 1Ma-Z51, https://doi.org/10.1190/geo2012-0063.1, 2013. a
Williams, S., Waterman, A., and Patterson, D.: The Roofline model offers
insight on how to improve the performance of software and hardware,
Communications of the ACM – A Direct Path to Dependable Software CACM Homepage archive, 52, 65–76, 2009. a
Witte, P., Louboutin, M., Lensink, K., Lange, M., Kukreja, N., Luporini, F.,
Gorman, G., and Herrmann, F. J.: Full-waveform inversion, Part 3:
Optimization, The Leading Edge, 37, 142–145, https://doi.org/10.1190/tle37020142.1, 2018.
a, b, c
Witte, P. A., Louboutin, M., Kukreja, N., Luporini, F., Lange, M., Gorman,
G. J., and Herrmann, F. J.: A large-scale framework for symbolic
implementations of seismic inversion algorithms in Julia,
available at: https://www.slim.eos.ubc.ca/Publications/Private/Submitted/20\%18/witte2018alf/witte2018alf.html
(last access: 24 March 2019),
Geophysics, 84, F57–F71, https://doi.org/10.1190/geo2018-0174.1, 2019. a, b
Yount, C.: Vector Folding: Improving Stencil Performance via Multi-dimensional
SIMD-vector Representation, in: 2015 IEEE 17th International Conference on
High Performance Computing and Communications, 2015 IEEE 7th International
Symposium on Cyberspace Safety and Security, and 2015 IEEE 12th International
Conference on Embedded Software and Systems, 24–26 August 2015, 865–870,
https://doi.org/10.1109/HPCC-CSS-ICESS.2015.27, 2015. a, b, c, d
Zhang, Y. and Mueller, F.: Auto-generation and Auto-tuning of 3D Stencil Codes
on GPU Clusters, in: Proceedings of the Tenth International Symposium on Code
Generation and Optimization, CGO '12, 155–164, ACM, New York, NY, USA,
https://doi.org/10.1145/2259016.2259037, 2012. a
Short summary
This paper presents Devito, a Python-based software. The aim of this software is to provide a high-level simple interface to users for the description and discretization of the mathematical definition of the physics. This research initially started as an attempt to improve research time, portability, and performance in exploration geophysics. We present the latest version of the software that is already making an impact in academics and industry.
This paper presents Devito, a Python-based software. The aim of this software is to provide a...