scholarly journals Causal networks of biosphere–atmosphere interactions

Author(s):  
Christopher Krich ◽  
Jakob Runge ◽  
Diego G. Miralles ◽  
Mirco Migliavacca ◽  
Oscar Perez-Priego ◽  
...  

Abstract. Local meteorological conditions and biospheric activity are tightly coupled. Understanding these links is an essential prerequisite for predicting the Earth system under climate change conditions. However, many empirical studies on the interaction between the biosphere and the atmosphere are based on correlative approaches that are not able to deduce causal paths, and only very few studies apply causal discovery methods. Here, we use a recently proposed causal graph discovery algorithm, which aims to reconstruct the causal dependency structure underlying a set of time series. We explore the potential of this method to infer temporal dependencies in biosphere-atmosphere interactions. Specifically we address the following questions: How do periodicity and heteroscedasticity influence causal detection rates, i.e. the detection of existing and non-existing links? How consistent are results for noise-contaminated data? Do results exhibit an increased information content that justifies the use of this causal-inference method? We explore the first question using artificial time series with well known dependencies that mimic real-world biosphere-atmosphere interactions. The two remaining questions are addressed jointly in two case studies utilizing observational data. Firstly, we analyse three replicated eddy covariance datasets from a Mediterranean ecosystem at half hourly time resolution allowing us to understand the impact of measurement uncertainties. Secondly, we analyse global NDVI time series (GIMMS 3g) along with gridded climate data to study large-scale climatic drivers of vegetation greenness. Overall, the results confirm the capacity of the causal discovery method to extract time-lagged linear dependencies under realistic settings. The violation of the method's assumptions increases the likelihood to detect false links. Nevertheless, we consistently identify interaction patterns in observational data. Our findings suggest that estimating a directed biosphere-atmosphere network at the ecosystem level can offer novel possibilities to unravel complex multi-directional interactions. Other than classical correlative approaches, our findings are constrained to a few meaningful set of relations which can be powerful insights for the evaluation of terrestrial ecosystem models.

2020 ◽  
Vol 17 (4) ◽  
pp. 1033-1061
Author(s):  
Christopher Krich ◽  
Jakob Runge ◽  
Diego G. Miralles ◽  
Mirco Migliavacca ◽  
Oscar Perez-Priego ◽  
...  

Abstract. The dynamics of biochemical processes in terrestrial ecosystems are tightly coupled to local meteorological conditions. Understanding these interactions is an essential prerequisite for predicting, e.g. the response of the terrestrial carbon cycle to climate change. However, many empirical studies in this field rely on correlative approaches and only very few studies apply causal discovery methods. Here we explore the potential for a recently proposed causal graph discovery algorithm to reconstruct the causal dependency structure underlying biosphere–atmosphere interactions. Using artificial time series with known dependencies that mimic real-world biosphere–atmosphere interactions we address the influence of non-stationarities, i.e. periodicity and heteroscedasticity, on the estimation of causal networks. We then investigate the interpretability of the method in two case studies. Firstly, we analyse three replicated eddy covariance datasets from a Mediterranean ecosystem. Secondly, we explore global Normalised Difference Vegetation Index time series (GIMMS 3g), along with gridded climate data to study large-scale climatic drivers of vegetation greenness. We compare the retrieved causal graphs to simple cross-correlation-based approaches to test whether causal graphs are considerably more informative. Overall, the results confirm the capacity of the causal discovery method to extract time-lagged linear dependencies under realistic settings. For example, we find a complete decoupling of the net ecosystem exchange from meteorological variability during summer in the Mediterranean ecosystem. However, cautious interpretations are needed, as the violation of the method's assumptions due to non-stationarities increases the likelihood to detect false links. Overall, estimating directed biosphere–atmosphere networks helps unravel complex multidirectional process interactions. Other than classical correlative approaches, our findings are constrained to a few meaningful sets of relations, which can be powerful insights for the evaluation of terrestrial ecosystem models.


2013 ◽  
Vol 10 (6) ◽  
pp. 7469-7516 ◽  
Author(s):  
M. T. Pham ◽  
W. J. Vanhaute ◽  
S. Vandenberghe ◽  
B. De Baets ◽  
N. E. C. Verhoest

Abstract. Of all natural disasters, the economic and environmental consequences of droughts are among the highest because of their longevity and widespread spatial extent. Because of their extreme behaviour, studying droughts generally requires long time series of historical climate data. Rainfall is a very important variable for calculating drought statistics, for quantifying historical droughts or for assessing the impact on other hydrological (e.g. water stage in rivers) or agricultural (e.g. irrigation requirements) variables. Unfortunately, time series of historical observations are often too short for such assessments. To circumvent this, one may rely on the synthetic rainfall time series from stochastic point process rainfall models, such as Bartlett–Lewis models. The present study investigates whether drought statistics are preserved when simulating rainfall with Bartlett–Lewis models. Therefore, a 105 yr 10 min rainfall time series obtained at Uccle, Belgium is used as test case. First, drought events were identified on the basis of the Effective Drought Index (EDI), and each event was characterized by two variables, i.e. drought duration (D) and drought severity (S). As both parameters are interdependent, a multivariate distribution function, which makes use of a copula, was fitted. Based on the copula, four types of drought return periods are calculated for observed as well as simulated droughts and are used to evaluate the ability of the rainfall models to simulate drought events with the appropriate characteristics. Overall, all Bartlett–Lewis type of models studied fail in preserving extreme drought statistics, which is attributed to the model structure and to the model stationarity caused by maintaining the same parameter set during the whole simulation period.


2020 ◽  
Vol 12 (19) ◽  
pp. 3207
Author(s):  
Ioannis Papoutsis ◽  
Charalampos Kontoes ◽  
Stavroula Alatza ◽  
Alexis Apostolakis ◽  
Constantinos Loupasakis

Advances in synthetic aperture radar (SAR) interferometry have enabled the seamless monitoring of the Earth’s crust deformation. The dense archive of the Sentinel-1 Copernicus mission provides unprecedented spatial and temporal coverage; however, time-series analysis of such big data volumes requires high computational efficiency. We present a parallelized-PSI (P-PSI), a novel, parallelized, and end-to-end processing chain for the fully automated assessment of line-of-sight ground velocities through persistent scatterer interferometry (PSI), tailored to scale to the vast multitemporal archive of Sentinel-1 data. P-PSI is designed to transparently access different and complementary Sentinel-1 repositories, and download the appropriate datasets for PSI. To make it efficient for large-scale applications, we re-engineered and parallelized interferogram creation and multitemporal interferometric processing, and introduced distributed implementations to best use computing cores and provide resourceful storage management. We propose a new algorithm to further enhance the processing efficiency, which establishes a non-uniform patch grid considering land use, based on the expected number of persistent scatterers. P-PSI achieves an overall speed-up by a factor of five for a full Sentinel-1 frame for processing in a 20-core server. The processing chain is tested on a large-scale project to calculate and monitor deformation patterns over the entire extent of the Greek territory—our own Interferometric SAR (InSAR) Greece project. Time-series InSAR analysis was performed on volumes of about 12 TB input data corresponding to more than 760 Single Look Complex Sentinel-1A and B images mostly covering mainland Greece in the period of 2015–2019. InSAR Greece provides detailed ground motion information on more than 12 million distinct locations, providing completely new insights into the impact of geophysical and anthropogenic activities at this geographic scale. This new information is critical to enhancing our understanding of the underlying mechanisms, providing valuable input into risk assessment models. We showcase this through the identification of various characteristic geohazard locations in Greece and discuss their criticality. The selected geohazard locations, among a thousand, cover a wide range of catastrophic events including landslides, land subsidence, and structural failures of various scales, ranging from a few hundredths of square meters up to the basin scale. The study enriches the large catalog of geophysical related phenomena maintained by the GeObservatory portal of the Center of Earth Observation Research and Satellite Remote Sensing BEYOND of the National Observatory of Athens for the opening of new knowledge to the wider scientific community.


2019 ◽  
Vol 12 (1) ◽  
pp. 45 ◽  
Author(s):  
Jiaqi Chen ◽  
Jiming Lv ◽  
Ning Li ◽  
Qingwei Wang ◽  
Jian Wang

There are a large number of lakes with beaded distribution in the semi-arid areas of the Inner Mongolian Plateau, and some of them have degraded or even disappeared during the past three decades. We studied the reasons of the disappearance of these lakes by determining the way of replenishment of these lakes and the impact of the natural-social environment of the basin, with the aim of saving these gradually disappearing lakes. Based on remote sensing image and hydrological analysis, this paper studied the recharge of Daihai Lake and Huangqihai Lake. The deep learning method was used to establish the time-series of lake evolution. The same method was combined with the innovative woodland and farmland extraction method to set up the time-series of ground classification composition in the basins. Using relevant survey data, combined with soil water infiltration test, water chemical, and isotopic signature analysis of various water bodies, we found that the Daihai Lake area is the largest in dry season and the smallest in rainy season and the other lake is not satisfied with this phenomenon. In addition, we calculated the specific recharge and consumption of the study basin. These experiments indicated that the exogenous groundwater is recharged directly through the faults at the bottom of Daihai Lake, while the exogenous groundwater is recharged in Huangqihai Lake through rivers indirectly. Large-scale exploitation of groundwater for agricultural irrigation and industrial production is the main cause of lake degradation. Reducing the extraction of groundwater for agricultural irrigation is an important measure to restore lake ecology.


2013 ◽  
Vol 17 (12) ◽  
pp. 5167-5183 ◽  
Author(s):  
M. T. Pham ◽  
W. J. Vanhaute ◽  
S. Vandenberghe ◽  
B. De Baets ◽  
N. E. C. Verhoest

Abstract. Of all natural disasters, the economic and environmental consequences of droughts are among the highest because of their longevity and widespread spatial extent. Because of their extreme behaviour, studying droughts generally requires long time series of historical climate data. Rainfall is a very important variable for calculating drought statistics, for quantifying historical droughts or for assessing the impact on other hydrological (e.g. water stage in rivers) or agricultural (e.g. irrigation requirements) variables. Unfortunately, time series of historical observations are often too short for such assessments. To circumvent this, one may rely on the synthetic rainfall time series from stochastic point process rainfall models, such as Bartlett–Lewis models. The present study investigates whether drought statistics are preserved when simulating rainfall with Bartlett–Lewis models. Therefore, a 105 yr 10 min rainfall time series obtained at Uccle, Belgium is used as a test case. First, drought events were identified on the basis of the Effective Drought Index (EDI), and each event was characterized by two variables, i.e. drought duration (D) and drought severity (S). As both parameters are interdependent, a multivariate distribution function, which makes use of a copula, was fitted. Based on the copula, four types of drought return periods are calculated for observed as well as simulated droughts and are used to evaluate the ability of the rainfall models to simulate drought events with the appropriate characteristics. Overall, all Bartlett–Lewis model types studied fail to preserve extreme drought statistics, which is attributed to the model structure and to the model stationarity caused by maintaining the same parameter set during the whole simulation period.


2020 ◽  
Vol 10 (1) ◽  
Author(s):  
Josué M. Polanco-Martínez ◽  
Javier Fernández-Macho ◽  
Martín Medina-Elizalde

AbstractThe wavelet local multiple correlation (WLMC) is introduced for the first time in the study of climate dynamics inferred from multivariate climate time series. To exemplify the use of WLMC with real climate data, we analyse Last Millennium (LM) relationships among several large-scale reconstructed climate variables characterizing North Atlantic: i.e. sea surface temperatures (SST) from the tropical cyclone main developmental region (MDR), the El Niño-Southern Oscillation (ENSO), the North Atlantic Multidecadal Oscillation (AMO), and tropical cyclone counts (TC). We examine the former three large-scale variables because they are known to influence North Atlantic tropical cyclone activity and because their underlying drivers are still under investigation. WLMC results obtained for these multivariate climate time series suggest that: (1) MDRSST and AMO show the highest correlation with each other and with respect to the TC record over the last millennium, and: (2) MDRSST is the dominant climate variable that explains TC temporal variability. WLMC results confirm that this method is able to capture the most fundamental information contained in multivariate climate time series and is suitable to investigate correlation among climate time series in a multivariate context.


2012 ◽  
Vol 9 (1) ◽  
pp. 175-214
Author(s):  
D. González-Zeas ◽  
L. Garrote ◽  
A. Iglesias ◽  
A. Sordo-Ward

Abstract. An important aspect to assess the impact of climate change on water availability is to have monthly time series representative of the current situation. In this context, a simple methodology is presented for application in large-scale studies in regions where a properly calibrated hydrologic model is not available, using the output variables simulated by regional climate models (RCMs) of the European project PRUDENCE under current climate conditions (period 1961–1990). The methodology compares different interpolation methods and alternatives to generate annual times series that minimize the bias with respect to observed values. The objective is to identify the best alternative to obtain bias-corrected, monthly runoff time series from the output of RCM simulations. This study uses information from 338 basins in Spain that cover the entire mainland territory and whose observed values of naturalised runoff have been estimated by the distributed hydrological model SIMPA. Four interpolation methods for downscaling runoff to the basin scale from 10 RCMs are compared with emphasis on the ability of each method to reproduce the observed behavior of this variable. The alternatives consider the use of the direct runoff of the RCMs and the mean annual runoff calculated using five functional forms of the aridity index, defined as the ratio between potential evaporation and precipitation. In addition, the comparison with respect to the global runoff reference of the UNH/GRDC dataset is evaluated, as a contrast of the "best estimator" of current runoff on a large scale. Results show that the bias is minimised using the direct original interpolation method and the best alternative for bias correction of the monthly direct runoff time series of RCMs is the UNH/GRDC dataset, although the formula proposed by Schreiber also gives good results.


MAUSAM ◽  
2021 ◽  
Vol 52 (1) ◽  
pp. 67-82
Author(s):  
J. R. KULKARNI ◽  
M. MUJUMDAR ◽  
S. P. GHARGE ◽  
V. SATYAN ◽  
G. B. PANT

Earlier investigations into the epochal behavior of fluctuations in All India Summer Monsoon Rainfall (AISMR) have indicated the existence of a Low Frequency Mode (LFM) in the 60-70 years range. One of the probable sources of this variability may be due to changes in solar irradiance. To investigate this, time series of 128-year solar irradiance data from 1871-1998 has been examined. The Wavelet Transform (WT) method is applied to extract the LFM from these time series, which show a very good correspondence. A case study has been carried out to test the sensitivity of AISMR to solar irradiance. The General Circulation Model (GCM) of the Center of Ocean-Land-Atmosphere (COLA) has been integrated in the control run (using the climatological value of solar constant i.e., 1365 Wm-2) and in the enhanced solar constant condition (enhanced by 10 Wm-2) for summer monsoon season of 1986. The study shows that the large scale atmospheric circulation over the Indian region, in the enhanced solar constant scenario is favorable to good monsoon activity. A conceptual model for the impact of solar irradiance on the AISMR at LFM is also suggested.


2019 ◽  
Vol 11 (14) ◽  
pp. 1665 ◽  
Author(s):  
Tianle He ◽  
Chuanjie Xie ◽  
Qingsheng Liu ◽  
Shiying Guan ◽  
Gaohuan Liu

Machine learning comprises a group of powerful state-of-the-art techniques for land cover classification and cropland identification. In this paper, we proposed and evaluated two models based on random forest (RF) and attention-based long short-term memory (A-LSTM) networks that can learn directly from the raw surface reflectance of remote sensing (RS) images for large-scale winter wheat identification in Huanghuaihai Region (North-Central China). We used a time series of Moderate Resolution Imaging Spectroradiometer (MODIS) images over one growing season and the corresponding winter wheat distribution map for the experiments. Each training sample was derived from the raw surface reflectance of MODIS time-series images. Both models achieved state-of-the-art performance in identifying winter wheat, and the F1 scores of RF and A-LSTM were 0.72 and 0.71, respectively. We also analyzed the impact of the pixel-mixing effect. Training with pure-mixed-pixel samples (the training set consists of pure and mixed cells and thus retains the original distribution of data) was more precise than training with only pure-pixel samples (the entire pixel area belongs to one class). We also analyzed the variable importance along the temporal series, and the data acquired in March or April contributed more than the data acquired at other times. Both models could predict winter wheat coverage in past years or in other regions with similar winter wheat growing seasons. The experiments in this paper showed the effectiveness and significance of our methods.


Sign in / Sign up

Export Citation Format

Share Document