Introduction

The Mongolian racerunner, Eremias argus, is a small lizard found throughout China, Mongolia, and Korean peninsula (Do et al. 2017; Orlova et al. 2019). In contrast to other Korean reptiles, E. argus only inhabits certain environments, such as coastal dunes and sandy banks along rivers and sea (Kim et al. 2011a). Given the limited habitat of E. argus, habitat disturbance can be fatal for members of this species. For example, species richness, species diversity, functional richness, functional dispersion, and the number of functional groups of lizards in the disturbed sites are lower than those in undisturbed sites (Berriozabal-Islas et al. 2017). This lizard is currently found in limited areas and there is an urgent need for its protection and management (Song et al. 2010, 2013). In South Korea, E. argus is an endangered species owing to its limited distribution and decreasing population size, and therefore, it is managed (Kim et al. 2011b). Research on E. argus is centered on population restoration by focusing on reproduction, home ranges, and food resources (Kim et al. 2010; Song et al. 2010; Jeong and Song 2010; Kim et al. 2012). Although there are several studies on restoring the population of E. argus, limited data are available on protecting and managing the habitat of this lizard in South Korea.

The Baramarye Coast includes one of the largest coastal sand dunes in Taean-gun, South Korea (Kim et al. 2011a; Park et al. 2013). Many of the coastal plants found in the coastal dunes act as both refuges and habitat for endangered E. argus (Kim et al. 2011a; Park et al. 2013). In 2008, a study on the distribution and population size of E. argus in Baramarye Coast was conducted (Song et al. 2013). The Taeanhaean National Park management designated this site as a natural protected area in 2009. The coastal dune restoration project removed concrete structures and sand trap installations from the area in 2013 and 2014 (Yu et al. 2016). As a result, coastal grasslands, which are the primary habitats of E. argus, increased by approximately 190%, from 8314 m2 in 2008 to 15,829 m2 in 2016 (Song et al. 2016). Furthermore, previous studies predicted continuous sand sedimentation in Baramarye Coast and Halmi Island due to wind direction and speed data (Park et al. 2013; Yu et al. 2016). The environmental conditions of Baramarye Coast are changing, and the subsequent effects on the different organisms present in the area, including E. argus, are unknown.

In this study, we sought to evaluate the hypothesis that the restoration of E. argus habitat, namely coastal dunes and grasslands, will affect the population expansion of Mongolian racerunner. The number of E. argus present in 2008 (Song et al. 2013) and 2016 was compared to evaluate the effect of habitat expansion on the lizard’s population size and determine the effects of restoration efforts. We predicted that the expansion of E. argus habitat, specifically the coastal grassland area, would lead to an increase in habitable area and consequently increase the size of E. argus population. Second, we analyzed the changes in E. argus habitat use and distribution due to habitat expansion. We predicted that the observed changes in habitat would increase the distribution of E. argus. This study not only identifies the effect of in situ conservation efforts on E. argus, but also provides important evidence for the future protection and management of organisms that are found in restricted habitats.

Materials & methods

Research area

Baramarye special protection zone (36.412417°N, 126.384440°E, WGS84, 3 m a.s.l.) is located at the southern end of Anmyeon-eup, Taean-gun. Baramarye Coast is characterized by well-developed sandy beaches and the presence of grassland patches. The lowest height from the sea level of the coast is normally 3–5 m, and the maximum height is 17 m (Park et al. 2013). A previous study reported 49 native and 8 non-native plants distributed in the Baramarye special protection zone, including the following communities: 5 arbors, 10 halophytes and sand dune plants, and 4 mixed plants (Park et al. 2013). Most plant communities form small patches and are usually distributed between 100 and 150 m from the coastline (Park et al. 2013). This zone is a part of the Taeanhaean National Park and therefore is well managed and preserved.

Field research

Considering the main activity season of E. argus, we conducted the field survey between April and June 2016 following the methodology of Song et al. (2013). To reduce the methodological difference between two studies in 2008 and 2016, we used exactly the same methods as in the previous study (Song et al. 2013): three researchers, capturing methods, survey time, and individual labeling. Usually, lizards immediately move as humans approach. We marked the initial discovery points of an individual animal with an iron stick to record the lizard’s original point. We used a portable GPS (Oregon 650, Garmin, USA) to obtain the coordinates of each location. To prevent repeated counting, all lizards were captured and kept in an ecorium until the completion of the field survey (Fig. 1).

Fig. 1
figure 1

Endangered Mongolian racerunner, Eremias argus, in Baramarye coast, South Korea

Individual tagging

Repeated capture of the same individual during a day could be a problem in estimating the population size; therefore, we used the tagging method. Typically, two tagging methods are used for the study of small lizards: inserting a passive integrated transponder microchip or partially clipping either a scale or toe (Langkilde and Shine 2006). Previous studies showed that toe clipping does not substantially affect lizard mobility, and thus we employed this method in the present study (Dodd 1993; Borges-Landáez and Shine 2003; Langkilde and Shine 2006). In a previous study, researchers clipped many toes for numbering (Song et al. 2016). In this study, however, we clipped only one joint of the longest fourth toe of the left foot to identify recaptured individuals. We released individual animals to the same location of their capture after tagging. We followed the guidelines for field research of live reptiles by Beaupre et al. (2004).

Population size estimation

We compared the total number of individuals, number of individuals per sampling month, and number of recaptured individuals in 2016 with those of a previous study (Song et al. 2013). To estimate the population size of E. argus in Baramarye Coast in 2016, we used the Peterson method (Begon 1979), and compared the result with those of the previous study (Song et al. 2013). The Peterson method is the most effective and simple method to estimate population size from two studies (Begon 1979). A few experimental assumptions are required for the Peterson method: 1) the population is closed, 2) all animals have the same chance of being caught in the first sample, 3) marking of individual animals does not affect their catchability, 4) animals do not lose their marks between the sampling periods, and 5) all marks are reported on discovery in the second sample (Begon 1979; Krebs 1998). We considered that the E. argus population was closed given that their habitat was limited to coastal dunes, and we found no E. argus individuals in other nearby regions. The number of researchers, duration of the investigation period, and sampling effort were the same in 2008 and 2016. We estimated the size of the E. argus population using the Peterson method with the following equations:

$$ N=\frac{M\left(C+1\right)}{\left(R+1\right)}\left(\mathrm{and}\right) $$
$$ Variance\ (N)=\frac{M^2\left(C+1\right)\left(C-R\right)}{{\left(R+1\right)}^2\left(R+2\right)}, $$

where, N refers to group size, M refers to the number of marked individuals in the first sample, R refers to the number of recaptured individuals in the second sample, and C refers to total number of individuals identified in the second sample. For population size estimation, we calculated N ± 1.96√(N) to obtain the 95% confidence interval.

Changes in the distribution of E. argus

To estimate the home range of E. argus in 2008 and 2016, we used the Kernel density (KD) method. It is well known that distinct boundaries such as cliffs, river, and ponds considerably affect the results of the Kernel method (Getz and Wilmers 2004; Getz et al. 2007). In our study site, there are no physical barriers to prevent the movement of lizards. Based on locations in which we captured E. argus in 2008 and 2016, the characteristics of their distribution and changes in their distribution were compared with the changes in their size of habitat area. For the KD analysis, the following options were set: area unit, square map unit; output value, density; and method, planar. We used ArcGIS ver 10.6 (ESRI, USA), and represented the result as 95% KD and 50% KD.

Results

During the field survey in the previous study (Song et al. 2013), 152 E. argus individuals were recorded (Fig. 2). In our study, we identified 147 E. argus individuals in April and 197 individuals in June, 2016 (Table 1). In the previous study, 15.1% (23 individuals) of E. argus were recaptured, while we recaptured 13.6% individuals (20 individuals; Table 1). We found no individual animals that were marked in 2008 during this study.

Fig. 2
figure 2

The location of Baramarye coast and the Halmi island (A). The spots (red: 2008, blue: 2016) indicated the points of discovery the Mongolian racerunners

Table 1 Number of captured Eremias argus from field surveys in Baramarye coast, South Korea. Recapture = recaptured lizards from second survey / lizards from first survey. *: Song et al. 2010

Population size estimation using the Peterson’s method predicted that the E. argus population in Baramarye Coast comprised 352–787 individuals in 2016. This estimation represents an increase by 126–137 (21.1–55.7%) lizards compared with that in the previous study (Song et al. 2013).

The distribution of the E. argus population was highly limited to dune grasslands near Baramarye Coast and northern Halmi Island in the 2008 study (Fig. 2). After habitat restoration efforts, the distribution of E. argus in the research site increased, and we found individuals in a larger area of dune grasslands near Baramarye Coast, northern Halmi Island, and the southern coast of Halmi Island in 2016 (Fig. 3).

Fig. 3
figure 3

The result of Kernel density analysis of the distribution of Eremias argus in Baramarye coast (a) in 2008, (b) in 2016. The black polygon indicates the area of significantly extended dune grassland after habitat restoration. Red and yellow polygons present 95% KDE and 50% KDE, respectively

The home range of E. argus was 12,092 m2 for 95% KD and 1329 m2 for 50% KD in 2008 (Fig. 3a). In 2016, after the restoration of coastal sand dune, the home range of 95% and 50% KD increased by 58,000 m2 (4.8-fold) and 4737 m2 (3.6-fold) compared with that in 2008 (Fig. 3b).

Discussion

In this study, we examined the effects of coastal sand dune restoration on the distribution and population of E. argus. We found evidence to support our hypothesis that the expansion of coastal dunes and grasslands can lead to an increase in the size of the E. argus population. Habitat expansion and restoration, such as the removal of threats within habitats, can lead to an increase in population (Towns 1994; Pike et al. 2011). In particular, species found only in certain habitats may be greatly affected by habitat restoration efforts (Amo et al. 2007). As a result of a habitat restoration project in Baramarye Coast that was carried out in 2013, the dunes and dune grasslands expanded, and may have resulted in an expansion in the distribution of E. argus. In other words, the changes in population size that coincide with habitat restoration, as identified in this study, present a meaningful case study for the effect of habitat restoration efforts on the presence of organisms.

Recently, the coastal region has high pressure from human actions such as tourism, building resorts, and public beach (Bramwell 2004; Nordstrom 2004). This pressure leads to overgrazing, deforestation, and exploitation and finally destroys the habitat of endemic species in coastal regions. However, the population of E. argus has been decreasing, because their habitats are mainly in coastal dunes vulnerable to human pressure (Song et al. 2010; Kim et al. 2011a; NIBR 2019). However, our results showed that the restoration and expansion of coastal sand dunes could recover the population of E. argus and even expand their home range. In other words, habitat restoration is an effective way of conserving species that inhabit a specific environment such as coastal dunes.

The introduction of invasive plant is another threat to E. argus populations. Recently, the abundance of invasive plant has been continually increasing in coastal dunes, resulting in negative effects on dune environments and the biodiversity present (Choi et al. 2006; Wouters et al. 2012). In general, the presence of invasive plants in coastal dunes is a result of topological changes and artificial introduction by visitors (Song and Cho 2007). Furthermore, the invasive plants change insect biodiversity and population dynamics in coastal dunes, which may also affect food availability for E. argus (Taean 2016a, b). Invasive plant removal and plant cover rate regulation are important for reptiles including E. argus, because open areas are necessary for basking (Kim et al. 2011a; Pike et al. 2011). Therefore, a consistent increase in invasive plants can have both direct and indirect effects on E. argus populations; thus, continuous management and removal of invasive species are required.

Vegetation cover of habitats is critical for E. argus (Kim et al. 2011a). The cover rate is related to the availability of moving and refuge spaces, degree of sun exposure, and availability of nearby food sources associated with dune vegetation (Kim et al. 2011a). In particular, the extent of light transmission by rice (Gramineae) plants makes the areas associated with these plants important refuges that protect from sun exposure; thus, these plants affect the presence and distribution of E. argus. However, the introduction of invasive plant could result in not only an increase in vegetation cover in coastal sand dunes but also reduce the available habitat of the lizards (Seo 2012).

In the 2008 study, all sampled E. argus individuals had their toes clipped (Song et al. 2013). In the current study, which was conducted in the same area, we did not find any previously marked E. argus individuals. It is possible that previously marked individuals died owing to various causes, such as disease, predation, and old age. A previous study estimated the life expectancy of E. argus to be 8–10 years (Kim et al. 2010), suggesting that 8 years is roughly equal to one generation. Generational changes are another important consideration for habitat or population restoration efforts.

In this study, we found that the expansion of habitat due to coastal dune restoration in Baramarye Coast can lead to an increase in the population size of E. argus, which is an endangered species. Our findings provide evidence that in situ conservation in the form of habitat expansion can be effective in preserving organisms that inhabit limited environments, such as coastal dunes.