Next Article in Journal
NF-kB Regulation and the Chaperone System Mediate Restorative Effects of the Probiotic Lactobacillus fermentum LF31 in the Small Intestine and Cerebellum of Mice with Ethanol-Induced Damage
Next Article in Special Issue
Competitive Advantage of Broussonetia papyrifera Growing in a Native Area as Suggested by Structural Diversity
Previous Article in Journal
Decoding the Behavior of a Queenless Colony Using Sound Signals
Previous Article in Special Issue
Effects of Clipping an Invasive Plant Species on the Growth of Planted Plants of Two Co-Occurring Species in a Greenhouse Study
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Spatial Patterns of Frangula alnus (Rosales: Rhamnaceae): Implications for Invasive Plant Management

by
Jennifer Greenleaf
1,
Roghaiyeh Karimzadeh
1,2 and
Yong-Lak Park
1,*
1
Division of Plant and Soil Sciences, West Virginia University, Morgantown, WV 26506, USA
2
Department of Plant Protection, Faculty of Agriculture, University of Tabriz, Tabriz 5166614888, Iran
*
Author to whom correspondence should be addressed.
Submission received: 30 September 2023 / Revised: 26 October 2023 / Accepted: 30 October 2023 / Published: 1 November 2023
(This article belongs to the Special Issue Biology, Ecology and Management of Invasive Alien Plants)

Abstract

:

Simple Summary

Glossy buckthorn (Frangula alnus) (Rosales: Rhamnaceae) is an invasive shrub from Europe, threatening native vegetation in open and disturbed habitats. To determine the spatial distribution patterns of glossy buckthorn, the density and presence of F. alnus were recorded in 1412 transect sample points located on four sites in the Allegheny National Forest, Pennsylvania, USA. Geostatistical analyses were used to determine spatial autocorrelation among individual F. alnus samples and to measure their spatial dependency. Spatial patterns were examined using spatial analyses by distance indices (SADIE), and buffer analyses were conducted to determine the relationships between proximity to unpaved forest roads and both F. alnus density and presence. The results of this study showed that the density and presence of F. alnus decreased as the distance from the forest road increased. We also found a significant spatial aggregation of F. alnus in all four study sites. These findings imply that spatially targeted management of F. alnus would be possible by locating the patches of F. alnus and applying control measures only where needed.

Abstract

Glossy buckthorn (Frangula alnus) (Rosales: Rhamnaceae) is an invasive shrub from Europe that has been invading North America for over a century and threatening native vegetation in open and disturbed habitats. The treatment of F. alnus is currently restricted to the roadside, suggesting any individual F. alnus residing within the forest would be left unmanaged and would continue to spread in the area. This research was conducted to determine the spatial patterns and relationship of F. alnus with forest roads. The presence and density of F. alnus at 1412 sample points were recorded on four sites in the Allegheny National Forest, Pennsylvania, USA. Buffer analyses were conducted along roads to determine the relationship between F. alnus density and proximity to forest roads. Geostatistics and spatial analysis by distance indices (SADIE) were used to characterize the spatial pattern of F. alnus. Results of this study showed that F. alnus was spatially aggregated and resided beyond forest roads. Both the density and presence of F. alnus decreased as the distance from the forest road increased. These results imply the potential for precision management of F. alnus by locating and managing only where F. alnus presents.

1. Introduction

Glossy buckthorn (Frangula alnus) (Rosales: Rhamnaceae) is an invasive shrub from Europe and has been invading North America for over a century [1]. In Europe, F. alnus was historically used to create superior charcoal for gunpowder production, is used as an antibacterial and antioxidant medicine [2], and plays a natural role in the ecosystem where it has coevolved with other wildlife [3]. In the late 1800s, it was intentionally introduced as an ornamental shrub to Canada. It then escaped cultivation and was first discovered growing wild in London, Canada in 1898 [4]. Since then, F. alnus has spread throughout Southwestern Canada and the Northeastern USA, proving itself to be an aggressive pest as it forms dense monospecific patches in disturbed areas, along roads, and throughout fields and clearings [5]. As of 2023, 29 states in the USA harbor this invasive plant [6].
Multiple factors contribute to the success of F. alnus in North America [7,8,9]. It appears to lack the natural herbivores that it experienced in Europe. Our observations of F. alnus populations in Belgium and Germany have noted that F. alnus are often small and sparse and that the leaves often show signs of pathogens and insect herbivory. In contrast, in the USA, F. alnus forms large, dense patches and the leaves appear to be largely undamaged. Known herbivores such as Zygina suavis (Hemiptera: Cicadellidae) and Gonopteryx rhamni (Lepidoptera: Pieridae) [10,11] could limit the growth and reproduction of F. alnus in its native habitats, but these specialist herbivores are not known to exist in North America.
One of the major regions in the USA where F. alnus invaded and established very successfully is the Allegheny National Forest, Pennsylvania, USA. The presence of F. alnus in the Allegheny National Forest was first reported in 1989. At that time, local forest management favored F. alnus due to timber harvests creating open canopies. Traditionally, detection of F. alnus is done by roadside surveys, stand exams, and surveys of areas marked to receive timber management. In 2007, the USDA Forest Service conducted an initial survey of F. alnus and a closely related species Rhamnus cathartica (Rosales: Rhamnaceae) in the Allegheny National Forest [12]. The survey area covered 1420 ha, and only 20 ha did not have either buckthorn species, and 55% of the survey area had a moderate or heavy infestation. Light feeding damage was also noted due to deer browsing, but it was also mentioned that this seemed to not affect invasion [12].
According to the Forest Service database entry for F. alnus in the Fire Effects Information System (https://www.fs.usda.gov/database/feis/plants/shrub/fraaln/all.html, accessed on 20 August 2023), the most effective method of control may be to target seedlings and to stop new invasions from occurring. Reducing disturbances is also recommended because F. alnus responds positively to openings and disturbances [13] although the complete eradication of F. alnus has historically not been successful [14,15]. According to Converse [16] and Schuster et al. [17], F. alnus may be managed by cutting, mowing, girdling, excavating, burning, and underplanting. Herbicide treatment of girdled stems and cut stumps, prescribed burns followed by foliar herbicide, and the removal of roots from cut areas also can be used to control F. alnus [18]. Previous studies [17,19] also suggested the complete removal of F. alnus for initial 2-year management to prevent its reinvasion. In the Allegheny National Forest, foliar herbicides (e.g., glyphosate), pre-emergent herbicides (e.g., sulfometuron methyl), and manual cutting or mastication are used to control F. alnus [20]. However, these practices of detection and treatment have been used for F. alnus only alongside forest roads where they are visually detectable, although F. alnus could be growing inside the forest far from the road. This indicates that any individual F. alnus residing within the forest would be left unmanaged and would continue to spread in the area. To maximize the management of F. alnus, it is vital to understand its spatial distribution patterns and how it is associated with roadways in the forest.
Determining the spatial distribution pattern of invasive plants is one of the fundamental measures providing important information about their population dynamics and dispersal ecology. Subsequently, visualizing the spatial distribution of the invasive plants on a map would be helpful for their population management [21]. Geostatistics is a statistical procedure that uses sample values and locations simultaneously to analyze spatial dependency and to estimate values at unsampled locations (i.e., interpolation) [22]. To quantify spatial dependency, the semivariance of paired samples is plotted against the lag or h, the separation distance between sample pairs (i.e., semivariogram) [22]. Geostatistics employs an interpolation method, called kriging, by using the spatial structure of data (i.e., semivariogram) to estimate values at unsampled locations [22]. Although semivariograms determine the spatial structure and spatial autocorrelation, they do not statistically determine the spatial distribution patterns such as aggregation, random, or uniform. Spatial analysis by distance indices (SADIE) is a methodology explicitly developed for determining spatial patterns by analyzing spatial count data [23]. In SADIE, the spatial pattern of samples is determined by calculating the index of aggregation (Ia) based on the distance to regularity [23], the minimum effort that the individuals in a sample would need to expend to move to an arrangement where there was an equal number in each sample unit. In addition, SADIE can be used to generate a clustering map, called a red-blue map, showing patches (i.e., spatial clusters of significantly higher count) or gaps (i.e., spatial clusters of significantly lower count). Detailed descriptions of this method and its applications can be found in published articles [23,24,25]. Geostatistics and SADIE together could be a novel way to determine the spatial structure of F. alnus, statistically test spatial aggregation, and generate interpolated maps of F. alnus distribution which can be used for developing a new strategy for management.
The objectives of this study were to determine (1) the relationship of F. alnus distribution with the proximity to forest roads and (2) the spatial distribution patterns of F. alnus using advanced geostatistical methods. Understanding the spatial distribution patterns and relationship of F. alnus with forest roads could increase both the detection and the successful treatment of F. alnus.

2. Materials and Methods

2.1. Study Sites

This study was conducted in the Allegheny National Forest near the borough of Ridgway, Pennsylvania, USA in 2021–2022. Four 2.5-ha sites with known infestations of F. alnus were selected for this study (Figure 1). The geocoordinates of the center of the four sites are as follows: site 1 (41.43704, −78.79137); site 2 (41.45674, −78.80247); site 3 (41.47875, −78.77288); site 4 (41.485216, −78.741848). The Allegheny National Forest is composed mostly of even-aged mature hardwood forests with closed canopies. Sites 1 and 2 consist of deciduous forests and unpaved forest roads. Site 3 is split down the middle by a powerline with an unpaved road and includes patches of grasslands and a deciduous forest. Site 4 is split by an unpaved road; on one half is a deciduous forest, and on the other a patch of coniferous trees and marshes.

2.2. Sampling Method

A total of 1412 sample points were placed every ca. 8 m on transects in four sites. Transects ran parallel to the longest borders of each site and also ran roughly parallel to roads within the sites. At each sample point, geo-coordinates and F. alnus data were recorded by using Avenza Maps version 4.2.2 (Avenza Systems Inc., Toronto, ON, Canada) and a GPS (Garmin GPSMAP® 64st, Garmin Ltd., Olathe, KS, USA).
Because many F. alnus inside the forest were small and short, it was impossible to take measurements of diameter at breast height (DBH). Therefore, we developed an index (called Buckthorn Index) that can include the measurement of small and short F. alnus equivalent to DBH. The Buckthorn Index (BI) was measured in the following manner. Within a circular area of 2.2 m2 at each sample point, all F. alnus individuals were given a value based on DBH and the number of F. alnus, but instead of measuring DBH directly, values were estimated by a point system. Each F. alnus within the 2.2-m2 area at each sample point that also reaches breast height was given one BI point for every 5 mm of stem diameter. Any F. alnus individuals that were too short to reach breast height were given a value of 1; this allowed for the inclusion of the many F. alnus saplings present in this study. The scores of BI values of all F. alnus found at each sample point were summed to represent the BI value for the sample point. Data from Avenza Maps was exported as a shapefile. Data from the Garmin GPS was exported to DNR GPS (Minnesota DNR, Saint Paul, MN, USA) and was then converted into a shapefile. Both shapefiles were combined in ArcGIS Pro Version 2.8.0 (Environmental Systems Research Institute, Redlands, CA, USA).

2.3. Buffer Analysis with GIS

The georeferenced BI data and site boundary shapefiles were uploaded to ArcGIS Pro. Shapefiles were created to designate all roads within and surrounding the four sites. Buffers were added site by site to the road shapefiles in increments of 5 m until all data points were within a buffer that represented a specific range of distance from the road (e.g., 0–5 m, 5–10 m, and so on). Regression analysis [26] was used to determine the relationship between BI values and the distance from the road. This buffer analysis was also done with the presence dataset (i.e., presence/absence of F. alnus at each sample point).

2.4. Geostatistical Analysis

To determine the spatial dependency of F. alnus distribution at each site, semivariograms were created for both the BI data and presence data. Semivariance was defined by the following equation [22]:
γ ( h ) = 1 2 N ( h ) i = 1 N ( h ) [ Z ( x i ) Z ( x i + h ) ] 2
where γ (h) is semivariance at distance between sample pairs h, N(h) is the number of sample pairs separated by lag distance h, and Z(xi) and Z(xi + h) are derived from sample values at point i. Active lag distances were chosen based on optimizing the lowest residual sum of squares (RSS) and highest R2 values. The uniform lag class distance intervals were between 7–9 m depending on the site. Three semivariogram parameters were used to describe the spatial structure of F. alnus distribution: range (A0, maximum distance at which two points are spatially correlated), sill (semivariance value at the range), partial sill (C, the difference between sill and nugget), and nugget (C0, y-intercept) [22].
The degree of spatial dependency (DSD) was measured by using the following formula [27]:
D S D = ( C C 0 + C ) × 100
Spatial dependency is considered weak when DSD < 25%, moderate when 26% < DSD < 75%, and strong when DSD > 76% [27]. The assumption that each lag class had more than 30 sample pairs was met, and eight analyses (4 sites, each with BI and presence values of F. alnus) were conducted.
Based on the assumption that forest roads would affect the distribution patterns of F. alnus, current detection and management have been focused on F. alnus only alongside forest roads. In such a case, the data would exhibit anisotropy (i.e., directional effect). Anisotropy was tested by constructing semivariograms in four directions (0, 45, 90, and 135°) with a 22.5° offset tolerance [22]. The anisotropy factor (i.e., the quotient between the minor and major axes of anisotropy) was used to determine the existence of directional effects. Where the anisotropy factor was <0.5, anisotropy would be used in the analyses [28].
Once semivariogram models were developed, the parameters were used for estimating BI and presence values at unsampled locations by using kriging (i.e., interpolation to create distribution maps). Cross-validation was used for assessing the precision of predictions and selecting the best model [28]. The models with mean prediction error close to zero, mean squared deviation ratio close to 1, and smallest root-mean-square prediction error were selected to generate distribution maps by using ArcGIS Pro. All geostatistical analyses were conducted using GS+ Version 9.0 (GammaDesign Software, Plainwell, MI, USA).

2.5. Spatial Analysis by Distance Indices (SADIE)

Spatial patterns of BI and the presence of F. alnus were investigated by calculating the index of aggregation (Ia) based on the distance to regularity. A value of overall Ia = 1 indicates a random spatial distribution; Ia > 1 indicates aggregation of observed counts; and Ia < 1 indicates a uniform spatial pattern. The statistical significance of spatial aggregation was determined by Pa values, the probability level of Ia. The Pa < 0.05 indicated that spatial aggregation was statistically significant. Clustering indices ( v i or v j ) were defined for each sample point to determine whether it belonged to a patch (significantly higher than the mean) or a gap (significantly lower than the mean). Overall spatial clustering was determined by v i ¯ (the mean value of v i of all patches) and v j ¯ (the mean value of v i of all gaps) values and their associated probabilities, P v i ¯ and P v j ¯ , respectively. Values of v i > 1.5 indicated that the sample point is within a patch and values of v j < −1.5 showed that the sample point is within a gap. Values of P v i ¯ < 0.05 and P v j ¯   < 0.05 indicated significant clustering into patches and gaps, respectively [23,29]. The analyses were conducted by using SADIEShell Version 2.0. The v i and v j values were used to produce red-blue maps and determine the locations of patches and gaps [24]. The maps were generated using the inverse distance weighted interpolation method in ArcGIS Pro.

3. Results

The average BI values of F. alnus in all four sites were 16.83 ± 4.28 (SEM), and maximum BI values at sites 1, 2, 3, and 4 were 82, 43, 109, and 105, respectively. Approximately 75% of sample points had at least one F. alnus and site 3 had the highest BI values (Table 1).

3.1. Analysis

The results of the buffer analysis and regression analysis showed that the presence of F. alnus (F = 51.29; d.f. = 1, 11; p < 0.001) and BI values (F = 49.75; d.f. = 1, 11; p < 0.001) decreased as the distance from the road increases (Figure 2). This indicates a higher density and presence of F. alnus closer to the forest road, but F. alnus still resided within the forest far from the forest road.

3.2. Geostatistical Analysis

Anisotropic semivariograms of both BI and presence values showed very low directional variations with an anisotropic factor > 0.5 for all datasets (Table S1). Therefore, we generated isotropic semivariograms for both BI and presence values. The semivariograms of presence values for F. alnus indicate that sites 1 and 4 fit exponential models, site 2 fit a linear model, and site 3 fit a nugget model. The ranges of presence values ranged from 4 to 204 m. DSD ranged from 0 to 87.94%; site 3 had a weak spatial dependency while sites 2 and 4 had a moderate spatial dependency, and site 1 had strong spatial dependency. The semivariograms of BI values (Table 2) indicate that site 1 fitted a spherical model while sites 2, 3, and 4 fit exponential models. These models indicate the presence of spatial dependency or spatial autocorrelation. The ranges of BI value ranged from 31 to 341 m, the distance of spatial dependency between two sample points. DSD ranged from 52.3–70.0%, indicating that all four sites exhibited moderate spatial dependency.
Interpolated maps of BI values by using kriging revealed that sites 1, 2, and 3 showed apparent positive associations of F. alnus density with proximity to the forest road (Figure 3). Only site 3 showed an apparent association with proximity to the forest road when the presence of F. alnus was mapped.

3.3. SADIE

The SADIE results of both BI and presence values on all four sites revealed significant spatial aggregations (i.e., Ia > 1 and Pa < 0.01) (Table 3). Eight red-blue maps derived from vi and vj values were generated using inverse distance weighted interpolation (Figure 4). No obvious correlation between BI values and the road was found in sites 1 and 4 while sites 2 and 3 showed obvious spatial associations of the presence of F. alnus with the roads.

4. Discussion

Although F. alnus has invaded North America since the late 1800s, this is the first study to determine the spatial patterns of F. alnus in North America. In Europe, the only study that examined the spatial structure of F. alnus distribution was conducted by Mosner et al. [30] who investigated the population genetic structure of F. alnus by relating patterns of genetic diversity with spatial distribution across Germany. Their spatial analyses showed significant spatial dependency or autocorrelation between individuals in an F. alnus population, but no relationship existed between environmental factors such as temperature and precipitation and the genetic distribution of F. alnus. The results of our study conducted at a much finer spatial scale showed that F. alnus exhibited spatial dependency and significant spatial aggregation, according to semivariograms and SADIE, respectively. In addition, our field observations noted that there were often thick patches of F. alnus growing in a thin line (approximately 1–3 m in thickness) on the edges of many forest roads and clearings (Figure 5). These linear patches of F. alnus are currently subject to control with mechanical control followed by herbicide application in the Allegheny National Forest where this study was conducted.
F. alnus in North America has a strong competitive advantage over native understory vegetation due to the lack of natural enemies, vegetative spread from underground runners, and tolerance to a wide range of soil conditions [10,11,31,32,33]. F. alnus aggressively outcompetes native regrowth by allelopathy and creates dense patches that inhibit any other understory plants from growing [5,32,33]. In addition, F. alnus maintains the unique ability to both leaf out before other deciduous trees and to retain its leaves longer into the fall than other deciduous trees [31]. These ecological characteristics equate to a longer growing season, ultimately enabling this invasive shrub to take over vital resources and crowd out native undergrowth with its vigorous maturation. F. alnus invades by taking advantage of gaps in forest canopies and along roadways. Burnham and Lee [34] found that F. alnus was 96 times more abundant in logged areas than in uncut plots and that forest disturbances benefit its growth and survival. Cunard and Lee [35] found that F. alnus loses its competitive advantage when surrounding shade-tolerant trees have greater basal areas and when sunlight is limited. They also suggested that F. alnus could be controlled by simply letting the forest advance to older successions; however, this would not work alongside roads where vegetation is regularly cleared to make room for vehicles in the forest. Therefore, we anticipated that F. alnus would exist in higher densities alongside the disturbed roadways, and falter under the forest canopy because F. alnus dispersion could be facilitated by vehicles and other machinery, or that birds may be utilizing the canopy gaps generated by roads to access F. alnus. The results of our study indicate that F. alnus density and its presence decrease as the distance from the road increases in the forest. Because the current management of F. alnus is limited to the side of the road where these invasive shrubs can be spotted from within a vehicle, the current method for F. alnus management may not be sufficient, given there are individual F. alnus outside the visible range from the forest road. We have found that many of them have successfully invaded under the tree canopies inside the forest.
Environmental factors such as edaphic and topographic factors and types of surrounding vegetation can influence the distribution pattern of F. alnus. Our study sites included habitats with both deciduous and coniferous trees, and we did not find their spatial relationship with F. alnus except forest roads. However, we noticed that F. alnus within the forest canopy and away from the road had fewer fruits, though many were >1 m tall. Future studies may seek to quantify this difference, with the importance being that if these individuals are not reproductive, they would be less of a threat and therefore less of a priority for treatment. It should be noted that the distribution of a plant is reliant not only on local environmental factors such as soil nutrient conditions and climate [7,9] but also on the species’ vegetative spread and seed dispersal methods. The seeds of F. alnus are typically bird-distributed, and it is reasonable to believe that a bird could travel between all four sites with the fruit; this would suggest correlations of both the genetics and the distribution of F. alnus between sites.
Based on our findings in this study, there are three major implications for F. alnus management. Firstly, although a higher density of F. alnus was found near forest roads, there still are patches of F. alnus inside the forest or hard-to-access areas where field surveyors generally do not check. Finding and managing such remote populations of F. alnus would prevent the spread and re-invasion of F. alnus after a successful control [19]. Current drone technologies including sensors and high-resolution satellite imagery could help find F. alnus located remotely or in hard-to-access areas [21]. Secondly, the strong spatial aggregation of F. alnus revealed from our study implies the possibility of spatially targeted invasive plant management. By detecting patches of F. alnus both roadside and inside the forest, control measures can be applied site-specifically. Such site-specific pest management generally increases pest management efficiency and reduces the number of herbicide inputs, application costs, and negative effects on surrounding wildlife [36]. Spatial distribution maps such as the ones generated in this study could be used to set management zones for site-specific pest management [36]. Lastly, the results of this research provide recommendations for how far from the road the control measures (such as herbicide, mechanical removal, or both) need to be applied to meet a certain goal under the allowable resources. We suggest treatments reach approximately 25, 60, and 70 m into the forest from the road to obtain goals of 50%, 90%, and 95% of F. alnus control, respectively, based on our regression analysis (Figure 2).

5. Conclusions

This is the first exploration of the spatial patterns of F. alnus on an individual plant scale in North America. Using advanced geostatistical methods (i.e., semivariograms and SADIE together), our study showed negative spatial associations of density and presence of F. alnus with the proximity to forest roads while the management of F. alnus is currently limited to what can be easily seen from the road. This information can improve the success and efficiency of F. alnus management. Monitoring efforts should consider that F. alnus is present inside the forest in areas that are not currently being reached by strict roadside treatments. Leaving these individuals untreated can lead to reinvasion in the area. This study also indicates that F. alnus is spatially aggregated, so we suggest spatially targeted detection and management of F. alnus.

Supplementary Materials

The following supporting information can be downloaded at: https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/biology12111393/s1, Table S1. Parameters of anisotropic semivariograms for Buckthorn Index (BI) and percent presence of F. alnus.

Author Contributions

Conceptualization, Y.-L.P.; methodology, J.G. and Y.-L.P.; software, J.G. and R.K.; validation, J.G. and Y.-L.P.; formal analysis, J.G. and R.K.; investigation, J.G., R.K. and Y.-L.P.; resources, Y.-L.P.; data curation, J.G., R.K. and Y.-L.P.; writing—original draft preparation, J.G.; writing—review and editing, J.G., R.K. and Y.-L.P.; visualization, J.G., R.K. and Y.-L.P.; supervision, Y.-L.P.; project administration, Y.-L.P.; funding acquisition, Y.-L.P. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the USDA NIFA AFRI Foundational and Applied Science (2021-67014-33757) and the Hatch Project (WVA00724) of the West Virginia University Agriculture and Forestry Experiment Station.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data will be available upon request.

Acknowledgments

We would like to thank Bella Hubbard, Shelby Meador, Vicki Kondo, and Kushal Naharki at West Virginia University for their help with the field survey. We also thank Richard Turcotte, Chris Hayes, Andrea Hille, Heather Smith, and Josh Hanson at the USDA Forest Service for their help with field coordination for the project.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

References

  1. Wafer, A.; Culley, T.M.; Stephens, K.; Stewart, J.R. Genetic comparison of introduced and native populations of common buckthorn (Rhamnus cathartica), a woody shrub introduced into North America from Europe. Invasive Plant Sci. Manag. 2020, 13, 68–75. [Google Scholar] [CrossRef]
  2. Vuletić, S.; Bekić, M.; Tomić, S.; Nikolić, B.; Cvetković, S.; Ganić, T.; Mitić-Ćulafić, D. Could alder buckthorn (Frangula alnus Mill) be a source of chemotherapeutics effective against hepato-and colorectal carcinoma? An in vitro study. Mutat. Res./Genet. Toxicol. Environ. Mutagen. 2023, 892, 503706. [Google Scholar] [CrossRef]
  3. Kremer, D.; Kosalec, I.; Locatelli, M.; Epifano, F.; Genovese, S.; Carlucci, G.; Končić, M.Z. Anthraquinone profiles, antioxidant and antimicrobial properties of Frangula rupestris (Scop.) Schur and Frangula alnus Mill. bark. Food Chem. 2012, 131, 1174–1180. [Google Scholar] [CrossRef]
  4. Catling, P.M.; Porebski, Z.S. The history of invasion and current status of glossy buckthorn, Rhamnus frangula, in southern Ontario. Can. Field Nat. 1994, 108, 305–310. [Google Scholar]
  5. Boettcher, T.; Gautam, S.; Cook, J. The impact of invasive buckthorn on ecosystem services and its potential for bioenergy production: A Review. J. Sustain. For. 2023, 42, 218–240. [Google Scholar] [CrossRef]
  6. The University of Georgia—Center for Invasive Species and Ecosystem Health. EDDMapS. Early Detection & Distribution Mapping System. Available online: https://www.eddmaps.org/distribution/viewmap.cfm?sub=5649 (accessed on 29 September 2023).
  7. Bibaud, H.A.; Ducey, M.J.; Simmons, M.J.; McCarthy, H.C.; Bibaud, R.N.; Lee, T.D. Factors limiting the success of invasive glossy buckthorn (Frangula alnus) in New Hampshire’s eastern white pine–hardwood forests. For. Ecol. Manag. 2020, 474, 118345. [Google Scholar] [CrossRef]
  8. Eisenhaure, S.E.; McCarthy, H.C.; Jenna, N.; Giguere, H.; Symonds, C.J.; Lee, T.D. Effects of turf, leaf litter, and soil compaction on emergence and establishment of invasive glossy buckthorn (Frangula alnus). For. Ecol. Manag. 2021, 484, 118933. [Google Scholar] [CrossRef]
  9. Custodio, T.; Houle, D.; Girard, F. Impact of environmental conditions on seed germination of glossy buckthorn (Frangula alnus (Mill)) in Eastern Canada. Forests 2023, 14, 1999. [Google Scholar] [CrossRef]
  10. Gassmann, A.; Tosevski, I.; Skinner, L. Use of native range surveys to determine the potential host range of arthropod herbivores for biological control of two related weed species, Rhamnus cathartica and Frangula alnus. Biol. Control 2008, 45, 11–20. [Google Scholar] [CrossRef]
  11. McKay, H. Egg-laying requirements of woodland butterflies; brimstones (Gonepteryx rhamni) and alder buckthorn (Frangula alnus). J. Appl. Ecol. 1991, 28, 731–743. [Google Scholar] [CrossRef]
  12. Steffan, T.; USDA Forest Service, Warren, PA, USA. Personal communication, 2022.
  13. Johnson, V.S.; Litvaitis, J.A.; Lee, T.D.; Frey, S.D. The role of spatial and temporal scale in colonization and spread of invasive shrubs in early successional habitats. For. Ecol. Manag. 2006, 228, 124–134. [Google Scholar] [CrossRef]
  14. Curtis, J.T. The Vegetation of Wisconsin: An Ordination of Plant Communities; University of Wisconsin Press: Madison, WI, USA, 1959. [Google Scholar]
  15. Taft, J.B.; Solecki, M.K. Vascular flora of the wetland and prairie communities of Gavin Bog and Prairie Nature Preserve, Lake County, Illinois. Rhodora 1990, 92, 142–165. [Google Scholar]
  16. Converse, C.K. Element Stewardship Abstract for Rhamnus cathartica, Rhamnus frangula (syn. Frangula alnus); The Nature Conservancy: Arlington, VA, USA, 1984. [Google Scholar]
  17. Schuster, M.J.; Wragg, P.D.; Roth, A.M.; Bockenstedt, P.; Frelich, L.E.; Reich, P.B. Using plants to control buckthorn (Rhamnus cathartica): Improved biotic resistance of forests through revegetation. Ecol. Eng. 2022, 182, 106730. [Google Scholar] [CrossRef]
  18. Heidorn, R. Vegetation management guideline: Exotic buckthorns—Common buckthorn (Rhamnus cathartica L.), glossy buckthorn (Rhamnus frangula L.), Dahurian buckthorn (Rhamnus davurica Pall.). Nat. Areas J. 1991, 11, 216–217. [Google Scholar]
  19. Schuster, M.J.; Wragg, P.D.; Roth, A.M.; Reich, P.B. No evidence of a long-lived seedbank in common buckthorn, Rhamnus cathartica L., within Minnesota deciduous forests. Biol. Invasions 2023, 25, 3043–3048. [Google Scholar] [CrossRef]
  20. Moore, A.; USDA Forest Service, Warren, PA, USA. Personal communication, 2022.
  21. Nininahazwe, F.; Varin, M.; Théau, J. Mapping common and glossy buckthorns (Frangula alnus and Rhamnus cathartica) using multi-date satellite imagery WorldView-3, GeoEye-1 and SPOT-7. Int. J. Digit. Earth 2023, 16, 31–42. [Google Scholar] [CrossRef]
  22. Isaaks, E.H.; Srivastava, R.M. Applied Geostatistics; Oxford University Press: New York, NY, USA, 1989; p. 561. [Google Scholar]
  23. Perry, J.N. Spatial analysis by distance indices. J. Anim. Ecol. 1995, 64, 303–314. [Google Scholar] [CrossRef]
  24. Perry, J.N.; Winder, L.; Holland, J.; Alston, R. Red–blue plots for detecting clusters in count data. Ecol. Lett. 1999, 2, 106–113. [Google Scholar] [CrossRef]
  25. Winder, L.; Alexander, C.; Griffiths, G.; Holland, J.; Woolley, C.; Perry, J. Twenty years and counting with SADIE: Spatial Analysis by Distance Indices software and review of its adoption and use. Rethink. Ecol. 2019, 4, 1–16. [Google Scholar] [CrossRef]
  26. SAS Institute. Using JMP Student Edition for Windows and Macintosh: The User’s Guide to Statistics with JMP Student Edition; SAS Institute: Cary, NC, USA, 2009. [Google Scholar]
  27. Grego, C.R.; Vieira, S.R.; Lourenção, A.L. Spatial distribution of Pseudaletia sequax Franclemlont in triticale under no-till management. Sci. Agric. 2006, 63, 321–327. [Google Scholar] [CrossRef]
  28. Goovaerts, P. Geostatistics for Natural Resources Evaluation; Oxford University Press: New York, NY, USA, 1997; p. 496. [Google Scholar]
  29. Ghahramani, M.; Karimzadeh, R.; Iranipour, S.; Sciarretta, A. Does harvesting affect the spatio-temporal signature of pests and natural enemies in alfalfa fields? Agronomy 2019, 9, 532. [Google Scholar] [CrossRef]
  30. Mosner, E.; Eimert, K.; Hüwe, U.; Ziegenhagen, B.; Janßen, A.; Leyer, I. Revisiting the provenance delineation of a widespread shrub, Frangula alnus—The role of spatial, temporal and environmental patterns. Tree Genet. Genom. 2017, 13, 63. [Google Scholar] [CrossRef]
  31. Labonté, J.; Drolet, G.; Sylvain, J.-D.; Thiffault, N.; Hébert, F.; Girard, F. Phenology-based mapping of an alien invasive species using time series of multispectral satellite data: A case-study with glossy buckthorn in Québec, Canada. Remote Sens. 2020, 12, 922. [Google Scholar] [CrossRef]
  32. Mays, A. Does Invasive Glossy Buckthorn (Frangula alnus) Have Allelopathic Effects on Native New England Wetland Margin Vegetation? Doctoral Dissertation, Mount Holyoke College, South Hadley, MA, USA, 2009. [Google Scholar]
  33. Godwin, H. Frangula alnus miller. J. Ecol. 1943, 31, 77–92. [Google Scholar] [CrossRef]
  34. Burnham, K.M.; Lee, T.D. Canopy gaps facilitate establishment, growth, and reproduction of invasive Frangula alnus in a Tsuga canadensis dominated forest. Biol. Invasions 2010, 12, 1509–1520. [Google Scholar] [CrossRef]
  35. Cunard, C.; Lee, T.D. Is patience a virtue? Succession, light, and the death of invasive glossy buckthorn (Frangula alnus). Biol. Invasions 2009, 11, 577–586. [Google Scholar] [CrossRef]
  36. Park, Y.L.; Krell, R.K.; Carroll, M. Theory, technology, and practice of site-specific insect pest management. J. Asia Pac. Entomol. 2007, 10, 89–101. [Google Scholar] [CrossRef]
Figure 1. The location and sample layout for each (2.5 ha) of the four research sites in the Allegheny National Forest, Pennsylvania, USA: site 1 (41.43704, −78.79137); site 2 (41.45674, −78.80247); site 3 (41.47875, −78.77288); site 4 (41.485216, −78.741848). Field surveys were conducted in 2021 and 2022. Yellow dots in each map indicate sample points where a measurement of F. alnus was taken, and white lines indicate unpaved forest roads.
Figure 1. The location and sample layout for each (2.5 ha) of the four research sites in the Allegheny National Forest, Pennsylvania, USA: site 1 (41.43704, −78.79137); site 2 (41.45674, −78.80247); site 3 (41.47875, −78.77288); site 4 (41.485216, −78.741848). Field surveys were conducted in 2021 and 2022. Yellow dots in each map indicate sample points where a measurement of F. alnus was taken, and white lines indicate unpaved forest roads.
Biology 12 01393 g001
Figure 2. Relationships of Buckthorn Index (BI) (a) and percent presence of F. alnus (b) with the distance from the forest road.
Figure 2. Relationships of Buckthorn Index (BI) (a) and percent presence of F. alnus (b) with the distance from the forest road.
Biology 12 01393 g002
Figure 3. Interpolation maps for the presence of F. alnus (%) and Buckthorn Index (BI). The maps were generated by using kriging. Gray lines indicate forest roads.
Figure 3. Interpolation maps for the presence of F. alnus (%) and Buckthorn Index (BI). The maps were generated by using kriging. Gray lines indicate forest roads.
Biology 12 01393 g003
Figure 4. Clustering maps for the presence of F. alnus (%) and Buckthorn Index (BI). Red and blue areas in cluster maps indicate patches and gaps, respectively. Blue areas represent “patches” where vi > 1.5 while red areas represent “gaps” where vj < −1.5. Gray lines indicate forest roads.
Figure 4. Clustering maps for the presence of F. alnus (%) and Buckthorn Index (BI). Red and blue areas in cluster maps indicate patches and gaps, respectively. Blue areas represent “patches” where vi > 1.5 while red areas represent “gaps” where vj < −1.5. Gray lines indicate forest roads.
Biology 12 01393 g004
Figure 5. Aerial view of linear infestations of F. alnus along the forest road. Yellow arrows indicate F. alnus. The aerial photos were taken by a drone (DJI Mavic 2 Pro, SZ DJI Technology Co., Ltd., Shenzhen, China).
Figure 5. Aerial view of linear infestations of F. alnus along the forest road. Yellow arrows indicate F. alnus. The aerial photos were taken by a drone (DJI Mavic 2 Pro, SZ DJI Technology Co., Ltd., Shenzhen, China).
Biology 12 01393 g005
Table 1. Summary of data collected for the presence and Buckthorn Index (BI) values of F. alnus.
Table 1. Summary of data collected for the presence and Buckthorn Index (BI) values of F. alnus.
SiteNumber of Sample PointsBI ± SEM% Samples with F. alnus (Presence)
145917.1 ± 0.0883.0
23947.2 ± 0.4567.8
328630.3 ± 1.5484.6
427312.8 ± 1.0664.5
Table 2. Parameters of geostatistical analysis for Buckthorn Index (BI) values and presence values. The degree of spatial dependency (DSD), expressed in percentage, was indicative of the spatial dependency between the sample points at each site. The range column expresses the maximum distance from the road in which F. alnus individuals are expected to be autocorrelated with each other. R2 and residual sums of squares (RSS) were both optimized to increase the accuracy of our calculations.
Table 2. Parameters of geostatistical analysis for Buckthorn Index (BI) values and presence values. The degree of spatial dependency (DSD), expressed in percentage, was indicative of the spatial dependency between the sample points at each site. The range column expresses the maximum distance from the road in which F. alnus individuals are expected to be autocorrelated with each other. R2 and residual sums of squares (RSS) were both optimized to increase the accuracy of our calculations.
VariableSiteModelNuggetSillPartial SillDSD (%)Range (m)R2RSS
Presence of F. alnus (%)1Exponential0.0170.1410.12487.943.70.4430.00113
2Nugget0.1820.2520.07027.78202.50.9430.00060
3Nugget0.1300.13000203.80.3090.00206
4Exponential0.1220.2460.12450.4131.80.9420.00106
Buckthorn Index (BI)1Spherical110.6369.1258.570.04155.00.993983
2Exponential61.1128.167.052.30340.80.858287
3Exponential284.0689.7405.758.8245.20.81142,461
4Exponential138.6340.4201.859.2830.90.8377860
Table 3. Parameters and associated statistics of SADIE. Ia represents the index of aggregation, and a value greater than one represents an aggregated spatial pattern. Pa indicates the significance of the respective Ia value, where a value of less than 0.05 is significant. vi and vj are the indices of clustering.
Table 3. Parameters and associated statistics of SADIE. Ia represents the index of aggregation, and a value greater than one represents an aggregated spatial pattern. Pa indicates the significance of the respective Ia value, where a value of less than 0.05 is significant. vi and vj are the indices of clustering.
VariableSiteIaPaMean vj *Mean vi **p (Mean vj)p (Mean vi)
Buckthorn Index (BI)15.180.0003−5.0654.579<0.0001<0.0001
23.490.0003−3.2132.912<0.0001<0.0001
31.740.0003−1.6491.7460.00150.0005
42.640.0003−2.8572.639<0.0001<0.0001
Presence (%)12.400.0003−2.4052.422<0.0001<0.0001
24.130.0003−4.0624.197<0.0001<0.0001
31.910.0003−1.8321.9440.0005<0.0001
41.830.0015−1.9041.7400.00080.0023
* vi = index of clustering for points in which the BI or presence value is above average for that site. ** vj = index of clustering for points in which the BI or presence value is below average for that site.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Greenleaf, J.; Karimzadeh, R.; Park, Y.-L. Spatial Patterns of Frangula alnus (Rosales: Rhamnaceae): Implications for Invasive Plant Management. Biology 2023, 12, 1393. https://0-doi-org.brum.beds.ac.uk/10.3390/biology12111393

AMA Style

Greenleaf J, Karimzadeh R, Park Y-L. Spatial Patterns of Frangula alnus (Rosales: Rhamnaceae): Implications for Invasive Plant Management. Biology. 2023; 12(11):1393. https://0-doi-org.brum.beds.ac.uk/10.3390/biology12111393

Chicago/Turabian Style

Greenleaf, Jennifer, Roghaiyeh Karimzadeh, and Yong-Lak Park. 2023. "Spatial Patterns of Frangula alnus (Rosales: Rhamnaceae): Implications for Invasive Plant Management" Biology 12, no. 11: 1393. https://0-doi-org.brum.beds.ac.uk/10.3390/biology12111393

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop