Climate-driven stock accounts of inland Green-Blue Ecosystems: Coupling the carbon and water cycles via Reverse Engineering and Geodetectors

Climate-driven stock accounts of inland Green-Blue Ecosystems: Coupling the carbon and water cycles via Reverse Engineering and Geodetectors | Research Square window.SnipcartSettings = { analytics: { enabled: false } }; (function() { var accessVector = localStorage.getItem('access_vector') || ''; window.dataLayer = window.dataLayer || []; if (accessVector) { window.dataLayer.push({ user: { profile: { profileInfo: { snid: accessVector } } } }); } })(); (function(w,d,s,l,i){w[l]=w[l]||[];w[l].push({'gtm.start':new Date().getTime(),event:'gtm.js'});var f=d.getElementsByTagName(s)[0],j=d.createElement(s),dl=l!='dataLayer'?'&l='+l:'';j.async=true;j.src='https://www.googletagmanager.com/gtm.js?id='+i+dl;f.parentNode.insertBefore(j,f);})(window,document,'script','dataLayer','GTM-K279D39R'); Browse Preprints In Review Journals COVID-19 Preprints AJE Video Bytes Research Tools Research Promotion AJE Professional Editing AJE Rubriq About Preprint Platform In Review Editorial Policies Our Team Advisory Board Help Center Sign In Submit a Preprint Cite Share Download PDF Research Article Climate-driven stock accounts of inland Green-Blue Ecosystems: Coupling the carbon and water cycles via Reverse Engineering and Geodetectors Bruna Almeida, Luís Monteiro, Pierre Scemama, Pedro Cabral This is a preprint; it has not been peer reviewed by a journal. https://doi.org/ 10.21203/rs.3.rs-4681296/v1 This work is licensed under a CC BY 4.0 License Status: Posted Version 1 posted You are reading this latest preprint version Abstract The critical role of inland Green-Blue Ecosystems (GBE) in delivering goods and services underscores the need to understand their relationships. This research investigates the impact of climate variables on GBE stock accounts by estimating Tree Cover Density (TCD) and Water & Wetness Probability Index (WWPI). Using supervised machine learning and factor analysis, we measured GBE extent and condition. Key predictors for Blue Ecosystems were topographic variables, while vegetation indices were crucial for Green Ecosystems. In 2018, 33% of the inland area was covered by forests, freshwater, and wetlands. Key climate-driven factors for forests included precipitation (0.65), aridity index (0.54), and evapotranspiration (0.44). For freshwater and wetlands, precipitation (0.69), aridity index (0.55), and elevation (0.42) were significant. This research enhances our understanding of how climate impacts GBE, influencing biomass density and water availability. It bridges socio-environmental science with engineering by integrating advanced modelling techniques, promoting ecosystem resilience and sustainability. Environmental Modelling Machine Learning Google Earth Engine Data fusion Copernicus Land Monitoring Services Figures Figure 1 Figure 2 Figure 3 Figure 4 Figure 5 Figure 6 Figure 7 Figure 8 Highlights Classification tasks measured GBE extent, regression estimated GBE condition Blue Ecosystems classification top three feature importance: DSM, NDWI and SLOPE Climate-driven stock accounts of GBE were analyzed through geodetectors The most important drivers for Green Ecosystems were PREC, AI, and PET The top three drivers for Blue Ecosystems were PREC, AI, and DEM 1. Introduction Understanding the terrestrial biosphere’s functioning requires an assessment of the magnitude and distribution of Green-Blue Ecosystems (GBE) (Crowther et al., 2015 ). Green Ecosystems (GE) support a considerable amount of global biodiversity, playing an important role in biogeochemical cycles, and providing several Ecosystem Services (ES) such as water quantity and quality regulation, food production, raw materials provisioning, genetic resources, and carbon storage and sequestration (Mulatu et al., 2017 ). Forest density influences soil and water retention rates, flood regulation, as well as competitive dynamics and habitat suitability for a variety of plant and animal species (Rijal et al., 2021 ). Blue Ecosystems (BE) provide the backbone of ecosystems, enabling agriculture, and influencing urban development decisions (Newbold et al., 2015 ). As forests are intricately connected to rainfall and water availability, changes in vegetation cover, as well as poorly planned and managed forests, have significant impacts on the carbon and water cycles, leading to pressures on the environment (Liquete et al., 2011 ). Extreme weather events and soil and water contamination are perceived as the top risk, directly impacting natural capital and biodiversity, and consequently decreasing human well-being (Lundqvist & Unver, 2018 ). Knowledge of the interactions between carbon and water cycles is vital for predicting ecosystem responses to climate change, as well as the influence on natural capital (Seddon et al., 2021 ). Despite significant advances, the coupling mechanisms of these cycles under socio-economic and climatic pressures remain poorly understood (Keith et al., 2021 ). Ecosystem Accounting (EA) encompass a statistical framework for quantifying and valuing the interactions between nature, society, and the economy, providing physical and monetary measurements of ecosystem assets (Bateman & Mace, 2020 ). In March 2021, the United Nations Statistical Commission adopted the System of Environmental-Economic Accounting Ecosystem-Accounting (SEEA-EA) as a baseline to inform policy development by countries. It represents a huge milestone in the history of accounting as “it moves beyond GDP and takes better account of biodiversity and ecosystems in national economic planning” (Edens et al., 2022 ). Thus, it is susceptible to shaping the future of the national statistic apparatus toward the inclusion of natural capital (Edens et al., 2022 ; Lange et al., 2022 ). Among the challenges still to be overcome is the substantial wide variety and amount of data required to produce comprehensive accounts (Bordt, 2018 ; Hein et al., 2020 ). Besides, extending locally generated ES models to other locations or scales has constraints due to the necessity for parameterization, calibration, and validation, which is sometimes hampered by a lack of ground truth data (Cord et al., 2017 ; Gosal et al., 2022 ). The opportunities associated with technological innovations include the potential for gathering, processing, analyzing and visualizing different types of data, and integrating socio-economic and environmental information (Fleming et al., 2022 ). Advancements in geospatial techniques, such as satellite imagery, geospatial mapping, drones, and sensors, are viewed as highly promising to improve decision-making on both individual and cross-sectoral levels facilitating the deployment of EA (Farrell et al., 2021 ). Besides, the increased computational power and the transition towards open-access data and open-source technologies are significantly contributing to enhancing natural capital accounting (del Río-Mena et al., 2023 ). Satellite Earth Observation (SEO) data and technologies, coupled with Geographic Information Systems (GIS), provide advantages such as synoptic and repeated coverage, historical spatial data analysis, and cost-effectiveness (Ramirez-Reyes et al., 2019 ). Therefore, where data is unavailable, imputation data techniques through statistics and Machine Learning (ML) are employed to fill in missing or incomplete data (Fleming et al., 2022 ). Additionally, data fusion is commonly applied, integrating information from diverse sources such as monitoring, statistics, modelling, or interviews (Braun et al., 2018 ). These approaches enable more frequent, and consistent assessment of ecosystem dynamics, regardless spatial scale of analysis (Hossain & Hashim, 2019 ). In developing models for quantifying natural stocks, accurate and timely data on the extent and status of ecosystems are critical and challenging to gather particularly due to the dynamic and intricate interactions between water/groundwater and vegetation (Ellison et al., 2017 ). Significant efforts have been made to unveil the potential of SEO in ES assessments (Willcock et al., 2018 ), such as the European Copernicus programme powered by the European Spatial Agency ( 2023 ). The platform yields free and open access to SEO imagery and datasets, providing insights into the status of natural resources not only for Europe but also on a global scale (Almeida & Cabral, 2023 ). Among the various products, it provides the Water & Wetness Probability Index (WWPI) and the Tree Cover Density (TCD) to assist in building policy-relevant ecological accounts. The WWPI is a biophysical indicator for freshwater and wetland ecosystems (Ludwig et al., 2019 ) utilized as a spatial proxy to support environmental mitigation, wetland protection, erosion control, flood monitoring, and streamflow regulation (Vargas et al., 2019 ). The dataset provides information on the status of BE reflecting their overall quality and conditions (Copernicus Programme, 2023 ). The TCD is an indicator describing forest ecosystem conditions through a continuous spectrum of crown cover information, population numbers, densities, and wood stocks (Crowther et al., 2015 ). It is acknowledged as an effective method for environmental analyses, forest management, supporting decision-making, and tracking changes in tree cover losses and gains (Chen et al., 2020 ). Therefore, these products have some limitations due to: a) the availability and quality of SEO data, b) the availability and quality of in-situ data, and c) the influence of climatic and topographic conditions (European Environment Agency, 2023 ). Satellite mission launch dates restrict SEO data availability, and the quality of imagery is influenced by the atmospheric effects on the light reflected off the Earth's surface and captured by the sensor. Another challenge is the lack of ground-truth data, which is expensive to acquire (Mairota et al., 2015 ). Some studies have demonstrated improvements with the integration of climate and topographic variables when assessing GBE. Crowther et al. ( 2015 ) demonstrated the strong associations between climatic variables and forest density in their research. Ludwig et al. ( 2019 ) combined multi-temporal optical imagery, topographic data, and spectral indices to build an automated wetland mapping detection, based on the WWPI. Madrigal-González et al. ( 2023 ) used climatic and topographic characteristics to investigate the relationship between tree density and water availability. Han et al (2021) estimated carbon stock in forests by considering temperature, precipitation, elevation, and soil type as factors. Alqadhi et al (2022) demonstrated the impact of landscape and topographic characteristics on ES capacity provisioning, considering: elevation, slope, aspect, drainage density, and precipitation. However, none of these studies used Reverse Engineering (RE) to reproduce TCD and WWPI taking into consideration climatic and topographic factors to estimate inland GBE at the national level. RE refers to the process of extracting knowledge or design information from products and subsequently reproducing or recreating them based on the acquired information (Nieves-Chinchilla et al., 2018 ). Through RE, it is possible to learn about the products, search for inconsistencies, or overall vulnerabilities, and assess whether there is a more efficient way to improve them (Wood, 2009 ). In this study, the WWPI and TCD will be modelled within Portugal's mainland in ArcGIS Pro (ESRI, 2023 ) using the Random Trees Regression algorithm, and spatial associations will be analyzed through geographical detectors utilizing the geodetectors package (version: 1.0–4) in R software (Wang et al., 2016 ). EA baselines (Kienast et al., 2009 ) will be followed to assess the climate-driven stock accounts of forests, water, and wetlands at the national level. The aim is to develop more integrated and updated versions of TCD and WWPI contributing to the broader scientific community by exploring the potential of SEO data and technologies combined with global climatic and topographic open-access data. Furthermore, this research will shed light on the interplay between aquatic and terrestrial ecosystems and foster innovative solutions for assessing the status and conditions of inland GBE. Advancing our understanding of the development of Copernicus products will deepen our knowledge of the replicability and applicability of such datasets in EA. 2. Materials and Methods 2.1. Study Area Portugal is located on the Iberian Peninsula in southwestern Europe, comprising a continental portion and two autonomous regions - the Azores and Madeira archipelagos. Portugal's mainland spans an approximate area of 90000 square km located between 36°57' − 42°9'N latitude and 6°12' − 9°30'W longitude in a climatic Mediterranean region (Beck et al., 2018 ) (Fig. 1 ). The environment changes importantly from north to south, with complex and diverse landscapes (Almeida & Cabral, 2021 ). Spatial and temporal variability is characterized by latitude, orography, and the effect of the Atlantic Ocean, which affects the water cycle primarily through oscillations in precipitation and temperature (Belo-Pereira et al., 2011 ). The territory is prone to climate variability, including droughts and desertification, like other southern European regions, notably in the southern sector, which receives the least precipitation and has the greatest temperatures (Oliveira et al., 2007 ). The climatic conditions, influencing vegetation patterns, contribute to a diverse array of Atlantic, European, and Mediterranean plant species within the country (Santos et al., 2023 ). Forest ecosystems predominantly consist of broad-leaved and coniferous forests, including Cork oak, Holm oak, Eucalyptus, Maritime pine, and Stone pine (Costa et al., 2022 ). This ecological diversity is evident across the country, with the northern regions characterized by flourishing broadleaf and coniferous forests like Eucalyptus, Maritime pine, and Stone pine. In contrast, the southern landscape exhibits arid-adapted Mediterranean species such as cork oak, holm oak, and olive trees, representing distinct climatic zones in the country (Fonseca et al., 2019 ). 2.2. Data Table 1 presents the data used in this study, describing their source, spatial resolution, and the role played by each dataset in the classification and regression tasks. Table 1 Description of datasets utilized in the research. Legend: C – Classification; R – Regression; GEE – Google Earth Engine. Data Source Pixel resolution Task Role Satellite imagery 2018 - Single bands European Space Agency (ESA) – GEE 10m C Feature Satellite imagery 2018 - Spectral indices European Space Agency (ESA) – GEE 10m C Feature ‍Global Digital Surface Model (DSM) (m) Earth Observation Research Center Japan Aerospace Exploration Agency (JAXA EORC) - GEE 30m C Feature Land Use and Land Cover Map (COS2018) 2018 (m) Directorate General for Territory - DGT ( https://snig.dgterritorio.gov.pt/ ) 10m C Target Digital Elevation Model (DEM) (m) & Slope (SLOPE) (%) NASA Shuttle Radar Topography Mission ( https://www.earthdata.nasa.gov/ ) 30m R Predictor Average Annual Precipitation 2018 (PREC) (mm) WorldClim ( https://worldclim.org/ ) 1km R Predictor Global Potential Evapotranspiration 2018 (PET) (mm) Consortium of Spatial Information - CSI ( http://www.cgiar-csi.org ) 1km R Predictor Global Aridity Index 2018 (AI) Consortium of Spatial Information -CSI ( http://www.cgiar-csi.org ) 1km R Predictor Water & Wetneess Probability Index (WWPI %) 2018 Copernicus Land Monitoring Services - CLMS ( https://land.copernicus.eu/en/products ) 10m R Response variable Tree Cover Density (TCD %) 2018 Copernicus Land Monitoring Services - CLMS ( https://land.copernicus.eu/en/products ) 10m R Response variable In the GEE platform, cloud-free composites of Sentinel-2 (S2) imagery encompassing the entire mainland of Portugal for 2018 were selected. All available spectral bands from S2 imagery, satellite-derived indices, and topographic features, were employed as classification predictors. These included the Normalized Difference Vegetation Index (NDVI)(Rouse et al., 1973 ) and the Enhanced Vegetation Index (EVI)(Huete et al., 1997 ) used to highlight and distinguish vegetation of different types and stages throughout their seasonal cycles; the Soil Adjusted Vegetation Index (SAVI)(Huete, 1988 ) to mitigate the effects of soil reflectance and improve vegetation response, and the NDWI (Normalized Difference Water Index) used to differentiate between water and soil, based on the optical characteristics of water ​(Ji et al., 2015 ). The S2-derived indices were calculated using the band arithmetic functions: (1), (2), (3), and (4): $$\:NDVI=\frac{B8-B4}{B8+B4}$$ 1 $$\:EVI=2.5*\frac{B8-B4}{B8+6*B4-7.5*B2+1}$$ 2 $$\:SAVI=\frac{B8-B4}{(B8+B4+L)*\left(1.0*L\right)}$$ 3 $$\:NDWI=\:\frac{B3-B8}{B3+B8}$$ 4 In the formulas, B refers to the spectral band, and L is the soil brightness correction that ranges between [0, 1]. In the SAVI calculation, L takes the value of 0.428, as suggested by Huete ( 1988 ). The spectral indices (1), (2) and (3) were used to classify forests, and (4) was employed in the freshwater and wetlands classification. Topographic variables, such as the Digital Surface Model (DSM) and Slope, were provided by the Japan Aerospace Exploration Agency (JAXA, 2023 ). These datasets were added with the bands and spectral indices as an 18-band data stack to perform the classification tasks. The classification training datasets included reference vectorial data representing GBE (forest, freshwater, and wetlands) were extracted from the Portuguese land use and land cover map for the year 2018 (COS2018), a national product developed by the Directorate General for Territory (DGT, 2018 ). Climatic variables, such as Precipitation (PREC), Potential Evapotranspiration (PET), Aridity Index (AI), and the Digital Elevation Model (DEM) and Slope (SLOPE) were further utilized as explanatory variables/factors in the regression tasks. Annual PREC was obtained from the WorldClim database (Fick & Hijmans, 2017 ); Global PET and the Global AI were obtained from the Consortium of Spatial Information, Global-AI and Global-PET Database (Zomer et al., 2022 ). Global PET and AI were derived from the WorldClim precipitation model and have the exact spatial resolution and temporal scales. The response variables (WWPI and TCD) are from the High-Resolution Layers portfolio of the Copernicus. The WWPI is a raster that depicts the presence of water and wet areas as an index ranging from 0 to 100% (Copernicus Programme, 2023 ). The WWPI workflow includes satellite data fusion, supervised ML classification and an unsupervised thresholding technique for water/wetness detection. Data fusion is based on optical data from S2 and Landsat 5/7/8 and radar data from Sentinel-1 multi-temporal imagery from 2012–2018. The classification task was performed with the RF algorithm, and post-processing included quality assurance, control procedures, visual improvements and error correction. The TCD is a vertical projection of tree crowns onto a horizontal Earth’s surface providing pixel-level information on tree cover density (%). In the TCD development, Sentinel-2A + B time series (Level-2A data) were employed to build a supervised classification model using an RF classifier to create a binary Tree Cover Mask (TCM) (non-tree covered areas/tree cover). Following that, a multiple linear regression was carried out inside the bounds of the TCM to estimate TCD values. 2.3. Methods Forward Engineering goes from concept to product, while RE does the opposite, applying techniques to create a blueprint of an existing product (Vallero, 2014 ). RE is frequently used to replicate existing products, test, or improve functionalities, and introduce new data (Nieves-Chinchilla et al., 2018 ). It is defined as the process of examining an object to determine its components and their interrelationships for further redesign and recreation (Wood, 2009 ). There are a few steps to follow when applying RE, such as defining the purpose, reviewing product specifications, identifying limitations and vulnerabilities, disassembling, analysing, redesigning, and reporting (Chikofsky & Cross II, 2002). The research development starts applying the RE approach to understand how TCD and WWPI were created, how they work, and how they could be improved, followed by employing the EA framework to quantify the physical measurements of inland GBE, through estimations of ecosystems extent (Fig. 2 a) and condition (Fig. 2 b). Ecosystem extent quantify the area per ecosystem type inside an accounting area. In contrast, ecosystem condition record the status of an ecosystem asset through indicators that reflect its condition (United Nations et al., 2021 ). The extent was estimated through classification tasks mapping GE and BE, while the condition of forest, freshwater and wetlands was predicted by regression tasks. Finally, there was the calculation of factor importance and interaction, and statistic tests between factors and within ecosystems, employing geodetectors. 2.3.1. Classification Task Two distinct classification tasks were conducted: one dedicated to creating the Forest Cover Mask (FCM) and another to creating the Freshwater and Wetlands Mask (FWM). The training data were extracted from the COS18 dataset, and used as samples for a Random Forest (RF) classification (Breiman, 2001 ), in Google Earth Engine (GEE) (Google, 2023 ). The scripts are provided in Tables S1 and S2 in the Supplementary Materials. Random point sampling was employed to provide training and validation data for the tasks. The RF ML algorithm evaluated the training features in terms of satellite-derived data, DSM, and slope. Additionally, to optimize the model and improve its accuracy, an analysis of feature importance and hyperparameter tuning was conducted for each classification. These analyses enabled the identification of the most relevant features and refining the RF algorithm, enhancing classification accuracy for each designated class. The classification result was then subject to an accuracy assessment procedure using the validation points by calculating different metrics, such as overall accuracy (OA), confusion matrix, producers’ accuracy (PA), users’ accuracy (UA), kappa coefficient, and F1 score. Lastly, post-processing techniques were applied to eliminate isolated pixels and noise. 2.3.2. Regression Task Predictive spatial ML regression was used to model the relationship between explanatory variables and response variables (WWPI and TCD) using ArcGIS Pro 2.9.0 software. The regression tasks were implemented using the tool Train Random Trees Regression Model . The ML spatial regression models for TCD and WWPI were built considering five explanatory variables (PREC, PET, AI, DEM, and SLOPE). The Percent Samples for Testing option was set to 20%, which means that one-fifth of the training sample, referred to as test location points, was used to quantify the error for interpolation in space. Three types of errors were measured: errors on training points, errors on test points, and errors on test location points. The maximum number of trees , maximum tree depth , and maximum number of samples were respectively set to 50, 30, and 100000. The maximum depth of each tree refers to the number of rules that each tree can generate to decide. It is critical to evaluate regression performance to understand how well the model is fitted and explained by independent variables (Hosmer & Lemeshow, 2000 ). The metrics utilized to detect bias and the amount of variance of the response variables were the coefficient of determination (R-squared) and regression error (Re). The tool outputs a table describing the importance of each predictor used in the model, as well as scatterplots of training data, test data, and test location data, and a regression definition file containing attribute information, statistics, and model performance. 2.3.3. Factor Analysis using Geodetectors Factor analysis was carried out to estimate the drivers of GBE using geographical detectors. Geodetector is a tool for investigating spatial stratified heterogeneity (SSH) through (1) measuring SSH of a variable Y; (2) testing the power of a determinant X of a dependent variable Y based on the consistency of their spatial distributions; and (3) investigating the interaction between two explanatory variables X1 and X2 (Wang et al., 2016 ). The strata of Y are a partition of Y by a categorical explanatory variable X (h(X)) (J. Yang et al., 2022 ). In this study, the strata refers to each ecosystem, hence the q-statistic was calculated for GE and BE. The q-statistic value is between [0, 1], then q = 0 indicates no relationship between Y and X; q = 1 shows that Y is governed by X. The tool is implemented by the geographical detector q -statistic (5): $$\:q=1-\left(\frac{1}{N{\sigma\:}^{2}}\right){\sum\:}_{h=1}^{a}{N}_{h}{\sigma\:}_{h}^{2}$$ 5 where N and \(\:{\sigma\:}^{2}\) represent the number of units and Y variance, respectively; the population Y is made up of a strata (h = 1, 2,..., a ), and \(\:{N}_{h}\) and \(\:{\sigma\:}_{h}^{2}\) represents the number of units and Y variance in stratum h, respectively. The factor analyses were calculated using the R geodetector package (Wang et al., 2016 ). The function returns the q statistic as well as the respective p-value. The ecological detector determines whether there are statistically significant differences between two risk factors X 1 ~ X 2 . The F statistic is utilized in this function to test the differences at a significant threshold of 0.05. The interaction detector determines whether the factors interact with Y. The risk detector computes the average values for each stratum of the explanatory variable (X) and shows if there are differences between strata. 3. Results 3.1. Forest, freshwater and wetlands mapping A binary classification task mapping forest/non-forest areas was performed at the national scale. Also, a multilabel classification task delineating freshwater and wetlands areas was conducted. Both classification maps are required for the calculation of ecosystem extent and are employed as masks in the subsequent step, the regression tasks (Fig. 3 ). The FCM and FWM classifications exhibit high accuracy, strong agreement, and reliability in their respective classifications. Table 2 presents the accuracy evaluation for forest classification (FCM) and freshwater and wetlands classification (FWM). Table 2 Accuracy evaluation. OA: Overall Accuracy; PA: Producer’s Accuracy; UA: User’s Accuracy. LULC OA (%) PA (%) UA (%) Kappa F1 Score GE 96.17 93.68 95.71 0.89 0.93 BE 94.66 95.83 96.67 0.94 0.94 The FWM classification achieved an overall accuracy (OA) of 94.66%, while the FCM classification achieved an even higher accuracy of 96.17%, indicating a high level of correctness in their respective classifications. For the producer’s accuracy (PA), the FCM classification demonstrated a result of 93.68%, while the FWM classification showed a slightly higher result (95.83%), indicating their ability to identify their respective categories accurately. As for user accuracy (UA), the FCM classification exhibited an accuracy of 95.71%, while the FWM classification demonstrated an accuracy of 96.67%, indicating the reliability of both models’ predictions. The forest cover classification exhibited a robust agreement, boasting a Kappa value of 0.89. The classification of freshwater and wetlands demonstrated an even higher level of agreement, with a Kappa value of 0.94. These results underscore both models’ strong consistency and reliability in accurately categorizing vegetation, freshwater, and wetlands. The F1 scores further reinforce the models’ effectiveness in correctly identifying these land cover categories. The FCM classification achieved a commendable F1 score of 0.93, while the FWM classification excelled with an impressive F1 score of 0.94. When looking at the mean percentage weight of feature importance for the water and wetlands classification at the national level (Fig. 4 ), the analysis provided insights into the features crucial for discriminating water bodies and wetland areas accurately. As such, among the 13 attributes assessed, the surface elevation (DSM) exhibited the highest significance (13.65%), followed by NDWI (9.70%) and Slope (9.30%). S2 Bands (8, 11, 12, 8A, 4, 6, 5 and 7) presented importance higher than 5%, while Bands (3 and 2) were less than 5%. Regarding the forest cover classification, a total of 15 attributes were considered as predictors. The feature analysis demonstrates that the spectral indices, SAVI (9.14%), NDVI (8.19%) and EVI (6.89%), showcased notable significance, but also S2 Bands 11 (8.70%), 3 (8.09%) and 12 (7.91%). The Short Wave Infrared (SWIR) bands, B11 and B12, capture variations in vegetation water content and the spongy mesophyll structure within vegetation canopies. Meanwhile, B3 is sensitive to healthy vegetation. The S2 Band 8A was the only feature contributing with less than 5%. The S2 Bands (7, 8, 6, 5, 2 and 4), DSM and Slope presented an importance weight of more than 5% but less than 6.89%. 3.2. FCD and FWPI predictions R-squared and residuals were used to assess the accuracy of FCD and FWPI predictions. The true values were compared to predicted values in terms of proportion of variance (R-squared), and residuals over targets. Figure 5 shows the scatterplot of predicted (FCD) over-reference values (TCD), and results for R-squared at training points (Tr), at test points (Te), and test location points (Tl), as well as the residual plot. The R-squared at Tr, Te and Tl, were respectively, 0.8746, 0.5456, and 0.3208. The R-squared for residuals is 0.4131 and the mean residual value is 0.04. These statistics were calculated using 100000 samples. The scatterplot of predictions (FWPI) over observed values (WWPI) and the residual plot for the FWPI model are shown in Fig. 6 . R-squared values at Tr, Te, and Tl were 0.8275, 0.6153, and 0.6555, respectively. The R-squared residual is 0.4765 and the mean residual value is 0.08. In the predictions of FWPI, the R-squared showcased superior results in Te and Tl. However, FCD performed better than FWPI, as R-squared at Tr is higher, and corresponds to 80% of samples. Table 3 presents the equations for each sample set (Tr, Te, and Tl) used to assess the models’ performance. Table 3 Performance evaluation by model (Tr: Training data at training locations, Te: Test data at training locations, and Tl: Test data at test locations). Model Residuals Tr Te Tl FCD \(\:y=0.32x-15.76\) \(\:y=0.74x+13.68\) \(\:y=0.54x+22.83\) \(\:y=0.32x+24.62\) FWPI \(\:y=0.36x-\:19.82\) \(\:y=0.71x+17.9\) \(\:y=0.55x+24.71\) \(\:y=0.6x+24.21\) The explanatory variables for both regression models were DEM, SLOPE, PREC, PET, and AI. Figure 7 describes the importance of each variable for FCD and FWPI. The predictors with the highest percentage values were more relevant to the model. All variables were shown to be important for both models, except the SLOPE. DEM and PREC were the factors that most contributed to the FWPI model, followed by PET and AI, and lastly, SLOPE. For the FCD model, PET and AI weighted the most, followed by PREC, DEM, and SLOPE. 3.3. Factor analysis using geodetectors The factor detector q-statistic indicates if explanatory variables stratify Y (GE and BE). Table 4 depicts the factor detector results. In the prediction of BE, PREC and AI presented q-values higher than 0.5; DEM and PET higher than 0.25; and SLOPE close to 0, meaning there was no association between FWPI predictions and SLOPE. For GE, PREC and AI show the highest degree of association with values higher than 0.5, followed by PET (0.44), SLOPE (0.13) and DEM (0.11). Table 4 Factor detector results for Blue Ecosystems (BE) and Green Ecosystems (GE) estimation. Y Statistics SLOPE DEM PET AI PREC BE q-statistic 0.06 0.42 0.29 0.55 0.69 BE p-value 1 9.04E-10 9.05E-10 9.35E-10 9.19E-10 GE q-statistic 0.13 0.11 0.44 0.54 0.65 GE p-value 1 9.50E-10 8.47E-10 9.66E-10 9.36E-10 The ecological detector determines if there are statistically significant differences between two covariates. This function utilizes the F-statistic to test for differences, employing a significance threshold of 0.05. In both GE and BE models, all variables were found to be significantly different. The interaction detector determines whether the factors interact with the covariate. All interactions between factors had q-values greater than individual factors, demonstrating that combining carbon and water cycle phenomenon-related factors increases the explanatory power of GBE (Fig. 8 ). The strongest interactions (greater than 0.75) were observed between AI and PREC (0.83), DEM and PREC (0.82), AI and PET, PET and PREC (0.78), and between DEM and PET (0.75). It was not verified that there were any interactions between AI with SLOPE, AI with DEM, and SLOPE with PREC, but SLOPE interacted with DEM (0.59) and PET (0.53). The risk detector computes the average values within ecosystems and determines if there are significant differences between the means of each stratum. Considering the country's extent, the mean value of GE corresponds to 51.47%. The mean value of BE referring to wetlands was 66.36% and 87.23% for freshwater ecosystems. All statistics were significantly different, considering a significant threshold of 0.05. 3.4. Inland GBE physical accounting: extent and condition In 2018, forests, freshwater bodies, and wetlands covered around 62% of the territory of Continental Portugal. Forests occupied a total area of 5287284.95 hectares covering 60.32% of the total country area, freshwater bodies filled an area of 128391.01 ha, constituting roughly 1.7%, and wetlands 0.33% occupying an area of 21792.51 ha. The condition of forests was assessed with the mean forest crown cover density (FCD), and the mean FWPI a proxy for the status of freshwater and wetlands. The GBE status assessment is calculated by pixel, meaning that, in an area of 100 square meters classified as forests, 51.47% of this area will account for aboveground biomass. In a pixel classified as freshwater and wetlands, respectively, 87.23% and 66.36% will effectively account for Blue Ecosystems. The area (hectares), proportional extent (%), status (%), and physical accounts (%) of forests, freshwater bodies and wetlands in 2018 at the national scale are presented in Table 5 . Table 5 Accounting for inland GBE extent and condition for Mainland Portugal in 2018. 2018 Forest Freshwater Wetlands Area (ha) 5287285 128391.01 21792.51 Proportional extent (%) 60.36 1.7 0.33 Status (%) 51.47 87.23 66.36 Physical accounts (%) 31.07 1.48 0.22 4. Discussion We proposed a novel methodology that relies on RE techniques to disassemble TCD and WWPI datasets from the Copernicus products portfolio and redesign and reproduce them introducing new data and methods to enhance usability. Besides, our study delivered a hands-on approach to estimating climate-driven stock accounts of forest, freshwater, and wetland ecosystems. The findings on the drivers of biomass density and water availability shed new light on the potential implications of land cover conversions and climate conditions on inland GBE. This information is critical, given that these ecosystems are biodiversity hotspots capable of supporting a wide range of livelihoods, and able to provide a variety of goods and services. 4.1. Cloud-based automated classification of inland GBE Open-access geospatial data streams and open-source software have made frequent, large-scale, high-resolution land cover mapping a reality (Arruda et al., 2021 )​. These recent breakthroughs in acquiring and processing big earth data have made large-scale, high-resolution land cover classifications more viable for a diverse range of user groups, organizations, and scholars (DeLancey et al., 2019 ). GEE has evolved in recent years as an appealing high-performance computing platform for the cloud-based processing of petabytes of earth data, offering computational capability for planetary-scale data processing, as well as the ability to apply well-known ML algorithms (Gorelick et al., 2017 ). The separation between GE and BE was implemented primarily to test the significant improvements in the accuracy and utility of such analyses (Koma et al., 2021 ). The fundamental motivation for segregating GE and BE into distinct categories arises from their inherently disparate spectral characteristics. Vegetation, characterized by its chlorophyll content and unique reflectance properties, displays pronounced spectral features in the green and near-infrared regions (Gitelson et al., 2003 ). Conversely, water bodies, including wetlands, exhibit spectral absorption patterns in the same regions (Khalid et al., 2021 ). These distinctive spectral properties lay the foundation for specialized classification approaches. The classification analysis has underscored the improvement in classification accuracy achieved through this separation. The Kappa coefficient in both classifications showed compelling results. This enhancement in accuracy is not merely an academic achievement but has practical implications for informed decision-making (Jagannathan & Divya, 2021 ). Beyond accuracy, the separation of GE and BE enables tailored analysis for specific applications. For instance, GE classification supports studies related to vegetation health, land-use planning, and forest management, while BE classification aids in monitoring water bodies, wetland conservation, and hydrological assessments (Kulawardhana, 2011 ). The F1 scores, indicative of the models' robustness in correctly identifying GE and BE categories, further emphasize the advantages of this separation (Alem & Kumar, 2022 ). GE achieved an F1 score of 0.93, while BE classification surpassed an F1 score of 0.94. These high F1 scores underscore the models' ability to effectively distinguish between these two categories. A stand-alone ML classification model for mapping GE and BE has significant implications for environmental assessment, empowering researchers to conduct focused analyses on ecosystems, and supporting decision-makers in crafting sustainable land management and conservation strategies (Chignell et al., 2018 ; Z. Yang et al., 2021 ). Furthermore, the outputs of these classification tasks, the FCM from GE mapping and FWM from the BE were utilized to delineate ecosystem assets and quantify physical measurements of inland GBE (United Nations et al., 2021 ). 4.2. From Reverse Engineering to satellite-based climate-driven models There is an emerging trend involving using SEO data to quantify ecosystem structure and functional traits which serve as more accurate indicators of ES than simplistic cover-based proxies (Schirpke et al., 2023). Their combination with regression models contributes to a comprehensive understanding of the social-ecological processes influencing the provision of ES and yields reliable data and results (C. Ramirez-Reyes et al., 2019 ) and thereby, provide cost-effective data to supply ecological accounts (Vargas et al., 2019 ). The WWPI and TCD have the potential to provide data to build ecological-economic accounts for GBE. Redesigning and reproducing these datasets focused on forest, freshwater and wetlands evidenced the interdependency of some ecosystem functions, such as biomass density, water yield and soil moisture water content, and climatic and topographic drivers, as shown by Mengist et al. ( 2021 ). Topography and climate are key elements in explaining changes in vegetation abundance at the local (i.e. pixel, plot, or parcel scale) and regional levels (Mpakairi et al., 2022 ). FCD predictions were better than FWPI when comparing the results of R-squared, regression errors and residuals. The independent variables considered in both models were topographic characteristics (SLOPE and DEM) and climatic factors (PREC, PET and AI). Using the same factors to estimate SSH within carbon and water cycles demonstrated the intrinsic relationships between aquatic and terrestrial ecosystems (Cantonati et al., 2020 ; Dwire et al., 2018 ; Hose et al., 2021 ; Huang et al., 2020 ; Kløve et al., 2014 ; Koit et al., 2021 ). The aboveground biomass density represented by FCD estimates varies spatially and temporally according to ecosystem types, altitude, latitude, climatic conditions, and water presence (Tabacchi et al., 2011 ). GE supports biodiversity and many livelihoods around lakes, perennial streams, and wetlands, which benefit from BE’s water and nutrient availability (Harris, 1984 ). Mpakairi et al. ( 2022 ) have shown that the density of healthy trees decreased due to limited water and nutrient availability, with water presence driving about 40% of biological processes in forests. Regarding the contribution of each variable to the models, all of them showed to be relevant, a part of SLOPE presenting very low importance to both predictions. DEM and PREC were the factors that most contributed to the FWPI model, followed by PET, AI, and SLOPE. For the FCD model, PET and AI contributed the most, followed by PREC, DEM, and SLOPE. The results are aligned with those presented in the studies of Crowther et al. ( 2015 ) and Madrigal-González et al. ( 2023 ), analysing the influence of climatic and topographic conditions in assessing GBE. These findings expanded our understanding of nature stock climate drivers, landscape behaviours, and increasing accuracy in ecosystem accounts. Therefore, unleashing the capacity to analyze the status and conditions of GBE helps to execute integrated natural resource management. The geographical detector analysis used to verify SSH enriched the knowledge of the drivers of GE and BE. The factor detectors indicate that both response variables are stratified by the exploratory variables, with more evidence of stratification shown by PREC, PET and AI. DEM showed a higher association with the BE than the GE model. Lastly, the physical accounts of inland GBE were assessed using the EA methodological framework. The national GBE extent account represented their magnitude in terms of area in hectares, while the GBE condition described their integrity and status (Czúcz et al., 2021 ). The aboveground biomass in the forest was measured by calculating the area classified as forest, and the mean estimated value of FCD was applied to assess GE condition. BE were evaluated by determining the extent of freshwater bodies and wetlands and the mean estimated values of FWPI for each ecosystem. The physical accounts evidenced that in 2018, 33% of the total inland area of Portugal was represented by forest, freshwater and wetlands ecosystems. The physical accounts offer means to understand how sustainable the production is through measurements of natural stocks and yields, providing baselines for the assessment of ecosystems’ capacity to deliver goods and services at multiscale (Keith et al., 2020 ). 4.3. Implications for future research There are challenges with reverse engineering, especially when it comes to digital products (Chikofsky & Cross II, 2002). Other issues outside understanding the dataset itself include determining how it was created when we have little or no access to inputs, intermediate outputs, or even the code, algorithms, and tools used in each task (Wood, 2009 ). Without clear documentation of the model architecture, hyperparameters, and training process, it may be difficult for others to reproduce methods or understand the model's behaviour, and how the model arrives at its predictions. The transparency through reporting the entire ML pipeline, including data preprocessing, model architecture, hyperparameters, and evaluation metrics enables others to reproduce and validate the research findings and improve the robustness and reliability of ML models, leading to better accuracy and more trustworthy results. Additionally, continuous monitoring and updating of models as new data becomes available can help maintain their performance over time. Nowadays, a wide range of technologies exists to aid in the forward/reverse engineering process, especially in the production of digital products. A good improvement would be to provide more information about these products' development and make version control available to the public. This practice helps in keeping track of dataset modifications, comparing several dataset versions, and revert to whichever previous version is required. Future research would further take these products apart to continue improving functionality, introducing new data and testing algorithms, as well as replicating to other areas in the globe. Additionally, it would be interesting to explore other products developed by Copernicus, which comprises a wide variety of products, such as Soil Water Index, Imperviousness Density and Impervious Built-up, Dominant Leaf Type and Forest Type, Coastal Zones and Protected Areas. These datasets are valuable resources for a wide range of geospatial services and environmental analysis. 5. Conclusion This research demonstrated the potential of the TCD and WWPI products provided by Copernicus as fundamental data for ecological-economic studies and presented an approach to overcome limitations and enhance their usability. Besides, it underscores the significance of climatic (PREC, PET, AT) and topographic (DEM and SLOPE) variables influencing the status and conditions of GBE, by analyzing spatial relationships in the TCD and WWPI predictions. RE proved to be a valuable tool for innovation and problem-solving while examining the product's functionality and behaviour to determine how it works and identifying any vulnerabilities, weaknesses, or areas for improvement. The outcomes demonstrate the synergy of combining global open-access data with open-source technologies in creating predictive models that offer actionable insights, thereby facilitating evidence-based decision-making across various environmental and societal challenges. Declarations Funding Open Access funding was provided by NOVA Information Management School NOVA University Lisbon. This study was funded by the research project MaSOT – Mapping Ecosystem Services from Earth Observations, funded by the Portuguese Science Foundation – FCT [EXPL/CTA-AMB/0165/2021], and by national funds through FCT (Fundação para a Ciência e a Tecnologia), under the project - UIDB/04152/2020 - Centro de Investigação em Gestão de Informação (MagIC)/NOVA IMS. Author Contribution B.A.: Conceptualization, Methodology, Investigation, Formal analysis, Writing - original draft, preparation, Writing - review and editing. L.M.: Conceptualization, Methodology, Investigation, Writing - review and editing. P.S.: Conceptualization, Methodology, Investigation, Writing - review and editing. P.C.: Conceptualization, Methodology, Investigation, Writing - review and editing. Acknowledgement This study was supported by the research project MaSOT – Mapping Ecosystem Services from Earth Observations, funded by the Portuguese Science Foundation – FCT [EXPL/CTA-AMB/0165/2021], and by national funds through FCT (Fundação para a Ciência e a Tecnologia), under the project - UIDB/04152/2020 - Centro de Investigação em Gestão de Informação (MagIC)/NOVA IMS. We are grateful to the Institut français du Portugal for fostering Franco-Portuguese collaborations through scientific periods in France (BIC-2023), and the IUEM for hosting. References Alem, A., & Kumar, S. (2022). Applied Artificial Intelligence End-to-End Convolutional Neural Network Feature Extraction for Remote Sensed Images Classification . https://doi.org/10.1080/08839514.2022.2137650 Almeida, B., & Cabral, P. (2021). Water yield modelling, sensitivity analysis and validation: A study for Portugal. ISPRS International Journal of Geo-Information , 10 (8). https://doi.org/10.3390/ijgi10080494 Almeida, B., & Cabral, P. (2023). Data-Driven Modelling of Freshwater Ecosystems: A Multiscale Framework Based on Global Geospatial Data. Proceedings of the 9th International Conference on Geographical Information Systems Theory, Applications and Management , 2023-April , 104–111. https://doi.org/10.5220/0012037800003473 Arruda, V. L. S., Piontekowski, V. J., Alencar, A., Pereira, R. S., & Matricardi, E. A. T. (2021). An alternative approach for mapping burn scars using Landsat imagery, Google Earth Engine, and Deep Learning in the Brazilian Savanna. Remote Sensing Applications: Society and Environment , 22 , 100472. https://doi.org/10.1016/j.rsase.2021.100472 Bateman, I. J., & Mace, G. M. (2020). The natural capital framework for sustainably efficient and equitable decision making. Nature Sustainability , 3 (10), 776–783. https://doi.org/10.1038/s41893-020-0552-3 Beck, H. E., Zimmermann, N. E., McVicar, T. R., Vergopolan, N., Berg, A., & Wood, E. F. (2018). Present and future Köppen-Geiger climate classification maps at 1-km resolution. Scientific Data , 5 (1), 180214. https://doi.org/10.1038/sdata.2018.214 Belo-Pereira, M., Dutra, E., & Viterbo, P. (2011). Evaluation of global precipitation data sets over the Iberian Peninsula. Journal of Geophysical Research: Atmospheres , 116 (D20). https://doi.org/https://doi.org/10.1029/2010JD015481 Bordt, M. (2018). Discourses in Ecosystem Accounting: A Survey of the Expert Community. Ecological Economics , 144 , 82–99. https://doi.org/10.1016/J.ECOLECON.2017.06.032 Braun, D., Damm, A., Hein, L., Petchey, O. L., & Schaepman, M. E. (2018). Spatio-temporal trends and trade-offs in ecosystem services: An Earth observation based assessment for Switzerland between 2004 and 2014. Ecological Indicators , 89 , 828–839. https://doi.org/10.1016/j.ecolind.2017.10.016 Breiman, L. (2001). Random forests. Machine Learning , 45 (1), 5–32. https://doi.org/10.1023/A:1010933404324 Cantonati, M., Stevens, L. E., Segadelli, S., Springer, A. E., Goldscheider, N., Celico, F., Filippini, M., Ogata, K., & Gargini, A. (2020). Ecohydrogeology: The interdisciplinary convergence needed to improve the study and stewardship of springs and other groundwater-dependent habitats, biota, and ecosystems. Ecological Indicators , 110 . https://doi.org/10.1016/j.ecolind.2019.105803 Chen, G., Singh, K. K., Lopez, J., & Zhou, Y. (2020). Tree canopy cover and carbon density are different proxy indicators for assessing the relationship between forest structure and urban socio-ecological conditions. Ecological Indicators , 113 , 106279. https://doi.org/10.1016/J.ECOLIND.2020.106279 Chignell, S. M., Luizza, M. W., Skach, S., Young, N. E., & Evangelista, P. H. (2018). An integrative modeling approach to mapping wetlands and riparian areas in a heterogeneous Rocky Mountain watershed. Remote Sensing in Ecology and Conservation , 4 (2), 150–165. https://doi.org/https://doi.org/10.1002/rse2.63 Chikofsky, E., & Cross II, J. H. (2002). Reverse Engineering. In Encyclopedia of Software Engineering . https://doi.org/https://doi.org/10.1002/0471028959.sof285 Copernicus Programme. (2023). CLC 2018 — Copernicus Land Monitoring Service . https://land.copernicus.eu/pan-european/corine-land-cover/clc2018 Cord, A. F., Brauman, K. A., Chaplin-Kramer, R., Huth, A., Ziv, G., & Seppelt, R. (2017). Priorities to Advance Monitoring of Ecosystem Services Using Earth Observation. Trends in Ecology & Evolution , 32 (6), 416–428. https://doi.org/10.1016/j.tree.2017.03.003 Costa, H., Benevides, P., Moreira, F. D., Moraes, D., & Caetano, M. (2022). Spatially Stratified and Multi-Stage Approach for National Land Cover Mapping Based on Sentinel-2 Data and Expert Knowledge. Remote Sensing 2022, Vol. 14, Page 1865 , 14 (8), 1865. https://doi.org/10.3390/RS14081865 Crowther, T. W., Glick, H. B., Covey, K. R., Bettigole, C., Maynard, D. S., Thomas, S. M., Smith, J. R., Hintler, G., Duguid, M. C., Amatulli, G., Tuanmu, M. N., Jetz, W., Salas, C., Stam, C., Piotto, D., Tavani, R., Green, S., Bruce, G., Williams, S. J., … Bradford, M. A. (2015). Mapping tree density at a global scale. Nature , 525 (7568), 201–205. https://doi.org/10.1038/nature14967 Czúcz, B., Keith, H., Driver, A., Jackson, B., Nicholson, E., & Maes, J. (2021). A common typology for ecosystem characteristics and ecosystem condition variables. One Ecosystem , 6 , 1–16. https://doi.org/10.3897/oneeco.6.e58218 del Río-Mena, T., Willemen, L., Vrieling, A., & Nelson, A. (2023). How remote sensing choices influence ecosystem services monitoring and evaluation results of ecological restoration interventions. Ecosystem Services , 64 . https://doi.org/10.1016/J.ECOSER.2023.101565 DeLancey, E. R., Kariyeva, J., Bried, J. T., & Hird, J. N. (2019). Large-scale probabilistic identification of boreal peatlands using Google Earth Engine, open-access satellite data, and machine learning. PLoS ONE , 14 (6), 1–23. https://doi.org/10.1371/journal.pone.0218165 DGT. (2018). Technical specifications of the land use and land cover map for Portugal Mainland . Dwire, K. A., Mellmann-Brown, S., & Gurrieri, J. T. (2018). Potential effects of climate change on riparian areas, wetlands, and groundwater-dependent ecosystems in the Blue Mountains, Oregon, USA. In Climate Services (Vol. 10, pp. 44–52). Elsevier B.V. https://doi.org/10.1016/j.cliser.2017.10.002 Edens, B., Maes, J., Hein, L., Obst, C., Siikamaki, J., Schenau, S., Javorsek, M., Chow, J., Chan, J. Y., Steurer, A., & Alfieri, A. (2022). Establishing the SEEA Ecosystem Accounting as a global standard. Ecosystem Services , 54 , 101413. https://doi.org/10.1016/J.ECOSER.2022.101413 Ellison, D., Morris, C. E., Locatelli, B., Sheil, D., Cohen, J., Murdiyarso, D., Gutierrez, V., Noordwijk, M. van, Creed, I. F., Pokorny, J., Gaveau, D., Spracklen, D. V., Tobella, A. B., Ilstedt, U., Teuling, A. J., Gebrehiwot, S. G., Sands, D. C., Muys, B., Verbist, B., … Sullivan, C. A. (2017). Trees, forests and water: Cool insights for a hot world. Global Environmental Change , 43 , 51–61. https://doi.org/10.1016/J.GLOENVCHA.2017.01.002 ESRI. (2023). ArcGIS Pro - ESRI (Environmental Systems Research Institute) . https://www.arcgis.com/index.html European Environment Agency. (2023). Copernicus Land Monitoring Service . European Spatial Agency. (2023). ESA - Space Science . https://www.esa.int/Science_Exploration/Space_Science Farrell, C., Coleman, L., Kelly-Quinn, M., Obst, C., Eigenraam, M., Norton, D., ODonoghue, C., Kinsella, S., Delargy, O., & Stout, J. (2021). Applying the System of Environmental Economic Accounting-Ecosystem Accounting (SEEA-EA) framework at catchment scale to develop ecosystem extent and condition accounts. One Ecosystem , 6 . https://doi.org/10.3897/oneeco.6.e65582 Fick, S. E., & Hijmans, R. J. (2017). WorldClim 2: new 1-km spatial resolution climate surfaces for global land areas. International Journal of Climatology , 37 (12), 4302–4315. https://doi.org/https://doi.org/10.1002/joc.5086 Fleming, A., O’Grady, A. P., Stitzlein, C., Ogilvy, S., Mendham, D., & Harrison, M. T. (2022). Improving acceptance of natural capital accounting in land use decision making: Barriers and opportunities. Ecological Economics , 200 , 107510. https://doi.org/10.1016/J.ECOLECON.2022.107510 Fonseca, F., de Figueiredo, T., Vilela, Â., Santos, R., de Carvalho, A. L., Almeida, E., & Nunes, L. (2019). Impact of tree species replacement on carbon stocks in a Mediterranean mountain area, NE Portugal. Forest Ecology and Management , 439 , 181–188. https://doi.org/https://doi.org/10.1016/j.foreco.2019.03.002 Gitelson, A. A., Gritz, Y., & Merzlyak, M. N. (2003). Relationships between leaf chlorophyll content and spectral reflectance and algorithms for non-destructive chlorophyll assessment in higher plant leaves. Journal of Plant Physiology , 160 (3), 271–282. https://doi.org/10.1078/0176-1617-00887 Google. (2023). Google Earth Engine . https://earthengine.google.com/ Gorelick, N., Hancher, M., Dixon, M., Ilyushchenko, S., Thau, D., & Moore, R. (2017). Google Earth Engine: Planetary-scale geospatial analysis for everyone. Remote Sensing of Environment , 202 , 18–27. https://doi.org/https://doi.org/10.1016/j.rse.2017.06.031 Gosal, A. S., Evans, P. M., Bullock, J. M., Redhead, J., Charlton, M. B., Cord, A. F., Johnson, A., & Ziv, G. (2022). Understanding the accuracy of modelled changes in freshwater provision over time. Science of The Total Environment , 833 , 155042. https://doi.org/10.1016/j.scitotenv.2022.155042 Harris, L. D. (1984). Island Biogeography Theory and the Preservation of Biotic Diversity . University of Chicago Press. https://doi.org/doi:10.7208/9780226219950 Hein, L., Bagstad, K. J., Obst, C., Edens, B., Schenau, S., Castillo, G., Soulard, F., Brown, C., Driver, A., Bordt, M., Steurer, A., Harris, R., & Caparrós, A. (2020). Progress in natural capital accounting for ecosystems. Science , 367 (6477), 514–515. https://doi.org/10.1126/science.aaz8901 Hose, G. C., Dabovic, J., Nelson, T., Pollitt, L., & Korbel, K. L. (2021). Groundwater Dependent Aquatic and Terrestrial Ecosystems. Reference Module in Earth Systems and Environmental Sciences . https://doi.org/10.1016/B978-0-12-819166-8.00118-3 Hosmer, D. W., & Lemeshow, S. (2000). Interpretation of the Fitted Logistic Regression Model. In Applied Logistic Regression (pp. 47–90). Wiley. https://doi.org/10.1002/0471722146.ch3 Hossain, M. S., & Hashim, M. (2019). Potential of Earth Observation (EO) technologies for seagrass ecosystem service assessments. International Journal of Applied Earth Observation and Geoinformation , 77 , 15–29. https://doi.org/10.1016/j.jag.2018.12.009 Huang, F., Chunyu, X., Zhang, D., Chen, X., & Ochoa, C. G. (2020). A framework to assess the impact of ecological water conveyance on groundwater-dependent terrestrial ecosystems in arid inland river basins. Science of the Total Environment , 709 . https://doi.org/10.1016/j.scitotenv.2019.136155 Huete, A. R. (1988). A soil-adjusted vegetation index (SAVI). Remote Sensing of Environment , 25 (3), 295–309. https://doi.org/10.1016/0034-4257(88)90106-X Huete, A. R., Liu, H. Q., & van Leeuwen, W. J. D. (1997). Use of vegetation indices in forested regions: Issues of linearity and saturation. International Geoscience and Remote Sensing Symposium (IGARSS) , 4 , 1966–1968. https://doi.org/10.1109/IGARSS.1997.609169 Jagannathan, J., & Divya, C. (2021). Deep learning for the prediction and classification of land use and land cover changes using deep convolutional neural network. Ecological Informatics , 65 , 101412. https://doi.org/10.1016/J.ECOINF.2021.101412 JAXA. (2023). ALOS World 3D - 30m . https://www.eorc.jaxa.jp/ALOS/en/dataset/aw3d30/aw3d30_e.htm Ji, L., Geng, X., Sun, K., Zhao, Y., & Gong, P. (2015). Target Detection Method for Water Mapping Using Landsat 8 OLI/TIRS Imagery. Water , 7 (2), 794–817. https://doi.org/10.3390/w7020794 Keith, H., Czúcz, B., Jackson, B., Driver, A., Nicholson, E., & Maes, J. (2020). A conceptual framework and practical structure for implementing ecosystem condition accounts. One Ecosystem , 5 , 1–54. https://doi.org/10.3897/oneeco.5.e58216 Keith, H., Vardon, M., Obst, C., Young, V., Houghton, R., & Mackey, B. (2021). Evaluating nature-based solutions for climate mitigation and conservation requires comprehensive carbon accounting. Science of The Total Environment , 769 , 144341. https://doi.org/10.1016/j.scitotenv.2020.144341 Khalid, H. W., Khalil, R. M. Z., & Qureshi, M. A. (2021). Evaluating spectral indices for water bodies extraction in western Tibetan Plateau. The Egyptian Journal of Remote Sensing and Space Science , 24 (3), 619–634. https://doi.org/10.1016/J.EJRS.2021.09.003 Kienast, F., Bolliger, J., Potschin, M., de Groot, R. S., Verburg, P. H., Heller, I., Wascher, D., & Haines-Young, R. (2009). Assessing Landscape Functions with Broad-Scale Environmental Data: Insights Gained from a Prototype Development for Europe. Environmental Management , 44 (6), 1099–1120. https://doi.org/10.1007/s00267-009-9384-7 Kløve, B., Ala-Aho, P., Bertrand, G., Gurdak, J. J., Kupfersberger, H., Kværner, J., Muotka, T., Mykrä, H., Preda, E., Rossi, P., Uvo, C. B., Velasco, E., & Pulido-Velazquez, M. (2014). Climate change impacts on groundwater and dependent ecosystems. Journal of Hydrology , 518 (PB), 250–266. https://doi.org/10.1016/J.JHYDROL.2013.06.037 Koit, O., Tarros, S., Pärn, J., Küttim, M., Abreldaal, P., Sisask, K., Vainu, M., Terasmaa, J., Retike, I., & Polikarpus, M. (2021). Contribution of local factors to the status of a groundwater dependent terrestrial ecosystem in the transboundary Gauja-Koiva River basin, North-Eastern Europe. Journal of Hydrology , 600 . https://doi.org/10.1016/j.jhydrol.2021.126656 Koma, Z., Seijmonsbergen, A. C., & Kissling, W. D. (2021). Classifying wetland-related land cover types and habitats using fine-scale lidar metrics derived from country-wide Airborne Laser Scanning. Remote Sensing in Ecology and Conservation , 7 (1), 80–96. https://doi.org/https://doi.org/10.1002/rse2.170 Kulawardhana, R. W. (2011). Remote sensing of vegetation: principles, techniques and applications. By Hamlyn G. Jones and Robin A Vaughan. Journal of Vegetation Science , 22 (6), 1151–1153. https://doi.org/https://doi.org/10.1111/j.1654-1103.2011.01319.x Lange, S., Campagne, C. S., Comte, A., Bank, E., Santos-Martín, F., Maes, J., & Burkhard, B. (2022). Progress on ecosystem accounting in Europe. Ecosystem Services , 57 , 101473. https://doi.org/10.1016/J.ECOSER.2022.101473 Liquete, C., Maes, J., Notte, A., & Bidoglio, G. (2011). Securing water as a resource for society: an ecosystem services perspective. Ecohydrology & Hydrobiology , 11 (3–4), 247–259. https://doi.org/10.2478/V10104-011-0044-1 Ludwig, C., Walli, A., Schleicher, C., Weichselbaum, J., & Riffler, M. (2019). A highly automated algorithm for wetland detection using multi-temporal optical satellite data. Remote Sensing of Environment , 224 , 333–351. https://doi.org/10.1016/J.RSE.2019.01.017 Lundqvist, J., & Unver, O. (2018). Alternative pathways to food security and nutrition – water predicaments and human behavior. Water Policy , 20 (5), 871–884. https://doi.org/10.2166/WP.2018.171 Madrigal-González, J., Calatayud, J., Ballesteros-Cánovas, J. A., Escudero, A., Cayuela, L., Marqués, L., Rueda, M., Ruiz-Benito, P., Herrero, A., Aponte, C., Sagardia, R., Plumptre, A. J., Dupire, S., Espinosa, C. I., Tutubalina, O. V., Myint, M., Pataro, L., López-Sáez, J., Macía, M. J., … Stoffel, M. (2023). Global patterns of tree density are contingent upon local determinants in the world’s natural forests. Communications Biology 2023 6:1 , 6 (1), 1–6. https://doi.org/10.1038/s42003-023-04419-8 Mairota, P., Cafarelli, B., Didham, R. K., Lovergine, F. P., Lucas, R. M., Nagendra, H., Rocchini, D., & Tarantino, C. (2015). Challenges and opportunities in harnessing satellite remote-sensing for biodiversity monitoring. Ecological Informatics , 30 , 207–214. https://doi.org/10.1016/j.ecoinf.2015.08.006 Mengist, W., Soromessa, T., & Feyisa, G. L. (2021). Monitoring Afromontane forest cover loss and the associated socio-ecological drivers in Kaffa biosphere reserve, Ethiopia. Trees, Forests and People , 6 . https://doi.org/10.1016/j.tfp.2021.100161 Mpakairi, K. S., Dube, T., Dondofema, F., & Dalu, T. (2022). Spatio–temporal variation of vegetation heterogeneity in groundwater dependent ecosystems within arid environments. Ecological Informatics , 69 , 101667. https://doi.org/10.1016/j.ecoinf.2022.101667 Mulatu, K. A., Mora, B., Kooistra, L., & Herold, M. (2017). Biodiversity monitoring in changing tropical forests: A review of approaches and new opportunities. In Remote Sensing (Vol. 9, Issue 10). MDPI AG. https://doi.org/10.3390/rs9101059 Newbold, T., Hudson, L. N., Hill, S. L. L., Contu, S., Lysenko, I., Senior, R. A., Börger, L., Bennett, D. J., Choimes, A., Collen, B., Day, J., De Palma, A., Díaz, S., Echeverria-Londoño, S., Edgar, M. J., Feldman, A., Garon, M., Harrison, M. L. K., Alhusseini, T., … Purvis, A. (2015). Global effects of land use on local terrestrial biodiversity. Nature , 520 (7545), 45–50. https://doi.org/10.1038/nature14324 Nieves-Chinchilla, J., Martínez, R., Farjas, M., Tubio-Pardavila, R., Cruz, D., & Gallego, M. (2018). Reverse engineering techniques to optimize facility location of satellite ground stations on building roofs. Automation in Construction , 90 , 156–165. https://doi.org/10.1016/J.AUTCON.2018.02.019 Oliveira, M., Novo, M. emilia, & Lobo-Ferreira, J.-P. (2007). Models to predict the impact of the climate changes on aquifer recharge. Water in Celtic Countries: Quantity, Quality and Climate Variability (Proceedings of the Fourth InterCeltic Colloquium on Hydrology and Management of Water Resources , 103–110. https://iahs.info/uploads/dms/13794.15-103-110-05-IC310-23-Oliveira.pdf Ramirez-Reyes, C., Brauman, K. A., Chaplin-Kramer, R., Galford, G. L., Adamo, S. B., Anderson, C. B., Anderson, C., Allington, G. R. H., Bagstad, K. J., Coe, M. T., Cord, A. F., Dee, L. E., Gould, R. K., Jain, M., Kowal, V. A., Muller-Karger, F. E., Norriss, J., Potapov, P., Qiu, J., … Wright, T. M. (2019). Reimagining the potential of Earth observations for ecosystem service assessments. Science of The Total Environment , 665 , 1053–1063. https://doi.org/10.1016/j.scitotenv.2019.02.150 Rijal, S., Rimal, B., Acharya, R. P., & Stork, N. E. (2021). Land use/land cover change and ecosystem services in the Bagmati River Basin, Nepal. Environmental Monitoring and Assessment , 193 (10), 651. https://doi.org/10.1007/s10661-021-09441-z Rouse, J., Haas, R., Schell J., & Deering D. (1973). Monitoring Vegetation Systems in the Great Plains with ERTS. Third ERTS Symposium, NASA . Santos, F. L. M., Couto, F. T., Dias, S. S., Ribeiro, N. A., & Salgado, R. (2023). Vegetation fuel characterization using machine learning approach over southern Portugal. Remote Sensing Applications: Society and Environment , 101017. https://doi.org/10.1016/J.RSASE.2023.101017 Seddon, N., Smith, A., Smith, P., Key, I., Chausson, A., Girardin, C., House, J., Srivastava, S., & Turner, B. (2021). Getting the message right on nature-based solutions to climate change. Global Change Biology , 27 (8), 1518–1546. https://doi.org/10.1111/GCB.15513 Tabacchi, G., Di Cosmo, L., & Gasparini, P. (2011). Aboveground tree volume and phytomass prediction equations for forest species in Italy. European Journal of Forest Research , 130 (6), 911–934. https://doi.org/10.1007/S10342-011-0481-9/TABLES/8 United Nations et al. (2021). System of Environmental-Economic Accounting-Ecosystem Accounting (SEEA-EA) . https://seea.un.org/ecosystem-accounting. Vallero, D. (2014). Grand Challenges. Fundamentals of Air Pollution , 953–961. https://doi.org/10.1016/B978-0-12-401733-7.00033-5 Vargas, L., Willemen, L., & Hein, L. (2019). Assessing the Capacity of Ecosystems to Supply Ecosystem Services Using Remote Sensing and An Ecosystem Accounting Approach. Environmental Management , 63 (1), 1–15. https://doi.org/10.1007/s00267-018-1110-x Wang, J.-F., Zhang, T.-L., & Fu, B.-J. (2016). A measure of spatial stratified heterogeneity. Ecological Indicators , 67 , 250–256. https://doi.org/10.1016/j.ecolind.2016.02.052 Willcock, S., Martínez-López, J., Hooftman, D. A. P., Bagstad, K. J., Balbi, S., Marzo, A., Prato, C., Sciandrello, S., Signorello, G., Voigt, B., Villa, F., Bullock, J. M., & Athanasiadis, I. N. (2018). Machine learning for ecosystem services. Ecosystem Services , 33 , 165–174. https://doi.org/10.1016/J.ECOSER.2018.04.004 Wood, W. H. (2009). Computational Representations of Function in Engineering Design. Philosophy of Technology and Engineering Sciences , 543–564. https://doi.org/10.1016/B978-0-444-51667-1.50024-0 Yang, J., Wang, J., Liao, X., Tao, H., & Li, Y. (2022). Chain modeling for the biogeochemical nexus of cadmium in soil-rice-human health system. Environment International , 167 , 107424. https://doi.org/10.1016/j.envint.2022.107424 Yang, Z., Bai, J., & Zhang, W. (2021). Mapping and assessment of wetland conditions by using remote sensing images and POI data. Ecological Indicators , 127 , 107485. https://doi.org/10.1016/J.ECOLIND.2021.107485 Zomer, R. J., Xu, J., & Trabucco, A. (2022). Version 3 of the Global Aridity Index and Potential Evapotranspiration Database. Scientific Data , 9 (1), 409. https://doi.org/10.1038/s41597-022-01493-1 Additional Declarations No competing interests reported. Supplementary Files SupplementaryMaterials.docx Cite Share Download PDF Status: Posted Version 1 posted You are reading this latest preprint version Research Square lets you share your work early, gain feedback from the community, and start making changes to your manuscript prior to peer review in a journal. As a division of Research Square Company, we’re committed to making research communication faster, fairer, and more useful. We do this by developing innovative software and high quality services for the global research community. Our growing team is made up of researchers and industry professionals working together to solve the most critical problems facing scientific publishing. Also discoverable on Platform About Our Team In Review Editorial Policies Advisory Board Help Center Resources Author Services Accessibility API Access RSS feed Manage Cookie Preferences © Research Square 2026 | ISSN 2693-5015 (online) Privacy Policy Terms of Service Do Not Sell My Personal Information {"props":{"pageProps":{"initialData":{"identity":"rs-4681296","acceptedTermsAndConditions":true,"allowDirectSubmit":true,"archivedVersions":[],"articleType":"Research Article","associatedPublications":[],"authors":[{"id":347547082,"identity":"d64048a5-c09e-49b1-ae5d-1e76014137d7","order_by":0,"name":"Bruna Almeida","email":"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAZAAAAAyAQMAAABI0h/eAAAABlBMVEX///8AAABVwtN+AAAACXBIWXMAAA7EAAAOxAGVKw4bAAABCElEQVRIiWNgGAWjYBACA3QBGX4QmVBAghYeyQaQFgxxfFoMDmAVRwBz6TNmD34w2Mnrzu4x/sDYdofH+PzqxA8PDBjk+cUOYNVi2ZdjbtjDkGy47c4ZMwnGtmc8ZjfebpYAOsxw5uwE7A47w2MmwcNwgHHbjRwzBsa2w0AtZzeAtCQY3MatRfIPwwF7oBaQww7zGM84u/kHIS3SQFsSgVoMJEBaDPh7t+G1xbKHrUxaxiA5edudY2USCecO80jc4N1mkWAggdMv5jzM2yTfVNjZbrvdvPnDh7LDcvz9Zzff/FFhI88vjV0L1HlALAHEYDUSEBKPchiAq+E/QITqUTAKRsEoGEkAAIiWWowUHAjsAAAAAElFTkSuQmCC","orcid":"","institution":"NOVA Information Management School (NOVA IMS), Universidade Nova de Lisboa","correspondingAuthor":true,"prefix":"","firstName":"Bruna","middleName":"","lastName":"Almeida","suffix":""},{"id":347547087,"identity":"7c0f3758-0681-4dde-bb8a-b78655824734","order_by":1,"name":"Luís Monteiro","email":"","orcid":"","institution":"NOVA Information Management School (NOVA IMS), Universidade Nova de Lisboa","correspondingAuthor":false,"prefix":"","firstName":"Luís","middleName":"","lastName":"Monteiro","suffix":""},{"id":347547089,"identity":"231535bd-5575-49a3-b2a7-958dd4a8e4af","order_by":2,"name":"Pierre Scemama","email":"","orcid":"","institution":"IFREMER, Univ Brest, CNRS, UMR AMURE","correspondingAuthor":false,"prefix":"","firstName":"Pierre","middleName":"","lastName":"Scemama","suffix":""},{"id":347547090,"identity":"93298b3d-74af-4c1c-a618-9c8d6eaa8171","order_by":3,"name":"Pedro Cabral","email":"","orcid":"","institution":"Nanjing University of Information Science and Technology","correspondingAuthor":false,"prefix":"","firstName":"Pedro","middleName":"","lastName":"Cabral","suffix":""}],"badges":[],"createdAt":"2024-07-03 15:01:58","currentVersionCode":1,"declarations":"","doi":"10.21203/rs.3.rs-4681296/v1","doiUrl":"https://doi.org/10.21203/rs.3.rs-4681296/v1","draftVersion":[],"editorialEvents":[],"editorialNote":"","failedWorkflow":false,"files":[{"id":63763856,"identity":"e4b995db-3ca7-4086-b074-966e06e24212","added_by":"auto","created_at":"2024-09-02 06:50:20","extension":"png","order_by":1,"title":"Figure 1","display":"","copyAsset":false,"role":"figure","size":336483,"visible":true,"origin":"","legend":"\u003cp\u003eGeographic location of the study area within Portugal’s mainland.\u003c/p\u003e","description":"","filename":"1.png","url":"https://assets-eu.researchsquare.com/files/rs-4681296/v1/91d65ba4afbb22df37be0d6d.png"},{"id":63765389,"identity":"2ca47c94-e798-4a63-a062-9e3c5c3fd058","added_by":"auto","created_at":"2024-09-02 07:06:18","extension":"png","order_by":2,"title":"Figure 2","display":"","copyAsset":false,"role":"figure","size":237079,"visible":true,"origin":"","legend":"\u003cp\u003eMethodology scheme depicting the estimations of ecosystem extent and condition.\u003c/p\u003e","description":"","filename":"2.png","url":"https://assets-eu.researchsquare.com/files/rs-4681296/v1/6c59106065d33e61a5b45ed0.png"},{"id":63763848,"identity":"ddbea554-8012-497b-824e-14ad62f883c2","added_by":"auto","created_at":"2024-09-02 06:50:18","extension":"png","order_by":3,"title":"Figure 3","display":"","copyAsset":false,"role":"figure","size":417700,"visible":true,"origin":"","legend":"\u003cp\u003eGBE mapping for Mainland Portugal.\u003c/p\u003e","description":"","filename":"3.png","url":"https://assets-eu.researchsquare.com/files/rs-4681296/v1/d6ffcd5e6967c98179bea512.png"},{"id":63763857,"identity":"60282595-41a0-444c-9b00-da0ce2d7c0ca","added_by":"auto","created_at":"2024-09-02 06:50:20","extension":"png","order_by":4,"title":"Figure 4","display":"","copyAsset":false,"role":"figure","size":30289,"visible":true,"origin":"","legend":"\u003cp\u003eMean % weight of feature Importance for Forest Cover Mask (FCM) and Freshwater \u0026amp; Wetlands Mask (FWM).\u003c/p\u003e","description":"","filename":"4.png","url":"https://assets-eu.researchsquare.com/files/rs-4681296/v1/c81623022740c72df7bcae84.png"},{"id":63764712,"identity":"41c8dd91-4ecc-4b55-a58b-9acb7bec040b","added_by":"auto","created_at":"2024-09-02 06:58:18","extension":"png","order_by":5,"title":"Figure 5","display":"","copyAsset":false,"role":"figure","size":179783,"visible":true,"origin":"","legend":"\u003cp\u003eScatterplot (top) and residual plot (bottom) of the Forest Cover Density (FCD) model.\u003c/p\u003e","description":"","filename":"5.png","url":"https://assets-eu.researchsquare.com/files/rs-4681296/v1/c9d10a02d13fde022af26a16.png"},{"id":63763854,"identity":"be60196a-d985-45a1-863a-e823f52eaec9","added_by":"auto","created_at":"2024-09-02 06:50:18","extension":"png","order_by":6,"title":"Figure 6","display":"","copyAsset":false,"role":"figure","size":195639,"visible":true,"origin":"","legend":"\u003cp\u003eScatterplot (top) and residual plot (bottom) of the Freshwater \u0026amp; Wetness Probability Index (FWPI) model.\u003c/p\u003e","description":"","filename":"6.png","url":"https://assets-eu.researchsquare.com/files/rs-4681296/v1/beea42f5b3347b28cbf59ef3.png"},{"id":63763849,"identity":"cd997fca-692e-4e0a-bd3d-2632e3c85a73","added_by":"auto","created_at":"2024-09-02 06:50:18","extension":"png","order_by":7,"title":"Figure 7","display":"","copyAsset":false,"role":"figure","size":23088,"visible":true,"origin":"","legend":"\u003cp\u003ePredictor importance for Freshwater \u0026amp; Wetness Probability Index (FWPI) and Forest Cover Density (FCD) regression models.\u003c/p\u003e","description":"","filename":"7.png","url":"https://assets-eu.researchsquare.com/files/rs-4681296/v1/c5495775398437034e2320bc.png"},{"id":63763852,"identity":"afb713f9-22c7-4f2a-86fc-4ed86fab3166","added_by":"auto","created_at":"2024-09-02 06:50:18","extension":"png","order_by":8,"title":"Figure 8","display":"","copyAsset":false,"role":"figure","size":20145,"visible":true,"origin":"","legend":"\u003cp\u003eGBE Factor interaction.\u003c/p\u003e","description":"","filename":"8.png","url":"https://assets-eu.researchsquare.com/files/rs-4681296/v1/c5b49b31352dc91f3ec0bb32.png"},{"id":66788261,"identity":"9c8eb455-c642-4236-b2f7-3f895060beab","added_by":"auto","created_at":"2024-10-16 13:17:12","extension":"pdf","order_by":0,"title":"","display":"","copyAsset":false,"role":"manuscript-pdf","size":2075293,"visible":true,"origin":"","legend":"","description":"","filename":"manuscript.pdf","url":"https://assets-eu.researchsquare.com/files/rs-4681296/v1/835f6f20-ca8f-489b-a782-180d02c32fc9.pdf"},{"id":63763853,"identity":"2be1bf56-a5a0-4602-bebf-b7ac6d388375","added_by":"auto","created_at":"2024-09-02 06:50:18","extension":"docx","order_by":1,"title":"","display":"","copyAsset":false,"role":"supplement","size":522685,"visible":true,"origin":"","legend":"","description":"","filename":"SupplementaryMaterials.docx","url":"https://assets-eu.researchsquare.com/files/rs-4681296/v1/9830efef74fa71b9fecbd660.docx"}],"financialInterests":"No competing interests reported.","formattedTitle":"Climate-driven stock accounts of inland Green-Blue Ecosystems: Coupling the carbon and water cycles via Reverse Engineering and Geodetectors","fulltext":[{"header":"Highlights","content":"\u003cul\u003e\n \u003cli\u003eClassification tasks measured GBE extent, regression estimated GBE condition\u003c/li\u003e\n \u003cli\u003eBlue Ecosystems classification top three feature importance: DSM, NDWI and SLOPE\u003c/li\u003e\n \u003cli\u003eClimate-driven stock accounts of GBE were analyzed through geodetectors\u003c/li\u003e\n \u003cli\u003eThe most important drivers for Green Ecosystems were PREC, AI, and PET\u0026nbsp;\u003c/li\u003e\n \u003cli\u003eThe top three drivers for Blue Ecosystems were PREC, AI, and DEM\u003cbr\u003e\u0026nbsp;\u003c/li\u003e\n\u003c/ul\u003e"},{"header":"1. Introduction","content":"\u003cp\u003eUnderstanding the terrestrial biosphere\u0026rsquo;s functioning requires an assessment of the magnitude and distribution of Green-Blue Ecosystems (GBE) (Crowther et al., \u003cspan citationid=\"CR18\" class=\"CitationRef\"\u003e2015\u003c/span\u003e). Green Ecosystems (GE) support a considerable amount of global biodiversity, playing an important role in biogeochemical cycles, and providing several Ecosystem Services (ES) such as water quantity and quality regulation, food production, raw materials provisioning, genetic resources, and carbon storage and sequestration (Mulatu et al., \u003cspan citationid=\"CR64\" class=\"CitationRef\"\u003e2017\u003c/span\u003e). Forest density influences soil and water retention rates, flood regulation, as well as competitive dynamics and habitat suitability for a variety of plant and animal species (Rijal et al., \u003cspan citationid=\"CR69\" class=\"CitationRef\"\u003e2021\u003c/span\u003e). Blue Ecosystems (BE) provide the backbone of ecosystems, enabling agriculture, and influencing urban development decisions (Newbold et al., \u003cspan citationid=\"CR65\" class=\"CitationRef\"\u003e2015\u003c/span\u003e). As forests are intricately connected to rainfall and water availability, changes in vegetation cover, as well as poorly planned and managed forests, have significant impacts on the carbon and water cycles, leading to pressures on the environment (Liquete et al., \u003cspan citationid=\"CR57\" class=\"CitationRef\"\u003e2011\u003c/span\u003e). Extreme weather events and soil and water contamination are perceived as the top risk, directly impacting natural capital and biodiversity, and consequently decreasing human well-being (Lundqvist \u0026amp; Unver, \u003cspan citationid=\"CR59\" class=\"CitationRef\"\u003e2018\u003c/span\u003e). Knowledge of the interactions between carbon and water cycles is vital for predicting ecosystem responses to climate change, as well as the influence on natural capital (Seddon et al., \u003cspan citationid=\"CR72\" class=\"CitationRef\"\u003e2021\u003c/span\u003e). Despite significant advances, the coupling mechanisms of these cycles under socio-economic and climatic pressures remain poorly understood (Keith et al., \u003cspan citationid=\"CR49\" class=\"CitationRef\"\u003e2021\u003c/span\u003e).\u003c/p\u003e \u003cp\u003eEcosystem Accounting (EA) encompass a statistical framework for quantifying and valuing the interactions between nature, society, and the economy, providing physical and monetary measurements of ecosystem assets (Bateman \u0026amp; Mace, \u003cspan citationid=\"CR5\" class=\"CitationRef\"\u003e2020\u003c/span\u003e). In March 2021, the United Nations Statistical Commission adopted the System of Environmental-Economic Accounting Ecosystem-Accounting (SEEA-EA) as a baseline to inform policy development by countries. It represents a huge milestone in the history of accounting as \u0026ldquo;it moves beyond GDP and takes better account of biodiversity and ecosystems in national economic planning\u0026rdquo; (Edens et al., \u003cspan citationid=\"CR24\" class=\"CitationRef\"\u003e2022\u003c/span\u003e). Thus, it is susceptible to shaping the future of the national statistic apparatus toward the inclusion of natural capital (Edens et al., \u003cspan citationid=\"CR24\" class=\"CitationRef\"\u003e2022\u003c/span\u003e; Lange et al., \u003cspan citationid=\"CR56\" class=\"CitationRef\"\u003e2022\u003c/span\u003e). Among the challenges still to be overcome is the substantial wide variety and amount of data required to produce comprehensive accounts (Bordt, \u003cspan citationid=\"CR8\" class=\"CitationRef\"\u003e2018\u003c/span\u003e; Hein et al., \u003cspan citationid=\"CR38\" class=\"CitationRef\"\u003e2020\u003c/span\u003e). Besides, extending locally generated ES models to other locations or scales has constraints due to the necessity for parameterization, calibration, and validation, which is sometimes hampered by a lack of ground truth data (Cord et al., \u003cspan citationid=\"CR16\" class=\"CitationRef\"\u003e2017\u003c/span\u003e; Gosal et al., \u003cspan citationid=\"CR36\" class=\"CitationRef\"\u003e2022\u003c/span\u003e).\u003c/p\u003e \u003cp\u003eThe opportunities associated with technological innovations include the potential for gathering, processing, analyzing and visualizing different types of data, and integrating socio-economic and environmental information (Fleming et al., \u003cspan citationid=\"CR31\" class=\"CitationRef\"\u003e2022\u003c/span\u003e). Advancements in geospatial techniques, such as satellite imagery, geospatial mapping, drones, and sensors, are viewed as highly promising to improve decision-making on both individual and cross-sectoral levels facilitating the deployment of EA (Farrell et al., \u003cspan citationid=\"CR29\" class=\"CitationRef\"\u003e2021\u003c/span\u003e). Besides, the increased computational power and the transition towards open-access data and open-source technologies are significantly contributing to enhancing natural capital accounting (del R\u0026iacute;o-Mena et al., \u003cspan citationid=\"CR20\" class=\"CitationRef\"\u003e2023\u003c/span\u003e).\u003c/p\u003e \u003cp\u003eSatellite Earth Observation (SEO) data and technologies, coupled with Geographic Information Systems (GIS), provide advantages such as synoptic and repeated coverage, historical spatial data analysis, and cost-effectiveness (Ramirez-Reyes et al., \u003cspan citationid=\"CR68\" class=\"CitationRef\"\u003e2019\u003c/span\u003e). Therefore, where data is unavailable, imputation data techniques through statistics and Machine Learning (ML) are employed to fill in missing or incomplete data (Fleming et al., \u003cspan citationid=\"CR31\" class=\"CitationRef\"\u003e2022\u003c/span\u003e). Additionally, data fusion is commonly applied, integrating information from diverse sources such as monitoring, statistics, modelling, or interviews (Braun et al., \u003cspan citationid=\"CR9\" class=\"CitationRef\"\u003e2018\u003c/span\u003e). These approaches enable more frequent, and consistent assessment of ecosystem dynamics, regardless spatial scale of analysis (Hossain \u0026amp; Hashim, \u003cspan citationid=\"CR41\" class=\"CitationRef\"\u003e2019\u003c/span\u003e).\u003c/p\u003e \u003cp\u003eIn developing models for quantifying natural stocks, accurate and timely data on the extent and status of ecosystems are critical and challenging to gather particularly due to the dynamic and intricate interactions between water/groundwater and vegetation (Ellison et al., \u003cspan citationid=\"CR25\" class=\"CitationRef\"\u003e2017\u003c/span\u003e). Significant efforts have been made to unveil the potential of SEO in ES assessments (Willcock et al., \u003cspan citationid=\"CR78\" class=\"CitationRef\"\u003e2018\u003c/span\u003e), such as the European Copernicus programme powered by the European Spatial Agency (\u003cspan citationid=\"CR28\" class=\"CitationRef\"\u003e2023\u003c/span\u003e). The platform yields free and open access to SEO imagery and datasets, providing insights into the status of natural resources not only for Europe but also on a global scale (Almeida \u0026amp; Cabral, \u003cspan citationid=\"CR3\" class=\"CitationRef\"\u003e2023\u003c/span\u003e). Among the various products, it provides the Water \u0026amp; Wetness Probability Index (WWPI) and the Tree Cover Density (TCD) to assist in building policy-relevant ecological accounts. The WWPI is a biophysical indicator for freshwater and wetland ecosystems (Ludwig et al., \u003cspan citationid=\"CR58\" class=\"CitationRef\"\u003e2019\u003c/span\u003e) utilized as a spatial proxy to support environmental mitigation, wetland protection, erosion control, flood monitoring, and streamflow regulation (Vargas et al., \u003cspan citationid=\"CR76\" class=\"CitationRef\"\u003e2019\u003c/span\u003e). The dataset provides information on the status of BE reflecting their overall quality and conditions (Copernicus Programme, \u003cspan citationid=\"CR15\" class=\"CitationRef\"\u003e2023\u003c/span\u003e). The TCD is an indicator describing forest ecosystem conditions through a continuous spectrum of crown cover information, population numbers, densities, and wood stocks (Crowther et al., \u003cspan citationid=\"CR18\" class=\"CitationRef\"\u003e2015\u003c/span\u003e). It is acknowledged as an effective method for environmental analyses, forest management, supporting decision-making, and tracking changes in tree cover losses and gains (Chen et al., \u003cspan citationid=\"CR12\" class=\"CitationRef\"\u003e2020\u003c/span\u003e). Therefore, these products have some limitations due to: a) the availability and quality of SEO data, b) the availability and quality of in-situ data, and c) the influence of climatic and topographic conditions (European Environment Agency, \u003cspan citationid=\"CR27\" class=\"CitationRef\"\u003e2023\u003c/span\u003e). Satellite mission launch dates restrict SEO data availability, and the quality of imagery is influenced by the atmospheric effects on the light reflected off the Earth's surface and captured by the sensor. Another challenge is the lack of ground-truth data, which is expensive to acquire (Mairota et al., \u003cspan citationid=\"CR61\" class=\"CitationRef\"\u003e2015\u003c/span\u003e).\u003c/p\u003e \u003cp\u003eSome studies have demonstrated improvements with the integration of climate and topographic variables when assessing GBE. Crowther et al. (\u003cspan citationid=\"CR18\" class=\"CitationRef\"\u003e2015\u003c/span\u003e) demonstrated the strong associations between climatic variables and forest density in their research. Ludwig et al. (\u003cspan citationid=\"CR58\" class=\"CitationRef\"\u003e2019\u003c/span\u003e) combined multi-temporal optical imagery, topographic data, and spectral indices to build an automated wetland mapping detection, based on the WWPI. Madrigal-Gonz\u0026aacute;lez et al. (\u003cspan citationid=\"CR60\" class=\"CitationRef\"\u003e2023\u003c/span\u003e) used climatic and topographic characteristics to investigate the relationship between tree density and water availability. Han et al (2021) estimated carbon stock in forests by considering temperature, precipitation, elevation, and soil type as factors. Alqadhi et al (2022) demonstrated the impact of landscape and topographic characteristics on ES capacity provisioning, considering: elevation, slope, aspect, drainage density, and precipitation. However, none of these studies used Reverse Engineering (RE) to reproduce TCD and WWPI taking into consideration climatic and topographic factors to estimate inland GBE at the national level. RE refers to the process of extracting knowledge or design information from products and subsequently reproducing or recreating them based on the acquired information (Nieves-Chinchilla et al., \u003cspan citationid=\"CR66\" class=\"CitationRef\"\u003e2018\u003c/span\u003e). Through RE, it is possible to learn about the products, search for inconsistencies, or overall vulnerabilities, and assess whether there is a more efficient way to improve them (Wood, \u003cspan citationid=\"CR79\" class=\"CitationRef\"\u003e2009\u003c/span\u003e).\u003c/p\u003e \u003cp\u003eIn this study, the WWPI and TCD will be modelled within Portugal's mainland in ArcGIS Pro (ESRI, \u003cspan citationid=\"CR26\" class=\"CitationRef\"\u003e2023\u003c/span\u003e) using the \u003cem\u003eRandom Trees Regression\u003c/em\u003e algorithm, and spatial associations will be analyzed through geographical detectors utilizing the geodetectors package (version: 1.0\u0026ndash;4) in R software (Wang et al., \u003cspan citationid=\"CR77\" class=\"CitationRef\"\u003e2016\u003c/span\u003e). EA baselines (Kienast et al., \u003cspan citationid=\"CR51\" class=\"CitationRef\"\u003e2009\u003c/span\u003e) will be followed to assess the climate-driven stock accounts of forests, water, and wetlands at the national level. The aim is to develop more integrated and updated versions of TCD and WWPI contributing to the broader scientific community by exploring the potential of SEO data and technologies combined with global climatic and topographic open-access data. Furthermore, this research will shed light on the interplay between aquatic and terrestrial ecosystems and foster innovative solutions for assessing the status and conditions of inland GBE. Advancing our understanding of the development of Copernicus products will deepen our knowledge of the replicability and applicability of such datasets in EA.\u003c/p\u003e"},{"header":"2. Materials and Methods","content":"\u003cdiv id=\"Sec3\" class=\"Section2\"\u003e \u003ch2\u003e2.1. Study Area\u003c/h2\u003e \u003cp\u003ePortugal is located on the Iberian Peninsula in southwestern Europe, comprising a continental portion and two autonomous regions - the Azores and Madeira archipelagos. Portugal's mainland spans an approximate area of 90000 square km located between 36\u0026deg;57' \u0026minus;\u0026thinsp;42\u0026deg;9'N latitude and 6\u0026deg;12' \u0026minus;\u0026thinsp;9\u0026deg;30'W longitude in a climatic Mediterranean region (Beck et al., \u003cspan citationid=\"CR6\" class=\"CitationRef\"\u003e2018\u003c/span\u003e) (Fig.\u0026nbsp;\u003cspan refid=\"Fig1\" class=\"InternalRef\"\u003e1\u003c/span\u003e).\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cp\u003eThe environment changes importantly from north to south, with complex and diverse landscapes (Almeida \u0026amp; Cabral, \u003cspan citationid=\"CR2\" class=\"CitationRef\"\u003e2021\u003c/span\u003e). Spatial and temporal variability is characterized by latitude, orography, and the effect of the Atlantic Ocean, which affects the water cycle primarily through oscillations in precipitation and temperature (Belo-Pereira et al., \u003cspan citationid=\"CR7\" class=\"CitationRef\"\u003e2011\u003c/span\u003e). The territory is prone to climate variability, including droughts and desertification, like other southern European regions, notably in the southern sector, which receives the least precipitation and has the greatest temperatures (Oliveira et al., \u003cspan citationid=\"CR67\" class=\"CitationRef\"\u003e2007\u003c/span\u003e).\u003c/p\u003e \u003cp\u003eThe climatic conditions, influencing vegetation patterns, contribute to a diverse array of Atlantic, European, and Mediterranean plant species within the country (Santos et al., \u003cspan citationid=\"CR71\" class=\"CitationRef\"\u003e2023\u003c/span\u003e). Forest ecosystems predominantly consist of broad-leaved and coniferous forests, including Cork oak, Holm oak, Eucalyptus, Maritime pine, and Stone pine (Costa et al., \u003cspan citationid=\"CR17\" class=\"CitationRef\"\u003e2022\u003c/span\u003e). This ecological diversity is evident across the country, with the northern regions characterized by flourishing broadleaf and coniferous forests like Eucalyptus, Maritime pine, and Stone pine. In contrast, the southern landscape exhibits arid-adapted Mediterranean species such as cork oak, holm oak, and olive trees, representing distinct climatic zones in the country (Fonseca et al., \u003cspan citationid=\"CR32\" class=\"CitationRef\"\u003e2019\u003c/span\u003e).\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec4\" class=\"Section2\"\u003e \u003ch2\u003e2.2. Data\u003c/h2\u003e \u003cp\u003eTable\u0026nbsp;\u003cspan refid=\"Tab1\" class=\"InternalRef\"\u003e1\u003c/span\u003e presents the data used in this study, describing their source, spatial resolution, and the role played by each dataset in the classification and regression tasks.\u003c/p\u003e \u003cp\u003e \u003cdiv class=\"gridtable\"\u003e\u003ctable float=\"Yes\" id=\"Tab1\" border=\"1\"\u003e \u003ccaption language=\"En\"\u003e \u003cdiv class=\"CaptionNumber\"\u003eTable 1\u003c/div\u003e \u003cdiv class=\"CaptionContent\"\u003e \u003cp\u003eDescription of datasets utilized in the research. Legend: C \u0026ndash; Classification; R \u0026ndash; Regression; GEE \u0026ndash; Google Earth Engine.\u003c/p\u003e \u003c/div\u003e \u003c/caption\u003e \u003ccolgroup cols=\"5\"\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c1\" colnum=\"1\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c2\" colnum=\"2\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c3\" colnum=\"3\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c4\" colnum=\"4\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c5\" colnum=\"5\"\u003e\u003c/div\u003e \u003cthead\u003e \u003ctr\u003e \u003cth align=\"left\" colname=\"c1\"\u003e \u003cp\u003eData\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c2\"\u003e \u003cp\u003eSource\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c3\"\u003e \u003cp\u003ePixel resolution\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c4\"\u003e \u003cp\u003eTask\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c5\"\u003e \u003cp\u003eRole\u003c/p\u003e \u003c/th\u003e \u003c/tr\u003e \u003c/thead\u003e \u003ctbody\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eSatellite imagery 2018 - Single bands\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eEuropean Space Agency (ESA) \u0026ndash; GEE\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e10m\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003eC\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003eFeature\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eSatellite imagery 2018 - Spectral indices\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eEuropean Space Agency (ESA) \u0026ndash; GEE\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e10m\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003eC\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003eFeature\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003e\u0026zwj;Global Digital Surface Model (DSM) (m)\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eEarth Observation Research Center Japan Aerospace Exploration Agency (JAXA EORC) - GEE\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e30m\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003eC\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003eFeature\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eLand Use and Land Cover Map (COS2018) 2018 (m)\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eDirectorate General for Territory - DGT (\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://snig.dgterritorio.gov.pt/\u003c/span\u003e\u003cspan address=\"https://snig.dgterritorio.gov.pt/\" targettype=\"URL\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e)\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e10m\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003eC\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003eTarget\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eDigital Elevation Model (DEM) (m) \u0026amp; Slope (SLOPE) (%)\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eNASA Shuttle Radar Topography Mission (\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://www.earthdata.nasa.gov/\u003c/span\u003e\u003cspan address=\"https://www.earthdata.nasa.gov/\" targettype=\"URL\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e)\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e30m\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003eR\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003ePredictor\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eAverage Annual Precipitation 2018 (PREC) (mm)\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eWorldClim (\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://worldclim.org/\u003c/span\u003e\u003cspan address=\"https://worldclim.org/\" targettype=\"URL\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e)\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e1km\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003eR\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003ePredictor\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eGlobal Potential Evapotranspiration 2018 (PET) (mm)\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eConsortium of Spatial Information - CSI (\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttp://www.cgiar-csi.org\u003c/span\u003e\u003cspan address=\"http://www.cgiar-csi.org\" targettype=\"URL\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e)\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e1km\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003eR\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003ePredictor\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eGlobal Aridity Index 2018 (AI)\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eConsortium of Spatial Information -CSI (\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttp://www.cgiar-csi.org\u003c/span\u003e\u003cspan address=\"http://www.cgiar-csi.org\" targettype=\"URL\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e)\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e1km\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003eR\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003ePredictor\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eWater \u0026amp; Wetneess Probability Index (WWPI %) 2018\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eCopernicus Land Monitoring Services - CLMS (\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://land.copernicus.eu/en/products\u003c/span\u003e\u003cspan address=\"https://land.copernicus.eu/en/products\" targettype=\"URL\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e)\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e10m\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003eR\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003eResponse variable\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eTree Cover Density (TCD %) 2018\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eCopernicus Land Monitoring Services - CLMS (\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003ehttps://land.copernicus.eu/en/products\u003c/span\u003e\u003cspan address=\"https://land.copernicus.eu/en/products\" targettype=\"URL\" class=\"RefTarget\"\u003e\u003c/span\u003e\u003c/span\u003e)\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e10m\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003eR\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003eResponse variable\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003c/tbody\u003e \u003c/colgroup\u003e \u003c/table\u003e\u003c/div\u003e \u003c/p\u003e \u003cp\u003eIn the GEE platform, cloud-free composites of Sentinel-2 (S2) imagery encompassing the entire mainland of Portugal for 2018 were selected. All available spectral bands from S2 imagery, satellite-derived indices, and topographic features, were employed as classification predictors. These included the Normalized Difference Vegetation Index (NDVI)(Rouse et al., \u003cspan citationid=\"CR70\" class=\"CitationRef\"\u003e1973\u003c/span\u003e) and the Enhanced Vegetation Index (EVI)(Huete et al., \u003cspan citationid=\"CR44\" class=\"CitationRef\"\u003e1997\u003c/span\u003e) used to highlight and distinguish vegetation of different types and stages throughout their seasonal cycles; the Soil Adjusted Vegetation Index (SAVI)(Huete, \u003cspan citationid=\"CR43\" class=\"CitationRef\"\u003e1988\u003c/span\u003e) to mitigate the effects of soil reflectance and improve vegetation response, and the NDWI (Normalized Difference Water Index) used to differentiate between water and soil, based on the optical characteristics of water ​(Ji et al., \u003cspan citationid=\"CR47\" class=\"CitationRef\"\u003e2015\u003c/span\u003e). The S2-derived indices were calculated using the band arithmetic functions: (1), (2), (3), and (4):\u003cdiv id=\"Equ1\" class=\"Equation\"\u003e\u003cdiv format=\"TEX\" class=\"mathdisplay\" id=\"FileID_Equ1\" name=\"EquationSource\"\u003e\n$$\\:NDVI=\\frac{B8-B4}{B8+B4}$$\u003c/div\u003e\u003cdiv class=\"EquationNumber\"\u003e1\u003c/div\u003e\u003c/div\u003e\u003cdiv id=\"Equ2\" class=\"Equation\"\u003e\u003cdiv format=\"TEX\" class=\"mathdisplay\" id=\"FileID_Equ2\" name=\"EquationSource\"\u003e\n$$\\:EVI=2.5*\\frac{B8-B4}{B8+6*B4-7.5*B2+1}$$\u003c/div\u003e\u003cdiv class=\"EquationNumber\"\u003e2\u003c/div\u003e\u003c/div\u003e\u003cdiv id=\"Equ3\" class=\"Equation\"\u003e\u003cdiv format=\"TEX\" class=\"mathdisplay\" id=\"FileID_Equ3\" name=\"EquationSource\"\u003e\n$$\\:SAVI=\\frac{B8-B4}{(B8+B4+L)*\\left(1.0*L\\right)}$$\u003c/div\u003e\u003cdiv class=\"EquationNumber\"\u003e3\u003c/div\u003e\u003c/div\u003e\u003cdiv id=\"Equ4\" class=\"Equation\"\u003e\u003cdiv format=\"TEX\" class=\"mathdisplay\" id=\"FileID_Equ4\" name=\"EquationSource\"\u003e\n$$\\:NDWI=\\:\\frac{B3-B8}{B3+B8}$$\u003c/div\u003e\u003cdiv class=\"EquationNumber\"\u003e4\u003c/div\u003e\u003c/div\u003e\u003c/p\u003e \u003cp\u003eIn the formulas, B refers to the spectral band, and L is the soil brightness correction that ranges between [0, 1]. In the SAVI calculation, L takes the value of 0.428, as suggested by Huete (\u003cspan citationid=\"CR43\" class=\"CitationRef\"\u003e1988\u003c/span\u003e). The spectral indices (1), (2) and (3) were used to classify forests, and (4) was employed in the freshwater and wetlands classification.\u003c/p\u003e \u003cp\u003eTopographic variables, such as the Digital Surface Model (DSM) and Slope, were provided by the Japan Aerospace Exploration Agency (JAXA, \u003cspan citationid=\"CR46\" class=\"CitationRef\"\u003e2023\u003c/span\u003e). These datasets were added with the bands and spectral indices as an 18-band data stack to perform the classification tasks. The classification training datasets included reference vectorial data representing GBE (forest, freshwater, and wetlands) were extracted from the Portuguese land use and land cover map for the year 2018 (COS2018), a national product developed by the Directorate General for Territory (DGT, \u003cspan citationid=\"CR22\" class=\"CitationRef\"\u003e2018\u003c/span\u003e).\u003c/p\u003e \u003cp\u003eClimatic variables, such as Precipitation (PREC), Potential Evapotranspiration (PET), Aridity Index (AI), and the Digital Elevation Model (DEM) and Slope (SLOPE) were further utilized as explanatory variables/factors in the regression tasks. Annual PREC was obtained from the WorldClim database (Fick \u0026amp; Hijmans, \u003cspan citationid=\"CR30\" class=\"CitationRef\"\u003e2017\u003c/span\u003e); Global PET and the Global AI were obtained from the Consortium of Spatial Information, Global-AI and Global-PET Database (Zomer et al., \u003cspan citationid=\"CR82\" class=\"CitationRef\"\u003e2022\u003c/span\u003e). Global PET and AI were derived from the WorldClim precipitation model and have the exact spatial resolution and temporal scales.\u003c/p\u003e \u003cp\u003eThe response variables (WWPI and TCD) are from the High-Resolution Layers portfolio of the Copernicus. The WWPI is a raster that depicts the presence of water and wet areas as an index ranging from 0 to 100% (Copernicus Programme, \u003cspan citationid=\"CR15\" class=\"CitationRef\"\u003e2023\u003c/span\u003e). The WWPI workflow includes satellite data fusion, supervised ML classification and an unsupervised thresholding technique for water/wetness detection. Data fusion is based on optical data from S2 and Landsat 5/7/8 and radar data from Sentinel-1 multi-temporal imagery from 2012\u0026ndash;2018. The classification task was performed with the RF algorithm, and post-processing included quality assurance, control procedures, visual improvements and error correction. The TCD is a vertical projection of tree crowns onto a horizontal Earth\u0026rsquo;s surface providing pixel-level information on tree cover density (%). In the TCD development, Sentinel-2A\u0026thinsp;+\u0026thinsp;B time series (Level-2A data) were employed to build a supervised classification model using an RF classifier to create a binary Tree Cover Mask (TCM) (non-tree covered areas/tree cover). Following that, a multiple linear regression was carried out inside the bounds of the TCM to estimate TCD values.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec5\" class=\"Section2\"\u003e \u003ch2\u003e2.3. Methods\u003c/h2\u003e \u003cp\u003eForward Engineering goes from concept to product, while RE does the opposite, applying techniques to create a blueprint of an existing product (Vallero, \u003cspan citationid=\"CR75\" class=\"CitationRef\"\u003e2014\u003c/span\u003e). RE is frequently used to replicate existing products, test, or improve functionalities, and introduce new data (Nieves-Chinchilla et al., \u003cspan citationid=\"CR66\" class=\"CitationRef\"\u003e2018\u003c/span\u003e). It is defined as the process of examining an object to determine its components and their interrelationships for further redesign and recreation (Wood, \u003cspan citationid=\"CR79\" class=\"CitationRef\"\u003e2009\u003c/span\u003e). There are a few steps to follow when applying RE, such as defining the purpose, reviewing product specifications, identifying limitations and vulnerabilities, disassembling, analysing, redesigning, and reporting (Chikofsky \u0026amp; Cross II, 2002).\u003c/p\u003e \u003cp\u003eThe research development starts applying the RE approach to understand how TCD and WWPI were created, how they work, and how they could be improved, followed by employing the EA framework to quantify the physical measurements of inland GBE, through estimations of ecosystems extent (Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003ea) and condition (Fig.\u0026nbsp;\u003cspan refid=\"Fig2\" class=\"InternalRef\"\u003e2\u003c/span\u003eb). Ecosystem extent quantify the area per ecosystem type inside an accounting area. In contrast, ecosystem condition record the status of an ecosystem asset through indicators that reflect its condition (United Nations et al., \u003cspan citationid=\"CR74\" class=\"CitationRef\"\u003e2021\u003c/span\u003e). The extent was estimated through classification tasks mapping GE and BE, while the condition of forest, freshwater and wetlands was predicted by regression tasks. Finally, there was the calculation of factor importance and interaction, and statistic tests between factors and within ecosystems, employing geodetectors.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cdiv id=\"Sec6\" class=\"Section3\"\u003e \u003ch2\u003e2.3.1. Classification Task\u003c/h2\u003e \u003cp\u003eTwo distinct classification tasks were conducted: one dedicated to creating the Forest Cover Mask (FCM) and another to creating the Freshwater and Wetlands Mask (FWM). The training data were extracted from the COS18 dataset, and used as samples for a Random Forest (RF) classification (Breiman, \u003cspan citationid=\"CR10\" class=\"CitationRef\"\u003e2001\u003c/span\u003e), in Google Earth Engine (GEE) (Google, \u003cspan citationid=\"CR34\" class=\"CitationRef\"\u003e2023\u003c/span\u003e). The scripts are provided in Tables S1 and S2 in the Supplementary Materials.\u003c/p\u003e \u003cp\u003eRandom point sampling was employed to provide training and validation data for the tasks. The RF ML algorithm evaluated the training features in terms of satellite-derived data, DSM, and slope. Additionally, to optimize the model and improve its accuracy, an analysis of feature importance and hyperparameter tuning was conducted for each classification. These analyses enabled the identification of the most relevant features and refining the RF algorithm, enhancing classification accuracy for each designated class.\u003c/p\u003e \u003cp\u003eThe classification result was then subject to an accuracy assessment procedure using the validation points by calculating different metrics, such as overall accuracy (OA), confusion matrix, producers\u0026rsquo; accuracy (PA), users\u0026rsquo; accuracy (UA), kappa coefficient, and F1 score. Lastly, post-processing techniques were applied to eliminate isolated pixels and noise.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec7\" class=\"Section3\"\u003e \u003ch2\u003e2.3.2. Regression Task\u003c/h2\u003e \u003cp\u003ePredictive spatial ML regression was used to model the relationship between explanatory variables and response variables (WWPI and TCD) using ArcGIS Pro 2.9.0 software. The regression tasks were implemented using the tool \u003cem\u003eTrain Random Trees Regression Model\u003c/em\u003e.\u003c/p\u003e \u003cp\u003eThe ML spatial regression models for TCD and WWPI were built considering five explanatory variables (PREC, PET, AI, DEM, and SLOPE). The \u003cem\u003ePercent Samples for Testing\u003c/em\u003e option was set to 20%, which means that one-fifth of the training sample, referred to as test location points, was used to quantify the error for interpolation in space. Three types of errors were measured: errors on training points, errors on test points, and errors on test location points. The \u003cem\u003emaximum number of trees\u003c/em\u003e, \u003cem\u003emaximum tree depth\u003c/em\u003e, and \u003cem\u003emaximum number of samples\u003c/em\u003e were respectively set to 50, 30, and 100000. The maximum depth of each tree refers to the number of rules that each tree can generate to decide.\u003c/p\u003e \u003cp\u003eIt is critical to evaluate regression performance to understand how well the model is fitted and explained by independent variables (Hosmer \u0026amp; Lemeshow, \u003cspan citationid=\"CR40\" class=\"CitationRef\"\u003e2000\u003c/span\u003e). The metrics utilized to detect bias and the amount of variance of the response variables were the coefficient of determination (R-squared) and regression error (Re). The tool outputs a table describing the importance of each predictor used in the model, as well as scatterplots of training data, test data, and test location data, and a regression definition file containing attribute information, statistics, and model performance.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec8\" class=\"Section3\"\u003e \u003ch2\u003e2.3.3. Factor Analysis using Geodetectors\u003c/h2\u003e \u003cp\u003eFactor analysis was carried out to estimate the drivers of GBE using geographical detectors. Geodetector is a tool for investigating spatial stratified heterogeneity (SSH) through (1) measuring SSH of a variable Y; (2) testing the power of a determinant X of a dependent variable Y based on the consistency of their spatial distributions; and (3) investigating the interaction between two explanatory variables X1 and X2 (Wang et al., \u003cspan citationid=\"CR77\" class=\"CitationRef\"\u003e2016\u003c/span\u003e). The strata of Y are a partition of Y by a categorical explanatory variable X (h(X)) (J. Yang et al., \u003cspan citationid=\"CR80\" class=\"CitationRef\"\u003e2022\u003c/span\u003e). In this study, the strata refers to each ecosystem, hence the q-statistic was calculated for GE and BE. The q-statistic value is between [0, 1], then q\u0026thinsp;=\u0026thinsp;0 indicates no relationship between Y and X; q\u0026thinsp;=\u0026thinsp;1 shows that Y is governed by X. The tool is implemented by the geographical detector \u003cem\u003eq\u003c/em\u003e-statistic (5):\u003cdiv id=\"Equ5\" class=\"Equation\"\u003e\u003cdiv format=\"TEX\" class=\"mathdisplay\" id=\"FileID_Equ5\" name=\"EquationSource\"\u003e\n$$\\:q=1-\\left(\\frac{1}{N{\\sigma\\:}^{2}}\\right){\\sum\\:}_{h=1}^{a}{N}_{h}{\\sigma\\:}_{h}^{2}$$\u003c/div\u003e\u003cdiv class=\"EquationNumber\"\u003e5\u003c/div\u003e\u003c/div\u003e\u003c/p\u003e \u003cp\u003ewhere N and \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:{\\sigma\\:}^{2}\\)\u003c/span\u003e\u003c/span\u003e represent the number of units and Y variance, respectively; the population Y is made up of \u003cem\u003ea\u003c/em\u003e strata (h\u0026thinsp;=\u0026thinsp;1, 2,..., \u003cem\u003ea\u003c/em\u003e), and \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:{N}_{h}\\)\u003c/span\u003e\u003c/span\u003e and \u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:{\\sigma\\:}_{h}^{2}\\)\u003c/span\u003e\u003c/span\u003e represents the number of units and Y variance in stratum h, respectively.\u003c/p\u003e \u003cp\u003eThe factor analyses were calculated using the R geodetector package (Wang et al., \u003cspan citationid=\"CR77\" class=\"CitationRef\"\u003e2016\u003c/span\u003e). The function returns the q statistic as well as the respective p-value. The ecological detector determines whether there are statistically significant differences between two risk factors \u003cem\u003eX\u003c/em\u003e\u003csub\u003e\u003cem\u003e1\u003c/em\u003e\u003c/sub\u003e\u0026thinsp;~\u0026thinsp;\u003cem\u003eX\u003c/em\u003e\u003csub\u003e\u003cem\u003e2\u003c/em\u003e\u003c/sub\u003e. The F statistic is utilized in this function to test the differences at a significant threshold of 0.05. The interaction detector determines whether the factors interact with Y. The risk detector computes the average values for each stratum of the explanatory variable (X) and shows if there are differences between strata.\u003c/p\u003e \u003c/div\u003e \u003c/div\u003e"},{"header":"3. Results","content":"\u003cdiv id=\"Sec10\" class=\"Section2\"\u003e \u003ch2\u003e3.1. Forest, freshwater and wetlands mapping\u003c/h2\u003e \u003cp\u003eA binary classification task mapping forest/non-forest areas was performed at the national scale. Also, a multilabel classification task delineating freshwater and wetlands areas was conducted. Both classification maps are required for the calculation of ecosystem extent and are employed as masks in the subsequent step, the regression tasks (Fig.\u0026nbsp;\u003cspan refid=\"Fig3\" class=\"InternalRef\"\u003e3\u003c/span\u003e).\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cp\u003eThe FCM and FWM classifications exhibit high accuracy, strong agreement, and reliability in their respective classifications. Table\u0026nbsp;\u003cspan refid=\"Tab2\" class=\"InternalRef\"\u003e2\u003c/span\u003e presents the accuracy evaluation for forest classification (FCM) and freshwater and wetlands classification (FWM).\u003c/p\u003e \u003cp\u003e \u003cdiv class=\"gridtable\"\u003e\u003ctable float=\"Yes\" id=\"Tab2\" border=\"1\"\u003e \u003ccaption language=\"En\"\u003e \u003cdiv class=\"CaptionNumber\"\u003eTable 2\u003c/div\u003e \u003cdiv class=\"CaptionContent\"\u003e \u003cp\u003eAccuracy evaluation. OA: Overall Accuracy; PA: Producer\u0026rsquo;s Accuracy; UA: User\u0026rsquo;s Accuracy.\u003c/p\u003e \u003c/div\u003e \u003c/caption\u003e \u003ccolgroup cols=\"6\"\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c1\" colnum=\"1\"\u003e\u003c/div\u003e \u003cdiv align=\"char\" char=\".\" class=\"colspec\" colname=\"c2\" colnum=\"2\"\u003e\u003c/div\u003e \u003cdiv align=\"char\" char=\".\" class=\"colspec\" colname=\"c3\" colnum=\"3\"\u003e\u003c/div\u003e \u003cdiv align=\"char\" char=\".\" class=\"colspec\" colname=\"c4\" colnum=\"4\"\u003e\u003c/div\u003e \u003cdiv align=\"char\" char=\".\" class=\"colspec\" colname=\"c5\" colnum=\"5\"\u003e\u003c/div\u003e \u003cdiv align=\"char\" char=\".\" class=\"colspec\" colname=\"c6\" colnum=\"6\"\u003e\u003c/div\u003e \u003cthead\u003e \u003ctr\u003e \u003cth align=\"left\" colname=\"c1\"\u003e \u003cp\u003eLULC\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c2\"\u003e \u003cp\u003eOA (%)\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c3\"\u003e \u003cp\u003ePA (%)\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c4\"\u003e \u003cp\u003eUA (%)\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c5\"\u003e \u003cp\u003eKappa\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c6\"\u003e \u003cp\u003eF1 Score\u003c/p\u003e \u003c/th\u003e \u003c/tr\u003e \u003c/thead\u003e \u003ctbody\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eGE\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c2\"\u003e \u003cp\u003e96.17\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e93.68\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e \u003cp\u003e95.71\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e \u003cp\u003e0.89\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c6\"\u003e \u003cp\u003e0.93\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eBE\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c2\"\u003e \u003cp\u003e94.66\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e95.83\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e \u003cp\u003e96.67\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c5\"\u003e \u003cp\u003e0.94\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c6\"\u003e \u003cp\u003e0.94\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003c/tbody\u003e \u003c/colgroup\u003e \u003c/table\u003e\u003c/div\u003e \u003c/p\u003e \u003cp\u003eThe FWM classification achieved an overall accuracy (OA) of 94.66%, while the FCM classification achieved an even higher accuracy of 96.17%, indicating a high level of correctness in their respective classifications. For the producer\u0026rsquo;s accuracy (PA), the FCM classification demonstrated a result of 93.68%, while the FWM classification showed a slightly higher result (95.83%), indicating their ability to identify their respective categories accurately. As for user accuracy (UA), the FCM classification exhibited an accuracy of 95.71%, while the FWM classification demonstrated an accuracy of 96.67%, indicating the reliability of both models\u0026rsquo; predictions.\u003c/p\u003e \u003cp\u003eThe forest cover classification exhibited a robust agreement, boasting a Kappa value of 0.89. The classification of freshwater and wetlands demonstrated an even higher level of agreement, with a Kappa value of 0.94. These results underscore both models\u0026rsquo; strong consistency and reliability in accurately categorizing vegetation, freshwater, and wetlands. The F1 scores further reinforce the models\u0026rsquo; effectiveness in correctly identifying these land cover categories. The FCM classification achieved a commendable F1 score of 0.93, while the FWM classification excelled with an impressive F1 score of 0.94.\u003c/p\u003e \u003cp\u003eWhen looking at the mean percentage weight of feature importance for the water and wetlands classification at the national level (Fig.\u0026nbsp;\u003cspan refid=\"Fig4\" class=\"InternalRef\"\u003e4\u003c/span\u003e), the analysis provided insights into the features crucial for discriminating water bodies and wetland areas accurately. As such, among the 13 attributes assessed, the surface elevation (DSM) exhibited the highest significance (13.65%), followed by NDWI (9.70%) and Slope (9.30%). S2 Bands (8, 11, 12, 8A, 4, 6, 5 and 7) presented importance higher than 5%, while Bands (3 and 2) were less than 5%.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cp\u003eRegarding the forest cover classification, a total of 15 attributes were considered as predictors. The feature analysis demonstrates that the spectral indices, SAVI (9.14%), NDVI (8.19%) and EVI (6.89%), showcased notable significance, but also S2 Bands 11 (8.70%), 3 (8.09%) and 12 (7.91%). The Short Wave Infrared (SWIR) bands, B11 and B12, capture variations in vegetation water content and the spongy mesophyll structure within vegetation canopies. Meanwhile, B3 is sensitive to healthy vegetation. The S2 Band 8A was the only feature contributing with less than 5%. The S2 Bands (7, 8, 6, 5, 2 and 4), DSM and Slope presented an importance weight of more than 5% but less than 6.89%.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec11\" class=\"Section2\"\u003e \u003ch2\u003e3.2. FCD and FWPI predictions\u003c/h2\u003e \u003cp\u003eR-squared and residuals were used to assess the accuracy of FCD and FWPI predictions. The true values were compared to predicted values in terms of proportion of variance (R-squared), and residuals over targets. Figure\u0026nbsp;\u003cspan refid=\"Fig5\" class=\"InternalRef\"\u003e5\u003c/span\u003e shows the scatterplot of predicted (FCD) over-reference values (TCD), and results for R-squared at training points (Tr), at test points (Te), and test location points (Tl), as well as the residual plot. The R-squared at Tr, Te and Tl, were respectively, 0.8746, 0.5456, and 0.3208. The R-squared for residuals is 0.4131 and the mean residual value is 0.04. These statistics were calculated using 100000 samples.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cp\u003eThe scatterplot of predictions (FWPI) over observed values (WWPI) and the residual plot for the FWPI model are shown in Fig.\u0026nbsp;\u003cspan refid=\"Fig6\" class=\"InternalRef\"\u003e6\u003c/span\u003e. R-squared values at Tr, Te, and Tl were 0.8275, 0.6153, and 0.6555, respectively. The R-squared residual is 0.4765 and the mean residual value is 0.08. In the predictions of FWPI, the R-squared showcased superior results in Te and Tl. However, FCD performed better than FWPI, as R-squared at Tr is higher, and corresponds to 80% of samples.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cp\u003eTable\u0026nbsp;\u003cspan refid=\"Tab3\" class=\"InternalRef\"\u003e3\u003c/span\u003e presents the equations for each sample set (Tr, Te, and Tl) used to assess the models\u0026rsquo; performance.\u003c/p\u003e \u003cp\u003e \u003cdiv class=\"gridtable\"\u003e\u003ctable float=\"Yes\" id=\"Tab3\" border=\"1\"\u003e \u003ccaption language=\"En\"\u003e \u003cdiv class=\"CaptionNumber\"\u003eTable 3\u003c/div\u003e \u003cdiv class=\"CaptionContent\"\u003e \u003cp\u003ePerformance evaluation by model (Tr: Training data at training locations, Te: Test data at training locations, and Tl: Test data at test locations).\u003c/p\u003e \u003c/div\u003e \u003c/caption\u003e \u003ccolgroup cols=\"5\"\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c1\" colnum=\"1\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c2\" colnum=\"2\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c3\" colnum=\"3\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c4\" colnum=\"4\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c5\" colnum=\"5\"\u003e\u003c/div\u003e \u003cthead\u003e \u003ctr\u003e \u003cth align=\"left\" colname=\"c1\"\u003e \u003cp\u003eModel\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c2\"\u003e \u003cp\u003eResiduals\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c3\"\u003e \u003cp\u003eTr\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c4\"\u003e \u003cp\u003eTe\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c5\"\u003e \u003cp\u003eTl\u003c/p\u003e \u003c/th\u003e \u003c/tr\u003e \u003c/thead\u003e \u003ctbody\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eFCD\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:y=0.32x-15.76\\)\u003c/span\u003e\u003c/span\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:y=0.74x+13.68\\)\u003c/span\u003e\u003c/span\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:y=0.54x+22.83\\)\u003c/span\u003e\u003c/span\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003e\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:y=0.32x+24.62\\)\u003c/span\u003e\u003c/span\u003e\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eFWPI\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:y=0.36x-\\:19.82\\)\u003c/span\u003e\u003c/span\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:y=0.71x+17.9\\)\u003c/span\u003e\u003c/span\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:y=0.55x+24.71\\)\u003c/span\u003e\u003c/span\u003e\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003e\u003cspan class=\"InlineEquation\"\u003e\u003cspan class=\"mathinline\"\u003e\\(\\:y=0.6x+24.21\\)\u003c/span\u003e\u003c/span\u003e\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003c/tbody\u003e \u003c/colgroup\u003e \u003c/table\u003e\u003c/div\u003e \u003c/p\u003e \u003cp\u003eThe explanatory variables for both regression models were DEM, SLOPE, PREC, PET, and AI. Figure\u0026nbsp;\u003cspan refid=\"Fig7\" class=\"InternalRef\"\u003e7\u003c/span\u003e describes the importance of each variable for FCD and FWPI. The predictors with the highest percentage values were more relevant to the model. All variables were shown to be important for both models, except the SLOPE. DEM and PREC were the factors that most contributed to the FWPI model, followed by PET and AI, and lastly, SLOPE. For the FCD model, PET and AI weighted the most, followed by PREC, DEM, and SLOPE.\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec12\" class=\"Section2\"\u003e \u003ch2\u003e3.3. Factor analysis using geodetectors\u003c/h2\u003e \u003cp\u003eThe factor detector q-statistic indicates if explanatory variables stratify Y (GE and BE). Table\u0026nbsp;\u003cspan refid=\"Tab4\" class=\"InternalRef\"\u003e4\u003c/span\u003e depicts the factor detector results. In the prediction of BE, PREC and AI presented q-values higher than 0.5; DEM and PET higher than 0.25; and SLOPE close to 0, meaning there was no association between FWPI predictions and SLOPE. For GE, PREC and AI show the highest degree of association with values higher than 0.5, followed by PET (0.44), SLOPE (0.13) and DEM (0.11).\u003c/p\u003e \u003cp\u003e \u003cdiv class=\"gridtable\"\u003e\u003ctable float=\"Yes\" id=\"Tab4\" border=\"1\"\u003e \u003ccaption language=\"En\"\u003e \u003cdiv class=\"CaptionNumber\"\u003eTable 4\u003c/div\u003e \u003cdiv class=\"CaptionContent\"\u003e \u003cp\u003eFactor detector results for Blue Ecosystems (BE) and Green Ecosystems (GE) estimation.\u003c/p\u003e \u003c/div\u003e \u003c/caption\u003e \u003ccolgroup cols=\"7\"\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c1\" colnum=\"1\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c2\" colnum=\"2\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c3\" colnum=\"3\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c4\" colnum=\"4\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c5\" colnum=\"5\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c6\" colnum=\"6\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c7\" colnum=\"7\"\u003e\u003c/div\u003e \u003cthead\u003e \u003ctr\u003e \u003cth align=\"left\" colname=\"c1\"\u003e \u003cp\u003eY\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c2\"\u003e \u003cp\u003eStatistics\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c3\"\u003e \u003cp\u003eSLOPE\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c4\"\u003e \u003cp\u003eDEM\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c5\"\u003e \u003cp\u003ePET\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c6\"\u003e \u003cp\u003eAI\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c7\"\u003e \u003cp\u003ePREC\u003c/p\u003e \u003c/th\u003e \u003c/tr\u003e \u003c/thead\u003e \u003ctbody\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eBE\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eq-statistic\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e0.06\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e0.42\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003e0.29\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e0.55\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c7\"\u003e \u003cp\u003e0.69\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eBE\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003ep-value\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e1\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e9.04E-10\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003e9.05E-10\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e9.35E-10\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c7\"\u003e \u003cp\u003e9.19E-10\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eGE\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003eq-statistic\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e0.13\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e0.11\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003e0.44\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e0.54\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c7\"\u003e \u003cp\u003e0.65\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eGE\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003ep-value\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c3\"\u003e \u003cp\u003e1\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c4\"\u003e \u003cp\u003e9.50E-10\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c5\"\u003e \u003cp\u003e8.47E-10\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c6\"\u003e \u003cp\u003e9.66E-10\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c7\"\u003e \u003cp\u003e9.36E-10\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003c/tbody\u003e \u003c/colgroup\u003e \u003c/table\u003e\u003c/div\u003e \u003c/p\u003e \u003cp\u003eThe ecological detector determines if there are statistically significant differences between two covariates. This function utilizes the F-statistic to test for differences, employing a significance threshold of 0.05. In both GE and BE models, all variables were found to be significantly different.\u003c/p\u003e \u003cp\u003eThe interaction detector determines whether the factors interact with the covariate. All interactions between factors had q-values greater than individual factors, demonstrating that combining carbon and water cycle phenomenon-related factors increases the explanatory power of GBE (Fig.\u0026nbsp;\u003cspan refid=\"Fig8\" class=\"InternalRef\"\u003e8\u003c/span\u003e). The strongest interactions (greater than 0.75) were observed between AI and PREC (0.83), DEM and PREC (0.82), AI and PET, PET and PREC (0.78), and between DEM and PET (0.75). It was not verified that there were any interactions between AI with SLOPE, AI with DEM, and SLOPE with PREC, but SLOPE interacted with DEM (0.59) and PET (0.53).\u003c/p\u003e \u003cp\u003e \u003c/p\u003e \u003cp\u003eThe risk detector computes the average values within ecosystems and determines if there are significant differences between the means of each stratum. Considering the country's extent, the mean value of GE corresponds to 51.47%. The mean value of BE referring to wetlands was 66.36% and 87.23% for freshwater ecosystems. All statistics were significantly different, considering a significant threshold of 0.05.\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec13\" class=\"Section2\"\u003e \u003ch2\u003e3.4. Inland GBE physical accounting: extent and condition\u003c/h2\u003e \u003cp\u003eIn 2018, forests, freshwater bodies, and wetlands covered around 62% of the territory of Continental Portugal. Forests occupied a total area of 5287284.95 hectares covering 60.32% of the total country area, freshwater bodies filled an area of 128391.01 ha, constituting roughly 1.7%, and wetlands 0.33% occupying an area of 21792.51 ha.\u003c/p\u003e \u003cp\u003eThe condition of forests was assessed with the mean forest crown cover density (FCD), and the mean FWPI a proxy for the status of freshwater and wetlands. The GBE status assessment is calculated by pixel, meaning that, in an area of 100 square meters classified as forests, 51.47% of this area will account for aboveground biomass. In a pixel classified as freshwater and wetlands, respectively, 87.23% and 66.36% will effectively account for Blue Ecosystems. The area (hectares), proportional extent (%), status (%), and physical accounts (%) of forests, freshwater bodies and wetlands in 2018 at the national scale are presented in Table\u0026nbsp;\u003cspan refid=\"Tab5\" class=\"InternalRef\"\u003e5\u003c/span\u003e.\u003c/p\u003e \u003cp\u003e \u003cdiv class=\"gridtable\"\u003e\u003ctable float=\"Yes\" id=\"Tab5\" border=\"1\"\u003e \u003ccaption language=\"En\"\u003e \u003cdiv class=\"CaptionNumber\"\u003eTable 5\u003c/div\u003e \u003cdiv class=\"CaptionContent\"\u003e \u003cp\u003eAccounting for inland GBE extent and condition for Mainland Portugal in 2018.\u003c/p\u003e \u003c/div\u003e \u003c/caption\u003e \u003ccolgroup cols=\"4\"\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c1\" colnum=\"1\"\u003e\u003c/div\u003e \u003cdiv align=\"left\" class=\"colspec\" colname=\"c2\" colnum=\"2\"\u003e\u003c/div\u003e \u003cdiv align=\"char\" char=\".\" class=\"colspec\" colname=\"c3\" colnum=\"3\"\u003e\u003c/div\u003e \u003cdiv align=\"char\" char=\".\" class=\"colspec\" colname=\"c4\" colnum=\"4\"\u003e\u003c/div\u003e \u003cthead\u003e \u003ctr\u003e \u003cth align=\"left\" colname=\"c1\"\u003e \u003cp\u003e2018\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c2\"\u003e \u003cp\u003eForest\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c3\"\u003e \u003cp\u003eFreshwater\u003c/p\u003e \u003c/th\u003e \u003cth align=\"left\" colname=\"c4\"\u003e \u003cp\u003eWetlands\u003c/p\u003e \u003c/th\u003e \u003c/tr\u003e \u003c/thead\u003e \u003ctbody\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eArea (ha)\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e5287285\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e128391.01\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e \u003cp\u003e21792.51\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eProportional extent (%)\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e60.36\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e1.7\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e \u003cp\u003e0.33\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003eStatus (%)\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e51.47\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e87.23\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e \u003cp\u003e66.36\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003ctr\u003e \u003ctd align=\"left\" colname=\"c1\"\u003e \u003cp\u003ePhysical accounts (%)\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"left\" colname=\"c2\"\u003e \u003cp\u003e31.07\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c3\"\u003e \u003cp\u003e1.48\u003c/p\u003e \u003c/td\u003e \u003ctd align=\"char\" char=\".\" colname=\"c4\"\u003e \u003cp\u003e0.22\u003c/p\u003e \u003c/td\u003e \u003c/tr\u003e \u003c/tbody\u003e \u003c/colgroup\u003e \u003c/table\u003e\u003c/div\u003e \u003c/p\u003e \u003c/div\u003e"},{"header":"4. Discussion","content":"\u003cp\u003eWe proposed a novel methodology that relies on RE techniques to disassemble TCD and WWPI datasets from the Copernicus products portfolio and redesign and reproduce them introducing new data and methods to enhance usability. Besides, our study delivered a hands-on approach to estimating climate-driven stock accounts of forest, freshwater, and wetland ecosystems. The findings on the drivers of biomass density and water availability shed new light on the potential implications of land cover conversions and climate conditions on inland GBE. This information is critical, given that these ecosystems are biodiversity hotspots capable of supporting a wide range of livelihoods, and able to provide a variety of goods and services.\u003c/p\u003e \u003cdiv id=\"Sec15\" class=\"Section2\"\u003e \u003ch2\u003e4.1. Cloud-based automated classification of inland GBE\u003c/h2\u003e \u003cp\u003eOpen-access geospatial data streams and open-source software have made frequent, large-scale, high-resolution land cover mapping a reality (Arruda et al., \u003cspan citationid=\"CR4\" class=\"CitationRef\"\u003e2021\u003c/span\u003e)​. These recent breakthroughs in acquiring and processing big earth data have made large-scale, high-resolution land cover classifications more viable for a diverse range of user groups, organizations, and scholars (DeLancey et al., \u003cspan citationid=\"CR21\" class=\"CitationRef\"\u003e2019\u003c/span\u003e). GEE has evolved in recent years as an appealing high-performance computing platform for the cloud-based processing of petabytes of earth data, offering computational capability for planetary-scale data processing, as well as the ability to apply well-known ML algorithms (Gorelick et al., \u003cspan citationid=\"CR35\" class=\"CitationRef\"\u003e2017\u003c/span\u003e).\u003c/p\u003e \u003cp\u003eThe separation between GE and BE was implemented primarily to test the significant improvements in the accuracy and utility of such analyses (Koma et al., \u003cspan citationid=\"CR54\" class=\"CitationRef\"\u003e2021\u003c/span\u003e). The fundamental motivation for segregating GE and BE into distinct categories arises from their inherently disparate spectral characteristics. Vegetation, characterized by its chlorophyll content and unique reflectance properties, displays pronounced spectral features in the green and near-infrared regions (Gitelson et al., \u003cspan citationid=\"CR33\" class=\"CitationRef\"\u003e2003\u003c/span\u003e). Conversely, water bodies, including wetlands, exhibit spectral absorption patterns in the same regions (Khalid et al., \u003cspan citationid=\"CR50\" class=\"CitationRef\"\u003e2021\u003c/span\u003e). These distinctive spectral properties lay the foundation for specialized classification approaches.\u003c/p\u003e \u003cp\u003eThe classification analysis has underscored the improvement in classification accuracy achieved through this separation. The Kappa coefficient in both classifications showed compelling results. This enhancement in accuracy is not merely an academic achievement but has practical implications for informed decision-making (Jagannathan \u0026amp; Divya, \u003cspan citationid=\"CR45\" class=\"CitationRef\"\u003e2021\u003c/span\u003e). Beyond accuracy, the separation of GE and BE enables tailored analysis for specific applications. For instance, GE classification supports studies related to vegetation health, land-use planning, and forest management, while BE classification aids in monitoring water bodies, wetland conservation, and hydrological assessments (Kulawardhana, \u003cspan citationid=\"CR55\" class=\"CitationRef\"\u003e2011\u003c/span\u003e).\u003c/p\u003e \u003cp\u003eThe F1 scores, indicative of the models' robustness in correctly identifying GE and BE categories, further emphasize the advantages of this separation (Alem \u0026amp; Kumar, \u003cspan citationid=\"CR1\" class=\"CitationRef\"\u003e2022\u003c/span\u003e). GE achieved an F1 score of 0.93, while BE classification surpassed an F1 score of 0.94. These high F1 scores underscore the models' ability to effectively distinguish between these two categories.\u003c/p\u003e \u003cp\u003eA stand-alone ML classification model for mapping GE and BE has significant implications for environmental assessment, empowering researchers to conduct focused analyses on ecosystems, and supporting decision-makers in crafting sustainable land management and conservation strategies (Chignell et al., \u003cspan citationid=\"CR13\" class=\"CitationRef\"\u003e2018\u003c/span\u003e; Z. Yang et al., \u003cspan citationid=\"CR81\" class=\"CitationRef\"\u003e2021\u003c/span\u003e). Furthermore, the outputs of these classification tasks, the FCM from GE mapping and FWM from the BE were utilized to delineate ecosystem assets and quantify physical measurements of inland GBE (United Nations et al., \u003cspan citationid=\"CR74\" class=\"CitationRef\"\u003e2021\u003c/span\u003e).\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec16\" class=\"Section2\"\u003e \u003ch2\u003e4.2. From Reverse Engineering to satellite-based climate-driven models\u003c/h2\u003e \u003cp\u003eThere is an emerging trend involving using SEO data to quantify ecosystem structure and functional traits which serve as more accurate indicators of ES than simplistic cover-based proxies (Schirpke et al., 2023). Their combination with regression models contributes to a comprehensive understanding of the social-ecological processes influencing the provision of ES and yields reliable data and results (C. Ramirez-Reyes et al., \u003cspan citationid=\"CR68\" class=\"CitationRef\"\u003e2019\u003c/span\u003e) and thereby, provide cost-effective data to supply ecological accounts (Vargas et al., \u003cspan citationid=\"CR76\" class=\"CitationRef\"\u003e2019\u003c/span\u003e).\u003c/p\u003e \u003cp\u003eThe WWPI and TCD have the potential to provide data to build ecological-economic accounts for GBE. Redesigning and reproducing these datasets focused on forest, freshwater and wetlands evidenced the interdependency of some ecosystem functions, such as biomass density, water yield and soil moisture water content, and climatic and topographic drivers, as shown by Mengist et al. (\u003cspan citationid=\"CR62\" class=\"CitationRef\"\u003e2021\u003c/span\u003e). Topography and climate are key elements in explaining changes in vegetation abundance at the local (i.e. pixel, plot, or parcel scale) and regional levels (Mpakairi et al., \u003cspan citationid=\"CR63\" class=\"CitationRef\"\u003e2022\u003c/span\u003e).\u003c/p\u003e \u003cp\u003eFCD predictions were better than FWPI when comparing the results of R-squared, regression errors and residuals. The independent variables considered in both models were topographic characteristics (SLOPE and DEM) and climatic factors (PREC, PET and AI). Using the same factors to estimate SSH within carbon and water cycles demonstrated the intrinsic relationships between aquatic and terrestrial ecosystems (Cantonati et al., \u003cspan citationid=\"CR11\" class=\"CitationRef\"\u003e2020\u003c/span\u003e; Dwire et al., \u003cspan citationid=\"CR23\" class=\"CitationRef\"\u003e2018\u003c/span\u003e; Hose et al., \u003cspan citationid=\"CR39\" class=\"CitationRef\"\u003e2021\u003c/span\u003e; Huang et al., \u003cspan citationid=\"CR42\" class=\"CitationRef\"\u003e2020\u003c/span\u003e; Kl\u0026oslash;ve et al., \u003cspan citationid=\"CR52\" class=\"CitationRef\"\u003e2014\u003c/span\u003e; Koit et al., \u003cspan citationid=\"CR53\" class=\"CitationRef\"\u003e2021\u003c/span\u003e). The aboveground biomass density represented by FCD estimates varies spatially and temporally according to ecosystem types, altitude, latitude, climatic conditions, and water presence (Tabacchi et al., \u003cspan citationid=\"CR73\" class=\"CitationRef\"\u003e2011\u003c/span\u003e). GE supports biodiversity and many livelihoods around lakes, perennial streams, and wetlands, which benefit from BE\u0026rsquo;s water and nutrient availability (Harris, \u003cspan citationid=\"CR37\" class=\"CitationRef\"\u003e1984\u003c/span\u003e). Mpakairi et al. (\u003cspan citationid=\"CR63\" class=\"CitationRef\"\u003e2022\u003c/span\u003e) have shown that the density of healthy trees decreased due to limited water and nutrient availability, with water presence driving about 40% of biological processes in forests.\u003c/p\u003e \u003cp\u003eRegarding the contribution of each variable to the models, all of them showed to be relevant, a part of SLOPE presenting very low importance to both predictions. DEM and PREC were the factors that most contributed to the FWPI model, followed by PET, AI, and SLOPE. For the FCD model, PET and AI contributed the most, followed by PREC, DEM, and SLOPE. The results are aligned with those presented in the studies of Crowther et al. (\u003cspan citationid=\"CR18\" class=\"CitationRef\"\u003e2015\u003c/span\u003e) and Madrigal-Gonz\u0026aacute;lez et al. (\u003cspan citationid=\"CR60\" class=\"CitationRef\"\u003e2023\u003c/span\u003e), analysing the influence of climatic and topographic conditions in assessing GBE. These findings expanded our understanding of nature stock climate drivers, landscape behaviours, and increasing accuracy in ecosystem accounts. Therefore, unleashing the capacity to analyze the status and conditions of GBE helps to execute integrated natural resource management.\u003c/p\u003e \u003cp\u003eThe geographical detector analysis used to verify SSH enriched the knowledge of the drivers of GE and BE. The factor detectors indicate that both response variables are stratified by the exploratory variables, with more evidence of stratification shown by PREC, PET and AI. DEM showed a higher association with the BE than the GE model. Lastly, the physical accounts of inland GBE were assessed using the EA methodological framework. The national GBE extent account represented their magnitude in terms of area in hectares, while the GBE condition described their integrity and status (Cz\u0026uacute;cz et al., \u003cspan citationid=\"CR19\" class=\"CitationRef\"\u003e2021\u003c/span\u003e). The aboveground biomass in the forest was measured by calculating the area classified as forest, and the mean estimated value of FCD was applied to assess GE condition. BE were evaluated by determining the extent of freshwater bodies and wetlands and the mean estimated values of FWPI for each ecosystem. The physical accounts evidenced that in 2018, 33% of the total inland area of Portugal was represented by forest, freshwater and wetlands ecosystems. The physical accounts offer means to understand how sustainable the production is through measurements of natural stocks and yields, providing baselines for the assessment of ecosystems\u0026rsquo; capacity to deliver goods and services at multiscale (Keith et al., \u003cspan citationid=\"CR48\" class=\"CitationRef\"\u003e2020\u003c/span\u003e).\u003c/p\u003e \u003c/div\u003e \u003cdiv id=\"Sec17\" class=\"Section2\"\u003e \u003ch2\u003e4.3. Implications for future research\u003c/h2\u003e \u003cp\u003eThere are challenges with reverse engineering, especially when it comes to digital products (Chikofsky \u0026amp; Cross II, 2002). Other issues outside understanding the dataset itself include determining how it was created when we have little or no access to inputs, intermediate outputs, or even the code, algorithms, and tools used in each task (Wood, \u003cspan citationid=\"CR79\" class=\"CitationRef\"\u003e2009\u003c/span\u003e). Without clear documentation of the model architecture, hyperparameters, and training process, it may be difficult for others to reproduce methods or understand the model's behaviour, and how the model arrives at its predictions. The transparency through reporting the entire ML pipeline, including data preprocessing, model architecture, hyperparameters, and evaluation metrics enables others to reproduce and validate the research findings and improve the robustness and reliability of ML models, leading to better accuracy and more trustworthy results. Additionally, continuous monitoring and updating of models as new data becomes available can help maintain their performance over time.\u003c/p\u003e \u003cp\u003eNowadays, a wide range of technologies exists to aid in the forward/reverse engineering process, especially in the production of digital products. A good improvement would be to provide more information about these products' development and make version control available to the public. This practice helps in keeping track of dataset modifications, comparing several dataset versions, and revert to whichever previous version is required. Future research would further take these products apart to continue improving functionality, introducing new data and testing algorithms, as well as replicating to other areas in the globe. Additionally, it would be interesting to explore other products developed by Copernicus, which comprises a wide variety of products, such as Soil Water Index, Imperviousness Density and Impervious Built-up, Dominant Leaf Type and Forest Type, Coastal Zones and Protected Areas. These datasets are valuable resources for a wide range of geospatial services and environmental analysis.\u003c/p\u003e \u003c/div\u003e"},{"header":"5. Conclusion","content":"\u003cp\u003eThis research demonstrated the potential of the TCD and WWPI products provided by Copernicus as fundamental data for ecological-economic studies and presented an approach to overcome limitations and enhance their usability. Besides, it underscores the significance of climatic (PREC, PET, AT) and topographic (DEM and SLOPE) variables influencing the status and conditions of GBE, by analyzing spatial relationships in the TCD and WWPI predictions. RE proved to be a valuable tool for innovation and problem-solving while examining the product's functionality and behaviour to determine how it works and identifying any vulnerabilities, weaknesses, or areas for improvement. The outcomes demonstrate the synergy of combining global open-access data with open-source technologies in creating predictive models that offer actionable insights, thereby facilitating evidence-based decision-making across various environmental and societal challenges.\u003c/p\u003e"},{"header":"Declarations","content":"\u003ch2\u003eFunding\u003c/h2\u003e \u003cp\u003eOpen Access funding was provided by NOVA Information Management School NOVA University Lisbon. This study was funded by the research project MaSOT \u0026ndash; Mapping Ecosystem Services from Earth Observations, funded by the Portuguese Science Foundation \u0026ndash; FCT [EXPL/CTA-AMB/0165/2021], and by national funds through FCT (Funda\u0026ccedil;\u0026atilde;o para a Ci\u0026ecirc;ncia e a Tecnologia), under the project - UIDB/04152/2020 - Centro de Investiga\u0026ccedil;\u0026atilde;o em Gest\u0026atilde;o de Informa\u0026ccedil;\u0026atilde;o (MagIC)/NOVA IMS.\u003c/p\u003e\u003ch2\u003eAuthor Contribution\u003c/h2\u003e\u003cp\u003eB.A.: Conceptualization, Methodology, Investigation, Formal analysis, Writing - original draft, preparation, Writing - review and editing. L.M.: Conceptualization, Methodology, Investigation, Writing - review and editing. P.S.: Conceptualization, Methodology, Investigation, Writing - review and editing. P.C.: Conceptualization, Methodology, Investigation, Writing - review and editing.\u003c/p\u003e\u003ch2\u003eAcknowledgement\u003c/h2\u003e\u003cp\u003eThis study was supported by the research project MaSOT \u0026ndash; Mapping Ecosystem Services from Earth Observations, funded by the Portuguese Science Foundation \u0026ndash; FCT [EXPL/CTA-AMB/0165/2021], and by national funds through FCT (Funda\u0026ccedil;\u0026atilde;o para a Ci\u0026ecirc;ncia e a Tecnologia), under the project - UIDB/04152/2020 - Centro de Investiga\u0026ccedil;\u0026atilde;o em Gest\u0026atilde;o de Informa\u0026ccedil;\u0026atilde;o (MagIC)/NOVA IMS. We are grateful to the Institut fran\u0026ccedil;ais du Portugal for fostering Franco-Portuguese collaborations through scientific periods in France (BIC-2023), and the IUEM for hosting.\u003c/p\u003e"},{"header":"References","content":"\u003col\u003e\n \u003cli\u003eAlem, A., \u0026amp; Kumar, S. (2022). \u003cem\u003eApplied Artificial Intelligence End-to-End Convolutional Neural Network Feature Extraction for Remote Sensed Images Classification\u003c/em\u003e. https://doi.org/10.1080/08839514.2022.2137650\u003c/li\u003e\n \u003cli\u003eAlmeida, B., \u0026amp; Cabral, P. (2021). Water yield modelling, sensitivity analysis and validation: A study for Portugal. \u003cem\u003eISPRS International Journal of Geo-Information\u003c/em\u003e, \u003cem\u003e10\u003c/em\u003e(8). https://doi.org/10.3390/ijgi10080494\u003c/li\u003e\n \u003cli\u003eAlmeida, B., \u0026amp; Cabral, P. (2023). Data-Driven Modelling of Freshwater Ecosystems: A Multiscale Framework Based on Global Geospatial Data. \u003cem\u003eProceedings of the 9th International Conference on Geographical Information Systems Theory, Applications and Management\u003c/em\u003e, \u003cem\u003e2023-April\u003c/em\u003e, 104\u0026ndash;111. https://doi.org/10.5220/0012037800003473\u003c/li\u003e\n \u003cli\u003eArruda, V. L. S., Piontekowski, V. J., Alencar, A., Pereira, R. S., \u0026amp; Matricardi, E. A. T. (2021). An alternative approach for mapping burn scars using Landsat imagery, Google Earth Engine, and Deep Learning in the Brazilian Savanna. \u003cem\u003eRemote Sensing Applications: Society and Environment\u003c/em\u003e, \u003cem\u003e22\u003c/em\u003e, 100472. https://doi.org/10.1016/j.rsase.2021.100472\u003c/li\u003e\n \u003cli\u003eBateman, I. J., \u0026amp; Mace, G. M. (2020). The natural capital framework for sustainably efficient and equitable decision making. \u003cem\u003eNature Sustainability\u003c/em\u003e, \u003cem\u003e3\u003c/em\u003e(10), 776\u0026ndash;783. https://doi.org/10.1038/s41893-020-0552-3\u003c/li\u003e\n \u003cli\u003eBeck, H. E., Zimmermann, N. E., McVicar, T. R., Vergopolan, N., Berg, A., \u0026amp; Wood, E. F. (2018). Present and future K\u0026ouml;ppen-Geiger climate classification maps at 1-km resolution. \u003cem\u003eScientific Data\u003c/em\u003e, \u003cem\u003e5\u003c/em\u003e(1), 180214. https://doi.org/10.1038/sdata.2018.214\u003c/li\u003e\n \u003cli\u003eBelo-Pereira, M., Dutra, E., \u0026amp; Viterbo, P. (2011). Evaluation of global precipitation data sets over the Iberian Peninsula. \u003cem\u003eJournal of Geophysical Research: Atmospheres\u003c/em\u003e, \u003cem\u003e116\u003c/em\u003e(D20). https://doi.org/https://doi.org/10.1029/2010JD015481\u003c/li\u003e\n \u003cli\u003eBordt, M. (2018). Discourses in Ecosystem Accounting: A Survey of the Expert Community. \u003cem\u003eEcological Economics\u003c/em\u003e, \u003cem\u003e144\u003c/em\u003e, 82\u0026ndash;99. https://doi.org/10.1016/J.ECOLECON.2017.06.032\u003c/li\u003e\n \u003cli\u003eBraun, D., Damm, A., Hein, L., Petchey, O. L., \u0026amp; Schaepman, M. E. (2018). Spatio-temporal trends and trade-offs in ecosystem services: An Earth observation based assessment for Switzerland between 2004 and 2014. \u003cem\u003eEcological Indicators\u003c/em\u003e, \u003cem\u003e89\u003c/em\u003e, 828\u0026ndash;839. https://doi.org/10.1016/j.ecolind.2017.10.016\u003c/li\u003e\n \u003cli\u003eBreiman, L. (2001). Random forests. \u003cem\u003eMachine Learning\u003c/em\u003e, \u003cem\u003e45\u003c/em\u003e(1), 5\u0026ndash;32. https://doi.org/10.1023/A:1010933404324\u003c/li\u003e\n \u003cli\u003eCantonati, M., Stevens, L. E., Segadelli, S., Springer, A. E., Goldscheider, N., Celico, F., Filippini, M., Ogata, K., \u0026amp; Gargini, A. (2020). Ecohydrogeology: The interdisciplinary convergence needed to improve the study and stewardship of springs and other groundwater-dependent habitats, biota, and ecosystems. \u003cem\u003eEcological Indicators\u003c/em\u003e, \u003cem\u003e110\u003c/em\u003e. https://doi.org/10.1016/j.ecolind.2019.105803\u003c/li\u003e\n \u003cli\u003eChen, G., Singh, K. K., Lopez, J., \u0026amp; Zhou, Y. (2020). Tree canopy cover and carbon density are different proxy indicators for assessing the relationship between forest structure and urban socio-ecological conditions. \u003cem\u003eEcological Indicators\u003c/em\u003e, \u003cem\u003e113\u003c/em\u003e, 106279. https://doi.org/10.1016/J.ECOLIND.2020.106279\u003c/li\u003e\n \u003cli\u003eChignell, S. M., Luizza, M. W., Skach, S., Young, N. E., \u0026amp; Evangelista, P. H. (2018). An integrative modeling approach to mapping wetlands and riparian areas in a heterogeneous Rocky Mountain watershed. \u003cem\u003eRemote Sensing in Ecology and Conservation\u003c/em\u003e, \u003cem\u003e4\u003c/em\u003e(2), 150\u0026ndash;165. https://doi.org/https://doi.org/10.1002/rse2.63\u003c/li\u003e\n \u003cli\u003eChikofsky, E., \u0026amp; Cross II, J. H. (2002). Reverse Engineering. In \u003cem\u003eEncyclopedia of Software Engineering\u003c/em\u003e. https://doi.org/https://doi.org/10.1002/0471028959.sof285\u003c/li\u003e\n \u003cli\u003eCopernicus Programme. (2023). \u003cem\u003eCLC 2018 \u0026mdash; Copernicus Land Monitoring Service\u003c/em\u003e. https://land.copernicus.eu/pan-european/corine-land-cover/clc2018\u003c/li\u003e\n \u003cli\u003eCord, A. F., Brauman, K. A., Chaplin-Kramer, R., Huth, A., Ziv, G., \u0026amp; Seppelt, R. (2017). Priorities to Advance Monitoring of Ecosystem Services Using Earth Observation. \u003cem\u003eTrends in Ecology \u0026amp; Evolution\u003c/em\u003e, \u003cem\u003e32\u003c/em\u003e(6), 416\u0026ndash;428. https://doi.org/10.1016/j.tree.2017.03.003\u003c/li\u003e\n \u003cli\u003eCosta, H., Benevides, P., Moreira, F. D., Moraes, D., \u0026amp; Caetano, M. (2022). Spatially Stratified and Multi-Stage Approach for National Land Cover Mapping Based on Sentinel-2 Data and Expert Knowledge. \u003cem\u003eRemote Sensing 2022, Vol. 14, Page 1865\u003c/em\u003e, \u003cem\u003e14\u003c/em\u003e(8), 1865. https://doi.org/10.3390/RS14081865\u003c/li\u003e\n \u003cli\u003eCrowther, T. W., Glick, H. B., Covey, K. R., Bettigole, C., Maynard, D. S., Thomas, S. M., Smith, J. R., Hintler, G., Duguid, M. C., Amatulli, G., Tuanmu, M. N., Jetz, W., Salas, C., Stam, C., Piotto, D., Tavani, R., Green, S., Bruce, G., Williams, S. J., \u0026hellip; Bradford, M. A. (2015). Mapping tree density at a global scale. \u003cem\u003eNature\u003c/em\u003e, \u003cem\u003e525\u003c/em\u003e(7568), 201\u0026ndash;205. https://doi.org/10.1038/nature14967\u003c/li\u003e\n \u003cli\u003eCz\u0026uacute;cz, B., Keith, H., Driver, A., Jackson, B., Nicholson, E., \u0026amp; Maes, J. (2021). A common typology for ecosystem characteristics and ecosystem condition variables. \u003cem\u003eOne Ecosystem\u003c/em\u003e, \u003cem\u003e6\u003c/em\u003e, 1\u0026ndash;16. https://doi.org/10.3897/oneeco.6.e58218\u003c/li\u003e\n \u003cli\u003edel R\u0026iacute;o-Mena, T., Willemen, L., Vrieling, A., \u0026amp; Nelson, A. (2023). How remote sensing choices influence ecosystem services monitoring and evaluation results of ecological restoration interventions. \u003cem\u003eEcosystem Services\u003c/em\u003e, \u003cem\u003e64\u003c/em\u003e. https://doi.org/10.1016/J.ECOSER.2023.101565\u003c/li\u003e\n \u003cli\u003eDeLancey, E. R., Kariyeva, J., Bried, J. T., \u0026amp; Hird, J. N. (2019). Large-scale probabilistic identification of boreal peatlands using Google Earth Engine, open-access satellite data, and machine learning. \u003cem\u003ePLoS ONE\u003c/em\u003e, \u003cem\u003e14\u003c/em\u003e(6), 1\u0026ndash;23. https://doi.org/10.1371/journal.pone.0218165\u003c/li\u003e\n \u003cli\u003eDGT. (2018). \u003cem\u003eTechnical specifications of the land use and land cover map for Portugal Mainland\u0026nbsp;\u003c/em\u003e.\u003c/li\u003e\n \u003cli\u003eDwire, K. A., Mellmann-Brown, S., \u0026amp; Gurrieri, J. T. (2018). Potential effects of climate change on riparian areas, wetlands, and groundwater-dependent ecosystems in the Blue Mountains, Oregon, USA. In \u003cem\u003eClimate Services\u003c/em\u003e (Vol. 10, pp. 44\u0026ndash;52). Elsevier B.V. https://doi.org/10.1016/j.cliser.2017.10.002\u003c/li\u003e\n \u003cli\u003eEdens, B., Maes, J., Hein, L., Obst, C., Siikamaki, J., Schenau, S., Javorsek, M., Chow, J., Chan, J. Y., Steurer, A., \u0026amp; Alfieri, A. (2022). Establishing the SEEA Ecosystem Accounting as a global standard. \u003cem\u003eEcosystem Services\u003c/em\u003e, \u003cem\u003e54\u003c/em\u003e, 101413. https://doi.org/10.1016/J.ECOSER.2022.101413\u003c/li\u003e\n \u003cli\u003eEllison, D., Morris, C. E., Locatelli, B., Sheil, D., Cohen, J., Murdiyarso, D., Gutierrez, V., Noordwijk, M. van, Creed, I. F., Pokorny, J., Gaveau, D., Spracklen, D. V., Tobella, A. B., Ilstedt, U., Teuling, A. J., Gebrehiwot, S. G., Sands, D. C., Muys, B., Verbist, B., \u0026hellip; Sullivan, C. A. (2017). Trees, forests and water: Cool insights for a hot world. \u003cem\u003eGlobal Environmental Change\u003c/em\u003e, \u003cem\u003e43\u003c/em\u003e, 51\u0026ndash;61. https://doi.org/10.1016/J.GLOENVCHA.2017.01.002\u003c/li\u003e\n \u003cli\u003eESRI. (2023). \u003cem\u003eArcGIS Pro - ESRI (Environmental Systems Research Institute)\u003c/em\u003e. https://www.arcgis.com/index.html\u003c/li\u003e\n \u003cli\u003eEuropean Environment Agency. (2023). \u003cem\u003eCopernicus Land Monitoring Service\u003c/em\u003e.\u003c/li\u003e\n \u003cli\u003eEuropean Spatial Agency. (2023). \u003cem\u003eESA - Space Science\u003c/em\u003e. https://www.esa.int/Science_Exploration/Space_Science\u003c/li\u003e\n \u003cli\u003eFarrell, C., Coleman, L., Kelly-Quinn, M., Obst, C., Eigenraam, M., Norton, D., ODonoghue, C., Kinsella, S., Delargy, O., \u0026amp; Stout, J. (2021). Applying the System of Environmental Economic Accounting-Ecosystem Accounting (SEEA-EA) framework at catchment scale to develop ecosystem extent and condition accounts. \u003cem\u003eOne Ecosystem\u003c/em\u003e, \u003cem\u003e6\u003c/em\u003e. https://doi.org/10.3897/oneeco.6.e65582\u003c/li\u003e\n \u003cli\u003eFick, S. E., \u0026amp; Hijmans, R. J. (2017). WorldClim 2: new 1-km spatial resolution climate surfaces for global land areas. \u003cem\u003eInternational Journal of Climatology\u003c/em\u003e, \u003cem\u003e37\u003c/em\u003e(12), 4302\u0026ndash;4315. https://doi.org/https://doi.org/10.1002/joc.5086\u003c/li\u003e\n \u003cli\u003eFleming, A., O\u0026rsquo;Grady, A. P., Stitzlein, C., Ogilvy, S., Mendham, D., \u0026amp; Harrison, M. T. (2022). Improving acceptance of natural capital accounting in land use decision making: Barriers and opportunities. \u003cem\u003eEcological Economics\u003c/em\u003e, \u003cem\u003e200\u003c/em\u003e, 107510. https://doi.org/10.1016/J.ECOLECON.2022.107510\u003c/li\u003e\n \u003cli\u003eFonseca, F., de Figueiredo, T., Vilela, \u0026Acirc;., Santos, R., de Carvalho, A. L., Almeida, E., \u0026amp; Nunes, L. (2019). Impact of tree species replacement on carbon stocks in a Mediterranean mountain area, NE Portugal. \u003cem\u003eForest Ecology and Management\u003c/em\u003e, \u003cem\u003e439\u003c/em\u003e, 181\u0026ndash;188. https://doi.org/https://doi.org/10.1016/j.foreco.2019.03.002\u003c/li\u003e\n \u003cli\u003eGitelson, A. A., Gritz, Y., \u0026amp; Merzlyak, M. N. (2003). Relationships between leaf chlorophyll content and spectral reflectance and algorithms for non-destructive chlorophyll assessment in higher plant leaves. \u003cem\u003eJournal of Plant Physiology\u003c/em\u003e, \u003cem\u003e160\u003c/em\u003e(3), 271\u0026ndash;282. https://doi.org/10.1078/0176-1617-00887\u003c/li\u003e\n \u003cli\u003eGoogle. (2023). \u003cem\u003eGoogle Earth Engine\u003c/em\u003e. https://earthengine.google.com/\u003c/li\u003e\n \u003cli\u003eGorelick, N., Hancher, M., Dixon, M., Ilyushchenko, S., Thau, D., \u0026amp; Moore, R. (2017). Google Earth Engine: Planetary-scale geospatial analysis for everyone. \u003cem\u003eRemote Sensing of Environment\u003c/em\u003e, \u003cem\u003e202\u003c/em\u003e, 18\u0026ndash;27. https://doi.org/https://doi.org/10.1016/j.rse.2017.06.031\u003c/li\u003e\n \u003cli\u003eGosal, A. S., Evans, P. M., Bullock, J. M., Redhead, J., Charlton, M. B., Cord, A. F., Johnson, A., \u0026amp; Ziv, G. (2022). Understanding the accuracy of modelled changes in freshwater provision over time. \u003cem\u003eScience of The Total Environment\u003c/em\u003e, \u003cem\u003e833\u003c/em\u003e, 155042. https://doi.org/10.1016/j.scitotenv.2022.155042\u003c/li\u003e\n \u003cli\u003eHarris, L. D. (1984). \u003cem\u003eIsland Biogeography Theory and the Preservation of Biotic Diversity\u003c/em\u003e. University of Chicago Press. https://doi.org/doi:10.7208/9780226219950\u003c/li\u003e\n \u003cli\u003eHein, L., Bagstad, K. J., Obst, C., Edens, B., Schenau, S., Castillo, G., Soulard, F., Brown, C., Driver, A., Bordt, M., Steurer, A., Harris, R., \u0026amp; Caparr\u0026oacute;s, A. (2020). Progress in natural capital accounting for ecosystems. \u003cem\u003eScience\u003c/em\u003e, \u003cem\u003e367\u003c/em\u003e(6477), 514\u0026ndash;515. https://doi.org/10.1126/science.aaz8901\u003c/li\u003e\n \u003cli\u003eHose, G. C., Dabovic, J., Nelson, T., Pollitt, L., \u0026amp; Korbel, K. L. (2021). Groundwater Dependent Aquatic and Terrestrial Ecosystems. \u003cem\u003eReference Module in Earth Systems and Environmental Sciences\u003c/em\u003e. https://doi.org/10.1016/B978-0-12-819166-8.00118-3\u003c/li\u003e\n \u003cli\u003eHosmer, D. W., \u0026amp; Lemeshow, S. (2000). Interpretation of the Fitted Logistic Regression Model. In \u003cem\u003eApplied Logistic Regression\u003c/em\u003e (pp. 47\u0026ndash;90). Wiley. https://doi.org/10.1002/0471722146.ch3\u003c/li\u003e\n \u003cli\u003eHossain, M. S., \u0026amp; Hashim, M. (2019). Potential of Earth Observation (EO) technologies for seagrass ecosystem service assessments. \u003cem\u003eInternational Journal of Applied Earth Observation and Geoinformation\u003c/em\u003e, \u003cem\u003e77\u003c/em\u003e, 15\u0026ndash;29. https://doi.org/10.1016/j.jag.2018.12.009\u003c/li\u003e\n \u003cli\u003eHuang, F., Chunyu, X., Zhang, D., Chen, X., \u0026amp; Ochoa, C. G. (2020). A framework to assess the impact of ecological water conveyance on groundwater-dependent terrestrial ecosystems in arid inland river basins. \u003cem\u003eScience of the Total Environment\u003c/em\u003e, \u003cem\u003e709\u003c/em\u003e. https://doi.org/10.1016/j.scitotenv.2019.136155\u003c/li\u003e\n \u003cli\u003eHuete, A. R. (1988). A soil-adjusted vegetation index (SAVI). \u003cem\u003eRemote Sensing of Environment\u003c/em\u003e, \u003cem\u003e25\u003c/em\u003e(3), 295\u0026ndash;309. https://doi.org/10.1016/0034-4257(88)90106-X\u003c/li\u003e\n \u003cli\u003eHuete, A. R., Liu, H. Q., \u0026amp; van Leeuwen, W. J. D. (1997). Use of vegetation indices in forested regions: Issues of linearity and saturation. \u003cem\u003eInternational Geoscience and Remote Sensing Symposium (IGARSS)\u003c/em\u003e, \u003cem\u003e4\u003c/em\u003e, 1966\u0026ndash;1968. https://doi.org/10.1109/IGARSS.1997.609169\u003c/li\u003e\n \u003cli\u003eJagannathan, J., \u0026amp; Divya, C. (2021). Deep learning for the prediction and classification of land use and land cover changes using deep convolutional neural network. \u003cem\u003eEcological Informatics\u003c/em\u003e, \u003cem\u003e65\u003c/em\u003e, 101412. https://doi.org/10.1016/J.ECOINF.2021.101412\u003c/li\u003e\n \u003cli\u003eJAXA. (2023). \u003cem\u003eALOS World 3D - 30m\u003c/em\u003e. https://www.eorc.jaxa.jp/ALOS/en/dataset/aw3d30/aw3d30_e.htm\u003c/li\u003e\n \u003cli\u003eJi, L., Geng, X., Sun, K., Zhao, Y., \u0026amp; Gong, P. (2015). Target Detection Method for Water Mapping Using Landsat 8 OLI/TIRS Imagery. \u003cem\u003eWater\u003c/em\u003e, \u003cem\u003e7\u003c/em\u003e(2), 794\u0026ndash;817. https://doi.org/10.3390/w7020794\u003c/li\u003e\n \u003cli\u003eKeith, H., Cz\u0026uacute;cz, B., Jackson, B., Driver, A., Nicholson, E., \u0026amp; Maes, J. (2020). A conceptual framework and practical structure for implementing ecosystem condition accounts. \u003cem\u003eOne Ecosystem\u003c/em\u003e, \u003cem\u003e5\u003c/em\u003e, 1\u0026ndash;54. https://doi.org/10.3897/oneeco.5.e58216\u003c/li\u003e\n \u003cli\u003eKeith, H., Vardon, M., Obst, C., Young, V., Houghton, R., \u0026amp; Mackey, B. (2021). Evaluating nature-based solutions for climate mitigation and conservation requires comprehensive carbon accounting. \u003cem\u003eScience of The Total Environment\u003c/em\u003e, \u003cem\u003e769\u003c/em\u003e, 144341. https://doi.org/10.1016/j.scitotenv.2020.144341\u003c/li\u003e\n \u003cli\u003eKhalid, H. W., Khalil, R. M. Z., \u0026amp; Qureshi, M. A. (2021). Evaluating spectral indices for water bodies extraction in western Tibetan Plateau. \u003cem\u003eThe Egyptian Journal of Remote Sensing and Space Science\u003c/em\u003e, \u003cem\u003e24\u003c/em\u003e(3), 619\u0026ndash;634. https://doi.org/10.1016/J.EJRS.2021.09.003\u003c/li\u003e\n \u003cli\u003eKienast, F., Bolliger, J., Potschin, M., de Groot, R. S., Verburg, P. H., Heller, I., Wascher, D., \u0026amp; Haines-Young, R. (2009). Assessing Landscape Functions with Broad-Scale Environmental Data: Insights Gained from a Prototype Development for Europe. \u003cem\u003eEnvironmental Management\u003c/em\u003e, \u003cem\u003e44\u003c/em\u003e(6), 1099\u0026ndash;1120. https://doi.org/10.1007/s00267-009-9384-7\u003c/li\u003e\n \u003cli\u003eKl\u0026oslash;ve, B., Ala-Aho, P., Bertrand, G., Gurdak, J. J., Kupfersberger, H., Kv\u0026aelig;rner, J., Muotka, T., Mykr\u0026auml;, H., Preda, E., Rossi, P., Uvo, C. B., Velasco, E., \u0026amp; Pulido-Velazquez, M. (2014). Climate change impacts on groundwater and dependent ecosystems. \u003cem\u003eJournal of Hydrology\u003c/em\u003e, \u003cem\u003e518\u003c/em\u003e(PB), 250\u0026ndash;266. https://doi.org/10.1016/J.JHYDROL.2013.06.037\u003c/li\u003e\n \u003cli\u003eKoit, O., Tarros, S., P\u0026auml;rn, J., K\u0026uuml;ttim, M., Abreldaal, P., Sisask, K., Vainu, M., Terasmaa, J., Retike, I., \u0026amp; Polikarpus, M. (2021). Contribution of local factors to the status of a groundwater dependent terrestrial ecosystem in the transboundary Gauja-Koiva River basin, North-Eastern Europe. \u003cem\u003eJournal of Hydrology\u003c/em\u003e, \u003cem\u003e600\u003c/em\u003e. https://doi.org/10.1016/j.jhydrol.2021.126656\u003c/li\u003e\n \u003cli\u003eKoma, Z., Seijmonsbergen, A. C., \u0026amp; Kissling, W. D. (2021). Classifying wetland-related land cover types and habitats using fine-scale lidar metrics derived from country-wide Airborne Laser Scanning. \u003cem\u003eRemote Sensing in Ecology and Conservation\u003c/em\u003e, \u003cem\u003e7\u003c/em\u003e(1), 80\u0026ndash;96. https://doi.org/https://doi.org/10.1002/rse2.170\u003c/li\u003e\n \u003cli\u003eKulawardhana, R. W. (2011). Remote sensing of vegetation: principles, techniques and applications. By Hamlyn G. Jones and Robin A Vaughan. \u003cem\u003eJournal of Vegetation Science\u003c/em\u003e, \u003cem\u003e22\u003c/em\u003e(6), 1151\u0026ndash;1153. https://doi.org/https://doi.org/10.1111/j.1654-1103.2011.01319.x\u003c/li\u003e\n \u003cli\u003eLange, S., Campagne, C. S., Comte, A., Bank, E., Santos-Mart\u0026iacute;n, F., Maes, J., \u0026amp; Burkhard, B. (2022). Progress on ecosystem accounting in Europe. \u003cem\u003eEcosystem Services\u003c/em\u003e, \u003cem\u003e57\u003c/em\u003e, 101473. https://doi.org/10.1016/J.ECOSER.2022.101473\u003c/li\u003e\n \u003cli\u003eLiquete, C., Maes, J., Notte, A., \u0026amp; Bidoglio, G. (2011). Securing water as a resource for society: an ecosystem services perspective. \u003cem\u003eEcohydrology \u0026amp; Hydrobiology\u003c/em\u003e, \u003cem\u003e11\u003c/em\u003e(3\u0026ndash;4), 247\u0026ndash;259. https://doi.org/10.2478/V10104-011-0044-1\u003c/li\u003e\n \u003cli\u003eLudwig, C., Walli, A., Schleicher, C., Weichselbaum, J., \u0026amp; Riffler, M. (2019). A highly automated algorithm for wetland detection using multi-temporal optical satellite data. \u003cem\u003eRemote Sensing of Environment\u003c/em\u003e, \u003cem\u003e224\u003c/em\u003e, 333\u0026ndash;351. https://doi.org/10.1016/J.RSE.2019.01.017\u003c/li\u003e\n \u003cli\u003eLundqvist, J., \u0026amp; Unver, O. (2018). Alternative pathways to food security and nutrition \u0026ndash; water predicaments and human behavior. \u003cem\u003eWater Policy\u003c/em\u003e, \u003cem\u003e20\u003c/em\u003e(5), 871\u0026ndash;884. https://doi.org/10.2166/WP.2018.171\u003c/li\u003e\n \u003cli\u003eMadrigal-Gonz\u0026aacute;lez, J., Calatayud, J., Ballesteros-C\u0026aacute;novas, J. A., Escudero, A., Cayuela, L., Marqu\u0026eacute;s, L., Rueda, M., Ruiz-Benito, P., Herrero, A., Aponte, C., Sagardia, R., Plumptre, A. J., Dupire, S., Espinosa, C. I., Tutubalina, O. V., Myint, M., Pataro, L., L\u0026oacute;pez-S\u0026aacute;ez, J., Mac\u0026iacute;a, M. J., \u0026hellip; Stoffel, M. (2023). Global patterns of tree density are contingent upon local determinants in the world\u0026rsquo;s natural forests. \u003cem\u003eCommunications Biology 2023 6:1\u003c/em\u003e, \u003cem\u003e6\u003c/em\u003e(1), 1\u0026ndash;6. https://doi.org/10.1038/s42003-023-04419-8\u003c/li\u003e\n \u003cli\u003eMairota, P., Cafarelli, B., Didham, R. K., Lovergine, F. P., Lucas, R. M., Nagendra, H., Rocchini, D., \u0026amp; Tarantino, C. (2015). Challenges and opportunities in harnessing satellite remote-sensing for biodiversity monitoring. \u003cem\u003eEcological Informatics\u003c/em\u003e, \u003cem\u003e30\u003c/em\u003e, 207\u0026ndash;214. https://doi.org/10.1016/j.ecoinf.2015.08.006\u003c/li\u003e\n \u003cli\u003eMengist, W., Soromessa, T., \u0026amp; Feyisa, G. L. (2021). Monitoring Afromontane forest cover loss and the associated socio-ecological drivers in Kaffa biosphere reserve, Ethiopia. \u003cem\u003eTrees, Forests and People\u003c/em\u003e, \u003cem\u003e6\u003c/em\u003e. https://doi.org/10.1016/j.tfp.2021.100161\u003c/li\u003e\n \u003cli\u003eMpakairi, K. S., Dube, T., Dondofema, F., \u0026amp; Dalu, T. (2022). Spatio\u0026ndash;temporal variation of vegetation heterogeneity in groundwater dependent ecosystems within arid environments. \u003cem\u003eEcological Informatics\u003c/em\u003e, \u003cem\u003e69\u003c/em\u003e, 101667. https://doi.org/10.1016/j.ecoinf.2022.101667\u003c/li\u003e\n \u003cli\u003eMulatu, K. A., Mora, B., Kooistra, L., \u0026amp; Herold, M. (2017). Biodiversity monitoring in changing tropical forests: A review of approaches and new opportunities. In \u003cem\u003eRemote Sensing\u003c/em\u003e (Vol. 9, Issue 10). MDPI AG. https://doi.org/10.3390/rs9101059\u003c/li\u003e\n \u003cli\u003eNewbold, T., Hudson, L. N., Hill, S. L. L., Contu, S., Lysenko, I., Senior, R. A., B\u0026ouml;rger, L., Bennett, D. J., Choimes, A., Collen, B., Day, J., De Palma, A., D\u0026iacute;az, S., Echeverria-Londo\u0026ntilde;o, S., Edgar, M. J., Feldman, A., Garon, M., Harrison, M. L. K., Alhusseini, T., \u0026hellip; Purvis, A. (2015). Global effects of land use on local terrestrial biodiversity. \u003cem\u003eNature\u003c/em\u003e, \u003cem\u003e520\u003c/em\u003e(7545), 45\u0026ndash;50. https://doi.org/10.1038/nature14324\u003c/li\u003e\n \u003cli\u003eNieves-Chinchilla, J., Mart\u0026iacute;nez, R., Farjas, M., Tubio-Pardavila, R., Cruz, D., \u0026amp; Gallego, M. (2018). Reverse engineering techniques to optimize facility location of satellite ground stations on building roofs. \u003cem\u003eAutomation in Construction\u003c/em\u003e, \u003cem\u003e90\u003c/em\u003e, 156\u0026ndash;165. https://doi.org/10.1016/J.AUTCON.2018.02.019\u003c/li\u003e\n \u003cli\u003eOliveira, M., Novo, M. emilia, \u0026amp; Lobo-Ferreira, J.-P. (2007). Models to predict the impact of the climate changes on aquifer recharge. \u003cem\u003eWater in Celtic Countries: Quantity, Quality and Climate Variability (Proceedings of the Fourth InterCeltic Colloquium on Hydrology and Management of Water Resources\u003c/em\u003e, 103\u0026ndash;110. https://iahs.info/uploads/dms/13794.15-103-110-05-IC310-23-Oliveira.pdf\u003c/li\u003e\n \u003cli\u003eRamirez-Reyes, C., Brauman, K. A., Chaplin-Kramer, R., Galford, G. L., Adamo, S. B., Anderson, C. B., Anderson, C., Allington, G. R. H., Bagstad, K. J., Coe, M. T., Cord, A. F., Dee, L. E., Gould, R. K., Jain, M., Kowal, V. A., Muller-Karger, F. E., Norriss, J., Potapov, P., Qiu, J., \u0026hellip; Wright, T. M. (2019). Reimagining the potential of Earth observations for ecosystem service assessments. \u003cem\u003eScience of The Total Environment\u003c/em\u003e, \u003cem\u003e665\u003c/em\u003e, 1053\u0026ndash;1063. https://doi.org/10.1016/j.scitotenv.2019.02.150\u003c/li\u003e\n \u003cli\u003eRijal, S., Rimal, B., Acharya, R. P., \u0026amp; Stork, N. E. (2021). Land use/land cover change and ecosystem services in the Bagmati River Basin, Nepal. \u003cem\u003eEnvironmental Monitoring and Assessment\u003c/em\u003e, \u003cem\u003e193\u003c/em\u003e(10), 651. https://doi.org/10.1007/s10661-021-09441-z\u003c/li\u003e\n \u003cli\u003eRouse, J., Haas, R., Schell J., \u0026amp; Deering D. (1973). \u003cem\u003eMonitoring Vegetation Systems in the Great Plains with ERTS. Third ERTS Symposium, NASA\u003c/em\u003e.\u003c/li\u003e\n \u003cli\u003eSantos, F. L. M., Couto, F. T., Dias, S. S., Ribeiro, N. A., \u0026amp; Salgado, R. (2023). Vegetation fuel characterization using machine learning approach over southern Portugal. \u003cem\u003eRemote Sensing Applications: Society and Environment\u003c/em\u003e, 101017. https://doi.org/10.1016/J.RSASE.2023.101017\u003c/li\u003e\n \u003cli\u003eSeddon, N., Smith, A., Smith, P., Key, I., Chausson, A., Girardin, C., House, J., Srivastava, S., \u0026amp; Turner, B. (2021). Getting the message right on nature-based solutions to climate change. \u003cem\u003eGlobal Change Biology\u003c/em\u003e, \u003cem\u003e27\u003c/em\u003e(8), 1518\u0026ndash;1546. https://doi.org/10.1111/GCB.15513\u003c/li\u003e\n \u003cli\u003eTabacchi, G., Di Cosmo, L., \u0026amp; Gasparini, P. (2011). Aboveground tree volume and phytomass prediction equations for forest species in Italy. \u003cem\u003eEuropean Journal of Forest Research\u003c/em\u003e, \u003cem\u003e130\u003c/em\u003e(6), 911\u0026ndash;934. https://doi.org/10.1007/S10342-011-0481-9/TABLES/8\u003c/li\u003e\n \u003cli\u003eUnited Nations et al. (2021). \u003cem\u003eSystem of Environmental-Economic Accounting-Ecosystem Accounting (SEEA-EA)\u003c/em\u003e. https://seea.un.org/ecosystem-accounting.\u003c/li\u003e\n \u003cli\u003eVallero, D. (2014). Grand Challenges. \u003cem\u003eFundamentals of Air Pollution\u003c/em\u003e, 953\u0026ndash;961. https://doi.org/10.1016/B978-0-12-401733-7.00033-5\u003c/li\u003e\n \u003cli\u003eVargas, L., Willemen, L., \u0026amp; Hein, L. (2019). Assessing the Capacity of Ecosystems to Supply Ecosystem Services Using Remote Sensing and An Ecosystem Accounting Approach. \u003cem\u003eEnvironmental Management\u003c/em\u003e, \u003cem\u003e63\u003c/em\u003e(1), 1\u0026ndash;15. https://doi.org/10.1007/s00267-018-1110-x\u003c/li\u003e\n \u003cli\u003eWang, J.-F., Zhang, T.-L., \u0026amp; Fu, B.-J. (2016). A measure of spatial stratified heterogeneity. \u003cem\u003eEcological Indicators\u003c/em\u003e, \u003cem\u003e67\u003c/em\u003e, 250\u0026ndash;256. https://doi.org/10.1016/j.ecolind.2016.02.052\u003c/li\u003e\n \u003cli\u003eWillcock, S., Mart\u0026iacute;nez-L\u0026oacute;pez, J., Hooftman, D. A. P., Bagstad, K. J., Balbi, S., Marzo, A., Prato, C., Sciandrello, S., Signorello, G., Voigt, B., Villa, F., Bullock, J. M., \u0026amp; Athanasiadis, I. N. (2018). Machine learning for ecosystem services. \u003cem\u003eEcosystem Services\u003c/em\u003e, \u003cem\u003e33\u003c/em\u003e, 165\u0026ndash;174. https://doi.org/10.1016/J.ECOSER.2018.04.004\u003c/li\u003e\n \u003cli\u003eWood, W. H. (2009). Computational Representations of Function in Engineering Design. \u003cem\u003ePhilosophy of Technology and Engineering Sciences\u003c/em\u003e, 543\u0026ndash;564. https://doi.org/10.1016/B978-0-444-51667-1.50024-0\u003c/li\u003e\n \u003cli\u003eYang, J., Wang, J., Liao, X., Tao, H., \u0026amp; Li, Y. (2022). Chain modeling for the biogeochemical nexus of cadmium in soil-rice-human health system. \u003cem\u003eEnvironment International\u003c/em\u003e, \u003cem\u003e167\u003c/em\u003e, 107424. https://doi.org/10.1016/j.envint.2022.107424\u003c/li\u003e\n \u003cli\u003eYang, Z., Bai, J., \u0026amp; Zhang, W. (2021). Mapping and assessment of wetland conditions by using remote sensing images and POI data. \u003cem\u003eEcological Indicators\u003c/em\u003e, \u003cem\u003e127\u003c/em\u003e, 107485. https://doi.org/10.1016/J.ECOLIND.2021.107485\u003c/li\u003e\n \u003cli\u003eZomer, R. J., Xu, J., \u0026amp; Trabucco, A. (2022). Version 3 of the Global Aridity Index and Potential Evapotranspiration Database. \u003cem\u003eScientific Data\u003c/em\u003e, \u003cem\u003e9\u003c/em\u003e(1), 409. https://doi.org/10.1038/s41597-022-01493-1\u003c/li\u003e\n\u003c/ol\u003e"}],"fulltextSource":"","fullText":"","funders":[],"hasAdminPriorityOnWorkflow":false,"hasManuscriptDocX":true,"hasOptedInToPreprint":true,"hasPassedJournalQc":"","hasAnyPriority":false,"hideJournal":true,"highlight":"","institution":"","isAcceptedByJournal":false,"isAuthorSuppliedPdf":false,"isDeskRejected":"","isHiddenFromSearch":false,"isInQc":false,"isInWorkflow":false,"isPdf":false,"isPdfUpToDate":true,"isWithdrawnOrRetracted":false,"journal":{"display":true,"email":"[email protected]","identity":"researchsquare","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":true,"externalIdentity":"","sideBox":"","snPcode":"","submissionUrl":"/submission","title":"Research Square","twitterHandle":"researchsquare","acdcEnabled":true,"dfaEnabled":false,"editorialSystem":"","reportingPortfolio":"","inReviewEnabled":false,"inReviewRevisionsEnabled":true},"keywords":"Environmental Modelling, Machine Learning, Google Earth Engine, Data fusion, Copernicus Land Monitoring Services","lastPublishedDoi":"10.21203/rs.3.rs-4681296/v1","lastPublishedDoiUrl":"https://doi.org/10.21203/rs.3.rs-4681296/v1","license":{"name":"CC BY 4.0","url":"https://creativecommons.org/licenses/by/4.0/"},"manuscriptAbstract":"\u003cp\u003eThe critical role of inland Green-Blue Ecosystems (GBE) in delivering goods and services underscores the need to understand their relationships. This research investigates the impact of climate variables on GBE stock accounts by estimating Tree Cover Density (TCD) and Water \u0026amp; Wetness Probability Index (WWPI). Using supervised machine learning and factor analysis, we measured GBE extent and condition. Key predictors for Blue Ecosystems were topographic variables, while vegetation indices were crucial for Green Ecosystems. In 2018, 33% of the inland area was covered by forests, freshwater, and wetlands. Key climate-driven factors for forests included precipitation (0.65), aridity index (0.54), and evapotranspiration (0.44). For freshwater and wetlands, precipitation (0.69), aridity index (0.55), and elevation (0.42) were significant. This research enhances our understanding of how climate impacts GBE, influencing biomass density and water availability. It bridges socio-environmental science with engineering by integrating advanced modelling techniques, promoting ecosystem resilience and sustainability.\u003c/p\u003e","manuscriptTitle":"Climate-driven stock accounts of inland Green-Blue Ecosystems: Coupling the carbon and water cycles via Reverse Engineering and Geodetectors","msid":"","msnumber":"","nonDraftVersions":[{"code":1,"date":"2024-09-02 06:50:13","doi":"10.21203/rs.3.rs-4681296/v1","editorialEvents":[{"type":"communityComments","content":0}],"status":"published","journal":{"display":true,"email":"[email protected]","identity":"researchsquare","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":true,"externalIdentity":"","sideBox":"","snPcode":"","submissionUrl":"/submission","title":"Research Square","twitterHandle":"researchsquare","acdcEnabled":true,"dfaEnabled":false,"editorialSystem":"","reportingPortfolio":"","inReviewEnabled":false,"inReviewRevisionsEnabled":true}}],"origin":"","ownerIdentity":"0e79206b-dd80-4ee8-a4ce-77ea22255b30","owner":[],"postedDate":"September 2nd, 2024","published":true,"recentEditorialEvents":[],"rejectedJournal":[],"revision":"","amendment":"","status":"posted","subjectAreas":[],"tags":[],"updatedAt":"2024-10-16T13:09:01+00:00","versionOfRecord":[],"versionCreatedAt":"2024-09-02 06:50:13","video":"","vorDoi":"","vorDoiUrl":"","workflowStages":[]},"version":"v1","identity":"rs-4681296","journalConfig":"researchsquare"},"__N_SSP":true},"page":"/article/[identity]/[[...version]]","query":{"redirect":"/article/rs-4681296","identity":"rs-4681296","version":["v1"]},"buildId":"qtupq5eGEP_6zYnWcrvyt","isFallback":false,"isExperimentalCompile":false,"dynamicIds":[84888],"gssp":true,"scriptLoader":[]}