Trees outside forests in agricultural landscapes: spatial distribution and impact on habitat connectivity for forest organisms

Various species of forest trees are commonly used for ornamental purposes and are therefore frequently found in nonforest ecosystems. They constitute an important component of the so-called trees outside forests (TOF). Not much is known, however, about the drivers of TOF spatial distribution either in urbanized or in agricultural landscapes since they are generally absent from forest inventories. The present study focused on the spatial distribution of TOF across agricultural landscapes and their potential role in the dispersal of a forest pest insect, the pine processionary moth, Thaumetopoea pityocampa (PPM). All the TOF belonging to the genera Pinus, Cedrus and Pseudotsuga were considered as potential hosts and inventoried within a 22 × 22 km study window. We fitted a nonstationary Poisson process to the empirical data and used the distance to the nearest building as a covariate. Both empirical and simulated data indicated that TOF associated to human artifacts/urbanized areas constituted the main source of landscape connectivity for the PPM in the open fields under study. Because they do not account for TOF, forest inventories dramatically underestimate landscape connectivity and provide an erroneous picture of the PPM habitat distribution. We conclude that TOF, especially the ornamental component, must be taken into account when it comes to understanding forest insect landscape dynamics or genetics. The omnipresence of TOF also suggests a potentially huge role in pest dispersal and invasive species expansion.


Introduction
Habitat degradation and fragmentation is considered to be one of the primary threats to biodiversity worldwide (Balmford et al. 2003). By disrupting dispersal, fragmentation alters a critical component of metapopulation dynamics and population persistence which strongly affects biological diversity (Wiens 2001;Hanski 2011). For that reason, species conservation issues are often viewed from a landscape perspective and this is one of the reasons why landscape ecology has become increasingly relevant for management and conservation of biological diversity (Gutzwiller 2002;Lindenmayer and Hobbs 2007;Wiens 2009;Hanski 2011). One widespread management option consists of developing green infrastructures and corridors to help connecting rural and urban or periurban territories in order to limit landscape fragmentation and facilitate species dispersal (Benedict and McMahon 2006;Clergeau and Blanc 2013). Trees and shrubs are major elements of these so-called green belts. Trees outside forests (TOF) are increasingly recognized as a prominent feature of either agricultural lands or built-up areas (Bellefontaine et al. 2001;Gutzwiller 2002;Konijnendijk et al. 2005;de Foresta et al. 2013). TOF are ''trees on land not defined as forest and other wooded land,'' a definition that is dependent on how forests and woodlands are defined themselves (Kleinn 2000;de Foresta et al. 2013). In agricultural landscapes, TOF are often referred to as scattered, isolated or paddock trees (Manning et al. 2006;Gibbons et al. 2008). In cities and towns, TOF are street and roadside trees, trees planted in parks, gardens and private yards (Tyrvainen et al. 2005;Wu 2007). TOF provide various ecosystem services such as control over soil erosion, nutrient and water cycling, biodiversity conservation or pest control (Plieninger et al. 2004;Lumsden and Bennett 2005;Manning et al. 2006). Urban woods and isolated trees also provide important social, aesthetic, and economic benefits (Tyrvainen et al. 2005;McDonnell et al. 2009).
A proportion of TOF are ornamental or amenity trees grown for decorative purposes in gardens and landscape design projects. These trees are native, naturalized or recently introduced exotic species. In recent decades, the exotic component of ornamental planting has increased as the international horticultural industry has developed and changed toward a global market. This evolution is associated with important consequences on the risk of emerging plant diseases (Dehnen-Schmutz et al. 2010;Ö öpik et al. 2013). However, while TOF benefits and services are increasingly acknowledged, it is worth noting that they also provide ecosystem disservices such as the production of allergens, potential damage of infrastructure and the introduction of new pests and pathogens (Roy et al. 2012;Bigsby et al. 2014;Dobbs et al. 2014).
Despite their ubiquity, little is known about ornamental TOF spatial distribution and the ecological and/or socio-economical processes that drive it at the landscape scale. Ornamental species are species reported as being frequently planted near houses and buildings  so that their spatial distribution is likely to be linked to human settlements but, to our knowledge, quantitative data describing such association are lacking. In addition, TOF (either ornamental or not) are generally absent from national forest inventories (Kleinn 2000;de Foresta et al. 2013) and this lack of data precludes a proper assessment of their contribution to landscape connectivity with regards to associated species.
In the present study, we tackle the question of TOF contribution to landscape connectivity with regards to a forest defoliating insect, the pine processionary moth (PPM) (Thaumetopoea pityocampa). The PPM is an important forest pest in the Mediterranean regions and its geographical range is currently expanding in relation to climate change (Roques 2015). Little is known about the landscape scale spread of the PPM, particularly in non-forest areas. We worked in large open-field agricultural landscapes of the Centre region of France where available forest inventories indicate a very low coverage of tree species suitable for our model insect (IFN 2010). In this recently colonized area, we thus expected a high degree of habitat fragmentation with regards to the PPM host species (referred to as host trees outside forests, HTOF) coverage. Using systematic inventory, we georeferenced all isolated PPM host trees either located on public or private lands within a large sampling window. We examined the possible relationship between human artifacts and HTOF distribution using data from the French institute IGN (Institut National de l'Information Géographique et Forestière) using point process modeling. Because exhaustive inventories are tedious and costly, spatial simulation may complete or enhance field data (O'Sullivan and Perry 2013). We explored how the resulting model could be used to simulate HTOF distribution over a large region of open-field where no HTOF inventory was done. Finally, we used landscape metrics to quantify HTOF contribution to landscape structure and particularly fragmentation with regard to the PPM.

Model species and study site
The PPM is the most important pine and cedar defoliator in southern Europe and northern Africa Landscape Ecol (Cielsa 2011). Pinus nigra, one of its preferred hosts, has been used for large-scale afforestation and ornamental plantations throughout France. The distribution of PPM is controlled by the minimum winter temperatures (Huchon and Demolin 1970), and the species range is currently expanding northward and in elevation due to climate change (Battisti et al. 2005;Robinet et al. 2007). Since the early 90s, the PPM's range has expanded towards north of France without being noticeably affected by the presence of large wide open-field areas such as the Beauce region (Roques et al. 2015a) despite limited female dispersal ability (Robinet et al. 2012) (Fig. 1). PPM dispersal in nonforest landscapes may strongly depend on scattered host trees although large-scale field data are lacking (Roques et al. 2015b). In order to unravel the potential role of HTOF as elements of landscape connectivity for the dispersal of PPM, we designed a field survey just north of the open-field region of Beauce in the centre of France (Fig. 1). This area is amongst Europe's most productive cereal regions and thus appeared to be a sensible choice given the scarcity of coniferous forest cover, the typical habitat of the PPM. Figure 1 shows the position of the Centre region of France (Fig. 1 left)

Field sampling
An exhaustive inventory of PPM HTOF was carried out in a 22 9 22 km = 484 km 2 area located in the north of the Centre region of France (Fig. 1). All trees belonging to the genera Pinus, Cedrus and Pseudotsuga were considered as potential hosts and are hereafter referred to as HTOF for host trees outside forests. The study site was located in the southern part of an ecoregion referred to as the Drouais-Thymerais region (Conseil régional Centre 2009) where the landscape mostly consisted of wide arable lands (cereal). All motorable roads and tracks were visited during autumn and winter 2009-2010. This period of the year was preferred because sighting and identification of PPM host trees (evergreen coniferous species) is easier when deciduous trees have lost their leaves (Rousselet et al. 2013). Every individual or small group of host trees (single trees, linear groups of trees, small woodlands) was observed by eye, and with binoculars when necessary, from the road and public land. The geographic coordinates were recorded using a Garmin TM GPS12 (or the location was mapped onto a georeferenced aerial photo using Pho-toExploreur BayoÓ in case of distant observation). In so doing, we obtained an exhaustive inventory of locations of PPM HTOF. The inventory data are available in Rousselet et al. (2015). Forest inventories and land use information available in the study area The French institute IGN (Institut National de l'information Géographique et Forestière http://www.ign.fr/) provides a national forest inventory (referred to as IFN for ''Institut National Forestier'' in French) of all the stands with an area [2.25 ha, i.e. approximately corresponding to a resolution of 150 9 150 m pixels (http://inventaire-forestier.ign.fr/carto/carto/afficher Carto/28) (Fig. 2). This institute also develops and maintains a database called BD Topo Ò IGN http:// professionnels.ign.fr/bdtopo, which is a vector database documenting both geographical and administrative entities at a 1:25,000 scale over the whole French territory.

Modelling HTOF spatial distribution
We modeled the HTOF spatial distribution in order to explore the link between tree spatial pattern and the distribution of buildings over the study area. The resulting model provided a convenient way of simulating HTOF distribution over unsampled areas. HTOF were modeled by fitting an inhomogeneous Poisson process, i.e. nonstationary Poisson process (Diggle 2003), upon the observed tree location within the experimental plot (Fig. 2a). The inhomogeneous Poisson process allows the introduction of covariates in the analysis of the observed point pattern by considering an intensity function k(x). In our case, the model was where k(x) denotes the intensity function at location x and z(x) stands for the covariate at location x. We used the log of 1? the distance to the nearest object corresponding to any sort of building as the covariate. The covariate values were computed from the data of the BDO TOPO 2.1 using the SIG software Quantum GIS (Quantum GIS Development Team 2014) (Fig. 2c).
Our rationale behind the use of a Poisson process was to mimic citizen behavior with regard to ornamental tree plantations. People plant ornamental species in the close vicinity of their house, for example in their yards or in front of their property  lands distributed all over the built-up areas, especially along public thoroughfares serving the buildings (Konijnendijk et al. 2005). We account for this behavior by using a nonstationary model where the nature of the covariate directly conveys the effect of building closeness. On the other hand, not all private or public owners would plant ornamental trees suitable for the PPM and this random aspect is accounted for through the Poisson process itself.
The fitted model was evaluated by performing residual analysis (Baddeley and Turner 2006) in which the point process is transformed into a Poisson process with uniform intensity on a given interval. The deviation of the transformed point process from a homogenous Poisson process indicates a possible model misfit. We used the Stoyan-Grabarnik diagnostic tool ''exponential energy marks'' (Stoyan and Grabarnik 1991) based on the reciprocal of the conditional intensity. Each point is given a weight depending on the conditional intensity of the fitted model and the sum of the weights for all the points divided by the area of the study plot has an expected value of 1 (Stoyan and Grabarnik 1991;. Model fitting was done using a maximum pseudolikelihood fitting algorithm implemented in the R statistical software (R Core Team 2014) within the package spatstat (Baddeley and Turner , 2006). Once the model linking HTOF to the covariate was fitted, it was used to perform stochastic simulations of the point process at hand within the sampling window (Møller and Waagepetersen 2004;. A total of 1000 simulations were used to examine the model performance with regards to a set of selected landscape metrics intended to characterize landscape features. The simulated point patterns were rasterized with a resolution equal to that of the IFN database i.e. pixels of 2.25 ha (150 9 150 m, Fig. 2) using the R package raster (Hijmans 2014). The resulting binary maps were used to compute the PLAND, PD and ED metrics (described below). The results were compared to the values computed for the empirical data set. Similarly, we simulated the HTOF distribution on the scale of the Beauce region (Fig. 1c).

Landscape configuration from the PPM viewpoint
We investigated the structure of landscapes derived from (i) the IFN raster database (ii) the HTOF database, and (iii) the IFN ? HTOF database by means of 3 key landscape metrics: percentage of landscape, patch density and edge density (Turner et al. 2001). The percentage of landscape (PLAND) was used to quantify the proportional abundance of the patches that corresponded to a suitable habitat for the PPM, i.e. corresponding to pixels where at least one individual host tree was present. Patch density (PD) described landscape structure by quantifying the degree of fragmentation of PPM host distribution. In addition, we used the edge density metric (ED) to quantify the shape and complexity of the habitat patches (Turner et al. 2001). This index appeared of interest given that the PPM is known to preferentially colonize forest stand edges (Samalens and Rossi 2011). Computations were done under R with the package SDMTools (VanDerWal et al. 2014).

Descriptive statistics in the inventory plot
According to the IGN, the exhaustively sampled area mainly consisted of an intensively managed agricultural landscape with only 12 % of forest and less than 1 % of pure or mixed stands with coniferous trees. We recorded a total of 3834 individuals or small groups of PPM HTOF in the survey window (Fig. 2a), which corresponded to a density of 7.9 occurrences per 100 ha. Figure 2 shows the observed HTOF (black crosses) (Fig. 2a) and the pixels (150 9 150 m) where the IFN inventory reported the presence of host trees (Fig. 2b). Of the 3834 HTOF, 285 (7.4 %) were separated from the nearest building by a distance \10 m; 1113 (29 %) were separated by a distance ranging from 10 to 20 m; and 627 (16.35 %) were separated by a distance ranging between 20 and 30 m (Fig. 2c, d). The spatial configuration of HTOF is depicted in Fig. 3. The average separating distance between trees and their closest neighbors was \50 m for more than 70 % HTOF (Fig. 3a). Figure 3b shows experimental data illustrating the flight abilities of female moths in the laboratory (data from Robinet et al. 2012). From this barplot, it can be seen that more than 54 % of female moths fly between 0 and 1 km; 16.21 % fly a distance between 1 and 2 km, while this proportion is 13.51 % for distances ranging from 2 to A B Fig. 3 Spatial configuration of HTOF. a Frequency distribution of the distance between each HTOF and its closest neighbor. b Frequency distribution of the flight capacities of female PPM (data from Robinet et al. 2012) Landscape Ecol 3 km. These flying abilities are compatible with the inter-HTOF distances reported above.
The percentage of landscape associated with PPM potential hosts appeared to be very limited when estimated on the basis of the IFN data (1.5 %), while it was more than 10 times higher when measured from observed HTOF (12.4 %) ( Table 1). In addition, both patch density and edge density were strongly underestimated when measured from the forest inventory alone (Table 1).

Inhomogeneous Poisson model and landscape distribution of HTOF
The coefficients of the fifth degree model were a = -9.40; b 1 = -18.77; b 2 = 16.30; b 3 = -4.99; b 4 = 0.64; b 5 = -0.029. The Stoyan-Grabarnik diagnostic tool divided by the plot area was 1.005526, which revealed a good fit. The values of the PLAND, PD and ED metrics for the empirical and simulated data set are given in Table 1. On average, the landscape patterns resulting from the model slightly overestimated the percentage of landscape (PLAND) where HTOF are present (16.8 %) as compared to the observed value (12.4 %). Similarly, patch density and edge density were slightly overestimated (Table 1). Figure 4a shows the binary maps of the HTOF distribution over the Beauce region available from the IFN (see also Fig. 1). Figure 4b gives an example of a simulated HTOF distribution over the same region. Table 2 gives the landscape metric values computed from the IFN data, the simulated HTOF distributions and the sum of these data sets. The coverage of trees suitable for the PPM given in traditional forest inventory (IFN) appeared to be markedly lower than the value estimated from the simulated distributions. Similarly, both the patch and the edge density were higher according to the model simulations.

Spatial distribution of TOF and the prominence of ornamental trees
The results reported here unambiguously indicate that most of the tree formations suitable for the PPM in our inventory area are TOF. While native broadleaved TOF located within or bordering the fields are rare and considered relictual in the intensively managed agricultural landscapes (Plieninger 2012), the PPM host tree species were not rare and appeared to be scattered across the landscape. Interestingly, no pine, cedar or Douglas fir (i.e. the host species of the PPM referred to as HTOF in this study) is native to the survey region even if natural regeneration sometimes occurs, especially in unused lands. The trees considered here are planted near human artifacts such as buildings (public land and residential, commercial and industrial areas) for decorative purposes. This leads to a strong spatial relationship between trees and buildings. As a consequence, 70 % of the inventoried tree formations were found in 4.5 % of the territory defined as urbanized areas in the BD CARTO database (BD Topo Ò IGN 2013). If HTOF are locally aggregated in the vicinity of buildings, the villages are evenly distributed in the study plot (Fig. 2) and this leads to a regular distribution tree patches throughout the landscape with an average inter-tree distance well below the flight capacity of the PPM. As a consequence, the landscape functional connectivity clearly appears high enough to allow the expansion of the moth across that open-field agricultural landscape. Such an even host tree Observed HTOF refers to observed host species in the field. Simulated HTOF refers to 1000 host tree distributions simulated in the inventory plot. IFN denotes the inventory data available from the IFN database. IFN ? HTOF represents the sum of the traditional inventory data plus a component standing for the HTOF contribution to landscape structure. Standard errors are given in parentheses Landscape Ecol Fig. 4 Binary mapping of the regional distribution of host trees suitable for the pine processionary moth (Thaumetopoea pityocampa) in the ecoregion of ''Beauce'' (see Fig. 1). a Data available for the French National Forest Inventory (IFN, currently grouped with the Institut National de l'Information Géographique et Forestière http//:inventaireforestier.ign.fr). b Example of simulated HTOF distribution. The raster resolution is of 2.25 ha, i.e. 150 9 150 m pixels in both maps distribution was previously unsuspected because our knowledge of coniferous species distribution in the region relied solely on the available forest inventory databases. Beyond the case of the PPM and its host trees, our results suggest that accounting for TOF in forest inventories would be important for many ecological studies especially when dealing with metapopulations or metacommunities (Hanski and Gaggiotti 2004;Holyoak et al. 2005). The present work focuses on a set of tree species mainly used for ornamental purposes. As such, the drivers of their dynamics fall largely within social forces such as consumers' choice and market fluctuations. In 2010, for instance, 475.2 million outside plants were purchased by the French people (FranceA-griMer 2011). Of these plants, 52.8 million were trees and shrubs including 15.1 million small fruit bushes, 4.2 million fruit trees and 33.5 million ornamental trees and shrubs. In the Beauce region, these ubiquitous trees and shrubs are planted in various types of environment depending on the buyer's location and cause lasting changes in local vegetation. It should be noted that some forest tree species are commonly used as ornamental trees, especially Pinus nigra, the preferred host of the PPM. In other words, citizens and local communities have a strong impact on host tree availability for the PPM through the choice of ornamental species. Addressing this point is beyond the scope of the present paper but it is worth noting that understanding the drivers of ornamental TOF dynamics would necessarily involve economic and social sciences.

Point process models for TOF
The model satisfactorily fitted the data, although simulated point patterns proved to slightly underestimate tree dispersion as indicated by the overestimation of PLAND, patch and edge density. A possible improvement would be to deepen the level of details taken into account when considering urban landscapes, for instance by distinguishing different types of buildings with regard to human occupation and function (farms, houses, public buildings…). In agricultural open-field landscapes like the Beauce region considered here, we suspect that ornamental trees are often planted in residential areas and in villages, but not at the same level in farmhouses and outbuildings. The pattern of ownership in urbanized areas may strongly affect the selection of ornamental species, their density and spatial arrangement in cities (Pham et al. 2011). Again, this goes beyond the scope of the present paper and a proper consideration of these aspects would require broadening the scope of our approach by including socioeconomic components.
The strength of point process modeling is that it allows simulating point patterns corresponding to the model fitted to observed data over an unsampled region (Diggle 2003;O'Sullivan and Perry 2013). This approach could be used to draw maps of TOF at large geographical scales where exhaustive data might be impossible to collect for obvious practical reasons. One very important limit of such an approach is that a model can only be used for simulating data in a landscape of similar nature and assuming that TOF generating processes at work are identical. Given this intrinsic limitation, simulations may provide usefuland cheap-information. TOF are generally assessed using stratified sampling and remote sensing-based approaches (Holmgren et al. 1994;Corona and Fattorini 2006;Tokola 2006;Lam et al. 2011;Baffetta et al. 2011;Lister et al. 2012;Meneguzzo et al. 2013;Westinga et al. 2013). A possible area of research would be to explore how point process models could be used to optimize field sampling and/or how point process simulation could be used to map TOF and complete the data available in forest inventories.
Landscape connectivity, pest spatial dynamics and management Data from the available forest inventory suggested that the Beauce ecoregion offers poor habitat resources for the PPM due to a very limited forest cover, especially in terms of conifers (IFN 2010) (Fig. 2a, b). Despite their very low density in comparison with forest stands, ornamental TOF associated with built-up areas provide numerous potential relays for the PPM females during the range expansion. Percentage of landscape, patch density and edge density all indicate that landscape permeability is strongly underestimated when only considering forest stands. Ornamental TOF therefore constitute an unsurprisingly, although rarely explicitly considered, important network corresponding to low tree densities. The resulting landscape connectivity with regards to associated invertebrate or pathogen species is thus of high functional importance. In the case of the PPM, given its flight capacities (Fig. 3), it appears that there is no barrier to spread of that forest insect in large open-field agricultural regions. More generally, models for pest dynamics may benefit from an explicit account for TOF as habitat resources and dispersal corridors. For example, taking TOF into consideration in landscape genetics may lead to a refined and more detailed perception of landscape constraints towards gene flows, hence a better understanding of species spatial ecology and spatial genetic variation (Kerdelhué et al. 2015;Roques et al. 2015b). Urban forests harbor a large number of insect pests also present in planted or natural forests, orchards or forest nurseries (Tello et al. 2005). In addition to native and naturalized pest species, urban trees are exposed to an increasing number of exotic pests and pathogens due to international trade (Smith et al. 2007;Tubby and Webber 2010;Kowarik 2011). Before spreading across large areas, most successful invasive species are firstly introduced in urbanized areas that act as hubs of international transport and then constitute major gateways by which species spread (Zippeter 2008;Colunga-Garcia et al. 2010). Urban areas are also hot spots of plant species richness, either for native or introduced species (Kowarik 2011) which, in turn, could be a facilitating factor of establishment success of insect pests. For example Colunga-Garcia et al. (2010) showed that the risk of invasion increased near to urbanized areas. Spread of invasive species can also be facilitated by green infrastructures that are intended to improve urban-rural gradients with regards to biodiversity management. So far, such ecosystem disservice has received less attention than beneficial effects (but see Säumel and Kowarik 2010). In the case of the PPM, P. nigra is amongst its preferred pine hosts and has massively been used as an ornamental tree in private and public lands. An obvious management option would be to discourage massive planting of P. nigra in areas frequented by public at risk like schools. Quantifying the contribution of different types of ownership (private yards, industrial and commercial areas managed by landscape gardeners, street trees or public green spaces managed by local public authorities) to the dynamics of ornamental tree plantings would help develop a better understanding of the dynamics of green infrastructures and their role in the spatial ecology of pests.