Testing of the SIBERIA landscape evolution model using the Tin Camp Creek,Northern Territory,Australia, field catchment

The SIBERIA landscape evolution model was used to simulate the geomorphic development of the Tin Camp Creek natural catchment over geological time. Measured hydrology, erosion and geomorphic data were used to calibrate the SIBERIA model, which was then used to make independent predictions of the landform geomorphology of the study site. The catchment, located in the Northern Territory, Australia is relatively untouched by Europeans so the hydrological and erosion processes that shaped the area can be assumed to be the same today as they have been in the past, subject to the caveats regarding long‐term climate fluctuation. A qualitative, or visual comparison between the natural and simulated catchments indicates that SIBERIA can match hillslope length and hillslope profile of the natural catchments. A comparison of geomorphic and hydrological statistics such as the hypsometric curve, width function, cumulative area distribution and area–slope relationship indicates that SIBERIA can model the geomorphology of the selected Tin Camp Creek catchments. Copyright 2002 © Environmental Research Institute of the Supervising Scientist, Commonwealth of Australia.     

The use of a landscape simulator in the validation of the SIBERIA landscape evolution model: transient landforms

SIBERIA is a physically based model for the geomorphic evolution of landforms. It is essential that the SIBERIA model be tested or validated against controlled landform development. Previous studies have demonstrated that SIBERIA is able to simulate declining equilibrium landforms and in this paper we examine SIBERIA's ability to simulate landforms as they evolve to their declining equilibrium form. These landscapes are termed transient landforms. Landscapes generated by SIBERIA were compared to those produced by a physical model (experimental model landforms) at stages of evolution. Comparison of the experimental landscapes with the simulated landscapes using total mass, hypsometric curve, width function, cumulative area distribution and area–slope demonstrate that SIBERIA can simulate the experimental model landscape during development (i.e. transient landscapes). Copyright © 2002 John Wiley & Sons, Ltd.     

Modelling of sediment dynamics in a laboratory‐scale experimental catchment

The variability of hillslope form and function is examined experimentally using a simple model catchment in which most landscape development parameters are either known or controlled. It is demonstrated that there is considerable variability in sediment output from similar catchments, subjected to the same hydrological processes, and for which the initial hillslope profiles are the same. The results demonstrate that, in the case of catchments with a linear initial hillslope profile, the sediment output is initially high but reduces through time, whereas for a concave initial profile the sediment output was smaller and relatively constant. Concave hillslope profiles also displayed reduced sediment output when compared with linear slopes with the same overall slope. Using this experimental model catchment data, the SIBERIA landscape evolution model was tested for its ability to predict temporal sediment transport. When calibrated for the rainfall and erodible material, SIBERIA is able to simulate mean temporal sediment output for the experimental catchment over a range of hillslope profiles and rainfall intensities. SIBERIA is also able to match the hillslope profile of the experimental catchments. The results of the study provide confidence in the ability of SIBERIA to predict temporal sediment output. The experimental and modelling data also demonstrate that, even with all geomorphic and hydrological variables being known and/or controlled, there is still a need for long‐term stream gauging to obtain reliable assessments of field catchment hydrology and sediment transport. Copyright © 2005 John Wiley & Sons, Ltd.     

Vegetation and cloud cover shape semi-arid carbonate landform development

We study the relationships among precipitation, vegetation, and morphological characteristics of watersheds draining either side of the Dhofar Mountains in southeastern Oman to understand the geomorphic signature of water availability in a semi-arid carbonate landscape. Water availability is expressed in terms of vegetation and cloud cover. The integral and the statistical moments of the hypsometric curve were used to determine whether hyper-arid, inland-draining watersheds are significantly different from seasonally wet watersheds on the coast side of the mountain range. We demonstrate that the vegetation and cloud cover are correlated, with locations with longer cloud periods also having a longer period with a vegetation canopy. The analysis shows that the hypsometric curve and its statistical moments capture the morphological difference between wet watersheds shaped by groundwater sapping and dry watersheds with fluvial morphology. Specifically, the curves exhibit two shapes: watersheds with more vegetation and cloud cover are characterized by higher convexity, and those with less vegetation and cloud cover are characterized by higher concavity. A variance analysis of cloud cover, vegetation, and hypsometric integral shows that they are significantly different between the wet and dry watersheds. The link between hydrology and morphology is not strong at the scale of a single watershed, but it is significant when the watersheds are aggregated in zones. The statistical moments of the hypsometric curve in the range of values of the integral and skewness show good separation between watersheds dominated by sapping and fluvial erosion processes. We can separate the watersheds draining the mountain range in two distinct groups on the basis of their bimodal hydrological and morphological characteristics. Our findings support other studies that hypothesize a trade-off from chemical- to mechanical-dominated denudation in carbonate terranes as precipitation decreases.     

Influence of bedrock structures on the spatial pattern of erosional landforms in small alpine catchments

Structural settings and lithological characteristics are traditionally assumed to influence the development of erosional landforms, such as gully networks and rock couloirs, in steep mountain rock basins. The structural control of erosion of two small alpine catchments of distinctive rock types is evaluated by comparing the correspondences between the orientations of their gullies and rock couloirs with (1) the sliding orientations of potential slope failures mechanisms, and (2) the orientation of the maximum joint frequency, this latter being considered as the direction exploited primarily by erosion and mass wasting processes. These characteristic orientations can be interpreted as structural weaknesses contributing to the initiation and propagation of erosion. The morphostructural analysis was performed using digital elevation models and field observations. The catchment comprised of magmatic intrusive rocks shows a clear structural control, mostly expressed through potential wedges failure. Such joint configurations have a particular geometry that encourages the development of gullies in hard rock, e.g. through enhanced gravitational and hydrological erosional processes. In the catchment underlain by sedimentary rocks, penetrative joints that act as structural weaknesses seem to be exploited by gullies and rock couloirs. However, the lithological setting and bedding configuration prominently control the development of erosional landforms, and influence not only the local pattern of geomorphic features, but the general morphology of the catchment. The orientations of the maximum joint frequency are clearly associated with the gully network, suggesting that its development is governed by anisotropy in rock strength. These two catchments are typical of bedrock‐dominated basins prone to intense processes of debris supply. This study suggests a quantitative approach for describing the relationship between bedrock jointing and geomorphic features geometry. Incorporation of bedrock structure can be relevant when studying processes governing the transfer of clastic material, for the assessment of sediment yields and in landforms evolution models. Copyright © 2012 John Wiley & Sons, Ltd.     

How do streamflow generation mechanisms affect watershed hypsometry?

Several studies have shown that the dominant streamflow generation mechanism in a river basin can leave distinct geomorphological signatures in basin topography. In particular, it has been suggested previously that basins generated by groundwater discharge tend to have a larger hypsometric integral than surface runoff basins because fluvial erosion is more focused in the valleys where groundwater discharge tends to occur. In this analysis, we aim to clarify this relationship by developing an alternative method to quantify the effects of streamflow generation mechanisms on basin hypsometry and by using a numerical model that can generate streamflow by different processes to evaluate the sensitivity of the results to the hydrological and geomorphological properties of the basin. The model results suggest that the hypsometric characteristics that are usually associated with groundwater discharge basins, such as a larger hypsometric integral, occur primarily when drainage networks are still advancing in the watershed. During later stages of development, an additional factor such as lithological controls or a distinct geomorphological process would be needed to preserve these features. The model results also show that the hypsometric effects are stronger when the parameters of the fluvial erosion process promote the influence of small discharge rates. Copyright © 2007 John Wiley & Sons, Ltd.     

Ephemeral landform development following rapid coastal uplift in the southern orogen of Taiwan

Newly emerged landscapes above sea level are characterized by rapidly evolving geomorphic systems where the initial fluvial pattern adapts to a former submarine topography. Such an early formed fluvial system establishes drainage basins and unstable landforms that characterize high topographic asymmetry which are prone to fast removal or reorganization. Transitional landscapes might form depositional systems as lakes or ponds that subsequently are incised, captured and incorporated into drainage basins. In this study we focus on the recently emerged Hengchun Peninsula to survey its paleoenvironment evolution. Three drillings performed in the Gangkou basin with fieldwork revealed several indicators that reconstructed stages of the landscape reorganization. The major finding shows an ephemeral large lake in the central part of the Hengchun Peninsula that was drained to the Pacific c. 6000 bp . The lake belonged to an ephemeral lakeland that was created after the emergence of the peninsula. Currently, several areas as relict landforms indicate this stage of topography evolution that through high rates of incision and subsequent captures, transforms into drainage basins. Furthermore, two drillings show brackish waters at the present estuary of the Gangkou basin. These two different paleoenvironments today build one system – Gangkou catchment. Long-term uplift rates show that a hanging wall of the Hengchun Fault plays a significant role in the creation of a lakeland by tilting the peninsula's surface. The tilt impacts on asymmetrical emergence of the peninsula and catchment development. Our study shows that a new geomorphic system might create depositional ephemeral landforms (lakes) that represent phases of early topography evolution after emergence above a sea level that are subjected to instantaneous rearrangement and evolves through large-scale phases before it reaches a topographic steady-state.     

A catchment scale evaluation of the SIBERIA and CAESAR landscape evolution models

Landscape evolution models provide a way to determine erosion rates and landscape stability over times scales from tens to thousands of years. The SIBERIA and CAESAR landscape evolution models both have the capability to simulate catchment–wide erosion and deposition over these time scales. They are both cellular, operate over a digital elevation model of the landscape, and represent fluvial and slope processes. However, they were initially developed to solve research questions at different time and space scales and subsequently the perspective, detail and process representation vary considerably between the models. Notably, CAESAR simulates individual events with a greater emphasis on fluvial processes whereas SIBERIA averages erosion rates across annual time scales. This paper describes how both models are applied to Tin Camp Creek, Northern Territory, Australia, where soil erosion rates have been closely monitored over the last 10 years. Results simulating 10 000 years of erosion are similar, yet also pick up subtle differences that indicate the relative strengths and weaknesses of the two models. The results from both the SIBERIA and CAESAR models compare well with independent field data determined for the site over different time scales. Representative hillslope cross‐sections are very similar between the models. Geomorphologically there was little difference between the modelled catchments after 1000 years but significant differences were revealed at longer simulation times. Importantly, both models show that they are sensitive to input parameters and that hydrology and erosion parameter derivation has long‐term implications for sediment transport prediction. Therefore selection of input parameters is critical. This study also provides a good example of how different models may be better suited to different applications or research questions. Copyright © 2010 John Wiley & Sons, Ltd and Commonwealth of Australia     

Tidal estuary width convergence: Theory and form in North Australian estuaries

In order to better understand the relations between tidal estuary shape and geomorphic processes, the width profiles of 79 tidal channels from within 30 estuaries in northern Australia have been extracted from LANDSAT 5 imagery using GIS. Statistics describing the shape and width convergence of individual channels and entire estuaries (which can contain several channels) are analysed along with proxies for the tidal range and fluvial inputs of the estuaries in question. The width profiles of most individual channels can be reasonably approximated with an exponential curve, and this is also true of the width profiles of estuaries. However, the shape of this exponential width profile is strongly related to the mouth width of the system. Channels and estuaries with larger mouths generally exhibit a more pronounced ‘funnel shape’ than those with narrower mouths, reflecting the hydrodynamic importance of the distance over which the channel or estuarine width converges. At the estuarine scale, this ‘convergence length’ also tends to be higher in estuaries which have larger catchments relative to their size. No clear relation between the estuarine width convergence length and tidal range could be discerned within the Northern Australian estuaries although, when these data are combined with data from other studies, a weak relationship emerges. Copyright © 2010 John Wiley & Sons, Ltd.     

Evaluating a new landform evolution model: A case study using a proposed mine rehabilitation landform

Landform evolution models are powerful tools for determining long-term erosional stability and denudation rates spanning geological timescales. SIBERIA, CAESAR and CHILD are examples of these model. The newly developed State Space Soil Production and Assessment Model (SSSPAM) coupled soilscape-landform evolution model has the ability to assess overall erosion rates of catchment scale landforms either using short-term precipitation events, variable precipitation or time-averaged precipitation (annual average). In addition, SSSPAM has the capability of developing the subsurface soil profile through weathering and armouring. In SSSPAM, physical processes of pedogenesis such as erosion and armouring, diffusion, sediment deposition and weathering are modelled using a state space matrix approach. In this article we simulate the short-term evolution (100 years) of a proposed post-mining landform using both SIBERIA and SSSPAM and compare the erosion and sediment output results. For the short-term simulations SSSPAM's armouring capability was disabled. The models were then used to simulate the evolution of the catchment for 10,000 years. Results demonstrate that the short-term SSSPAM simulation results compare well with the results from the established landform evolution model SIBERIA. The long-term armouring disabled SSSPAM simulations produces simulated erosion rates comparable with SIBERIA simulations both of which are similar to upper limit of field measured denudation rates. The SSSPAM simulation using armouring demonstrated that armouring reduced the erosion rate of the catchment by a factor of 4 which is comparable with the lower limit of field measured denudation rates. This observation emphasizes the importance of armouring in long-term evolution of landforms. Soil profile cross-sections developed from the same results show that SSSPAM can also reproduce subsurface soil evolution and stratification and spatial variability of soil profile characteristics typically observed in the field.     

Orogen‐scale drainage network evolution and response to erodibility changes: insights from numerical experiments

The continuous feedbacks among tectonics, surface processes, and climate are reflected in the distribution of catchments on active mountain ranges. Previous studies have shown a regularity of valley spacing across mountain ranges worldwide, but the origin of this geomorphological feature is currently not well known. In this work, we use a landscape evolution model to investigate the process of fluvial network organization and the evolution of regular ridge‐and‐valley patterns on simulated mountain ranges. In particular, we investigate the behavior of such patterns when subjected to a perturbation in landscape processes from a previous steady state, resulting from a sudden variation in the pattern of bedrock erodibility, from homogeneous to a gradient. We analyze the time evolution of the mean ratio λ' between the linear spacing of adjacent valleys and the half width of the mountain range. We show how a valley spacing ratio of ~0.5 is first achieved at steady state under uniform bedrock erodibility. After applying the gradient of bedrock erodibility across the landscape, we observe that λ' first increases and then decreases to a new steady‐state value that is smaller than the original value. A detailed analysis of the simulations, through observations of surface ‘snapshots’ at repeated time intervals, allows to gain some insight into the mechanisms governing this fluvial network reorganization process, driven by the migration of the main divide toward the side characterized by lower bedrock erodibility. On both sides of the range the new steady‐state valley spacing is obtained through mechanisms of catchment reorganization and competition between adjacent fluvial networks. In particular, catchment reorganization is characterized by the growth of smaller catchments between shrinking larger catchments on the side with lower erodibility, and the growth of larger catchments on the side with higher erodibility. Copyright © 2014 John Wiley & Sons, Ltd.     

Hydrologic response to valley‐scale structure in alpine headwaters

Few systematic studies of valley‐scale geomorphic drivers of streamflow regimes in complex alpine headwaters have compared response between catchments. As a result, little guidance is available for regional‐scale hydrological research and monitoring efforts that include assessments of ecosystem function. Physical parameters such as slope, elevation range, drainage area and bedrock geology are often used to stratify differences in streamflow response between sampling sites within an ecoregion. However, these metrics do not take into account geomorphic controls on streamflow specific to glaciated mountain headwaters. The coarse‐grained nature of depositional features in alpine catchments suggests that these landforms have little water storage capacity because hillslope runoff moves rapidly just beneath the rock mantle before emerging in fluvial networks. However, recent studies show that a range of depositional features, including talus slopes, protalus ramparts and ‘rock‐ice’ features may have more storage capacity than previously thought. To better evaluate potential differences in streamflow response among basins with extensive coarse depositional features and those without, we examined the relationships between streamflow discharge, stable isotopes, water temperature and the amplitude of the diurnal signal at five basin outlets. We also quantified the percentages of colluvial channel length measured along the stepped longitudinal profile. Colluvial channels, characterized by the presence of surficial, coarse‐grained depositional features, presented sediment‐rich, transport‐limited morphologies that appeared to have a cumulative effect on the timing and volume of flow downstream. Measurements taken from colluvial channels flowing through depositional landforms showed median recession constants (K_r) of 0.9–0.95, δ¹⁸O values of ≥?14.5 and summer diurnal amplitudes ≤0.8 as compared with more typical surface water recession constant values of 0.7, δ¹⁸O ≤ ?13.5 and diurnal amplitudes >2.0. Our results demonstrated strong associations between the percentage of colluvial channel length within a catchment and moderated streamflow regimes, water temperatures, diurnal signals and depleted δ¹⁸O related to groundwater influx. Copyright © 2014 John Wiley & Sons, Ltd.     

Application of a confluence‐based sediment‐fingerprinting approach to a dynamic sedimentary catchment,New Zealand

Fine sediment is a dynamic component of the fluvial system, contributing to the physical form, chemistry and ecological health of a river. It is important to understand rates and patterns of sediment delivery, transport and deposition. Sediment fingerprinting is a means of directly determining sediment sources via their geochemical properties, but it faces challenges in discriminating sources within larger catchments. In this research, sediment fingerprinting was applied to major river confluences in the Manawatu catchment as a broad‐scale application to characterizing sub‐catchment sediment contributions for a sedimentary catchment dominated by agriculture. Stepwise discriminant function analysis and principal component analysis of bulk geochemical concentrations and geochemical indicators were used to investigate sub‐catchment geochemical signatures. Each confluence displayed a unique array of geochemical variables suited for discrimination. Geochemical variation in upstream sediment samples was likely a result of the varying geological source compositions. The Tiraumea sub‐catchment provided the dominant signature at the major confluence with the Upper Manawatu and Mangatainoka sub‐catchments. Subsequent downstream confluences are dominated by the upstream geochemical signatures from the main stem of Manawatu River. Variability in the downstream geochemical signature is likely due to incomplete mixing caused in part by channel configuration. Results from this exploratory investigation indicate that numerous geochemical elements have the ability to differentiate fine sediment sources using a broad‐scale confluence‐based approach and suggest there is enough geochemical variation throughout a large sedimentary catchment for a full sediment fingerprint model. Combining powerful statistical procedures with other geochemical analyses is critical to understanding the processes or spatial patterns responsible for sediment signature variation within this type of catchment. Copyright © 2015 John Wiley & Sons, Ltd.     

Fluvial carbon flux from headwater peatland streams: significance of particulate carbon flux

The extensive blanket peatlands of the UK uplands account for almost half of total national terrestrial carbon storage. However, much of the blanket peat is severely eroded so that the contemporary role of the peatland system in carbon sequestration is compromised by losses of organic carbon in dissolved (DOC) and particulate (POC) form in the fluvial system. This paper presents the first detailed assessment of dissolved and organic carbon losses from a severely eroded headwater peatland (River Ashop, South Pennines, UK). Total annual fluvial organic carbon losses range from 29–106 Mg C km,^‐2 decreasing from the headwaters to the main catchment outlet. In contrast to less eroded systems fluvial organic carbon flux is dominated by POC. POC:DOC ratios decrease from values of 4 in the headwaters to close to unity at the catchment outlet. These results demonstrate the importance of eroding headwater sites as sources of POC to the fluvial system. Comparison with a range of catchment characteristics reveals that drainage density is the best predictor of POC:DOC but there is scatter in the relation in the headwaters. Steep declines in specific POC yield from headwater catchments are consistent with storage of POC within the fluvial system. Key to the significance of fluvial carbon flux in greenhouse gas budgets is understanding the fate of fluvial carbon. Further work on the fate of POC and the role of floodplains in fluvial carbon cycling is urgently required. Copyright © 2012 John Wiley & Sons, Ltd.     

Quantifying flows along hydrological pathways by applying a new filtering algorithm in conjunction with master recession curve analysis

Quantifying the proportion of the river hydrograph derived from the different hydrological pathways is essential for understanding the behaviour of a catchment. This paper describes a new approach using the output from master recession curve analysis to inform a new algorithm based on the Lyne and Hollick ‘one‐parameter’ signal analysis filtering algorithm. This approach was applied to six catchments (including two subcatchments of these) in Ireland. The conceptual model for each catchment consists of four main flow pathways: overland flow, interflow, shallow groundwater and deep groundwater. The results were compared with those of the master recession curve analysis, a recharge coefficient approach developed in Ireland and the semi‐distributed, lumped and deterministic hydrological model Nedbør‐Afstrømings‐Model. The new algorithm removes the ‘free variable’ aspect that is typically associated with filtering algorithms and provides a means of estimating the contribution of each pathway that is consistent with the results of hydrograph separation in catchments that are dominated by quick response pathways. These types of catchments are underlain by poorly productive aquifers that are not capable of providing large baseflows in the river. Such aquifers underlie over 73% of Ireland, ensuring that this new algorithm is applicable in the majority of catchments in Ireland and potentially in those catchments internationally that are strongly influenced by the quick‐responding hydrological pathways. Copyright © 2013 John Wiley & Sons, Ltd.     

STATISTICAL RELATIONS BETWEEN GEOMORPHIC PARAMETERS- A CASE STUDY OF THE YUNNAN REACH OF THE LANCANGJIANG RIVER IN SOUTHWESTERN CHINA

The properties of rivers and their catchments can be expressed by statistical relationships between geomorphic parameters. These statistical relationships may reveal some inherent differences in geomorphic evolution for different reaches or different order tributaries of a river basin. A case study was undertaken of the Yunnan reach of the Langcangjiang River. The catchment area, channel length and gradient of the first-, second- and third-order tributaries all with catchment areas larger than 100 km2 in the Yunnan reach were the main geomorphic parameters evaluated. The correlation between catchment area and channel length as well as between catchment area and channel gradient with respect to the total tributaries, different reach tributaries, and different order tributaries were revealed using statistical methods. In general, the channel length as a function of catchment area, was best expressed by a quadratic function where channel length increases with increasing catchment area (half parabola), while the channel gradient as a function of catchment area is best expressed by an exponential decay function. Comparison of the best-fit formulas revealed the following phenomena: the lower Yunnan reach tributaries and the first-order tributaries have a dominant effect on geomorphic parameters of the total tributaries. In addition, the statistical relationships indicate that the river geomorphic system in the upper and lower Yunnan reaches evolved differently. This study method used to differentiate river characteristics by determining statistical relationships between geomorphic parameters may be extended to other rivers and their catchments.     

Identification of calanco,a badland landform in the northern Apennines,Italy

Calanco (plural, calanchi) is a term widely used in the northern Apennines, Italy, to define a type of badland formed in clayey bedrock. However, no precise geomorphological definition of calanco has been established and a variety of map symbols are used to indicate the presence of calanco landforms. With the aim of developing an improved approach to identifying calanchi, a group of experienced surveyors identified 24 catchments with calanco characteristics among 67 catchments located between Bologna and Faenza in the northern Apennines. The morphology of each catchment was classified using traditional quantitative geomorphic approaches including fieldwork, map interpretation, hypsometric curve construction and computation of the annual sediment yield. Consideration of the parameters produced by these approaches indicated that none was capable of representing the presence of calanchi unequivocally and the basins were grouped into five classes on the basis of number and type of calanco criteria that they met. A characteristic of calanchi that is evident on topographic maps is crenulation of the contour lines and in this study a new topographic parameter was developed to represent the degree of contour crenulation. This parameter, LO/LF, is defined as the ratio of the actual length of a contour line (LO) to the length of the same line smoothed by an algorithm based on a moving average (LF). Calculated values of LO/LF ranged from 1·05 to 1·38. To test whether high values of the contour crenulation parameter were associated with calanchi, LO/LF values were added to other criteria for the five classes of catchment. Class 1 catchments, consisting of 14 of the 24 calanchi catchments identified in the field, displayed all of the criteria defining calanchi, and were characterized by the highest values of LO/LF (mean value 1·27 ± 0·15). It is proposed, therefore, that the contour crenulation ratio (LO/LF) may be useful in identifying the calanco landform. Copyright © 2000 John Wiley & Sons, Ltd.