The Louisiana Amphibian Monitoring Program from 1997 to 2017: Results, analyses, and lessons learned

To determine trends in either frog distribution or abundance in the State of Louisiana, we reviewed and analyzed frog call data from the Louisiana Amphibian Monitoring Program (LAMP). The data were collected between 1997 and 2017 using North American Amphibian Monitoring Program protocols. Louisiana was divided into three survey regions for administration and analysis: the Florida Parishes, and 2 areas west of the Florida parishes called North and South. Fifty-four routes were surveyed with over 12,792 stops and 1,066 hours of observation. Observers heard 26 species of the 31 species reported to be in Louisiana. Three of the species not heard were natives with ranges that did not overlap with survey routes. The other two species were introduced species, the Rio Grande Chirping Frog (Eleutherodactylus cystignathoides) and the Cuban Treefrog (Osteopilus septentrionalis). Both seem to be limited to urban areas with little to no route coverage. The 15 most commonly occurring species were examined in detail using the percentage of stops at which they observed along a given survey and their call indices. Most species exhibited a multimodal, concave, or convex pattern of abundance over a 15-year period. Among LAMP survey regions, none of the species had synchronous population trends. Only one group of species, winter callers, regularly co-occur. Based on the species lists, the North region could be seen as a subset of the South. However, based on relative abundance, the North was more similar to Florida parishes for both the winter and summer survey runs. Our analyses demonstrate that long-term monitoring (10 years or more) may be necessary to determine population and occupancy trends, and that frog species may have different local demographic patterns across large geographic areas.


Introduction
In the 1990's, concerns over perceived declines in North American amphibian populations led the members of the Declining Amphibian Population Task Force to recommend the implementation of a statistically defensible amphibian monitoring program [1]. In 1997 the North American Amphibian Monitoring Program (NAAMP) was established to answer the question, 'Are these perceived declines real?" To address this question, NAAMP utilized a network of largely volunteer observers who monitored calling frogs along roadside routes, similar to the Breeding Bird Survey [2]. From 1997 to 2015 the U.S. Geological Survey (USGS) coordinated NAAMP. The USGS developed the monitoring protocol, provided the states with randomized route starting points, and stored data provided by the states. The states recruited and trained agency personnel and volunteers to conduct the surveys. In 2015, the USGS terminated the NAAMP. The Louisiana Amphibian Monitoring Program (LAMP) is Louisiana's state 'chapter' of NAAMP. The purpose of the LAMP is to determine if there are changes in frog distribution and abundance in the state of Louisiana over time. The first LAMP surveys were conducted in 1997, and the program is still active as of 2021. Fifty-nine routes were set up across the state, 54 of which were surveyed one or more times from 1997 through 2017. Using data collected with NAAMP protocols we want to determine if there were changes in frog distribution and abundance in Louisiana over the observation period. We had three broad questions: (1) 'Can we use call data to detect trends in frog populations?' (2) 'Are there trends in frog populations?' and (3) 'If there are trends, are these associated with frog communities as a whole or with individual species, or both? ' We wanted to look at the above question at two hierarchical levels: for frog communities as a whole and for individual species. Finally, we wanted to use LAMP data to describe Louisiana's anuran communities. Our community-and species-level questions were as follows: Community Level Questions: 1. Are there significant changes in species richness over time for a given route or in a given region?
2. Are there significant changes in species call index (as a proxy for abundance) over time for a given route or in a given region? That is, independent of species richness, are more or fewer frogs calling?
3. Do some species co-occur and can they be considered as a community or as indicators of a particular habitat?

Do the different LAMP regions have different community compositions?
Species Level Questions: 1. Which species were detected during the surveys, where, and when?
2. Which species were the most observed and which have the highest abundance (call index) when observed?
3. Did any species show a change in frequency of observations or abundances (call index) for a given region, route, or time of year?

Methods
LAMP data used in this study can be downloaded from the USGS's ScienceBase data server [3].

Survey routes
Thirty-seven (37) of the 64 parishes (i.e., 'counties') in Louisiana have one or more survey routes. Starting with locations and directions randomly selected by the NAAMP program, LAMP volunteers set up survey routes with 10 stops each. Once established, survey routes and stops can only be changed for reasons of safety or compliance with protocols. S1 Table lists LAMP Route names, NAAMP numbers, and the latitude and longitude of first stop. Fig 1 is a state map showing LAMP route starting points. The first surveys were conducted in 1997. The total number of stops observed varied between routes for two reasons. Not all routes were surveyed every year, or completely surveyed in a given year. A subset of 54 of the 59 established routes were deemed to have sufficient data for analysis ( Table 1). The state was divided into three regions for administrative and analysis purposes based on what was known about frog and toad distributions and regional phenology at the time the survey was initiated. The regions are the Florida Parishes, North and South (Fig 1). The Florida Parishes region (henceforth, 'Florida region') is north and east of the Mississippi River and north of Lake Pontchartrain and Lake Maurepas. The area of the state west of the Mississippi River is divided between North and South at latitude 31˚41' N. This line reflects a geographical break in the randomly assigned route locations between north and south Louisiana, and the 16˚C average annual temperature line [4], where north of that line the winters are colder, longer, and frogs start calling later in the year. The 40-day survey windows (or 'runs') for the South and Florida regions are the same; the survey windows for the North begin about 30 days later than those in the South due to their colder climate ( Table 2).
The LAMP followed the NAAMP survey protocols [5]. Each survey was initiated at least ½ hour after sunset and before midnight. The surveys were run three times a year during  Table 2). The Florida region is comprised of the 'Florida Parishes', locations east or north of the Mississippi River. The areas west of the Mississippi River are divided at approximately 31˚41'N latitude. See S1 Table for  The base map is from ESRI online (for a description of the copyright information see: http://goto.arcgisonline.com/maps/World_Topo_Map), and is in turn based on USGS topo data, and is not copyrighted. The base map is a state topo map with parish (i.e., county) boundaries and larger geographic features labeled. ArcInfo was used to superimpose on the public domain base map the locations of NAAMP route survey routes and then to generate a TIFF file of the new map. The LAMP regional boundary lines were overlaid on the TIFF file using MS PowerPoint software and then saved as a new TIFF file. The coordinates used are listed in S1 Table. https://doi.org/10.1371/journal.pone.0257869.t001

PLOS ONE
The Louisiana Amphibian Monitoring Program: 1997Program: -2017 windows meant to capture 'winter', 'spring' or 'summer' calling frog species. The observers did not need permits or approvals to conduct their surveys. All routes and stops were on public roads along right-of-ways, and didn't require landowner permissions. No permits were required by either state or federal agencies and no endangered or protected species were interacted with. Volunteers only passively listened for frogs and did not come into contact with them nor call to them in order to get a response. These observation windows were called Runs '1', '2', and '3' respectively. Each stop was considered a separate observation and ideally there were a total of 30 observations per route, per year. The frog call observers listened for 5-minutes at each stop, logged frog species heard calling, and assessed an index of calling activity for each species. Using the stratified by habitat protocol [5] stops were spaced a minimum of 0.8 km apart so that the same individuals were not heard at neighboring stops. As a quality control measure, observers were required to take an online test on their ability to identify frog calls and rank calling activity [6]. In addition to the frog calling data, observers reported the time of day the survey was conducted and recent weather conditions. Data sheets from survey runs were forwarded to the Louisiana state coordinator (J. Boundy) who performed quality control and entered the data into the national database.

Route data selection
The routes used depended on the analysis ( Table 1). All stops with one or more species calling were used for analyses of species richness, phenology, and average call intensity (54 routes met these conditions). All stops with two or more species calling were used for co-occurrence analyses (45 routes met these conditions). For trend analysis using general linear models, routes had to have data sets from 8 or more different years, and the surveys needed to span 20 or more years (21 routes met these conditions). For General Additive Model (GAM) analysis route data had to have data sets from 8 or more different years, surveys spanning 20 or more years, and no gaps in the data of more than 5 years (19 routes met these conditions). The data were examined at different hierarchical levels in both space and time. In all cases time (years) was the independent variable. When all three runs per year were combined data are referred to as "route" or "region". "Region" data included data from all routes in a given region. When data were examined by using all the routes in a given region for a given run, data are referred to as "region-run" data (3 regions, 3 runs per year, 9 possible combinations). When data were examined looking at individual routes separated by run, it is referred to as "route-run" data (54 routes, 3 runs, 162 possible combinations). Individual species trends were examined by both region-run and route-run (171 possible combinations).

Species selected for analysis
All species and all observations were used to determine when species called and where they were observed. A species required a minimum of 200 observations for trend analyses. For a

PLOS ONE
The Louisiana Amphibian Monitoring Program: 1997-2017 species route-run trend analysis the species needed to be detected over 8 or more years on the route-run.

Route and region: Species richness and average call index
Species richness for each stop was the number of species heard. Species richness for a routerun was the number of different species observed over all stops during a run. The average species richness for a run-route was the average species richness of all stops along the route, including stops where no frogs were observed. We used a linear model to estimate the change (slope) in route species richness over the 1997-2017 observation period for all the routes by run for each region. The average call index (ACI) for a route-run was the average of all call indices of all species observed. Stops where no frogs were heard calling were not used in calculating the average. Species richness and ACIs were calculated for both the regions and routes, and by run. We used a general linear model to estimate the change (slope) in ACI over the 1997-2017 observation period.

Regional assessment of species co-occurrence
A Two-Way Species Indicator Analysis (TWINSPAN) in PC Ord software [7] was used to determine if groups of species formed 'communities' of co-associations [8,9]. In this analysis each stop was treated as a separate observation and only observations with two or more species were used. The species used for this analysis also needed a minimum of 10 observations from routes with 5 or more years of observation for that species.

Relative abundance analysis
A Pielou's Evenness Index (J) is a measure that relates species richness (i.e., the number of other categories) to their relative abundances [10]. J is the ratio of the measured Shannon diversity index (Eq 1) to perfectly even abundances (Eq 2). If J equals 1, all species are present in equal abundance. The smaller J, the more skewed the observations.
for i = 1 to R; where p i is the relative abundance of species or categories 'i', and R is the total number of species or categories where H max = ln(R) J was calculated for the species observations statewide and by region-run. J was also used to compare the evenness of the sampling effort between regions using the routes as categories and number of stops as observations.

Species: Percent observation and average call index
The data were analyzed to determine if there were trends in: how often species were observed; species richness, or the species call abundance index. All runs in which a given species was observed were used to calculate that species' percent observation for that run. For example, if a species was observed on three different stops percent observation for that run was 0.3. This analysis was done for regions by run, and for routes by run.
A species average call index (sACI) was the average of all non-zero call indices of that species on a route-run or region-run. A routes-run ACI was the average for all species calling on a given region and run. A region-run ACI was the average for all species calling on a given route and run. We used PROC REG in SAS 9.4 [11] for linear model analysis to look for trends in route-run and region-run ACIs. We use PROC GAM in SAS to look for species sACI trends by route-run and region-run.

Generalized Additive Model (GAM) analysis
When we ran simple linear regressions for route-run and region-run species richness and ACI, the assumptions of normality and homogeneity were violated. Since those could not be fixed with a transformation, we ran a generalized additive model (GAM) using the generalized cross validation approach, whereby the estimated parameters are chosen via a generalized cross validation [12]. GAMs are used when there is an expectation that behavior along the dependent variable may not be linear or may change over time. The GAM procedure fits generalized additive models as defined by Hastie and Tibshirani [13]. PROC GAM in SAS 9.4 was used using non-parametric regression and smoothing [11].
We ran GAM analysis for the species with 200 or more observations (15). An individual route-run used in this analysis had to have at least 8 stops visited. A species needed a minimum of 8 or more route-run observations for GAM analysis of percent observations. Route-run combinations found to have significant trends then had their GAMs classified as one of the following trend patterns: increasing, decreasing, concave, convex, or multimodal.

Route statistics
The number of routes surveyed, and the number of runs completed varied from year to year. Survey data for a given route might be incomplete with missed stops, runs, or even consecutive years missing. Additionally, even when a route-run combination was completed, there might be one or more stops where no frogs were heard. Finally, there was with a significant decline in the number of routes surveyed starting in 2016 (Fig 2).
Between 1997 and 2017, 75 observers surveyed 54 routes and 12,792 stops. Sampling was uneven between the regions (Table 1) with the South having the largest number of routes (37) and stops observed (8,772), followed by North (12 Routes/ 2,540 stops) and Florida (5 routes/ 1,480 stops). Between 1999 and 2015 the number of stops varied between 500 and 700 per year. Species richness at a given stop varied between 0 (i.e., no species heard at the stop) to 11 (Fig 3).
If we treat each route as a category (i.e., as a 'species') and the number of stops observed as abundances, we can use J to compare the evenness of sampling between the different LAMP regions. We can then ask the question "Were some regions more reliant on heavily sampled routes for their results than other regions?" The J for the three regions using the number of stops visited per route was 0.961, 0.844 and 0.930 for Florida, North, and South respectively. Thus, while the number of routes surveyed was different in each region; the sampling effort for the routes within a region was comparable between regions. and 3.0. Statewide, the frequency a species was observed was not a good predictor of average calling intensity (Fig 4). The taxonomy used here follows the Integrated Taxonomic Information System [14].
The frequency a given species was observed when plotted against its rank abundance was best modeled as an exponential function (Fig 5). J for the state overall was 0.30. Species could be grouped by when they were heard calling as either winter callers, spring callers, or summer callers (Table 4). Cajun Chorus Frogs (Pseudacris fouquettei), Spring Peepers (P. crucifer), L. sphenocephalus, Crawfish Frogs (L. areolatus) and Pickerel Frogs (L palustris) were winter (e.g., Run 1) callers. Two toad species were primarily spring (e.g., Run 2) callers: American Toads (Anaxyrus americanus) and Oak Toads (Anaxyrus quercicus). Southern (Acris gryllus) and Northern Cricket Frogs (Acris crepitans) were almost evenly divided between spring and summer, and the rest of the species were observed calling 62% or more of the time during the summer (e.g., Run 3). One summer species, S. holbrookii, was never heard during the spring Run.
All species that called in winter also called in spring or summer, while 2/3 of species classified as Spring or Summer callers did not call in winter. Because of the above observations in addition to statewide trend analyses, some trends were analyzed by region, or region and run. Since there were only two spring calling species, we combined the spring and summer runs for some analyses (see below).
There were 5 species reported to be in Louisiana that were not observed: 3 native species, the Dusky Gopher Frog (L. sevosa) Ornate Chorus Frog (P. ornata) and Strecker's Chorus Frog (P. streckeri); and 2 exotic species, the Rio Grande Chirping Frog (Eleutherodactylus cystignathoides), and the Cuban Treefrog (Osteopilus septentrionalis. Distribution maps [15,16] indicate that the documented locations of the L. sevosa, P. ornata and P. streckeri were not in the vicinity of LAMP routes. Also, L. sevosa and P. ornata are considered extirpated from the state, not having been seen since the 1960's [16]. E. cystignathoides has been documented in urban areas around Baton Rouge, Lafayette, Lake Charles, Alexandria, and Shreveport [17] since the early 1990's (J. Boundy, pers. obs.). O. septentrionalis has been documented in the cities of Baton Rouge, Lafayette and New Orleans since 2018 [18]. Neither species was observed on any survey route as of 2017.
Region: Florida, South, and North had 23, 22, and 17 species in their surveys respectively ( Table 5). The most frequently observed species in the South was L. sphenocephalus and for Florida and North, P. crucifer. Both are winter calling species.
When we examined abundance by region and run ( Table 6) the most commonly observed species for the regions as a whole were winter callers. For the spring and summer runs, the most common species were Green Treefrogs (Dryophytes cinereus) in the Florida region, Cope's Gray Treefrog (D. chrysoscelis) in the North region, and A. crepitans in the South region. The Pielou evenness index, J, for species observations was more even for Florida (0.8255) than for North (0.7935) and South (0.7940). Winter runs were consistently more uneven, meaning their observations were dominated by fewer species.

Frequency observed and call intensity
Region analysis. The ACI varied year to year ( Table 7). The only region-run combination with an R-squared greater than 0.3 was Florida-Run 2. It appeared to have a sustained decrease in average call intensity after 2004 (Fig 6).  Table 8). The Rsquared values ranged from 0.2357 to 0.8445.
Species richness. Regional trends. North-Run 1 and the South-Run 3 had significant positive changes in mean species richness (Table 9). No other region-run combination had significant changes.
Route-run trends. Twenty-one routes had sufficient data to examine species richness. Using a p-value of 0.05 or less we found 16 route-runs with significance (Table 10). Five (5) routerun combinations had significant decreases in species richness and 11 had increases. Interestingly, in one case (Rayne), there was a decrease in Run 2 but an increase in Run 3.

Two Way Indicator Species Analysis (TWINSPAN)
TWINSPAN was conducted separately for each LAMP region (Fig 7). The first group to separate out for all three regions were winter calling species (Run 1). Beyond the winter calling species no groupings were consistent between regions.

PLOS ONE
The Louisiana Amphibian Monitoring Program: 1997-2017

PLOS ONE
The Louisiana Amphibian Monitoring Program: 1997-2017

Generalized Additive Model (GAM) analysis
We examined the 15 most common species along their possible route-run combinations over time to determine if they had changes in percent observation or their call index (as a proxy for abundance). Plots of the 12 species with 10 or more years' worth of data and significant GAMs (p � 0.05) were examined for model type (S1 File). There were 39 significant GAM models classified: multimodal = 32, convex = 6, concave = 9, and linearly increasing = 2 (Table 11). Seven (7) models showed an overall increase, and 1 an overall decrease over the observation period. There were 9 regional models, and the remaining 30 models were route specific. At 12 each, North and Florida had relatively more models than South region's 15 models. This may be due to the fact that many South routes had less than 10 years' worth of survey data or breaks of 5 or more years between surveys, whereas the routes in the other two regions tended to have more complete data sets. At 10, the Bronze Frog (L. clamitans) had the greatest number of significant GAM models.

Community level questions
As noted above, because different species tended to call during different times of year, the seasonal survey runs were examined independently.

PLOS ONE
The Louisiana Amphibian Monitoring Program: 1997Program: -2017 When regressing the number of observations against rank abundance, the exponential model over predicted expected observations of the most common species, however its performance is much better than the alternatives models (logarithmic, linear, and polynomial), which do much worse with respect to predicting the abundance of less common species.
(1) Were there significant changes in species richness over time on a given route or region? No, with the following exceptions. By Region: Only two region-run combinations showed a statistically significant (0.05 alpha level) change in species richness over time, North-Winter and South-Summer (Table 9). Both were small and positive. By Route: The majority of routerun combinations in all regions had no significant change in species richness. Of the 16 that did, 5 routes-runs showed declines in species and 11 had increases (Table 10).
(2) Are there significant changes in species call index (as a proxy for community abundance) over time on a given route or region?
No. With the few exceptions noted in the Results section, there were no changes in ACI for the majority of region-run or route combinations examined.
(3) Do some species co-occur together? Yes, for winter callers. No for spring and summer callers. TWINSPAN showed that winter calling species formed a group in all regions of the state. An additional species, L. areolatus,

PLOS ONE
The Louisiana Amphibian Monitoring Program: 1997-2017 was present in the South region. Patterns with respect to other species were not apparent in TWINSPAN. Part of this may be due to sampling bias. More than twice as many routes in the South than for Florida and North combined were surveyed and when statewide data were pooled the groupings largely mirrored what the South TWINSPAN showed and did not provide insight for the state as a whole. Therefore, we only present the separate regional   Table 7. Only Florida-Run 2 had a significant slope. https://doi.org/10.1371/journal.pone.0257869.g006 TWINSPAN results (Fig 7). However, when the TWINSPAN results for different regions were compared, no other groups besides winter callers were consistent between regions. This may be in part due to the way differences in species composition and abundances lined up between the regions (see below). (4) Are there differences between the LAMP regions in community composition? Yes. Run 1 had lower J indices than the scores for Runs 2 and 3 for all regions. Indeed, there was a core of 3 winter calling species (L. sphenocephalus, P. fouquettei, and P. crucifer) that made up more than 95% of the observations in all three regions.
Differences in the number of routes sampled and the number of runs completed make direct quantitative comparisons difficult. However, we can make a few observations. Florida had 9 species that were not found in either the South or North regions. The South and North It should be noted that the species that did set the regions apart were relatively rare ( Table 6). The uniquely Florida species only accounted for 2.16% of the relative abundances for that region, and the species that were not found in Florida only accounted for 0.48% and 0.05% of the relative abundances in South and North respectively.
When we compare the most abundant species however, a different picture emerges (Table 6). Of the 6 most common species in each region, North shares 5 of its 6 with Florida, while South shares only 3 each with Florida and North. When we look at the highest ranked species (those that make up 59% or more of the observations), all 4 of North's species are shared with Florida, but only 2 of its 4 species are shared with South. So, while comparisons of the species list would group South and North, comparisons of actual observations would group North and Florida. A comparison of the region-runs finds the same pattern, whereby the lists of the most observed species in North and Florida, by rank or relative abundance, are more similar to each other than either is to South. This is true for both Run 1 and Runs 2 and 3 combined.
In summary, the regions and runs do differ. Winter calling species held up as a community in and of themselves both in composition and relative abundances. If we were simply comparing species lists, the North would be seen as a subset of the South. When comparing abundances however, North and Florida are more similar to each other than either is to South both regionally and seasonally.

Species level questions
1. Which species were detected during the surveys, when and where?
Twenty-six species were detected (Table 4). Five were classified as winter callers, 2 spring callers, 2 were equally divided between spring and summer, and the rest (17) were summer callers. No species were found outside of the regions for which they were already known. Three native species were not observed, but the survey routes did not overlap with their known distributions and two are considered extirpated since the 1960's. Two recently introduced species, O. septentrionalis and E. cystignathoides, were not observed along any survey route. This may be due to fact that their distributions may be in areas without survey routes.
2. Did any species show a change in the frequency of observations or abundances for a given region, route or run? Yes they did, but no consistent pattern emerged for any given species across all routes in a given region or statewide. The sACI for a species along a given route-run was not related to how often a species was observed along the same route (Fig 4). The rarest species in the LAMP data, S. holbrookii, had the highest sACI. Eleven (11) of the 15 most common species had at least one GAM with significant changes in how often they were observed along a route over time (Table 11). All but one model had at least one change in slope sign from positive to negative or negative to positive. Seven (7) models showed an overall increase and 1 a decrease in percent observation over time. The majority of GAM models of percentage of stops where a species was observed on a route-run were not significant, and those that were in general had multimodal patterns.

What are the most common species?
Lithobates sphenocephalus and P. crucifer were the most commonly observed species. The most observed species statewide was L. sphenocephalus. When examined by region the most common species in the South was L. sphenocephalus and the most common species in Florida and North were P. crucifer. Both are winter callers. The species with the highest sACI, S. holbrookii, was also the species with the fewest observations suggesting that for some species, the call index may not be a good measure of abundance.

Previous studies
Villena et al. [19] used a portion of the LAMP data as part of a larger study of frog population trends in the Southeast. Their study used fewer species and covered fewer years than the data used for this study. They used occupancy modeling to document changes in their occupancy.
In Louisiana only two toad species, Fowler's Toad (Anaxyrus fowleri) and Southern Toad (Anaxyrus terrestris) were found to have significant changes in their occupancy slopes (positive for Fowler's, negative for Southern). The other species were found to have slopes with confidence intervals that overlapped with 0 (no change). In our analyses, we found A. fowleri to have a multimodal abundance pattern. We did not conduct an analysis of the A. terrestris data because there were not enough observations (N = 45) to meet our inclusion criteria (>200 observations). In this study species patterns varied between regions. A. fowleri made up less than 4% of the observations for Florida and South but 11.86% of the observations for North. The only significant GAM for A. fowleri was for Rocky Branch (North). It included only 6 years and was multimodal.
In this study, 14 species were found to have significant GAM models with respect to changes in observations for 71 different route-run combinations ( Table 8). The relevance of this is the majority of species had multiple peaks or dips in their populations over the study period and that short observation windows of less than 10 years are probably inadequate to document long-term trends.
A study in the Atchafalaya Basin of Louisiana [20] found that all 12 species detected in the study experienced declines in occupancy during their 5-year observation period (2002)(2003)(2004)(2005)(2006). Our GAM analysis of percent observations along a route (Table 11) suggests that their observations may not be generalizable for the rest of Louisiana, and that their observation window may have been too short to assess long-term species status.

Conclusions
Taken together, our results lead to three broad conclusions. (1) Can call monitoring detect trends in frog populations? Yes. We were able to follow trends in abundance and distribution for some routes for 20 years. During this time populations often went through one or more cycles of changes in abundance. The length of these cycles demonstrated the importance of long-term (10+ years) monitoring. Observation periods of less than 10 years, or with long multiyear gaps, are likely to be too short and/or incomplete to determine if a species is experiencing long-term declines or increases in abundance or occupancy. (2) Are frog abundances in Louisiana declining? Perhaps in some locations, but overall no. In some cases abundances increased. Our analyses detected no long-term trends in either the percent observation along a given route or region, or ACI along a given route or region. We did not detect long-term species richness trends in any of the three LAMP survey regions. (3) Are their regional trends in frog abundance, and are frog populations in a given region synchronous? With a couple exceptions, no to both questions. The routes were, in general, not synchronized with respect to their patterns of species richness, abundance (call index), and species occupancy (percent presence) and different regions had different dominant species.
New threats to frog species have emerged over the last 10 years. O. septentrionalis, has recently invaded the cities of New Orleans and Lafayette, and have been expanding their range from urban areas [21]. They are predators on other species of frogs [21]. A second potential emerging threat are Giant Apple Snails (Pomacea maculata), which are potential frog egg predators [22]. Pomacea maculata have been expanding their range throughout the southern parishes of Louisiana [23] and are often found in many of the same habitats that frogs occupy. Continued monitoring my help inform whether these or other unknown threats are impacting frog populations. Monitoring provides an important baseline for future assessments. The LAMP data provides a baseline that can be used to determine if there are significant changes in frog populations or distributions caused by these or other threats.
Supporting information S1 Table. LAMP routes. Route name, NAAMP number, latitude and longitude of first stop. (DOCX) S1 File. S1-S21 Figs. GAM plots of the percentage of stops a species was observed calling along a given route-run versus year. (ZIP)