Fluvial incision and tectonic uplift across the Himalayas of central Nepal

The pattern of fluvial incision across the Himalayas of central Nepal is estimated from the distribution of Holocene and Pleistocene terraces and from the geometry of modern channels along major rivers draining across the range. The terraces provide good constraints on incision rates across the Himalayan frontal folds (Sub-Himalaya or Siwaliks Hills) where rivers are forced to cut down into rising anticlines and have abandoned numerous strath terraces. Farther north and upstream, in the Lesser Himalaya, prominent fill terraces were deposited, probably during the late Pleistocene, and were subsequently incised. The amount of bedrock incision beneath the fill deposits is generally small, suggesting a slow rate of fluvial incision in the Lesser Himalaya. The terrace record is lost in the high range where the rivers are cutting steep gorges. To complement the terrace study, fluvial incision was also estimated from the modern channel geometries using an estimate of the shear stress exerted by the flowing water at the bottom of the channel as a proxy for river incision rate. This approach allows quantification of the effect of variations in channel slope, width, and discharge on the incision rate of a river; the determination of incision rates requires an additional lithological calibration. The two approaches are shown to yield consistent results when applied to the same reach or if incision profiles along nearby parallel reaches are compared. In the Sub-Himalaya, river incision is rapid, with values up to 10–15 mm/yr. It does not exceed a few millimeters per year in the Lesser Himalaya, and rises abruptly at the front of the high range to reach values of ∼4–8 mm/yr within a 50-km-wide zone that coincides with the position of the highest Himalayan peaks. Sediment yield derived from the measurement of suspended load in Himalayan rivers suggests that fluvial incision drives hillslope denudation of the landscape at the scale of the whole range. The observed pattern of erosion is found to closely mimic uplift as predicted by a mechanical model taking into account erosion and slip along the flat-ramp-flat geometry of the Main Himalayan Thrust fault. The morphology of the range reflects a dynamic equilibrium between present-day tectonics and surface processes. The sharp relief together with the high uplift rates in the Higher Himalaya reflects thrusting over the midcrustal ramp rather than the isostatic response to reincision of the Tibetan Plateau driven by late Cenozoic climate change, or late Miocene reactivation of the Main Central Thrust

Active folding of fluvial terraces across the Siwaliks Hills, Himalayas of central Nepal

We analyze geomorphic evidence of recent crustal deformation in the sub-Himalaya of central Nepal, south of the Kathmandu Basin. The Main Frontal Thrust fault (MFT), which marks the southern edge of the sub-Himalayan fold belt, is the only active structure in that area. Active fault bend folding at the MFT is quantified from structural geology and fluvial terraces along the Bagmati and Bakeya Rivers. Two major and two minor strath terraces are recognized and dated to be 9.2, 2.2, and 6.2, 3.7 calibrated (cal) kyr old, respectively. Rock uplift of up to 1.5 cm/yr is derived from river incision, accounting for sedimentation in the Gangetic plain and channel geometry changes. Rock uplift profiles are found to correlate with bedding dip angles, as expected in fault bend folding. It implies that thrusting along the MFT has absorbed 21 ± 1.5 mm/yr of N-S shortening on average over the Holocene period. The ±1.5 mm/yr defines the 68% confidence interval and accounts for uncertainties in age, elevation measurements, initial geometry of the deformed terraces, and seismic cycle. At the longitude of Kathmandu, localized thrusting along the Main Frontal Thrust fault must absorb most of the shortening across the Himalaya. By contrast, microseismicity and geodetic monitoring over the last decade suggest that interseismic strain is accumulating beneath the High Himalaya, 50–100 km north of the active fold zone, where the Main Himalayan Thrust (MHT) fault roots into a ductile décollement beneath southern Tibet. In the interseismic period the MHT is locked, and elastic deformation accumulates until being released by large (M_w > 8) earthquakes. These earthquakes break the MHT up to the near surface at the front of the Himalayan foothills and result in incremental activation of the MFT

Investigation of the relationships between basin morphology, tectonic uplift, and denudation from the study of an active fold belt in the Siwalik Hills, central Nepal

The present study investigates correlations between an extensive range of geomorphic properties that can be estimated from a digital elevation model and the uplift rate on geological timescales. The analysis focuses on an area in the Siwalik Hills (central Nepal), where lithology and climate can be considered as uniform. This area undergoes rapid tectonic uplift at rates of up to 15 mm yr^(−1), which are derived from the geometric pattern of a fault-bend model of fold growth. The selected geomorphic properties can be divided in two categories, depending on whether or not the vertical dimension is taken into account. None of the planar properties are significantly correlated to uplift rate, unlike those that include the vertical dimension, such as the mean elevation of basins, hypsometric curve, and hypsometric integral, and relief defined by the amplitude factor of length scaling analysis. Correlation between relief and uplift rate is observed for all length scales of topography shorter than 600 m, which suggests that all orders of the streams are able to adjust to the tectonic signal. Simple mass balance considerations imply that the average elevation is only 10% of surface uplift, suggesting that a dynamic equilibrium has been reached quite rapidly. Using a simple two-process model for erosion, we find that fairly high diffusion coefficients (order of 10 m^2 yr^(−1)) and efficient transport of the material by rivers are required. This unusually high value for mass diffusivity at small length scales may be obtained by either a very efficient linear diffusion or by landsliding. Actually, both processes may be active, which appears likely given the nature of the unconsolidated substratum and the favorable climatic conditions. Local relief in the study area may therefore be used to predict either uplift or denudation, but the prediction is calibrated only for that specific climatic and lithologic conditions and cannot be systematically applied to other contexts

Persistence of full glacial conditions in the central Pacific until 15,000 years ago

The magnitude of atmospheric cooling during the Last Glacial Maximum and the timing of the transition into the current interglacial period remain poorly constrained in tropical regions, partly because of a lack of suitable climate records. Glacial moraines provide a method of reconstructing past temperatures, but they are relatively rare in the tropics. Here we present a reconstruction of atmospheric temperatures in the central Pacific during the last deglaciation on the basis of cosmogenic ^3He ages of moraines and numerical modelling of the ice cap on Mauna Kea volcano, Hawaii—the only highland in the central Pacific on which moraines that formed during the last glacial period are preserved. Our reconstruction indicates that the Last Glacial Maximum occurred between 19,000 and 16,000 years ago in this region and that temperatures at high elevations were about 7 °C lower than today during this interval. Glacial retreat began about 16,000 years ago, but temperatures were still about 6.5 °C lower than today until 15,000 years ago. When combined with estimates of sea surface temperatures in the central Pacific Ocean, our reconstruction indicates that the lapse rate during the Last Glacial Maximum was higher than at present, which is consistent with the proposal that the atmosphere was drier at that time. Furthermore, the persistence of full glacial conditions until 15,000 years ago is consistent with the relatively late and abrupt transition to warmer temperatures in Greenland5, indicating that there may have been an atmospheric teleconnection between the central Pacific and North Atlantic regions during the last deglaciation

Interseismic Strain Accumulation on the Himalayan Crustal Ramp (Nepal)

The Departement of Mines and Geology has been monitoring the seismicity of the Central Himalayas of Nepal since 1985. Intense microseismicity and frequent medium‐size earthquakes (mL<4) tend to cluster beneath the topographic front of the Higher Himalaya. This 10–20km deep seismicity also correlates with a zone of localized uplift that has been evidenced from geodetic data. Both microseismic and geodetic data indicate strain accumulation on a mid‐crustal ramp that had been previously inferred from geological and geophysical evidence. This ramp connects a flat decollement under the Lesser and Sub‐Himalaya with a deeper decollement under the Higher Himalaya, and probably acts as a geometric asperity where strain and stress build up during the interseismic period. The large Himalayan earthquakes could nucleate there and probably activate the whole flat‐and‐ramp system up to the blind thrusts of the Sub‐Himalaya

Testing fluvial erosion models using the transient response of bedrock rivers to tectonic forcing in the Apennines, Italy

The transient response of bedrock rivers to a drop in base level can be used to discriminate between competing fluvial erosion models. However, some recent studies of bedrock erosion conclude that transient river long profiles can be approximately characterized by a transport‐limited erosion model, while other authors suggest that a detachment‐limited model best explains their field data. The difference is thought to be due to the relative volume of sediment being fluxed through the fluvial system. Using a pragmatic approach, we address this debate by testing the ability of end‐member fluvial erosion models to reproduce the well‐documented evolution of three catchments in the central Apennines (Italy) which have been perturbed to various extents by an independently constrained increase in relative uplift rate. The transport‐limited model is unable to account for the catchments’response to the increase in uplift rate, consistent with the observed low rates of sediment supply to the channels. Instead, a detachment‐limited model with a threshold corresponding to the field‐derived median grain size of the sediment plus a slope‐dependent channel width satisfactorily reproduces the overall convex long profiles along the studied rivers. Importantly, we find that the prefactor in the hydraulic scaling relationship is uplift dependent, leading to landscapes responding faster the higher the uplift rate, consistent with field observations. We conclude that a slope‐ dependent channel width and an entrainment/erosion threshold are necessary ingredients when modeling landscape evolution or mapping the distribution of fluvial erosion rates in areas where the rate of sediment supply to channels is low

Kinematics of fault-related folding derived from a sandbox experiment

We analyze the kinematics of fault tip folding at the front of a fold-and-thrust wedge using a sandbox experiment. The analog model consists of sand layers intercalated with low-friction glass bead layers, deposited in a glass-sided experimental device and with a total thickness h = 4.8 cm. A computerized mobile backstop induces progressive horizontal shortening of the sand layers and therefore thrust fault propagation. Active deformation at the tip of the forward propagating basal décollement is monitored along the cross section with a high-resolution CCD camera, and the displacement field between pairs of images is measured from the optical flow technique. In the early stage, when cumulative shortening is less than about h/10, slip along the décollement tapers gradually to zero and the displacement gradient is absorbed by distributed deformation of the overlying medium. In this stage of detachment tip folding, horizontal displacements decrease linearly with distance toward the foreland. Vertical displacements reflect a nearly symmetrical mode of folding, with displacements varying linearly between relatively well defined axial surfaces. When the cumulative slip on the décollement exceeds about h/10, deformation tends to localize on a few discrete shear bands at the front of the system, until shortening exceeds h/8 and deformation gets fully localized on a single emergent frontal ramp. The fault geometry subsequently evolves to a sigmoid shape and the hanging wall deforms by simple shear as it overthrusts the flat ramp system. As long as strain localization is not fully established, the sand layers experience a combination of thickening and horizontal shortening, which induces gradual limb rotation. The observed kinematics can be reduced to simple analytical expressions that can be used to restore fault tip folds, relate finite deformation to incremental folding, and derive shortening rates from deformed geomorphic markers or growth strata

A Comparison of ACE Measurements of Galactic Cosmic-Ray Abundances and Energy Spectra for Two Successive Solar Minima

Using current solar minimum measurements from the Cosmic Ray Isotope Spectrometer (CRIS) onboard the Advanced Composition Explorer (ACE), we report the observed elemental abundances and energy spectra for C, O, Si, and Fe in the energy range of ~50-500MeV/nucleon. These measurements are compared to prior CRIS observations for the 1997-98 solar minimum period, when the solar magnetic field was of the opposite polarity. By April 2009, the current solar minimum intensities for each of the four elements has surpassed the peak intensities observed during the previous solar minimum. We also examine the correlation of the galactic cosmic ray intensities with the tilt angle of the heliospheric current sheet. These studies will be important for understanding solar modulation processes, such as changes in drift effects during a solar cycle

A novel topographic parameterization scheme indicates that martian gullies display the signature of liquid water

Martian gullies resemble gullies carved by water on Earth, yet are thought to have formed in an extremely cold (2-driven processes. That this argument persists demonstrates the limitations of morphological interpretations made from 2D images, especially when similar-looking landforms can form by very different processes. To overcome this we have devised a parameterization scheme, based on statistical discriminant analysis and hydrological terrain analysis of meter-scale digital topography data, which can distinguish between dry and wet surface processes acting on a landscape. Applying this approach to new meter-scale topographic datasets of Earth, the Moon and Mars, we demonstrate that martian gullied slopes are dissimilar to dry, gullied slopes on Earth and the Moon, but are similar to both terrestrial debris flows and fluvial gullies. We conclude that liquid water was integral to the process by which martian gullies formed. Finally, our work shows that quantitative 3D analyses of landscape have great potential as a tool in planetary science, enabling remote assessment of processes acting on planetary surfaces

Controls and limits on bedrock channel geometry

Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/94837/1/jgrf728.pd