Lara Wallberg, Laura Suarez‐Gutierrez, Daniela Matei, Wolfgang A. Müller
Abstract. The internal variability of European summer temperatures has been linked to various mechanisms on seasonal to sub- and multi-decadal timescales. We find that sub-decadal timescales dominate summer temperature variability over large parts of the continent and determine mechanisms controlling extremely warm summers on sub-decadal timescales. We show that the sub-decadal warm phases of bandpass-filtered European summer temperatures, hereinafter referred to as extremely warm European summers, are related to a strengthening of the North Atlantic Ocean subtropical gyre, an increase in meridional heat transport, and an accumulation of ocean heat content in the North Atlantic several years prior to the extreme summer. This ocean warming affects the ocean–atmosphere heat fluxes, leading to a weakening and northward displacement of the jet stream and increased probability of occurrence of high-pressure systems over Scandinavia. Thus, our findings link the occurrence of extremely warm European summers to the accumulation of heat in the North Atlantic Ocean and provide the potential to improve the predictability of extremely warm summers several years ahead, which is of great societal interest.
{"title":"Extremely warm European summers preceded by sub-decadal North Atlantic ocean heat accumulation","authors":"Lara Wallberg, Laura Suarez‐Gutierrez, Daniela Matei, Wolfgang A. Müller","doi":"10.5194/esd-15-1-2024","DOIUrl":"https://doi.org/10.5194/esd-15-1-2024","url":null,"abstract":"Abstract. The internal variability of European summer temperatures has been linked to various mechanisms on seasonal to sub- and multi-decadal timescales. We find that sub-decadal timescales dominate summer temperature variability over large parts of the continent and determine mechanisms controlling extremely warm summers on sub-decadal timescales. We show that the sub-decadal warm phases of bandpass-filtered European summer temperatures, hereinafter referred to as extremely warm European summers, are related to a strengthening of the North Atlantic Ocean subtropical gyre, an increase in meridional heat transport, and an accumulation of ocean heat content in the North Atlantic several years prior to the extreme summer. This ocean warming affects the ocean–atmosphere heat fluxes, leading to a weakening and northward displacement of the jet stream and increased probability of occurrence of high-pressure systems over Scandinavia. Thus, our findings link the occurrence of extremely warm European summers to the accumulation of heat in the North Atlantic Ocean and provide the potential to improve the predictability of extremely warm summers several years ahead, which is of great societal interest.\u0000","PeriodicalId":48931,"journal":{"name":"Earth System Dynamics","volume":"59 25","pages":""},"PeriodicalIF":7.3,"publicationDate":"2024-01-09","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"139441563","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2023-12-22DOI: 10.5194/esd-14-1363-2023
Yinglin Tian, Deyu Zhong, Sarosh Alam Ghausi, Guangqian Wang, A. Kleidon
Abstract. A dominant term in the surface energy balance and central to global warming is downwelling longwave radiation (Rld). It is influenced by radiative properties of the atmospheric column, in particular by greenhouse gases, water vapor, clouds, and differences in atmospheric heat storage. We use the semi-empirical equation derived by Brutsaert (1975) to identify the leading terms responsible for the spatial–temporal climatological variations in Rld. This equation requires only near-surface observations of air temperature and humidity. We first evaluated this equation and its extension by Crawford and Duchon (1999) with observations from FLUXNET, the NASA-CERES dataset, and the ERA5 reanalysis. We found a strong spatiotemporal correlation between estimated Rld and the datasets above, with r2 ranging from 0.87 to 0.98 across the datasets for clear-sky and all-sky conditions. We then used the equations to show that changes in lower atmospheric heat storage explain more than 95 % and around 73 % of diurnal range and seasonal variations in Rld, respectively, with the regional contribution decreasing with latitude. Seasonal changes in the emissivity of the atmosphere play a second role, which is controlled by anomalies in cloud cover at high latitudes but dominated by water vapor changes at midlatitudes and subtropics, especially over monsoon regions. We also found that as aridity increases over the region, the contributions from changes in emissivity and lower atmospheric heat storage tend to offset each other (−40 and 20–30 W m−2, respectively), explaining the relatively small decrease in Rld with aridity (−(10–20) W m−2). These equations thus provide a solid physical basis for understanding the spatiotemporal variability of surface downwelling longwave radiation. This should help us to better understand and interpret climatological changes, such as those associated with extreme events and global warming.
{"title":"Understanding variations in downwelling longwave radiation using Brutsaert's equation","authors":"Yinglin Tian, Deyu Zhong, Sarosh Alam Ghausi, Guangqian Wang, A. Kleidon","doi":"10.5194/esd-14-1363-2023","DOIUrl":"https://doi.org/10.5194/esd-14-1363-2023","url":null,"abstract":"Abstract. A dominant term in the surface energy balance and central to global warming is downwelling longwave radiation (Rld). It is influenced by radiative properties of the atmospheric column, in particular by greenhouse gases, water vapor, clouds, and differences in atmospheric heat storage. We use the semi-empirical equation derived by Brutsaert (1975) to identify the leading terms responsible for the spatial–temporal climatological variations in Rld. This equation requires only near-surface observations of air temperature and humidity. We first evaluated this equation and its extension by Crawford and Duchon (1999) with observations from FLUXNET, the NASA-CERES dataset, and the ERA5 reanalysis. We found a strong spatiotemporal correlation between estimated Rld and the datasets above, with r2 ranging from 0.87 to 0.98 across the datasets for clear-sky and all-sky conditions. We then used the equations to show that changes in lower atmospheric heat storage explain more than 95 % and around 73 % of diurnal range and seasonal variations in Rld, respectively, with the regional contribution decreasing with latitude. Seasonal changes in the emissivity of the atmosphere play a second role, which is controlled by anomalies in cloud cover at high latitudes but dominated by water vapor changes at midlatitudes and subtropics, especially over monsoon regions. We also found that as aridity increases over the region, the contributions from changes in emissivity and lower atmospheric heat storage tend to offset each other (−40 and 20–30 W m−2, respectively), explaining the relatively small decrease in Rld with aridity (−(10–20) W m−2). These equations thus provide a solid physical basis for understanding the spatiotemporal variability of surface downwelling longwave radiation. This should help us to better understand and interpret climatological changes, such as those associated with extreme events and global warming.\u0000","PeriodicalId":48931,"journal":{"name":"Earth System Dynamics","volume":"5 2","pages":""},"PeriodicalIF":7.3,"publicationDate":"2023-12-22","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"138994193","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2023-12-21DOI: 10.5194/esd-14-1333-2023
Y. Quilcaille, L. Gudmundsson, S. Seneviratne
Abstract. Climate emulators are models calibrated on Earth system models (ESMs) to replicate their behavior. Thanks to their low computational cost, these tools are becoming increasingly important to accelerate the exploration of emission scenarios and the coupling of climate information to other models. However, the emulation of regional climate extremes and water cycle variables has remained challenging. The MESMER emulator was recently expanded to represent regional temperature extremes in the new “MESMER-X” version, which is targeted at impact-related variables, including extremes. This paper presents a further expansion of MESMER-X to represent indices related to fire weather and soil moisture. Given a trajectory of global mean temperature, the extended emulator generates spatially resolved realizations for the seasonal average of the Canadian Fire Weather Index (FWI), the number of days with extreme fire weather, the annual average of the soil moisture, and the annual minimum of the monthly average soil moisture. For each ESM, the emulations mimic the statistical distributions and the spatial patterns of these indicators. For each of the four variables considered, we evaluate the performances of the emulations by calculating how much their quantiles deviate from those of the ESMs. Given how it performs over a large range of annual indicators, we argue that this framework can be expanded to further variables. Overall, the now expanded MESMER-X emulator can emulate several climate variables, including climate extremes and soil moisture availability, and is a useful tool for the exploration of regional climate changes and their impacts.
{"title":"Extending MESMER-X: a spatially resolved Earth system model emulator for fire weather and soil moisture","authors":"Y. Quilcaille, L. Gudmundsson, S. Seneviratne","doi":"10.5194/esd-14-1333-2023","DOIUrl":"https://doi.org/10.5194/esd-14-1333-2023","url":null,"abstract":"Abstract. Climate emulators are models calibrated on Earth system models (ESMs) to replicate their behavior. Thanks to their low computational cost, these tools are becoming increasingly important to accelerate the exploration of emission scenarios and the coupling of climate information to other models. However, the emulation of regional climate extremes and water cycle variables has remained challenging. The MESMER emulator was recently expanded to represent regional temperature extremes in the new “MESMER-X” version, which is targeted at impact-related variables, including extremes. This paper presents a further expansion of MESMER-X to represent indices related to fire weather and soil moisture. Given a trajectory of global mean temperature, the extended emulator generates spatially resolved realizations for the seasonal average of the Canadian Fire Weather Index (FWI), the number of days with extreme fire weather, the annual average of the soil moisture, and the annual minimum of the monthly average soil moisture. For each ESM, the emulations mimic the statistical distributions and the spatial patterns of these indicators. For each of the four variables considered, we evaluate the performances of the emulations by calculating how much their quantiles deviate from those of the ESMs. Given how it performs over a large range of annual indicators, we argue that this framework can be expanded to further variables. Overall, the now expanded MESMER-X emulator can emulate several climate variables, including climate extremes and soil moisture availability, and is a useful tool for the exploration of regional climate changes and their impacts.\u0000","PeriodicalId":48931,"journal":{"name":"Earth System Dynamics","volume":"56 36","pages":""},"PeriodicalIF":7.3,"publicationDate":"2023-12-21","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"138949352","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2023-12-20DOI: 10.5194/esd-14-1317-2023
Chunhua Shen, John C. Moore, H. Kuswanto, Liyun Zhao
Abstract. The Indonesia Throughflow (ITF) is the only low-latitude channel between the Pacific and Indian oceans, and its variability has important effects on global climate and biogeochemical cycles. Climate models consistently predict a decline in ITF transport under global warming, but it has not yet been examined under solar geoengineering scenarios. We use standard parameterized methods for estimating the ITF – the Amended Island Rule and buoyancy forcing – to investigate the ITF under the SSP2-4.5 and SSP5-8.5 greenhouse gas scenarios and the geoengineering experiments G6solar and G6sulfur, which reduce net global mean radiative forcing from SSP5-8.5 levels to SSP2-4.5 levels using solar dimming and sulfate aerosol injection strategies, respectively. Six-model ensemble-mean projections for 2080–2100 show reductions of 19 % under the G6solar scenario and 28 % under the G6sulfur scenario relative to the historical (1980–2014) ITF, which should be compared with reductions of 23 % and 27 % under SSP2-4.5 and SSP5-8.5. Despite standard deviations amounting to 5 %–8 % for each scenario, all scenarios are significantly different from each other (p<0.05) when the whole 2020–2100 simulation period is considered. Thus, significant weakening of the ITF occurs under all scenarios, but G6solar more closely approximates SSP2-4.5 than G6sulfur does. In contrast with the other three scenarios, which show only reductions in forcing due to ocean upwelling, the G6sulfur experiment shows a large reduction in ocean surface wind stress forcing accounting for 47 % (38 %–65 % across the model range) of the decline in wind + upwelling-driven ITF transport. There are also reductions in deep-sea upwelling in extratropical western boundary currents.
{"title":"The Indonesian Throughflow circulation under solar geoengineering","authors":"Chunhua Shen, John C. Moore, H. Kuswanto, Liyun Zhao","doi":"10.5194/esd-14-1317-2023","DOIUrl":"https://doi.org/10.5194/esd-14-1317-2023","url":null,"abstract":"Abstract. The Indonesia Throughflow (ITF) is the only low-latitude channel between the Pacific and Indian oceans, and its variability has important effects on global climate and biogeochemical cycles. Climate models consistently predict a decline in ITF transport under global warming, but it has not yet been examined under solar geoengineering scenarios. We use standard parameterized methods for estimating the ITF – the Amended Island Rule and buoyancy forcing – to investigate the ITF under the SSP2-4.5 and SSP5-8.5 greenhouse gas scenarios and the geoengineering experiments G6solar and G6sulfur, which reduce net global mean radiative forcing from SSP5-8.5 levels to SSP2-4.5 levels using solar dimming and sulfate aerosol injection strategies, respectively. Six-model ensemble-mean projections for 2080–2100 show reductions of 19 % under the G6solar scenario and 28 % under the G6sulfur scenario relative to the historical (1980–2014) ITF, which should be compared with reductions of 23 % and 27 % under SSP2-4.5 and SSP5-8.5. Despite standard deviations amounting to 5 %–8 % for each scenario, all scenarios are significantly different from each other (p<0.05) when the whole 2020–2100 simulation period is considered. Thus, significant weakening of the ITF occurs under all scenarios, but G6solar more closely approximates SSP2-4.5 than G6sulfur does. In contrast with the other three scenarios, which show only reductions in forcing due to ocean upwelling, the G6sulfur experiment shows a large reduction in ocean surface wind stress forcing accounting for 47 % (38 %–65 % across the model range) of the decline in wind + upwelling-driven ITF transport. There are also reductions in deep-sea upwelling in extratropical western boundary currents.","PeriodicalId":48931,"journal":{"name":"Earth System Dynamics","volume":"331 ","pages":""},"PeriodicalIF":7.3,"publicationDate":"2023-12-20","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"139170630","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2023-12-13DOI: 10.5194/esd-14-1295-2023
L. de Mora, R. Swaminathan, Richard P. Allan, Jerry C. Blackford, D. Kelley, Phil Harris, Christopher Jones, Colin G. Jones, S. Liddicoat, R. Parker, T. Quaife, Jeremy Walton, A. Yool
Abstract. We show that the distribution of anthropogenic carbon between the atmosphere, land surface, and ocean differs with the choice of projection scenario even for identical changes in mean global surface temperature. Warming thresholds occur later in lower-CO2-emissions scenarios and with less carbon in the three main reservoirs than in higher-CO2-emissions scenarios. At 2 ∘C of warming, the mean carbon allocation differs by up to 63 PgC between scenarios, which is equivalent to approximately 6 years of the current global total emissions. At the same warming level, higher-CO2-concentration scenarios have a lower combined ocean and land carbon allocation fraction of the total carbon compared to lower-CO2-concentration scenarios. The warming response to CO2, quantified as the equilibrium climate sensitivity, ECS, directly impacts the global warming level exceedance year and hence the carbon allocation. Low-ECS models have more carbon than high-ECS models at a given warming level because the warming threshold occurs later, allowing more emissions to accumulate. These results are important for carbon budgets and mitigation strategies as they impact how much carbon the ocean and land surface could absorb at a given warming level. Carbon budgeting will be key to reducing the impacts of anthropogenic climate change, and these findings could have critical consequences for policies aimed at reaching net zero.
{"title":"Scenario choice impacts carbon allocation projection at global warming levels","authors":"L. de Mora, R. Swaminathan, Richard P. Allan, Jerry C. Blackford, D. Kelley, Phil Harris, Christopher Jones, Colin G. Jones, S. Liddicoat, R. Parker, T. Quaife, Jeremy Walton, A. Yool","doi":"10.5194/esd-14-1295-2023","DOIUrl":"https://doi.org/10.5194/esd-14-1295-2023","url":null,"abstract":"Abstract. We show that the distribution of anthropogenic carbon between the atmosphere, land surface, and ocean differs with the choice of projection scenario even for identical changes in mean global surface temperature. Warming thresholds occur later in lower-CO2-emissions scenarios and with less carbon in the three main reservoirs than in higher-CO2-emissions scenarios. At 2 ∘C of warming, the mean carbon allocation differs by up to 63 PgC between scenarios, which is equivalent to approximately 6 years of the current global total emissions. At the same warming level, higher-CO2-concentration scenarios have a lower combined ocean and land carbon allocation fraction of the total carbon compared to lower-CO2-concentration scenarios. The warming response to CO2, quantified as the equilibrium climate sensitivity, ECS, directly impacts the global warming level exceedance year and hence the carbon allocation. Low-ECS models have more carbon than high-ECS models at a given warming level because the warming threshold occurs later, allowing more emissions to accumulate. These results are important for carbon budgets and mitigation strategies as they impact how much carbon the ocean and land surface could absorb at a given warming level. Carbon budgeting will be key to reducing the impacts of anthropogenic climate change, and these findings could have critical consequences for policies aimed at reaching net zero.\u0000","PeriodicalId":48931,"journal":{"name":"Earth System Dynamics","volume":"42 4","pages":""},"PeriodicalIF":7.3,"publicationDate":"2023-12-13","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"138976369","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2023-12-12DOI: 10.5194/esd-14-1277-2023
Takahito Mitsui, Matteo Willeit, Niklas Boers
Abstract. The glacial–interglacial cycles of the Quaternary exhibit 41 kyr periodicity before the Mid-Pleistocene Transition (MPT) around 1.2–0.8 Myr ago and ∼ 100 kyr periodicity after that. From the viewpoint of dynamical systems, proposed mechanisms generating these periodicities are broadly divided into two types: (i) nonlinear forced responses of a mono- or multi-stable climate system to the astronomical forcing or (ii) synchronization of internal self-sustained oscillations to the astronomical forcing. In this study, we investigate the dynamics of glacial cycles simulated by the Earth system model of intermediate complexity CLIMBER-2 with a fully interactive carbon cycle, which reproduces the MPT under gradual changes in volcanic-CO2 degassing and regolith cover. We report that, in this model, the dominant frequency of glacial cycles is set in line with the principle of synchronization. It is found that the model exhibits self-sustained oscillations in the absence of astronomical forcing. Before the MPT, glacial cycles synchronize to the 41 kyr obliquity cycles because the self-sustained oscillations have periodicity relatively close to 41 kyr. After the MPT the timescale of internal oscillations becomes too long to follow every 41 kyr obliquity cycle, and the oscillations synchronize to the 100 kyr eccentricity cycles that modulate the amplitude of climatic precession. The latter synchronization occurs with the help of the 41 kyr obliquity forcing, which enables some terminations and glaciations to occur robustly at their right timing. We term this phenomenon vibration-enhanced synchronization because of its similarity to the noise-enhanced synchronization known in nonlinear science. While we interpret the dominant periodicities of glacial cycles as the result of synchronization of internal self-sustained oscillations to the astronomical forcing, the Quaternary glacial cycles show facets of both synchronization and forced response.
{"title":"Synchronization phenomena observed in glacial–interglacial cycles simulated in an Earth system model of intermediate complexity","authors":"Takahito Mitsui, Matteo Willeit, Niklas Boers","doi":"10.5194/esd-14-1277-2023","DOIUrl":"https://doi.org/10.5194/esd-14-1277-2023","url":null,"abstract":"Abstract. The glacial–interglacial cycles of the Quaternary exhibit 41 kyr periodicity before the Mid-Pleistocene Transition (MPT) around 1.2–0.8 Myr ago and ∼ 100 kyr periodicity after that. From the viewpoint of dynamical systems, proposed mechanisms generating these periodicities are broadly divided into two types: (i) nonlinear forced responses of a mono- or multi-stable climate system to the astronomical forcing or (ii) synchronization of internal self-sustained oscillations to the astronomical forcing. In this study, we investigate the dynamics of glacial cycles simulated by the Earth system model of intermediate complexity CLIMBER-2 with a fully interactive carbon cycle, which reproduces the MPT under gradual changes in volcanic-CO2 degassing and regolith cover. We report that, in this model, the dominant frequency of glacial cycles is set in line with the principle of synchronization. It is found that the model exhibits self-sustained oscillations in the absence of astronomical forcing. Before the MPT, glacial cycles synchronize to the 41 kyr obliquity cycles because the self-sustained oscillations have periodicity relatively close to 41 kyr. After the MPT the timescale of internal oscillations becomes too long to follow every 41 kyr obliquity cycle, and the oscillations synchronize to the 100 kyr eccentricity cycles that modulate the amplitude of climatic precession. The latter synchronization occurs with the help of the 41 kyr obliquity forcing, which enables some terminations and glaciations to occur robustly at their right timing. We term this phenomenon vibration-enhanced synchronization because of its similarity to the noise-enhanced synchronization known in nonlinear science. While we interpret the dominant periodicities of glacial cycles as the result of synchronization of internal self-sustained oscillations to the astronomical forcing, the Quaternary glacial cycles show facets of both synchronization and forced response.\u0000","PeriodicalId":48931,"journal":{"name":"Earth System Dynamics","volume":"5 6","pages":""},"PeriodicalIF":7.3,"publicationDate":"2023-12-12","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"139008455","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2023-12-05DOI: 10.5194/esd-14-1261-2023
M. Paprota, W. Sulisz
Abstract. A semi-analytical solution to an advection–diffusion equation is coupled with a nonlinear wavemaker model to investigate the effect of strong nonlinearity on wave-induced mixing. The comparisons with weakly nonlinear model predictions reveal that in the case of waves of higher steepness, enhanced mixing affects the subsurface layer of the water column. A fully nonlinear model captures the neglected higher-order terms from a weakly nonlinear solution and provides a reliable estimation of the time-mean velocity field. The corrected wave-induced mass-transport velocity leads to improved estimates of subsurface mixing intensity and ocean surface temperature.
{"title":"The effect of strong nonlinearity on wave-induced vertical mixing","authors":"M. Paprota, W. Sulisz","doi":"10.5194/esd-14-1261-2023","DOIUrl":"https://doi.org/10.5194/esd-14-1261-2023","url":null,"abstract":"Abstract. A semi-analytical solution to an advection–diffusion equation is coupled with a nonlinear wavemaker model to investigate the effect of strong nonlinearity on wave-induced mixing. The comparisons with weakly nonlinear model predictions reveal that in the case of waves of higher steepness, enhanced mixing affects the subsurface layer of the water column. A fully nonlinear model captures the neglected higher-order terms from a weakly nonlinear solution and provides a reliable estimation of the time-mean velocity field. The corrected wave-induced mass-transport velocity leads to improved estimates of subsurface mixing intensity and ocean surface temperature.\u0000","PeriodicalId":48931,"journal":{"name":"Earth System Dynamics","volume":"132 9","pages":""},"PeriodicalIF":7.3,"publicationDate":"2023-12-05","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"138599248","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2023-11-29DOI: 10.5194/esd-14-1239-2023
F. van Oorschot, R. J. van der Ent, M. Hrachowitz, Emanuele Di Carlo, F. Catalano, S. Boussetta, G. Balsamo, Andrea Alessandri
Abstract. Vegetation largely controls land surface–atmosphere interactions. Although vegetation is highly dynamic across spatial and temporal scales, most land surface models currently used for reanalyses and near-term climate predictions do not adequately represent these dynamics. This causes deficiencies in the variability of modeled water and energy states and fluxes from the land surface. In this study we evaluated the effects of integrating spatially and temporally varying land cover and vegetation characteristics derived from satellite observations on modeled evaporation and soil moisture in the Hydrology Tiled ECMWF Scheme for Surface Exchanges over Land (HTESSEL) land surface model. Specifically, we integrated interannually varying land cover from the European Space Agency Climate Change Initiative and interannually varying leaf area index (LAI) from the Copernicus Global Land Services (CGLS). Additionally, satellite data on the fraction of green vegetation cover (FCover) from CGLS were used to formulate and integrate a spatially and temporally varying effective vegetation cover parameterization. The effects of these three implementations on model evaporation fluxes and soil moisture were analyzed using historical offline (land-only) model experiments at the global scale, and model performances were quantified with global observational products of evaporation (E) and near-surface soil moisture (SMs). The interannually varying land cover consistently altered the evaporation and soil moisture in regions with major land cover changes. The interannually varying LAI considerably improved the correlation of SMs and E with respect to the reference data, with the largest improvements in semiarid regions with predominantly low vegetation during the dry season. These improvements are related to the activation of soil moisture–evaporation feedbacks during vegetation-water-stressed periods with interannually varying LAI in combination with interannually varying effective vegetation cover, defined as an exponential function of LAI. The further improved effective vegetation cover parameterization consistently reduced the errors of model effective vegetation cover, and it regionally improved SMs and E. Overall, our study demonstrated that the enhanced vegetation variability consistently improved the near-surface soil moisture and evaporation variability, but the availability of reliable global observational data remains a limitation for complete understanding of the model response. To further explain the improvements found, we developed an interpretation framework for how the model development activates feedbacks between soil moisture, vegetation, and evaporation during vegetation water stress periods.
摘要植被在很大程度上控制着地表与大气的相互作用。虽然植被在空间和时间尺度上具有高度动态性,但目前用于再分析和近期气候预测的大多数陆表模式并不能充分反映这些动态。这就造成了模拟的水和能量状态以及来自地表的通量的变异性存在缺陷。在这项研究中,我们评估了在水文平铺 ECMWF 陆地表面交换方案(HTESSEL)陆地表面模式中整合卫星观测得出的时空变化的土地覆盖和植被特征对模型蒸发和土壤湿度的影响。具体而言,我们整合了欧洲空间局气候变化倡议(Climate Change Initiative)的年际间变化土地覆盖和哥白尼全球陆地服务(CGLS)的年际间变化叶面积指数(LAI)。此外,我们还利用哥白尼全球陆地服务的绿色植被覆盖率(FCover)卫星数据,制定并整合了时空变化的有效植被覆盖参数。利用全球尺度的历史离线(纯陆地)模型实验分析了这三种实现方式对模型蒸发通量和土壤水分的影响,并利用全球蒸发(E)和近地表土壤水分(SMs)观测产品对模型性能进行了量化。年际变化的土地覆被持续改变了土地覆被发生重大变化地区的蒸发量和土壤水分。与参考数据相比,年际变化的 LAI 显著改善了 SMs 和 E 的相关性,在旱季植被主要较低的半干旱地区改善最大。这些改进与年际变化的 LAI 结合年际变化的有效植被覆盖度(定义为 LAI 的指数函数)在植被缺水期激活土壤水分-蒸发反馈有关。总之,我们的研究表明,植被可变性的增强持续改善了近地表土壤水分和蒸发的可变性,但可靠的全球观测数据的可用性仍然是全面理解模型响应的限制因素。为了进一步解释所发现的改进,我们建立了一个解释框架,说明在植被水分胁迫期,模型的发展如何激活土壤水分、植被和蒸发之间的反馈。
{"title":"Interannual land cover and vegetation variability based on remote sensing data in the HTESSEL land surface model: implementation and effects on simulated water dynamics","authors":"F. van Oorschot, R. J. van der Ent, M. Hrachowitz, Emanuele Di Carlo, F. Catalano, S. Boussetta, G. Balsamo, Andrea Alessandri","doi":"10.5194/esd-14-1239-2023","DOIUrl":"https://doi.org/10.5194/esd-14-1239-2023","url":null,"abstract":"Abstract. Vegetation largely controls land surface–atmosphere interactions. Although vegetation is highly dynamic across spatial and temporal scales, most land surface models currently used for reanalyses and near-term climate predictions do not adequately represent these dynamics. This causes deficiencies in the variability of modeled water and energy states and fluxes from the land surface. In this study we evaluated the effects of integrating spatially and temporally varying land cover and vegetation characteristics derived from satellite observations on modeled evaporation and soil moisture in the Hydrology Tiled ECMWF Scheme for Surface Exchanges over Land (HTESSEL) land surface model. Specifically, we integrated interannually varying land cover from the European Space Agency Climate Change Initiative and interannually varying leaf area index (LAI) from the Copernicus Global Land Services (CGLS). Additionally, satellite data on the fraction of green vegetation cover (FCover) from CGLS were used to formulate and integrate a spatially and temporally varying effective vegetation cover parameterization. The effects of these three implementations on model evaporation fluxes and soil moisture were analyzed using historical offline (land-only) model experiments at the global scale, and model performances were quantified with global observational products of evaporation (E) and near-surface soil moisture (SMs). The interannually varying land cover consistently altered the evaporation and soil moisture in regions with major land cover changes. The interannually varying LAI considerably improved the correlation of SMs and E with respect to the reference data, with the largest improvements in semiarid regions with predominantly low vegetation during the dry season. These improvements are related to the activation of soil moisture–evaporation feedbacks during vegetation-water-stressed periods with interannually varying LAI in combination with interannually varying effective vegetation cover, defined as an exponential function of LAI. The further improved effective vegetation cover parameterization consistently reduced the errors of model effective vegetation cover, and it regionally improved SMs and E. Overall, our study demonstrated that the enhanced vegetation variability consistently improved the near-surface soil moisture and evaporation variability, but the availability of reliable global observational data remains a limitation for complete understanding of the model response. To further explain the improvements found, we developed an interpretation framework for how the model development activates feedbacks between soil moisture, vegetation, and evaporation during vegetation water stress periods.","PeriodicalId":48931,"journal":{"name":"Earth System Dynamics","volume":"20 1","pages":""},"PeriodicalIF":7.3,"publicationDate":"2023-11-29","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"139211131","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2023-11-27DOI: 10.5194/esd-14-1211-2023
Chenwei Xiao, S. Zaehle, Hui Yang, J. Wigneron, Christiane Schmullius, Ana Bastos
Abstract. Drought events are projected to become more severe and frequent across many regions in the future, but their impacts will likely differ among ecosystems depending on their ability to maintain functioning during droughts, i.e., ecosystem resistance. Plant species have diverse strategies to cope with drought. As a result, divergent responses of different vegetation types for similar levels of drought severity have been observed. It remains unclear whether such divergence can be explained by different drought duration; co-occurring compounding effects, e.g., heat stress or memory effects; management practices; etc. Here, we provide a global synthesis of vegetation resistance to drought and heat using different proxies for vegetation condition, namely the vegetation optical depth (L-VOD) data from the ESA Soil Moisture and Ocean Salinity (SMOS) passive L-band microwave mission and enhanced vegetation index (EVI) and kernel normalized difference vegetation index (kNDVI) from NASA MODIS. Due to its longer wavelength, L-VOD has the advantage over more commonly used vegetation indices (such as kNDVI, EVI) in that it provides different information on vegetation structure and biomass and suffers from less saturation over dense forests. We apply a linear model accounting for drought and temperature effects to characterize ecosystem resistance by their sensitivity to drought duration and temperature anomalies. We analyze how ecosystem resistance varies with land cover across the globe and investigate the potential effects of forest management and crop irrigation. We compare estimates of ecosystem resistance to drought and heat as retrieved from L-VOD, kNDVI, and EVI products. We find that regions with higher forest fraction show stronger ecosystem resistance to extreme droughts than cropland for all three vegetation proxies. L-VOD indicates that primary forests tend to be more resistant to drought events than secondary forests when controlling for the differences in background climate, but this cannot be detected in EVI and kNDVI. The difference is possibly related to EVI and kNDVI saturation in dense forests. In tropical primary evergreen broadleaf forests, old-growth trees tend to be more resistant to drought than young trees from L-VOD and kNDVI. Irrigation increases the drought resistance of cropland substantially. Forest harvest decreases the drought resistance of forests. Our results suggest that ecosystem resistance can be better monitored using L-VOD in dense forests and highlight the role of forest cover, forest management, and irrigation in determining ecosystem resistance to droughts.
摘要预计未来许多地区的干旱事件将变得更加严重和频繁,但其影响可能因生态系统在干旱期间维持功能的能力(即生态系统抵抗力)而有所不同。植物物种应对干旱的策略多种多样。因此,人们观察到不同植被类型对类似严重干旱程度的不同反应。目前还不清楚这种差异是否可以用不同的干旱持续时间、同时出现的复合效应(如热胁迫或记忆效应)、管理方法等来解释。在此,我们利用不同的植被状况代用指标,即来自欧空局土壤水分和海洋盐度(SMOS)被动 L 波段微波任务的植被光学深度(L-VOD)数据,以及来自美国宇航局 MODIS 的增强植被指数(EVI)和核归一化差异植被指数(kNDVI),对植被的抗旱和抗热能力进行了全球综合分析。由于波长较长,L-VOD 与更常用的植被指数(如 kNDVI、EVI)相比具有优势,它能提供不同的植被结构和生物量信息,并且在密林中的饱和度较低。我们采用一个考虑干旱和温度影响的线性模型,通过生态系统对干旱持续时间和温度异常的敏感性来描述生态系统抗性。我们分析了生态系统抵抗力如何随全球土地覆盖而变化,并研究了森林管理和作物灌溉的潜在影响。我们比较了从 L-VOD、kNDVI 和 EVI 产品中获取的生态系统对干旱和高温的抵抗力估计值。我们发现,就所有三种植被代用指标而言,森林覆盖率较高的地区比耕地对极端干旱的生态系统抵抗力更强。L-VOD 表明,在控制背景气候差异的情况下,原始森林往往比次生林更能抵御干旱事件,但 EVI 和 kNDVI 无法检测到这一点。这种差异可能与密林中的 EVI 和 kNDVI 饱和度有关。在热带原始常绿阔叶林中,从 L-VOD 和 kNDVI 来看,老树往往比幼树更抗旱。灌溉可大幅提高耕地的抗旱性。森林采伐降低了森林的抗旱性。我们的研究结果表明,在密林中使用 L-VOD 可以更好地监测生态系统的抗旱性,并强调了森林覆盖率、森林管理和灌溉在决定生态系统抗旱性中的作用。
{"title":"Land cover and management effects on ecosystem resistance to drought stress","authors":"Chenwei Xiao, S. Zaehle, Hui Yang, J. Wigneron, Christiane Schmullius, Ana Bastos","doi":"10.5194/esd-14-1211-2023","DOIUrl":"https://doi.org/10.5194/esd-14-1211-2023","url":null,"abstract":"Abstract. Drought events are projected to become more severe and frequent across many regions in the future, but their impacts will likely differ among ecosystems depending on their ability to maintain functioning during droughts, i.e., ecosystem resistance. Plant species have diverse strategies to cope with drought. As a result, divergent responses of different vegetation types for similar levels of drought severity have been observed. It remains unclear whether such divergence can be explained by different drought duration; co-occurring compounding effects, e.g., heat stress or memory effects; management practices; etc. Here, we provide a global synthesis of vegetation resistance to drought and heat using different proxies for vegetation condition, namely the vegetation optical depth (L-VOD) data from the ESA Soil Moisture and Ocean Salinity (SMOS) passive L-band microwave mission and enhanced vegetation index (EVI) and kernel normalized difference vegetation index (kNDVI) from NASA MODIS. Due to its longer wavelength, L-VOD has the advantage over more commonly used vegetation indices (such as kNDVI, EVI) in that it provides different information on vegetation structure and biomass and suffers from less saturation over dense forests. We apply a linear model accounting for drought and temperature effects to characterize ecosystem resistance by their sensitivity to drought duration and temperature anomalies. We analyze how ecosystem resistance varies with land cover across the globe and investigate the potential effects of forest management and crop irrigation. We compare estimates of ecosystem resistance to drought and heat as retrieved from L-VOD, kNDVI, and EVI products. We find that regions with higher forest fraction show stronger ecosystem resistance to extreme droughts than cropland for all three vegetation proxies. L-VOD indicates that primary forests tend to be more resistant to drought events than secondary forests when controlling for the differences in background climate, but this cannot be detected in EVI and kNDVI. The difference is possibly related to EVI and kNDVI saturation in dense forests. In tropical primary evergreen broadleaf forests, old-growth trees tend to be more resistant to drought than young trees from L-VOD and kNDVI. Irrigation increases the drought resistance of cropland substantially. Forest harvest decreases the drought resistance of forests. Our results suggest that ecosystem resistance can be better monitored using L-VOD in dense forests and highlight the role of forest cover, forest management, and irrigation in determining ecosystem resistance to droughts.","PeriodicalId":48931,"journal":{"name":"Earth System Dynamics","volume":"15 1","pages":""},"PeriodicalIF":7.3,"publicationDate":"2023-11-27","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"139230030","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Pub Date : 2023-11-17DOI: 10.5194/esd-14-1183-2023
Daniel Boateng, S. Mutz, Armelle Ballian, M. Meijers, K. Methner, S. Botsyun, A. Mulch, T. Ehlers
Abstract. This study presents the simulated response of regional climate and the oxygen isotopic composition of precipitation (δ18Op) to different along-strike topographic evolution scenarios. These simulations are conducted to determine if the previously hypothesized diachronous surface uplift in the Western and Eastern Alps would produce δ18Op signals in the geologic record that are sufficiently large and distinct to be detected using stable isotope paleoaltimetry. We present a series of topographic sensitivity experiments conducted with the water-isotope-tracking atmospheric general circulation model (GCM) ECHAM5-wiso. The topographic scenarios are created from the variation of two free parameters, (1) the elevation of the Western–Central Alps and (2) the elevation of the Eastern Alps. The results indicate Δδ18Op values (i.e., the difference between δ18Op values at the low- and high-elevation sites) of up to −8 ‰ along the strike of the Alps for the diachronous uplift scenarios, primarily due to changes in orographic precipitation and adiabatic lapse rate driven localized changes in near-surface variables. These simulated magnitudes of Δδ18Op values suggest that the expected isotopic signal would be significant enough to be preserved and measured in geologic archives. Moreover, the simulated slight δ18Op differences of 1 ‰–2 ‰ across the low-elevation sites support the use of the δ–δ paleoaltimetry approach and highlight the importance of sampling far-field low-elevation sites to differentiate between the different surface uplift scenarios. The elevation-dependent rate of change in δ18Op (“isotopic lapse rate”) varies depending on the topographic configuration and the extent of the surface uplift. Most of the changes are significant (e.g., −1.04 ‰ km−1 change with slope error of ±0.09 ‰ km−1), while others were within the range of the statistical uncertainties (e.g., −0.15 ‰ km−1 change with slope error of ±0.13 ‰ km−1). The results also highlight the plausible changes in atmospheric circulation patterns and associated changes in moisture transport pathways in response to changes in the topography of the Alps. These large-scale atmospheric dynamics changes can complicate the underlying assumption of stable isotope paleoaltimetry and therefore require integration with paleoclimate modeling to ensure accurate reconstruction of the paleoelevation of the Alps.
{"title":"The effects of diachronous surface uplift of the European Alps on regional climate and the oxygen isotopic composition of precipitation","authors":"Daniel Boateng, S. Mutz, Armelle Ballian, M. Meijers, K. Methner, S. Botsyun, A. Mulch, T. Ehlers","doi":"10.5194/esd-14-1183-2023","DOIUrl":"https://doi.org/10.5194/esd-14-1183-2023","url":null,"abstract":"Abstract. This study presents the simulated response of regional climate and the oxygen isotopic composition of precipitation (δ18Op) to different along-strike topographic evolution scenarios. These simulations are conducted to determine if the previously hypothesized diachronous surface uplift in the Western and Eastern Alps would produce δ18Op signals in the geologic record that are sufficiently large and distinct to be detected using stable isotope paleoaltimetry. We present a series of topographic sensitivity experiments conducted with the water-isotope-tracking atmospheric general circulation model (GCM) ECHAM5-wiso. The topographic scenarios are created from the variation of two free parameters, (1) the elevation of the Western–Central Alps and (2) the elevation of the Eastern Alps. The results indicate Δδ18Op values (i.e., the difference between δ18Op values at the low- and high-elevation sites) of up to −8 ‰ along the strike of the Alps for the diachronous uplift scenarios, primarily due to changes in orographic precipitation and adiabatic lapse rate driven localized changes in near-surface variables. These simulated magnitudes of Δδ18Op values suggest that the expected isotopic signal would be significant enough to be preserved and measured in geologic archives. Moreover, the simulated slight δ18Op differences of 1 ‰–2 ‰ across the low-elevation sites support the use of the δ–δ paleoaltimetry approach and highlight the importance of sampling far-field low-elevation sites to differentiate between the different surface uplift scenarios. The elevation-dependent rate of change in δ18Op (“isotopic lapse rate”) varies depending on the topographic configuration and the extent of the surface uplift. Most of the changes are significant (e.g., −1.04 ‰ km−1 change with slope error of ±0.09 ‰ km−1), while others were within the range of the statistical uncertainties (e.g., −0.15 ‰ km−1 change with slope error of ±0.13 ‰ km−1). The results also highlight the plausible changes in atmospheric circulation patterns and associated changes in moisture transport pathways in response to changes in the topography of the Alps. These large-scale atmospheric dynamics changes can complicate the underlying assumption of stable isotope paleoaltimetry and therefore require integration with paleoclimate modeling to ensure accurate reconstruction of the paleoelevation of the Alps.","PeriodicalId":48931,"journal":{"name":"Earth System Dynamics","volume":"504 2","pages":""},"PeriodicalIF":7.3,"publicationDate":"2023-11-17","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"139264031","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":2,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}