IoT Monitoring of Urban Tree Ecosystem Services: Possibilities and Challenges

: Urban green infrastructure plays an increasingly signiﬁcant role in sustainable urban development planning as it provides important regulating and cultural ecosystem services. Monitoring of such dynamic and complex systems requires technological solutions which provide easy data collection, processing, and utilization at a ﬀ ordable costs. To meet these challenges a pilot study was conducted using a network of wireless, low cost, and multiparameter monitoring devices, which operate using Internet of Things (IoT) technology, to provide real-time monitoring of regulatory ecosystem services in the form of meaningful indicators for both human health and environmental policies. The pilot study was set in a green area situated in the center of Moscow, which is exposed to the heat island e ﬀ ect as well as high levels of anthropogenic pressure. Sixteen IoT devices were installed on individual trees to monitor their ecophysiological parameters from 1 July to 31 November 2019 with a time resolution of 1.5 h. These parameters were used as input variables to quantify indicators of ecosystem services related to climate, air quality, and water regulation. Our results showed that the average tree in the study area during the investigated period reduced extreme heat by 2 ◦ C via shading, cooled the surrounding area by transferring 2167 ± 181 KWh of incoming solar energy into latent heat, transpired 137 ± 49 mm of water, sequestered 8.61 ± 1.25 kg of atmospheric carbon, and removed 5.3 ± 0.8 kg of particulate matter (PM 10 ). The values of the monitored processes varied spatially and temporally when considering di ﬀ erent tree species (up to ﬁve to ten times), local environmental conditions, and seasonal weather. Thus, it is important to use real-time monitoring data to deepen understandings of the processes of urban forests. There is a new opportunity of applying IoT technology not only to measure trees functionality through ﬂuxes presents the results of a pilot study conducted using a network of and monitoring to monitor This study aimed

This work presents the results of a pilot study which was conducted in a green area situated in the center of Moscow using a network of wireless, low cost, and multiparameter monitoring devices TreeTalkers (TT+) [36], to monitor single-tree ecophysiological parameters. This study aimed at testing the feasibility of the proposed technological solution to provide real-time monitoring of regulatory ecosystem services, reported in the form of meaningful indicators for potential end-users and relevant for both human health and environmental policy targets. This is particularly relevant since the New Moscow Development Project, adopted in 2012, is continuing to expand UGI as result of active urbanization on an area of more than 1500 km 2 , the impact of which on soils and ecosystems services is already visible [37]. As such, its monitoring and operational management is greatly required.

Study Site and Network Setup
With a population of over 12 million people, Moscow is the largest metropolis in Europe. The urban vegetation is made of species characteristic of the south taiga zone [38], mainly conserved in urban parks. It also includes introduced genera such as Tilia, Acer, Salix, and other ornamental forms of trees and shrubs. According to the Köppen climate classification, the climate in this area is a humid continental climate. The study site is a greened public place (green area) created in 1948, located in the Bolotnaya square in central Moscow 600 m to the south of the Kremlin, on Balchug Island. The area is under the higher influence of urban heat (mean annual air temperature is +1.8 • C than the rural reference- [39,40]) with a high level of anthropogenic pressure. The area covers about 4.5 hectares (Figure 1). TreeTalkers (TT+) [36], to monitor single-tree ecophysiological parameters. This study aimed at testing the feasibility of the proposed technological solution to provide real-time monitoring of regulatory ecosystem services, reported in the form of meaningful indicators for potential end-users and relevant for both human health and environmental policy targets. This is particularly relevant since the New Moscow Development Project, adopted in 2012, is continuing to expand UGI as result of active urbanization on an area of more than 1500 km 2 , the impact of which on soils and ecosystems services is already visible [37]. As such, its monitoring and operational management is greatly required.

Study Site and Network Setup
With a population of over 12 million people, Moscow is the largest metropolis in Europe. The urban vegetation is made of species characteristic of the south taiga zone [38], mainly conserved in urban parks. It also includes introduced genera such as Tilia, Acer, Salix, and other ornamental forms of trees and shrubs. According to the Köppen climate classification, the climate in this area is a humid continental climate. The study site is a greened public place (green area) created in 1948, located in the Bolotnaya square in central Moscow 600 m to the south of the Kremlin, on Balchug Island. The area is under the higher influence of urban heat (mean annual air temperature is +1.8 °C than the rural reference- [39,40]) with a high level of anthropogenic pressure. The area covers about 4.5 hectares ( Figure 1). The nearby Balchug meteorological station (700 meters distance) provided meteorological data at a 3 h frequency. The daily mean temperature varied between −6.7 °C in February and +19.2 °C in July, and the amount of precipitation was about 700 mm/year, mainly in the summer months.
TreeTalker+ devices (TT+) are microprocessor-based IoT platforms built around the ATMEGA328p (Atmel Corp) chip, equipped with a LoRa transceiver for radio transmission to a central The nearby Balchug meteorological station (700 meters distance) provided meteorological data at a 3 h frequency. The daily mean temperature varied between −6.7 • C in February and +19.2 • C in July, and the amount of precipitation was about 700 mm/year, mainly in the summer months.
TreeTalker+ devices (TT+) are microprocessor-based IoT platforms built around the ATMEGA328p (Atmel Corp) chip, equipped with a LoRa transceiver for radio transmission to a central gateway, which collects the individual tree data and send it to the cloud using GSM/GPRS technology. The TT+ were first developed and tested by our team [36,41,42] and later, based on our design, industrialized by the Nature4 Benefit Corporation (www.nature4.org) based in Italy, which is a start-up company dedicated to environmental monitoring. The TT+ sensors are able to measure (1) the sap-flow density, using the transient thermal dissipation method based on a heating/cooling cycle of 10 min every 1.5 h [43]; (2) the light transmission spectra through the canopy in 12 spectral bands, using two spectrometers (VIS and NIR); (3) diameter growth with an optical IR-pulsed device; and (4) stem position and oscillation in three axes with an on-board accelerometer. In addition, air temperature and humidity are recorded at the single-tree level. Specifications and pictures are listed in Table 1 and Figure 2, respectively. gateway, which collects the individual tree data and send it to the cloud using GSM/GPRS technology. The TT+ were first developed and tested by our team [36,41,42] and later, based on our design, industrialized by the Nature4 Benefit Corporation (www.nature4.org) based in Italy, which is a startup company dedicated to environmental monitoring. The TT+ sensors are able to measure (1) the sapflow density, using the transient thermal dissipation method based on a heating/cooling cycle of 10 min every 1.5 h [43]; (2) the light transmission spectra through the canopy in 12 spectral bands, using two spectrometers (VIS and NIR); (3) diameter growth with an optical IR-pulsed device; and (4) stem position and oscillation in three axes with an on-board accelerometer. In addition, air temperature and humidity are recorded at the single-tree level. Specifications and pictures are listed in Table 1 and Figure  2, respectively.  Sixteen TreeTalker+ (TT+) devices were installed on corresponding trees of four different tree species, five on Acer platanoides (average diameter at breast height (DBH) 38.7 cm), three on Betula pendula (average DBH 21.8 cm), three on Larix sibirica (average DBH 32.1 cm), and five on Tilia cordata (average DBH 34.1 cm). At the onset of installation, all trees were characterized by height, diameter, age, standing type (on the edge of the square or inside it), and VTA (visual tree assessment) [44]. The investigated trees' VTA scores ranged from 1 to 3; from a scale of 1 (healthy condition) to 7 (severe decline), (see Appendix A). A reference device TT+ (TT-R) was mounted outside of the tree canopies Sixteen TreeTalker+ (TT+) devices were installed on corresponding trees of four different tree species, five on Acer platanoides (average diameter at breast height (DBH) 38.7 cm), three on Betula pendula (average DBH 21.8 cm), three on Larix sibirica (average DBH 32.1 cm), and five on Tilia cordata (average DBH 34.1 cm). At the onset of installation, all trees were characterized by height, diameter, age, standing type (on the edge of the square or inside it), and VTA (visual tree assessment) [44]. The investigated trees' VTA scores ranged from 1 to 3; from a scale of 1 (healthy condition) to 7 (severe decline), (see Appendix A). A reference device TT+ (TT-R) was mounted outside of the tree canopies to collect climate data and incoming solar radiation reference spectra. On individual trees, the devices were placed at a height of 3 m from the ground on the north side of the trunk, and the solar-powered batteries, on the south side. The 3 m height was chosen to reduce the risk of damage or theft of devices. All measurements were conducted from 01 July till 31 November 2019 with 1.5 h temporal resolution.

Choice of ES Indicators
A wide range of existing ES indicators have been investigated in previous studies [11,45]. For the present study, the ES indicators were chosen on the basis of the possibility to be directly estimated by our measured parameters. The ES indicators and the relative measured variables, including algorithm references are shown on Table 2. Direct measurements (like air temperature or relative humidity) gave us an opportunity to calculate ESI without any additional assumptions. However, all the ESI that are labeled as "indirect" required us to introduce some adjustment factors or assumptions. One of the most frequently-used "indirect" predictors of ES is tree leaf area index (LAI) since many ESI can be estimated by the amount of leaf area [11,45]. For the purpose of this paper we monitored LAI in real-time by using the two on-board spectrometers and we present particulate adsorption (calculated based on iTree model's equations [46]) as an example of usefulness of LAI indicator indirect use.
where ∆V is tree stem volume (m 3 ); BCEF is the biomass conversion and expansion factor, obtained using this formula BCEF = BEF * ρ, where BEF is the biomass expansion factor to obtain the total above-ground biomass and ρ is the density of the wood (kg/m 3 ); R is the ratio of below-ground to above-ground biomass; and CF is the carbon fraction of the dry mass (conventionally equal to 0.5 as suggested by the IPCC [69]. The BCEF and R were taken from literature according to species and age of the tree [70] and CF (biomass conversion into carbon) was taken as 0.5. Trunk volume (V) was calculated using the height measured directly in field and basal area increments (BAI according to [71,72]). where w is the trunk diameter expansion measured in real-time with the TreeTalker+ IR distance sensor. Stem dendrometric coefficients for different species were taken into account for the estimation of the trunk shape in the final calculation of biomass increment [69,70]. Uncertainties in carbon sequestration estimates arising from measurements errors and the uncertainty associated to the parameters retrieved from literature were assessed performing error propagation analysis by means of the "errors" R package [73]. Errors associated to tree height, stem circumference at breast height, and stem radial growth measurements were considered as 0.5 and 0.01 m, respectively. Stem radial growth error was the accuracy of the TreeTalker+ IR distance sensor (Table 1) while uncertainties in BCEF and R were retrieved from [70].

Climate Regulation via Air Temperature Control
We used direct measurements of temperature changes based on the difference between data from thermo-hygrometers of individual TT, measuring climate parameters at 3 m height under the crown space, and TT-R (reference outside station). Temperature measurements were taken by the TT device with a measuring cycle of 1.5 h as for all the other parameters.

Water Fluxes and Energy Consumption through Transpiration
The transpiration rate of whole plants is closely approximated by the sap-flow rate in the main stem or trunk. We implemented the thermal dissipation method in the TT+ platform, first developed by Granier [74]. Granier's original method is based on the assessment of the amount of heat dissipated from a heated probe in relation to a reference probe in proximity (about 10 cm in our case), with continuous heating. To reduce energy consumption, the method was modified according to Do et al. [43,75] and instead of continuous heating, cycles of cooling (80 min) and heating (10 min) were used. The heated and reference probes of TT+ were installed into the tree trunk at a vertical distance of 10 cm. The probes had a diameter of 2.5 mm each and were installed at the depth of 2.5 cm. The TT+ heat dissipation probe was installed at the 3.5 m height of the trunk and was well protected from direct sun heating by the canopy and the northernly orientation. The sap wood area for each tree was assessed by using literature data on its relation to tree diameter [76][77][78]. Sap flow was calculated with the assumption that the whole trunk sapwood area was conducting water. It was assumed that daily transpiration was equal to daily sum of sap flow.
The energy absorbed by the tree for transpiration was calculated based on the equation: where λ is the energy spent in transpiration (the latent heat for vaporization of water = 2264.705 KJ/Kg) and T is the transpiration, which is calculated by sap-flow density and the estimated sapwood area. We have considered that the potential runoff mitigation of urban trees as ecosystem service indicators can be shown as the ratio of transpiration (T) to precipitation (P): where both terms are calculated using mm. Tree-level transpiration has been recalculated per unit of crown area. The equation could be, in principle, implemented with evapotranspiration, with evaporation estimated by simple meteorological models and canopy interception through LAI, but in the present paper we simplify the approach to derive a runoff mitigation indicator, represented by R p .

LAI
There are many methods and protocols to estimate LAI [79,80]. Light transmission through the canopy as porous media can be treated according to Beer's Law [81]. In this way LAI can be estimated by the extinction of photosynthetic light radiation through the canopy [82]. Photosynthetically-active radiation was measured outside and below the canopy (with TT-R and TT+ spectrometers, respectively). Since the light is also blocked by the woody components of canopy (i.e., branches and twigs), the extinction of the light profile gives the PAI (plant area index) as follows: PAI consists of wood area index (WAI) and leaf area index (LAI). Assuming that WAI is constant throughout the vegetation period and LAI = 0 after defoliation (second part of October and November), WAI for each tree was calculated as the mean PAI of November. Light extinction coefficient k was calculated per each species utilizing LAI light measurements at the foliage peak with a hemispherical photo done with Kodak PIXPRO SP360 and Hemisfer software [83].

Particulate Adsorption
Dry deposition of solid particles on the canopy was calculated according to i-Tree Eco dry deposition model [46]: where C is the PM 10 concentration (g m −3 ), V d the velocity of deposition (m s −1 ), LAI PM10 the leaf area index for pollutant deposition (m 2 m −2 ), and V dPM10max , V dPM10avg , and V dPM10min are maximum, average, and minimum deposition velocity for PM 10 -0.0064 ms −1 , 0.0025 ms −1 , and 0.01 ms −1 , respectively. Values were multiplied by canopy area to show adsorption per individual tree. Pollutant concentration was obtained from nearby open access PM 10 sensors via the sensor.community web portal in Moscow (https://sensor.community).

Data Processing
Data collection with TreeTalker+ devices was organized according to the following scheme. All types of devices (TT+ and TT-R) made measurement every 90 min, stored data in internal memory, and then according to a predefined time window transmitted data to the gateway (TT-cloud) device via low-power wide-area network (LoRa) wireless networking protocol. The TT-cloud device is a gateway, whose purpose is to collect data from all TT devices on site, store it and then transmit it to an online database via WiFi or wireless mobile networks.
All remote data were collected and processed using the R computing language [84]. Field data was organized in a table and added to computation during early stages of processing. All measurement which deviated by more than three standard deviations from the weekly mean was filtered out. Filtered data was linearly interpolated. Most of data gaps were less than four measurements (6 h) and the maximum gap was 44 measurements (66 h) due to water-damaged electronics after heavy rain. Data were filled with data from trees with the closest parameters (species, trunk diameter, height, canopy size, and position on site). Statistical analysis of results (correlation analysis) was also performed in R environment.

Carbon Sequestration
As the diameter increment is an important component of carbon storage across the season, relative biomass accumulation was calculated using Equation (1). The growth dynamics show biomass increase until the end of September (around day 260 of the year) due to the warm autumn ( Figure 3).

Carbon Sequestration
As the diameter increment is an important component of carbon storage across the season, relative biomass accumulation was calculated using Equation (1). The growth dynamics show biomass increase until the end of September (around day 260 of the year) due to the warm autumn ( Figure 3). Betula pendula growth rates ranged from 2.19 to 2.79 kg C per tree along the season and they represent the lowest values among the investigated trees. Acer platanoides growth rates ranged from 8.35 to 15.3 kg C per tree along the season. Tilia cordata growth rates were slightly lower than those for Acer platanoides, ranging from 1.76 to 6.27 kg C, showing a marked individual variation. Larix sibirica trees showed similar growth rates with values ranging from 4.95 to 7.66 kg C per tree. In total the average amount of carbon accumulated in the tree biomass for the investigated species was 11.58 ± 1.47 kg C (Acer platanoides, ± standard error), 2.5 ± 0.15 kg C (Betula pendula), 6.03 ± 1.31 kg C (Larix sibirica), and 10.89 ± 2.89 kg C (Tilia cordata). The overall uncertainty in the biomass carbon sequestration estimates arising from the propagation of measurements and parameter uncertainties, was on average 3.6% ranging from a minimum of 2.3% (Larix spp.) to a maximum of 5.6% (Betula spp.). Betula pendula growth rates ranged from 2.19 to 2.79 kg C per tree along the season and they represent the lowest values among the investigated trees. Acer platanoides growth rates ranged from 8.35 to 15.3 kg C per tree along the season. Tilia cordata growth rates were slightly lower than those for Acer platanoides, ranging from 1.76 to 6.27 kg C, showing a marked individual variation. Larix sibirica trees showed similar growth rates with values ranging from 4.95 to 7.66 kg C per tree. In total the average amount of carbon accumulated in the tree biomass for the investigated species was 11.58 ± 1.47 kg C (Acer platanoides, ± standard error), 2.5 ± 0.15 kg C (Betula pendula), 6.03 ± 1.31 kg C (Larix sibirica), and 10.89 ± 2.89 kg C (Tilia cordata). The overall uncertainty in the biomass carbon sequestration estimates arising from the propagation of measurements and parameter uncertainties, was on average 3.6% ranging from a minimum of 2.3% (Larix spp.) to a maximum of 5.6% (Betula spp.). In general carbon sequestration related to resultant biomass increment significantly was correlated mostly to tree stem diameter, exhibiting a Pearson r coefficient of 0.87, followed by canopy size and tree transpiration with mutually similar, yet much weaker, relations (r = 0.39 and 0.36, respectively).
General knowledge of urban tree-carbon dynamics, including the balance of growth, mortality, and planting rates is quite data-limited [47]. Although the carbon stock density of urban green areas is generally smaller than that of forests, the potential storage capacity is still considerably higher [85,86]. Growth rates of urban trees may be accelerated by the heat island effect, as a result of increased temperature, longer growing season, and potentially higher N deposition [87]. Conversion of diameter growth in biomass and hence carbon sequestration is very much related to the BEF coefficients which we derived by [70], specifically for the Russian environment. One of the most comprehensive studies on urban tree-carbon sequestration is one across a wide range of US cities [47]. These authors showed an average annual net carbon uptake, per tree, of about 0.226 kg C m −2 year −1 . Unfortunately, they did not analyze individual species performance but only the overall carbon accumulated by the unit of green area. Scaling larch trees carbon sequestration per unit of crown area we obtain an average value of 0.14 ± 0.04 kg C m −2 , which is the closest result in relation to Nowak et al. [47], but for broad leaved species this parameter was more than two times higher (with a maximum of 0.48 ± 0.14 kg C m −2 over Forests 2020, 11, 775 9 of 25 the study period for Acer platanoides). For Tilia cordata of the same age, Moser et al. [88] reported a mean annual carbon uptake of 4.58 kg C per tree (2.48-7.12) which is lower compared with our estimate of 10.89 ± 2.89 kg C (4.3-17.79).
However it compares well with our mean annual growth increment of 4.5 ± 1.18, which was determined using the same approach using total stock and tree age. Using a continuous dendrometer with larch species of about the same age (60 years), Deslauriers et al. [89,90] found an annual radial increment of 1.8 mm which is comparable with our estimates of 2.07 ± 0.33 mm. Using a large national inventory for birch trees in Finland, Repola et al. [91], found an annual radial increment of 1.1 mm, while we observed an growth of 2.79 ± 0.91 mm. However, we have to consider the difference between forests and isolated urban trees. This also showed the importance of considering the difference between the current and the mean increment in future carbon sequestration analyses. Although it is difficult to compare our data with existing annual growth rates, due to the lack of full season data coverage (about four out of six months) our data show a reasonable consistency with the annual mean increment of the same trees. More importantly they confirm the possibility of monitoring the carbon sequestration of urban trees in real-time.

Cooling Effect
The local climate control from an ecosystem-services perspective is the mitigation of extreme temperatures and providing a more comfortable urban microclimate. In Figure 4 the diurnal difference mean changes of temperature between the reference station and the space under the tree crown is presented as the mean monthly day, for the investigated species. During the day, the temperature difference was maximal, peaking early afternoon in July, August, and September. This dynamic was reversed in October. During the day trees were cooling the surrounding air via shading and evaporative cooling, showing an effect of up to about 2 • C degrees with the external temperatures (Appendix B). In October, all the species showed a warming effect at midday, of the same order. An opposite behavior is recorded for nighttime periods where, usually during summer months, trees were slightly warmer than the surrounding air. In October they showed a cooling effect during the night.  Another important effect of trees is the mitigation of climate extremes. For this purpose, we have estimated the daily differences between maximum and minimum temperatures of the reference station and the individual tree recording ( Figure 5). While temperature amplitudes outside of the canopy (the black line) reached a maximum of about 10 °C degrees in August, under the canopy this amplitude was 3 °C degrees lower. All the species showed similar dynamics through the investigated Another important effect of trees is the mitigation of climate extremes. For this purpose, we have estimated the daily differences between maximum and minimum temperatures of the reference station and the individual tree recording ( Figure 5). While temperature amplitudes outside of the canopy (the black line) reached a maximum of about 10 • C degrees in August, under the canopy this amplitude was 3 • C degrees lower. All the species showed similar dynamics through the investigated period, but under the larch, temperature extremes were more like the reference station both during summer and particularly after autumn leaf drop. These findings correspond well to [60], who explained observed seasonal dynamics of temperature reduction rates by differences in transpiration and to [92], who also found that within the canopy radius of 4.5 m of Acer platanoides or Tilia cordata trees, daytime air temperature decreased up to 3.5 °C during August in comparison to the unshaded surrounding area. Our results agree with a meta-analysis study, which showed relations between individual tree characteristics and daily/seasonal temperature reduction dynamics in different climate and urban conditions [92]. Air temperature reduction by urban green infrastructure of about 0.5-2.5 °C degrees was also shown in several papers, which used computational modeling [52,93] and satellite based data for land-surface temperature [94,95].

Run-Off Mitigation and Energy Consumption by Trees via Transpiration
Ensembles of transpiration daily patterns show the typical diurnal behavior modulated by changing environmental parameters, with a variation in the onset and the peak of the transpiration cycle across the individuals, particularly in Acer platanoides and Tilia cordata, showing the variation of conditions to sun exposure (Appendix C). Sap flux during night for all individuals was negligible, while during the day it increased up to 3 lh −1 for Betula pendula trees and 5 lh −1 for several Tilia cordata and Acer platanoides trees in July. Maximal values of daily flux were detected for all species in July, and reached 8 lh −1 for Larix sibirica, 7.6 lh −1 for Acer platanoides, 5.9 lh −1 for Tilia cordata, and 3.8 lh −1 for Betula pendula, while average daily (24 h) flux for those species in July was 2.90 ± 0.08, 2.85 ± 0.08, 2.18 ± 0.05, and 1.74 ± 0.04 lh −1 , respectively. In following months, sap flux values gradually decreased for all species. These findings correspond well to [60], who explained observed seasonal dynamics of temperature reduction rates by differences in transpiration and to [92], who also found that within the canopy radius of 4.5 m of Acer platanoides or Tilia cordata trees, daytime air temperature decreased up to 3.5 • C during August in comparison to the unshaded surrounding area. Our results agree with a meta-analysis study, which showed relations between individual tree characteristics and daily/seasonal temperature reduction dynamics in different climate and urban conditions [92]. Air temperature reduction by urban green infrastructure of about 0.5-2.5 • C degrees was also shown in several papers, which used computational modeling [52,93] and satellite based data for land-surface temperature [94,95].

Run-Off Mitigation and Energy Consumption by Trees via Transpiration
Ensembles of transpiration daily patterns show the typical diurnal behavior modulated by changing environmental parameters, with a variation in the onset and the peak of the transpiration cycle across the individuals, particularly in Acer platanoides and Tilia cordata, showing the variation of conditions to sun exposure (Appendix C). Sap flux during night for all individuals was negligible, while during the day it increased up to 3 lh −1 for Betula pendula trees and 5 lh −1 for several Tilia cordata and Acer platanoides trees in July. Maximal values of daily flux were detected for all species in July, and reached 8 lh −1 for Larix sibirica, 7.6 lh −1 for Acer platanoides, 5.9 lh −1 for Tilia cordata, and 3.8 lh −1 for Betula pendula, while average daily (24 h) flux for those species in July was 2.90 ± 0.08, 2.85 ± 0.08, 2.18 ± 0.05, and 1.74 ± 0.04 lh −1 , respectively. In following months, sap flux values gradually decreased for all species.
For the purpose of the current study we report the amount of cumulated transpired water in relation to seasonal rainfall. The units are expressed in mm of water, where for each of the trees, transpiration rates have been converted to mm by using the tree canopy area. The main purpose is to show, from an ecosystem-services perspective, the role of trees to complement (or support) the drainage systems by actively lowering the soil water content, thus increasing its marginal water storage capacity and mitigating the volume of stormwater entering the drainage channels. Our results show that Acer platanoides and Tilia cordata trees consumed 107 ± 39 mm (± standard error) and 122 ± 32 mm of water, respectively, through the investigated period. Larix sibirica trees reduced 66 ± 40 mm of water, while young Betula pendula trees reduced 240 ± 53 mm (Appendix A, Figure 6). In general, this cumulative process was mostly linear with some differences in rates that were found to be associated with different VTA scores. On the other hand, there was no significant response to the heavy rains. According to our data we can say that up to 130 ± 28% (± standard error) of precipitation can be transpired per square meter of canopy area for Betula pendula. Meanwhile mature broadleaf trees such as Tilia cordata and Acer platanoides transpired 66 ± 17% and 58 ± 21% of precipitated water during the study period. The only investigated coniferous species showed minimal values yielding 36 ± 21% of precipitated water transpired. The high values of Betula pendula, although showing a lower rate of transpiration per tree, can be explained by a greater efficiency per canopy area. This also shows how the birch trees can very efficiently utilize soil moisture storage. drainage systems by actively lowering the soil water content, thus increasing its marginal water storage capacity and mitigating the volume of stormwater entering the drainage channels. Our results show that Acer platanoides and Tilia cordata trees consumed 107 ± 39 mm (± standard error) and 122 ± 32 mm of water, respectively, through the investigated period. Larix sibirica trees reduced 66 ± 40 mm of water, while young Betula pendula trees reduced 240 ± 53 mm (Appendix A, Figure 6). In general, this cumulative process was mostly linear with some differences in rates that were found to be associated with different VTA scores. On the other hand, there was no significant response to the heavy rains. According to our data we can say that up to 130 ± 28% (± standard error) of precipitation can be transpired per square meter of canopy area for Betula pendula. Meanwhile mature broadleaf trees such as Tilia cordata and Acer platanoides transpired 66 ± 17% and 58 ± 21% of precipitated water during the study period. The only investigated coniferous species showed minimal values yielding 36 ± 21% of precipitated water transpired. The high values of Betula pendula, although showing a lower rate of transpiration per tree, can be explained by a greater efficiency per canopy area. This also shows how the birch trees can very efficiently utilize soil moisture storage. For the investigated period, the ability of individual trees to transpire precipitated water volume around them positively correlated with the size of those species canopies. This supports the findings of similar works [88,96]. Taking into account all the trees within the study location, the total annual plot rainfall (excluding any extreme rainfall events) was approximately equal to total transpiration, what was shown for boreal urban trees of other municipalities [61,63]. Such information obtained in real-time could considerably contribute to cities' stormwater management [97][98][99]. On the other For the investigated period, the ability of individual trees to transpire precipitated water volume around them positively correlated with the size of those species canopies. This supports the findings of similar works [88,96]. Taking into account all the trees within the study location, the total annual plot rainfall (excluding any extreme rainfall events) was approximately equal to total transpiration, what was shown for boreal urban trees of other municipalities [61,63]. Such information obtained in real-time could considerably contribute to cities' stormwater management [97][98][99]. On the other hand, it is widely discussed that rain interception by leaf buffering during heavy rains also contributes to run-off mitigation [66,100,101], which is based on leaf area index, as discussed further in the text.
Using the energy balance equation (Equation (3)) it was possible to estimate the amount of absorbed radiation that is subtracted from the environment due to transpiration. The diurnal graph shows the increase of the transfer latent heat (L) during the daytime when transpiration starts following the sunrise (Figure 7). The daily range of adsorbed energy during July and August was 0.5-2.2 kWh for Acer platanoides, Larix sibirica, and Tilia cordata trees, while for Betula pendula it was a bit lower (0.2-1.5 kWh). Already in September there was a decrease in absolute numbers below 1.2 kWh for all species, but the range for Acer platanoides was 0.6-1.  Using the energy balance equation (Equation (3)) it was possible to estimate the amount of absorbed radiation that is subtracted from the environment due to transpiration. The diurnal graph shows the increase of the transfer latent heat (L) during the daytime when transpiration starts following the sunrise (Figure 7). The daily range of adsorbed energy during July and August was 0.5-2.2 kWh for Acer platanoides, Larix sibirica, and Tilia cordata trees, while for Betula pendula it was a bit lower (0.2-1.5 kWh). Already in September there was a decrease in absolute numbers below 1.2 kWh for all species, but the range for Acer platanoides was 0.6-1.  Energy latent heat fluxes varied between tree species, individual trees, and seasons. We found that main variability between individual trees was due to the difference in stem size and canopy area. While there are few publications investigating the energy balance of boreal urban trees, we showed comparable measurements for Tilia cordata in summer months to another study [88].  Energy latent heat fluxes varied between tree species, individual trees, and seasons. We found that main variability between individual trees was due to the difference in stem size and canopy area. While there are few publications investigating the energy balance of boreal urban trees, we showed comparable measurements for Tilia cordata in summer months to another study [88].

LAI as a Proxy for ES Provision
The dynamics of LAI are a good proxy for several types of ES, such as wind velocity reduction, noise reduction, pollution regulation, and erosion protection via leaves acting as a buffer (see references from Table 2). From our radiation sensor we can calculate the overall PAI, that is, a sum of the leaves and wood (bark) light interception effects. Thus, after comparing periods with and without leaves we could derive both leaf and wood area index. These periods are clearly visible in Figure 8 where we can easily distinguish the time of defoliation which lasted for one week in the first days of October. Average LAI for all trees (Appendix A) was 3.77 with less variation in Betula pendula (3.5-3.7) and Larix sibirica (3.6-3.9) and higher variation for Acer platanoides   Despite the clear pattern shown in Figure 8, LAI, estimated by the spectrometer, was influenced by the angle of view of the spectrometer, which is ±20° and positioned in the north side of the trunk, thus recording only one section of the crown. For a better absolute determination it is necessary to compare the sensor with more precise techniques (e.g., camera with fish-eye lens or other LAI meters) which are essential for calibration and validation [80,102]. That said, the absolute values (3)(4) are comparable to most of the papers [80,82,103] in similar species. While the strong correlation between LAI and DBH and also between LAI and tree height was shown previously for Canadian and Bavarian boreal forests [103,104] our results did not allow us to confirm this correlation. This lack of correlation is possibly explained by the small number of trees in consideration and the urban conditions of our study site [105]. However the real-time dynamics of LAI and WAI is the advantage of our device which is important for plant physiological processes and models [106]. Despite the clear pattern shown in Figure 8, LAI, estimated by the spectrometer, was influenced by the angle of view of the spectrometer, which is ±20 • and positioned in the north side of the trunk, thus recording only one section of the crown. For a better absolute determination it is necessary to compare the sensor with more precise techniques (e.g., camera with fish-eye lens or other LAI meters) which are essential for calibration and validation [80,102]. That said, the absolute values (3)(4) are comparable to most of the papers [80,82,103] in similar species. While the strong correlation between LAI and DBH and also between LAI and tree height was shown previously for Canadian and Bavarian boreal forests [103,104] our results did not allow us to confirm this correlation. This lack of correlation is possibly explained by the small number of trees in consideration and the urban conditions of our study site [105]. However the real-time dynamics of LAI and WAI is the advantage of our device which is important for plant physiological processes and models [106].

Particulate Adsorption
Particulate adsorption is influenced by two main parameters -PAI and particulate concentration in air, according to the model described in Methods. The adsorption dynamic during the investigated period mostly reflected the changes in particulate matter concentration, which showed peaks at the end of July and beginning of August, and again in September and late November. While the reasons for these pollution dynamics were not our focus, we can say that the average absorption rates differed among the investigated trees, mostly due to their LAI. The lowest adsorption rates were shown by Betula pendula trees with an average of 9.3 g per day, and the highest were Acer platanoides and Larix sibirica trees with 51.2 and 51.4 g per day, respectively. Tilia cordata trees were in the middle with 25.5 g per day per tree (Appendix E). The highest adsorption peaks (in terms of maximum deposition velocity) were up to 1200 g (Larix sibirica) during several days with high concentrations in the air. This parameter showed drastic variation across time. The cumulative particulate mass adsorption over the season was 7.7 ± 1.5 kg for Acer platanoides (± standard error), 7.7 ± 2.1 kg for Larix sibirica, 3.8 ± 0.2 kg for Tilia cordata, and 1.4 ± 0.1 kg for Betula pendula (Figure 9, Appendix A).
Saebø et al. [68] showed that Betula pendula improves air quality by particulate adsorption much better than Acer platanoides and Tilia cordata, however tree size should be taken into account as a factor to interpret the results of our study. It is well known that healthy large trees remove about 60 times more pollution annually than healthy small trees [67]. Additionally, leaf wax or leaf hair density and topography also have great influence [68,107]. In our case we could not include leaf morphology or topography since our estimation was based on LAI. However the total adsorption by tree seems comparable with several works yielding an average of 0.5-5 g m −2 rates of adsorption in different cities [18,108,109], which resulted in 10-200 g of particulate matter adsorption daily per tree. Thus, our results from the high-traffic Moscow center look comparable.
trees were in the middle with 25.5 g per day per tree (Appendix E). The highest adsorption peaks (in terms of maximum deposition velocity) were up to 1200 g (Larix sibirica) during several days with high concentrations in the air. This parameter showed drastic variation across time. The cumulative particulate mass adsorption over the season was 7.7 ± 1.5 kg for Acer platanoides (± standard error), 7.7 ± 2.1 kg for Larix sibirica, 3.8 ± 0.2 kg for Tilia cordata, and 1.4 ± 0.1 kg for Betula pendula (Figure 9, Appendix A). Saebø et al. [68] showed that Betula pendula improves air quality by particulate adsorption much better than Acer platanoides and Tilia cordata, however tree size should be taken into account as a factor to interpret the results of our study. It is well known that healthy large trees remove about 60 times more pollution annually than healthy small trees [67]. Additionally, leaf wax or leaf hair density and topography also have great influence [68,107]. In our case we could not include leaf morphology or topography since our estimation was based on LAI. However the total adsorption by tree seems comparable with several works yielding an average of 0.5-5 g m −2 rates of adsorption in different cities [18,108,109], which resulted in 10-200g of particulate matter adsorption daily per tree. Thus, our results from the high-traffic Moscow center look comparable.

Conclusions
According to our findings we can summarize ecosystem-services estimation for an individual tree as averaging 8.61 ± 1.25 kg (± standard error) of carbon stored, 137 ± 49 mm of water transpired, 2167 ± 181 kWh spent for microclimate regulation, and 5309 ± 808 g of PM10 adsorbed per the investigated period (July-November 2019). These numbers could be easily transformed to monetary

Conclusions
According to our findings we can summarize ecosystem-services estimation for an individual tree as averaging 8.61 ± 1.25 kg (± standard error) of carbon stored, 137 ± 49 mm of water transpired, 2167 ± 181 kWh spent for microclimate regulation, and 5309 ± 808 g of PM 10 adsorbed per the investigated period (July-November 2019). These numbers could be easily transformed to monetary values with the use of local prices for each of these services [110]. For Moscow, it would amount to about $150 per tree during the study period, mostly due to energy and particulate adsorption.
There are several approaches to provide ES information for the green infrastructure in urban areas. However, most of the inventory approaches, even when based on high resolution imaging, are limited by the temporal resolution which sometimes is important for detecting an early onset of ES decline. Our results show that an IoT tree-level network, using individual tree physiology sensing devices, such as TreeTalker+, or other similar devices, can be used in principle for monitoring urban green infrastructure ES in real-time. Furthermore, for some of the ES indicators, (e.g., water and cooling effects), they are most often based on models with indirectly-derived parameters [111]. Having real-time and individual tree data can improve our predictions and urban green infrastructure planning. There are several advantages for increasing the granularity of ES monitoring, particularly in that individual trees can be managed with a greater accuracy. The cost of monitoring is therefore critical for IoT expansion in green infrastructure monitoring. In recent years technological development and low cost microprocessors, traditionally used in automation and industry processes (Industry 4.0), are creating new opportunities for their expansion in environmental monitoring, which we could define as a Nature 4.0 transformation [36]. Typically, the total cost of a cluster of 20 trees is about 6000 euros, including the gateway, which corresponds to about 300 euros per tree. A single gateway can host up to 48 trees and the cost could decline with increasing monitored trees.
However, there are limitations and improvements to be considered in future work. First of all, the power consumption of the TreeTalker+ devices, used in the current work, is still a big limitation. Batteries need to be replaced every 1-1.5 months which requires a considerable human resource input. New batteries are being developed with much larger capacity that in principle could extend the battery life duration. In terms of improvement, a new IR sensor for distance sensing of canopy temperature could be very useful for improving the energy balance estimation and cooling effects. In particular, the installation of an anemometer will provide additional data on wind speed in the canopy, which influences the delivery of several ES. In addition, simple PM 2.5-10 optical devices can be included in the processor platform to get useful data on air quality using trees as monitoring stations. Further studies need to be conducted, but a noise sensor and microphone could also be included with the aim to provide useful information on the noise pollution and "soundscape" quality generated by trees in parks [35,112] and also to evaluate associated biodiversity with the help of recorded bird songs [34,113]. Nevertheless, the technical development of sensors along with people engagement to citizen science will be inevitable [32,114,115], thus it will be important to adapt them to the task of monitoring those parameters that are important for urban planning decisions [116].
Among the indicators presented in the article, perhaps not all of them can be directly used for practical purposes. For example, wind speed as well as air temperature and humidity under the canopy of city trees can be presented "as is" for the general public. However, in order to provide clear information on the quality of the urban environment associated with green infrastructure for public consumption, specific scales of air quality, microclimate comfort, and noise pollution levels should be developed. On the other hand, for spatial planning tasks, it will be useful to create an urban tree database on annual or seasonal indicators of ecosystem services provided by tree species at their specific age, height, and condition. This could be very useful for operational management of urban green infrastructures. [117]. In addition, it is also necessary to take into account disservices associated with urban trees such weakened and diseased trees falling on cars, infrastructure and buildings, and the allergic reaction of people to tree pollen [118,119]. These parameters should also be continuously monitored and reported in real-time for rapid response or timely prevention. Funding: Monitoring and modeling ecosystem services of urban trees was supported by Russian Scientific Foundation Project # 19-77-300-12. Data processing and analysis was supported by "RUDN University program 5-100". Appendix A Table A1. Tree description and summary of ecosystem services produced by each tree per investigated period (July-November 2019).