Hierarchic species-area relationships and the management of forest habitat islands in intensive farmland

Habitat loss and fragmentation due to land use changes are major threats to biodiversity in forest ecosystems, and they are expected to have important impacts on many taxa and at various spatial scales. Species richness and area relationships (SARs) have been used to assess species diversity patterns and drivers, and thereby in the establishment of conservation and management strategies. Here we propose a hierarchical approach to achieve deeper insights on SARs in small forest islets in intensive farmland and to address the impacts of decreasing naturalness on such relationships. In the intensive dairy landscapes of Northwest Portugal, where small forest stands (dominated by pines, eucalypts or both) represent semi-natural habitat islands, 50 small forest stands were selected and surveyed for vascular plant diversity. A hierarchical analytical framework was devised to determine species richness and inter- and intra-patch SARs for the whole set of forest patches (general patterns) and for each type of forest (specific patterns). Differences in SARs for distinct groups were also tested by considering subsets of species (native, alien, woody, and herbaceous). Overall, values for species richness were confirmed to be different between forest patches exhibiting different levels of naturalness. Whereas higher values of plant diversity were found in pine stands, higher values for alien species were observed in eucalypt stands. Total area of forest (inter-patch SAR) was found not to have a significant impact on species richness for any of the targeted groups of species. However, significant intra-patch SARs were obtained for all groups of species and forest types. A hierarchical approach was successfully applied to scrutinise SARs along a gradient of forest naturalness in intensively managed landscapes. Dominant canopy tree and management intensity were found to reflect differently on distinct species groups as well as to compensate for increasing stand area, buffering SARs among patches, but not within patches. Thus, the maintenance of small semi-natural patches dominated by pines, under extensive practices of forest management, will promote native plant diversity while at the same time contributing to limit the expansion of problematic alien invasive species.

Slopes of Avian Species-Area Relationships, Human Population Density, and Environmental Factors

There is increasing interest in how humans influence spatial patterns in biodiversity. One of the most frequently noted and marked of these patterns is the increase in species richness with area, the species–area relationship (SAR). SARs are used for a number of conservation purposes, including predicting extinction rates, setting conservation targets, and identifying biodiversity hotspots. Such applications can be improved by a detailed understanding of the factors promoting spatial variation in the slope of SARs, which is currently the subject of a vigorous debate. Moreover, very few studies have considered the anthropogenic influences on the slopes of SARs; this is particularly surprising given that in much of the world areas with high human population density are typically those with a high number of species, which generates conservation conflicts. Here we determine correlates of spatial variation in the slopes of species–area relationships, using the British avifauna as a case study. Whilst we focus on human population density, a widely used index of human activities, we also take into account (1) the rate of increase in habitat heterogeneity with increasing area, which is frequently proposed to drive SARs, (2) environmental energy availability, which may influence SARs by affecting species occupancy patterns, and (3) species richness. We consider environmental variables measured at both local (10 km × 10 km) and regional (290 km × 290 km) spatial grains, but find that the former consistently provides a better fit to the data. In our case study, the effect of species richness on the slope SARs appears to be scale dependent, being negative at local scales but positive at regional scales. In univariate tests, the slope of the SAR correlates negatively with human population density and environmental energy availability, and positively with the rate of increase in habitat heterogeneity. We conducted two sets of multiple regression analyses, with and without species richness as a predictor. When species richness is included it exerts a dominant effect, but when it is excluded temperature has the dominant effect on the slope of the SAR, and the effects of other predictors are marginal.

Sampling-Design Effects on Properties of Species-Area Relationships - A Case Study from Estonian Dry Grassland Communities

Despite widespread use of species-area relationships (SARs), dispute remains over the most representative SAR model. Using data of small-scale SARs of Estonian dry grassland communities, we address three questions: (1) Which model describes these SARs best when known artifacts are excluded? (2) How do deviating sampling procedures (marginal instead of central position of the smaller plots in relation to the largest plot; single values instead of average values; randomly located subplots instead of nested subplots) influence the properties of the SARs? (3) Are those effects likely to bias the selection of the best model? Our general dataset consisted of 16 series of nested-plots (1 cm(2)-100 m(2), any-part system), each of which comprised five series of subplots located in the four corners and the centre of the 100-m(2) plot. Data for the three pairs of compared sampling designs were generated from this dataset by subsampling. Five function types (power, quadratic power, logarithmic, Michaelis-Menten, Lomolino) were fitted with non-linear regression. In some of the communities, we found extremely high species densities (including bryophytes and lichens), namely up to eight species in 1 cm(2) and up to 140 species in 100 m(2), which appear to be the highest documented values on these scales. For SARs constructed from nested-plot average-value data, the regular power function generally was the best model, closely followed by the quadratic power function, while the logarithmic and Michaelis-Menten functions performed poorly throughout. However, the relative fit of the latter two models increased significantly relative to the respective best model when the single-value or random-sampling method was applied, however, the power function normally remained far superior. These results confirm the hypothesis that both single-value and random-sampling approaches cause artifacts by increasing stochasticity in the data, which can lead to the selection of inappropriate models.

Cichlid species-area relationships are shaped by adaptive radiations that scale with area

A positive relationship between species richness and island size is thought to emerge from an equilibrium between immigration and extinction rates, but the influence of species diversification on the form of this relationship is poorly understood. Here, we show that within-lake adaptive radiation strongly modifies the species-area relationship for African cichlid fishes. The total number of species derived from in situ speciation increases with lake size, resulting in faunas orders of magnitude higher in species richness than faunas assembled by immigration alone. Multivariate models provide evidence for added influence of lake depth on the species-area relationship. Diversity of clades representing within-lake radiations show responses to lake area, depth and energy consistent with limitation by these factors, suggesting that ecological factors influence the species richness of radiating clades within these ecosystems. Together, these processes produce lake fish faunas with highly variable composition, but with diversities that are well predicted by environmental variables.

Out of Amazonia again and again: episodic crossing of the Andes promotes diversification in a lowland forest flycatcher

Most Neotropical lowland forest taxa occur exclusively on one side of the Andes despite the availability of appropriate habitat on both sides. Almost all molecular phylogenies and phylogenetic analyses of species assemblages (i.e. area cladograms) have supported the hypothesis that Andean uplift during the Late Pliocene created a vicariant barrier affecting lowland lineages in the region. However, a few widespread plant and animal species occurring in lowland forests on both sides of the Andes challenge the generality of this hypothesis. To understand the role of the Andes in the history of such organisms, we reconstructed the phylogeographic history of a widespread Neotropical flycatcher (Mionectes oleagineus) in the context of the other four species in the genus. A molecular phylogeny based on nuclear and mitochondrial sequences unambiguously showed an early basal split between montane and lowland Mionectes. The phylogeographic reconstruction of lowland taxa revealed a complex history, with multiple cases in which geographically proximate populations do not represent sister lineages. Specifically, three populations of M. oleagineus west of the Andes do not comprise a monophyletic clade; instead, each represents an independent lineage with origins east of the Andes. Divergence time estimates suggest that at least two cross-Andean dispersal events post-date Andean uplift.

Transferência de oxigênio em reatores de leito fluidizado com circulação em tubos concêntricos em meios bifásico e trifásico com variação da relação entre diâmetros

A eficiência do reator de leito fluidizado com circulação em tubos concêntricos depende das condições hidrodinâmicas que influem na transferência de oxigênio ao biofilme. Este trabalho investigou a influência da relação entre diâmetros dos tubos e da concentração de meio suporte (areia), sobre o coeficiente global de transferência de oxigênio (K La). Os ensaios - em reatores de 2,6 m de altura, com diâmetro externo de 250 mm e internos de 100, 125, 150 e 200 mm - empregaram vazões de ar até 2.500 L.h-1 e concentrações de até 150 g.L-1 de areia. O K La aumentou ligeiramente com 30 g.L-1 e diminuiu para concentrações maiores, confirmando relatos da literatura em condições semelhantes. Um modelo para K La em meio bifásico foi ajustado para as diversas relações ensaiadas entre a área externa e a interna, postulando-se uma redução na razão entre a transferência na fase líquida e o diâmetro da bolha com o aumento da vazão de ar.

Unraveling the drivers of community dissimilarity and species extinction in fragmented landscapes

Communities in fragmented landscapes are often assumed to be structured by species extinction due to habitat loss, which has led to extensive use of the species-area relationship (SAR) in fragmentation studies. However, the use of the SAR presupposes that habitat loss leads species to extinction but does not allow for extinction to be offset by colonization of disturbed-habitat specialists. Moreover, the use of SAR assumes that species richness is a good proxy of community changes in fragmented landscapes. Here, we assessed how communities dwelling in fragmented landscapes are influenced by habitat loss at multiple scales; then we estimated the ability of models ruled by SAR and by species turnover in successfully predicting changes in community composition, and asked whether species richness is indeed an informative community metric. To address these issues, we used a data set consisting of 140 bird species sampled in 65 patches, from six landscapes with different proportions of forest cover in the Atlantic Forest of Brazil. We compared empirical patterns against simulations of over 8 million communities structured by different magnitudes of the power-law SAR and with species-specific rules to assign species to sites. Empirical results showed that, while bird community composition was strongly influenced by habitat loss at the patch and landscape scale, species richness remained largely unaffected. Modeling results revealed that the compositional changes observed in the Atlantic Forest bird metacommunity were only matched by models with either unrealistic magnitudes of the SAR or by models ruled by species turnover, akin to what would be observed along natural gradients. We show that, in the presence of such compositional turnover, species richness is poorly correlated with species extinction, and z values of the SAR strongly underestimate the effects of habitat loss. We suggest that the observed compositional changes are driven by each species reaching its individual extinction threshold: either a threshold of forest cover for species that disappear with habitat loss, or of matrix cover for species that benefit from habitat loss.

Biogeographische Beziehungen zwischen den Alpen, dem Kaukasus und den asiatischen Hochgebirgen

Hintergrund: Biogeographische Verbindungen der Alpen nach Asien wurden schon im 19. Jahrhundert durch floristische Vergleiche entdeckt. Ziele: In den vorliegenden drei Artikeln wird untersucht, ob diese Verbindung entlang einer nördlichen, arktisch-borealen Route oder entlang einer südlichen Route über dazwischen liegende Gebirge wie den Kaukasus bestand. Methoden: Molekulare Phylogenien von Epimedium (Berberidaceae) und der Tribus der Hyoscyameae (Solanaceae), deren Vertreter disjunkt in den Alpen, dem Kaukasus und in asiatischen Gebirgen vorkommen, wurden mit nukleären und plastidären Markern erstellt und bio¬geographisch ausgewertet. Zur Datierung der Diversifizierungsereignisse und Arealbil¬dung diente eine molekulare Uhr. Ein aktueller floristischer Vergleich von Gattungen, die in den Alpen, den asiatischen Gebirgen und in Gebieten entlang der potentiellen Routen vor¬kom¬men, wurde unternommen und ausgewertet. Daran anschließend wurde nach molekular-phylogenetischer Literatur für interessante Gruppen aus diesem Ver¬gleich recherchiert, um diese biogeogra¬phisch bezüglich der Fragestellung zu interpretieren. Ergebnisse: Von 429 Gattungen, die in den Alpen und im Himalaya vorkommen, wachsen 218 entlang der nördlichen und der südlichen Route. 203 kommen nur entlang der süd¬lichen und drei nur entlang der nördlichen Route vor. Fünf kommen nur in den Alpen und im Himalaya vor. Epimedium, Scopolia/Physochlaina aus den Hyoscyameae und Primula sect. Auricula sind Beispiele für eine nördliche biogeographische Verbindung zwischen den Alpen und Asien. Diese bestand zumindest in den ersten beiden Fällen wahrscheinlich aus einem durchgängigen Laubwaldgürtel, der durch das abkühlende und trockener werdende Klima im Pliozän und Pleistozän fragmentiert wurde und heute nicht mehr besteht. Es handelt sich also um ein Vikarianzmuster und weniger um eine Migrationsroute. Die größere Diver¬sität dieser Guppen, die in Asien beobachtet werden kann ist sekundär entstanden. Atropa aus den Hyoscyameae, Brachypodium, die Subtribus Loliinae aus den Poaceae, Bupleurum und Doronicum sind Beispiele für eine südliche Verbindung. Hier liegen die Diversitätszentren im Mediterraneum und die Vorkommen im Osten sind abgeleitet. Schlussfolgerungen: Die Datengrundlage aus der Literaturrecherche ist nicht sehr breit und die meisten Phylogenien waren für die Fragestellung nicht aussagekräftig. Die histori¬schen Erkenntnisse und die Ideen über den Zusammenhang der Alpenflora mit der asiatischer Gebirge werden grundsätzlich bestärkt. Die Flora der Alpen enthält ein Element das über eine nördliche Verbindung noch lange mit asiatischen Gebirgen in Kontakt stand und ein südliches Element, das aus dem Mittelmeer¬raum bzw. aus SW Asien stammt. Die Beispiele für eine nördliche Verbindung sprechen allerdings nicht für eine Migration von Asien nach Europa sondern zeigen ein Vikarianzmuster auf.

Dry grasslands of NW Bulgarian mountains: first insights into diversity, ecology and syntaxonomy

We present the data of the 3rdresearch expedition of the European Dry Grasslands Group (EDGG), which was conducted in 2011 in two contrasting areas in NW Bulgarian mountains. The aim was to collect plot data for comparing Bulgarian dry grasslands with those of other parts of Europe in terms of syntaxonomy and biodiversity. We sampled 15 nested-plot series (0.0001–100 m²) and 68 normal plots(10 m²) covering the full variety of dry grassland types occurring in the Vratsa area (Balkan Mts.) and the Koprivshtitsa area (Sredna Gora Mt.). In the plots all vascular plants, terricolous non-vascular plants and a set of soil and other environmental parameters were determined. By applying modified TWIN-SPAN, we distinguished 10 floristically well characterised vegetation types at the association level. After comparison with the regional and European literature, we propose to place them within three classes and five orders: Festuco-Brometea with the orders Stipo pulcherrimae-Festucetalia pallentis (xerophilous dry grasslands of base-rich rocks; alliance Saturejion montanae), Brachypodietalia pinnate (meso-xeric, basiphilous grasslands; alliances Cirsio-Brachypodion pinnate and Chyrsopogono grylli-Danthonion calycinae),Calluno-Ulicetea with the order Nardetalia stricae (lowland to montane Nardus swards; alliance Violion caninae), and Koelerio-Corynephoretea with the orders Sedo-Scleranthetalia (open communities of skeleton-rich, acidic soils; alliance unclear) and Trifolioarvensis-Festucetalia ovinae(closed, meso-xeric, acidophilous grasslands; alliance Armerio rumelicae-Potentillion). The Violion caninae with the association Festuco rubrae-Genistelletum sagittalisis reported from Bulgaria for the first time, while the two occurring Koelerio-Corynephoretea communities are described as new associations (Cetrario aculeatae-Plantaginetum radicatae, Plantagini radicatae-Agrostietum capillaris). According to DCA the main floristic gradient was largely determined by soil conditions, differentiating the Festuco-Brometea communities on soils with high pH and high humus content from the Koelerio-Corynephoretea communities on acidic, humus-poor soils, while the Calluno-Ulicetea stands are the connecting link. At 10 m² Festuco-Brometea and Calluno-Ulicetea stands were richer in species across all investigated taxa and in vascular plants than Koelerio-Corynephoretea stands; the latter were richest in lichen species, while bryophyte richness did not differ significantly among syntaxa. Among the Bulgarian classes, the species-area relationships tended to be steepest in the Festuco-Brometea (i.e. highest beta diversity), but both alpha and beta diversity clearly fell behind the Festuco-Brometea communities in the Transylvanian Plateau, Romania, located less than 500 km north of the study region. Overall, our study contributes to a more adequate placement of the Bulgarian dry grasslands in the European syntaxonomic system and provides valuable data for large-scale analyses of biodiversity patterns

Sampling sites, sampled birds and plants of young larch plantations in Tokachi district, eastern Hokkaido, northern Japan, 2011

Models and data used to describe species-area relationships confound sampling with ecological process as they fail to acknowledge that estimates of species richness arise due to sampling. This compromises our ability to make ecological inferences from and about species-area relationships. We develop and illustrate hierarchical community models of abundance and frequency to estimate species richness. The models we propose separate sampling from ecological processes by explicitly accounting for the fact that sampled patches are seldom completely covered by sampling plots and that individuals present in the sampling plots are imperfectly detected. We propose a multispecies abundance model in which community assembly is treated as the summation of an ensemble of species-level Poisson processes and estimate patch-level species richness as a derived parameter. We use sampling process models appropriate for specific survey methods. We propose a multispecies frequency model that treats the number of plots in which a species occurs as a binomial process. We illustrate these models using data collected in surveys of early-successional bird species and plants in young forest plantation patches. Results indicate that only mature forest plant species deviated from the constant density hypothesis, but the null model suggested that the deviations were too small to alter the form of species-area relationships. Nevertheless, results from simulations clearly show that the aggregate pattern of individual species density-area relationships and occurrence probability-area relationships can alter the form of species-area relationships. The plant community model estimated that only half of the species present in the regional species pool were encountered during the survey. The modeling framework we propose explicitly accounts for sampling processes so that ecological processes can be examined free of sampling artefacts. Our modeling approach is extensible and could be applied to a variety of study designs and allows the inclusion of additional environmental covariates.

Estimate of the total ice volume of Svalbard glaciers and their potential contribution to sea-level rise = Estimación del volumen de hielo de Svalbard y su contribución potencial a la subida del nivel del mar

El objetivo final de las investigaciones recogidas en esta tesis doctoral es la estimación del volumen de hielo total de los ms de 1600 glaciares de Svalbard, en el Ártico, y, con ello, su contribución potencial a la subida del nivel medio del mar en un escenario de calentamiento global. Los cálculos más exactos del volumen de un glaciar se efectúan a partir de medidas del espesor de hielo obtenidas con georradar. Sin embargo, estas medidas no son viables para conjuntos grandes de glaciares, debido al coste, dificultades logísticas y tiempo requerido por ellas, especialmente en las regiones polares o de montaña. Frente a ello, la determinación de áreas de glaciares a partir de imágenes de satélite sí es viable a escalas global y regional, por lo que las relaciones de escala volumen-área constituyen el mecanismo más adecuado para las estimaciones de volúmenes globales y regionales, como las realizadas para Svalbard en esta tesis. Como parte del trabajo de tesis, hemos elaborado un inventario de los glaciares de Svalbard en los que se han efectuado radioecosondeos, y hemos realizado los cálculos del volumen de hielo de más de 80 cuencas glaciares de Svalbard a partir de datos de georradar. Estos volúmenes han sido utilizados para calibrar las relaciones volumen-área desarrolladas en la tesis. Los datos de georradar han sido obtenidos en diversas campañas llevadas a cabo por grupos de investigación internacionales, gran parte de ellas lideradas por el Grupo de Simulación Numérica en Ciencias e Ingeniería de la Universidad Politécnica de Madrid, del que forman parte la doctoranda y los directores de tesis. Además, se ha desarrollado una metodología para la estimación del error en el cálculo de volumen, que aporta una novedosa técnica de cálculo del error de interpolación para conjuntos de datos del tipo de los obtenidos con perfiles de georradar, que presentan distribuciones espaciales con unos patrones muy característicos pero con una densidad de datos muy irregular. Hemos obtenido en este trabajo de tesis relaciones de escala específicas para los glaciares de Svalbard, explorando la sensibilidad de los parámetros a diferentes morfologías glaciares, e incorporando nuevas variables. En particular, hemos efectuado experimentos orientados a verificar si las relaciones de escala obtenidas caracterizando los glaciares individuales por su tamaño, pendiente o forma implican diferencias significativas en el volumen total estimado para los glaciares de Svalbard, y si esta partición implica algún patrón significativo en los parámetros de las relaciones de escala. Nuestros resultados indican que, para un valor constante del factor multiplicativo de la relacin de escala, el exponente que afecta al área en la relación volumen-área decrece según aumentan la pendiente y el factor de forma, mientras que las clasificaciones basadas en tamaño no muestran un patrón significativo. Esto significa que los glaciares con mayores pendientes y de tipo circo son menos sensibles a los cambios de área. Además, los volúmenes de la población total de los glaciares de Svalbard calculados con fraccionamiento en grupos por tamaño y pendiente son un 1-4% menores que los obtenidas usando la totalidad de glaciares sin fraccionamiento en grupos, mientras que los volúmenes calculados fraccionando por forma son un 3-5% mayores. También realizamos experimentos multivariable para obtener estimaciones óptimas del volumen total mediante una combinación de distintos predictores. Nuestros resultados muestran que un modelo potencial simple volumen-área explica el 98.6% de la varianza. Sólo el predictor longitud del glaciar proporciona significación estadística cuando se usa además del área del glaciar, aunque el coeficiente de determinación disminuye en comparación con el modelo más simple V-A. El predictor intervalo de altitud no proporciona información adicional cuando se usa además del área del glaciar. Nuestras estimaciones del volumen de la totalidad de glaciares de Svalbard usando las diferentes relaciones de escala obtenidas en esta tesis oscilan entre 6890 y 8106 km3, con errores relativos del orden de 6.6-8.1%. El valor medio de nuestras estimaciones, que puede ser considerado como nuestra mejor estimación del volumen, es de 7.504 km3. En términos de equivalente en nivel del mar (SLE), nuestras estimaciones corresponden a una subida potencial del nivel del mar de 17-20 mm SLE, promediando 19_2 mm SLE, donde el error corresponde al error en volumen antes indicado. En comparación, las estimaciones usando las relaciones V-A de otros autores son de 13-26 mm SLE, promediando 20 _ 2 mm SLE, donde el error representa la desviación estándar de las distintas estimaciones. ABSTRACT The final aim of the research involved in this doctoral thesis is the estimation of the total ice volume of the more than 1600 glaciers of Svalbard, in the Arctic region, and thus their potential contribution to sea-level rise under a global warming scenario. The most accurate calculations of glacier volumes are those based on ice-thicknesses measured by groundpenetrating radar (GPR). However, such measurements are not viable for very large sets of glaciers, due to their cost, logistic difficulties and time requirements, especially in polar or mountain regions. On the contrary, the calculation of glacier areas from satellite images is perfectly viable at global and regional scales, so the volume-area scaling relationships are the most useful tool to determine glacier volumes at global and regional scales, as done for Svalbard in this PhD thesis. As part of the PhD work, we have compiled an inventory of the radio-echo sounded glaciers in Svalbard, and we have performed the volume calculations for more than 80 glacier basins in Svalbard from GPR data. These volumes have been used to calibrate the volume-area relationships derived in this dissertation. Such GPR data have been obtained during fieldwork campaigns carried out by international teams, often lead by the Group of Numerical Simulation in Science and Engineering of the Technical University of Madrid, to which the PhD candidate and her supervisors belong. Furthermore, we have developed a methodology to estimate the error in the volume calculation, which includes a novel technique to calculate the interpolation error for data sets of the type produced by GPR profiling, which show very characteristic data distribution patterns but with very irregular data density. We have derived in this dissertation scaling relationships specific for Svalbard glaciers, exploring the sensitivity of the scaling parameters to different glacier morphologies and adding new variables. In particular, we did experiments aimed to verify whether scaling relationships obtained through characterization of individual glacier shape, slope and size imply significant differences in the estimated volume of the total population of Svalbard glaciers, and whether this partitioning implies any noticeable pattern in the scaling relationship parameters. Our results indicate that, for a fixed value of the factor in the scaling relationship, the exponent of the area in the volume-area relationship decreases as slope and shape increase, whereas size-based classifications do not reveal any clear trend. This means that steep slopes and cirque-type glaciers are less sensitive to changes in glacier area. Moreover, the volumes of the total population of Svalbard glaciers calculated according to partitioning in subgroups by size and slope are smaller (by 1-4%) than that obtained considering all glaciers without partitioning into subgroups, whereas the volumes calculated according to partitioning in subgroups by shape are 3-5% larger. We also did multivariate experiments attempting to optimally predict the volume of Svalbard glaciers from a combination of different predictors. Our results show that a simple power-type V-A model explains 98.6% of the variance. Only the predictor glacier length provides statistical significance when used in addition to the predictor glacier area, though the coefficient of determination decreases as compared with the simpler V-A model. The predictor elevation range did not provide any additional information when used in addition to glacier area. Our estimates of the volume of the entire population of Svalbard glaciers using the different scaling relationships that we have derived along this thesis range within 6890-8106 km3, with estimated relative errors in total volume of the order of 6.6-8.1% The average value of all of our estimates, which could be used as a best estimate for the volume, is 7,504 km3. In terms of sea-level equivalent (SLE), our volume estimates correspond to a potential contribution to sea-level rise within 17-20 mm SLE, averaging 19 _ 2 mm SLE, where the quoted error corresponds to our estimated relative error in volume. For comparison, the estimates using the V-A scaling relations found in the literature range within 13-26 mm SLE, averaging 20 _ 2 mm SLE, where the quoted error represents the standard deviation of the different estimates.

A theoretical exploration of catchment-scale sediment delivery

This paper proposes a theoretical explanation of the variations of the sediment delivery ratio (SDR) versus catchment area relationships and the complex patterns in the behavior of sediment transfer processes at catchment scale. Taking into account the effects of erosion source types, deposition, and hydrological controls, we propose a simple conceptual model that consists of two linear stores arranged in series: a hillslope store that addresses transport to the nearest streams and a channel store that addresses sediment routing in the channel network. The model identifies four dimensionless scaling factors, which enable us to analyze a variety of effects on SDR estimation, including (1) interacting processes of erosion sources and deposition, (2) different temporal averaging windows, and (3) catchment runoff response. We show that the interactions between storm duration and hillslope/channel travel times are the major controls of peak-value-based sediment delivery and its spatial variations. The interplay between depositional timescales and the travel/residence times determines the spatial variations of total-volume-based SDR. In practical terms this parsimonious, minimal complexity model could provide a sound physical basis for diagnosing catchment to catchment variability of sediment transport if the proposed scaling factors can be quantified using climatic and catchment properties.