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.
Model description paper
| Highlight paper
| 
27 Mar 2019
Model description paper | Highlight paper |
| 27 Mar 2019
| 27 Mar 2019
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
Impact of ocean vertical-mixing parameterization on Arctic sea ice and upper-ocean properties using the NEMO-SI3 model
Bridging the gap: a new module for human water use in the Community Earth System Model version 2.2.1
A new lightning scheme in the Canadian Atmospheric Model (CanAM5.1): implementation, evaluation, and projections of lightning and fire in future climates
Methane dynamics in the Baltic Sea: investigating concentration, flux, and isotopic composition patterns using the coupled physical–biogeochemical model BALTSEM-CH4 v1.0
Split-explicit external mode solver in the finite volume sea ice–ocean model FESOM2
Applying double cropping and interactive irrigation in the North China Plain using WRF4.5
The sea ice component of GC5: coupling SI3 to HadGEM3 using conductive fluxes
CICE on a C-grid: new momentum, stress, and transport schemes for CICEv6.5
HyPhAICC v1.0: a hybrid physics–AI approach for probability fields advection shown through an application to cloud cover nowcasting
CICERO Simple Climate Model (CICERO-SCM v1.1.1) – an improved simple climate model with a parameter calibration tool
Development of a plant carbon–nitrogen interface coupling framework in a coupled biophysical-ecosystem–biogeochemical model (SSiB5/TRIFFID/DayCent-SOM v1.0)
Dynamical Madden–Julian Oscillation forecasts using an ensemble subseasonal-to-seasonal forecast system of the IAP-CAS model
Implementation of a brittle sea ice rheology in an Eulerian, finite-difference, C-grid modeling framework: impact on the simulated deformation of sea ice in the Arctic
HSW-V v1.0: localized injections of interactive volcanic aerosols and their climate impacts in a simple general circulation model
A 3D-Var assimilation scheme for vertical velocity with CMA-MESO v5.0
Updating the radiation infrastructure in MESSy (based on MESSy version 2.55)
An urban module coupled with the Variable Infiltration Capacity model to improve hydrothermal simulations in urban systems
Bayesian hierarchical model for bias-correcting climate models
Evaluation of the coupling of EMACv2.55 to the land surface and vegetation model JSBACHv4
Reduced floating-point precision in regional climate simulations: an ensemble-based statistical verification
TorchClim v1.0: a deep-learning plugin for climate model physics
Linking global terrestrial and ocean biogeochemistry with process-based, coupled freshwater algae–nutrient–solid dynamics in LM3-FANSY v1.0
Validating a microphysical prognostic stratospheric aerosol implementation in E3SMv2 using observations after the Mount Pinatubo eruption
Architectural Insights and Training Methodology Optimization of Pangu-Weather
Implementing detailed nucleation predictions in the Earth system model EC-Earth3.3.4: sulfuric acid–ammonia nucleation
Modeling biochar effects on soil organic carbon on croplands in a microbial decomposition model (MIMICS-BC_v1.0)
Hector V3.2.0: functionality and performance of a reduced-complexity climate model
Evaluation of CMIP6 model simulations of PM2.5 and its components over China
Robust handling of extremes in quantile mapping – "Murder your darlings"
Assessment of a tiling energy budget approach in a land surface model, ORCHIDEE-MICT (r8205)
Multivariate adjustment of drizzle bias using machine learning in European climate projections
Development and evaluation of the interactive Model for Air Pollution and Land Ecosystems (iMAPLE) version 1.0
A perspective on the next generation of Earth system model scenarios: towards representative emission pathways (REPs)
Evaluating downscaled products with expected hydroclimatic co-variances
Software sustainability of global impact models
CARIB12: A Regional Community Earth System Model / Modular Ocean Model 6 Configuration of the Caribbean Sea
Parallel SnowModel (v1.0): a parallel implementation of a distributed snow-evolution modeling system (SnowModel)
LB-SCAM: a learning-based method for efficient large-scale sensitivity analysis and tuning of the Single Column Atmosphere Model (SCAM)
Quantifying the impact of SST feedback frequency on Madden–Julian oscillation simulations
Systematic and objective evaluation of Earth system models: PCMDI Metrics Package (PMP) version 3
A revised model of global silicate weathering considering the influence of vegetation cover on erosion rate
Evaluation of global fire simulations in CMIP6 Earth system models
A radiative–convective model computing precipitation with the maximum entropy production hypothesis
Design, evaluation and future projections of the NARCliM2.0 CORDEX-CMIP6 Australasia regional climate ensemble
Introducing the MESMER-M-TPv0.1.0 module: Spatially Explicit Earth System Model Emulation for Monthly Precipitation and Temperature
Leveraging regional mesh refinement to simulate future climate projections for California using the Simplified Convection-Permitting E3SM Atmosphere Model Version 0
Machine learning parameterization of the multi-scale Kain–Fritsch (MSKF) convection scheme and stable simulation coupled in the Weather Research and Forecasting (WRF) model using WRF–ML v1.0
A computationally light-weight model for ensemble forecasting of environmental hazard: General TAMSAT-ALERT v1.2.1
Impacts of spatial heterogeneity of anthropogenic aerosol emissions in a regionally refined global aerosol–climate model
cfr (v2024.1.26): a Python package for climate field reconstruction
Sofia Allende, Anne Marie Treguier, Camille Lique, Clément de Boyer Montégut, François Massonnet, Thierry Fichefet, and Antoine Barthélemy
Geosci. Model Dev., 17, 7445–7466, https://doi.org/10.5194/gmd-17-7445-2024, https://doi.org/10.5194/gmd-17-7445-2024, 2024
Short summary
Short summary
We study the parameters of the turbulent-kinetic-energy mixed-layer-penetration scheme in the NEMO model with regard to sea-ice-covered regions of the Arctic Ocean. This evaluation reveals the impact of these parameters on mixed-layer depth, sea surface temperature and salinity, and ocean stratification. Our findings demonstrate significant impacts on sea ice thickness and sea ice concentration, emphasizing the need for accurately representing ocean mixing to understand Arctic climate dynamics.
Sabin I. Taranu, David M. Lawrence, Yoshihide Wada, Ting Tang, Erik Kluzek, Sam Rabin, Yi Yao, Steven J. De Hertog, Inne Vanderkelen, and Wim Thiery
Geosci. Model Dev., 17, 7365–7399, https://doi.org/10.5194/gmd-17-7365-2024, https://doi.org/10.5194/gmd-17-7365-2024, 2024
Short summary
Short summary
In this study, we improved a climate model by adding the representation of water use sectors such as domestic, industry, and agriculture. This new feature helps us understand how water is used and supplied in various areas. We tested our model from 1971 to 2010 and found that it accurately identifies areas with water scarcity. By modelling the competition between sectors when water availability is limited, the model helps estimate the intensity and extent of individual sectors' water shortages.
Cynthia Whaley, Montana Etten-Bohm, Courtney Schumacher, Ayodeji Akingunola, Vivek Arora, Jason Cole, Michael Lazare, David Plummer, Knut von Salzen, and Barbara Winter
Geosci. Model Dev., 17, 7141–7155, https://doi.org/10.5194/gmd-17-7141-2024, https://doi.org/10.5194/gmd-17-7141-2024, 2024
Short summary
Short summary
This paper describes how lightning was added as a process in the Canadian Earth System Model in order to interactively respond to climate changes. As lightning is an important cause of global wildfires, this new model development allows for more realistic projections of how wildfires may change in the future, responding to a changing climate.
Erik Gustafsson, Bo G. Gustafsson, Martijn Hermans, Christoph Humborg, and Christian Stranne
Geosci. Model Dev., 17, 7157–7179, https://doi.org/10.5194/gmd-17-7157-2024, https://doi.org/10.5194/gmd-17-7157-2024, 2024
Short summary
Short summary
Methane (CH4) cycling in the Baltic Proper is studied through model simulations, enabling a first estimate of key CH4 fluxes. A preliminary budget identifies benthic CH4 release as the dominant source and two main sinks: CH4 oxidation in the water (92 % of sinks) and outgassing to the atmosphere (8 % of sinks). This study addresses CH4 emissions from coastal seas and is a first step toward understanding the relative importance of open-water outgassing compared with local coastal hotspots.
Tridib Banerjee, Patrick Scholz, Sergey Danilov, Knut Klingbeil, and Dmitry Sidorenko
Geosci. Model Dev., 17, 7051–7065, https://doi.org/10.5194/gmd-17-7051-2024, https://doi.org/10.5194/gmd-17-7051-2024, 2024
Short summary
Short summary
In this paper we propose a new alternative to one of the functionalities of the sea ice model FESOM2. The alternative we propose allows the model to capture and simulate fast changes in quantities like sea surface elevation more accurately. We also demonstrate that the new alternative is faster and more adept at taking advantages of highly parallelized computing infrastructure. We therefore show that this new alternative is a great addition to the sea ice model FESOM2.
Yuwen Fan, Zhao Yang, Min-Hui Lo, Jina Hur, and Eun-Soon Im
Geosci. Model Dev., 17, 6929–6947, https://doi.org/10.5194/gmd-17-6929-2024, https://doi.org/10.5194/gmd-17-6929-2024, 2024
Short summary
Short summary
Irrigated agriculture in the North China Plain (NCP) has a significant impact on the local climate. To better understand this impact, we developed a specialized model specifically for the NCP region. This model allows us to simulate the double-cropping vegetation and the dynamic irrigation practices that are commonly employed in the NCP. This model shows improved performance in capturing the general crop growth, such as crop stages, biomass, crop yield, and vegetation greenness.
Ed Blockley, Emma Fiedler, Jeff Ridley, Luke Roberts, Alex West, Dan Copsey, Daniel Feltham, Tim Graham, David Livings, Clement Rousset, David Schroeder, and Martin Vancoppenolle
Geosci. Model Dev., 17, 6799–6817, https://doi.org/10.5194/gmd-17-6799-2024, https://doi.org/10.5194/gmd-17-6799-2024, 2024
Short summary
Short summary
This paper documents the sea ice model component of the latest Met Office coupled model configuration, which will be used as the physical basis for UK contributions to CMIP7. Documentation of science options used in the configuration are given along with a brief model evaluation. This is the first UK configuration to use NEMO’s new SI3 sea ice model. We provide details on how SI3 was adapted to work with Met Office coupling methodology and documentation of coupling processes in the model.
Jean-François Lemieux, William H. Lipscomb, Anthony Craig, David A. Bailey, Elizabeth C. Hunke, Philippe Blain, Till A. S. Rasmussen, Mats Bentsen, Frédéric Dupont, David Hebert, and Richard Allard
Geosci. Model Dev., 17, 6703–6724, https://doi.org/10.5194/gmd-17-6703-2024, https://doi.org/10.5194/gmd-17-6703-2024, 2024
Short summary
Short summary
We present the latest version of the CICE model. It solves equations that describe the dynamics and the growth and melt of sea ice. To do so, the domain is divided into grid cells and variables are positioned at specific locations in the cells. A new implementation (C-grid) is presented, with the velocity located on cell edges. Compared to the previous B-grid, the C-grid allows for a natural coupling with some oceanic and atmospheric models. It also allows for ice transport in narrow channels.
Rachid El Montassir, Olivier Pannekoucke, and Corentin Lapeyre
Geosci. Model Dev., 17, 6657–6681, https://doi.org/10.5194/gmd-17-6657-2024, https://doi.org/10.5194/gmd-17-6657-2024, 2024
Short summary
Short summary
This study introduces a novel approach that combines physics and artificial intelligence (AI) for improved cloud cover forecasting. This approach outperforms traditional deep learning (DL) methods in producing realistic and physically consistent results while requiring less training data. This architecture provides a promising solution to overcome the limitations of classical AI methods and contributes to open up new possibilities for combining physical knowledge with deep learning models.
Marit Sandstad, Borgar Aamaas, Ane Nordlie Johansen, Marianne Tronstad Lund, Glen Philip Peters, Bjørn Hallvard Samset, Benjamin Mark Sanderson, and Ragnhild Bieltvedt Skeie
Geosci. Model Dev., 17, 6589–6625, https://doi.org/10.5194/gmd-17-6589-2024, https://doi.org/10.5194/gmd-17-6589-2024, 2024
Short summary
Short summary
The CICERO-SCM has existed as a Fortran model since 1999 that calculates the radiative forcing and concentrations from emissions and is an upwelling diffusion energy balance model of the ocean that calculates temperature change. In this paper, we describe an updated version ported to Python and publicly available at https://github.com/ciceroOslo/ciceroscm (https://doi.org/10.5281/zenodo.10548720). This version contains functionality for parallel runs and automatic calibration.
Zheng Xiang, Yongkang Xue, Weidong Guo, Melannie D. Hartman, Ye Liu, and William J. Parton
Geosci. Model Dev., 17, 6437–6464, https://doi.org/10.5194/gmd-17-6437-2024, https://doi.org/10.5194/gmd-17-6437-2024, 2024
Short summary
Short summary
A process-based plant carbon (C)–nitrogen (N) interface coupling framework has been developed which mainly focuses on plant resistance and N-limitation effects on photosynthesis, plant respiration, and plant phenology. A dynamic C / N ratio is introduced to represent plant resistance and self-adjustment. The framework has been implemented in a coupled biophysical-ecosystem–biogeochemical model, and testing results show a general improvement in simulating plant properties with this framework.
Yangke Liu, Qing Bao, Bian He, Xiaofei Wu, Jing Yang, Yimin Liu, Guoxiong Wu, Tao Zhu, Siyuan Zhou, Yao Tang, Ankang Qu, Yalan Fan, Anling Liu, Dandan Chen, Zhaoming Luo, Xing Hu, and Tongwen Wu
Geosci. Model Dev., 17, 6249–6275, https://doi.org/10.5194/gmd-17-6249-2024, https://doi.org/10.5194/gmd-17-6249-2024, 2024
Short summary
Short summary
We give an overview of the Institute of Atmospheric Physics–Chinese Academy of Sciences subseasonal-to-seasonal ensemble forecasting system and Madden–Julian Oscillation forecast evaluation of the system. Compared to other S2S models, the IAP-CAS model has its benefits but also biases, i.e., underdispersive ensemble, overestimated amplitude, and faster propagation speed when forecasting MJO. We provide a reason for these biases and prospects for further improvement of this system in the future.
Laurent Brodeau, Pierre Rampal, Einar Ólason, and Véronique Dansereau
Geosci. Model Dev., 17, 6051–6082, https://doi.org/10.5194/gmd-17-6051-2024, https://doi.org/10.5194/gmd-17-6051-2024, 2024
Short summary
Short summary
A new brittle sea ice rheology, BBM, has been implemented into the sea ice component of NEMO. We describe how a new spatial discretization framework was introduced to achieve this. A set of idealized and realistic ocean and sea ice simulations of the Arctic have been performed using BBM and the standard viscous–plastic rheology of NEMO. When compared to satellite data, our simulations show that our implementation of BBM leads to a fairly good representation of sea ice deformations.
Joseph P. Hollowed, Christiane Jablonowski, Hunter Y. Brown, Benjamin R. Hillman, Diana L. Bull, and Joseph L. Hart
Geosci. Model Dev., 17, 5913–5938, https://doi.org/10.5194/gmd-17-5913-2024, https://doi.org/10.5194/gmd-17-5913-2024, 2024
Short summary
Short summary
Large volcanic eruptions deposit material in the upper atmosphere, which is capable of altering temperature and wind patterns of Earth's atmosphere for subsequent years. This research describes a new method of simulating these effects in an idealized, efficient atmospheric model. A volcanic eruption of sulfur dioxide is described with a simplified set of physical rules, which eventually cools the planetary surface. This model has been designed as a test bed for climate attribution studies.
Hong Li, Yi Yang, Jian Sun, Yuan Jiang, Ruhui Gan, and Qian Xie
Geosci. Model Dev., 17, 5883–5896, https://doi.org/10.5194/gmd-17-5883-2024, https://doi.org/10.5194/gmd-17-5883-2024, 2024
Short summary
Short summary
Vertical atmospheric motions play a vital role in convective-scale precipitation forecasts by connecting atmospheric dynamics with cloud development. A three-dimensional variational vertical velocity assimilation scheme is developed within the high-resolution CMA-MESO model, utilizing the adiabatic Richardson equation as the observation operator. A 10 d continuous run and an individual case study demonstrate improved forecasts, confirming the scheme's effectiveness.
Matthias Nützel, Laura Stecher, Patrick Jöckel, Franziska Winterstein, Martin Dameris, Michael Ponater, Phoebe Graf, and Markus Kunze
Geosci. Model Dev., 17, 5821–5849, https://doi.org/10.5194/gmd-17-5821-2024, https://doi.org/10.5194/gmd-17-5821-2024, 2024
Short summary
Short summary
We extended the infrastructure of our modelling system to enable the use of an additional radiation scheme. After calibrating the model setups to the old and the new radiation scheme, we find that the simulation with the new scheme shows considerable improvements, e.g. concerning the cold-point temperature and stratospheric water vapour. Furthermore, perturbations of radiative fluxes associated with greenhouse gas changes, e.g. of methane, tend to be improved when the new scheme is employed.
Yibing Wang, Xianhong Xie, Bowen Zhu, Arken Tursun, Fuxiao Jiang, Yao Liu, Dawei Peng, and Buyun Zheng
Geosci. Model Dev., 17, 5803–5819, https://doi.org/10.5194/gmd-17-5803-2024, https://doi.org/10.5194/gmd-17-5803-2024, 2024
Short summary
Short summary
Urban expansion intensifies challenges like urban heat and urban dry islands. To address this, we developed an urban module, VIC-urban, in the Variable Infiltration Capacity (VIC) model. Tested in Beijing, VIC-urban accurately simulated turbulent heat fluxes, runoff, and land surface temperature. We provide a reliable tool for large-scale simulations considering urban environment and a systematic urban modelling framework within VIC, offering crucial insights for urban planners and designers.
Jeremy Carter, Erick A. Chacón-Montalván, and Amber Leeson
Geosci. Model Dev., 17, 5733–5757, https://doi.org/10.5194/gmd-17-5733-2024, https://doi.org/10.5194/gmd-17-5733-2024, 2024
Short summary
Short summary
Climate models are essential tools in the study of climate change and its wide-ranging impacts on life on Earth. However, the output is often afflicted with some bias. In this paper, a novel model is developed to predict and correct bias in the output of climate models. The model captures uncertainty in the correction and explicitly models underlying spatial correlation between points. These features are of key importance for climate change impact assessments and resulting decision-making.
Anna Martin, Veronika Gayler, Benedikt Steil, Klaus Klingmüller, Patrick Jöckel, Holger Tost, Jos Lelieveld, and Andrea Pozzer
Geosci. Model Dev., 17, 5705–5732, https://doi.org/10.5194/gmd-17-5705-2024, https://doi.org/10.5194/gmd-17-5705-2024, 2024
Short summary
Short summary
The study evaluates the land surface and vegetation model JSBACHv4 as a replacement for the simplified submodel SURFACE in EMAC. JSBACH mitigates earlier problems of soil dryness, which are critical for vegetation modelling. When analysed using different datasets, the coupled model shows strong correlations of key variables, such as land surface temperature, surface albedo and radiation flux. The versatility of the model increases significantly, while the overall performance does not degrade.
Hugo Banderier, Christian Zeman, David Leutwyler, Stefan Rüdisühli, and Christoph Schär
Geosci. Model Dev., 17, 5573–5586, https://doi.org/10.5194/gmd-17-5573-2024, https://doi.org/10.5194/gmd-17-5573-2024, 2024
Short summary
Short summary
We investigate the effects of reduced-precision arithmetic in a state-of-the-art regional climate model by studying the results of 10-year-long simulations. After this time, the results of the reduced precision and the standard implementation are hardly different. This should encourage the use of reduced precision in climate models to exploit the speedup and memory savings it brings. The methodology used in this work can help researchers verify reduced-precision implementations of their model.
David Fuchs, Steven C. Sherwood, Abhnil Prasad, Kirill Trapeznikov, and Jim Gimlett
Geosci. Model Dev., 17, 5459–5475, https://doi.org/10.5194/gmd-17-5459-2024, https://doi.org/10.5194/gmd-17-5459-2024, 2024
Short summary
Short summary
Machine learning (ML) of unresolved processes offers many new possibilities for improving weather and climate models, but integrating ML into the models has been an engineering challenge, and there are performance issues. We present a new software plugin for this integration, TorchClim, that is scalable and flexible and thereby allows a new level of experimentation with the ML approach. We also provide guidance on ML training and demonstrate a skillful hybrid ML atmosphere model.
Minjin Lee, Charles A. Stock, John P. Dunne, and Elena Shevliakova
Geosci. Model Dev., 17, 5191–5224, https://doi.org/10.5194/gmd-17-5191-2024, https://doi.org/10.5194/gmd-17-5191-2024, 2024
Short summary
Short summary
Modeling global freshwater solid and nutrient loads, in both magnitude and form, is imperative for understanding emerging eutrophication problems. Such efforts, however, have been challenged by the difficulty of balancing details of freshwater biogeochemical processes with limited knowledge, input, and validation datasets. Here we develop a global freshwater model that resolves intertwined algae, solid, and nutrient dynamics and provide performance assessment against measurement-based estimates.
Hunter York Brown, Benjamin Wagman, Diana Bull, Kara Peterson, Benjamin Hillman, Xiaohong Liu, Ziming Ke, and Lin Lin
Geosci. Model Dev., 17, 5087–5121, https://doi.org/10.5194/gmd-17-5087-2024, https://doi.org/10.5194/gmd-17-5087-2024, 2024
Short summary
Short summary
Explosive volcanic eruptions lead to long-lived, microscopic particles in the upper atmosphere which act to cool the Earth's surface by reflecting the Sun's light back to space. We include and test this process in a global climate model, E3SM. E3SM is tested against satellite and balloon observations of the 1991 eruption of Mt. Pinatubo, showing that with these particles in the model we reasonably recreate Pinatubo and its global effects. We also explore how particle size leads to these effects.
Deifilia Aurora To, Julian Quinting, Gholam Ali Hoshyaripour, Markus Götz, Achim Streit, and Charlotte Debus
EGUsphere, https://doi.org/10.5194/egusphere-2024-1714, https://doi.org/10.5194/egusphere-2024-1714, 2024
Short summary
Short summary
Pangu-Weather is a breakthrough machine learning model in medium-range weather forecasting that considers three-dimensional atmospheric information. We show that using a simpler 2D framework improves robustness, speeds up training, and reduces computational needs by 20–30%. We introduce a training procedure that varies the importance of atmospheric variables over time to speed up training convergence. Decreasing computational demand increases accessibility of training and working with the model.
Carl Svenhag, Moa K. Sporre, Tinja Olenius, Daniel Yazgi, Sara M. Blichner, Lars P. Nieradzik, and Pontus Roldin
Geosci. Model Dev., 17, 4923–4942, https://doi.org/10.5194/gmd-17-4923-2024, https://doi.org/10.5194/gmd-17-4923-2024, 2024
Short summary
Short summary
Our research shows the importance of modeling new particle formation (NPF) and growth of particles in the atmosphere on a global scale, as they influence the outcomes of clouds and our climate. With the global model EC-Earth3 we show that using a new method for NPF modeling, which includes new detailed processes with NH3 and H2SO4, significantly impacts the number of particles in the air and clouds and changes the radiation balance of the same magnitude as anthropogenic greenhouse emissions.
Mengjie Han, Qing Zhao, Xili Wang, Ying-Ping Wang, Philippe Ciais, Haicheng Zhang, Daniel S. Goll, Lei Zhu, Zhe Zhao, Zhixuan Guo, Chen Wang, Wei Zhuang, Fengchang Wu, and Wei Li
Geosci. Model Dev., 17, 4871–4890, https://doi.org/10.5194/gmd-17-4871-2024, https://doi.org/10.5194/gmd-17-4871-2024, 2024
Short summary
Short summary
The impact of biochar (BC) on soil organic carbon (SOC) dynamics is not represented in most land carbon models used for assessing land-based climate change mitigation. Our study develops a BC model that incorporates our current understanding of BC effects on SOC based on a soil carbon model (MIMICS). The BC model can reproduce the SOC changes after adding BC, providing a useful tool to couple dynamic land models to evaluate the effectiveness of BC application for CO2 removal from the atmosphere.
Kalyn Dorheim, Skylar Gering, Robert Gieseke, Corinne Hartin, Leeya Pressburger, Alexey N. Shiklomanov, Steven J. Smith, Claudia Tebaldi, Dawn L. Woodard, and Ben Bond-Lamberty
Geosci. Model Dev., 17, 4855–4869, https://doi.org/10.5194/gmd-17-4855-2024, https://doi.org/10.5194/gmd-17-4855-2024, 2024
Short summary
Short summary
Hector is an easy-to-use, global climate–carbon cycle model. With its quick run time, Hector can provide climate information from a run in a fraction of a second. Hector models on a global and annual basis. Here, we present an updated version of the model, Hector V3. In this paper, we document Hector’s new features. Hector V3 is capable of reproducing historical observations, and its future temperature projections are consistent with those of more complex models.
Fangxuan Ren, Jintai Lin, Chenghao Xu, Jamiu A. Adeniran, Jingxu Wang, Randall V. Martin, Aaron van Donkelaar, Melanie S. Hammer, Larry W. Horowitz, Steven T. Turnock, Naga Oshima, Jie Zhang, Susanne Bauer, Kostas Tsigaridis, Øyvind Seland, Pierre Nabat, David Neubauer, Gary Strand, Twan van Noije, Philippe Le Sager, and Toshihiko Takemura
Geosci. Model Dev., 17, 4821–4836, https://doi.org/10.5194/gmd-17-4821-2024, https://doi.org/10.5194/gmd-17-4821-2024, 2024
Short summary
Short summary
We evaluate the performance of 14 CMIP6 ESMs in simulating total PM2.5 and its 5 components over China during 2000–2014. PM2.5 and its components are underestimated in almost all models, except that black carbon (BC) and sulfate are overestimated in two models, respectively. The underestimation is the largest for organic carbon (OC) and the smallest for BC. Models reproduce the observed spatial pattern for OC, sulfate, nitrate and ammonium well, yet the agreement is poorer for BC.
Peter Berg, Thomas Bosshard, Denica Bozhinova, Lars Bärring, Joakim Löw, Carolina Nilsson, Gustav Strandberg, Johan Södling, Johan Thuresson, Renate Wilcke, and Wei Yang
Geosci. Model Dev. Discuss., https://doi.org/10.5194/gmd-2024-98, https://doi.org/10.5194/gmd-2024-98, 2024
Revised manuscript accepted for GMD
Short summary
Short summary
When bias adjusting climate model data using quantile mapping, one needs to prescribe what to do at the tails of the distribution, where a larger range of data is likely encountered outside the calibration period. The end result is highly dependent on the method used, and we show that one needs to exclude data in the calibration range to activate the extrapolation functionality also in that time period, else there will be discontinuities in the timeseries.
Yi Xi, Chunjing Qiu, Yuan Zhang, Dan Zhu, Shushi Peng, Gustaf Hugelius, Jinfeng Chang, Elodie Salmon, and Philippe Ciais
Geosci. Model Dev., 17, 4727–4754, https://doi.org/10.5194/gmd-17-4727-2024, https://doi.org/10.5194/gmd-17-4727-2024, 2024
Short summary
Short summary
The ORCHIDEE-MICT model can simulate the carbon cycle and hydrology at a sub-grid scale but energy budgets only at a grid scale. This paper assessed the implementation of a multi-tiling energy budget approach in ORCHIDEE-MICT and found warmer surface and soil temperatures, higher soil moisture, and more soil organic carbon across the Northern Hemisphere compared with the original version.
Georgia Lazoglou, Theo Economou, Christina Anagnostopoulou, George Zittis, Anna Tzyrkalli, Pantelis Georgiades, and Jos Lelieveld
Geosci. Model Dev., 17, 4689–4703, https://doi.org/10.5194/gmd-17-4689-2024, https://doi.org/10.5194/gmd-17-4689-2024, 2024
Short summary
Short summary
This study focuses on the important issue of the drizzle bias effect in regional climate models, described by an over-prediction of the number of rainy days while underestimating associated precipitation amounts. For this purpose, two distinct methodologies are applied and rigorously evaluated. These results are encouraging for using the multivariate machine learning method random forest to increase the accuracy of climate models concerning the projection of the number of wet days.
Xu Yue, Hao Zhou, Chenguang Tian, Yimian Ma, Yihan Hu, Cheng Gong, Hui Zheng, and Hong Liao
Geosci. Model Dev., 17, 4621–4642, https://doi.org/10.5194/gmd-17-4621-2024, https://doi.org/10.5194/gmd-17-4621-2024, 2024
Short summary
Short summary
We develop the interactive Model for Air Pollution and Land Ecosystems (iMAPLE). The model considers the full coupling between carbon and water cycles, dynamic fire emissions, wetland methane emissions, biogenic volatile organic compound emissions, and trait-based ozone vegetation damage. Evaluations show that iMAPLE is a useful tool for the study of the interactions among climate, chemistry, and ecosystems.
Malte Meinshausen, Carl-Friedrich Schleussner, Kathleen Beyer, Greg Bodeker, Olivier Boucher, Josep G. Canadell, John S. Daniel, Aïda Diongue-Niang, Fatima Driouech, Erich Fischer, Piers Forster, Michael Grose, Gerrit Hansen, Zeke Hausfather, Tatiana Ilyina, Jarmo S. Kikstra, Joyce Kimutai, Andrew D. King, June-Yi Lee, Chris Lennard, Tabea Lissner, Alexander Nauels, Glen P. Peters, Anna Pirani, Gian-Kasper Plattner, Hans Pörtner, Joeri Rogelj, Maisa Rojas, Joyashree Roy, Bjørn H. Samset, Benjamin M. Sanderson, Roland Séférian, Sonia Seneviratne, Christopher J. Smith, Sophie Szopa, Adelle Thomas, Diana Urge-Vorsatz, Guus J. M. Velders, Tokuta Yokohata, Tilo Ziehn, and Zebedee Nicholls
Geosci. Model Dev., 17, 4533–4559, https://doi.org/10.5194/gmd-17-4533-2024, https://doi.org/10.5194/gmd-17-4533-2024, 2024
Short summary
Short summary
The scientific community is considering new scenarios to succeed RCPs and SSPs for the next generation of Earth system model runs to project future climate change. To contribute to that effort, we reflect on relevant policy and scientific research questions and suggest categories for representative emission pathways. These categories are tailored to the Paris Agreement long-term temperature goal, high-risk outcomes in the absence of further climate policy and worlds “that could have been”.
Seung H. Baek, Paul A. Ullrich, Bo Dong, and Jiwoo Lee
EGUsphere, https://doi.org/10.5194/egusphere-2024-1456, https://doi.org/10.5194/egusphere-2024-1456, 2024
Short summary
Short summary
We evaluate downscaled products by examining locally relevant covariances during convective and frontal precipitation events. Common statistical downscaling techniques preserve expected covariances during convective precipitation. However, they dampen future intensification of frontal precipitation captured in global climate models and dynamical downscaling. This suggests statistical downscaling may not fully resolve non-stationary hydrologic processes as compared to dynamical downscaling.
Emmanuel Nyenah, Petra Döll, Daniel S. Katz, and Robert Reinecke
Geosci. Model Dev. Discuss., https://doi.org/10.5194/gmd-2024-97, https://doi.org/10.5194/gmd-2024-97, 2024
Revised manuscript accepted for GMD
Short summary
Short summary
Research software is crucial for scientific progress but is often developed by scientists with limited training, time, and funding, leading to software that is hard to understand, (re)use, modify, and maintain. Our study across 10 research sectors highlights strengths in version control, open-source licensing, and documentation while emphasizing the need for containerization and code quality. Recommendations include workshops, code quality metrics, funding, and adherence to FAIR standards.
Giovanni G. Seijo-Ellis, Donata Giglio, Gustavo M. Marques, and Frank O. Bryan
EGUsphere, https://doi.org/10.5194/egusphere-2024-1378, https://doi.org/10.5194/egusphere-2024-1378, 2024
Short summary
Short summary
A CESM/MOM6 regional configuration of the Caribbean Sea was developed as a response to the rising need of high-resolution models for climate impact studies. The configuration is validated for the period of 2000–2020 and improves significant errors in a low resolution model. Oceanic properties are well represented. Patterns of freshwater associated with the Amazon river are well captured and the mean flows across the multiple passages in the Caribbean Sea agree with observations.
Ross Mower, Ethan D. Gutmann, Glen E. Liston, Jessica Lundquist, and Soren Rasmussen
Geosci. Model Dev., 17, 4135–4154, https://doi.org/10.5194/gmd-17-4135-2024, https://doi.org/10.5194/gmd-17-4135-2024, 2024
Short summary
Short summary
Higher-resolution model simulations are better at capturing winter snowpack changes across space and time. However, increasing resolution also increases the computational requirements. This work provides an overview of changes made to a distributed snow-evolution modeling system (SnowModel) to allow it to leverage high-performance computing resources. Continental simulations that were previously estimated to take 120 d can now be performed in 5 h.
Jiaxu Guo, Juepeng Zheng, Yidan Xu, Haohuan Fu, Wei Xue, Lanning Wang, Lin Gan, Ping Gao, Wubing Wan, Xianwei Wu, Zhitao Zhang, Liang Hu, Gaochao Xu, and Xilong Che
Geosci. Model Dev., 17, 3975–3992, https://doi.org/10.5194/gmd-17-3975-2024, https://doi.org/10.5194/gmd-17-3975-2024, 2024
Short summary
Short summary
To enhance the efficiency of experiments using SCAM, we train a learning-based surrogate model to facilitate large-scale sensitivity analysis and tuning of combinations of multiple parameters. Employing a hybrid method, we investigate the joint sensitivity of multi-parameter combinations across typical cases, identifying the most sensitive three-parameter combination out of 11. Subsequently, we conduct a tuning process aimed at reducing output errors in these cases.
Yung-Yao Lan, Huang-Hsiung Hsu, and Wan-Ling Tseng
Geosci. Model Dev., 17, 3897–3918, https://doi.org/10.5194/gmd-17-3897-2024, https://doi.org/10.5194/gmd-17-3897-2024, 2024
Short summary
Short summary
This study uses the CAM5–SIT coupled model to investigate the effects of SST feedback frequency on the MJO simulations with intervals at 30 min, 1, 3, 6, 12, 18, 24, and 30 d. The simulations become increasingly unrealistic as the frequency of the SST feedback decreases. Our results suggest that more spontaneous air--sea interaction (e.g., ocean response within 3 d in this study) with high vertical resolution in the ocean model is key to the realistic simulation of the MJO.
Jiwoo Lee, Peter J. Gleckler, Min-Seop Ahn, Ana Ordonez, Paul A. Ullrich, Kenneth R. Sperber, Karl E. Taylor, Yann Y. Planton, Eric Guilyardi, Paul Durack, Celine Bonfils, Mark D. Zelinka, Li-Wei Chao, Bo Dong, Charles Doutriaux, Chengzhu Zhang, Tom Vo, Jason Boutte, Michael F. Wehner, Angeline G. Pendergrass, Daehyun Kim, Zeyu Xue, Andrew T. Wittenberg, and John Krasting
Geosci. Model Dev., 17, 3919–3948, https://doi.org/10.5194/gmd-17-3919-2024, https://doi.org/10.5194/gmd-17-3919-2024, 2024
Short summary
Short summary
We introduce an open-source software, the PCMDI Metrics Package (PMP), developed for a comprehensive comparison of Earth system models (ESMs) with real-world observations. Using diverse metrics evaluating climatology, variability, and extremes simulated in thousands of simulations from the Coupled Model Intercomparison Project (CMIP), PMP aids in benchmarking model improvements across generations. PMP also enables efficient tracking of performance evolutions during ESM developments.
Haoyue Zuo, Yonggang Liu, Gaojun Li, Zhifang Xu, Liang Zhao, Zhengtang Guo, and Yongyun Hu
Geosci. Model Dev., 17, 3949–3974, https://doi.org/10.5194/gmd-17-3949-2024, https://doi.org/10.5194/gmd-17-3949-2024, 2024
Short summary
Short summary
Compared to the silicate weathering fluxes measured at large river basins, the current models tend to systematically overestimate the fluxes over the tropical region, which leads to an overestimation of the global total weathering flux. The most possible cause of such bias is found to be the overestimation of tropical surface erosion, which indicates that the tropical vegetation likely slows down physical erosion significantly. We propose a way of taking this effect into account in models.
Fang Li, Xiang Song, Sandy P. Harrison, Jennifer R. Marlon, Zhongda Lin, L. Ruby Leung, Jörg Schwinger, Virginie Marécal, Shiyu Wang, Daniel S. Ward, Xiao Dong, Hanna Lee, Lars Nieradzik, Sam S. Rabin, and Roland Séférian
Geosci. Model Dev. Discuss., https://doi.org/10.5194/gmd-2024-85, https://doi.org/10.5194/gmd-2024-85, 2024
Revised manuscript accepted for GMD
Short summary
Short summary
This study provides the first comprehensive assessment of historical fire simulations from 19 CMIP6 ESMs. Most models reproduce global total, spatial pattern, seasonality, and regional historical changes well, but fail to simulate the recent decline in global burned area and underestimate the fire sensitivity to wet-dry conditions. They addressed three critical issues in CMIP5. We present targeted guidance for fire scheme development and methodologies to generate reliable fire projections.
Quentin Pikeroen, Didier Paillard, and Karine Watrin
Geosci. Model Dev., 17, 3801–3814, https://doi.org/10.5194/gmd-17-3801-2024, https://doi.org/10.5194/gmd-17-3801-2024, 2024
Short summary
Short summary
All accurate climate models use equations with poorly defined parameters, where knobs for the parameters are turned to fit the observations. This process is called tuning. In this article, we use another paradigm. We use a thermodynamic hypothesis, the maximum entropy production, to compute temperatures, energy fluxes, and precipitation, where tuning is impossible. For now, the 1D vertical model is used for a tropical atmosphere. The correct order of magnitude of precipitation is computed.
Giovanni Di Virgilio, Jason Evans, Fei Ji, Eugene Tam, Jatin Kala, Julia Andrys, Christopher Thomas, Dipayan Choudhury, Carlos Rocha, Stephen White, Yue Li, Moutassem El Rafei, Rishav Goyal, Matthew Riley, and Jyothi Lingala
Geosci. Model Dev. Discuss., https://doi.org/10.5194/gmd-2024-87, https://doi.org/10.5194/gmd-2024-87, 2024
Revised manuscript accepted for GMD
Short summary
Short summary
We introduce new climate models that simulate Australia’s future climate at regional scales, including at an unprecedented resolution of 4 km for 1950–2100. We describe the model design process used to create these new climate models. We show how the new models perform relative to previous-generation models, and compare their climate projections. This work is of national and international relevance as it can help guide climate model design and the use and interpretation of climate projections.
Sarah Schöngart, Lukas Gudmundsson, Mathias Hauser, Peter Pfleiderer, Quentin Lejeune, Shruti Nath, Sonia Isabelle Seneviratne, and Carl-Friedrich Schleußner
EGUsphere, https://doi.org/10.5194/egusphere-2024-278, https://doi.org/10.5194/egusphere-2024-278, 2024
Short summary
Short summary
Precipitation and temperature are two of the most impact-relevant climatic variables. Their joint distribution largely determines the division into climate regimes. Yet, projecting precipitation and temperature data under different emission scenarios relies on complex models that are computationally expensive. In this study, we propose a method that allows to generate monthly means of local precipitation and temperature at low computational costs.
Jishi Zhang, Peter Bogenschutz, Qi Tang, Philip Cameron-smith, and Chengzhu Zhang
Geosci. Model Dev., 17, 3687–3731, https://doi.org/10.5194/gmd-17-3687-2024, https://doi.org/10.5194/gmd-17-3687-2024, 2024
Short summary
Short summary
We developed a regionally refined climate model that allows resolved convection and performed a 20-year projection to the end of the century. The model has a resolution of 3.25 km in California, which allows us to predict climate with unprecedented accuracy, and a resolution of 100 km for the rest of the globe to achieve efficient, self-consistent simulations. The model produces superior results in reproducing climate patterns over California that typical modern climate models cannot resolve.
Xiaohui Zhong, Xing Yu, and Hao Li
Geosci. Model Dev., 17, 3667–3685, https://doi.org/10.5194/gmd-17-3667-2024, https://doi.org/10.5194/gmd-17-3667-2024, 2024
Short summary
Short summary
In order to forecast localized warm-sector rainfall in the south China region, numerical weather prediction models are being run with finer grid spacing. The conventional convection parameterization (CP) performs poorly in the gray zone, necessitating the development of a scale-aware scheme. We propose a machine learning (ML) model to replace the scale-aware CP scheme. Evaluation against the original CP scheme has shown that the ML-based CP scheme can provide accurate and reliable predictions.
Emily Black, John Ellis, and Ross Maidment
Geosci. Model Dev. Discuss., https://doi.org/10.5194/gmd-2024-75, https://doi.org/10.5194/gmd-2024-75, 2024
Revised manuscript accepted for GMD
Short summary
Short summary
We present General TAMSAT-ALERT: a computationally lightweight and versatile tool for generating ensemble forecasts from time series data. General TAMSAT-ALERT is capable of combining multiple streams of monitoring and forecasting data into probabilistic hazard assessments. As such, it complements existing systems and enhances their utility for actionable hazard assessment.
Taufiq Hassan, Kai Zhang, Jianfeng Li, Balwinder Singh, Shixuan Zhang, Hailong Wang, and Po-Lun Ma
Geosci. Model Dev., 17, 3507–3532, https://doi.org/10.5194/gmd-17-3507-2024, https://doi.org/10.5194/gmd-17-3507-2024, 2024
Short summary
Short summary
Anthropogenic aerosol emissions are an essential part of global aerosol models. Significant errors can exist from the loss of emission heterogeneity. We introduced an emission treatment that significantly improved aerosol emission heterogeneity in high-resolution model simulations, with improvements in simulated aerosol surface concentrations. The emission treatment will provide a more accurate representation of aerosol emissions and their effects on climate.
Feng Zhu, Julien Emile-Geay, Gregory J. Hakim, Dominique Guillot, Deborah Khider, Robert Tardif, and Walter A. Perkins
Geosci. Model Dev., 17, 3409–3431, https://doi.org/10.5194/gmd-17-3409-2024, https://doi.org/10.5194/gmd-17-3409-2024, 2024
Short summary
Short summary
Climate field reconstruction encompasses methods that estimate the evolution of climate in space and time based on natural archives. It is useful to investigate climate variations and validate climate models, but its implementation and use can be difficult for non-experts. This paper introduces a user-friendly Python package called cfr to make these methods more accessible, thanks to the computational and visualization tools that facilitate efficient and reproducible research on past climates.
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...
