Influence of Climate and Management on Patterns of Taxonomic and Functional Diversity of Recreational Park Vegetation

Recreational urban parks support diverse assemblages of plants that through their functions, contribute beneficial services to billions of individuals throughout the world. Drivers of vegetation-derived services in parks are complex, as climate and park management interact with the functioning of multiple species of vegetation types. Yet, informal observations suggest that recreational parks are constructed consistently to specific principles of landscape design. Here we ask: what are the patterns of functional traits and vegetation diversity in cities of varying climate in the United States, and how do these patterns result in a consistent typology of recreational park? We hypothesized that increased aridity would exclude species not adapted to warm, dry climates, thereby reducing local, or alpha, taxonomic diversity and shifting community composition. However, a similar preference of park managers in the United States for suites of service-based functional traits leads to similarity of mean values of services traits in recreational parks among cities, regardless of climate differences. We tested this hypothesis by surveying lawn species, comprised of herbaceous turf and spontaneous plants, and woody species in fifteen recreational parks across Baltimore MD, Riverside CA, and Palm Springs CA, three cities that contain multiple parks but differ in regional climate. With increasing aridity, taxonomic alpha diversity decreased and plant physiology shifted, yet no differences were observed among most service-based functional traits. Among the cities surveyed, no significant differences were observed in functional dispersion of woody and spontaneous species or most service-based traits. Taxonomic composition differed in each city for all vegetation types, while suites of service-based traits differed between Baltimore and the two more arid cities of Riverside and Palm Springs. Our results suggest that across the United States, service-based functional traits are consistent, even when arising from unique compositions and abundances of species in recreational parks. We interpret these results as an interaction between climate and the preferences of recreation park managers for services, creating a pattern of vegetation diversity where taxonomic alpha and beta diversity vary among regions while specific suites of services remain available.


INTRODUCTION
Urban vegetation comprises novel ecosystems by bringing together native and non-native plant species of no natural analog (Hobbs et al., 2014;Aronson et al., 2017). In cities across climates and cultures, one finds novel vegetation arrangements within recreational parks (Threlfall et al., 2016;Weems, 2016;Talal and Santelmann, 2019). Recreational parks are a public form of green infrastructure that provides services to individuals through access to recreation, cooler temperatures, and improvement to overall health (Bolund and Hunhammar, 1999;Ayala-Azcárraga et al., 2019). These urban greenspaces contain multiple species of trees planted to provide esthetic benefits and shade as well as large expanses of green turf for recreation and gathering (Tinsley et al., 2002;Pataki et al., 2013;Talal and Santelmann, 2019). However, specific services provided to individuals and neighborhoods by planted vegetation may be contingent on the diversity of species and vegetation function (Larson K. L. et al., 2016;Vieira et al., 2018).

Taxonomic and Functional Diversity in Recreational Park Vegetation
Examining patterns of taxonomic and functional diversity among communities allows us to understand the drivers behind local species assembly and functioning (Anderson et al., 2011;Johnson et al., 2015). What drives urban diversity is a blend of biophysical and social factors (Padullés Cubino et al., 2018). Climate and management preferences can determine the expression of plant functional traits in parks, from the physiological of leaves to esthetic characteristics preferred by people (McCarthy et al., 2011;Pataki et al., 2013;Avolio et al., 2018). Vegetation physiological function has also been linked to the provisioning of ecosystem service-based traits Swan et al., 2016). The influence of climate and management preference can vary depending on the type of vegetation studied (Nielsen et al., 2014). A metaanalysis of species richness in urban parks found out of ten studies focused on vascular plants, seven studies restricted results to only woody species (Nielsen et al., 2014). Other recent studies in urban green systems focus only on trees (Avolio et al., 2015;Gillespie et al., 2017;Kendal et al., 2018), or herbaceous turf and weedy species (Knapp et al., 2012;Wheeler et al., 2017;Padullés Cubino et al., 2018). To understand the causes of vegetation diversity in an urban plant community, we need to know how regional (climate) and local (management) influence different vegetation types in a specific urban community (Aronson et al., 2016), such as recreational parks.
Building on the framework of environmental filtering of communities (Ackerly, 2003;Spasojevic and Suding, 2012;Kraft et al., 2015), extreme aridity reduces physiological niche space, thus reducing taxonomic diversity and shifting physiological function of park vegetation (Spasojevic et al., 2014;Kraft et al., 2015;Pearse et al., 2018). Cities of extreme climates generally have lower species richness of urban trees and lawn vegetation compared to milder climates Padullés Cubino et al., 2018). Within lawn communities, these patterns may arise due to differences in irrigation practices across cities, as well as differences in turf varietals (Kendal et al., 2018;Pincetl et al., 2019). Cities in more arid climates, such as Los Angeles and Phoenix, irrigate their lawns, while more mesic cities have both irrigated and non-irrigated lawns (Wheeler et al., 2017). Warmer climates also allow land managers to plant species of turf which use the C4 photosynthetic pathway, which can maintain productivity under high temperatures and aridity, while C3 turf species are more physiologically restricted, making them a poor choice for hot and dry cities environments (Wherley and Sinclair, 2009;Beard, 2013). Furthermore, lawn communities are often colonized by spontaneously growing herbaceous species that has a greater range of functional properties compared to turf, potentially allowing more spontaneous species to pass through environmental filters that limit turf species (Robinson and Lundholm, 2012).

Services Preferred by Recreational Park Management
While taxonomic composition and physiological functions are strongly driven by climate, similar preferences for suites of traits related to specific service-based functions among urban residents may decrease functional diversity of urban plant communities (Fukami et al., 2005;Larson K. L. et al., 2016). The service-based traits chosen for recreational parks reflect how park managers design park vegetation, where the park tree composition can provide services such as shade and visual esthetics not found in the local native communities . Lawns provide similar services across cities, such as vibrant greenness and resiliency to active use (Larson K. L. et al., 2016). However, the effect of preference in parks on service-based traits can differ depending on the plant community and type of trait. In warmer cities, there can be an increase in service-based traits relating to showy flowers and fruits, however across a gradient of aridity, differences in tree heights or leaf areas have not been identified Pearse et al., 2018).
Service-based functional traits can have distributions less limited by climate compared to ecophysiological traits. For example, esthetic flower production is a valued service trait not limited by climate due to seasonally available flowering species in most regions (Konijnendijk et al., 2005;Threlfall et al., 2016;Goodness, 2018). However, esthetic fall foliage is prevalent in broadleaf deciduous species, which are likely to be found mostly in biomes like temperate deciduous forests. Preferences allow for similarity of esthetic flower traits across climates while having fewer effects on traits of esthetic fall foliage. Alternatively, spontaneous vegetation is not actively planted, making its dispersal more limited than horticultural vegetation (De Wet and Harlan, 1975). After undergoing selection by dispersal and environmental filters, the selected traits of spontaneous species in recreational parks are those that can overcome mowing, weeding, and physical dispersal barriers. Overall, these at times independent variations of taxonomic, ecophysiological, and service-based diversity implies multiple drivers may determine the pattern of diversity in urban park vegetation.
While visually one can see that recreational parks appear similar in form to one another regarding their design, the question remains who controls that typology and the diversity found within? Urban park tree diversity has been found to respond to different physical and social drivers compared to residential yards and street trees (Kendal et al., 2012). What vegetation is found in recreational parks is ultimately a function of the park managers who decide what is planted. A recent study based in Portland OR found the reasons behind vegetation design in recreational parks are often a function of economics, maintenance of the park environment from human disturbance, and a desire to provide beneficial services for visitors (Talal and Santelmann, 2020). Similar management desires for economy of maintenance were found in surveys of park managers in Hong Kong (Chan et al., 2014). The differences between park visitor preferences and manager preferences can often be due to budget constraints (Baur et al., 2013), where park visitors may desire an increase in colorful foliage and more flowers, managers are looking for ways to improve maintenance on a restricted budget (Talal and Santelmann, 2020). While park managers may have similar goals in different countries, the species diversity found in parks is often determined by national and local horticultural trends Nielsen et al., 2014;Roman et al., 2018), as well as regional climate .
To better understand how climate and management preferences influence the distributions of recreation park diversity, we ask: how do park biodiversity, community composition, and values of service-based traits differ among three mid-sized United States cities of varying climates? We predict that alpha diversity will vary among recreational parks in arid and mesic cities due to the physiological tolerances of vegetation, and that each city will have a unique taxonomic beta diversity. Functional diversity and service provisioning, however, will not differ across climates owing to the preference of park managers for a specific typology of vegetation composition and provided services within recreational parks. To resolve the gap spanning multiple vegetation types within the same community, we focus on three different vegetation types, including horticultural woody species and lawn communities comprised of turfgrass and spontaneous herbaceous vegetation. By quantifying patterns in taxonomic and functional diversity and their relationship to services provisioning in recreational parks, we describe how broad-scale climatic variables interact with the preferences of recreational park managers to influence park vegetation diversity and function.

Study System
We assessed climate and park managers preference as drivers of park vegetation biodiversity by sampling 15 recreational parks located in three cities within the United States: Baltimore MD, Riverside CA, and Palm Springs CA (Figure 1 and Table 1). Cities were selected to represent a gradient of aridity (measured as the vapor pressure deficit, VPD). Baltimore (mean summer VPD: 0.89 kPa) is in a temperate humid climate characterized by the regular precipitation events throughout the growing season. Riverside (mean summer VPD: 2.7 kPa) is in a Mediterranean climate of intermediate aridity and Palm Springs (mean summer VPD: 6.2 kPa) is in a desert climate of extreme aridity, and both cities are characterized by seasonal precipitation during winter and early spring with very little precipitation during summer and fall. Within Riverside and Palm Springs, urban vegetation health is maintained through regular irrigation throughout the year. In each city, we selected five recreational parks. We focused only on highly managed recreational parks, rather than natural parks or parks with remnant natural areas, to highlight effects of park managers' preferences for recreation. Recreational parks are characterized by large lawn areas with multiple species of planted, maintained woody species. Lawns, along with spontaneous species found within the lawn, are the only herbaceous species present in the park and no area is unmanaged. To test whether species-area relationships influence diversity (Nielsen et al., 2014;van der Maarel and Franklin, 2015), our parks represented a range of areas, from 10,000 to 300,000 m 2 . Park size was determined via ArcGIS 10.6 (ESRI, Redlands, CA, United States) polygon delineation ( Table 1). Our field sampling protocol was designed to collect all taxonomic data for the five parks within each region over 7 days of field work in each city. To assess climate as a driver of physiological function, parks were sampled during mid-afternoon in peak summer months. Our experimental protocol is designed to be transferrable to other cities to characterize recreational park vegetation can within a approximately 5 days.

Functional Trait Analysis
Our research protocol is designed for maximum repeatability, even by researchers without access to physiology wet labs. To that end we collected functional data in silico, using trait databases and primary literature. To incorporate functional diversity metrics into our analysis of recreational park species, we used the BIEN database (Maitner et al., 2018), the Global Wood Density Database (Chave et al., 2009) and primary literature to assemble average ecophysiological (EP) traits for each sampled species ( Table 2). Following Westoby's leaf-height-seed plant ecology strategy scheme, we included traits that relate to a holistic FIGURE 1 | Study area of cities of interest. Five recreational parks were sampled in each city. All parks are representative of a recreational park typology, comprising expansive lawns and individually planted trees. Images are sourced from Google Earth Pro. ecological strategy: Specific Leaf Area, Seed Dry Mass, and Height (Westoby, 1998). To incorporate traits that relate to hydraulic architecture in stems, we included wood density as an additional trait (Swenson and Enquist, 2007;Chave et al., 2009) for woody species. To analyze influences of park manager preference, we included a series of service-based traits for woody species (ES). These traits are based on desired attributes found in nursery stock . We included traits reflecting shade and fruit provisioning, as well as flower, fruit, foliage, and bark esthetics. While there are many service-based traits found in horticultural vegetation, such as carbon capture (Pataki et al., 2006), pollution and dust mitigation (Wang et al., 2019), and microclimate regulation (Shiflett et al., 2017), we chose to focus on traits that were more associated with attributes preferred by urban residents and urban land managers Larson K. L. et al., 2016). We used horticultural sources to determine ES traits for each woody species ( Table 2). ES traits were analyzed as presence/absence, where if a horticulture record listed a trait as associated with the species, that trait was given a value of 1 and if no record of the traits was found, we recorded a value of 0. We only used species that had records for all trait factors for subsequent analyses. As there was an unequal number of species containing all ES and all EP traits, we separated woody species into two categories to capture all species with a full spectrum of service-based traits, and all species with a full spectrum of EP traits.
We analyzed herbaceous species as a lawn, encompassing all herbaceous species, as well as separated by functional groups of turf species or spontaneous species. As turf and spontaneous species were all herbaceous, we did not include wood density as an EP trait and focused only on those traits in the Leaf-Height-Seed strategy. ES traits were not calculated for turf or spontaneous species. Turf is selected based on its ability to persist; thus if turf is present a singular selectable service is being provided (Christians et al., 2016). Spontaneous species are not preferentially selected for planting. We filled gaps in trait data through searches of primary sources ( Table 2).

Statistical Analysis
We conducted all statistical analyses in R (R Core Team. 2019). For each vegetation type we created a site by species matrix, each site being one sampled recreational park. Taxonomic alpha diversity of vegetation types in parks was represented by the Shannon-Wiener index calculated in the vegan package in R (Oksanen et al., 2008). Statistical significance was determined using one-way ANOVA calculated for each cover type; differences among cities were determined with a Tukey's HSD post hoc test. Taxonomic beta diversity between any two parks was calculated using Brey-Curtis distance, and visualized through Principle Coordinates Analysis (PCoA). Within-city beta diversity was defined as the Euclidean distance to centroid, and significance was determined through a PERMANOVA, while differences in Park area was determined through ArcGIS, percent impermeable surface was determined from the 2016 National Land Cover Database (Homer et al., 2012), and adjacent land use was determined through county zoning maps of each city.
composition were determined using a pairwise PERMANOVA in the RVaideMemorie package in R (Maxime, 2019). To examine differences in beta diversity, we calculated dissimilarity matrices between all parks in the study, using both taxonomic and functional characteristics. Functional beta diversity was calculated using Gowers distance. Any trait distribution that did not fit a normal distribution was log-transformed and all trait data were z-transformed before any multivariate analysis was completed. Using PERMANOVA, we tested for compositional shifts in taxonomic and functional diversity for all vegetation types, and visualized results by the PCoA. The dissimilarity matrices allowed us to calculate homogeneity of variance through PERMANOVA and we identified compositional shifts in taxonomic and functional diversity. Functional alpha diversity was defined as the functional dispersion (FDis) of each park. FDis is an abundanceweighted measure of the distribution of functional traits in a community and quantifies the range of functional strategies within a community (Laliberte and Legendre, 2010). To quantify individual functional metrics, we calculated community weighted trait means (CWM) for each trait. CWM values provide a singletrait metric that is weighted by species abundances within a sampled community (Zuo et al., 2016). FDis and CWM analyses were completed using the FD package in R (Laliberté et al., 2015).

Taxonomic and Functional Alpha Diversity
Taxonomic and functional alpha diversity was first analyzed to determine if park area acted as a significant driver. Park area was not correlated with either taxonomic or functional alpha diversity. We then conducted subsequent analyses investigating city-wide climate as a primary driver of diversity metrics.
For both woody and turf vegetation, taxonomic alpha diversity, measured as Shannon-Weiner diversity, was lower in recreational parks of Palm Springs compared to the more mesic city of Baltimore (Woody species EP traits, p = 0.011; woody species ES traits, p = 0.021; turf species, p = 0.008). Differences in Shannon diversity of woody vegetation between Riverside and either Baltimore or Palm Springs were not detected. Differences in lawn taxonomic diversity were observed between the mesic Baltimore and the two drier cities of Riverside (p = 0.050) and Palm Springs (p = 0.010). However, when isolating turf  (Schmid and Brenzel, 2001) 5. (Maitner et al., 2018) 6. (Kröber et al., 2012) 7.  8. (Johnson and Gerhold, 2001) Trait values are determined from cited sources. and spontaneous species, the turf in Palm Springs was different compared to the other two cities (Palm Springs -Riverside p = 0.024, Palm Springs -Baltimore p = 0.008) and no differences were found when comparing spontaneous species assemblages (p = 0.115) (Figure 2A). Functional alpha diversity, FDis, of both EP and ES traits in woody species were not significantly different among the three cities (EP p = 0.487; ES p = 0.659) (Figure 2A). Likewise, FDis of all lawn species and only spontaneous species were not significantly different across cities (p = 0.247). Within turf species, there was a significant difference in city-scale FDis means (p = 0.034). No significant difference of FDis means was found between Baltimore (mean FDis = 0.206) and Riverside (mean FDis = 0.200), however Palm Springs (mean FDis = 0.047) values were significantly lower than Baltimore (p = 0.050) (Figure 2C).

Functional Traits
Woody EP traits all showed shifts in CWM across cities. The three traits associated with the Leaf-Height-Seed strategy all decreased FIGURE 2 | Metrics of alpha diversity for each vegetation type (mean and standard deviation). Asterix refer to significant differences between groups tested through ANOVA with a post hoc Tukey HSD test (*0.001 < p < 0.05). Herbaceous species are both spontaneous and turf species (A) Taxonomi.c alpha diversity was measured by Shannon-Weiner Diversity, (B) Functional alpha diversity (FDis) of woody species, (C) FDis of herbaceous type species. FDis was measured by Functional Dispersion, an abundance weighted metric of functional trait ranges.
from Baltimore to Riverside and Palm Springs (SLA: Baltimore -Riverside p < 0.001, Baltimore -Palm Springs p < 0.001; Height: Baltimore -Palm Springs p = 0.002; Seed Dry Mass: Baltimore -Riverside p < 0.001, Baltimore -Palm Springs p < 0.001). Wood density increased from Baltimore to Riverside and Palm Springs (Baltimore -Riverside p = 0.031, Baltimore -Palm Springs p = 0.001) (Figure 3A). Overall values of woody ES traits display little dissimilarity among cities. Shade and fruit provisioning, and FIGURE 3 | Community weighted means (CWM) of woody functional traits (mean and standard deviation). Asterix refer to significant differences between groups tested through ANOVA with a post hoc Tukey HSD test (*0.001 < p < 0.05, ***p < 0.0001). (A) Traits related to ecophysiological characteristics; (B) Traits related to service-based functional trait characteristics. esthetic traits of flowers, fruit and bark were not different across our study cities. In each city, the highest value for any servicebased trait was that of shade provisioning. The lone service functional trait which exhibited differences among cities was that of esthetic fall foliage; more trees associated with the fall foliage trait were found in Baltimore compared to Riverside and Palm Springs (Baltimore -Riverside p < 0.001, Baltimore -Palm Springs p < 0.001) (Figure 3B).
Within lawn species, only seed dry mass differed across cities, with Palm Springs having significantly lower values FIGURE 4 | CWM of ecophysiological traits of herbaceous, turf, and spontaneous species (mean and standard deviation). Asterix refer to significant differences between groups tested through ANOVA with a post hoc Tukey HSD test ( * p < 0.05, * * p < 0.001, * * * p < 0.0001). Observe that y-axis varies for each trait type.
than Baltimore (p = 0.031). Similar results were observed in spontaneous species, where seed dry mass was higher in Baltimore compared to Riverside and Palm Springs (Baltimore-Riverside, p = 0.016; Baltimore-Palm Springs, p = 0.012). Turf species in Riverside showed larger specific leaf area compared to Baltimore (p = 0.021), while vegetation height in all herbaceous vegetation groups was not different among cities (Figure 4).

Composition / Beta Diversity
All cities displayed a significant difference in taxonomic composition from each other except in lawn species composition between Riverside and Palm Springs (Figure 5A, p values displayed in Table 3). In a pairwise PERMANOVA, we found no difference between all lawn and only spontaneous functional composition among cities; woody ES and EP and turf composition were significantly different in Baltimore compared to both Riverside and Palm Springs but not between the two more arid cities (Figure 5 and Table 3B). No differences were found in the distance to centroid of taxonomic compositions among the cities or vegetation types. Regarding functional composition, only woody physiological composition showed a significant difference in variance around the centroid between parks in Baltimore and Palm Springs (Table 3A). Bold font represents a significant p-value. ES Wood refers to woody species which also have service-based functional traits recorded. EP Wood refers to woody species which also have ecophysiological traits recorded.

Regional Climate as a Driver of Park Diversity
The decrease of taxonomic alpha diversity with increasing aridity is consistent with our hypothesis of climate filtering. In arid cities like Palm Springs, the number of woody and turf grass species that can withstand physiological stress is limited compared to the horticultural species found in more mesic cities (Pearse et al., 2018; Figure 2). In recreational parks, taxonomic diversity is further filtered by park planners selecting species that provide specific desired services (Talal and Santelmann, 2019). The increased selection pressure of park management combined with extreme aridity could explain why we found a decrease in diversity while other studies found a positive correlation between species diversity of all urban trees and warmer climates . Conclusions based on broad sampling of urban vegetation differ when confining studies to a specific type of green infrastructures, such as residential lawns or recreational parks as in this study, where management practices regarding planting and watering are more specialized (Pearse et al., 2018). As increasingly arid climates appear to reduce species diversity within recreational parks, differences in climate also result in unique taxonomic and functional composition in recreational parks across cities (Figure 5). All vegetation types in recreation parks have unique compositions among our study cities. We hypothesized that to achieve desired services, the taxonomic composition will vary with minimal effect on the values of service-based traits. Unique taxonomic compositions can vary across cities due to climatedriven changes in the availability of urban vegetation , temporally through variability of nursery offerings , or depending on local consumer preferences (Conway and Vander Vecht, 2015;Roman et al., 2018). While many species found in recreational parks are not native to the region (Wheeler et al., 2017;Talal and Santelmann, 2019), the lower Shannon diversity in Palm Springs and the difference in composition among all cities, explain how horticultural vegetation physiologically responds to the native climate.
For recreational parks to maintain healthy vegetation, planted vegetation must be able to functionally acclimate to regional climate. Lower specific leaf area is often connected to adaptations that confer tolerance to drought, such as reduced wilting and ability to maintain photosynthesis in extremely arid environments (Poorter et al., 2009;De Micco and Aronne, 2012). The higher wood density seen in arid urban park trees can lower cavitation risk, and mortality caused by extreme aridity may be better avoided by these species than by species of lower wood densities (Wright et al., 2004;Savi et al., 2015). However, resistance to drought is likely a co-benefit as park managers are more concerned with planting for ease of maintenance and economic reasons, and high wood density is associated with long lived trees which would need to fewer replacements (Díaz et al., 2015;Talal and Santelmann, 2020). This ecological tradeoff between stress tolerance and growth is also found in the trait of seed dry mass. Species adapted to arid regions generally have seeds with lower dry mass to prevent desiccation before precipitation brings on germination (Westoby, 1998). Recreational park tree species found in arid regions exhibit adaptations to reduce heat and drought stress by having lower specific leaf area, seed dry mass, and higher wood density. Variation in functional traits ( Figure 3A) among cities will result in park visitors in Baltimore experiencing a different functional composition than in Riverside and Palm Springs ( Figure 5B).
The photosynthetic pathway of C4 turfgrasses allows for continuous transpiration during extreme heat and aridity, and in a well-irrigated Palm Springs park there is plenty of water for turf to maintain productivity. In a mild city like Baltimore, there are better conditions for C3 turf to highly productive yearround. However, while the physiological traits of the C3 turf in Baltimore and the C4 turf in Palm Springs were similar, the range of trait values was smaller in Palm Springs. This result is indicative of fewer varieties of C4 turf grass for planting compared to C3 varieties in nurseries throughout the country (Trammell et al., 2019).
Spontaneous species are not subject to the same environmental or preferential drivers as turf species (Niinemets and Peñuelas, 2008;Knapp et al., 2012). Yet, many varieties of spontaneous species persist in cities and are of similar functional type (Wheeler et al., 2017), resulting in similar physiological traits and taxonomic diversity. Spontaneous species maintained distinct compositions among cities, while the turf species were less dissimilar ( Figure 5A). However, unlike turf composition, spontaneous functional composition was not significantly different across cities. For spontaneous species to be found within a recreational park lawn, the species must disperse and germinate within the climate of the city, and then establish despite management and recreational activities enacted on lawns (Abu-Dieyeh and Watson, 2005;Anderson and Minor, 2019). Adaptations to regional climate may explain the variety in lawn taxonomy that we observed, while the adaptation to intense human impacts could explain the lack of difference in functional composition.

Management Preference as a Driver of Park Diversity
Depending on the specific functional trait, management preference or climate can be a greater influence on park plants. The functional responses specific leaf area in turf and seed mass in spontaneous species are potential evidence of aridity creating regions of unique taxonomic and functional compositions (Figure 4). Regarding suites of service-based traits, Baltimore's climate allows for a greater abundance of species with the highly valued trait of esthetic fall foliage, as evidenced by the Maryland Department of Natural Resources' weekly fall updates tracking the changing of the leaves (Maryland Department of Natural Resources 2019). This unique difference in a single climate driven trait can drive significant changes in total service-based trait composition in recreational parks ( Figure 5B). Similarly, while fruit production is not a preferential service-based trait by park managers, we did find a greater abundance of fruit trees in the Palm Springs and Riverside regions, which could be indicative to the legacy of citrus in southern Californian agriculture (Farmer, 2013).
In recreational parks, woody species are generally planted to provide valued service traits of shade, and esthetic appreciation Avolio et al., 2018) and turf is planted to provide the service of vegetative greenness and recreational play areas (Ayala-Azcárraga et al., 2019). Park trees are spaced far enough apart for esthetics to be appreciated while providing copious but fragmented shaded areas, creating a physical arrangement of vegetation where similar service traits are available across cities (Goodness, 2018;Talal and Santelmann, 2019). While regional aridity is correlated with park vegetation physiology, the preference for specific values of ES traits is generally similar across urban regions. Similarity of servicebased traits in recreational parks allows park visitors, who value certain services, to have similar experiences in cities regardless of regional horticulture. While our study cities exhibited a unique functional composition of woody EP and ES traits, we found no differences in FDis or compositional variation ( Figure 2B and Table 3A). We suggest this pattern represents an example of park managers' preference for a similar arrangement of function in recreation parks across cities, despite a regionally varying species diversity. Parks in Palm Springs achieve a comparable range of functional breadth to milder cities while exhibiting a unique taxonomic and function composition and a significantly smaller taxonomic diversity (Figures 2, 3, 5).
One of the defining characteristics of woody species in recreational parks is that they are not a monoculture; there is generally a variety of tree species (Jim and Chen, 2008;Nielsen et al., 2014). Moreover, surveys of urban residents have identified specific highly valued service-based traits in trees, such as shade, beauty, and fruit production . Variation in the desires of urban residents result in a large variety of species that can provide these specific services Talal and Santelmann, 2020). Therefore, a wide variety of woody species will result in a similar FDis of physiological traits, even when the distribution of species is different. While parks may provide varying amounts of services, the similarity of dispersion of service-based traits in cities with differing Shannon diversity implies management practices that select for a specific range of service functions in a park.
Woody species may not be a monoculture, but our results show that planted turf acts like a plantation. We interpret this result as evidence that the highest value service trait is that it remains green and usable for recreation and is easily maintained (Christians et al., 2016;Larson K. L. et al., 2016). If a minimal number of species available can achieve this goal, there is little need to expand planting to species with unique functions. Maintenance regimes can be standardized to the similar taxa, which fits with the stated desires of park managers to minimize maintenance costs (Chan et al., 2014;Talal and Santelmann, 2020). There are other varieties of native grasses or sedges that are viable in regions of high aridity, yet these varieties are not commonly cultivated as turf.

Urban Form as a Driver of Diversity
Interestingly, we did not find significant relationships between park area and alpha diversity, either taxonomic or functional. This is in contradiction to the recent review of diversity drivers in urban parks (Nielsen et al., 2014). The discrepancy between our study and current literature could arise from our focus on urban recreational parks, where the entire park area is actively managed only for recreation (Weems, 2016). Conversely, the review by Nielsen et al. (2014) includes parks managed for recreation, agriculture, and natural areas. Including multiple varieties of urban park incorporates multiple habitats as well, which influence species-area relationships.
While we used observations of urban park vegetation as evidence of management preferences, there are other management interactions outside of vegetation preference that can influence the diversity in a community. Fertilization, pesticide application, prioritizing play equipment, and access to local nurseries can all lead to variations in taxonomic and functional vegetation diversity (Kjelgren and Clark, 1993;Politi Bertoncini et al., 2012;Chan et al., 2014;Cavender-Bares et al., 2020). Urban soil profiles can be both heterogeneous within and among cities (Crum et al., 2016;Herrmann et al., 2018). Regional climate and human facilitation are major filters leading to taxonomic and functional diversity (Aronson et al., 2016;Pearse et al., 2018), however, future work would be inclined to explore other potential narrower urban filters on recreational park diversity.

CONCLUSION
Climate and management preferences both play key roles in determining recreational park structure and composition, through driving differences in taxonomic alpha and beta diversity while maintaining similarity in the value and distribution of service-based traits. We show that regional climate drivers affect the taxonomic diversity and composition of each city's parks. Furthermore, the number of functional strategies also reflects a stabilization of FDis among cities. By showing how Shannon diversity of woody and turf species both responds to climatic shifts, FDis of woody species and turf species diverge in this regard, we can infer woody species FDis is influenced by park management preference while turf has a strong limitation to extreme climates. Integrating woody, turf, and spontaneous vegetation with multiple metrics of diversity allows for these results and answers a call for work to incorporate multiple functional types within a singular study (Nielsen et al., 2014). To develop a comprehensive diversity framework for entire cities, future work should incorporate more cities to represent more regional climates and differences in local horticultural preference.
Following the paradigms of urban ecology "in", "of " and "for" the city, our hypothesis required a synthesis of both climate drivers and the influence of park managers' preference to understand the patterns of recreational park diversity (Pickett et al., 2016). Our study synthesizes a biotic (ecology "in" the city) and a social-ecological (ecology "of " the city) influencer of diversity and our results provide a rationale for why specific vegetation types and functions are more influenced by either climate or preference. Our results show the influence of management preference guiding diversity in recreational parks, where service-based traits and FDis (excepting turf) does not change while taxonomic alpha diversity decreases into arid regions. Using our approach, we can identify functional cobenefits that could enhance the selection of park vegetation to provide climate resiliency along with traditional park management. The baseline days of extreme heat in Palm Springs are projected to increase from 135 to 179 days >35 • C by the end of the century (Sun et al., 2015), creating opportunities for the current stewardship for service-based function to shift in these extreme cities as the availability of viable vegetation decreases. By incorporating the results of this study with local urban park planning are cities at risk of extreme heat, we can move to practicing ecology "for" the city as well.

DATA AVAILABILITY STATEMENT
The datasets generated for this study can be found in the doi: 10.6086/D1FT1R.

AUTHOR CONTRIBUTIONS
PI, DB, CS, and GJ designed the study. PI, DB, and MR collected all the data for the study. PI and DB completed the statistical analysis of the data. PI wrote the manuscript. All authors contributed to the editing of the manuscript and approved the submitted version.

FUNDING
This work was supported by NSF CBE -1444758 and CNH -1924288.