Spatio-Temporal Changes in Structure for a Mediterranean Urban Forest : Santiago , Chile 2002 to 2014

There is little information on how urban forest ecosystems in South America and Mediterranean climates change across both space and time. This study statistically and spatially analyzed the spatio-temporal dynamics of Santiago, Chile’s urban forest using tree and plot-level data from permanent plots from 2002 to 2014. We found mortality, ingrowth, and tree cover remained stable over the analysis period and similar patterns were observed for basal area (BA) and biomass. However, tree cover increased, and was greater in the highest socioeconomic stratum neighborhoods while it dropped in the medium and low strata. Growth rates for the five most common tree species averaged from 0.12 to 0.36 cm ̈ year ́1. Spatially, tree biomass and BA were greater in the affluent, northeastern sections of the city and in southwest peri-urban areas. Conversely, less affluent central, northwest, and southern areas showed temporal losses in BA and biomass. Overall, we found that Santiago’s urban forest follows similar patterns as in other parts of the world; affluent areas tend to have more and better managed urban forests than poorer areas, and changes are primarily influenced by social and ecological drivers. Nonetheless, care is warranted when comparing urban forest structural metrics measured with similar sampling-monitoring approaches across ecologically disparate regions and biomes.


Introduction
Urban forests are characterized by unique soils and urban morphologies, heterogeneous vegetation structure and composition, and novel assemblages of native and exotic tree species [1][2][3].The spatial and temporal characteristics of urban forest structure are driven by biophysical factors such as topography, climate, biogeochemical cycles, and disturbances such as drought [2,4,5].They are also affected by socioeconomic factors such as management and planning regimes, people's preferences, and socio-political budgets and directives [6][7][8].As all these drivers alter the structure and composition of an urban forest, so too are ecological processes affected and subsequently the provision of ecosystem services [6].A few studies have examined urban forests, their functions, and mortality with some examples from Mediterranean climate urban forests in California (United States, US), Spain, and Italy [9][10][11][12].However, other than these few studies, there is little information on how and why urban forest structure in South America, the global south, and Mediterranean climates changes across space and time.
Urban forests are a key component in cities and provider of ecosystem services as they influence the well-being of urban inhabitants [11,13].They are able to store and sequester carbon, regulate hydrologic cycles, ameliorate climate, remove air pollutants, provide habitat for fauna and space for recreation and spiritual enjoyment, among other services [13,14].However, multi-scale economic and socio-ecological drivers can alter urban forest structure, thus understanding the spatio-temporal dynamics of urban forest structure and composition can provide insights on their management and planning [15].
Tree planting preferences and management decisions by both private individuals and communities can influence urban forest cover and density, growth, mortality, and distribution [6].Similarly, species composition, age diversity, condition, site characteristics, and socio-ecological disturbances have an effect on growth of the urban forest [2,4,7,8], affecting several ecosystem processes, disservices and services.The type of land use and building density (i.e., urban morphology) and socioeconomic status in particular, affect the structure and distribution of the urban forest, leading to changes in biomass with consequences for carbon storage and sequestration [5,6,10,12].Growth, mortality, and regeneration in both urban and natural forests depend on environmental factors such as rainfall, temperature, and soil conditions; however, urban tree growth and biomass have been shown to be often greater in urban than natural landscapes [16].Other anthropogenic factors affecting growth and mortality include exposure to pollution, artificial irrigation, and vegetation management and maintenance practices [7,8,16,17].
Irrigation becomes extremely influential in Mediterranean semi-arid and arid cities where frequent droughts and water scarcity during growing periods can cause reduced growth, basal area, and biomass, and can increase mortality rates [15].These structural effects can be measured in growth rates, biomass, basal area, and mortality rates, and are key for analyzing temporal patterns of carbon sequestration and storage [17], tree wood waste biomass [5], and overall effects of climate on urban forest structure [6].Research on urban tree growth, biomass, and basal area is scarce and, to our knowledge, limited to a number of species mainly from North America [10,16,18] and a few studies from other regions [19][20][21].Further information on drivers of urban tree growth and mortality can improve decision-making on tree selection, urban forest carbon accounting, and benefit estimates [10,17,22,23].
Temporal studies of urban forest ecosystems are not common and are mainly based on coarse resolution information derived from satellite imagery [7,[24][25][26][27].Most research for monitoring the urban forest has been based on land use/land cover change analyses, with several examples from the US [5,7,16], Europe [24], and Asia [25], with few examples from Latin America [26] and Africa [27].Research based on inventories and field data are particularly scarce, aside from monitoring of planted street trees [10,20] and citation therein], and a few examples from the subtropics [5,8,16,28] and humid temperate areas of North America [7] exist.However, most of these mentioned studies are located in temperate and industrialized regions, with a scarce few studies from elsewhere [2,15] such as those in South America and from cities in Mediterranean climates.The use of permanent urban forest plots for monitoring the structure and composition of the urban forest is increasingly being used for not only recognizing necessary changes in planning and management goals, but also for distinguishing the most effective practices that maximize the provision of ecosystem services [4,7].Nevertheless, such research and information is rare for urban forests in South America and in Mediterranean climates, such as south and western Australia [15], southern Europe [20,23], and southern Africa [27].
The aim of this study is to analyze the spatio-temporal dynamics of a South American urban forest with a Mediterranean climate using tree and plot-level data from permanent monitoring plots.
Additionally, to our knowledge, such an approach would make this study one of the first of its kind in South America.Our specific study objectives are to analyze temporal changes from 2002 to 2014 and spatial differences in urban forest structural characteristics from Santiago, Chile including: (1) overall tree population growth, mortality, and ingrowth; (2) basal area, biomass, and tree cover change dynamics; (3) correlates of structural change; and (4) spatial patterns across the study area.We propose to address these research questions using statistical and spatial analyses of field data collected during 2002 and 2014.This type of study can be used to better understand the effects of urbanization and land use change on the structure of not only Latin American but Mediterranean climate urban forests as well.Findings can also be compared to similar plot based studies from other urban forests.Similarly, structural information, such as growth rates and biomass change, is key towards assessing ecosystem services and disservices that are frequently being used in urban planning, sustainability and climate change initiatives, and land management decisions.

Study Area and Field Sampling
Santiago, Chile is located in the middle of the Chilean Mediterranean climate zone (33 ˝27 1 S-70 ˝41 1 W).According to the Köppen climate classification systems, Santiago has a cool, semi-arid climate with warm, dry, hot, summers (November to March; [29]).Temperatures vary from an average of 20 ˝C in January to 8 ˝C in June-July with an annual average of 14.4 ˝C, while mean rainfall is 312.5 mm per year ( [29]).The city lies in the center of a valley surrounded by a coastal mountain range to the west and the foothills of the Andes to the east.The elevation varies from 400 to 900 m with an average of 540 m.In 2002 its population was 5.3 million inhabitants with a population density of 10,000 people/km 2 and 55,700 ha of built up area, which increased to 61,679 ha by 2009 [30].The study area of 967 km 2 is within the Santiago Metropolitan area and its 2014 population of approximately 7 million inhabitants, and encompasses multiple land use, land covers, and tenures.The eastern higher elevation portion of the study area is located in Andean piedmont shrublands while the western portion, once an Acacia spp.and grass dominated alluvial plain, has now mostly been altered to agricultural and urban land covers [12,30].
A total of 200 stratified random 400 m 2 plots, originally measured during January and February 2002 using criteria outlined in [12], were allocated across all of the Santiago Metropolitan area's 36 comunas and an additional four unincorporated ones (i.e., Colina, Lampa, Puente Alto, and San Bernardo).A comuna in Chile is a geographically and administratively delineated municipality.Of the 200 originally measured plots in the study area, 192 were relocated and re-measured during October 2014 to January 2015 (Figure 1).
Re-measured plots were assigned land use classes based on existing land use and land cover types originally defined in 2002 ( [12]; Table 1).We also classified the 2002 pervious surface cover types on each plot into percent: maintained grass, herbaceous vegetation, and bare soil.Percent impervious surface covers were also measured and include: cement, buildings, and paved surfaces.Plots were also post-stratified into low, medium, or high socioeconomic strata based on the Asociación Chilena de Empresas de Investigación de Mercado and an approach described in [12]; socioeconomic strata being defined as a combination of average annual income, education, vehicle ownership, and house services (e.g., fixed telephone).Due to the low number of individual trees per species found on the plots, we grouped tree species into five tree growth form types based on Lawrence et al.'s [16] classes: conifers, broadleaf-evergreen, broadleaf-deciduous, shrubs, and palms.

Plot Matching and Data Analysis
During 2014, plot locations were relocated using 2002 measured field reference data (e.g., Geographical Information System plot centroid coordinates, aerial photographs, plot center photos, and plot sketches).Once true plot center was identified using azimuth and distance to permanent reference object measurements (e.g., utility poles, manhole covers, building corners), we identified trees, palms, and shrubs originally measured in 2002.A tree, palm, or shrub is defined as a woody plant or palm with a diameter at breast height (DBH) of 2.5 cm or greater at 1.5 m above the ground surface.Plot data in 2002 recorded each tree's distance and direction to the plot center.This same information was available in 2014 to locate and identify individuals that were measured in 2002.If species matched and DBH was similar to the 2002 measurement, or wherever diameter was indicated to have been originally measured in the case of forked or deformed stems; trees were re-measured.Missing trees that were originally measured in 2002, but were not located in 2014, were further investigated to determine if stumps were present to confirm mortality or removal.Trees in 2014 that were not present in 2002 were also recorded.See [16,28] for details on specific plot and tree matching approaches.
Urban forest structure changes were analyzed using tree mortality and ingrowth variables as previously done in studies from North America [16,28].Specifically, we defined mortality as 2002

Plot Matching and Data Analysis
During 2014, plot locations were relocated using 2002 measured field reference data (e.g., Geographical Information System plot centroid coordinates, aerial photographs, plot center photos, and plot sketches).Once true plot center was identified using azimuth and distance to permanent reference object measurements (e.g., utility poles, manhole covers, building corners), we identified trees, palms, and shrubs originally measured in 2002.A tree, palm, or shrub is defined as a woody plant or palm with a diameter at breast height (DBH) of 2.5 cm or greater at 1.5 m above the ground surface.Plot data in 2002 recorded each tree's distance and direction to the plot center.This same information was available in 2014 to locate and identify individuals that were measured in 2002.If species matched and DBH was similar to the 2002 measurement, or wherever diameter was indicated to have been originally measured in the case of forked or deformed stems; trees were re-measured.Missing trees that were originally measured in 2002, but were not located in 2014, were further investigated to determine if stumps were present to confirm mortality or removal.Trees in 2014 that were not present in 2002 were also recorded.See [16,28] for details on specific plot and tree matching approaches.
Urban forest structure changes were analyzed using tree mortality and ingrowth variables as previously done in studies from North America [16,28].Specifically, we defined mortality as 2002 trees that could not be matched to 2014 trees.As such, mortality does not distinguish between tree removal due to maintenance or land clearing activities.Ingrowth was defined as the presence of a new tree in 2014, not measured in 2002 that could be the result of planting (i.e., existing trees that grew into the 2.5 DBH criteria), or natural regeneration.Average annual tree diameter growth (∆DBH; cm¨year ´1) was estimated for matched trees by calculating the change in DBH divided by the total number of days between measurements and then annualized using 365 days per year [16].Plot-level basal area change (∆BA; m 2 ¨year ´1) was also estimated using average annualized tree diameter growth for all trees on a plot [28].
We estimated individual tree biomass (kg) using aboveground allometric biomass equations from [19] and from the GlobAllomeTree database [31].If species-specific equations were not available, we used equations from the same genus, family, or tree type (i.e., conifer and hardwoods) following [32] and [33]'s approaches.Plot-level biomass was obtained as the sum of individual tree biomass on a specific plot; temporal changes in biomass (∆Biomass kg¨year ´1) were estimated by subtracting 2014 biomass from 2002 and dividing by the time since measurements.We did not include new biomass from 2002 due to ingrowth or biomass from 2014 lost to mortality.Plot-level basal area and biomass estimates were then converted to a per hectare basis for subsequent analyses.

Statistical Analysis
We tested for differences in DBH, BA, and biomass change at the species-level for both the five most frequent tree species as well as tree classification groups based on their form.We also analyzed differences in 2014 basal area and biomass and 2002 to 2014 ∆DBH, ∆BA, and ∆Biomass by comuna, socioeconomic strata, soil surface cover types, and land use classification.The 2014 biomass data was transformed using a logarithmic function for subsequent statistical analysis.We also analyzed for statistical differences in tree cover change as well as mortality and ingrowth according to socioeconomic strata, soil surface cover, and land use.
To test for statistical differences in mortality and ingrowth according to tree form, we used a five-sample test for equality of proportions without continuity correction using the R procedure prop.test[34] and alternative two-sided hypotheses.To better identify the plot-level surface cover factors, or correlates, that can possibly be driving the changes in structure, we determined the correlation between surface cover types and DBH, BA, and biomass (B).Specifically, we calculated Spearman correlation coefficients between 2002 to 2014 change in basal area and biomass with 2014 surface cover data using the cor procedure.Residuals were evaluated using Q-Q plots to assess their distribution and statistical differences were determined using Analyses of Variance with the aov function.This function was also used to test for differences in plot-level tree cover among strata while a paired t-test was used to determine city-wide tree cover differences during the analysis period.
Spatial correlations were determined using the Spatial Dependence: Weighting Schemes, Statistics and Model (spdep), and Data Analysis in Ecology (pgirmess) packages with the R correlog function.Finally, to better visually assess spatial patterns in BA and biomass we mapped urban tree biomass using spatial interpolation with an Inverse Distance Weighting (IDW) method using QGIS 2.10 and a distance coefficient of 0.05.Spatial autocorrelation was determined using Moran's I with the correlog procedure in the pgirmess package.All analyses were done using R version 3.1.3[34].

Results
Of the original 200 plots from 2002, only eight could not be re-measured due to lack of access and permission.Given Santiago's semi-arid environment, high building density, and low tree cover; our overall sample size and number of matched trees (Table 2) was low compared to other urban forest ecosystem studies from the subtropics [16,28,33] and more humid temperate areas [7].We found that 2014 tree biomass (kg¨ha ´1) was greater in the northeastern sections of the study area (Figure 2) and significant differences were found among socioeconomic strata (p < 0.01; Table 3).In general, wealthier comunas (i.e., the northeastern section of the study area) had greater tree biomass than lower income ones.We observed that city-wide plot-level tree cover was not significantly different (p = 0.45) during the analysis period, with a value of 18.2% (+/´1.6%) in 2002 that dropped to 16.6% (+/´1.6%) in 2014.Tree cover in the high socioeconomic stratum was significantly greater than the other two strata and increased during the period of analysis from 17.3% to 23.1%.Conversely, tree cover dropped from 16.8% to 13.1% in the medium strata, and from 20.1% to 12.8% in the low strata (Table 3).We also observed a high variability in ∆BA and ∆Biomass, but no significant differences were found across socioeconomic strata (p = 0.4, Table 3) or land uses (p = 0.71, Table 4).Overall, we found no statistically significant differences in both ∆BA and ∆Biomass, and we observed almost a static biomass accumulation and BA increases during the analyzed period (Table 4).Overall, biomass (p = 0.42) and BA (p = 0.37) were not significantly different across land uses (Table 4).However, commercial/industrial areas did present lower average biomass and BA values than other land use types.

Spatial Differences in Tree Biomass, Cover, and Basal Area
We found that 2014 tree biomass (kg•ha −1 ) was greater in the northeastern sections of the study area (Figure 2) and significant differences were found among socioeconomic strata (p < 0.01; Table 3).In general, wealthier comunas (i.e., the northeastern section of the study area) had greater tree biomass than lower income ones.We observed that city-wide plot-level tree cover was not significantly different (p = 0.45) during the analysis period, with a value of 18.2% (+/− 1.6%) in 2002 that dropped to 16.6% (+/− 1.6%) in 2014.Tree cover in the high socioeconomic stratum was significantly greater than the other two strata and increased during the period of analysis from 17.3% to 23.1%.Conversely, tree cover dropped from 16.8% to 13.1% in the medium strata, and from 20.1% to 12.8% in the low strata (Table 3).We also observed a high variability in ΔBA and ΔBiomass, but no significant differences were found across socioeconomic strata (p = 0.4, Table 3) or land uses (p = 0.71, Table 4).Overall, we found no statistically significant differences in both ΔBA and ΔBiomass, and we observed almost a static biomass accumulation and BA increases during the analyzed period (Table 4).Overall, biomass (p = 0.42) and BA (p = 0.37) were not significantly different across land uses (Table 4).However, commercial/industrial areas did present lower average biomass and BA values than other land use types.

Changes in DBH, BA, and Biomass According to Species and Tree Form
The five most frequent tree species represented approximately 25% of all our matched trees, and their overall growth rate by species ranged from 0.12 to 0.36 cm¨year ´1 (Table 5).These species also had an increase in biomass between 0.46 to 6.47 kg¨year ´1 during the analyzed period; with R. pseudoacacia having the highest growth rates and P. ceracifera the greatest biomass change (Table 5).Actual ∆DBH (p < 0.0001) and ∆Biomass (p = 0.03) during the analysis period were significantly different among the five most frequent tree species (Table 5).Growth diameter rates (∆DBH) were greater for conifers, while conifers and broadleaf deciduous trees exhibited the greatest change in biomass (Table 6).The large standard errors found for conifers and palms are likely due to a very small sample size, stem shrinkage or swelling, and the palm's stem sheath and fronds that makes DBH re-measurements unreliable [4] (Table 6).

Surface Cover Correlates
We assessed the use of plot-level surface cover as an indicator for urban forest structural changes and found that both BA and biomass were poorly correlated with all measured surface covers (Table 7).Increases in cover indicate lower BA and biomass, however, low building cover was not indicative of BA and Biomass.We only found that cement coverage was correlated to 2002 Basal Area (r = 0.43) and 2002 and 2014 Biomass (r = 0.41 and r = 0.36, respectively; Table 7).

Spatio-Temporal Changes and Patterns in Biomass and Basal Area
Similar to the previous statistical analyses, our spatial analyses also showed greater tree biomass towards the northeastern section of the city and a cluster of plots with greater biomass in the peri-urban, southwest comunas of the study area (Figures 2 and 3).These same plots had greater basal area.Spatial interpolation by IDW of ∆Biomass also exhibited a similar pattern to the plot-level biomass distribution (Figures 2 and 4).Plots and areas in the central, northwest, and southern sections of Santiago showed losses in both ∆BA and ∆Biomass (Figure 3). Figure 4 maps biomass change but should be interpreted carefully since the IDW spatial interpolation based on the plot influence distance for urban areas has high uncertainty, given the non-continuous heterogeneous character of cities.Similarly, we found some significantly high spatial autocorrelation according to Moran's I statistic in plots that were less than 4 km in proximity relative to each other.

Discussion
Although no clear evidence can be observed in the 2002 and 2014 Landsat imagery (Figure 1), the Santiago Metropolitan Area has undergone noticeable urbanization in the form of infill and land use changes since 2002 [30].Rapid economic development and recent large-scale infrastructure projects such as new transportation rights of way and hubs, tunnels, and building-housing projects [35] has negatively affected the structure of Santiago's urban forest.These changes have occurred for most of the city with the exception of the most affluent areas, the northeast and few areas in the oldest sections of the city center [35].Santiago's semi-arid climate and relative drought related water scarcity explains its low tree cover; however, even small changes can have noticeable effects in terms of the urban forest structure and subsequent functions and services.The tree scarcity in semi-arid Santiago resulted in different patterns than previous studies that used similar field sampling approaches from wet and humid urban forests in temperate and subtropical North America (e.g., Syracuse, USA [7]; Gainesville, USA [16]; Orlando, USA [33] and San Juan, Puerto Rico [28]); making this study unique.
These socio-ecological and urbanization effects were observed in the results from both our statistical and spatial analyses.This is particularly evident from the differences in tree cover, losses in biomass throughout the study area, and BA changes among socioeconomic strata (Table 3).More affluent comunas are mostly located in the northeastern section of the city and, because of greater resources, have more tree cover, larger trees (i.e., greater BA), more abundant trees, and therefore, greater ∆Biomass (0.29 Mg¨ha ´1¨year ´1) [6,15,36].Conversely, plots on the medium and lower socioeconomic strata had losses in overall tree cover and negative ∆Biomass.These plots generally correspond to areas with higher building densities and comunas with lesser economic resources ( [12]; Figure 4).Spatially, plots in the northeastern section of the city had greater tree biomass; however, in the far southwest section of study area we also found high biomass areas corresponding to agricultural land from the peri-urban section (Figure 2).Rainfall gradients across the study areas could also be playing a role in structural differences as upper elevation plots in the Andean piedmont receive more rain, but management and maintenance regimes (i.e., irrigation and fertilization) could be masking the effects of precipitation on tree growth and mortality.
Poor correlations between surface cover and BA and biomass might be a result of having small sample sizes and a number of plots with both high cement cover and high biomass and others with high grass cover and low biomass.High cement cover and high biomass are typical of large street trees near building areas.Using the sum of the different vegetation cover (horizontal) and height (vertical) estimates might increase the correlations, as observed by [33].We noted that cover types are regularly used as indicators of urbanization, planting space, and urban forest structure and function [7,16], but given our semi-arid and high building densities, this might have resulted in different patterns.
We acknowledge that our sample size, in terms of plot density and tree numbers, was low relative to similar urban forest studies using fixed area plots [7,16,28,33].Factors-not analyzed in this study-such as climate change and socio-ecological dynamics-could have also driven these changes.Additionally, our use of forest grown allometric equations on urban trees can often lead to errors of up to 40% in biomass estimates [5].Similarly, the sample size for conifers and palms was very low, seasonal stem shrinkage and swelling is likely [4,16], and this was also confounded by the difficulties in measuring palms using measurement techniques developed for single stem temperate trees [4,16,28].However, we note that all these cited studies had the same limitations.Indeed, given the context of our study area's: size, socio-political dynamics, and limited infrastructure that characterizes cities from developing countries, we feel that our results do provide a better understanding of urban forests outside the frequently studied areas of the United States and Canada, giving more insights to understudied areas such as Latin America [1,13,32], Mediterranean climates [12,28,32], and Australia [3,15].
When comparing urban forest structure results such as annual mortality to other studies using inventory data, we can distinguish the effect of climate and management.Mortality of urban trees was apparently more related to management than to climate; Santiago showed much lower mortality rates (3.0%) than subtropical San Juan, Puerto Rico [28], with values between 30%-40%.Conversely a subtropical city in the United States, Gainesville [16] had an average mortality rate between 10%-19%.Looking at growth rates, climate appears to be a stronger driver than management; values for subtropical cities of San Juan and Gainesville are closer to each other (0.7-1.0 cm¨year ´1), while Santiago had lower average values between 0.1 to 0.4 cm¨year ´1.In terms of the temporal changes to urban forest structure, there are no similarities to relevant studies [7,28].
Studies that use permanent random plots often report number of trees per unit area-as opposed to basal area-as a measure of tree density [7]; however, multi-stemmed trees with a shrub form or secondary tropical forests and mangroves can confound comparisons, thus our preference for basal area.Values for tree density can vary from 222-328 trees per ha in San Juan in 2001 and 2010, respectively [28].While there is an average of 34 trees per ha in a temperate city such as Syracuse, US [7], in Santiago, a Mediterranean shrubland biome, this value reached 64 trees per ha [12].Meanwhile, there was an overall annual net loss of approximately four trees per ha in subtropical Gainesville [16].Tucker-Lima et al. report a basal area of 4.6 m 2 ¨ha ´1 in San Juan, which is very similar to our 2014 estimate of 4.8 m 2 ¨ha ´1.Changes in basal area in Santiago of 0.1 m 2 ¨ha ´1 were, however, much lower than in San Juan (1.0 m 2 ¨ha ´1) [28].
The growing body of literature using similar sized, long-term monitoring plots as utilized in this study can be used to compare trends for urban forests across different regions of the globe [5,6,16,28,33,36].The characteristics of urban forest growth, mortality, and the effects of site characteristic such as irrigation, and ecological disturbance on these, could be analyzed against field measured plot and site correlates [5,16,20,21].However, the presence of palms, multi-stemmed tall shrubs, and size criteria for tree-shrub differentiation across different biomes can confound some of these comparisons.Care is also warranted when comparing shrubland dominated biomes such as Santiago, to dense subtropical secondary and mangrove forests or temperate forests.Similarly, because of the difficulty in sampling large, heterogeneous, urbanized areas, sampling intensities can be low and result in larger uncertainties.This is particularly true in Santiago's densely built, semi-arid, urban context, which also resulted in a reduced sample of measured and matched trees.
Our study provides one of the few comparative insights into how a South American and Mediterranean urban forest changes across space and time.Future research could analyze spatio-temporal changes in urban forest composition and its subsequent effects such as the spatial dynamics of ecosystem service provision and disservice hotspots [11,13,21,22,[24][25][26]. Plot level structure-function information could also be used to test land management and planning scenario effects on the demand and supply of services [4,7,14,33].Quantitative analyses could also determine the socio-ecological causal factors behind, and the drivers related to, changes in urban forests, their processes, and services [1,3,6,23,36].More basic research could use these monitoring sites for better understanding floral diversity and the occurrence of invasives and understory dynamics in urbanized forests, as well as for developing allometric equations [36][37][38][39].Permanent plots could also facilitate dendrochronological analysis of different tree species and their growth and mortality as affected by climate change, pollution, maintenance practices, and community preferences.

Conclusions
Overall, we found that urban forest mortality, ingrowth, and tree cover in the greater Santiago area remained stable during 2002 to 2014.We also noticed slight losses in basal area and biomass change across the study area.However, there were some noticeable trends during the analysis period in that tree cover increased and was greater in the highest socioeconomic stratum; however, it decreased in the medium and low strata.Similarly, the less affluent central, northwest, and southern plots, in particular, exhibited losses in BA and biomass from 2002 to 2014.
As previously mentioned, other studies have used plot-level data to analyze changes in an urban forest but, to our knowledge, most of these studies are from North America.Here, we describe for the first time the spatio-temporal dynamics of a South American urban forest in a Mediterranean climate.Overall, we found that these urban forests follow similar trends as in other parts of the world.For example, affluent areas tend to have more and better condition urban forests than poorer areas.Thus, barring substantial ecological disturbance (e.g., storms, drought, pests, or urbanization), changes in urban forest structure will primarily depend on human management (e.g., maintenance and irrigation) and biophysical factors (i.e., growing space and climate).Findings could also be used to identify which tree species perform better in terms of growth and basal area.Accordingly, this information could be used to identify tree functional traits that are most associated with ecosystem service provision (e.g., carbon offsets) and disservice minimization (e.g., allergenic tree locations).
However, Latin America and the rest of the developing world are generally characterized by unplanned land uses and marked socioeconomic inequities.Thus, studies such as ours can be used for targeting specific sites for improved management and setting monitoring and evaluation standards for municipalities.Conversely, areas that are maintaining their urban forest and standards-despite limited budgets-could also be identified.Most importantly, issues of environmental justice and resilience to the effects of climate and socio-political changes could be better addressed with this long-term data.Indeed, the growing body of literature calls for the development of an available global urban forest monitoring database that could be used for research by the global scientific community.

Figure 1 .
Figure 1.Plot locations within the Santiago de Chile study area in 2002 (left) and 2014 (right).Plots located on high, medium, and low socioeconomic strata are represented in red, gray, and blue circles, respectively.Vegetated land covers are green in color, bare soil areas are in pink, while orange and brown depict sparsely vegetated areas.High to low urban density areas are dark to light magenta, respectively, and dark magenta, linear features are major transportation rights of way.Linear light magenta areas in the southeast and far east in the Andean foothills are a river flood plain and barren rock, respectively.Note: images are from Landsat-7 RGB 742 taken on February 2002 and January 2014.

Figure 1 .
Figure 1.Plot locations within the Santiago de Chile study area in 2002 (left) and 2014 (right).Plots located on high, medium, and low socioeconomic strata are represented in red, gray, and blue circles, respectively.Vegetated land covers are green in color, bare soil areas are in pink, while orange and brown depict sparsely vegetated areas.High to low urban density areas are dark to light magenta, respectively, and dark magenta, linear features are major transportation rights of way.Linear light magenta areas in the southeast and far east in the Andean foothills are a river flood plain and barren rock, respectively.Note: images are from Landsat-7 RGB 742 taken on February 2002 and January 2014.

Figure 2 .
Figure 2. Aboveground tree biomass distribution in Santiago, Chile in 2014.Larger sized circles represent plots with higher biomass estimates.

Figure 2 .
Figure 2. Aboveground tree biomass distribution in Santiago, Chile in 2014.Larger sized circles represent plots with higher biomass estimates.

Figure 3 .
Figure 3. Plot-level urban tree biomass change in Santiago, Chile from 2002 to 2014 based on 192 remeasured plots.

Figure 4 .
Figure 4. Urban tree biomass change prediction map for Santiago, Chile using Inverse Distance Weighting interpolation using the re-measured plots.

Figure 3 .
Figure 3. Plot-level urban tree biomass change in Santiago, Chile from 2002 to 2014 based on 192 re-measured plots.

Figure 3 .
Figure 3. Plot-level urban tree biomass change in Santiago, Chile from 2002 to 2014 based on 192 remeasured plots.

Figure 4 .
Figure 4. Urban tree biomass change prediction map for Santiago, Chile using Inverse Distance Weighting interpolation using the 192 re-measured plots.

Figure 4 .
Figure 4. Urban tree biomass change prediction map for Santiago, Chile using Inverse Distance Weighting interpolation using the 192 re-measured plots.

Table 1 .
2002 land use/cover classes analyzed in Santiago, Chile.

Table 1 .
2002 land use/cover classes analyzed in Santiago, Chile.

Table 2
below provides an overview of the basal area, ingrowth, and mortality in Santiago, Chile's urban forest from 2002 to 2014.

Table 2 .
Basal area and average annual urban forest percent mortality and ingrowth for different tree forms in Santiago, Chile during 2002 to 2014.

Table 2 .
Basal area and average annual urban forest percent mortality and ingrowth for different tree forms in Santiago, Chile during 2002 to 2014.

Table 3 .
Mean aboveground tree biomass in 2014, change in aboveground tree biomass (∆Biomass), and mean percent tree cover change (∆Tree Cover%) according to socioeconomic strata during 2002 to 2014 in Santiago, Chile.

Table 4 .
Mean tree biomass, change in tree biomass (∆Biomass) during 2002 to 2014, mean basal area, and change in basal area (∆BA) during 2002 to 2014 in Santiago, Chile.
Note: SE is Standard error.

Table 5 .
Growth diameter rate (∆DBH), change in biomass (∆Biomass) and basal area for the five most frequent tree species in Santiago, Chile during 2002 to 2014.

Table 6 .
Annual average change in growth diameter rate (∆DBH) and aboveground biomass for different tree forms in Santiago Chile during 2002 to 2014.
* Standard errors were much too large to report value.