Species-area Relationship Research Articles

Touristification and retail diversity: theory and evidence

This study explores the impact of touristification on retail diversity in neighbourhoods, presenting a species-area relationship model to demonstrate how tourism-driven changes can negatively affect local residents by displacing unique local retailers. Empirical regression analysis, using differences-in-difference (DID) models and neighbourhood clustering to define retail areas in Hong Kong, found that tourist-concentrated areas have fewer diverse high-street retailers. The result argues that the trend of touristification, which often leads to a decrease in retail diversity, could undermine the resilience of a tourism destination. This is evidenced by the noticeable downturn in the retail industry seen in many previously over-touristed cities after the pandemic. The findings prompt further research into tourism’s role in changing retail landscapes, particularly in touristified neighbourhoods.

Extending Species-Area Relationships Into the Realm of Ecoacoustics: The Soundscape-Area Relationship.

The rise in species richness with area is one of the few ironclad ecological relationships. Yet, little is known about the spatial scaling of alternative dimensions of diversity. Here, we provide empirical evidence for a relationship between the richness of acoustic traits emanating from a landscape, or soundscape richness, and island area, which we term the SoundScape-Area Relationship (SSAR). We show a positive relationship between the gamma soundscape richness and island area. This relationship breaks down at the smallest spatial scales, indicating a small-island effect. Moreover, we demonstrate a positive spatial scaling of the plot-scale alpha soundscape richness, but not the beta soundscape turnover, suggesting a direct effect of species on acoustic trait diversity. We conclude that the general scaling of biodiversity can be extended into the realm of ecoacoustics, implying soundscape metrics are sensitive to fundamental ecological patterns and useful in disentangling their complex mechanistic drivers.

Isolation limits spring pollination in a UK fragmented landscape

Animal-mediated pollination determines the reproductive success of most flowering plants; this process however can be disrupted by environmental degradation, with habitat loss and fragmentation highlighted as a top driver of pollination deficits. Despite being a pervasive stressor worldwide, we still have rather limited empirical evidence on its effects on pollination services, especially for early spring pollination syndromes. We investigate this using a potted plant phytometry experiment in which we placed English Bluebell (Hyacinthoides non-scripta)—a species largely pollinated in spring—into a fragmented woodland habitat. We selected 51 woodland patches which varied in both size and distance from each other and placed 153 pots of bluebell plants in the patches for c.4 weeks to measure pollination. The woodlands were located in a matrix of grassland, the latter being of low plant species richness and overall in the patches, woodland plants showed a positive species-area relationship. We collected traits on bluebell reproduction, these included the number, size, quality of seeds, the number of seed capsules and the number of flowers that failed to set any seeds. We found that seed traits responded differently to patch area and isolation. Patch isolation negatively affected the number of seeds and capsules, whilst it did not affect the size and quality of seeds. Patch area had no effect on any traits, suggesting that patch area might not necessarily be a factor that affects pollination in this species. The number of flowers that failed to set seed was unaffected by either patch area or isolation. Our study suggests that woodland fragmentation impacts the pollination of understory spring flowering plants. Our results highlight the use of multiple traits of phytometer plants to evaluate pollination and the importance of connectivity in maintaining pollination services in small-fragmented landscapes.

Assessing the potential of species loss caused by deforestation in a mature subtropical broadleaf forest in central China

Deforestation is a major type of land use change to accommodate growing population, especially in developing countries. The risk of diversity loss due to habitat loss can be estimated using the species-area relationship based on abundance of each species. However, deforestation often occurs before there is any understanding of the impact of deforestation on tree diversity. Here, we assessed the potential effect of forest habitat destruction on the loss of species richness in a mature subtropical broadleaf forest in central China. We surveyed and constructed the species-area relationship for 54 400 m2 plots, and simulated habitat loss scenarios by randomly and aggregately sampling plots. Rank-abundance of the 21 tree species was best fitted by the Zipf-Mandelbrot model, and our sample size was sufficient by the criterion of Hill numbers at orders q = 0, 1, and 2. We found that the number of species lost due to habitat loss was well predicted by the random placement species loss-area loss curve, and was lower than that due to aggregated habitat destruction by less than one species. The probability of losing one species reached 40% when losing 16 plots by aggregated sampling, 10 plots fewer than that by random sampling. Moreover, the probability of losing two species was 10–22 % higher by aggregately sampling than that by randomly sampling when losing 17 – 34 plots (0.68 – 1.36 ha). Considering that aggregated deforestation is common in reality, the results imply that the number of tree species lost due to deforestation could be higher than the theoretical estimation. Our study suggests the importance of assessing the impact of deforestation on tree diversity before selective logging in subtropical forests.

Landscape, Human Disturbance, and Climate Factors Drive the Species Richness of Alien Invasive Plants on Subtropical Islands.

Invasive alien plants (IAPs) pose a significant threat to island biodiversity and severely impact ecosystems. Understanding the species-area relationship and environmental determinants of growth forms for IAP species on subtropical islands is crucial for establishing an IAP's early warning mechanism, enhancing island ecological management, and protecting the ecosystems of Fujian and other subtropical islands. The study identified significant species-area relationships for IAPs and different life-form plants (trees, shrubs, and herbs), with slopes of 0.27, 0.16, 0.15, and 0.24, respectively. The small island effect does not apply to all species. Isolation has little effect on species richness, and the IAPs on Fujian islands do not conform to the isolation effect in island biogeography. Landscape factors are the main determinants of IAPs and different life-form species richness, with area, shape index, and perimeter-area ratio being the three primary landscape factors. These environmental factors are closely related to habitat heterogeneity. Besides landscape factors, different life forms respond differently to environmental factors. Climate drives the species richness distribution of shrubs and herbs, while trees are mainly influenced by human activities. Overall, landscape, human disturbance, and climate jointly drive the distribution of IAPs, with landscape factors being the most significant.

Fern Species-Area Relationship in Urban Anthropogenic Islands in Slawi, Tegal, Central Java

In anthropogenic islands as urban parks, the fern species richness and composition may be determined ecologically by the quality of habitat, including area greenness, or biogeographically by area size. As the development of the theory of island biogeography also includes man-made parks, it is feasible to test whether area-species relationship applies in these urban parks, and is more pronouncedly evident compared to another ecological factor, such as NDVI. Total species number and composition of 8 urban parks in a kecamatan in Tegal Regency were collected and arranged in clustering methods to understand the similarity between parks. The similarity analysis result is important for the management of the parks in Slawi. The species richness data is subsequently tested using Pearson correlation and regression against NDVI and area sizes. The relation between NDVI and species richness is non-significant (p=0.058), while area size and species richness is significant (p=0.003). This signifies that the urban fern species richness is determined by area as biogeographical factor, compared to NDVI as ecological factor. This result is important for the purpose of designing and managing urban parks as evidently size is important in the effort of attracting native biodiversity into urban parks and in turn enhancing the well-being of urban population.

Diversity of Ladybird Beetles (Coleoptera: Coccinellidae) in Tenerife and La Gomera (Canary Islands): The Role of Size and Other Island Characteristics.

This paper provides new data on the ladybird beetles (Coccinellidae) from two islands in the Canary archipelago: Tenerife, the largest island, and La Gomera, the second smallest. As they clearly differ in size but are similar in location and geological age, they are a suitable model for testing the species-area relationship. Our study shows that, in line with this main assumption of the theory of island biogeography, clearly more species occur on a large island (Tenerife) than on a small one (La Gomera). The field surveys documented the occurrence of 35 ladybird species on Tenerife (including 5 not previously reported from this island) and of 20 species on La Gomera (2 species new to the island). Coelopterus sp. collected on Tenerife (a single female that could not be identified to species) is the first record of this genus for the whole Canary Islands. Taking our data and previously published records into account, 47 species of Coccinellidae are known to occur on Tenerife and 26 species on La Gomera. Tenerife has by far the richest ladybird fauna of all the Canary Islands (the next in line, Gran Canaria, has 41 recorded species), but it also has the highest number of non-native ladybird species. All of the ten non-native species recorded in the Canary Islands are found on Tenerife, and for most of them, Tenerife was the island of their first appearance in the archipelago. This island, much more distant from the mainland than the other relatively large islands (Fuerteventura, Lanzarote), appears to be the main recipient of ladybirds immigrating to the Canary Islands. Tenerife can play this role probably because of its great habitat diversity and altitude variation, as well as intensive tourism and trade-related transport.

Scale dependence and mechanisms of grazing‐induced biodiversity changes depend on herbivore type in semiarid grasslands

Abstract Grasslands support a variety of herbivores that profoundly impact biodiversity and ecosystem functioning at local and regional scales. Understanding how different herbivores influence plant diversity across multiple scales is crucial for effective biodiversity conservation. However, most studies have focussed on the effects of grazing intensity on plant diversity, neglecting the impacts and associated mechanisms of different herbivores across scales. Based on a 7‐year grazing experiment with 16 plots in Inner Mongolia's typical steppe, we compared the scale‐dependent effects of different herbivores (cattle, sheep and goats) on plant diversity at moderate grazing intensity. We used the species–area relationship (SAR) and linear mixed models to analyse the changes in species richness across multiple scales (1, 16, 32, 64, 128, 256, 512 and 1024 m2). To understand SAR slope changes, we examined four potential mechanisms: overall species richness, total number of individuals, species abundance distribution and species aggregation. We found that grazing had scale effects on plant diversity, as indicated by a significant decrease in the SAR slope. This decrease in slope suggests that the negative impacts of grazing were more pronounced at larger scales. Additionally, the scale effects of different herbivores on plant diversity varied. Cattle had stronger positive impacts on diversity than sheep and goats at both small (1 m2) and large (1024 m2) scales. Conversely, sheep and goats had stronger negative effects on diversity than did the control at the large scale (1024 m2). The decrease in the SAR slope was mainly caused by changes in overall species richness and species abundance distribution. This was primarily because cattle had a stronger impact on dominant species, whereas sheep and goats had a stronger impact on rare species. Synthesis: Our results demonstrate that even at the same intensity, the grazing effects on species diversity vary with the study scale, herbivore type and their interaction. This finding has important implications for the conservation of grassland biodiversity given that grazing typically occurs at larger scales than in studies with quadrats and involves various herbivores with distinct effects on plant diversity.

Habitat quality drives the species–area relationship of plants and soil microbes in an ocean archipelago

Habitat quality drives the species–area relationship of plants and soil microbes in an ocean archipelago

Biodiversity on the line: life cycle impact assessment of power lines on birds and mammals in Norway

The global shift towards renewable energy plays an important role in fighting climate change. To facilitate the global growth of renewable energy production, the expansion of the electric grid becomes inevitable. Yet further construction of power lines poses a risk to biodiversity. Power lines traverse natural habitats and can lead to habitat conversion, fragmentation, and loss. Moreover, due to collisions and electrocutions, power lines kill hundreds of millions of birds each year. These impacts, however, have so far not been incorporated into decision-making. Life cycle assessment (LCA) is a widely used framework to compare environmental impacts and support decision-makers in planning and promoting sustainable strategies. We adapted existing life cycle impact assessment (LCIA) models to quantify the three main impacts of power lines on biodiversity: collision, electrocution, and habitat conversion and fragmentation. Our models incorporated species-area and species-habitat relationships to assess the effects of power lines on the diversity of birds and non-flying mammals in Norway, as the country is currently committed to reducing its emissions by expanding its renewable energy capacity. Overall, habitat conversion and fragmentation had the highest impact across the three impact categories, particularly affecting mammal richness. Furthermore, distribution lines often affected species richness more than transmission lines. The effect of the three impact pathways varied among different species groups, highlighting the vulnerability of certain species to habitat change, collision, and electrocution. Integrating LCIA models that quantify the impacts of power lines on biodiversity into LCA can support decision-makers with tools to promote the development of the electric grid without overlooking its effect on species richness. In addition, our models set the stage for a comprehensive assessment of the effects of electricity generation and transmission on biodiversity.

Power-law scaling in intratumoral microbiota of colorectal cancer

It has recently been proposed that the study of microbial dynamics in humans may gain insights from island biogeographical theory. Here, we test whether the diversity of the intratumoral microbiota of colorectal cancer tumors (CRC) follows a power law with tumor size akin to the island species-area relationship. We confirm a direct correlation between the quantity of Amplicon Sequence Variants (ASVs) within CRC tumors and tumor sizes, following a (log)power model, explaining 47% of the variation. Understanding the processes involved, potentially through the analogy of tumors and islands, may ultimately contribute to future clinical and therapeutic strategies.

New data on spiders (Arachnida, Araneae) from the islands of the Strait of Sicily (Southern Italy) with taxonomic notes on Poecilochroa loricata Kritscher, 1996 (Araneae, Gnaphosidae) and eight new records for Europe

We present a checklist of spiders inhabiting the Pelagie archipelago and Pantelleria island, in the Strait of Sicily (Southern Italy). Data were compiled from both literature sources and unpublished materials stored in museum collections. In total, we report new data on 100 species, bringing the total number of species documented for the islands of the Strait of Sicily to 148. Among these, 8 are new for Europe and 9 new for Italy. The island of Lampedusa hosts the highest number of species (107), followed by Pantelleria (63) and Linosa (25). The most represented families are Gnaphosidae in Lampedusa and Linosa, while Salticidae is dominant in Pantelleria. Most of the species present on the islands of the Strait of Sicily have a Holarctic distribution. Several rare species are recorded, including Haplodrassus crassipes (Lucas, 1846), Palliduphantes labilis (Simon, 1913) and Xysticus promiscuus O. Pickard-Cambridge, 1876. In addition, the female of Poecilochroa loricata Kritscher, 1996 (Gnaphosidae) is here described and illustrated for the first time. As revealed by the study of the species-area relationship, the islands of the Strait of Sicily host a higher number of species compared to other small islands in Italy. Accordingly, species richness recorded in Lampedusa and Linosa is higher than expected, while further investigations are suggested to increase knowledge of the local spider diversity.

Accounting for imperfect detection when estimating species-area relationships and beta-diversity.

Ecologists have historically quantified fundamental biodiversity patterns, including species-area relationships (SARs) and beta diversity, using observed species counts. However, imperfect detection may often bias derived community metrics and subsequent community models. Although several statistical methods claim to correct for imperfect detection, their performance in species-area and β-diversity research remains unproven. We examine inaccuracies in the estimation of SARs and β-diversity parameters that emerge from imperfect detection, and whether such errors can be mitigated using a non-parametric diversity estimator (iNEXT.3D) and Multi-Species Occupancy Models (MSOMs). We simulated 28,350 sampling regimes of 2835 fragmented communities, varying the mean and standard deviation of species detection probabilities, and the number of sampling repetitions. We then quantified the bias, accuracy, and precision of derived estimates of model coefficients for SARs and the effects of patch area on β-diversity (pairwise Sørensen similarity). Imperfect detection biased estimates of all evaluated parameters, particularly when mean detection probabilities were low, and there were few sampling repetitions. Observed counts consistently underestimated species richness and SAR z-values, and overestimated SAR c-values; iNEXT.3D and MSOMs only partially resolved these biases. iNEXT.3D provided the best estimates of SAR z-values, although MSOM estimates were generally comparable. All three methods accurately estimated pairwise Sørensen similarity in most circumstances, but only MSOMs provided unbiased estimates of the coefficients of models examining covariate effects on β-diversity. Even when using iNEXT.3D or MSOMs, imperfect detection consistently caused biases in SAR coefficient estimates, calling into question the robustness of previous SAR studies. Furthermore, the inability of observed counts and iNEXT.3D to estimate β-diversity model coefficients resulted from a systematic, area-related bias in Sørensen similarity estimates. Importantly, MSOMs corrected for these biases in β-diversity assessments, even in suboptimal scenarios. Nonetheless, as estimator performance consistently improved with increasing sampling repetitions, the importance of appropriate sampling effort cannot be understated.

Secondary production and biomass in mussel assemblages relate to species richness and stream size but not life history

AbstractIncreases in species richness with habitat area (species–area relationship, or SAR) and increases in ecosystem function with species richness (biodiversity–ecosystem functioning, or BEF) are widely studied ecological patterns. Incorporating functional trait analysis into assemblage datasets may help clarify interpretations of SAR and BEF relationships in natural ecological systems. For example, life history theory can be used to make predictions about what species are most important in generating ecosystem function given a certain set of environmental conditions. We used quantitative assemblage data for freshwater mussels at nine sites in western Alabama, USA, to test for SAR and BEF relationships. At each site, we calculated species richness, mussel assemblage density, and two fundamental metrics of ecosystem function: biomass and secondary production. We also tested whether the proportional biomass and production contributions from species belonging to each of three life history strategies—opportunistic strategists adapted to unstable or frequently disturbed habitats, periodic strategists adapted to habitats subject to predictable large‐scale disturbances, and equilibrium strategists adapted to stable habitats—varied longitudinally with stream drainage area, a proxy for habitat area. Species richness increased with stream size (SAR), and both biomass and production increased with species richness (BEF) and mussel density. There were few longitudinal changes in the proportional contributions of the different life history strategy classifications that we used, but the invasive clam Corbicula fluminea contributed proportionally more biomass and production at sites that had smaller drainage areas. This study provides further evidence for a clear longitudinal SAR in stream‐dwelling taxa. It also suggests BEF relationships for biomass and secondary production in natural assemblages but underscores the importance of assemblage density in BEF studies that use observational field data. Variation in proportional biomass and production contributions by different life history strategies was likely limited by the size of the stream size gradient in our study, as contributions were uniformly high for species with life history traits better adapted to stable and productive habitats such as mid‐sized rivers with low or predictable hydrologic disturbance frequencies. This highlights the need to understand how organisms' functional traits govern their relationships to the environment at different scales.

High levels of species' extirpation in an urban environment-A case study from Berlin, Germany, covering 1700-2023.

Species loss is highly scale-dependent, following the species-area relationship. We analysed spatio-temporal patterns of species' extirpation on a multitaxonomic level using Berlin, the capital city of Germany. Berlin is one of the largest cities in Europe and has experienced a strong urbanisation trend since the late nineteenth century. We expected species' extirpation to be exceptionally high due to the long history of urbanisation. Analysing 37 regional Red Lists of Threatened Plants, Animals and Fungi of Berlin (covering 9498 species), we found that 16% of species were extirpated, a rate 5.9 times higher than at the German scale and 47.1 times higher than at the European scale. Species' extirpation in Berlin is comparable to that of another German city with a similarly broad taxonomic coverage, but much higher than in regional areas with less human impact. The documentation of species' extirpation started in the eighteenth century and is well documented for the nineteenth and twentieth centuries. We found an average annual extirpation of 3.6 species in the nineteenth century, 9.6 species in the twentieth century and the same number of extirpated species as in the nineteenth century were documented in the twenty-first century, despite the much shorter time period. Our results showed that species' extirpation is higher at small than on large spatial scales, and might be negatively influenced by urbanisation, with different effects on different taxonomic groups and habitats. Over time, we found that species' extirpation is highest during periods of high human alterations and is negatively affected by the number of people living in the city. But, there is still a lack of data to decouple the size of the area and the human impact of urbanisation. However, cities might be suitable systems for studying species' extirpation processes due to their small scale and human impact.

Landscape configuration can flip species-area relationships in dynamic meta-food-webs.

Spatial and trophic processes profoundly influence biodiversity, yet ecological theories often treat them independently. The theory of island biogeography and related theories on metacommunities predict higher species richness with increasing area across islands or habitat patches. In contrast, food-web theory explores the effects of traits and network structure on coexistence within local communities. Exploring the mechanisms by which landscape configurations interact with food-web dynamics in shaping metacommunities is important for our understanding of biodiversity. Here, we use a meta-food-web model to explore the role of landscape configuration in determining species richness and show that when habitat patches are interconnected by dispersal, more species can persist on smaller islands than predicted by classical theory. When patch sizes are spatially aggregated, this effect flattens the slope of the species-area relationship. Surprisingly, when landscapes have random patch-size distributions, the slope of the species-area relationships can even flip and become negative. This could be explained by higher biomass densities of lower trophic levels that then support species occupying higher trophic levels, which only persist on small and well-connected patches. This highlights the importance of simultaneously considering landscape configuration and local food-web dynamics to understand drivers of species-area relationships in metacommunities.This article is part of the theme issue 'Diversity-dependence of dispersal: interspecific interactions determine spatial dynamics'.

Climatic and biogeographic factors show contrasted effects on continental and volcanic ISARs

Abstract Aim The study of islands biotas is fundamental to understand biodiversity patterns and process, both on evolutionary and ecological time scales. This study explores the influence of biogeographic and climatic factors on island species–area relationships (ISARs) in volcanic and continental islands, aiming to detect differences in slopes and intercepts between these island origins. Methods Data from 5049 vegetation plots on 58 Central Mediterranean Sea islands were collected from various sources. Islands were categorized as volcanic or continental based on their geological origin. area, isolation, maximum elevation, eccentricity, mean annual temperature and annual precipitation were calculated for each island. By using a moving window, we created groups of islands, and ISARs were fitted to each group using the Arrhenius power functions. Linear models and a permutation test were employed to examine how ISAR model parameters varied with the above-mentioned variables. Results While ISAR intercept values did not differ between island origins, volcanic islands showed higher ISAR slopes than continental islands. Whereas increasing island isolation and mean annual temperature increased ISAR intercepts on continental islands groups, it decreased on volcanic ones. Increasing annual precipitation decreased ISAR intercept on continental islands groups and increased on volcanic islands groups. Increasing island isolation, eccentricity and mean annual temperature increased the ISAR slope on volcanic islands groups while increasing annual precipitation decreased it. Increasing maximum elevation increased ISAR slope on continental islands groups. Main conclusions Our study provides evidence about the importance of island geological origin in determining the observed ISAR patterns. Biogeographic and climatic factors are pivotal in shaping species richness patterns on islands, exerting varying influences on both volcanic and continental islands.

Distribution of native plants and birds on 27 remote Tuamotu atolls (South Pacific Ocean): a contribution to the Island species-area relationship

ABSTRACT Island species-area relationship (ISAR) is the most documented pattern in island biogeography. Different hypotheses were advanced to explain this pattern. In this study, I selected 27 remote Tuamotu atolls in the South Pacific Ocean with complete surveys of native species richness of birds and vascular plants to test the influence of four abiotic predictors on species richness (atoll emerged area, habitat diversity, mean elevation, and number of islets). Linear regressions were used to assess the relative influence of predictors on native species richness while stepwise regression was then used to identify the best model. Atoll area was a significant predictor to explain native bird and plant species richness, attesting ISAR on the remote surveyed atolls. Stepwise model demonstrated that both habitat diversity and atoll area explained bird species richness, whereas atoll area and mean elevation were the best predictors for native plants. These results suggest that ISAR can be related to different hypotheses, depending on the taxon studied. Among hypotheses, the simple “target-area” hypothesis was a suitable framework to explain ISAR of native birds, while the “disturbance” hypothesis was relevant to support ISAR of native plants observed on the atolls.

Effect of Plateau pika on Soil Microbial Assembly Process and Co-Occurrence Patterns in the Alpine Meadow Ecosystem.

Burrowing animals are a critical driver of terrestrial ecosystem functioning, but we know little about their effects on soil microbiomes. Here, we evaluated the effect of burrowing animals on microbial assembly processes and co-occurrence patterns using soil microbiota from a group of habitats disturbed by Plateau pikas (Ochtona curzoniae). Pika disturbance had different impacts on bacterial and fungal communities. Fungal diversity generally increased with patch area, whereas bacterial diversity decreased. These strikingly different species-area relationships were closely associated with their community assembly mechanisms. The loss of bacterial diversity on larger patches was largely driven by deterministic processes, mainly due to the decline of nutrient supply (e.g., organic C, inorganic N). In contrast, fungal distribution was driven primarily by stochastic processes that dispersal limitation contributed to their higher fungal diversity on lager patches. A bacterial co-occurrence network exhibited a positive relationship of nodes and linkage numbers with patch area, and the fungal network presented a positive modularity-area relationship, suggesting that bacteria tended to form a closer association community under pika disturbance, while fungi tended to construct a higher modularity network. Our results suggest that pikas affects the microbial assembly process and co-occurrence patterns in alpine environments, thereby enhancing the current understanding of microbial biogeography under natural disturbances.

EDNA Metabarcoding Reveals the Species-Area Relationship of Amphibians on the Zhoushan Archipelago.

The species-area relationship is important for understanding species diversity patterns at spatial scales, but few studies have examined the relationship using environmental DNA (eDNA) techniques. We investigated amphibian diversity on 21 islands of the Zhoushan Archipelago and nearby mainland areas in China using the combination of eDNA metabarcoding and the traditional line transect method (TLTM) and identified the species-area relationship for amphibians on the islands. The mean detection probability of eDNA is 0.54, while the mean detection probability of TLTM is 0.24. The eDNA metabarcoding detected eight amphibian species on the islands and nine species in the mainland areas, compared with seven species on the islands and nine species in the mainland areas that were identified by TLTM. Amphibian richness on the islands increased with island area and habitat diversity. The species-area relationship for amphibians in the archipelago was formulated as the power function (S = 0.47A0.21) or exponential function (S = 2.59 + 2.41 (logA)). Our results suggested that eDNA metabarcoding is more sensitive for the detection of amphibian species. The combined use of eDNA metabarcoding and the traditional line transect method may optimize the survey results for amphibians.