Pub Date : 2024-02-29DOI: 10.5194/gmd-17-1813-2024
A. Savita, J. Kjellsson, Robin Pilch Kedzierski, M. Latif, Tabea Rahm, Sebastian Wahl, Wonsun Park
Abstract. We examine the impact of horizontal resolution and model time step on the climate of the OpenIFS version 43r3 atmospheric general circulation model. A series of simulations for the period 1979–2019 are conducted with various horizontal resolutions (i.e. ∼100, ∼50, and ∼25 km) while maintaining the same time step (i.e. 15 min) and using different time steps (i.e. 60, 30, and 15 min) at 100 km horizontal resolution. We find that the surface zonal wind bias is significantly reduced over certain regions such as the Southern Ocean and the Northern Hemisphere mid-latitudes and in tropical and subtropical regions at a high horizontal resolution (i.e. ∼25 km). Similar improvement is evident too when using a coarse-resolution model (∼100 km) with a smaller time step (i.e. 30 and 15 min). We also find improvements in Rossby wave amplitude and phase speed, as well as in weather regime patterns, when a smaller time step or higher horizontal resolution is used. The improvement in the wind bias when using the shorter time step is mostly due to an increase in shallow and mid-level convection that enhances vertical mixing in the lower troposphere. The enhanced mixing allows frictional effects to influence a deeper layer and reduces wind and wind speed throughout the troposphere. However, precipitation biases generally increase with higher horizontal resolutions or smaller time steps, whereas the surface air temperature bias exhibits a small improvement over North America and the eastern Eurasian continent. We argue that the bias improvement in the highest-horizontal-resolution (i.e. ∼25 km) configuration benefits from a combination of both the enhanced horizontal resolution and the shorter time step. In summary, we demonstrate that, by reducing the time step in the coarse-resolution (∼100 km) OpenIFS model, one can alleviate some climate biases at a lower cost than by increasing the horizontal resolution.
{"title":"Assessment of climate biases in OpenIFS version 43r3 across model horizontal resolutions and time steps","authors":"A. Savita, J. Kjellsson, Robin Pilch Kedzierski, M. Latif, Tabea Rahm, Sebastian Wahl, Wonsun Park","doi":"10.5194/gmd-17-1813-2024","DOIUrl":"https://doi.org/10.5194/gmd-17-1813-2024","url":null,"abstract":"Abstract. We examine the impact of horizontal resolution and model time step on the climate of the OpenIFS version 43r3 atmospheric general circulation model. A series of simulations for the period 1979–2019 are conducted with various horizontal resolutions (i.e. ∼100, ∼50, and ∼25 km) while maintaining the same time step (i.e. 15 min) and using different time steps (i.e. 60, 30, and 15 min) at 100 km horizontal resolution. We find that the surface zonal wind bias is significantly reduced over certain regions such as the Southern Ocean and the Northern Hemisphere mid-latitudes and in tropical and subtropical regions at a high horizontal resolution (i.e. ∼25 km). Similar improvement is evident too when using a coarse-resolution model (∼100 km) with a smaller time step (i.e. 30 and 15 min). We also find improvements in Rossby wave amplitude and phase speed, as well as in weather regime patterns, when a smaller time step or higher horizontal resolution is used. The improvement in the wind bias when using the shorter time step is mostly due to an increase in shallow and mid-level convection that enhances vertical mixing in the lower troposphere. The enhanced mixing allows frictional effects to influence a deeper layer and reduces wind and wind speed throughout the troposphere. However, precipitation biases generally increase with higher horizontal resolutions or smaller time steps, whereas the surface air temperature bias exhibits a small improvement over North America and the eastern Eurasian continent. We argue that the bias improvement in the highest-horizontal-resolution (i.e. ∼25 km) configuration benefits from a combination of both the enhanced horizontal resolution and the shorter time step. In summary, we demonstrate that, by reducing the time step in the coarse-resolution (∼100 km) OpenIFS model, one can alleviate some climate biases at a lower cost than by increasing the horizontal resolution.\u0000","PeriodicalId":12799,"journal":{"name":"Geoscientific Model Development","volume":null,"pages":null},"PeriodicalIF":5.1,"publicationDate":"2024-02-29","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140413816","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2024-02-29DOI: 10.5194/gmd-17-1789-2024
Kees Nederhoff, M. van Ormondt, Jay Veeramony, A. V. van Dongeren, J. Antolínez, T. Leijnse, D. Roelvink
Abstract. Tropical-cyclone impacts can have devastating effects on the population, infrastructure, and natural habitats. However, predicting these impacts is difficult due to the inherent uncertainties in the storm track and intensity. In addition, due to computational constraints, both the relevant ocean physics and the uncertainties in meteorological forcing are only partly accounted for. This paper presents a new method, called the Tropical Cyclone Forecasting Framework (TC-FF), to probabilistically forecast compound flooding induced by tropical cyclones, considering uncertainties in track, forward speed, and wind speed and/or intensity. The open-source method accounts for all major relevant physical drivers, including tide, surge, and rainfall, and considers TC uncertainties through Gaussian error distributions and autoregressive techniques. The tool creates temporally and spatially varying wind fields to force a computationally efficient compound-flood model, allowing for the computation of probabilistic wind and flood hazard maps for any oceanic basin in the world as it does not require detailed information on the distribution of historical errors. A comparison of TC-FF and JTWC operational ensembles, both based on DeMaria et al. (2009), revealed minor differences of <10 %, suggesting that TC-FF can be employed as an alternative, for example, in data-scarce environments. The method was applied to Cyclone Idai in Mozambique. The underlying physical model showed reliable skill in terms of tidal propagation, reproducing the storm surge generation during landfall and flooding near the city of Beira (success index of 0.59). The method was successfully applied to forecasting the impact of Idai with different lead times. The case study analyzed needed at least 200 ensemble members to get reliable water levels and flood results 3 d before landfall (<1 % flood probability error and <20 cm sampling errors). Results showed the sensitivity of forecasting, especially with increasing lead times, highlighting the importance of accounting for cyclone variability in decision-making and risk management.
{"title":"Accounting for uncertainties in forecasting tropical-cyclone-induced compound flooding","authors":"Kees Nederhoff, M. van Ormondt, Jay Veeramony, A. V. van Dongeren, J. Antolínez, T. Leijnse, D. Roelvink","doi":"10.5194/gmd-17-1789-2024","DOIUrl":"https://doi.org/10.5194/gmd-17-1789-2024","url":null,"abstract":"Abstract. Tropical-cyclone impacts can have devastating effects on the population, infrastructure, and natural habitats. However, predicting these impacts is difficult due to the inherent uncertainties in the storm track and intensity. In addition, due to computational constraints, both the relevant ocean physics and the uncertainties in meteorological forcing are only partly accounted for. This paper presents a new method, called the Tropical Cyclone Forecasting Framework (TC-FF), to probabilistically forecast compound flooding induced by tropical cyclones, considering uncertainties in track, forward speed, and wind speed and/or intensity. The open-source method accounts for all major relevant physical drivers, including tide, surge, and rainfall, and considers TC uncertainties through Gaussian error distributions and autoregressive techniques. The tool creates temporally and spatially varying wind fields to force a computationally efficient compound-flood model, allowing for the computation of probabilistic wind and flood hazard maps for any oceanic basin in the world as it does not require detailed information on the distribution of historical errors. A comparison of TC-FF and JTWC operational ensembles, both based on DeMaria et al. (2009), revealed minor differences of <10 %, suggesting that TC-FF can be employed as an alternative, for example, in data-scarce environments. The method was applied to Cyclone Idai in Mozambique. The underlying physical model showed reliable skill in terms of tidal propagation, reproducing the storm surge generation during landfall and flooding near the city of Beira (success index of 0.59). The method was successfully applied to forecasting the impact of Idai with different lead times. The case study analyzed needed at least 200 ensemble members to get reliable water levels and flood results 3 d before landfall (<1 % flood probability error and <20 cm sampling errors). Results showed the sensitivity of forecasting, especially with increasing lead times, highlighting the importance of accounting for cyclone variability in decision-making and risk management.\u0000","PeriodicalId":12799,"journal":{"name":"Geoscientific Model Development","volume":null,"pages":null},"PeriodicalIF":5.1,"publicationDate":"2024-02-29","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140408807","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2024-02-29DOI: 10.5194/gmd-17-1831-2024
Ngoc B. Trinh, M. Herrmann, C. Ulses, P. Marsaleix, Thomas Duhaut, Thai To Duy, C. Estournel, R. K. Shearman
Abstract. The South China Sea throughflow (SCSTF) connects the South China Sea (SCS) with neighboring seas and oceans, transferring surface water of the global thermohaline circulation between the Pacific and Indian oceans. A configuration of the SYMPHONIE ocean model at high resolution (4 km) and including an explicit representation of tides is implemented over this region, and a simulation is analyzed over 2010–2018. Comparisons with in situ and satellite data and other available simulations at coarser resolution show the good performance of the model and the relevance of the high resolution for reproducing the spatial and temporal variability of the characteristics of surface dynamics and water masses over the SCS. The added value of an online computation of each term of the water, heat, and salt SCS budgets (surface, lateral oceanic and river fluxes, and internal variations) is also quantitatively demonstrated: important discards are obtained with offline computation, with relative biases of ∼40 % for lateral oceanic inflows and outflows. The SCS water volume budget, including the SCSTF, is analyzed at climatological and seasonal scales. The SCS receives on average a 4.5 Sv yearly water volume input, mainly from the Luzon Strait. It laterally releases this water to neighboring seas, mainly to the Sulu Sea through Mindoro Strait (49 %), to the East China Sea via Taiwan Strait (28 %), and to the Java Sea through Karimata Strait (22 %). The seasonal variability of this water volume budget is driven by lateral interocean exchanges. Surface interocean exchanges, especially at Luzon Strait, are all driven by monsoon winds that favor winter southwestward flows and summer northeastward surface flows. Exchanges through Luzon Strait deep layers show a stable sandwiched structure with vertically alternating inflows and outflows. Last, differences in flux estimates induced by the use of a high-resolution model vs. a low-resolution model are quantified.
{"title":"New insights into the South China Sea throughflow and water budget seasonal cycle: evaluation and analysis of a high-resolution configuration of the ocean model SYMPHONIE version 2.4","authors":"Ngoc B. Trinh, M. Herrmann, C. Ulses, P. Marsaleix, Thomas Duhaut, Thai To Duy, C. Estournel, R. K. Shearman","doi":"10.5194/gmd-17-1831-2024","DOIUrl":"https://doi.org/10.5194/gmd-17-1831-2024","url":null,"abstract":"Abstract. The South China Sea throughflow (SCSTF) connects the South China Sea (SCS) with neighboring seas and oceans, transferring surface water of the global thermohaline circulation between the Pacific and Indian oceans. A configuration of the SYMPHONIE ocean model at high resolution (4 km) and including an explicit representation of tides is implemented over this region, and a simulation is analyzed over 2010–2018. Comparisons with in situ and satellite data and other available simulations at coarser resolution show the good performance of the model and the relevance of the high resolution for reproducing the spatial and temporal variability of the characteristics of surface dynamics and water masses over the SCS. The added value of an online computation of each term of the water, heat, and salt SCS budgets (surface, lateral oceanic and river fluxes, and internal variations) is also quantitatively demonstrated: important discards are obtained with offline computation, with relative biases of ∼40 % for lateral oceanic inflows and outflows. The SCS water volume budget, including the SCSTF, is analyzed at climatological and seasonal scales. The SCS receives on average a 4.5 Sv yearly water volume input, mainly from the Luzon Strait. It laterally releases this water to neighboring seas, mainly to the Sulu Sea through Mindoro Strait (49 %), to the East China Sea via Taiwan Strait (28 %), and to the Java Sea through Karimata Strait (22 %). The seasonal variability of this water volume budget is driven by lateral interocean exchanges. Surface interocean exchanges, especially at Luzon Strait, are all driven by monsoon winds that favor winter southwestward flows and summer northeastward surface flows. Exchanges through Luzon Strait deep layers show a stable sandwiched structure with vertically alternating inflows and outflows. Last, differences in flux estimates induced by the use of a high-resolution model vs. a low-resolution model are quantified.\u0000","PeriodicalId":12799,"journal":{"name":"Geoscientific Model Development","volume":null,"pages":null},"PeriodicalIF":5.1,"publicationDate":"2024-02-29","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140412076","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2024-02-28DOI: 10.5194/gmd-17-1749-2024
Louis Thiry, Long Li, Guillaume Roullet, Etienne Mémin
Abstract. This paper presents MQGeometry, a multi-layer quasi-geostrophic (QG) equation solver for non-rectangular geometries. We advect the potential vorticity (PV) with finite volumes to ensure global PV conservation using a staggered discretization of the PV and stream function (SF). Thanks to this staggering, the PV is defined inside the domain, removing the need to define the PV on the domain boundary. We compute PV fluxes with upwind-biased interpolations whose implicit dissipation replaces the usual explicit (hyper-)viscous dissipation. The discretization presented here does not require tuning of any additional parameter, e.g., additional eddy viscosity. We solve the QG elliptic equation with a fast discrete sine transform spectral solver on rectangular geometry. We extend this fast solver to non-rectangular geometries using the capacitance matrix method. Subsequently, we validate our solver on a vortex-shear instability test case in a circular domain, on a vortex–wall interaction test case, and on an idealized wind-driven double-gyre configuration in an octagonal domain at an eddy-permitting resolution. Finally, we release a concise, efficient, and auto-differentiable PyTorch implementation of our method to facilitate future developments on this new discretization, e.g., machine-learning parameterization or data-assimilation techniques.
{"title":"MQGeometry-1.0: a multi-layer quasi-geostrophic solver on non-rectangular geometries","authors":"Louis Thiry, Long Li, Guillaume Roullet, Etienne Mémin","doi":"10.5194/gmd-17-1749-2024","DOIUrl":"https://doi.org/10.5194/gmd-17-1749-2024","url":null,"abstract":"Abstract. This paper presents MQGeometry, a multi-layer quasi-geostrophic (QG) equation solver for non-rectangular geometries. We advect the potential vorticity (PV) with finite volumes to ensure global PV conservation using a staggered discretization of the PV and stream function (SF). Thanks to this staggering, the PV is defined inside the domain, removing the need to define the PV on the domain boundary. We compute PV fluxes with upwind-biased interpolations whose implicit dissipation replaces the usual explicit (hyper-)viscous dissipation. The discretization presented here does not require tuning of any additional parameter, e.g., additional eddy viscosity. We solve the QG elliptic equation with a fast discrete sine transform spectral solver on rectangular geometry. We extend this fast solver to non-rectangular geometries using the capacitance matrix method. Subsequently, we validate our solver on a vortex-shear instability test case in a circular domain, on a vortex–wall interaction test case, and on an idealized wind-driven double-gyre configuration in an octagonal domain at an eddy-permitting resolution. Finally, we release a concise, efficient, and auto-differentiable PyTorch implementation of our method to facilitate future developments on this new discretization, e.g., machine-learning parameterization or data-assimilation techniques.\u0000","PeriodicalId":12799,"journal":{"name":"Geoscientific Model Development","volume":null,"pages":null},"PeriodicalIF":5.1,"publicationDate":"2024-02-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140421506","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2024-02-28DOI: 10.5194/gmd-17-1765-2024
Marlene Klockmann, Udo von Toussaint, E. Zorita
Abstract. We present a new framework for the reconstruction of climate indices based on proxy data such as tree rings. The framework is based on the supervised learning method Gaussian Process Regression (GPR) and aims at preserving the amplitude of past climate variability. It can adequately handle noise-contaminated proxies and variable proxy availability over time. To this end, the GPR is formulated in a modified input space, termed here embedding space. We test the new framework for the reconstruction of the Atlantic multi-decadal variability (AMV) in a controlled environment with pseudo-proxies derived from coupled climate-model simulations. In this test environment, the GPR outperforms benchmark reconstructions based on multi-linear principal component regression. On AMV-relevant timescales, i.e. multi-decadal, the GPR is able to reconstruct the true amplitude of variability even if the proxies contain a realistic non-climatic noise signal and become sparser back in time. Thus, we conclude that the embedded GPR framework is a highly promising tool for climate-index reconstructions.
{"title":"Towards variance-conserving reconstructions of climate indices with Gaussian process regression in an embedding space","authors":"Marlene Klockmann, Udo von Toussaint, E. Zorita","doi":"10.5194/gmd-17-1765-2024","DOIUrl":"https://doi.org/10.5194/gmd-17-1765-2024","url":null,"abstract":"Abstract. We present a new framework for the reconstruction of climate indices based on proxy data such as tree rings. The framework is based on the supervised learning method Gaussian Process Regression (GPR) and aims at preserving the amplitude of past climate variability. It can adequately handle noise-contaminated proxies and variable proxy availability over time. To this end, the GPR is formulated in a modified input space, termed here embedding space. We test the new framework for the reconstruction of the Atlantic multi-decadal variability (AMV) in a controlled environment with pseudo-proxies derived from coupled climate-model simulations. In this test environment, the GPR outperforms benchmark reconstructions based on multi-linear principal component regression. On AMV-relevant timescales, i.e. multi-decadal, the GPR is able to reconstruct the true amplitude of variability even if the proxies contain a realistic non-climatic noise signal and become sparser back in time. Thus, we conclude that the embedded GPR framework is a highly promising tool for climate-index reconstructions.\u0000","PeriodicalId":12799,"journal":{"name":"Geoscientific Model Development","volume":null,"pages":null},"PeriodicalIF":5.1,"publicationDate":"2024-02-28","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140418225","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2024-02-27DOI: 10.5194/gmd-17-1729-2024
Aaron A. Naidoo-Bagwell, F. Monteiro, K. Hendry, Scott Burgan, Jamie D. Wilson, B. Ward, A. Ridgwell, D. Conley
Abstract. We extend the ecological component (ECOGEM) of the carbon-centric Grid-Enabled Integrated Earth system model (cGEnIE) to include a diatom functional group. ECOGEM represents plankton community dynamics via a spectrum of ecophysiological traits originally based on size and plankton food web (phyto- and zooplankton; EcoGEnIE 1.0), which we developed here to account for a diatom functional group (EcoGEnIE 1.1). We tuned EcoGEnIE 1.1, exploring a range of ecophysiological parameter values specific to phytoplankton, including diatom growth and survival (18 parameters over 550 runs) to achieve best fits to observations of diatom biogeography and size class distribution as well as to global ocean nutrient and dissolved oxygen distributions. This, in conjunction with a previously developed representation of opal dissolution and an updated representation of the ocean iron cycle in the water column, resulted in an improved distribution of dissolved oxygen in the water column relative to the previous EcoGEnIE 1.0, with global export production (7.4 Gt C yr−1) now closer to previous estimates. Simulated diatom biogeography is characterised by larger size classes dominating at high latitudes, notably in the Southern Ocean, and smaller size classes dominating at lower latitudes. Overall, diatom biological productivity accounts for ∼20 % of global carbon biomass in the model, with diatoms outcompeting other phytoplankton functional groups when dissolved silica is available due to their faster maximum photosynthetic rates and reduced palatability to grazers. Adding a diatom functional group provides the cGEnIE Earth system model with an extended capability to explore ecological dynamics and their influence on ocean biogeochemistry.
{"title":"A diatom extension to the cGEnIE Earth system model – EcoGEnIE 1.1","authors":"Aaron A. Naidoo-Bagwell, F. Monteiro, K. Hendry, Scott Burgan, Jamie D. Wilson, B. Ward, A. Ridgwell, D. Conley","doi":"10.5194/gmd-17-1729-2024","DOIUrl":"https://doi.org/10.5194/gmd-17-1729-2024","url":null,"abstract":"Abstract. We extend the ecological component (ECOGEM) of the carbon-centric Grid-Enabled Integrated Earth system model (cGEnIE) to include a diatom functional group. ECOGEM represents plankton community dynamics via a spectrum of ecophysiological traits originally based on size and plankton food web (phyto- and zooplankton; EcoGEnIE 1.0), which we developed here to account for a diatom functional group (EcoGEnIE 1.1). We tuned EcoGEnIE 1.1, exploring a range of ecophysiological parameter values specific to phytoplankton, including diatom growth and survival (18 parameters over 550 runs) to achieve best fits to observations of diatom biogeography and size class distribution as well as to global ocean nutrient and dissolved oxygen distributions. This, in conjunction with a previously developed representation of opal dissolution and an updated representation of the ocean iron cycle in the water column, resulted in an improved distribution of dissolved oxygen in the water column relative to the previous EcoGEnIE 1.0, with global export production (7.4 Gt C yr−1) now closer to previous estimates. Simulated diatom biogeography is characterised by larger size classes dominating at high latitudes, notably in the Southern Ocean, and smaller size classes dominating at lower latitudes. Overall, diatom biological productivity accounts for ∼20 % of global carbon biomass in the model, with diatoms outcompeting other phytoplankton functional groups when dissolved silica is available due to their faster maximum photosynthetic rates and reduced palatability to grazers. Adding a diatom functional group provides the cGEnIE Earth system model with an extended capability to explore ecological dynamics and their influence on ocean biogeochemistry.\u0000","PeriodicalId":12799,"journal":{"name":"Geoscientific Model Development","volume":null,"pages":null},"PeriodicalIF":5.1,"publicationDate":"2024-02-27","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140425166","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2024-02-26DOI: 10.5194/gmd-17-1709-2024
M. Butzin, Ying Ye, Christoph Völker, Özgür Gürses, J. Hauck, Peter Köhler
Abstract. In this paper we describe the implementation of the carbon isotopes 13C and 14C (radiocarbon) into the marine biogeochemistry model REcoM3. The implementation is tested in long-term equilibrium simulations where REcoM3 is coupled with the ocean general circulation model FESOM2.1, applying a low-resolution configuration and idealized climate forcing. Focusing on the carbon-isotopic composition of dissolved inorganic carbon (δ13CDIC and Δ14CDIC), our model results are largely consistent with reconstructions for the pre-anthropogenic period. Our simulations also exhibit discrepancies, e.g. in upwelling regions and the interior of the North Pacific. Some of these differences are due to the limitations of our ocean circulation model setup, which results in a rather shallow meridional overturning circulation. We additionally study the accuracy of two simplified modelling approaches for dissolved inorganic 14C, which are faster (15 % and about a factor of five, respectively) than the complete consideration of the marine radiocarbon cycle. The accuracy of both simplified approaches is better than 5 %, which should be sufficient for most studies of Δ14CDIC.
{"title":"Carbon isotopes in the marine biogeochemistry model FESOM2.1-REcoM3","authors":"M. Butzin, Ying Ye, Christoph Völker, Özgür Gürses, J. Hauck, Peter Köhler","doi":"10.5194/gmd-17-1709-2024","DOIUrl":"https://doi.org/10.5194/gmd-17-1709-2024","url":null,"abstract":"Abstract. In this paper we describe the implementation of the carbon isotopes 13C and 14C (radiocarbon) into the marine biogeochemistry model REcoM3. The implementation is tested in long-term equilibrium simulations where REcoM3 is coupled with the ocean general circulation model FESOM2.1, applying a low-resolution configuration and idealized climate forcing. Focusing on the carbon-isotopic composition of dissolved inorganic carbon (δ13CDIC and Δ14CDIC), our model results are largely consistent with reconstructions for the pre-anthropogenic period. Our simulations also exhibit discrepancies, e.g. in upwelling regions and the interior of the North Pacific. Some of these differences are due to the limitations of our ocean circulation model setup, which results in a rather shallow meridional overturning circulation. We additionally study the accuracy of two simplified modelling approaches for dissolved inorganic 14C, which are faster (15 % and about a factor of five, respectively) than the complete consideration of the marine radiocarbon cycle. The accuracy of both simplified approaches is better than 5 %, which should be sufficient for most studies of Δ14CDIC.\u0000","PeriodicalId":12799,"journal":{"name":"Geoscientific Model Development","volume":null,"pages":null},"PeriodicalIF":5.1,"publicationDate":"2024-02-26","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140430469","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2024-02-26DOI: 10.5194/gmd-17-1667-2024
Ferdinand Briegel, Jonas Wehrle, D. Schindler, A. Christen
Abstract. As the frequency and intensity of heatwaves will continue to increase in the future, accurate and high-resolution mapping and forecasting of human outdoor thermal comfort in urban environments are of great importance. This study presents a machine-learning-based outdoor thermal comfort model with a good trade-off between computational cost, complexity, and accuracy compared to common numerical urban climate models. The machine learning approach is basically an emulation of different numerical urban climate models. The final model consists of four submodels that predict air temperature, relative humidity, wind speed, and mean radiant temperature based on meteorological forcing and geospatial data on building forms, land cover, and vegetation. These variables are then combined into a thermal index (universal thermal climate index – UTCI). All four submodel predictions and the final model output are evaluated using street-level measurements from a dense urban sensor network in Freiburg, Germany. The final model has a mean absolute error of 2.3 K. Based on a city-wide simulation for Freiburg, we demonstrate that the model is fast and versatile enough to simulate multiple years at hourly time steps to predict street-level UTCI at 1 m spatial resolution for an entire city. Simulations indicate that neighbourhood-averaged thermal comfort conditions vary widely between neighbourhoods, even if they are attributed to the same local climate zones, for example, due to differences in age and degree of urban vegetation. Simulations also show contrasting differences in the location of hotspots during the day and at night.
{"title":"High-resolution multi-scaling of outdoor human thermal comfort and its intra-urban variability based on machine learning","authors":"Ferdinand Briegel, Jonas Wehrle, D. Schindler, A. Christen","doi":"10.5194/gmd-17-1667-2024","DOIUrl":"https://doi.org/10.5194/gmd-17-1667-2024","url":null,"abstract":"Abstract. As the frequency and intensity of heatwaves will continue to increase in the future, accurate and high-resolution mapping and forecasting of human outdoor thermal comfort in urban environments are of great importance. This study presents a machine-learning-based outdoor thermal comfort model with a good trade-off between computational cost, complexity, and accuracy compared to common numerical urban climate models. The machine learning approach is basically an emulation of different numerical urban climate models. The final model consists of four submodels that predict air temperature, relative humidity, wind speed, and mean radiant temperature based on meteorological forcing and geospatial data on building forms, land cover, and vegetation. These variables are then combined into a thermal index (universal thermal climate index – UTCI). All four submodel predictions and the final model output are evaluated using street-level measurements from a dense urban sensor network in Freiburg, Germany. The final model has a mean absolute error of 2.3 K. Based on a city-wide simulation for Freiburg, we demonstrate that the model is fast and versatile enough to simulate multiple years at hourly time steps to predict street-level UTCI at 1 m spatial resolution for an entire city. Simulations indicate that neighbourhood-averaged thermal comfort conditions vary widely between neighbourhoods, even if they are attributed to the same local climate zones, for example, due to differences in age and degree of urban vegetation. Simulations also show contrasting differences in the location of hotspots during the day and at night.\u0000","PeriodicalId":12799,"journal":{"name":"Geoscientific Model Development","volume":null,"pages":null},"PeriodicalIF":5.1,"publicationDate":"2024-02-26","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140429878","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2024-02-26DOI: 10.5194/gmd-17-1689-2024
Sven Karsten, H. Radtke, M. Gröger, H. T. M. Ho-Hagemann, Hossein Mashayekh, T. Neumann, H. E. M. Meier
Abstract. In this article the development of a high-resolution Earth System Model (ESM) for the Baltic Sea region is described. In contrast to conventional coupling approaches, the presented model features an additional (technical) component, the flux calculator, which calculates fluxes between the model components on a common exchange grid. This approach naturally ensures conservation of exchanged quantities, a locally consistent treatment of the fluxes, and facilitates interchanging model components in a straightforward manner. The main purpose of this model is to downscale global reanalysis or climate model data to the Baltic Sea region as typically, global model grids are too coarse to resolve the region of interest sufficiently. The regional ESM consists of the Modular Ocean Model 5 (MOM5) for the ocean and the COSMO model in CLimate Mode (CCLM, version 5.0_clm3) for the atmosphere. The bi-directional ocean–atmosphere coupling allows for a realistic air–sea feedback that outperforms the traditional approach of using uncoupled standalone models, as typically pursued with the EURO-CORDEX protocol. In order to address marine environmental problems (e.g., eutrophication and oxygen depletion), the ocean model is internally coupled with the marine biogeochemistry model, ERGOM, set up for the Baltic Sea's hydrographic conditions. The regional ESM can be used for various scientific questions such as climate sensitivity experiments, reconstruction of ocean dynamics, study of past climates, and natural variability, as well as investigation of ocean–atmosphere interactions. Therefore, it can serve for a better understanding of natural processes via attribution experiments that relate observed changes to mechanistic causes.
{"title":"Flux coupling approach on an exchange grid for the IOW Earth System Model (version 1.04.00) of the Baltic Sea region","authors":"Sven Karsten, H. Radtke, M. Gröger, H. T. M. Ho-Hagemann, Hossein Mashayekh, T. Neumann, H. E. M. Meier","doi":"10.5194/gmd-17-1689-2024","DOIUrl":"https://doi.org/10.5194/gmd-17-1689-2024","url":null,"abstract":"Abstract. In this article the development of a high-resolution Earth System Model (ESM) for the Baltic Sea region is described. In contrast to conventional coupling approaches, the presented model features an additional (technical) component, the flux calculator, which calculates fluxes between the model components on a common exchange grid. This approach naturally ensures conservation of exchanged quantities, a locally consistent treatment of the fluxes, and facilitates interchanging model components in a straightforward manner. The main purpose of this model is to downscale global reanalysis or climate model data to the Baltic Sea region as typically, global model grids are too coarse to resolve the region of interest sufficiently. The regional ESM consists of the Modular Ocean Model 5 (MOM5) for the ocean and the COSMO model in CLimate Mode (CCLM, version 5.0_clm3) for the atmosphere. The bi-directional ocean–atmosphere coupling allows for a realistic air–sea feedback that outperforms the traditional approach of using uncoupled standalone models, as typically pursued with the EURO-CORDEX protocol. In order to address marine environmental problems (e.g., eutrophication and oxygen depletion), the ocean model is internally coupled with the marine biogeochemistry model, ERGOM, set up for the Baltic Sea's hydrographic conditions. The regional ESM can be used for various scientific questions such as climate sensitivity experiments, reconstruction of ocean dynamics, study of past climates, and natural variability, as well as investigation of ocean–atmosphere interactions. Therefore, it can serve for a better understanding of natural processes via attribution experiments that relate observed changes to mechanistic causes.\u0000","PeriodicalId":12799,"journal":{"name":"Geoscientific Model Development","volume":null,"pages":null},"PeriodicalIF":5.1,"publicationDate":"2024-02-26","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140428698","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2024-02-26DOI: 10.5194/gmd-17-1651-2024
Zheqi Shen, Yihao Chen, Xiaojing Li, Xunshu Song
Abstract. This study investigates parameter estimation (PE) to enhance climate forecasts of a coupled general circulation model by adjusting the background vertical diffusivity coefficients in its ocean component. These parameters were initially identified through sensitivity experiments and subsequently estimated by assimilating the sea surface temperature and temperature–salinity profiles. This study expands the coupled data assimilation system of the Community Earth System Model (CESM) and the ensemble adjustment Kalman filter (EAKF) to enable parameter estimation. PE experiments were performed to establish balanced initial states and adjusted parameters for forecasting the El Niño–Southern Oscillation (ENSO). Comparing the model states between the PE experiment and a state estimation (SE) experiment revealed that PE can significantly reduce the uncertainty of these parameters and improve the quality of analysis. The forecasts obtained from PE and SE experiments further validate that PE has the potential to improve the forecast skill for ENSO.
摘要本研究调查了参数估计(PE),通过调整海洋部分的背景垂直扩散系数来增强耦合大气环流模式的气候预测。这些参数最初是通过敏感性实验确定的,随后通过同化海面温度和温度-盐度剖面进行估算。这项研究扩展了群落地球系统模式(CESM)和集合调整卡尔曼滤波器(EAKF)的耦合数据同化系统,以实现参数估计。为建立平衡的初始状态和预报厄尔尼诺-南方涛动(ENSO)的调整参数,进行了 PE 实验。通过比较 PE 实验和状态估计(SE)实验之间的模型状态,发现 PE 可以显著降低这些参数的不确定性,提高分析质量。PE 和 SE 实验得出的预测结果进一步验证了 PE 有可能提高厄尔尼诺/南方涛动的预测技能。
{"title":"Parameter estimation for ocean background vertical diffusivity coefficients in the Community Earth System Model (v1.2.1) and its impact on El Niño–Southern Oscillation forecasts","authors":"Zheqi Shen, Yihao Chen, Xiaojing Li, Xunshu Song","doi":"10.5194/gmd-17-1651-2024","DOIUrl":"https://doi.org/10.5194/gmd-17-1651-2024","url":null,"abstract":"Abstract. This study investigates parameter estimation (PE) to enhance climate forecasts of a coupled general circulation model by adjusting the background vertical diffusivity coefficients in its ocean component. These parameters were initially identified through sensitivity experiments and subsequently estimated by assimilating the sea surface temperature and temperature–salinity profiles. This study expands the coupled data assimilation system of the Community Earth System Model (CESM) and the ensemble adjustment Kalman filter (EAKF) to enable parameter estimation. PE experiments were performed to establish balanced initial states and adjusted parameters for forecasting the El Niño–Southern Oscillation (ENSO). Comparing the model states between the PE experiment and a state estimation (SE) experiment revealed that PE can significantly reduce the uncertainty of these parameters and improve the quality of analysis. The forecasts obtained from PE and SE experiments further validate that PE has the potential to improve the forecast skill for ENSO.\u0000","PeriodicalId":12799,"journal":{"name":"Geoscientific Model Development","volume":null,"pages":null},"PeriodicalIF":5.1,"publicationDate":"2024-02-26","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"140428880","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}