1. Introduction
Climate change presents one of the most significant challenges for the global environment and society at large today [
1]. It impacts ecosystems and their biodiversity by stressing and sometimes even threatening current habitats due to heat and water stress. It can also cause health-related effects and socio-economic impacts induced by increasing heat waves, droughts, and flooding events [
2]. Although climate change effects have a global impact, these challenges must be addressed locally to build and foster resilience, be better prepared, and manage accompanying risks [
3]. Hence, urban and municipal administrations play a crucial role in implementing effective measures to protect their cities and their citizens, especially the more vulnerable population groups, against these threats. For instance, those who are vulnerable include children, the elderly, socioeconomically disadvantaged, disabled, or underinsured individuals, or those with certain medical conditions [
4]. Nevertheless, constrained resources at the municipal level, such as minimal or a lack of financial means, ability, and expertise, often impede progress. This is especially true for smaller, less well-off municipalities to advance and improve adaptative planning and climate change mitigation measures [
5].
The European Copernicus program, the earth observation branch of the European Union space program, celebrated its 25th anniversary in June 2023. Initially known as the Global Monitoring for Environment and Security Programme (GMES), Copernicus was introduced in 1998 with the goal of supplying environmental data to support a diverse range of fields, such as urban planning, agriculture, disaster relief, and climate change. Hence, the program aims to not only integrate and provide satellite data but also non-space data to provide insights from earth observation. There are several thematic platforms. One of them is the Copernicus Climate Change Service (C3S), which provides information and data specifically related to climate indicators. Others, for instance, focus on land, marine, air, or atmospheric monitoring [
6]. The platform offers a wide range of earth observation data, such as environmental and ecological projections and climate monitoring, stored in the Climate Data Store (CDS) with different spatial resolutions [
7]. Although some climate indicators and gridded products may have coarser resolution and may not capture all urban-scale details, they are still valuable. For example, these can be biophysical or climate indicators and parameters, which help to understand vegetation health responses to climate variations [
8]. Extensive data processing is needed and essential to gain new insights from the results and help decision-makers formulate new or better climate change adaptation and mitigation strategies or measures. Finer resolution data and a more detailed scale analysis shall also be conducted to better understand and learn from the data what happens at a local scale [
9].
Monitoring vegetation dynamics is useful for climate adaptation purposes, especially at the municipal level, considering vegetation provides a range of ecological benefits, such as reducing the urban heat island effect. Numerous studies have been conducted and emphasise the intrinsic relationship between climate and vegetation [
7,
8,
9], for example, by monitoring vegetation response to weather conditions [
10].
The climate envelope model is used by scientists to calculate scenarios and derive new insights and knowledge. The models describe the relationships between species occurrences and bioclimate variables. The derived results of the model may indicate where plant species can thrive under specific climate conditions. It may also help identify regions and plant species that are prone to being more vulnerable to climatic changes [
11,
12]. Despite their limitations, such as often being investigated under equilibrium conditions that do not account for competition, dispersal, or nutrient supply [
9], climate envelopes are quite useful in understanding vegetation responses to climate change overall.
This study aims to investigate the health status of the vegetation and its correlation with climate conditions in the respective study area, the City of Constance, at Lake Constance, in southern Germany (Chapter 2). We developed a systematic approach and a simple tool for monitoring vegetation changes and health status on a smaller urban scale with a finer resolution compared to the coarse climate data from the Climate Data Store. We downscaled leaf area index (LAI), a fraction of absorbed photosynthetically active radiation (FAPAR), and bioclimate indicators in coarse resolution data using finer resolution data obtained from satellite images and local authorities, including Sentinel-2, Landsat 8, and a digital elevation model (DEM), to generate vegetation health and climatic parameters [
10,
13,
14]. Using the climate envelope model, we used the result as the input for modelling vegetation–climate relationships. The processing steps were carried out through the model builder in ArcGIS Pro 3.2, which can be transformed into toolboxes and a series of Python codes [
15]. This approach can provide knowledge and tools to support municipal decision-makers in identifying vulnerable locations and vegetation types for effective climate adaptation strategies, thereby enhancing local climate resilience.
4. Results
4.1. Vegetation Health Proxies
Vegetation compositions in the City of Constance exhibit considerable diversity. In this study, we assessed vegetation health by examining satellite-derived biophysical parameters, specifically leaf area index (LAI) and fraction of absorbed photosynthetically active radiation (FAPAR).
4.1.1. Vegetation Indices, LAI, and FAPAR
The LAI and FAPAR products from CDS with 300 m resolution were insufficient for application in urban vegetation studies. Therefore, we used several downscaling methods to generate finer-resolution data, and an empirical approach was used to estimate LAI and FAPAR at a higher resolution. It adapted relationships from general linear regression equations derived from field data collection [
24].
Vegetation indices generated from Sentinel-2 at a resolution of 10 m were computed, representing various aspects of vegetation, including the MSI, NDVI, SAVI, PSRI, and RNDVI. These indices were derived from specific bands commonly associated with capturing vegetation characteristics based on their reflectance properties [
31]. NDVI and SAVI showed strong correlations with leaf area index (LAI) and fraction of absorbed photosynthetically active radiation (FAPAR), reaching 0.71 and 0.74, respectively. RNDVI exhibited slightly lower correlations with LAI and FAPAR, at 0.69 and 0.70, respectively. In contrast, MSI and PSRI showed negative correlations, with MSI having a correlation of 0.66 for LAI and 0.67 for FAPAR, while PSRI had the lowest correlation values, with 0.52 for LAI and 0.50 for FAPAR.
4.1.2. Linear Model
Linear regression equations were developed to predict LAI and FAPAR using combinations of these vegetation indices (
Table 5 and
Table 6). The most optimal regression equations have high R-squared values and low AIC values. In this study, explanatory regression was used to investigate the consistency of vegetation indices in predicting LAI and FAPAR. NDVI, SAVI, and PSRI showed significant and stable influences on LAI and FAPAR, with NDVI and SAVI demonstrating 100% positive relationships, while PSRI showed 91.67% positive and 8.33% negative relationships. The NDVI, PSRI, SAVI, and PSRI combinations produced the most optimal regression equations with high R-squared and low AIC values.
Considering the study area’s characteristics with varied vegetation cover and tree cover, the NDVI and PSRI equations were chosen. However, if the study area is primarily covered by low vegetation, SAVI is recommended for use. The vegetation index of SAVI attempts to minimise the influence of soil brightness in areas of low vegetation cover using a soil brightness correction factor.
The mathematical equations derived from 2018 data were used to estimate LAI and FAPAR in 2022, assuming no significant landscape changes in the City of Constance within the four years. The coefficient determination (R-squared) for the NDVI and PSRI equations was approximately 0.78, indicating the model’s accuracy based on CDS data as a reference.
The equations used for generating LAI and FAPAR were
The comparison of the originally downloaded LAI and FAPAR from CDS and the generated linear model LAI and FAPAR is shown in
Figure 3.
The performance of the linear model in estimating LAI and FAPAR was compared to the SNAP software results for the same year. The R-squared values for the linear model and SNAP-derived LAI were 0.65 and 0.92 for FAPAR. SNAP-derived LAI tended to produce higher values than the linear model, which is consistent with previous studies [
32,
33]. However, for the study area’s vegetation, which primarily consists of tree-covered forest and urban areas, both linear models and SNAP-derived LAI and FAPAR are feasible.
4.1.3. Vegetation Health Classification
The classification value for vegetation health status is presented in
Table 7, and the map of vegetation health classification is shown in
Figure 4. Similar studies using different datasets to estimate LAI have demonstrated comparable values for corresponding types of vegetation cover. The classification outcomes indicate the presence of negative LAI and FAPAR values across all vegetation types. Therefore, to understand the causes, we overlayed the results with land cover data. The overlay analysis showed that these areas were frequently found at mixed and misclassified pixels.
Healthy trees were predominantly found in forested areas with a dense concentration of trees, especially in urban forests. Stressed trees were commonly found in street areas and close to built-up regions, with dead branches as indicators of stress symptoms. Warmer temperatures can enhance the carbon assimilation rate, leading to enlarged canopy cover in trees. However, water deficits can cause defoliation, especially in urban areas where rising temperatures may result in early summer defoliation.
While stressed trees in urban areas may face challenges, their presence is crucial in mitigating urban heat and creating a cooling effect. Understanding the correlation between vegetation and climate is essential, especially in urban areas, as urban vegetation can mitigate the negative impacts of climate change while being vulnerable to increasing temperatures and drought events.
In other vegetations, the difference between healthy and stressed vegetation is based on the amount of vegetation cover. Stressed vegetation was predominantly found in grasslands with high soil reflectance and unplanted cropland areas, even during peak summer when satellite images were acquired (July).
4.2. Vegetation and Climate Relations
Various bioclimate parameters were utilised across different vegetation types to identify optimal climate conditions conducive to the healthy growth of each vegetation type. This investigation aimed to establish the correlation between local-scale climate conditions and vegetation.
4.2.1. Vegetation Climate Model
The vegetation climate model identifies optimal climatic conditions for vegetation growth and adaptation to extreme climate changes, as shown in
Figure 5a. It defines climatic boundaries for trees and other vegetation in the City of Constance, assuming they will not grow if local climate variables exceed those defining its envelope. The range value of the bioclimate envelope is shown in
Table 8.
This model predicts the natural location range of vegetation to grow in healthy and stressed conditions. It assumes vegetation can grow well within its predicted natural conditions. However, specific adjustments may be needed outside of these locations. The model was combined with LAI- and FAPAR-derived vegetation health status to determine the probability of healthy and stressed vegetation, as shown in
Figure 5b. The model achieved high r-squared values for training (0.97) and validation (0.84) data, utilising six downscaled bioclimatic indicators as explanatory variables.
The bioclimate envelope model outlines the potential occurrence range for each vegetation type based on six variables’ maximum and minimum values. Minimum precipitation during the driest month (BIO14) represents the maximum drought vegetation can withstand, with trees and other vegetation requiring 0.02 mm. The upper limit of precipitation indicates the drought level inducing dormancy. Trees show better acclimatisation ability under changing climates (higher BIO14) than other vegetation types. The mean temperature during the warmest quarter (BIO10) is relatively similar among trees and other vegetation, with trees having a slightly wider range. Warmer temperatures within an optimal range stimulate photosynthetic activities. The lower limit represents the minimum temperature requirement for growth. The annual moisture index is useful for predicting vegetation types and forest areas. The minimum temperature of the coldest month (BIO06) can help identify anomalies in cold temperatures that may impact vegetation, while the annual temperature range (BIO07) indicates the potential effects of extreme temperatures on vegetation. BIO10 and BIO14 indicate the ability of that vegetation to withstand the average warmest period and the driest month.
4.2.2. Model Validation
The climate envelope model was validated by comparing it with existing vegetation from ESA WorldCover 2020, resulting in an overall accuracy of 86.7%. The result was summarised in the confusion matrix, as shown in
Table 9.
The kappa value, representing the level of agreement beyond chance, was 0.74, which is slightly lower than the overall accuracy. The model utilised a bioclimatic envelope concept, associating various climate aspects with species occurrences to estimate suitable conditions for vegetation. Producer accuracy for vegetation types shown in the climate envelope and ESA WorldCover ranged from 70.1% to 75.7%, indicating the percentage of reference pixels classified correctly. User accuracy varied for trees and other vegetation: 83.1% for other vegetation, 93.4% for trees, and 98.6% for non-vegetation.
4.3. Integrate Workflows into a Toolbox
This study utilised a range of global climatic data from CDS and other resources for climate monitoring, integrating different data sources to improve the results at the municipal level. Multiple processes of data acquisition and processing were employed in this study. The main process involved sequences of geoprocessing performed using model builders from ArcGIS Pro combined with a Python notebook. The model builder is an automated tool that connects data and available tools in ArcGIS Pro to execute workflows efficiently.
We created a toolbox called vegetation health, containing four toolsets representing different geoprocessing steps: vegetation indices (
Figure 6), biophysical processors (
Figure 7), vegetation health (
Figure 8), downscale climate indicators (
Figure 9a), and climate envelopes (
Figure 9b).
The biophysical processor toolset generated LAI and FAPAR using a model builder with a raster calculator and equations derived from general linear regression. The chosen vegetation indices (NDVI and PSRI) better predicted LAI and FAPAR than others.
The vegetation health toolset consisted of four model builders to classify vegetation health status based on LAI and FAPAR values. The process involved separating each vegetation type, identifying their health status, combining all types into a single map, and classifying vegetation health using LAI and FAPAR.
The downscale climate indicators toolset downscaled bioclimate indicators using EBK regression prediction with coarse and fine resolution rasters as explanatory variables, producing bioclimate indicators in finer resolution. The climate envelope toolset executed forest-based classification and regression, building a vegetation climate model called the climate envelope. It used vegetation-type point samples and bioclimate indicators in finer resolution as explanatory variables, providing a vegetation map showing the probability of healthy and stressed vegetation growth.
5. Discussion
5.1. Consideration of Approach and Interpretation of Vegetation Climate Model
Climate factors play a crucial role in influencing vegetation greenness, which is one of the key indicators of vegetation health. We used greenness and senescence indexes for assessing vegetation health in this study. The vegetation response to climate variation and adaptability is complex yet challenging to accurately simulate. Most studies assumed that vegetation has a fixed response pattern to climate change. Despite its limitations, the vegetation–climate relations model is important to understand how it impacts the environment.
Understanding how vegetation responds to local weather changes is crucial in microclimate. Microclimate refers to localised variations in heat and water moisture levels near the earth’s surface, leading to temperature and humidity differences compared to the surrounding areas. This local atmospheric condition can be influenced by a range of factors, such as energy absorption, shading, and wind speeds, which either trap or remove heat and moisture. The variations surrounding vegetation could potentially influence vegetation health, with healthier vegetation located near denser vegetation and stressed vegetation in isolated areas. Vegetation growing close to dense vegetation benefits from shading and moisture. In contrast, isolated vegetation surrounded by non-vegetation areas tends to experience more stress. It was shown in the vegetation health classification from the linear models of LAI and FAPAR that stressed vegetation is often found farther away from other vegetated areas.
This study combined the big data platform and local data to provide an adequate municipal-level model on a finer resolution scale. We performed multiple processes to downscale existing coarse data using remote sensing and local data. Statistical downscaling methods were employed, representing a flexible and straightforward approach to enhance the data of coarse resolution. Notably, though our models can be used for monitoring vegetation health, they cannot thoroughly describe the relationship between the climate and vegetation, especially for non-trees, which consist of various types of vegetation, such as wetland-sensitive biomes to temperature and with the highest interannual variability.
The climate envelope model was valuable for assessing suitable tree and other vegetation locations based on current climate conditions. For example, areas close to parks can be cooler during daylight hours compared to rooftops due to the cooling effect of transpiration. Sparse foliage areas exhibit higher temperatures due to less evaporation than areas covered by dense foliage. Besides the climate factors, elevated atmospheric CO2 concentration, varying nitrogen deposition rates, land use, and other anthropic factors could also influence vegetation health, which may bring a greater potential for vegetation change due to the more complex factors. We have not considered these additional factors in the model in this study. Further, we can explore the more complex social-ecological systems by inputting more natural and anthropogenic variables and coupling them with specific physical processes based on advanced modelling to better understand the complex relationship between vegetation and environments.
5.2. Additional Management Considerations
We produced linear model equations and a vegetation–climate model that exhibits potential applicability to municipalities sharing similarities with the City of Constance. This contribution extends beyond the specific study area, offering a transferable framework for regions exhibiting comparable characteristics. Additionally, the incorporation of a toolbox into our methodology serves as an effective means to introduce these models to municipal stakeholders lacking expertise in spatial data analysis. This approach enhances accessibility, facilitating wider utilisation and implementation in diverse municipal contexts.
Vegetation climate models at the municipal level offer valuable contributions to the conservation and management of urban ecosystems. The bioclimate envelope models derived from this study have numerous sources of uncertainty, including method choices, data used, and climate dynamics. Despite its limitations, the analysis provides useful tools for assessing climate impacts on urban vegetation, especially the information on the potential location for vegetation growth. This information empowers decision-makers to explore existing ecosystems within these environments, understand their function, and gather essential information for adaptation measures and planning. Integrating non-climatic factors and adaptive capacity information enhances the potential for conducting comprehensive ecological climate change vulnerability assessment in urban environments, presenting a crucial step towards effective urban ecosystem conservation.