# Estimation of Urban Tree Chlorophyll Content and Leaf Area Index Using Sentinel-2 Images and 3D Radiative Transfer Model Inversion

Théo Le Saint, Jean Nabucet, Laurence Hubert-Moy, Karine Adeline

► To cite this version:

Théo Le Saint, Jean Nabucet, Laurence Hubert-Moy, Karine Adeline. Estimation of Urban Tree Chlorophyll Content and Leaf Area Index Using Sentinel-2 Images and 3D Radiative Transfer Model Inversion. Remote Sensing, 2024, 16 (20), pp.3867. 10.3390/rs16203867. hal-04879418

HAL Id: hal-04879418

https://hal.science/hal-04879418v1

Submitted on 10 Jan 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Article

# Estimation of Urban Tree Chlorophyll Content and Leaf Area Index Using Sentinel-2 Images and 3D Radiative Transfer Model Inversion

Théo Le Saint<sup>1,2,*</sup>, Jean Nabucet<sup>1</sup>, Laurence Hubert-Moy<sup>1</sup> and Karine Adeline<sup>2 ID</sup>

<figure><img src="image_3.png" /></figure>

check for
updates

**Citation:** Le Saint, T.; Nabucet, J.; Hubert-Moy, L.; Adeline, K. Estimation of Urban Tree Chlorophyll Content and Leaf Area Index Using Sentinel-2 Images and 3D Radiative Transfer Model Inversion. *Remote Sens.* 2024, *16*, 3867. https://doi.org/10.3390/rs16203867

Academic Editors: Paolo Santi, Amin Anjomshoaa and Priyanka Nadia DeSouza

Received: 6 September 2024

Revised: 12 October 2024

Accepted: 16 October 2024

Published: 18 October 2024

<figure><img src="image_4.png" /></figure>

**Copyright:** © 2024 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

**Abstract:** Urban trees play an important role in mitigating effects of climate change and provide essential ecosystem services. However, the urban environment can stress trees, requiring the use of effective monitoring methods to assess their health and functionality. The objective of this study, which focused on four deciduous tree species in Rennes, France, was to evaluate the ability of hybrid inversion models to estimate leaf chlorophyll content (LCC), leaf area index (LAI), and canopy chlorophyll content (CCC) of urban trees using eight Sentinel-2 (S2) images acquired in 2021. Simulations were performed using the 3D radiative transfer model DART, and the hybrid inversion models were developed using machine-learning regression algorithms (random forest (RF) and gaussian process regression). Model performance was assessed using in situ measurements, and relations between satellite data and in situ measurements were investigated using spatial allocation (SA) methods at the pixel and tree scales. The influence of including environment features (EFs) as model inputs was also assessed. The results indicated that random forest models that included EFs and used the pixel-scale SA method were the most accurate with $R^2$ values of 0.33, 0.29, and 0.46 for LCC, LAI, and CCC, respectively, with notable variability among species.

**Keywords:** DART; random forest regression; gaussian process regression; vegetation traits; satellite time series

## 1. Introduction

In the realm of climate-change mitigation, urban trees play a central role in the strategies of land-use planning because they provide a variety of essential ecosystem services [1,2], such as temperature regulation by providing shade and evapotranspiration [3], carbon storage [4], and biodiversity preservation [5]. As a result, the presence of vegetation in urban areas is crucial for human health and the overall quality of life [6]. However, the urban environment contains many stress factors that can slow tree development [7]. Urban trees can experience mechanical and chemical disturbances due to planting conditions (proximity to buildings, edaphic conditions), and atmospheric conditions such as temperature and air quality, and light pollution at night. While urban trees contribute significantly to urban well-being, they must contend with a limiting environment that can damage their health. Several studies demonstrate notable impacts of these stress factors: urban trees have a shorter lifespan than rural trees [8], isolated and street trees experience more stress than those in parks [9], and young trees have a higher mortality rate in urban areas [10].

In this context, effectively monitoring the functional and mechanical condition of urban trees at a large scale is challenging. Traditional field inventories of urban trees are usually limited to a few trees and are performed annually or even every 10 years because they require large amounts of human resources and time. Conversely, satellite

<footnote n="1">1 UMR 6554 CNRS, LETG, University of Rennes, Place du Recteur Henri Le Moal, 35000 Rennes, France; jean.nabucet@univ-rennes2.fr (J.N.); laurence.moy@univ-rennes2.fr (L.H.-M.)</footnote>

<footnote n="2">2 DOTA, ONERA, Université de Toulouse, 31055 Toulouse, France; karine.adeline@onera.fr</footnote>

<footnote>* Correspondence: theo.le-saint@univ-rennes2.fr</footnote>
remote sensing can collect information on urban vegetation across entire cities with both
high spatial and temporal resolutions. Over the past decade, the use of satellite remote
sensing in urban areas for vegetation-related purposes has grown steadily [11,12]. Most
studies have focused on mapping large urban green spaces and classifying the species
in them. Relatively fewer studies have concentrated on smaller vegetation patches such
as isolated or aligned street trees, despite a large amount of cumulative coverage of such
patches. Furthermore, few studies have investigated the functional status of urban trees to
assess health or estimate biomass, instead focussing mainly on large green infrastructures.
Vegetation status can be characterized using proxies, which can be estimated using models
based on satellite images. The most commonly used vegetation traits are (i) leaf chlorophyll
content (LCC), often expressed in µg of chlorophyll a and b per cm² of leaf; (ii) leaf area
index (LAI), which corresponds to the sum of leaf areas of the tree crown relative to the
tree crown projection at ground area expressed in m²/m²; and (iii) canopy chlorophyll
content (CCC), expressed in µg/cm², which is the product of LCC and LAI. LCC is closely
linked to photosynthetic activity [13,14], biomass production [15,16], and levels of air
and soil pollution [17]. Because chlorophyll pigments are sensitive to changes in external
conditions, LCC is a good indicator of environmental stress and changes in temperature
and humidity [18]. At the canopy scale, LAI can be used to identify tree phenological
stages [19] and their photosynthetic potential [20]. CCC provides a more holistic and
accurate representation of chlorophyll content at the canopy scale than leaf-scale estimates.
These three vegetation traits can be estimated using satellite hyperspectral data due to
the latter’s wide spectral range (0.4–2.5 µm) [21]. However, current satellite hyperspectral
images, such as EnMAP and PRISMA images, have spatial and temporal resolutions that
are too coarse (30 m and 27–29 days at nadir) for the study of urban trees. While high or
very high spatial resolution (0.5–5 m) multispectral satellite imagery such as WorldView,
RapidEye, Pleiades, QuickBird, and PlanetScope are increasingly used to study vegetation
in urban environments, these data are not freely accessible conversely to Sentinel-2 (S2)
imagery. S2 stands out for its high spatial resolution (10 m for visible and near infrared
and 20 m for red edge (RE) (0.67–0.76 µm)/short-wave infrared), relatively large number
of spectral bands (10 bands at 10–20 m resolution), and temporal resolution suited to
monitoring intra-annual vegetation dynamics [22].

The methods used most often to estimate vegetation traits using S2 data can be
classified into two main types [23]: empirical or physical. Empirical methods are based on
regression models that derive empirical relationships between spectral bands or vegetation
indices and the variable of interest (e.g., LCC, LAI). Statistical models are widely used
due to their accessibility and ease of application. However, models developed for specific
locations or sensor configurations face limitations when applied in other contexts, such
as different vegetation types, sensor configurations, or data acquisition times. These
limitations arise from the need for large and heterogeneous in situ datasets for model
training and validation.

Physical methods are based on the use of radiative transfer models (RTM) to generate
synthetic spectral datasets or Look-Up Tables (LUT) mimicking the simulation of remote
sensing images. RTM models simulate the physical processes of electromagnetic radiation
in a medium. In remote sensing, these models are used to generate synthetic images based
on the optical and geometric properties of the scene and according to the characteristics
of a given sensor. They can be divided into two subcategories: LUT-based models and
hybrid models. LUT-based methods rely on the minimization of a cost function computed
between the simulated (from LUT) and measured spectra for each pixel of the remote
sensing image. Post-processing can include the use of adapted spectral intervals or the
use of vegetation indices. The best fit should provide the most appropriate estimated
variable value found in the LUT. The hybrid models rely on machine learning algorithms
(MLRA) to train a model on the LUT before being applied on the remote sensing images to
derive inversion maps of vegetation traits. In average, they are increasingly used, since
they exploit the potential of physically based methods combined with the flexibility and
computational efficiency of MLRA [23]. Moreover, predicting variables using lookup-table-based approaches may be less accurate, as different combinations of the input parameters can produce the same simulated spectra [24], especially in urban environments that contain many light interactions and a variety of objects and materials.

Hybrid inversion methods are well-suited for complex systems or scenarios in which traditional inversion methods may be insufficient. They can represent non-linear relationships, multi-scale processes, and interactions among components more effectively, making them suitable for studying intricate environmental phenomena such as those in urban contexts. By using hybrid inversion models with S2 data, many studies have succeeded in estimating vegetation traits, such as LCC in open-canopy conifer forests [25], LCC and LAI in deciduous broadleaf forests [26], and LCC, LAI, and CCC in mixed mountain forests [27]. In the same way, hybrid methods accurately estimated LAI and LCC in cropland [28,29].

These studies used the PROSPECT model [30,31] to simulate leaf reflectance and either SAIL [32] or INFORM [33] to simulate canopy reflectance. SAIL simulates the canopy as a horizontal turbid medium, whose properties depend on a simple representation of the canopy through the LAI and the leaf angle distribution. PROSAIL, which corresponds to the combined PROSPECT and SAIL models [34], has been extensively used for LAI and LCC retrieval of crops [35,36] and forests [37] due to its relative simplicity and accuracy. However, PROSAIL cannot accurately simulate the interactions occurring in complex canopies and environments in three dimensions such as urban canopies. INFORM can simulate more complex canopy structures [26] but cannot represent 3D scenes analogous to an urban environment. Choosing an RTM adapted to the context and able to model the complexity of interactions among the multiple scene components is essential for estimating vegetation traits in urban environments. The discrete anisotropic radiative transfer model (DART) [38,39] has many advantages for modelling three-dimensional (3D) urban scenes [40]. For example, DART has been used with S2 to retrieve the albedo of urban canopies [41] and the LCC and LAI of the sparse canopy of olive trees [42]. The latter study has similarities with this study, since the canopy of sparse olive plants is also heterogenous with a high variability in the optical properties of the ground, leading to a high proportion of S2 mixed pixels.

To our knowledge, no study has used S2 images to estimate vegetation traits of trees in urban areas. DART can simulate an urban scene, including the diversity of materials, the layout of buildings, the structure of the canopy, and its optical properties. These elements are essential for modelling a pseudo-realistic scene that reflects the heterogeneity of mixed pixels, which are common in urban remote sensing. In this study, we used a global approach to model urban trees that covers many urban configurations based on the definition of local climate zones (LCZ). The LCZ classification system provides a characterization of the physical nature of urban morphology, with LCZs defined as areas of uniform surface cover, structure, materials, and human activity that span hundreds of meters to several kilometers [43]. Originally, this classification system was dedicated to the study of urban climate and heat islands, but it is increasingly used as a reference framework for studying urban vegetation [44].

Among the MLRAs used with hybrid inversion approaches, random forest regression (RFR) [45] and Gaussian process regression (GPR) [46] are widely used. A study that compared MLRAs used with hybrid inversion showed that RFR and GPR slightly outperformed artificial neural network and support vector machines [27]. In the hybrid inversion studies cited, MLRAs were always trained with spectral features or derived features (vegetation indices). From images simulated by DART, pixel-level features can be estimated, such as the proportion of shadow, underlying vegetation, or canopy cover. These features can be used as training features and can improve model performance, as they interact with spectral features.
Finally, an important point when estimating vegetation traits using remote sensing is the approach used to combine field and satellite observations. Because field datasets are essential for validating predictions, in situ measurement protocols must be adapted to the spatial resolution of the images to ensure spatial consistency between the two data sources. For heterogeneous canopies that cover a large area, the most common method is to perform in situ measurements at the plot scale and extract the pixel values that correspond spatially to the plot. In this study, we used two spatial allocation (SA) methods, one at tree-scale and the other at pixel-scale. The tree-scale method simply assigns the vegetation traits of the dominant tree (i.e., the one with the largest crown area in the pixel) to the pixel. The pixel-scale method is more comprehensive and considered the percentage of tree(s) at the pixel scale by weighting the vegetation traits of the tree(s) by their crown area in the pixel. The two methods were evaluated and used to build the LUTs (i.e., to match the DART input parameters (LAI and LCC) and simulated pixel values) and for the validation step (i.e., to match in situ measurements and real S2 pixels). The aim of the study was to evaluate the performance of hybrid inversion models for estimating LCC, LAI, and CCC of urban trees from S2 images. Specifically, three components were investigated to determine the best strategy for estimating vegetation traits: (i) the contribution of environmental features (EFs) to training models, (ii) the spatial allocation (SA) method used to match vegetation traits and pixel values, and (iii) the MLRA used (RFR or GPR) to estimate vegetation traits.

## 2. Materials and Methods

### 2.1. *Study Sites*

The study sites were located in Rennes, in northwestern France (48°10′N, 1°68′W) (Figure 1a,b). Rennes is a medium-sized city of 222,485 inhabitants, with a population density of 4414 people/km² (INSEE, 2020) [47]. Rennes is located in an oceanic climate zone, with a projected increase in mean annual temperature from 12.4 °C to 12.9 °C and 13.2 °C for representative concentration pathway scenarios 4.5 and 8.5, respectively, from 1991–2020 to 2031–2060 (Haut Conseil pour le Climat en Bretagne 2023) [48]. Rennes managed 130,000 trees in 2023 in public areas, including trees in parks, ornamental trees, and street trees. This study focused on four of the five most common deciduous genera in Rennes—oak (*Quercus*), maple (*Acer*), plane (*Platanus*), and ash (*Fraxinus*)—as deciduous trees are the dominant trees in Rennes. We chose to focus on trees with a crown diameter of at least 10 m, to ensure that the proportion of tree canopy in the Sentinel-2 pixel was sufficiently high. Therefore, we excluded linden trees (Tilia), the fifth most common deciduous genus, which is often pruned into a rectangular shape, which restricts the extent of the crown as seen from the zenith. We selected the most common species for each genus: *Platanus acerifolia* (PL), *Quercus rubra* (QR), *Acer platanoides* (AC), and *Fraxinus excelsior* (FR). A total of 117 trees were monitored in four monospecific alignments composed of 29 PL, 29 AC, 30 QR, and 29 FR (Figure 1c–f).
<figure><img src="image_5.png" /><figcaption>**Figure 1.** (a) Location of the study area in northwestern France (red point); (b) Location of the four study sites of tree alignments in the city of Rennes (source: Rennes Métropole's 2021 orthophoto-graph); (c) *Quercus rubra*; (d) *Platanus acerifolia*; (e) *Acer platanoides*; (f) *Fraxinus excelsior* (background: Google Map 3D). Tree-row boundaries in (c–f) are displayed with red lines.</figcaption></figure>

## 2.2. Methodological Framework

Tree vegetation traits were estimated from S2 images in four steps (Figure 2): (1) collection of real data to create a validation dataset that included field data, S2 images, and ancillary data; (2) simulation of S2 images based on a design of experiment and DART modelling/simulations; (3) configuration of the real and simulated datasets, which included extraction of spectral features and EFs and matching of field and satellite data; and
(4) estimation of vegetation traits using RFR and GPR, which included feature selection and hyper-parametrization, model training, and model validation.

<figure><img src="image_6.png" /><figcaption>**Figure 2.** Flowchart developed to estimate vegetation traits. LCC: leaf chlorophyll content; LAI: leaf area index; CCC: canopy chlorophyll content.</figcaption></figure>
## 2.3. Real Data

### 2.3.1. Field Data

LCC and leaf area density (LAD) were measured in situ for the 117 trees on eight dates during the growing season (Table 1) according to the following protocol:

* LCC was measured using a Dualex leaf-clip (FORCE A, Orsay, France). Two leaves were collected in each cardinal direction on the edge of the crown and as high up as possible using a lopper. Two Dualex measurements were taken per leaf. Mean LCC per tree, which equaled the mean of all 16 Dualex readings on a given date, was calculated for the eight dates. We used the Dualex device rather than the widely used SPAD and CCM-200 chlorophyll meters, as it responds linearly to increasing chlorophyll content rather than non-linearly. An equation developed for dicot species was used to retrieve LCC from the Dualex reading [49]:

$$LCC = -4.84 + 1.24 \cdot \text{Dualex} \quad (1)$$

* LAD was measured using a canopy analyzer (LAI-2200, LiCor, Lincoln, NE, USA). The measurement protocol was adapted according to the user manual and the protocol of Wei et al. [50].

<table><caption>**Table 1.** Dates of field measurements and of Sentinel-2 images and presence (P) or absence (A) of corresponding leaf canopy cover (LCC) and leaf area density (LAD) measurements. AC: *Acer platanoides*; FR: *Fraxinus excelsior*; PL: *Platanus acerifolia*; QR: *Quercus rubra*.</caption><thead><tr><th colspan="2">Dates</th><th rowspan="3">Image Lag Time<br>(Days)</th><th colspan="6">LCC Measurement</th><th colspan="6">LAD Measurement</th></tr><tr><th rowspan="2">Field</th><th rowspan="2">Sentinel-2</th><th colspan="4">Species</th><th colspan="4">Species</th></tr><tr><th>AC</th><th>FR</th><th>PL</th><th>QR</th><th>AC</th><th>FR</th><th>PL</th><th>QR</th></tr></thead><tbody><tr><td>27 April 2021</td><td>23 April 2021</td><td>-4</td><td>A</td><td>A</td><td>A</td><td>A</td><td>P</td><td>A</td><td>P</td><td>A</td><td>P</td><td>P</td></tr><tr><td>11 May 2021</td><td>6 May 2021</td><td>-5</td><td>P</td><td>P</td><td>P</td><td>P</td><td>P</td><td>P</td><td>P</td><td>P</td><td>P</td><td>P</td></tr><tr><td>2 June 2021</td><td>31 May 2021</td><td>-2</td><td>P</td><td>P</td><td>P</td><td>P</td><td>P</td><td>P</td><td>P</td><td>P</td><td>P</td><td>P</td></tr><tr><td>23 June 2021</td><td>15 June 2021</td><td>-8</td><td>P</td><td>P</td><td>P</td><td>P</td><td>P</td><td>P</td><td>P</td><td>P</td><td>P</td><td>P</td></tr><tr><td>21 July 2021</td><td>20 July 2021</td><td>-1</td><td>P</td><td>P</td><td>P</td><td>P</td><td>P</td><td>P</td><td>P</td><td>P</td><td>P</td><td>P</td></tr><tr><td>17 August 2021</td><td>14 August 2021</td><td>-3</td><td>P</td><td>P</td><td>P</td><td>P</td><td>P</td><td>P</td><td>P</td><td>P</td><td>P</td><td>P</td></tr><tr><td>1 September 2021</td><td>5 September 2021</td><td>4</td><td>P</td><td>P</td><td>P</td><td>P</td><td>P</td><td>P</td><td>P</td><td>P</td><td>P</td><td>P</td></tr><tr><td>20 September 2021</td><td>13 September 2021</td><td>-7</td><td>P</td><td>P</td><td>P</td><td>P</td><td>A</td><td>P</td><td>P</td><td>P</td><td>P</td><td>P</td></tr></tbody></table>

Data are available at https://zenodo.org/records/12751353 (accessed on 17 July 2024); field survey and protocols for measuring LCC and LAD are described by Le Saint et al. 2024 [51].

### 2.3.2. Sentinel-2 Data

The eight cloud-free S2 L2A images acquired closest to the field measurement dates were selected (Table 1) and downloaded from the Copernicus Browser [52]. The mean and longest time lags between the dates of field measurement and S2 images were 4.2 and 8 days, respectively. These images, which were already corrected for atmospheric effects using the Sen2Cor algorithm, were co-registered with a precision of < 0.8 pixels using the COREGIS method [53]. Only the bands at 10 and 20 m spatial resolution were retained, and the 20 m bands were resampled to 10 m using the nearest-neighbour method.

### 2.3.3. Ancillary Data

The ancillary data included several open-source layers used to derive EFs (Table 2). A digital terrain model (DTM) and a digital surface model (DSM) were used to derive the proportion of shadow (Pshadow). They were also used to calculate LAI by providing tree-crown volume, whereas tree-crown extent was delineated from orthophotographs.
The tree-crown extent was also used to derive canopy cover (CCP) at the pixel scale. The grass-extent layer was used to derive the grass proportion (Pgrass) at the pixel scale.

<table><caption>**Table 2.** Description, resolution, and altimetric precision (for raster data), and sources of ancillary data.</caption><thead><tr><th>Name</th><th>Description</th><th>Resolution/Altimetric Precision</th><th>Source</th></tr></thead><tbody><tr><td>Digital terrain model</td><td>Spatial raster representing the elevation of the surface of bare Earth, free of natural and built features</td><td>0.5 m/0.2 m</td><td>Opendata Rennes Métropole</td></tr><tr><td>Digital surface model</td><td>Spatial raster representing the elevation of the surface of bare Earth, including natural and built features</td><td>0.5 m/0.2 m</td><td>Opendata Rennes Métropole</td></tr><tr><td>Orthophotographs</td><td>Optical visible orthophotographs of Rennes Métropole acquired in 2021</td><td>0.05 m/-</td><td>Opendata Rennes Métropole</td></tr><tr><td>Tree-crown extent</td><td>Spatial vector of tree-crown extent</td><td>-</td><td>Manual digitization of orthophotographs</td></tr><tr><td>Grass extent</td><td>Spatial vector of grass extent</td><td>-</td><td>OpenStreetMap ("landuse" key and "grass" value)</td></tr></tbody></table>

## 2.4. Simulated Data

### 2.4.1. 3D Urban Scene Modelling

The 3D urban scenes including trees were defined using the local climate zone (LCZ) taxonomy [43,54]. Four LCZ types were selected from the most dominant LCZ built-cover types encountered in European cities with more than 100,000 inhabitants [55]: LCZ 2 (Compact Midrise), LCZ 5 (Open Midrise), LCZ 6 (Open Low-rise), and LCZ 8 (Large Low-rise). Three LCZ parameters were used to build contrasting urban scenes: height of roughness elements, which is the mean height of the buildings; aspect ratio, calculated as the height of roughness elements divided by the width of the urban canyon; and fraction of area in buildings, which is the percentage of the ground surface occupied by buildings. Their values (Table 3) did not vary during the simulations.

<table><caption>**Table 3.** Local climate zone (LCZ) parameters used for Compact Midrise (LCZ 2), Open Midrise (LCZ 5), Open Low-rise (LCZ 6), and Large Low-rise (LCZ 8). The height of roughness elements is the mean height of the buildings, the aspect ratio is calculated as the height of roughness elements divided by the width of the urban canyon, and the fraction of area in buildings is the percentage of the ground surface occupied by buildings.</caption><thead><tr><th>Parameter</th><th>LCZ 2</th><th>LCZ 5</th><th>LCZ 6</th><th>LCZ 8</th></tr></thead><tbody><tr><td>Height of roughness elements [m]</td><td>18.0</td><td>18.0</td><td>6.5</td><td>6.5</td></tr><tr><td>Aspect ratio</td><td>1.375</td><td>0.525</td><td>0.525</td><td>0.200</td></tr><tr><td>Fraction of area in buildings [%]</td><td>55</td><td>30</td><td>30</td><td>40</td></tr></tbody></table>

Only deciduous tree species were modelled. We defined a 3D tree model (reference tree) with a tree height of 15 m, crown height of 11 m, and trunk diameter of 0.4 m, based on statistical analysis of the urban tree database of the Lyon Metropolitan Area (France) (https://data.grandlyon.com/jeux-de-donnees/arbres-alignement-metropole-lyon/info, accessed on 17 July 2024). The 3D tree was modelled with an ellipsoid crown and a cylindrical trunk with branches. This simplistic representation was already used for other studies in urban and forest contexts [56]. To consider alignment heterogeneity, we created an alignment of five trees that contained two profiles of tree functional traits. Specifically, the central tree and the two trees at the ends of the alignment had one profile of optical and biophysical traits (profile A), while the other two trees had the other profile
(profile B) (Figures 3 and 4). The 3D urban scenes were generated using R, and the final 3D objects (.obj format) were then input into DART.

<figure><img src="image_7.png" /><figcaption>**Figure 3.** Three-dimensional zenith and perspective views of urban scenes defined for the four local climate zones (LCZs): Compact Midrise (LCZ 2), Open Midrise (LCZ 5), Open Low-rise (LCZ 6), and Large Low-rise (LCZ 8).</figcaption></figure>

<figure><img src="image_8.png" /><figcaption>**Figure 4.** Illustration of the four extraction windows used to extract spectral features at 10 m resolution and the corresponding crown area. Delta (Δ) represents the offset between the window centre and the tree centroid. Theta (θ) represents the offset angle between the window centre and the alignment axis relative to the tree centroid. CCP represents the resulting canopy cover percentage when tree diameter was set to 10 m. (A) and (B) correspond to the two different tree profiles in the alignment.</figcaption></figure>

### 2.4.2. DART Parametrization and Simulation

Although DART includes dozens of tunable parameters, 32 parameters were used to simulate S2 images (Table 4). See Appendix A for details of these parameters. They included 12 tree-endogenous and 20 tree-exogenous parameters that have direct or indirect influence, respectively, on simulated tree canopy reflectance. Whereas some parameters (15)
were fixed, others (17) varied according to predefined values/ranges and distributions. For the purpose of generalization, the min and max bounds for the leaf traits were determined with reference to the literature [57,58].

<table><caption>**Table 4.** DART parametrization, describing the section of the DART graphical user interface and the parameter's name, category, type ((F)ixed or (V)ariable), and value (if F) or range (if V).</caption><thead><tr><th>DART Section</th><th>Parameter Name</th><th>Category</th><th>Type</th><th>Values and Range</th></tr></thead><tbody><tr><td>Global settings</td><td>Light propagation mode</td><td>exogenous</td><td>F</td><td>Bi-directional (DART-Lux)</td></tr><tr><td rowspan="4">Sensor settings</td><td>Spectral bands</td><td>exogenous</td><td>F</td><td>According to Sentinel-2 sensor</td></tr><tr><td>Zenithal angle</td><td>exogenous</td><td>F</td><td>2.8 [°]</td></tr><tr><td>Azimuth angle</td><td>exogenous</td><td>F</td><td>182 [°]</td></tr><tr><td>Spatial resolution</td><td>exogenous</td><td>F</td><td>1 [m]</td></tr><tr><td rowspan="3">Direction input parameter</td><td>Hour</td><td>exogenous</td><td>F</td><td>11:07 UTC</td></tr><tr><td>Day</td><td>exogenous</td><td>F</td><td>Day 15 of each month</td></tr><tr><td>Month</td><td>exogenous</td><td>V</td><td>From March to November</td></tr><tr><td rowspan="2">Atmosphere</td><td>Atmosphere model</td><td>exogenous</td><td>F</td><td>USSTD76</td></tr><tr><td>Aerosol properties</td><td>exogenous</td><td>F</td><td>Urban Type<br>Aerosol optical depth = 1</td></tr><tr><td rowspan="4">Scene optical properties</td><td>Roof</td><td>exogenous</td><td>V</td><td>See Appendix A</td></tr><tr><td>Wall</td><td>exogenous</td><td>V</td><td>See Appendix A</td></tr><tr><td>Impervious ground</td><td>exogenous</td><td>V</td><td>See Appendix A</td></tr><tr><td>Pervious ground</td><td>exogenous</td><td>V</td><td>See Appendix A</td></tr><tr><td rowspan="3">Earth scene</td><td>Dimensions</td><td>exogenous</td><td>F</td><td>100 m × 100 [m]</td></tr><tr><td>Latitude</td><td>exogenous</td><td>F</td><td>48.1° N</td></tr><tr><td>Longitude</td><td>exogenous</td><td>F</td><td>1.68° W</td></tr><tr><td rowspan="4">Tree planting conditions</td><td>Distance to nearest building</td><td>exogenous</td><td>V</td><td>LCZ2 and LCZ6: 5.0–6.5 [m]<br>LCZ5 and LCZ8: 6–16 [m]</td></tr><tr><td>Tree exposure</td><td>exogenous</td><td>V</td><td>Shady side or sunny side</td></tr><tr><td>Street orientation <sup>1</sup></td><td>exogenous</td><td>V</td><td>0, 45, 90 and 135 [°]</td></tr><tr><td>Percentage of grass on the ground <sup>2</sup></td><td>exogenous</td><td>V</td><td>0–100 [%]</td></tr><tr><td rowspan="4">Tree</td><td>Tree-crown diameter</td><td>endogenous</td><td>V</td><td>10 and 12 [m]</td></tr><tr><td>Other geometric parameters</td><td>endogenous</td><td>F</td><td>See Appendix A</td></tr><tr><td>Leaf angle distribution</td><td>endogenous</td><td>V</td><td>Plagiophile and planophile</td></tr><tr><td>Leaf area density (LAD)</td><td>endogenous</td><td>V</td><td>0.1 and 1.2 [m<sup>2</sup>/m<sup>3</sup>]</td></tr><tr><td rowspan="8">Leaf</td><td>Clumping factor</td><td>endogenous</td><td>V</td><td>0–50 [%]</td></tr><tr><td>Structure coefficient (N)</td><td>endogenous</td><td>V</td><td>1.1–2.3 [arbitrary unit]</td></tr><tr><td>Leaf chlorophylls content (Cab)</td><td>endogenous</td><td>V</td><td>5–60 [µg/cm<sup>2</sup>]</td></tr><tr><td>Carotenoid content (Car)</td><td>endogenous</td><td>V</td><td>2.5–25 [µg/cm<sup>2</sup>]</td></tr><tr><td>Brown pigment</td><td>endogenous</td><td>F</td><td>0 [arbitrary unit]</td></tr><tr><td>Anthocyanin</td><td>endogenous</td><td>F</td><td>0 [µg/cm<sup>2</sup>]</td></tr><tr><td>Equivalent water thickness</td><td>endogenous</td><td>V</td><td>0.004–0.024 [cm]</td></tr><tr><td>Dry matter content</td><td>endogenous</td><td>V</td><td>0.002–0.014 [g/cm<sup>2</sup>]</td></tr></tbody></table>

DART inputs were created from the variable parameters (Table 4) using a design of experiment (DOE) computational approach based on a Latin hypercube sampling (Python v3.9 library OpenTurns v1.34 [59]). A total of 80,000 simulations were performed (i.e., 20,000 per LCZ). Thus, 80,000 coordinate points [X1, ..., Xp] were generated, where p is the number of parameters, to cover the p-dimensional space as well as possible. For continuous variable parameters (type = V in Table 4), values were generated using a uniform distribution.

<footnote n="1">¹ The street orientation was set at an anti-clockwise angle, with the following orientations: 0°: west/east, 45°: south-west/north-east, 90°: north/south, and 135° south-east/north-west. ² The area not covered by grass was assumed to be mineral.</footnote>
DART was configured to generate orthorectified images with a size of 100 m × 100 m, and a spatial resolution of 1 m. The images were generated according to the S2 spectral and geometrical (view angles) characteristics. Specifically, only the 10 m (VIS and PIR) and 20 m (RE, SWIR) bands were simulated, and central wavelengths and bandwidths were defined based on characteristics of the S2 MSI sensor. Zenithal and azimuthal angles were set to the mean S2 view angles: 3° (with zenith = 0°) and 180° (with north = 0°), respectively. As the synthetic S2 images were initially simulated at 1 m resolution, we used a method based on extraction windows to extract pixel spectral values, considering the 10 m and 20 m spatial resolutions of S2 (see Section 2.4.3). DART version 5.8.10v1259 was used.

### 2.4.3. Spatial Window Extraction

The simulated images were generated at a spatial resolution of 1 m, which represented interactions among elements of the 3D scenes more precisely. We then aggregated spectral values at 10 m (B02, B03, B04, and B08) and 20 m resolution (B05, B06, B07, B8A, B11, and B12) using four pair-extraction windows of 10 × 10 and 20 × 20 m (Figure 4): one tree-centred window and three windows shifted by 5 m at 0°, 45°, and 90° with respect to the centre of the tree-canopy in the direction of the tree alignment. The 20 × 20 m windows were configured so that their lower left quarter corresponded to the associated 10 × 10 m window. This extraction-window method was applied to consider S2 pixel heterogeneity (mixed pixels), geolocation uncertainty, and urban conditions. This approach allowed a wide range of pixel compositions to be considered, with different proportions of canopy and shadow in the pixel, as well as pixels that overlapped two trees. In terms of computational efficiency, it allowed a diverse set of pixels to be obtained from a single simulated image, while decreasing the total number of simulations and simplifying parametrization.

## 2.5. Dataset Configuration

### 2.5.1. Environmental Feature Extraction

Contextual information at the pixel scale was included as three EFs in the RFR and GPR models:

1. Green proportion (Pgreen): the percentage of the pixel covered by underlying vegetation, which strongly influences reflectance when the tree canopy does not cover the entire pixel and/or the tree has low LAI;

2. Shadow proportion (Pshadow): the percentage of the pixel covered by shadow, which also strongly influences reflectance, particularly of one near-infrared S2 band (B08) with solar angle decrease;

3. Canopy cover pixel (CCP): the percentage of the pixel covered by tree canopy, which indicates the percentage of pixel purity.

For the simulated dataset (SDS), Pgreen had already been defined as an input parameter for DART (Table 4) and thus was known at the pixel scale for each simulated image. Pshadow was calculated by DART at the pixel scale for each simulated image using a direct insolation mask. CCP was calculated at the pixel scale by intersecting the four extraction windows previously defined and the crown extent simulated by DART (Figure 4). Depending on the window and crown diameter, CCP ranged from 39% (window 4 with a 10 m-diameter crown) to 95% (window 1 with a 12 m-diameter crown).

For the real dataset (RDS), EFs were calculated using different methods than those used for the SDS. First, we considered a spatial reference grid (S2 grid) that corresponded to S2 images at 10 m resolution, then EFs were calculated as follows:

1. Pgreen was calculated at the pixel scale by intersecting the S2 grid and the grass-extent layer;

2. Pshadow was calculated using a raster layer of potential direct incoming solar radiation (kWh/m²) at 1 m resolution derived from the DSM using the Potential Incoming Solar Radiation algorithm (Terrain Analysis; Lighting, Visibility) in SAGA software v7.8.2 [60]. Pixels with > 0 kWh/m² were classified as non-shadow (0), while those
with 0 kWh/m² were classified as shadow. This layer was then aggregated to the S2 grid resolution (10 m);

3. For CCP calculation, the S2 grid and the tree-crown extent were spatially intersected. Finally, CCP was calculated at the pixel scale based on this spatial intersection.

### 2.5.2. Spectral Feature Extraction

Spectral features were extracted on specific spatial locations depending on the dataset. For the RDS, pixels with CCP < 35% were removed, as this value was close to the minimum CCP of the SDS (39%); the remaining pixels were considered pixels of interest. Spectral features were then extracted from the pixels of interest for the RDS and in the windows for SDS, and 18 vegetation indices correlated with LCC and LAI according to the literature were calculated for both datasets (Table 5).

<table><caption>**Table 5.** Vegetation indices used in the study. Abbrev.: abbreviation, Ref.: reference.</caption><thead><tr><th>Index</th><th>Abbrev.</th><th>Equation (with S2 Band Names)</th><th>Ref.</th></tr></thead><tbody><tr><td>Red-green-blue vegetation index</td><td>RGBVI</td><td>RGBVI = <sup>B03<sup>2</sup>-(B02-B04)</sup>&frasl;<sub>B03<sup>2</sup>+(B02-B04)</sub></td><td>[61]</td></tr><tr><td>Green leaf index</td><td>GLI</td><td>GLI = <sup>2·B03-B04-B02</sup>&frasl;<sub>2·B03+B04+B02</sub></td><td>[62]</td></tr><tr><td>Normalized green-blue difference index</td><td>NGBDI</td><td>NGBDI = <sup>B02-B03</sup>&frasl;<sub>B02+B03</sub></td><td>[63]</td></tr><tr><td>Structure insensitive pigment index</td><td>SIPI</td><td>SIPI = <sup>B08-B02</sup>&frasl;<sub>B08-B04</sub></td><td>[64,65]</td></tr><tr><td>Normalized difference vegetation index</td><td>NDVI</td><td>NDVI = <sup>B08-B04</sup>&frasl;<sub>B08+B04</sub></td><td>[66]</td></tr><tr><td>Atmospherically resistant vegetation index</td><td>ARVI</td><td>ARVI = (<sup>B08-(B04-1*(B02-B04))</sup>&frasl;<sub>B08+(B04-1*(B02-B04))</sub>)</td><td>[67]</td></tr><tr><td>Enhanced vegetation index</td><td>EVI</td><td>EVI = 2.5·<sup>B08-B04</sup>&frasl;<sub>B08+B04</sub></td><td>[68,69]</td></tr><tr><td>Optimized soil adjusted vegetation index</td><td>OSAVI</td><td>OSAVI = (1 + 1.16)·(B08 - B04)/(B08 + B04 + 0.16)</td><td>[70]</td></tr><tr><td>Modified chlorophyll absorption in reflectance index 2</td><td>MCARI2</td><td>MCARI2 → See reference</td><td>[71]</td></tr><tr><td>Red-edge normalized difference vegetation index</td><td>NDVIRE</td><td>NDVIRE = <sup>B08-B05</sup>&frasl;<sub>B08+B05</sub></td><td>[72]</td></tr><tr><td>Sentinel-2 LAI index</td><td>SELI</td><td>SELI = <sup>B8A-B05</sup>&frasl;<sub>B8A+B05</sub></td><td>[73]</td></tr><tr><td>Mixed leaf area index vegetation index</td><td>MixLAIVI</td><td>MixLAIVI = <sup>B04+B11<sup>2</sup></sup>&frasl;<sub>2</sub></td><td>[74]</td></tr><tr><td>Transformed chlorophyll absorption in reflectance index</td><td>TCARI</td><td>TCARI = 3·[(B05 - B04) - 0.2·(B05 - B03)·(B05/B04)]</td><td>[75]</td></tr><tr><td>TCARI/OSAVI</td><td>CCII</td><td>CCII = <sup>TCARI</sup>&frasl;<sub>OSAVI</sub></td><td>[71,75,76]</td></tr><tr><td>Sentinel-2-based triangular vegetation index</td><td>STVI</td><td>STVI → See equation in [77]</td><td>[77]</td></tr><tr><td>Greenness component of Sentinel-2 tasseled cap transformation</td><td>TCT_G</td><td>TCT_G → See equation in [78]</td><td>[78]</td></tr><tr><td>Brightness component of Sentinel-2 tasseled cap transformation</td><td>TCT_B</td><td>TCT_B → See equation in [78]</td><td>[78]</td></tr><tr><td>Wetness component of Sentinel-2 tasseled cap transformation</td><td>TCT_W</td><td>TCT_W → See equation in [78]</td><td>[78]</td></tr></tbody></table>

### 2.5.3. Spatial Allocation for Field-Satellite Data Matching

To match field and satellite data, we applied two spatial allocation (SA) methods that intersected pixels and crown areas: one at the tree scale and one at the pixel scale. The former simply allocated the field measurements of LCC, LAI, and CCC of the dominant tree in a given pixel to that pixel (LAItree, LCCtree, and CCCtree, respectively). The pixel-scale method considered the percentage of trees at the pixel scale by weighting the field measurements of LCC, LAI, and CCC of the tree(s) by their crown area in the pixel (LAIpix, LCCpix, and CCCpix, respectively).
Both methods were applied to the SDS and RDS. Before doing so, tree LAD (m²/m³) was converted into LAI (m²/m²), as the in situ LAI-2200 measurements and DART input values for a given tree were LAD. This step required knowing geometric properties of the trees, which were derived from the auxiliary data for real trees or from DART input parameters for simulated trees (Table 6). This step was based on two simplifying assumptions: (i) crown height equals 2/3 of the tree height (i.e., live crown ratio [79–81], and (ii) the tree crown can be considered an ellipsoid, which can be used to calculate its volume.

<table><caption>**Table 6.** Tree geometric properties and equations used to calculate leaf area index (LAI) from leaf area density (LAD).</caption><thead><tr><th>Property</th><th>Name</th><th>Unit</th><th>Equation</th><th>Description</th></tr></thead><tbody><tr><td>Tree height</td><td>H<sub>tree</sub></td><td>m</td><td>-</td><td>Tree height, given by DSM-DTM, according to the crown centroid or corresponding to DART input</td></tr><tr><td>Tree crown height</td><td>H<sub>crown</sub></td><td>m</td><td>H<sub>tree</sub> · <sup>2</sup>&frasl;<sub>3</sub></td><td>Crown height, assumed to equal 2/3 of the tree height</td></tr><tr><td>Tree crown ellipsoid semi-axis (b, c)</td><td>b, c</td><td>m</td><td>-</td><td>Semi-axis of the tree crown (considered as and ellipsoid). Calculated from 2D crown-delineation polygons or corresponding to DART input</td></tr><tr><td>Tree crown volume</td><td>V<sub>crown</sub></td><td>m<sup>3</sup></td><td><sup>4</sup>&frasl;<sub>3</sub>π · <sup>H<sub>crown</sub></sup>&frasl;<sub>2</sub> · b · c</td><td>Crown volume, calculated as that of an ellipsoid</td></tr><tr><td>Tree crown area</td><td>A<sub>crown</sub></td><td>m<sup>2</sup></td><td>-</td><td>Area projected onto the ground by the crown</td></tr><tr><td>Tree LAD</td><td>LAD<sub>tree</sub></td><td>m<sup>2</sup>/m<sup>3</sup></td><td>-</td><td>Leaf area density, measured with LAI-2200 or corresponding to DART input</td></tr><tr><td>Tree LAI</td><td>LAI<sub>tree</sub></td><td>m<sup>2</sup>/m<sup>2</sup></td><td>LAD<sub>tree</sub> · <sup>V<sub>crown</sub></sup>&frasl;<sub>A<sub>crown</sub></sub></td><td>Leaf area index of the tree</td></tr></tbody></table>

Both SA methods intersected crown area and the extraction window for SDS or the S2 grid for RDS to calculate five properties (Table 7): pixel area, intersection area between a given tree and a pixel (when several trees intersected a pixel, the intersection area of each tree was calculated), total canopy area of the pixel, percentage of canopy cover of the pixel, and the percentage of canopy in the pixel for a given tree.

<table><caption>**Table 7.** Properties calculated from intersecting crown area and the extraction window for the simulated dataset or the Sentinel-2 grid for the real dataset.</caption><thead><tr><th>Property</th><th>Name</th><th>Unit</th><th>Equation</th><th>Description</th></tr></thead><tbody><tr><td>Pixel area</td><td>A<sub>pix</sub></td><td>m<sup>2</sup></td><td>-</td><td>Constant area for a pixel of 10 m resolution (100 m<sup>2</sup>)</td></tr><tr><td>Intersection area tree-pixel</td><td>Ainter<sub>tree</sub></td><td>m<sup>2</sup></td><td>-</td><td>Intersection area between a given tree and a given pixel. A tree can overlap several pixels and vice-versa.</td></tr><tr><td>Total canopy area</td><td>TCA</td><td>m<sup>2</sup></td><td>Σ<sub>i=1</sub><sup>n</sup> Ainter<sub>tree<sub>i</sub></sub></td><td>Total canopy area for a given pixel</td></tr><tr><td>Canopy cover pixel</td><td>CCP</td><td>%</td><td>TCA/A<sub>pix</sub></td><td>Percentage of canopy cover in the pixel</td></tr><tr><td>Percentage of canopy cover</td><td>PCCP<sub>tree</sub></td><td>%</td><td>Ainter<sub>tree</sub>/TCA</td><td>Percentage of canopy in the pixel for a given tree</td></tr></tbody></table>

For the tree-scale method, LCCtree equalled the value measured in situ or the DART input parameter, LAI<sub>tree</sub> was calculated from LAD (Table 6), and CCC<sub>tree</sub> equalled LAI<sub>tree</sub> × LCC<sub>tree</sub> (Table 8). For the pixel-scale method, the vegetation traits for each pixel were calculated using equations (Table 8).
<table><caption>**Table 8.** Equations used to calculate vegetation traits at the tree scale and pixel scale.</caption><thead><tr><th>Trait</th><th>Name</th><th>Unit</th><th>Equation</th><th>Description</th></tr></thead><tbody><tr><td>LAI tree</td><td>LAI<sub>tree</sub></td><td>m<sup>2</sup>/m<sup>2</sup></td><td>LAD<sub>tree</sub> · <sup>V<sub>crown</sub></sup>/<sub>A<sub>crown</sub></sub></td><td>Leaf area index for the tree</td></tr><tr><td>LCC tree</td><td>LCC<sub>tree</sub></td><td>µg/cm<sup>2</sup></td><td>-</td><td>Dualex leaf-clip reading for a given tree or corresponding to DART input.</td></tr><tr><td>CCC tree</td><td>CCC<sub>tree</sub></td><td>µg/m<sup>2</sup></td><td>LAI<sub>tree</sub> · LCC<sub>tree</sub></td><td>Canopy chlorophyll content</td></tr><tr><td>LAI pixel</td><td>LAI<sub>pix</sub></td><td>m<sup>2</sup>/m<sup>2</sup></td><td>Σ<sub>i=1</sub><sup>n</sup> LAI<sub>tree<sub>i</sub></sub> · Ainter<sub>tree<sub>i</sub></sub></td><td>Weighted sum of LAI<sub>tree</sub> in the pixel. n equals the number of trees intersected with the pixel.</td></tr><tr><td>LCC pixel</td><td>LCC<sub>pix</sub></td><td>µg/cm<sup>2</sup></td><td>Σ<sub>i=1</sub><sup>n</sup> LCC<sub>tree<sub>i</sub></sub> · PCCP<sub>tree<sub>i</sub></sub></td><td>Leaf chlorophyll content at the pixel scale</td></tr><tr><td>CCC pixel</td><td>CCC<sub>pix</sub></td><td>µg/m<sup>2</sup></td><td>Σ<sub>i=1</sub><sup>n</sup> CC<sub>tree<sub>i</sub></sub> · Ainter<sub>tree<sub>i</sub></sub></td><td>Weighted sum of CCC<sub>tree</sub> in the pixel. n equals the number of trees intersected with the pixel.</td></tr></tbody></table>

## 2.6. Machine Learning Regression Algorithms: Building Strategy and Training

For each of the three vegetation traits (LCC, LAI, and CCC), we built a total of 24 regression models depending on different configurations defined by three components: (i) the type of MLRA (RFR or GPR), (ii) the exclusion or inclusion of EFs (models trained using only spectral features or combining these spectral features with the three EFs (CCP, Pgreen, and Pshadow)), and (iii) the SA method (tree scale or pixel scale).

RFR and GPR have been widely used to estimate LAI and LCC from S2 data using hybrid inversion approaches with PROSAIL [29,82], INFORM [27], and DART [42] models. Non-parametric MLRAs such as RFR and GPR provide faster and more accurate estimates than simple parametric MLRAs do [23]. Their ability to capture non-linear relationships and interactions among variables makes them suitable for simulating complex datasets. RFR is an ensemble-learning algorithm that combines the predictions of several decision trees built from bootstrap data samples. This method, which can capture complex relationships between explanatory and interaction features, addresses overfitting problems effectively while maintaining high prediction accuracy. GPR uses Bayesian inference and a pre-set covariance kernel function that is optimized to fit the data. The kernels provide the expected correlation between different observations. In this study, the SDS and RDS included spectral, derived (vegetation indices), and contextual features. We used a radial basis function as the kernel to train GPR.

To identify an optimal set of input spectral features to use to train the models, a Pearson correlation matrix was calculated, and, when two features were strongly correlated (i.e., r > 0.95), the feature with the largest mean absolute correlation was removed. This step aimed to improve model interpretability and performance by removing redundant features.

Next, for each target variable and both SA methods (i.e., LCCtree, LCCpix, LAItree, LAIpix, CCCtree, and CCCpix), recursive feature elimination (RFE) was applied, which iteratively removed less informative features by calculating variable importance based on random forests. Backward selection of features was based on the root mean square error (RMSE), which was calculated using a cross-validation procedure with 10 folds and two repeats. RFE was applied to each target variable, and the feature set with the lowest RMSE (Figure 5) was used for model training. For RFR models, hyperparameter tuning was performed using a grid-search approach with an RMSE cost function. Two parameters were used for the grid search: the number of trees {100, 200, ..., 1000} and the number of random features selected at each division {n/4, n/3, n/2} (where n is the number of training features). The parameter combination with the lowest RMSE was used for model training. Overall, 24 models were trained and evaluated (three target variables × eight configurations per target variable).
<figure><img src="image_9.png" /><figcaption>**Figure 5. Ranking of variable importance according to RFE for each vegetation trait and spatial allocation. LCC: Leaf chlorophyll content; LAI: Leaf area index; CCC: Canopy chlorophyll content.**</figcaption></figure>

## 2.7. Validation

All models were cross-validated using a k-fold strategy, with k = 10. Once models were trained with the SDS, they were applied to the RDS and validated using field measurements in the RDS. The minimum, mean, and maximum values of the vegetation traits included the RDS, which were determined based on the two SA methods and in situ measurements, are presented in Table 9.
<table><caption>**Table 9.** Minimum, mean, and maximum values for the three vegetation traits (leaf chlorophyll content, leaf area index, and canopy chlorophyll content) based on the two spatial allocation methods (at tree and pixel scales) and in situ measurements. Min: minimum; Max: Maximum.</caption><thead><tr><th rowspan="2">Vegetation Traits</th><th colspan="3">Tree-Scale</th><th colspan="3">Pixel-Scale</th></tr><tr><th>Min</th><th>Mean</th><th>Max</th><th>Min</th><th>Mean</th><th>Max</th></tr></thead><tbody><tr><td>Leaf chlorophyll content [µg/cm²]</td><td>11.4</td><td>28.9</td><td>49.3</td><td>11.6</td><td>28.9</td><td>48.3</td></tr><tr><td>Leaf area index [m²/m²]</td><td>0.94</td><td>3.78</td><td>8.24</td><td>0.01</td><td>2.36</td><td>7.92</td></tr><tr><td>Canopy chlorophyll content [µg/m²]</td><td>11</td><td>111</td><td>279</td><td>1</td><td>71</td><td>289</td></tr></tbody></table>

Five metrics were used to assess the performance of the regression models:

* Coefficient of determination (R²): The coefficient of determination measures the proportion of the variance in the dependent variable that is explained by the regression model. It ranges from 0 to 1, where 1 indicates a perfect fit of the model to the observed data.

$$R^2(\text{Obs, Pred}) = 1 - \frac{\sum_{i=1}^{n} (\text{Pred}_i - \text{Obs}_i)^2}{\sum_{i=1}^{n} (\text{Pred}_i - \overline{\text{Pred}_i})^2} \quad (2)$$

* Root mean squared error (RMSE): RMSE is a measure of model accuracy that calculates the square root of the mean of the squares of the differences between predictions and observations, thus indicating the mean deviation between them.

$$\text{RMSE}(\text{Obs, Pred}) = \sqrt{\frac{1}{n} \sum_{i=1}^{n} (\text{Pred}_i - \text{Obs}_i)^2} \quad (3)$$

* Symmetric mean absolute percentage error (SMAPE): SMAPE equals the percentage difference between predictions and observations while accounting for their scales. In this study, it provided an intuitive interpretation and allowed model performances for the three vegetation traits to be compared.

$$\text{SMAPE}(\text{Obs, Pred}) = \frac{100}{n} \sum_{i=1}^{n} \frac{|\text{Pred}_i - \text{Obs}_i|}{(|\text{Obs}_i| + |\text{Pred}_i|) \cdot 0.5} \quad (4)$$

* Bias: Bias equals the mean difference between predictions and observations, which indicates a model's tendency to overestimate or underestimate the dependent variable.

$$\text{Bias}(\text{Obs, Pred}) = \frac{1}{n} \sum_{i=1}^{n} (\text{Obs}_i - \text{Pred}_i) \quad (5)$$

* Bias standard deviation (BSD): The standard deviation of the bias measures the distribution of prediction errors around the mean bias, which indicates the variability in differences between predictions and observations.

$$\text{BSD}(\text{Obs, Pred}) = \sqrt{\frac{1}{n} \sum_{i=1}^{n} [(\text{Obs}_i - \text{Pred}_i) - \text{Bias}]^2} \quad (6)$$

For Equations (2)–(6), Obs is the observed value of the target variable, Pred is the predicted value of the target variable, and n is the total number of samples in the RDS.

Finally, temporal metrics were calculated to assess the temporal consistency between estimated and observed values. Temporal consistency ensures that observed changes or trends in predicted vegetation traits reflect real phenomena rather than artefacts introduced by inconsistencies in data acquisition or processing [83]. This is particularly important in the context of vegetation monitoring, when trends and/or temporal anomalies derived
from image time-series are most often more relevant than information obtained from single-date observations [22,84].

The LCC, LAI, and CCC time series were estimated for each pixel using the best models (according to the metrics used). For each pixel, the estimated series were compared to the observed series. Three temporal similarity metrics available in the R package *TSclust* [85] were used:

* $d_{EUCL}$, based on the Euclidean distance between values observed at the same points in time

* $d_{COR}$, based on Pearson correlation between the two series [86]

* $d_{CORT}$, based on temporal correlation between the two series, used to include both conventional measures for the proximity of observations and temporal correlation to estimate the proximity behaviour and dynamic [87].

## 3. Results

### 3.1. Accuracy of LCC, LAI, and CCC Estimation Using the Simulated Dataset

Values of metrics calculated for all configurations are shown in Table 10. Cross-validation scatterplots are shown in Appendix B.

<table><caption>**Table 10.** Performance metrics of leaf chlorophyll content (LCC), leaf area index (LAI), and canopy chlorophyll content (CCC) obtained by cross validation using the simulated dataset. SA: spatial allocation; MLRA: machine learning regression algorithm, RMSE: root mean square error; SMAPE: symmetric mean absolute percentage error, BSD: bias standard deviation; RMSE is expressed in µg/cm² for LCC and CCC, m²/m² for LAI.</caption><thead><tr><th>Target Variable</th><th>SA Method</th><th>Environment Features</th><th>MLRA</th><th>R<sup>2</sup></th><th>RMSE</th><th>SMAPE</th><th>BIAS</th><th>BSD</th></tr></thead><tbody><tr><td rowspan="8">LCC</td><td>Tree</td><td>Yes</td><td>GPR</td><td>0.68</td><td>5.95</td><td>16%</td><td>0.55</td><td>5.93</td></tr><tr><td>Tree</td><td>Yes</td><td>RFR</td><td>0.71</td><td>5.74</td><td>16%</td><td>0.74</td><td>5.69</td></tr><tr><td>Tree</td><td>No</td><td>GPR</td><td>0.67</td><td>6.08</td><td>17%</td><td>0.55</td><td>6.06</td></tr><tr><td>Tree</td><td>No</td><td>RFR</td><td>0.69</td><td>5.87</td><td>16%</td><td>0.69</td><td>5.83</td></tr><tr><td>Pixel</td><td>Yes</td><td>GPR</td><td>0.79</td><td>4.38</td><td>12%</td><td>0.25</td><td>4.38</td></tr><tr><td>Pixel</td><td>Yes</td><td>RFR</td><td>0.82</td><td>4.16</td><td>11%</td><td>0.46</td><td>4.14</td></tr><tr><td>Pixel</td><td>No</td><td>GPR</td><td>0.77</td><td>4.57</td><td>13%</td><td>0.28</td><td>4.56</td></tr><tr><td>Pixel</td><td>No</td><td>RFR</td><td>0.80</td><td>4.34</td><td>12%</td><td>0.43</td><td>4.32</td></tr><tr><td rowspan="8">LAI</td><td>Tree</td><td>Yes</td><td>GPR</td><td>0.29</td><td>1.80</td><td>37%</td><td>0.36</td><td>1.76</td></tr><tr><td>Tree</td><td>Yes</td><td>RFR</td><td>0.35</td><td>1.77</td><td>36%</td><td>0.49</td><td>1.7</td></tr><tr><td>Tree</td><td>No</td><td>GPR</td><td>0.13</td><td>2.01</td><td>43%</td><td>0.47</td><td>1.96</td></tr><tr><td>Tree</td><td>No</td><td>RFR</td><td>0.17</td><td>1.98</td><td>41%</td><td>0.56</td><td>1.90</td></tr><tr><td>Pixel</td><td>Yes</td><td>GPR</td><td>0.51</td><td>1.2</td><td>31%</td><td>0.21</td><td>1.18</td></tr><tr><td>Pixel</td><td>Yes</td><td>RFR</td><td>0.56</td><td>1.15</td><td>30%</td><td>0.27</td><td>1.12</td></tr><tr><td>Pixel</td><td>No</td><td>GPR</td><td>0.37</td><td>1.36</td><td>35%</td><td>0.27</td><td>1.33</td></tr><tr><td>Pixel</td><td>No</td><td>RFR</td><td>0.39</td><td>1.36</td><td>34%</td><td>0.30</td><td>1.32</td></tr><tr><td rowspan="8">CCC</td><td>Tree</td><td>Yes</td><td>GPR</td><td>0.48</td><td>58.58</td><td>37%</td><td>12.14</td><td>57.31</td></tr><tr><td>Tree</td><td>Yes</td><td>RFR</td><td>0.52</td><td>58.73</td><td>37%</td><td>16.16</td><td>56.47</td></tr><tr><td>Tree</td><td>No</td><td>GPR</td><td>0.40</td><td>62.93</td><td>42%</td><td>14.19</td><td>61.31</td></tr><tr><td>Tree</td><td>No</td><td>RFR</td><td>0.44</td><td>62.13</td><td>41%</td><td>16.47</td><td>59.91</td></tr><tr><td>Pixel</td><td>Yes</td><td>GPR</td><td>0.64</td><td>36.88</td><td>30%</td><td>6.24</td><td>36.35</td></tr><tr><td>Pixel</td><td>Yes</td><td>RFR</td><td>0.68</td><td>36.39</td><td>29%</td><td>8.59</td><td>35.37</td></tr><tr><td>Pixel</td><td>No</td><td>GPR</td><td>0.56</td><td>40.94</td><td>34%</td><td>7.78</td><td>40.20</td></tr><tr><td>Pixel</td><td>No</td><td>RFR</td><td>0.58</td><td>40.84</td><td>33%</td><td>9.09</td><td>39.82</td></tr></tbody></table>

For LCC estimation, R² ranged from 0.68–0.82. SA influenced LCC accuracy the most, with R² of 0.67–0.71 for tree-scale SA and 0.77–0.82 for pixel-scale SA. The two MLRAs performed similarly, with slightly higher R² and lower RMSE for RFR (R² = 0.80–0.82 and RMSE = 4.16–4.34 µg/cm²) than for GPR (R² = 0.77–0.79 and RMSE = 4.38–4.57 µg/cm²). Nevertheless, GPR had a smaller bias than RFR. Inclusion of EFs did not influence LCC estimation significantly, with models that did so having only slightly better evaluation
metrics. Overall, the best model for estimating LCC was RFR with pixel-scale SA and inclusion of EFs, which had high accuracy during cross-validation.

For LAI estimation, R² ranged from 0.13–0.56. To a greater extent than that for LCC models, SA influenced accuracy the most, with R² of 0.13–0.35 for tree-scale SA and 0.37–0.52 for pixel-scale SA. The type of MLRA had a pattern similar to that for LCC models, with RFR having slightly higher R², lower RMSE, but smaller bias than those for GPR. Unlike for LCC models, inclusion of EFs influenced LAI estimation strongly, leading the RFR model with pixel-scale SA and inclusion of EFs to have R² 0.17 higher than that of the same configuration but without EFs.

CCC estimation followed the same pattern as those for LCC and LAI estimation. The best model for estimating CCC was RFR with pixel-scale SA and inclusion of EFs (R² = 0.68 and RMSE = 36.39 µg/m²).

## 3.2. Accuracy of LCC, LAI, and CCC Estimation Using the Real Dataset

### 3.2.1. Overall Accuracy Assessment and Consistency between Cross-Validation and Validation Using the Real Dataset

For LCC accuracy, R² ranged from 0.24–0.33 (Table 11). SA influenced LCC accuracy slightly, with R² of 0.24–0.3 for tree-scale SA and 0.24–0.33 for pixel-scale SA. RFR outperformed GPR (SMAPE = 16% and 21%, respectively). Inclusion of EFs did not influence LCC estimation significantly, with models that did so having only a small decrease in bias. Similar to the results obtained using cross-validation, the LCC model with the best performance was RFR with pixel-scale SA and inclusion of EFs.

<table><caption>**Table 11.** Evaluation metrics used to assess accuracy of leaf chlorophyll content (LCC), leaf area index (LAI), and canopy chlorophyll content (CCC) using the real dataset. SA: spatial allocation; MLRA: machine learning regression algorithm, RMSE: root mean square error; SMAPE: symmetric mean absolute percentage error, BSD: bias standard deviation; RMSE is expressed in µg/cm² for LCC and CCC, m²/m² for LAI.</caption><thead><tr><th>Target Variable</th><th>SA Method</th><th>Environment Features</th><th>MLRA</th><th>R<sup>2</sup></th><th>RMSE</th><th>SMAPE</th><th>BIAS</th><th>BSD</th></tr></thead><tbody><tr><td rowspan="6">LCC</td><td>Tree</td><td>Yes</td><td>GPR</td><td>0.24</td><td>7.66</td><td>22%</td><td>-1.34</td><td>7.55</td></tr><tr><td>Tree</td><td>Yes</td><td>RFR</td><td>0.30</td><td>5.83</td><td>16%</td><td>-1.57</td><td>5.61</td></tr><tr><td>Tree</td><td>No</td><td>GPR</td><td>0.24</td><td>7.48</td><td>22%</td><td>-1.26</td><td>7.38</td></tr><tr><td>Tree</td><td>No</td><td>RFR</td><td>0.29</td><td>6.01</td><td>17%</td><td>-1.87</td><td>5.72</td></tr><tr><td>Pixel</td><td>Yes</td><td>GPR</td><td>0.27</td><td>7.46</td><td>21%</td><td>-1.78</td><td>7.25</td></tr><tr><td>Pixel</td><td>Yes</td><td>RFR</td><td>0.33</td><td>5.64</td><td>16%</td><td>-1.87</td><td>5.32</td></tr><tr><td rowspan="6">LAI</td><td>Pixel</td><td>No</td><td>GPR</td><td>0.24</td><td>7.45</td><td>21%</td><td>-1.71</td><td>7.25</td></tr><tr><td>Pixel</td><td>No</td><td>RFR</td><td>0.32</td><td>5.87</td><td>17%</td><td>-2.29</td><td>5.40</td></tr><tr><td>Tree</td><td>Yes</td><td>GPR</td><td>0.02</td><td>1.94</td><td>43%</td><td>-0.28</td><td>1.92</td></tr><tr><td>Tree</td><td>Yes</td><td>RFR</td><td>0.03</td><td>1.59</td><td>36%</td><td>-0.43</td><td>1.54</td></tr><tr><td>Tree</td><td>No</td><td>GPR</td><td>0.04</td><td>1.54</td><td>35%</td><td>-0.29</td><td>1.51</td></tr><tr><td>Tree</td><td>No</td><td>RFR</td><td>0.04</td><td>1.43</td><td>33%</td><td>-0.28</td><td>1.40</td></tr><tr><td rowspan="6">CCC</td><td>Pixel</td><td>Yes</td><td>GPR</td><td>0.12</td><td>1.63</td><td>57%</td><td>-0.87</td><td>1.38</td></tr><tr><td>Pixel</td><td>Yes</td><td>RFR</td><td>0.29</td><td>1.18</td><td>47%</td><td>-0.58</td><td>1.03</td></tr><tr><td>Pixel</td><td>No</td><td>GPR</td><td>0.25</td><td>1.47</td><td>53%</td><td>-0.90</td><td>1.17</td></tr><tr><td>Pixel</td><td>No</td><td>RFR</td><td>0.22</td><td>1.30</td><td>49%</td><td>-0.69</td><td>1.1</td></tr><tr><td>Tree</td><td>Yes</td><td>GPR</td><td>0.12</td><td>57.11</td><td>43%</td><td>-10.94</td><td>56.08</td></tr><tr><td>Tree</td><td>Yes</td><td>RFR</td><td>0.13</td><td>49.22</td><td>39%</td><td>-12.85</td><td>47.54</td></tr><tr><td rowspan="6">CCC</td><td>Tree</td><td>No</td><td>GPR</td><td>0.16</td><td>49.77</td><td>38%</td><td>-10.86</td><td>48.59</td></tr><tr><td>Tree</td><td>No</td><td>RFR</td><td>0.13</td><td>48.84</td><td>38%</td><td>-13.23</td><td>47.03</td></tr><tr><td>Pixel</td><td>Yes</td><td>GPR</td><td>0.27</td><td>51.00</td><td>58%</td><td>-29.01</td><td>41.97</td></tr><tr><td>Pixel</td><td>Yes</td><td>RFR</td><td>0.46</td><td>36.44</td><td>49%</td><td>-18.79</td><td>31.23</td></tr><tr><td>Pixel</td><td>No</td><td>GPR</td><td>0.31</td><td>48.99</td><td>56%</td><td>-29.73</td><td>38.96</td></tr><tr><td>Pixel</td><td>No</td><td>RFR</td><td>0.32</td><td>42.94</td><td>54%</td><td>-24.26</td><td>35.45</td></tr></tbody></table>
For LAI accuracy, SA had a significant influence: the four models with tree-scale SA had R² close to 0, unlike those with pixel-scale SA, for which R² ranged from 0.12–0.29. For the MLRAs, GPR did not achieve SMAPE < 50%, unlike RFR, which had SMAPE of 47% and 49%. Unlike for the LCC models, inclusion of EFs improved LAI accuracy. Comparing the two RFR models with the pixel-scale SA, the one that included EFs had higher accuracy than the one that excluded them, with a higher R² (0.29 vs. 0.22, respectively), lower RMSE (1.18 vs. 1.30 m²/m², respectively) and lower bias (–0.58 vs. –0.69, respectively).

For CCC accuracy, SA also had a significant influence: models with pixel-scale SA had much higher R² (0.27–0.46) than those with tree-scale SA (0.12–0.16). Similarly, RMSE was lower with pixel-scale SA (36–51 µg/m²) than with tree-scale SA (49–57 µg/m²), but the bias was much larger with pixel-scale SA. The RFR model with inclusion of EFs outperformed the other models according to all metrics.

### 3.2.2. Accuracy of LCC, LAI, and CCC Estimation by Tree Species

Using the best model for each target variable, we estimated the target values by tree species. The model used to estimate LCC showed high tree-species dependence, with R² ranging from 0.11 for PL to 0.45 for QR (Table 12).

<table><caption>**Table 12.** Evaluation metrics for assessing model accuracy by tree species using the real dataset. Only the best model configurations are shown (i.e., MLRA = RFR, environmental features included = yes; spatial allocation = pixel scale). SA: spatial allocation; MLRA: machine learning regression algorithm, RMSE: root mean square error; SMAPE: symmetric mean absolute percentage error, BSD: bias standard deviation; RMSE is expressed in µg/cm² for LCC and CCC, m²/m² for LAI.</caption><thead><tr><th>Target Variable</th><th>Species</th><th>R<sup>2</sup></th><th>RMSE</th><th>SMAPE</th><th>BIAS</th><th>BSD</th></tr></thead><tbody><tr><td rowspan="4">LCC</td><td>AC</td><td>0.26</td><td>5.48</td><td>0.16</td><td>-2.45</td><td>4.91</td></tr><tr><td>FR</td><td>0.41</td><td>5.45</td><td>0.14</td><td>-2.07</td><td>5.06</td></tr><tr><td>PL</td><td>0.11</td><td>5.66</td><td>0.16</td><td>-0.36</td><td>5.66</td></tr><tr><td>QR</td><td>0.45</td><td>5.99</td><td>0.17</td><td>-2.81</td><td>5.29</td></tr><tr><td rowspan="4">LAI</td><td>AC</td><td>0.17</td><td>0.95</td><td>0.36</td><td>-0.42</td><td>0.85</td></tr><tr><td>FR</td><td>0.36</td><td>1.26</td><td>0.47</td><td>-0.61</td><td>1.10</td></tr><tr><td>PL</td><td>0.03</td><td>1.34</td><td>0.65</td><td>-0.98</td><td>0.92</td></tr><tr><td>QR</td><td>0.43</td><td>1.15</td><td>0.37</td><td>-0.18</td><td>1.14</td></tr><tr><td rowspan="4">CCC</td><td>AC</td><td>0.28</td><td>31.07</td><td>0.42</td><td>-17.75</td><td>25.53</td></tr><tr><td>FR</td><td>0.5</td><td>41.05</td><td>0.48</td><td>-21.54</td><td>35.03</td></tr><tr><td>PL</td><td>0.1</td><td>37.74</td><td>0.64</td><td>-26.42</td><td>26.99</td></tr><tr><td>QR</td><td>0.55</td><td>36.26</td><td>0.38</td><td>-6.96</td><td>35.66</td></tr></tbody></table>

FR and QR followed a similar pattern for LCC (Figure 6), with overestimation from low to median LCC (10–30 µg/cm²). Conversely, LCC patterns for AC and PL were more scattered. We also found a phenological trend, with LCC lower for early dates (i.e., day of the year (DOY) 126 and 166) than for later dates.
<figure><img src="image_10.png" /><figcaption>**Figure 6.** Scatterplot of estimated vs. observed leaf chlorophyll content (LCC) (spatial allocation = pixel scale, machine-learning regression algorithm = RFR, environmental features included = yes) for the four tree species: *Acer platanoides* (AC), *Fraxinus excelsior* (FR), *Platanus acerifolia* (PL), *Quercus rubra* (QR). DOY: day of year of the Sentinel-2 image.</figcaption></figure>

The model used to estimate LAI was also species-dependent, with R² ranging from near-zero for PL (0.03) to moderate for QR (0.43). A similar pattern was found for the four species, with overestimation of low LAI (0–2), and few LAI < 2 being estimated. The best performances were observed for QR and FR, but with a positive bias (i.e., overestimation) for LAI of 0–2.5 and a negative bias for LAI of 2.5–6.0 (Figure 7).
<figure><img src="image_11.png" /><figcaption>**Figure 7.** Scatterplot of estimated vs. observed leaf area index (LAI) (spatial allocation = pixel scale, machine-learning regression algorithm = RFR, environmental features included = yes) for the four tree species: *Acer platanoides* (AC), *Fraxinus excelsior* (FR), *Platanus acerifolia* (PL), and *Quercus rubra* (QR). DOY: day of year of the Sentinel-2 image.</figcaption></figure>

CCC estimation followed a pattern similar to that of LAI, with low R² for PL and AC (0.10 and 0.28, respectively) and clustering of estimated values (Figure 8). Performances were better for QR and FR (R² of 0.55 and 0.50, respectively), which followed a distinctive pattern with overestimated low CCC and underestimated high CCC, and an inflection point around 100 µg/m².
<figure><img src="image_12.png" /><figcaption>**Figure 8.** Scatterplot of estimated vs. observed canopy chlorophyll content (CCC) (spatial al- location = pixel scale, machine-learning regression algorithm = RFR, environmental features in- cluded = yes) for the four tree species: *Acer platanoides* (AC), *Fraxinus excelsior* (FR), *Platanus acerifolia* (PL), and *Quercus rubra* (QR). DOY: day of the year (DOY) of the Sentinel-2 image.</figcaption></figure>

Overall, for the three target variables, PL and AC followed a similar pattern, while FR
and QR followed another similar pattern.

3.3. *Consistency between Estimated and Observed Time Series*

For each target variable and species, we examined the mean time series (Figure 9)
and distribution of temporal metric values (Figure 10) to assess the temporal consistency
between estimated and observed series. High adequacy was observed between estimated and observed LCC time series according to d<sub>EUCL</sub> and d<sub>CORT</sub> (Figure 9). There were no significant differences among the four species. In particular, the observed LCC time series decreased from DOY 166 to 201 and then increased from DOY 201 to 226. This pattern was also found, to a lesser extent, in the estimated LCC time series for AC and QR. For LAI, the estimated time series differed much more from the observed time series. Temporal trends were observed for FR (d<sub>CORT</sub> = 0.04) and QR (d<sub>CORT</sub> = 0.03). AC and PL showed differences between the estimated and observed time series. Values of d<sub>EUCL</sub> were particularly high for PL. The estimated and observed CCC time series showed high adequacy, with similar temporal trends for AC and FR. However, a strong positive bias in the estimated CCC time series was observed for AC, FR, and PL (d<sub>EUCL</sub> = 7.41, 9.52, and 10.08, respectively).

<figure><img src="image_13.png" /><figcaption>**Figure 9.** Observed (blue) and estimated (red) mean time series of LCCpix, LAIpix, and CCCpix by tree species (AC, FR, PL, and QR). LCCpix: leaf chlorophyll content at pixel scale; LAIpix: leaf area index at pixel scale; CCCpix: canopy chlorophyll content at pixel scale.</figcaption></figure>
<figure><img src="image_14.png" /><figcaption>**Figure 10.** Distribution of similarity metrics for each target variable by tree species (QR, PL, FR, and AC) and metric ($d_{\text{EUCL}}$, $d_{\text{COR}}$, and $d_{\text{CORT}}$). The curves represent the relative shape of the distributions. QR: *Quercus rubra*; PL: *Platanus acerifolia*; FR: *Fraxinus excelsior*; and AC: *Acer platanoides*. $d_{\text{EUCL}}$: Euclidean distance; $d_{\text{COR}}$: Pearson correlation distance; $d_{\text{CORT}}$: temporal correlation distance.</figcaption></figure>

## 4. Discussion

Results highlighted two main points: (i) evaluating accuracy using cross-validation revealed that (i) pixel-scale SA estimated all target variables better than tree-scale SA, (ii) including EFs in model training improved model performance (especially for LAI and CCC), and (iii) RFR gave slightly better estimates than GPR. Comparing the three target variables, the models had good accuracy for LCC estimation and moderate accuracy for CCC and LAI estimation (SMAPE = 11%, 29%, and 30%, respectively); (ii) the highest $R^2$ obtained using the RDS was 0.46 for CCC, 0.33 for LCC, and 0.29 for LAI. We reached the same conclusions using the RDS as when using SDS: for all three target variables, the best configuration for the model was RFR with pixel-scale SA and inclusion of EFs. However, the models had lower overall accuracy when using the RDS than when using SDS.

### 4.1. Overall Performance of Vegetation Trait Estimation

The best performances for LCC estimation in the study were $R^2$ of 0.33 with high variability among species, $R^2$ ranging from 0.11 (PL) to 0.45 (QR), and RMSE of 5.64 µg/cm². Although LCC, LAI, and CCC of urban trees had never been estimated using S2 images,
LCC has been estimated using a non-species-specific model for urban trees based on hyperspectral images at 1 m spatial resolution in an urban area that contained PL, with R² of 0.77 [88]. Another study that used species-specific models in an urban area based on hyperspectral images at 2 m spatial resolution [89] had R² of 0.33 and 0.62, for AC and PL, respectively.

The combination of S2 data with hybrid inversion models based on PROSAIL has been successfully used to estimate LAI and LCC with high accuracy. A study conducted on tree plantations showed R² values above 0.8 for LAI estimation [36], and another study including seven vegetation types (crops, forests, wetlands, ...) showed an R² value of 0.62 for LCC estimation [37]. However, these models did not take into account the complexity of urban environments.

Among studies that estimated vegetation traits using S2 imagery, the most similar case studies were based on forest canopies or sparsely wooded areas and used hybrid inversion methods based on DART or INFORM. Comparing the best model in the present study (RFR with pixel-scale SA and inclusion of EFs) to INFORM inversion for mixed forest [27], we found similar validation R² for LCC estimation (0.33 vs. 0.34, respectively). For LAI and CCC estimation, however, the present study's models had lower R² than those of the other study: 0.29 vs. 0.47 and 0.46 vs. 0.65, respectively. However, in both studies, model validation using simulated datasets showed the same pattern, with very good performance for LCC, good performance for CCC, and moderate performance for LAI. Another study based on DART and based on an agricultural area with olive orchards had good performances for LAI (RMSE = 0.58 m²/m²) and LCC (RMSE = 2.5 µg/cm²) [42]. These good performances can be explained by the high detail in 3D models of the trees and prior knowledge of soil optical properties extracted from S2 pure pixels, which was not possible in the present study because it is nearly impossible to extract pure pixels of soils in urban areas from images with a spatial resolution of 10 m, and because soil composition varies greatly among pixels in urban areas.

Validation using the real dataset did not show high accuracy for estimation of vegetation traits. However, while high estimation accuracy may be needed for certain eco-physiological applications in agricultural areas, moderate accuracy may be sufficient for monitoring purposes in urban areas. There was high temporal consistency between observed and estimated vegetation traits, particularly according to the dCORT metric, which highlighted similar trends for estimated and observed values. The temporal trend in vegetation traits may be more informative than the absolute values of these traits when studying the phenology and response of vegetation to stress events [90].

## 4.2. DART Modelling

The DART model has been used in studies that examined various aspects of the urban environment, such as retrieving spectral signatures [91,92], assessing radiative budgets [41], investigating temperature patterns [93], and studying impacts of 3D structures on emissivity [94]. Nevertheless, to the best of our knowledge, no inversion models of vegetation traits of urban trees have been developed using DART or other RTMs. While physical approaches are inherently general and applicable in many contexts, they are limited by competing demands of model realism and inversion feasibility. A simplified approach to complex physical processes is usually required to constrain inversion in a remote sensing context [95,96]. Unlike some approaches used, mainly for calibration, that involve direct simulations with high accuracy and definition of 3D scenes, this study's approach is more global and aims to cover a wide range of urban configurations and model pseudo-realistic urban scenes. Thus, we made several simplifications, notably to the 3D tree model. Tree structure determined the proportion of woody elements in the crown, as branch shape, length, and density can strongly influence reflectance, especially when LAI is low [50], and thus the accuracy of LAI estimation. The present study's two simplifying assumptions (crown height equalled 2/3 of tree height and the crown was an ellipsoid) enabled us first to acquire a validation dataset using LAI-2200 measurements and a DSM—
without requiring terrestrial LIDAR data, which are much more difficult to acquire—and
then to aggregate these values at the pixel scale. However, the live crown ratio and crown
shape may be species-specific [79], which may explain the overall moderate performance
of LAI estimation and the variability in performances among tree species.

4.3. *Environment Features*

One key point of this study was the use of contextual information for simulation
and model training. The LCZ classification system was used as a framework for the
construction of pseudo-realistic urban scenes. While the four LCZ types were selected
among the most dominant LCZ built-cover types encountered in European cities with
more than 100,000 inhabitants, these LCZ types did not cover all the urban environments
of our study site. In particular, the dense, high-rise urban centre corresponding to LCZ
type 1 (“Compact high-rise”) was excluded. However, it does not affect the transferability
of the models, because tree vegetation is less prevalent in this LCZ type than in the four
LCZ types selected. Nevertheless, many tree-exogenous parameters were used for the
DART simulations, including differences in LCZ type, material type, contextual parameter,
and underlying vegetation. However, using a large number of parameters can increase
the probability of encountering ill-posed problems, when nearly identical spectra may
correspond to different combinations of model input parameters [97,98]. To mitigate
these issues and improve model performance, many studies developed regularization
strategies, which incorporate supplementary information and adopt constrained inversion
methods [99–102]. However, these strategies often rely on adjusting RTM input parameters
to better match the magnitude of in situ vegetation traits or contextual features such as
soil spectra. While these approaches can improve performance, they may decrease the
generality of the resulting inversion models.

In the present study, contextual input information was summarized into three EFs
at the pixel scale (CCP, Pshadow, and Pgreen). The strategy was to consider using these
features to train the models (as they interact with spectral features) rather than to filter or
restrict the simulated data so that their characteristics match those of the RDS. However,
building the model in this way required obtaining real data that corresponded to these
features to apply the inversion model. Vegetation cover and DSM data are now widely
available at the scale of major European cities via Open Data city portals, OpenStreetMap,
the Copernicus portal, and Google Earth Engine, including a highly accurate, Europe-wide
DSM [103].

The models that included EFs performed better for all three vegetation traits, especially LAI. In the machine-learning approach, EFs were considered interaction features, as they were not directly correlated with vegetation traits, although they strongly influenced spectral features. One advantage of RFR is its ability to exploit interaction effects between the features effectively [104]. Including these three EFs allowed a wider range of pixels of interest to be considered (first for the validation dataset, and then for mapping), without having to filter pure pixels, pixels in shadow, or pixels with underlying vegetation, which makes this approach particularly efficient for urban areas.

4.4. *Spatial Allocation*

The models’ performances showed the importance of the SA method used. Corre-
spondence between the spatial coverage of the objects of interest (i.e., trees) and satellite
images is essential for model validation, particularly in an urban environment in which the
dimensions of the objects of interest are similar to those of the pixel. Correspondence is not
a problem when the canopy studied is spatially homogeneous and in situ measurements
can be performed in a homogeneous way. In many studies, in situ measurements are
made at the plot scale, using a square that is usually larger than the resolution of the
satellite image. For example, using S2 imagery, Ali et al. [27] and Brown et al. [26] took
measurements in plots with 30 and 40 m sides, respectively. In sparse canopy, methods
need to be adapted to match satellite and in situ data. For example, Makhloufi et al. [42]
measured the LAI of isolated trees and then related the measurements to the area of a plot to estimate total LAI at the plot scale. We followed a similar approach, which considered the presence of several trees and their canopy cover in the pixel. This approach was able to address details of urban vegetation layout, such as street tree alignments. Although overall performance was moderate, the fact that pixel-scale models outperformed tree-scale models demonstrated the benefits of this approach.

**5. Conclusions**

This study examined an innovative approach for estimating vegetation traits LCC, LAI, and CCC of urban trees using S2 imagery and a hybrid inversion method based on the DART model. Combining DART with RFR and GPR effectively captured the complexity and heterogeneity of urban environments. One major contribution of this study is the inclusion of EFs such as CCP, Pshadow, and Pgreen in the regression models, which significantly improved model performance, particularly for estimating LAI and CCC. The study also showed that pixel-scale SA provided more accurate results for all vegetation traits, highlighting the need for precise spatial matching between in situ measurements and satellite data in urban areas. Species-specific analysis revealed variability in model performance, with certain species showing higher accuracy. Despite moderate overall performance, the models showed strong temporal consistency, particularly for LCC, making them valuable for monitoring intra-annual urban tree dynamics. This research highlights the potential of freely available Sentinel-2 imagery for monitoring and managing urban trees. Future research will focus on species-specific modelling, extending the method to different urban contexts, and exploring the inclusion of additional environment variables to further improve model performance and applicability.

**Author Contributions:** T.L.S.: Conceptualization, visualization, methodology, investigation, writing—original draft; J.N.: Conceptualization, methodology, investigation, project administration, funding acquisition, supervision; L.H.-M.: Supervision, writing—review and editing; K.A.: Conceptualization, methodology, investigation, supervision, writing—review and editing. All authors have read and agreed to the published version of the manuscript.

**Funding:** This study was supported by Rennes Métropole and the Association Nationale de la Recherche et de la Technologie (ASTRESS project and grant no. 2021/0301).

**Data Availability Statement:** Field survey data for LCC and LAI measurement are available at https://zenodo.org/records/12751353 (accessed on 17 July 2024).

**Conflicts of Interest:** The authors declare no conflicts of interest.

**Appendix A. Detailed Description of DART Input Parameters**

Appendix A.1. *Tree-Exogenous Parameters*

Appendix A.1.1. Sensor Settings, Direction Input Parameter, and Atmosphere

The illumination of the scene is important, as it determines the distance of the shadow cast by the buildings (and the energy balance of the scene). These conditions change throughout the year, characterized by a change in the solar incidence angles (zenith and azimuth). To reproduce this change, the date can be set in DART (which will determine the solar angles for a given latitude and longitude) and then changed. We considered dates from 15 March to 15 November, which coincides with the vegetative period of the trees. Concerning parameters for the atmosphere, urban areas have a specific atmosphere because the concentration of housing and industry produces large amounts of emissions that influence the composition of the atmosphere and aerosols. Images acquired in urban areas are therefore influenced more by absorption and scattering, mainly in the visible and near-infrared wavelengths. The DART model can use the USSTD76 atmospheric model [105] to reproduce the atmospheric influence on the signal. Optical properties of aerosols can also be set; this study used an urban aerosol with a constant multiplier of one for the optical depth of aerosols.
## Appendix A.1.2. Spectral Library

One characteristic of the urban environment is the variety of materials, most of which are human-made or mineral, such as rough stone (e.g., sandstone, limestone), granite, asphalt, cement, concrete, brick, shingle, slate, red tile, and glass, or even metals, such as steel and zinc. Each material has specific spectral signatures due to its intrinsic physical and chemical composition that will influence the land-surface radiation balance in the scene, as well as the reflectance of the pixel analyzed. Due to the geometric configuration of the four LCZs selected and the simplified 3D representations of the buildings, and to reproduce plausible contexts that represented a broad range of materials, four classes were considered: roofs, walls, impermeable ground, and permeable ground. Each of these classes was assigned potential reflectance spectra from a DART spectral library with different levels of detail. A total of 18 spectra were retained (Figure A1).

<figure><img src="image_15.png" /><figcaption>**Figure A1.** Spectral library showing the four spectral categories for (a) pervious ground (n = 2), (b) impervious ground (n = 4), (c) walls (n = 6) and (d) roofs (n = 6). Each spectrum is shown with its reflectance value [%] as a function of wavelength [µm]. The vertical bands correspond to the central wavelength and width of Sentinel-2 spectral bands.</figcaption></figure>
### Appendix A.1.3. Earth Scene and Tree Planting Context

The scene was set in Rennes, France, to match the real Sentinel-2 data. The scene dimensions were set to 100 × 100 m based on the scale of LCZs, and to ensure that electromagnetic interactions in the tree's environment could actually occur. The conditions under which a tree is planted are a key factor in its development. It is important to distinguish two elements: the absolute planting conditions, which include the type of ground (e.g., soil, grass, human-made materials) at the base of the tree, and the environmental conditions (e.g., the presence of buildings, the street). The former can be classified as restricted (a ca. 1 m² square of bare soil or grass surrounded by impermeable materials), linear (trees planted on a strip of bare soil or grass surrounded by impermeable materials), or open (trees located on extended grass or in a park) [106]. In the simulations, these conditions were reflected by the ratio of permeable to impermeable ground defined with the LCZs (Table 1). The latter can be described using three parameters: the distance from the tree to the nearest building, street orientation, and exposure (planted on the sunny or shady side of the street). These three parameters are likely to strongly influence how shade will impact the tree.

### Appendix A.2. Tree-Endogenous Parameters

#### Appendix A.2.1. Tree Structural Parameters

Tree structural parameters were determined by analysing the tree row database of the Lyon metropolitan area, which is one of the most complete databases for France in Open-Data (https://data.grandlyon.com/jeux-de-donnees/arbres-alignement-metropole-lyon/info; accessed on 17 July 2024). The database contains 50,000 records with the following geometric variables: crown height, crown diameter, total height, and trunk diameter. To build a reference profile, the database was limited to the five most common genera—*Tilia*, *Acer*, *Platanus*, *Fraxinus*, and *Quercus*—which agrees with urban tree databases for other European cities, and to crown diameters of 8–14 m. With these limits, mean values were calculated for tree structural parameters (Table A1). The tree was further modelled with an ellipsoid crown and a cylindrical trunk without branches.

<table><caption>**Table A1.** Structural parameters of the tree modelled in DART.</caption><thead><tr><td>Parameter</td><td>Value [m]</td></tr></thead><tbody><tr><td>Tree height</td><td>15</td></tr><tr><td>Trunk height under the crown</td><td>4</td></tr><tr><td>Trunk height in the crown</td><td>6</td></tr><tr><td>Trunk diameter</td><td>0.4</td></tr></tbody></table>

Three parameters were variable: leaf area density (LAD), the percentage of crown holes, and leaf-angle distribution. The leaf area index (LAI) was calculated as the total leaf area of a canopy per unit area of ground, which is suitable for a tree canopy with a homogeneous height such as that of dense forests. For a single tree, LAI can be defined using the LAD, which equals the total leaf area per unit volume of canopy (i.e., the tree crown) in m²/m³. Next, a clumping factor indicates the degree of spatial aggregation of leaves in the canopy. DART can calculate the percentage of crown volume occupied by leaves. If 100% of the crown volume is occupied, the LAD is homogeneous at the crown scale and the clumping factor equals 0, but if only 50% of the crown volume is occupied, the LAD is not homogeneous at the crown scale, and leaves are artificially aggregated for a given initial LAD. We used a crown hole percentage from 0–50% as a proxy to simulate the clumping factor. Finally, the leaf angle distribution refers to the statistical distribution of the angular orientation of leaves at the scale of a tree (e.g., planophilic, plagiophilic, extremophilic, uniform). Several methods exist to estimate it, including several predefined mathematical functions. From the database available in [107], which provides leaf-inclination angles for
temperate and boreal broadleaf woody species, the two most common distributions (i.e., planophile and plagiophile) for deciduous tree species were extracted.

## Appendix A.2.2. Leaf Parameters

The DART model is coupled with the PROSPECT model. In this study, the PROSPECT-D version was used [31], as it can simulate the reflectance and transmittance spectra of a leaf in the spectral range of 0.4–2.5 µm (Figure A2). To simulate these spectra, the model uses seven parameters: the structure coefficient (N), chlorophyll content (Cab), carotenoid content (Car), anthocyanin content, brown pigments, equivalent water thickness (EWT), and dry matter (LMA). The simulated spectrum is then assigned to the optical properties of the tree leaves in DART.

<figure><img src="image_16.png" /></figure>

**Figure A2.** Example of two leaf spectra simulated by the PROSPECT model. Spectrum 1: N = 2.3, Cab = 60, Car = 25, EWT = 0.024 and LMA = 0.018. Spectrum 2: N = 1.7, Cab = 32.5, Car = 13.75, EWT = 0.014, and LMA = 0.008. Spectra 1 and 2 correspond to maximum and median values used in the virtual experimental design, respectively. The vertical bands correspond to the central wavelength and width of Sentinel-2 spectral bands. N = structure coefficient; Cab = chlorophyll (a+b) content; Car = carotenoids; EWT = equivalent water thickness; LMA = leaf mass area.
Appendix B. Cross-Validation Scatterplots

<figure><img src="image_17.png" /><figcaption>**Figure A3.** Cross-validation scatterplot for LCC according to the SA method (tree (a) or pixel (b)), MLRA (GPR or RFR) and inclusion of EF (with or without). LCC: leaf chlorophyll content; SA: spatial allocation; MLRA: machine learning regression algorithm; GPR: gaussian process regression; RFR: random forest regression; EF: environment feature.</figcaption></figure>
<figure><img src="image_18.png" /><figcaption>**Figure A4.** Cross-validation scatterplot for LAI according to the SA method (tree (a) or pixel (b)), MLRA (GPR or RFR) and inclusion of EF (with or without). LAI: leaf area index; SA: spatial allocation; MLRA: machine learning regression algorithm; GPR: gaussian process regression; RFR: random forest regression; EF: environment feature.</figcaption></figure>
<figure><img src="image_19.png" /><figcaption>**Figure A5.** Cross-validation scatterplot for CCC according to SA method (tree (a) or pixel (b)), MLRA (GPR or RFR) and inclusion of EF (with or without). CCC: canopy chlorophyll content; SA: spatial allocation; MLRA: machine learning regression algorithm; GPR: gaussian process regression; RFR: random forest regression; EF: environment feature.</figcaption></figure>
## References

1. Xu, F.; Yan, J.; Heremans, S.; Somers, B. Pan-European Urban Green Space Dynamics: A View from Space between 1990 and 2015. *Landsc. Urban Plan.* **2022**, *226*, 104477. [CrossRef]

2. Bolund, P.; Hunhammar, S. Ecosystem Services in Urban Areas. *Ecol. Econ.* **1999**, *29*, 293–301. [CrossRef]

3. Andersson-Sköld, Y.; Thorsson, S.; Rayner, D.; Lindberg, F.; Janhäll, S.; Jonsson, A.; Moback, U.; Bergman, R.; Granberg, M. An Integrated Method for Assessing Climate-Related Risks and Adaptation Alternatives in Urban Areas. *Clim. Risk Manag.* **2015**, *7*, 31–50. [CrossRef]

4. Nowak, D.J.; Crane, D.E. Carbon Storage and Sequestration by Urban Trees in the USA. *Environ. Pollut.* **2002**, *116*, 381–389. [CrossRef]

5. Andersson, E.; Barthel, S.; Ahrné, K. Measuring Social-Ecological Dynamics Behind the Generation of Ecosystem Services. *Ecol. Appl.* **2007**, *17*, 1267–1278. [CrossRef]

6. Wolf, K.L.; Lam, S.T.; McKeen, J.K.; Richardson, G.R.A.; van den Bosch, M.; Bardekjian, A.C. Urban Trees and Human Health: A Scoping Review. *Int. J. Environ. Res. Public Health* **2020**, *17*, 4371. [CrossRef]

7. Czaja, M.; Kołton, A.; Muras, P. The Complex Issue of Urban Trees—Stress Factor Accumulation and Ecological Service Possibilities. *Forests* **2020**, *11*, 932. [CrossRef]

8. Sæbø, A.; Borzan, Ž.; Ducatillion, C.; Hatzistathis, A.; Lagerström, T.; Supuka, J.; García-Valdecantos, J.L.; Rego, F.; Van Slycken, J. The Selection of Plant Materials for Street Trees, Park Trees and Urban Woodland. In *Urban Forests and Trees: A Reference Book*; Konijnendijk, C., Nilsson, K., Randrup, T., Schipperijn, J., Eds.; Springer: Berlin/Heidelberg, Germany, 2005; pp. 257–280, ISBN 978-3-540-27684-5.

9. Ma, B.; Hauer, R.J.; Östberg, J.; Koeser, A.K.; Wei, H.; Xu, C. A Global Basis of Urban Tree Inventories: What Comes First the Inventory or the Program. *Urban For. Urban Green.* **2021**, *60*, 127087. [CrossRef]

10. Hilbert, D.; Roman, L.; Koeser, A.; Vogt, J.; van Doorn, N. Urban Tree Mortality: A Literature Review. *Arboric. Urban For. AUF* **2019**, *45*, 167–200. [CrossRef]

11. Neyns, R.; Canters, F. Mapping of Urban Vegetation with High-Resolution Remote Sensing: A Review. *Remote Sens.* **2022**, *14*, 1031. [CrossRef]

12. García-Pardo, K.A.; Moreno-Rangel, D.; Domínguez-Amarillo, S.; García-Chávez, J.R. Remote Sensing for the Assessment of Ecosystem Services Provided by Urban Vegetation: A Review of the Methods Applied. *Urban For. Urban Green.* **2022**, *74*, 127636. [CrossRef]

13. Mattila, H.; Valev, D.; Havurinne, V.; Khorobrykh, S.; Virtanen, O.; Antinluoma, M.; Mishra, K.B.; Tyystjärvi, E. Degradation of Chlorophyll and Synthesis of Flavonols during Autumn Senescence—the Story Told by Individual Leaves. *AoB Plants* **2018**, *10*, ply028. [CrossRef] [PubMed]

14. Croft, H.; Chen, J.M.; Luo, X.; Bartlett, P.; Chen, B.; Staebler, R.M. Leaf Chlorophyll Content as a Proxy for Leaf Photosynthetic Capacity. *Glob. Change Biol.* **2017**, *23*, 3513–3524. [CrossRef] [PubMed]

15. Liu, C.; Liu, Y.; Lu, Y.; Liao, Y.; Nie, J.; Yuan, X.; Chen, F. Use of a Leaf Chlorophyll Content Index to Improve the Prediction of Above-Ground Biomass and Productivity. *PeerJ* **2019**, *6*, e6240. [CrossRef]

16. Filella, I.; Penuelas, J. The Red Edge Position and Shape as Indicators of Plant Chlorophyll Content, Biomass and Hydric Status. *Int. J. Remote Sens.* **1994**, *15*, 1459–1470. [CrossRef]

17. Ramos-Montaño, C. Vehicular Emissions Effect on the Physiology and Health Status of Five Tree Species in a Bogotá, Colombia Urban Forest. *Rev. Biol. Trop.* **2020**, *68*, 1001–1015. [CrossRef]

18. Talebzadeh, F.; Valeo, C. Evaluating the Effects of Environmental Stress on Leaf Chlorophyll Content as an Index for Tree Health. *IOP Conf. Ser. Earth Environ. Sci.* **2022**, *1006*, 012007. [CrossRef]

19. Peñuelas, J.; Rutishauser, T.; Filella, I. Phenology Feedbacks on Climate Change. *Science* **2009**, *324*, 887–888. [CrossRef]

20. Duncan, W.G. Leaf Angles, Leaf Area, and Canopy Photosynthesis. *1. Crop Sci.* **1971**, *11*, 482–485. [CrossRef]

21. Halme, E.; Pellikka, P.; Möttus, M. Utility of Hyperspectral Compared to Multispectral Remote Sensing Data in Estimating Forest Biomass and Structure Variables in Finnish Boreal Forest. *Int. J. Appl. Earth Obs. Geoinf.* **2019**, *83*, 101942. [CrossRef]

22. Granero-Belinchon, C.; Adeline, K.; Briottet, X. Impact of the Number of Dates and Their Sampling on a NDVI Time Series Reconstruction Methodology to Monitor Urban Trees with Venus Satellite. *Int. J. Appl. Earth Obs. Geoinf.* **2021**, *95*, 102257. [CrossRef]

23. Verrelst, J.; Rivera, J.P.; Veroustraete, F.; Muñoz-Marí, J.; Clevers, J.G.P.W.; Camps-Valls, G.; Moreno, J. Experimental Sentinel-2 LAI Estimation Using Parametric, Non-Parametric and Physical Retrieval Methods—A Comparison. *ISPRS J. Photogramm. Remote Sens.* **2015**, *108*, 260–272. [CrossRef]

24. Combal, B.; Baret, F.; Weiss, M.; Trubuil, A.; Macé, D.; Pragnère, A.; Myneni, R.; Knyazikhin, Y.; Wang, L. Retrieval of Canopy Biophysical Variables from Bidirectional Reflectance: Using Prior Information to Solve the Ill-Posed Inverse Problem. *Remote Sens. Environ.* **2003**, *84*, 1–15. [CrossRef]

25. Zarco-Tejada, P.J.; Hornero, A.; Beck, P.S.A.; Kattenborn, T.; Kempeneers, P.; Hernández-Clemente, R. Chlorophyll Content Estimation in an Open-Canopy Conifer Forest with Sentinel-2A and Hyperspectral Imagery in the Context of Forest Decline. *Remote Sens. Environ.* **2019**, *223*, 320–335. [CrossRef]
26. Brown, L.A.; Ogutu, B.O.; Dash, J. Estimating Forest Leaf Area Index and Canopy Chlorophyll Content with Sentinel-2: An Evaluation of Two Hybrid Retrieval Algorithms. *Remote Sens. 2019*, *11*, 1752. [CrossRef]

27. Ali, A.M.; Darvishzadeh, R.; Skidmore, A.; Gara, T.W.; Heurich, M. Machine Learning Methods' Performance in Radiative Transfer Model Inversion to Retrieve Plant Traits from Sentinel-2 Data of a Mixed Mountain Forest. *Int. J. Digit. Earth 2021*, *14*, 106–120. [CrossRef]

28. Amin, E.; Verrelst, J.; Rivera-Caicedo, J.P.; Pipia, L.; Ruiz-Verdú, A.; Moreno, J. Prototyping Sentinel-2 Green LAI and Brown LAI Products for Cropland Monitoring. *Remote Sens. Environ.* **2021**, *255*, 112168. [CrossRef]

29. de Sá, N.C.; Baratchi, M.; Hauser, L.T.; van Bodegom, P. Exploring the Impact of Noise on Hybrid Inversion of PROSAIL RTM on Sentinel-2 Data. *Remote Sens.* **2021**, *13*, 648. [CrossRef]

30. Jacquemoud, S.; Baret, F. PROSPECT: A Model of Leaf Optical Properties Spectra. *Remote Sens. Environ.* **1990**, *34*, 75–91. [CrossRef]

31. Féret, J.-B.; Gitelson, A.A.; Noble, S.D.; Jacquemoud, S. PROSPECT-D: Towards Modeling Leaf Optical Properties through a Complete Lifecycle. *Remote Sens. Environ.* **2017**, *193*, 204–215. [CrossRef]

32. Verhoef, W. Light Scattering by Leaf Layers with Application to Canopy Reflectance Modeling: The SAIL Model. *Remote Sens. Environ.* **1984**, *16*, 125–141. [CrossRef]

33. Atzberger, C. Development of an Invertible Forest Reflectance Model: The INFOR-Model. In *A Decade of Trans-European Remote Sensing Cooperation, Proceedings of the 20th EARSeL Symposium, Dresden, Germany, 14–16 June 2000*; University of Trier: Trier, Germany, 2000; Volume 14, pp. 39–44.

34. Jacquemoud, S.; Verhoef, W.; Baret, F.; Bacour, C.; Zarco-Tejada, P.J.; Asner, G.P.; François, C.; Ustin, S.L. PROSPECT + SAIL Models: A Review of Use for Vegetation Characterization. *Remote Sens. Environ.* **2009**, *113*, S56–S66. [CrossRef]

35. Darvishzadeh, R.; Wang, T.; Skidmore, A.; Vrieling, A.; O’Connor, B.; Gara, T.W.; Ens, B.J.; Paganini, M. Analysis of Sentinel-2 and RapidEye for Retrieval of Leaf Area Index in a Saltmarsh Using a Radiative Transfer Model. *Remote Sens.* **2019**, *11*, 671. [CrossRef]

36. Wan, L.; Ryu, Y.; Dechant, B.; Lee, J.; Zhong, Z.; Feng, H. Improving Retrieval of Leaf Chlorophyll Content from Sentinel-2 and Landsat-7/8 Imagery by Correcting for Canopy Structural Effects. *Remote Sens. Environ.* **2024**, *304*, 114048. [CrossRef]

37. Sinha, S.K.; Padalia, H.; Dasgupta, A.; Verrelst, J.; Rivera, J.P. Estimation of Leaf Area Index Using PROSAIL Based LUT Inversion, MLRA-GPR and Empirical Models: Case Study of Tropical Deciduous Forest Plantation, North India. *Int. J. Appl. Earth Obs. Geoinf.* **2020**, *86*, 102027. [CrossRef]

38. Gastellu-Etchegorry, J.-P.; Lauret, N.; Yin, T.; Landier, L.; Kallel, A.; Malenovsky, Z.; Bitar, A.A.; Aval, J.; Benhmida, S.; Qi, J.; et al. DART: Recent Advances in Remote Sensing Data Modeling with Atmosphere, Polarization, and Chlorophyll Fluorescence. *IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens.* **2017**, *10*, 2640–2649. [CrossRef]

39. Wang, Y.; Kallel, A.; Yang, X.; Regaieg, O.; Lauret, N.; Guilleux, J.; Chavanon, E.; Gastellu-Etchegorry, J.-P. DART-Lux: An Unbiased and Rapid Monte Carlo Radiative Transfer Method for Simulating Remote Sensing Images. *Remote Sens. Environ.* **2022**, *274*, 112973. [CrossRef]

40. Zhen, Z.; Benromdhane, N.; Kallel, A.; Wang, Y.; Regaieg, O.; Boitard, P.; Landier, L.; Chavanon, E.; Lauret, N.; Guilleux, J.; et al. DART: A 3D Radiative Transfer Model for Urban Studies. In *Proceedings of the 2023 Joint Urban Remote Sensing Event (JURSE)*, Heraklion, Greece, 17–19 May 2023; pp. 1–4.

41. Landier, L.; Gastellu-Etchegorry, J.P.; Al Bitar, A.; Chavanon, E.; Lauret, N.; Feigenwinter, C.; Mitraka, Z.; Chrysoulakis, N. Calibration of Urban Canopies Albedo and 3D Shortwave Radiative Budget Using Remote-Sensing Data and the DART Model. *Eur. J. Remote Sens.* **2018**, *51*, 739–753. [CrossRef]

42. Makhloufi, A.; Kallel, A.; Chaker, R.; Gastellu-Etchegorry, J.-P. Retrieval of Olive Tree Biophysical Properties from Sentinel-2 Time Series Based on Physical Modelling and Machine Learning Technique. *Int. J. Remote Sens.* **2021**, *42*, 8542–8571. [CrossRef]

43. Stewart, I.D.; Oke, T.R. Local Climate Zones for Urban Temperature Studies. *Bull. Am. Meteorol. Soc.* **2012**, *93*, 1879–1900. [CrossRef]

44. Aslam, A.; Rana, I.A. The Use of Local Climate Zones in the Urban Environment: A Systematic Review of Data Sources, Methods, and Themes. *Urban Clim.* **2022**, *42*, 101120. [CrossRef]

45. Breiman, L. Random Forests. *Mach. Learn.* **2001**, *45*, 5–32. [CrossRef]

46. Rasmussen, C.E. Gaussian Processes in Machine Learning. In *Advanced Lectures on Machine Learning: ML Summer Schools 2003, Canberra, Australia, 2–14 February 2003, Tübingen, Germany, 4–16 August 2003*; Revised Lectures; Bousquet, O., von Luxburg, U., Rätsch, G., Eds.; Springer: Berlin/Heidelberg, Germany, 2004; pp. 63–71, ISBN 978-3-540-28650-9.

47. INSEE (French National Institute of Statistics and Economic Studies). *Population Census 2020*; INSEE: Montrouge, France, 2022.

48. Haut Conseil pour le Climat en Bretagne. Le changement Climatique En Bretagne. *Bulletin*, 2023; pp. 10–11.

49. Casa, R.; Castaldi, F.; Pascucci, S.; Pignatti, S. Chlorophyll Estimation in Field Crops: An Assessment of Handheld Leaf Meters and Spectral Reflectance Measurements. *J. Agric. Sci.* **2015**, *153*, 876–890. [CrossRef]

50. Wei, S.; Yin, T.; Dissegna, M.A.; Whittle, A.J.; Ow, G.L.F.; Yusof, M.L.M.; Lauret, N.; Gastellu-Etchegorry, J.-P. An Assessment Study of Three Indirect Methods for Estimating Leaf Area Density and Leaf Area Index of Individual Trees. *Agric. For. Meteorol.* **2020**, *292–293*, 108101. [CrossRef]

51. Sain, T.L.; Nabucet, J.; Sulmon, C.; Pellen, J.; Adeline, K.; Hubert-moy, L. A Spatio-Temporal Dataset for Ecophysiological Monitoring of Urban Trees. *Data Brief* **2024**, *57*, 111010. [CrossRef]

52. Copernicus Browser. Available online: https://browser.dataspace.copernicus.eu/ (accessed on 7 October 2024).
53. Stumpf, A.; Michéa, D.; Malet, J.-P. Improved Co-Registration of Sentinel-2 and Landsat-8 Imagery for Earth Surface Motion Measurements. *Remote Sens.* **2018**, *10*, 160. [CrossRef]

54. Zhao, C.; Weng, Q.; Wang, Y.; Hu, Z.; Wu, C. Use of Local Climate Zones to Assess the Spatiotemporal Variations of Urban Vegetation Phenology in Austin, Texas, USA. *GIScience Remote Sens.* **2022**, *59*, 393–409. [CrossRef]

55. Demuzere, M.; Bechtel, B.; Middel, A.; Mills, G. Mapping Europe into Local Climate Zones. *PLOS ONE* **2019**, *14*, e0214474. [CrossRef]

56. Gascon, F.; Gastellu-Etchegorry, J.-P.; Lefevre-Fonollosa, M.-J.; Dufrene, E. Retrieval of Forest Biophysical Variables by Inverting a 3-D Radiative Transfer Model and Using High and Very High Resolution Imagery. *Int. J. Remote Sens.* **2004**, *25*, 5601–5616. [CrossRef]

57. Miraglio, T.; Adeline, K.; Huesca, M.; Ustin, S.; Briottet, X. Joint Use of PROSAIL and DART for Fast LUT Building: Application to Gap Fraction and Leaf Biochemistry Estimations over Sparse Oak Stands. *Remote Sens.* **2020**, *12*, 2925. [CrossRef]

58. Feret, J.-B.; François, C.; Asner, G.P.; Gitelson, A.A.; Martin, R.E.; Bidel, L.P.R.; Ustin, S.L.; le Maire, G.; Jacquemoud, S. PROSPECT-4 and 5: Advances in the Leaf Optical Properties Model Separating Photosynthetic Pigments. *Remote Sens. Environ.* **2008**, *112*, 3030–3043. [CrossRef]

59. Baudin, M.; Dutfoy, A.; Iooss, B.; Popelin, A.-L. OpenTURNS: An Industrial Software for Uncertainty Quantification in Simulation. In *Handbook of Uncertainty Quantification*; Ghanem, R., Higdon, D., Owhadi, H., Eds.; Springer International Publishing: Cham, Switzerland, 2017; pp. 2001–2038, ISBN 978-3-319-12385-1.

60. Conrad, O.; Bechtel, B.; Bock, M.; Dietrich, H.; Fischer, E.; Gerlitz, L.; Wehberg, J.; Wichmann, V.; Böhner, J. System for Automated Geoscientific Analyses (SAGA) v. 2.1.4. *Geosci. Model Dev.* **2015**, *8*, 1991–2007. [CrossRef]

61. Bendig, J.; Yu, K.; Aasen, H.; Bolten, A.; Bennertz, S.; Broscheit, J.; Gnyp, M.L.; Bareth, G. Combining UAV-Based Plant Height from Crop Surface Models, Visible, and near Infrared Vegetation Indices for Biomass Monitoring in Barley. *Int. J. Appl. Earth Obs. Geoinf.* **2015**, *39*, 79–87. [CrossRef]

62. Louhaichi, M.; Borman, M.M.; Johnson, D.E. Spatially Located Platform and Aerial Photography for Documentation of Grazing Impacts on Wheat. *Geocarto Int.* **2001**, *16*, 65–70. [CrossRef]

63. Woebbecke, D.M.; Meyer, G.E.; Bargen, V.; Mortensen, D.A. Color Indices for Weed Identification Under Various Soil, Residue, and Lighting Conditions. *Trans. ASAE* **1995**, *38*, 259–269. [CrossRef]

64. Penuelas, J.; Baret, F.; Filella, I. Semi-Empirical Indices to Assess Carotenoids/Chlorophyll a Ratio from Leaf Spectral Reflectance. *Photosynthetica* **1995**, *31*, 221.

65. Peñuelas, J.; Inoue, Y. Reflectance Indices Indicative of Changes in Water and Pigment Contents of Peanut and Wheat Leaves. *Photosynthetica* **1999**, *36*, 355–360. [CrossRef]

66. Rouse, J.W.; Haas, R.H.; Schell, J.A.; Deering, D.W. Monitoring Vegetation Systems in the Great Plains with ERTS. In *Proceedings of the NASA Scientific and Technical publications, Greenbelt, MD, USA*, 1 January 1974.

67. Kaufman, Y.J.; Tanre, D. Atmospherically Resistant Vegetation Index (ARVI) for EOS-MODIS. *IEEE Trans. Geosci. Remote Sens.* **1992**, *30*, 261–270. [CrossRef]

68. Liu, H.Q.; Huete, A. A Feedback Based Modification of the NDVI to Minimize Canopy Background and Atmospheric Noise. *IEEE Trans. Geosci. Remote Sens.* **1995**, *33*, 457–465. [CrossRef]

69. Son, N.-T.; Chen, C.-F.; Chen, C.-R.; Guo, H.-Y. Classification of Multitemporal Sentinel-2 Data for Field-Level Monitoring of Rice Cropping Practices in Taiwan. *Adv. Space Res.* **2020**, *65*, 1910–1921. [CrossRef]

70. Rondeaux, G.; Steven, M.; Baret, F. Optimization of Soil-Adjusted Vegetation Indices. *Remote Sens. Environ.* **1996**, *55*, 95–107. [CrossRef]

71. Haboudane, D.; Miller, J.R.; Pattey, E.; Zarco-Tejada, P.J.; Strachan, I.B. Hyperspectral Vegetation Indices and Novel Algorithms for Predicting Green LAI of Crop Canopies: Modeling and Validation in the Context of Precision Agriculture. *Remote Sens. Environ.* **2004**, *90*, 337–352. [CrossRef]

72. Gitelson, A.; Merzlyak, M.N. Quantitative Estimation of Chlorophyll-a Using Reflectance Spectra: Experiments with Autumn Chestnut and Maple Leaves. *J. Photochem. Photobiol. B Biol.* **1994**, *22*, 247–252. [CrossRef]

73. Pasqualotto, N.; Delegido, J.; Van Wittenberghe, S.; Rinaldi, M.; Moreno, J. Multi-Crop Green LAI Estimation with a New Simple Sentinel-2 LAI Index (SeLI). *Sensors* **2019**, *19*, 904. [CrossRef] [PubMed]

74. Lymburner, L. Estimation of Canopy-Average Surface-Specific Leaf Area Using Landsat TM Data. *Photogramm. Eng. Remote Sens.* **2000**, *66*, 183–192.

75. Haboudane, D.; Miller, J.R.; Tremblay, N.; Zarco-Tejada, P.J.; Dextraze, L. Integrated Narrow-Band Vegetation Indices for Prediction of Crop Chlorophyll Content for Application to Precision Agriculture. *Remote Sens. Environ.* **2002**, *81*, 416–426. [CrossRef]

76. Wu, C.; Niu, Z.; Tang, Q.; Huang, W. Estimating Chlorophyll Content from Hyperspectral Vegetation Indices: Modeling and Validation. *Agric. For. Meteorol.* **2008**, *148*, 1230–1241. [CrossRef]

77. Qian, B.; Ye, H.; Huang, W.; Xie, Q.; Pan, Y.; Xing, N.; Ren, Y.; Guo, A.; Jiao, Q.; Lan, Y. A Sentinel-2-Based Triangular Vegetation Index for Chlorophyll Content Estimation. *Agric. For. Meteorol.* **2022**, *322*, 109000. [CrossRef]

78. Shi, T.; Xu, H. Derivation of Tasseled Cap Transformation Coefficients for Sentinel-2 MSI At-Sensor Reflectance Data. *IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens.* **2019**, *12*, 4038–4048. [CrossRef]
79. Fang, F.; Im, J.; Lee, J.; Kim, K. An Improved Tree Crown Delineation Method Based on Live Crown Ratios from Airborne LiDAR Data. *GIScience Remote Sens.* **2016**, *53*, 402–419. [CrossRef]

80. Dyer, M.E.; Burkhart, H.E. Compatible Crown Ratio and Crown Height Models. *Can. J. For. Res.* **1987**, *17*, 572–574. [CrossRef]

81. Holdaway, M.R. Modeling Tree Crown Ratio. *For. Chron.* **1986**, *62*, 451–455. [CrossRef]

82. Guo, A.; Ye, H.; Li, G.; Zhang, B.; Huang, W.; Jiao, Q.; Qian, B.; Luo, P. Evaluation of Hybrid Models for Maize Chlorophyll Retrieval Using Medium- and High-Spatial-Resolution Satellite Images. *Remote Sens.* **2023**, *15*, 1784. [CrossRef]

83. Garrigues, S.; Lacaze, R.; Baret, F.; Morisette, J.T.; Weiss, M.; Nickeson, J.E.; Fernandes, R.; Plummer, S.; Shabanov, N.V.; Myneni, R.B.; et al. Validation and Intercomparison of Global Leaf Area Index Products Derived from Remote Sensing Data. *J. Geophys. Res. Biogeosciences* **2008**, *113*. [CrossRef]

84. Misra, G.; Cawkwell, F.; Wingler, A. Status of Phenological Research Using Sentinel-2 Data: A Review. *Remote Sens.* **2020**, *12*, 2760. [CrossRef]

85. Montero, P.; Vilar, J.A. TSclust: An R Package for Time Series Clustering. *J. Stat. Softw.* **2015**, *62*, 1–43. [CrossRef]

86. Golay, X.; Kollias, S.; Stoll, G.; Meier, D.; Valavanis, A.; Boesiger, P. A New Correlation-Based Fuzzy Logic Clustering Algorithm for FMRI. *Magn. Reson. Med.* **1998**, *40*, 249–260. [CrossRef]

87. Douzal, A.; Nagabhushan, P. Adaptive Dissimilarity Index for Measuring Time Series Proximity. *Adv. Data Anal. Classif.* **2007**, *1*, 5–21. [CrossRef]

88. Delegido, J.; Van Wittenberghe, S.; Verrelst, J.; Ortiz, V.; Veroustraete, F.; Valcke, R.; Samson, R.; Rivera, J.P.; Tenjo, C.; Moreno, J. Chlorophyll Content Mapping of Urban Vegetation in the City of Valencia Based on the Hyperspectral NAOC Index. *Ecol. Indic.* **2014**, *40*, 34–42. [CrossRef]

89. Degerickx, J.; Roberts, D.A.; McFadden, J.P.; Hermy, M.; Somers, B. Urban Tree Health Assessment Using Airborne Hyperspectral and LiDAR Imagery. *Int. J. Appl. Earth Obs. Geoinf.* **2018**, *73*, 26–38. [CrossRef]

90. Wu, K.; Chen, J.; Yang, H.; Yang, Y.; Hu, Z. Spatiotemporal Variations in the Sensitivity of Vegetation Growth to Typical Climate Factors on the Qinghai–Tibet Plateau. *Remote Sens.* **2023**, *15*, 2355. [CrossRef]

91. Zhen, Z.; Gastellu-Etchegorry, J.-P.; Chen, S.; Yin, T.; Chavanon, E.; Lauret, N.; Guilleux, J. Quantitative Analysis of DART Calibration Accuracy for Retrieving Spectral Signatures Over Urban Area. *IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens.* **2021**, *14*, 10057–10068. [CrossRef]

92. Adeline, K.R.M.; Paparoditis, N.; Briottet, X.; Gastellu-Etchegorry, J.-P. Material Reflectance Retrieval in Urban Tree Shadows with Physics-Based Empirical Atmospheric Correction. In *Proceedings of the Joint Urban Remote Sensing Event 2013*, Sao Paulo, Brazil, 21–23 April 2013; pp. 279–283.

93. Dissegna, M.A.; Yin, T.; Wu, H.; Lauret, N.; Wei, S.; Gastellu-Etchegorry, J.-P.; Grêt-Regamey, A. Modeling Mean Radiant Temperature Distribution in Urban Landscapes Using DART. *Remote Sens.* **2021**, *13*, 1443. [CrossRef]

94. Zhang, Q.; Wang, D.; Gastellu-Etchegorry, J.-P.; Yang, J.; Qian, Y. Impact of 3-D Structures on Directional Effective Emissivity in Urban Areas Based on DART Model. *Build. Environ.* **2023**, *239*, 110410. [CrossRef]

95. Houborg, R.; McCabe, M.; Cescatti, A.; Gao, F.; Schull, M.; Gitelson, A. Joint Leaf Chlorophyll Content and Leaf Area Index Retrieval from Landsat Data Using a Regularized Model Inversion System (REGFLEC). *Remote Sens. Environ.* **2015**, *159*, 203–221. [CrossRef]

96. Houborg, R.; McCabe, M.F. A Hybrid Training Approach for Leaf Area Index Estimation via Cubist and Random Forests Machine-Learning. *ISPRS J. Photogramm. Remote Sens.* **2018**, *135*, 173–188. [CrossRef]

97. Zurita-Milla, R.; Laurent, V.C.E.; van Gijsel, J.A.E. Visualizing the Ill-Posedness of the Inversion of a Canopy Radiative Transfer Model: A Case Study for Sentinel-2. *Int. J. Appl. Earth Obs. Geoinf.* **2015**, *43*, 7–18. [CrossRef]

98. Baret, F.; Buis, S. Estimating Canopy Characteristics from Remote Sensing Observations: Review of Methods and Associated Problems. In *Advances in Land Remote Sensing: System, Modeling, Inversion and Application*; Liang, S., Ed.; Springer: Dordrecht, The Netherlands, 2008; pp. 173–201, ISBN 978-1-4020-6450-0.

99. Schiefer, F.; Schmidtlein, S.; Kattenborn, T. The Retrieval of Plant Functional Traits from Canopy Spectra through RTM-Inversions and Statistical Models Are Both Critically Affected by Plant Phenology. *Ecol. Indic.* **2021**, *121*, 107062. [CrossRef]

100. Fernández-Guisuraga, J.M.; Suárez-Seoane, S.; Quintano, C.; Fernández-Manso, A.; Calvo, L. Comparison of Physical-Based Models to Measure Forest Resilience to Fire as a Function of Burn Severity. *Remote Sens.* **2022**, *14*, 5138. [CrossRef]

101. Atzberger, C.; Richter, K. Spatially Constrained Inversion of Radiative Transfer Models for Improved LAI Mapping from Future Sentinel-2 Imagery. *Remote Sens. Environ.* **2012**, *120*, 208–218. [CrossRef]

102. Rivera, J.P.; Verrelst, J.; Leonenko, G.; Moreno, J. Multiple Cost Functions and Regularization Options for Improved Retrieval of Leaf Chlorophyll Content and LAI through Inversion of the PROSAIL Model. *Remote Sens.* **2013**, *5*, 3280–3304. [CrossRef]

103. Liu, S.; Brandt, M.; Nord-Larsen, T.; Chave, J.; Reiner, F.; Lang, N.; Tong, X.; Ciais, P.; Igel, C.; Pascual, A.; et al. The Overlooked Contribution of Trees Outside Forests to Tree Cover and Woody Biomass across Europe. *Sci. Adv.* **2023**, *9*, eadh4097. [CrossRef] [PubMed]

104. Wright, M.N.; Ziegler, A.; König, I.R. Do Little Interactions Get Lost in Dark Random Forests? *BMC Bioinform.* **2016**, *17*, 145. [CrossRef] [PubMed]

105. Atmosphere, U.S. *US Standard Atmosphere*; National Oceanic and Atmospheric Administration: Silver Spring, MD, USA, 1976.
106. Yu, K.; Van Geel, M.; Ceulemans, T.; Geerts, W.; Ramos, M.M.; Sousa, N.; Castro, P.M.L.; Kastendeuch, P.; Najjar, G.; Ameglio, T.; et al. Foliar Optical Traits Indicate That Sealed Planting Conditions Negatively Affect Urban Tree Health. *Ecol. Indic.* **2018**, *95*, 895–906. [CrossRef]

107. Chianucci, F.; Pisek, J.; Raabe, K.; Marchino, L.; Ferrara, C.; Corona, P. A Dataset of Leaf Inclination Angles for Temperate and Boreal Broadleaf Woody Species. *Ann. For. Sci.* **2018**, *75*, 50. [CrossRef]

**Disclaimer/Publisher's Note:** The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.