scispace - formally typeset
Search or ask a question
Journal ArticleDOI

Aridity is expressed in river topography globally.

TL;DR: A global dataset of river longitudinal profiles is presented and it is shown that river profiles become straighter with increasing aridity and numerical modelling suggests that this can be explained by rainfall–runoff regimes in different climate zones.
Abstract: It has long been suggested that climate shapes land surface topography through interactions between rainfall, runoff and erosion in drainage basins1,2,3,4. The longitudinal profile of a river (elevation versus distance downstream) is a key morphological attribute that reflects the history of drainage basin evolution, so its form should be diagnostic of the regional expression of climate and its interaction with the land surface5,6,7,8,9. However, both detecting climatic signatures in longitudinal profiles and deciphering the climatic mechanisms of their development have been challenging, owing to the lack of relevant global data and to the variable effects of tectonics, lithology, land surface properties and human activities10,11. Here we present a global dataset of 333,502 river longitudinal profiles, and use it to explore differences in overall profile shape (concavity) across climate zones. We show that river profiles are systematically straighter with increasing aridity. Through simple numerical modelling, we demonstrate that these global patterns in longitudinal profile shape can be explained by hydrological controls that reflect rainfall–runoff regimes in different climate zones. The most important of these is the downstream rate of change in streamflow, independent of the area of the drainage basin. Our results illustrate that river topography expresses a signature of aridity, suggesting that climate is a first-order control on the evolution of the drainage basin.

Summary (2 min read)

METHODS

  • Having confirmation from two climatic indices (K-G and AI), which are computed in distinct ways (e.g., AI represents the balance between PET and P), gives us confidence that the authors have captured real climate influences on long profile development.
  • This was performed using an algorithm 34 which minimizes the topographic change required to ensure all DEM cells flow to the DEM base level.
  • This is the channel that was extracted in their analysis for this area and which is included in GLoPro.
  • Given their source in SRTM data, the extracted profiles represent the water surface profile for perennial rivers and the bed topography profile for ephemeral rivers.

Normalized Concavity Index (NCI).

  • The authors define the endpoints of the longitudinal profile (L0, E0) and (Ln, En) where L is distance downstream, E is elevation, and where the subscripts 0 and n indicate the most upstream and downstream points, respectively.
  • Then, at each measured point along the profile, the vertical offset between the river profile and the fitted straight line is calculated as EL -YL.
  • The authors then calculate the median value of all offsets, normalized by the total topographic relief along the profile (E0 -En) to enable comparison across scales (Extended Data Fig. 1b ).

NCI = median[(EL-YL)/(E0-En)]

  • (1) There have been previous concavity indices developed in the literature, such as Stream Concavity Index (SCI) 7 , Concavity Index (θ) 40 , and Chi (χ) transformation 41 .
  • SCI, for example, calculates the area between channel elevation and the straight line connecting the endpoints of channel, similar to NCI.
  • SCI is sensitive to local variations along the profile (e.g., knickpoints) and requires smoothing.
  • Since their goal was to explore conditions where the relationship between area and channel discharge are weak for complete river profiles, the authors opted for a different metric.

as an example).

  • The river extraction methods and concavity calculation result in an internally consistent NCI dataset.
  • The impact of channel head location on NCI is minimal because only the longest river of each basin or sub-basin was analyzed (not smaller tributaries).
  • The authors confirmed that NCI for extracted rivers in GLoPro are not correlated with key river metrics, such as river length, gradient, relief, or basin area (Extended Data Fig. 4 ).
  • Therefore, the authors were confident in using it to compare rivers of different sizes and across climate zones.

2. riverid:

  • The unique name given to each river record in GLoPro.
  • Comprises the K-G climate zone that the river is within and a unique alphanumeric string.
  • Used to identify the associated data for the river recorded in rivers.

SELECT elevation, length FROM profiles WHERE riverid like 'Aw_75_river_72';

  • Note that due to the size of the profiles table, queries can take a few minutes to complete.
  • To learn more about using SQL databases in a research context, the authors recommend the training materials provided by Software Carpentry: http://swcarpentry.github.io/sql-novice-survey.

Kernel density estimation (KDE).

  • In several figures in the paper, the authors present plots generated based on kernel density estimation (KDE).
  • Statistical differences of the NCI distributions were analyzed using the Kolmogorov-Smirnov test (K-S test) between distribution pairs across climate zones.
  • K-S test is a nonparametric test for checking whether two continuous, one-dimensional data samples, X1 and X2, come from the same distribution.
  • The authors simulated variations in downstream discharge and their impact on long profile evolution.
  • Since other model parameters can also affect long profile concavities, the authors conducted sensitivity analyses to discharge (Qmax), median grain size , tectonic uplift, and base level change.

Calculation of α values from real rivers.

  • To develop a real-world understanding of α and its variation in different climate zones, the authors downloaded multidecadal mean daily streamflow data for rivers from the US Geological Survey's National Water Information System (https://waterdata.usgs.gov/nwis).
  • For each main K-G climate zone, the authors selected 5 rivers, spanning a range of river lengths, with at least three gauging stations along the same river (a total of 20 rivers), ensuring via Google Earth satellite imagery that there are no obvious anthropogenic factors that could influence the downstream variation in discharge.
  • Then the authors extracted α for each power law fit from equation (3) (Extended Data Table 2 ).
  • The exponent for arid channels is closest to zero for small floods and increases slightly for higher flood recurrence intervals.
  • The analysis of α values was not exhaustive.

Did you find this useful? Give us your feedback

Content maybe subject to copyright    Report

This is an Open Access document downloaded from ORCA, Cardiff University's institutional
repository: http://orca.cf.ac.uk/124261/
This is the author’s version of a work that was submitted to / accepted for publication.
Citation for final published version:
Chen, Shiuan-An, Michaelides, Katerina, Grieve, Stuart W. D. and Singer, Michael Bliss 2019.
Aridity is expressed in river topography globally. Nature 573 , pp. 573-577. 10.1038/s41586-019-
1558-8 file
Publishers page: http://dx.doi.org/10.1038/s41586-019-1558-8 <http://dx.doi.org/10.1038/s41586-
019-1558-8>
Please note:
Changes made as a result of publishing processes such as copy-editing, formatting and page
numbers may not be reflected in this version. For the definitive version of this publication, please
refer to the published source. You are advised to consult the publisher’s version if you wish to cite
this paper.
This version is being made available in accordance with publisher policies. See
http://orca.cf.ac.uk/policies.html for usage policies. Copyright and moral rights for publications
made available in ORCA are retained by the copyright holders.

Aridity is expressed in river topography globally 1
2
Shiuan-An Chen
1
*, Katerina Michaelides
1,2
, Stuart W. D. Grieve
3
and Michael Bliss Singer
2,4,5
3
4
1
School of Geographical Sciences, University of Bristol, Bristol, BS8 1SS, UK.
2
Earth Research Institute, 5
University of California Santa Barbara, Santa Barbara, California 91306, USA.
3
University College 6
London, London, WC1E 6BT, UK.
4
School of Earth and Ocean Sciences, Cardiff University, Cardiff, CF10 7
3AT, UK.
5
Water Research Institute, Cardiff University, Cardiff, CF10 3AX, UK. 8
*e-mail: sc16970@bristol.ac.uk 9
10
It has long been suggested that climate shapes land surface topography, through interactions between 11
rainfall, runoff, and erosion in drainage basins
1-4
. The longitudinal profile of a river (elevation versus 12
distance downstream) is a key morphological attribute that reflects the history of drainage basin 13
evolution, so its form should be diagnostic of the regional expression of climate and its interaction with 14
the land surface
5-9
. However, both detecting climatic signatures in longitudinal profiles and 15
deciphering the climatic mechanisms of their development have been challenging due to the lack of 16
relevant data across the globe, and due to the variable effects of tectonics, lithology, land-surface 17
properties, and humans
10,11
. Here we present a global dataset of river longitudinal profiles (n = 18
333,502), and use it to explore differences in overall profile shape (concavity) across climate zones. 19
We show that river profiles are systematically straighter with increasing aridity. Through simple 20
numerical modeling, we demonstrate that these global patterns in longitudinal profile shape can be 21
explained by hydrological controls that reflect rainfall-runoff regimes in different climate zones. The 22
most important of these is the downstream rate-of-change in streamflow independent of drainage 23
basin area. Our results illustrate that river topography inherits a signature of aridity, suggesting that 24
climate is a first-order control on drainage basin evolution. 25

Conventional theory presents river longitudinal profiles (long profiles) as having a generally concave-up 26
shape, with knickpoints and other fluctuations expressing the interactions of several independent variables: 27
climate, tectonics, lithology, and human impacts
11-13
. This characteristic shape of long profiles has been 28
interpreted to arise due to downstream flow increase with drainage area, which erodes the riverbed, 29
transports sediment from upstream to downstream, and produces fining profiles in bed material grain 30
size
13,14
. However, there are long profiles with overall concavity much closer to zero (straighter) than the 31
typical concave-up profile shape
15-17
, yet there is limited understanding of the global distribution of long 32
profile concavities and their relation to climate. Stream power incision theory states that channel erosion is 33
intrinsically tied to an assumed relationship between river discharge (Q) and drainage area (Q~A
c
). Based 34
on this theory, an expression has been derived that links supply-limited river long profile concavity to the 35
exponent c
18
, illustrating that profiles will be concave up for c > 0, straight for c = 0, and convex for c < 0, 36
and a similar dependency of profile concavity on the Q-A relationship has been derived for 37
transport-limited fluvial systems
19
. Previous work has largely emphasized long profile concavity for cases 38
where c > 0, despite evidence that c in many river basins, especially in drylands, may vary flood to flood 39
between negative, zero, and positive values
8,17,20
. Of particular interest here is to ascertain whether the 40
climatic expression within river channel hydrology may be a first-order control on long profile shape, and 41
whether its climatic signature is preserved across the globe. 42
A river experiences a cascade: from climate to hydrology to erosion, which evolves its long profile. 43
Therefore, the climatic expression within streamflow should be a first-order control on long profile shape
6-8
. 44
Numerical analysis of profile shape responses to a distribution of flow events above the threshold for 45
bedrock incision has demonstrated part of this dependency
5,8,21
. However, there is limited global evidence 46
of how the hydrologic expression of climate affects long profiles, across a wide range of climate zones. 47
Climate determines the precipitation regime within a region. In turn, the precipitation regime controls the 48
rate and frequency of water supply to the land surface, a proportion of which generates runoff over drainage 49
basins, subject to losses by infiltration and evapotranspiration. Flow in rivers occurs when runoff reaches 50

the channel, with notable baseflow contributions from groundwater and subsurface drainage in humid 51
regions and potential for prolonged periods of no flow in arid channels. The flow of water within a river is 52
a key driver of landscape evolution, through the corresponding downstream force exerted on the stream bed, 53
the associated channel erosion, and the expression of local river incision at each elevation position along 54
the long profile. Therefore, we propose that the climate-streamflow relationship exerts a strong control on 55
long profiles. 56
Cimate is expressed differently in the downstream rate-of-change in streamflow between arid and humid 57
endmember rivers. In arid climates, streamflow tends to decrease downstream in all but extreme floods
22
58
for two main reasons: 1) Low annual rainfall, limited areal coverage of rainstorms, and short duration of 59
rainfall events generates partial area runoff
23
. This results in a small proportion of basin tributaries 60
contributing streamflow to the mainstem for limited periods of time. 2) Rivers are typically ephemeral (no 61
permanent flow)
24
, so channels lose water through dry, porous beds (transmission losses
22
) because water 62
tables are well below the channel
25
. Thus, the commonly assumed power law relationship between 63
streamflow and drainage area (with positive exponent c) breaks down
20
such that the long-term average 64
value of c may be negative, positive, or zero. In contrast, humid channels have perennial flow (all year 65
round), supported by baseflow from groundwater, and they accumulate flow from adjoining tributaries, 66
producing downstream increases in discharge
13
(positive c). We intuit that there is a spectrum of prevailing 67
downstream changes in streamflow across the globe based on the regional expression of climate within 68
discharge regimes (e.g., dryland hydrology, mountain front orography
5
), rather than simply on drainage 69
basin area. Given the obvious link between streamflow and riverbed erosion, we hypothesize that climatic 70
signatures are imprinted within river long profiles, superimposed upon other exogenous controls. In other 71
words, we expect a great deal of scatter typical of environmental data, but we hypothesize that climate will 72
reveal itself as a first-order control on long profile shape. 73
To test this hypothesis, we produced a new and unprecedented database of Global Longitudinal Profiles
74
(GLoPro) of rivers between 60°N and 56°S (Fig.1) extracted from NASA’s 30-m Shuttle Radar 75

Topography Mission Digital Elevation Model (SRTM-DEM)
26
. The profiles were extracted using 76
LSDTopoTools
27
, software with advanced capabilities in topographic analysis, employing a conservative 77
threshold for upstream drainage area and an algorithm of downstream flow accumulation, both of which 78
reduce the likelihood of Type 1 errors (Methods). For each profile we computed the Normalized Concavity 79
Index (NCI), a metric computed based solely on profile geometry (Methods; Extended Data Fig.1) that 80
allows for standardized comparisons of river profile shapes across the globe. The NCI is negative if the 81
profile is concave-up, zero if the profile is straight, and positive if the profile is convex-up. 82
We categorized each profile in GLoPro using the Köppen-Geiger (K-G) climate classification
28
and the 83
quantitative Aridity Index (AI = Precipitation/Potential Evapotranspiration)
29
, to investigate relationships 84
between climate and river long profile shape and to test whether the expression of aridity is detectable in 85
NCI. K-G is based on temperature and precipitation thresholds, emphasizing vegetation response to climate. 86
AI is a scale that represents the balance between precipitation and evaporative demand, and it declines with 87
aridity. Here we addressed the null hypothesis that there are no differences in NCI between climate 88
categories. We did not censor GLoPro for any other natural or anthropogenic factors, and it includes both 89
bedrock and alluvial rivers. We do not make any assumptions about whether the profiles are steady-state 90
(equilibrium) or transient, but we assumed that climate categories in K-G and AI have not changed 91
significantly over the timescales of long profile development (Methods). 92
The global distribution of NCI values does not suggest any strong geographic biases, although there are 93
clear concentrations of convex (Southern Siberia), concave (SE Asia), and nearly straight (Arabian 94
peninsula) rivers (Fig.1). NCI distributions of different climate classes (Fig.2a) overlap and display great 95
breadth, reflecting the large sample size and the many interacting independent variables (climate, tectonics, 96
lithology, and human factors) that affect drainage basin development. Nevertheless, statistically significant 97
differences between distributions are evident (Extended Data Fig.5). Comparing the four main K-G climate 98
zones, all NCI distributions are negatively skewed, revealing that river long profiles are generally 99
concave-up (Fig.2a). However, compared to the other three main climate zones (Tropical, Temperate, and 100

Citations
More filters
Journal ArticleDOI
TL;DR: The availability of planetary-scale geospatial datasets that can support the study of water-related disasters in the Anthropocene is rapidly growing as mentioned in this paper, and they review 124 global and free datasets allowing...
Abstract: The availability of planetary-scale geospatial datasets that can support the study of water-related disasters in the Anthropocene is rapidly growing. We review 124 global and free datasets allowing ...

43 citations

Journal ArticleDOI
TL;DR: A new framework for connecting flood changes to longitudinal variability in river conveyance, precipitation climatology, flows and sediment connectivity is proposed and the results show how this system of interacting units in the atmospheric, hydrologic and geomorphological realm function as a nonlinear filter that fundamentally alters the frequency of flood events.
Abstract: An improved understanding of changes in flood hazard and the underlying driving mechanisms is critical for predicting future changes for better adaptation strategies. While recent increases in flooding across the world have been partly attributed to a range of atmospheric or landscape drivers, one often-forgotten driver of changes in flood properties is the variability of river conveyance capacity. This paper proposes a new framework for connecting flood changes to longitudinal variability in river conveyance, precipitation climatology, flows and sediment connectivity. We present a first step, based on a regional analysis, towards a longer-term research effort that is required to decipher the circular causality between floods and rivers. The results show how this system of interacting units in the atmospheric, hydrologic and geomorphological realm function as a nonlinear filter that fundamentally alters the frequency of flood events. To revise and refine our estimation of future flood risk, this work highlights that multidriver attribution studies are needed, that include boundary conditions such as underlying climate, water and sediment connectivity, and explicit estimations of river conveyance properties.

41 citations

Journal ArticleDOI
TL;DR: In this paper, an hourly potential evapotranspiration (PET) dataset (hPET) was developed for the global land surface at 0.1° spatial resolution, based on output from the recently developed ERA5-Land reanalysis dataset, over the period 1981 to present.
Abstract: Challenges exist for assessing the impacts of climate and climate change on the hydrological cycle on local and regional scales, and in turn on water resources, food, energy, and natural hazards. Potential evapotranspiration (PET) represents atmospheric demand for water, which is required at high spatial and temporal resolutions to compute actual evapotranspiration and thus close the water balance near the land surface for many such applications, but there are currently no available high-resolution datasets of PET. Here we develop an hourly PET dataset (hPET) for the global land surface at 0.1° spatial resolution, based on output from the recently developed ERA5-Land reanalysis dataset, over the period 1981 to present. We show how hPET compares to other available global PET datasets, over common spatiotemporal resolutions and time frames, with respect to spatial patterns of climatology and seasonal variations for selected humid and arid locations across the globe. We provide the data for users to employ for multiple applications to explore diurnal and seasonal variations in evaporative demand for water.

34 citations

Journal ArticleDOI

31 citations


Cites background from "Aridity is expressed in river topog..."

  • ...These examples, along with similar recent work quantifying global fluvial geomorphic patterns (e.g., Chen et al., 2019; Frasson, Pavelsky, et al., 2019; Lin et al., 2020), suggest that RS is coming of age in its ability to provide global‐scale data that honors local differences in rivers....

    [...]

References
More filters
Journal ArticleDOI
TL;DR: In this article, the authors analyzed longitudinal profiles of bedrock streams in central Kentucky and of coastal plain streams in southeast Texas to determine the extent to which they exhibit smoothly concave profiles and to relate profile convexities to environmental controls.

78 citations

Journal ArticleDOI
TL;DR: A preliminary survey of some areas of the Central and Northern Namib Desert was undertaken in order to identify geomorphic and sedimentary features that could potentially provide information on past climatic change in the region.

77 citations

Journal ArticleDOI
TL;DR: In this paper, an analytic solution for non-steady state peak discharge was proposed to investigate how changes in storm duration alter hydrograph characteristics and the shape of the resulting equilibrium landscape, showing that decreasing relative storm duration results in downstreamdecreasing equilibrium channel concavity and decreasing valley density.
Abstract: [1] Landscape evolution models that calculate erosion as a function of discharge usually assume steady state runoff conditions and compute discharge as a power law function of the contributing area. This assumption is appropriate for small catchments and for regions in which the climate is characterized by long-lasting rainfall events. With larger catchments or shorter storms, however, the travel time of a water particle from the divide to the catchment outlet is typically longer than the characteristic storm duration. Hence a hydrologic steady state cannot be reached. This paper offers an analytic solution for nonsteady state peak discharge and investigates how changes in storm duration alter hydrograph characteristics and the shape of the resulting equilibrium landscape. An asymptotic function is used to incorporate the effect of storm duration into a landscape evolution model. We demonstrate with numerical simulations and analytical results that decreasing relative storm duration results in downstream-decreasing equilibrium channel concavity and decreasing valley density. It is also shown that this nonsteady state runoff system is less stable than the system of steady state runoff processes. Under certain conditions, no fixed dynamic equilibrium state exists. Furthermore, analysis of the shape of runoff hydrographs reveals feedback mechanisms manifested between hydrograph shape, erosion, and morphogenesis, demonstrating a compensatory mechanism between external precipitation input and internal runoff production.

76 citations

Journal ArticleDOI
TL;DR: In this article, the authors examine how a gradient of mean annual precipitation is expressed in the topography of low- to medium-relief, sediment-mantled, water-limited ecosystems.
Abstract: [1] Drainage density has long been observed to vary among climates, a relationship often attributable to differences in the erosivity of runoff and resistivity of vegetation. There is also evidence, though much less, that relief and channel concavity also vary with climate. The biophysical chain of events that connect climate to these topographic expressions, however, deserves greater attention. Using a numerical landscape evolution model, we examine how a gradient of mean annual precipitation is expressed in the topography of low- to medium-relief, sediment-mantled, water-limited ecosystems. We find equilibrium landscapes have the lowest drainage density, greatest relief, and lowest concavity at intermediate levels of rainfall. This climatic threshold represents the transition from vegetation-dominated sediment flux in drier climates to runoff-dominated flux with more precipitation. In drier climates, marginal increases in precipitation manifest themselves primarily as increases in vegetation while runoff is relatively constant; this acts to suppress sediment transport and decreases the drainage density. In wetter environments, marginal increases in precipitation lead to increases primarily in runoff while vegetation cover remains relatively constant; this results in increased sediment transport and increases the drainage density. The location of the transition depends at least in part on plant community structure and composition. The modeling study illustrates the complexities inherent in biogeomorphic systems but also that a simplified conceptual model of landscape evolution may indeed be sufficient to understand the large-scale patterns. The study also illustrates the opportunities offered by approaching questions of landform development from an ecohydrological perspective.

70 citations

Journal ArticleDOI
TL;DR: In this paper, the authors analyzed nearly one million digitally mapped river junctions throughout the contiguous United States and showed that branching angles vary systematically with climatic aridity, averaging roughly 45° in the driest places.
Abstract: Branching stream networks are a ubiquitous feature of the Earth's surface, but the processes that shape them, and their dependence on the climate in which they grow, remain poorly understood. Research has mainly focused on climatic controls of channel incision rates, while the climatic influence on planform geometry has often been overlooked. Here we analyze nearly one million digitally mapped river junctions throughout the contiguous United States and show that branching angles vary systematically with climatic aridity. In arid landscapes, which are thought to be dominated by surface runoff erosion, junction angles average roughly 45° in the driest places. Branching angles are systematically wider in humid regions, averaging roughly 72°, which is the theoretically predicted angle for network growth in a diffusive field such as groundwater seepage. The correlation of mean junction angle with aridity is stronger than with topographic gradient, downstream concavity, or other geometric factors that have been proposed as controls of junction angles. Thus, it may be possible to identify channelization processes from stream network geometry in relict landscapes, such as those on Mars.

69 citations

Trending Questions (2)
What is the relationship between topography and drainage?

The provided paper does not directly address the relationship between topography and drainage. It focuses on the relationship between river topography and aridity.

How does the status of river basins affect local and global climate patterns?

The provided paper does not directly address how the status of river basins affects local and global climate patterns.