Rapid crustal uplift in Patagonia due to enhanced ice loss

A vertical crustal uplift rate of 39 mm yr^? 1 is measured between 2003 and 2006 using Global Positioning System (GPS) measurements at the northeastern edge of the Southern Patagonia Icefield (SPI). This is the largest present-day glacial isostatic rate ever recorded. The combination of SPI's rapid melting and the unique regional slab-window tectonics that promotes a relatively low viscosity, is central to our interpretation of the observations. The two effects lead to a strong interaction of short relaxation times with ice loads that change on a comparable time scale. The profile of GPS observations link ice loss to the soft viscoelastic isostatic flow response over the time scale of the Little Ice Age (LIA), including ice loss in the period of observation. The agreement of the results with our model predictions strongly suggests the large crustal uplift in Patagonia is due an accelerated glacier wasting since the termination of the LIA and that the effective regional mantle viscosity is near 4.0–8.0 × 10¹⁸ Pa s. A century-long diminution of the icefields, at rates that are about 1/4 – 1/2 the contemporary loss rates, is consistent with multidecadal-scale temperature trends estimated for the past 50–100 years and is, in fact, a key feature in any model capable of explaining the uplift data.     

Streamflow variation due to glacier melting and climate change in upstream Heihe River Basin,Northwest China

Streamflow simulation is often challenging in mountainous watersheds because of incomplete hydrological models, irregular topography, immeasurable snowpack or glacier, and low data resolution. In this study, a semi-distributed conceptual hydrological model (SWAT-Soil Water Assessment Tool) coupled with a glacier melting algorithm was applied to investigate the sensitivity of streamflow to climatic and glacial changes in the upstream Heihe River Basin. The glacier mass balance was calculated at daily time-step using a distributed temperature-index melting and accumulation algorithm embedded in the SWAT model. Specifically, the model was calibrated and validated using daily streamflow data measured at Yingluoxia Hydrological Station and decadal ice volume changes derived from survey maps and remote sensing images between 1960 and 2010. This study highlights the effects of glacier melting on streamflow and their future changes in the mountainous watersheds. We simulate the contribution of glacier melting to streamflow change under different scenarios of climate changes in terms of temperature and precipitation dynamics. The rising temperature positively contributed to streamflow due to the increase of snowmelt and glacier melting. The rising precipitation directly contributes to streamflow and it contributed more to streamflow than the rising temperature. The results show that glacial meltwater has contributed about 3.25 billion m³ to streamflow during 1960–2010. However, the depth of runoff within the watershed increased by about 2.3 mm due to the release of water from glacial storage to supply the intensified evapotranspiration and infiltration. The simulation results indicate that the glacier made about 8.9% contribution to streamflow in 2010. The research approach used in this study is feasible to estimate the glacial contribution to streamflow in other similar mountainous watersheds elsewhere.     

Latitudinal patterns of export production recorded in surface sediments of the Chilean Patagonian fjords (41–55°S) as a response to water column productivity

The Chilean Patagonian fjords region (41–56°S) is characterized by highly complex geomorphology and hydrographic conditions, and strong seasonal and latitudinal patterns in precipitation, freshwater discharge, glacier coverage, and light regime; all of these directly affect biological production in the water column. In this study, we compiled published and new information on water column properties (primary production, nutrients) and surface sediment characteristics (biogenic opal, organic carbon, molar C/N, bulk sedimentary δ¹³C_org) from the Chilean Patagonian fjords between 41°S and 55°S, describing herein the latitudinal pattern of water column productivity and its imprint in the underlying sediments. Based on information collected at 188 water column and 118 sediment sampling sites, we grouped the Chilean fjords into four main zones: Inner Sea of Chiloé (41° to ～44°S), Northern Patagonia (44° to ～47°S), Central Patagonia (48–51°S), and Southern Patagonia (Magellan Strait region between 52° and 55°S). Primary production in the Chilean Patagonian fjords was the highest in spring–summer, reflecting the seasonal pattern of water column productivity. A clear north–south latitudinal pattern in primary production was observed, with the highest average spring and summer estimates in the Inner Sea of Chiloé (2427 and 5860 mg C m^?2 d^?1) and Northern Patagonia (1667 and 2616 mg C m^?2 d^?1). This pattern was closely related to the higher availability of nutrients, greater solar radiation, and extended photoperiod during the productive season in these two zones. The lowest spring value was found in Caleta Tortel, Central Patagonia (91 mg C m^?2 d^?1), a site heavily influenced by glacier meltwater and river discharge loaded with glacial sediments. Biogenic opal, an important constituent of the Chilean fjord surface sediments (Si_OPAL ～1–13%), reproduced the general north–south pattern of primary production and was directly related to water column silicic acid concentrations. Surface sediments were also rich in organic carbon content and the highest values corresponded to locations far away from glacier influence, sites within fjords, and/or semi-enclosed and protected basins, reflecting both autochthonous (water column productivity) and allochthonous sources (contribution of terrestrial organic matter from fluvial input to the fjords). A gradient was observed from the more oceanic sites to the fjord heads (west–east) in terms of bulk sedimentary δ¹³C_org and C/N ratios; the more depleted (δ¹³C_org ?26‰) and higher C/N (23) values corresponded to areas close to rivers and glaciers. A comparison of the Chilean Patagonian fjords with other fjord systems in the world revealed high variability in primary production for all fjord systems as well as similar surface sediment geochemistry due to the mixing of marine and terrestrial organic carbon.     

Glacier changes and their impacts on the discharge in the past half-century in Tekes watershed,Central Asia

The glacier is an important and stable water supply in Central Asia. Monitoring the change of glacier and understanding the impacts of glacier change on river discharge are critical to predict the downstream water availability change in future. Glacier changes were discussed and their impacts on river discharge were evaluated by hydrological modeling with a distributed hydrological model SWAT under two land use and land cover scenarios (1970 and 2007) in Tekes watershed, the most important source of water discharge to the Ili River. Compared to the glacier area of 1511 km² in 1970s it decreased by 332 km² in 2007, which resulted in the contribution the discharge from precipitation in the glacier area to the average annual discharge of the watershed changing from 9.8% in the period 1966–1975 to 7.8% in the period 2000–2008. In the month scale, with the decrease of glacier area, the distribution of the contribution of monthly discharge from precipitation in the glacier area to the total of the watershed changed from bimodal pattern to unimodal pattern. By linking a hydrological model to remote sensing image analysis and Chinese glacier inventories to determine glacier area change our approach in quantifying the impacts of glacier changes on hydrology at different scales, will provide quantitative information for stakeholders in making decisions for water resource management.     

The radiocarbon reservoir effect in proglacial lakes: Examples from Antarctica

We studied the apparently old radiocarbon ages from lakes in the dry valleys of Antarctica. The radiocarbon reservoir effect in these lakes results from two components: the inherited age and the residence age. The inherited age is derived from input of old carbon, primarily from subsurface melt of adjacent glaciers. The residence age comes from in situ aging of lake water in an environment sealed from the atmosphere. Our results indicate that surface melt of glaciers introduces little ancient carbon to the lake system, because of rapid gas equilibration with the atmosphere. Subsurface melt in lakes with large glacier cross-sectional areas at the grounding line, however, can contribute a significant amount of ancient carbon, leading to lake-bottom reservoir effects in excess of ∼ 2700 yr. This value can increase to ∼ 20,000 yr immediately at the grounding line. In most lakes, however, surface melt far exceeds that from the subsurface and dilutes the effect of ancient carbon, making the inherited age relatively low. Residence ages generally are on the order of a few thousand years, but can be as much as ∼ 10,000 yr. Because a residence age is reset when the lake loses its ice cover and is exposed to wind-driven mixing, its magnitude can provide important information about lake history.     

Stability of massive ground ice bodies in University Valley,McMurdo Dry Valleys of Antarctica: Using stable O–H isotope as tracers of sublimation in hyper-arid regions

To date, studies of the stability of subsurface ice in the McMurdo Dry Valleys of Antarctica have been mainly based on climate-based vapor diffusion models. In University Valley (1800 m), a small glacier is found at the base of the head of the valley, and adjacent to the glacier, a buried body of massive ice was uncovered beneath 20–40 cm of loose cryotic sediments and sandstone boulders. This study assesses the origin and stability of the buried body of massive ice by measuring the geochemistry and stable O–H isotope composition of the ice and applies a sublimation and molecular diffusion model that accounts for the observed trends. The results indicate that the buried massive ice body represents an extension of the adjacent glacier that was buried by a rock avalanche during a cold climate period. The contrasting δ¹⁸O profiles and regression slope values between the uppermost 6 cm of the buried massive ice (upward convex δ¹⁸O profile and S_D-18O = 5.1) and that below it (progressive increase in δ¹⁸O and S_D-18O = 6.4) suggest independent post-depositional processes affected the isotope composition of the ice. The upward convex δ¹⁸O profile in the uppermost 6 cm is consistent with the ice undergoing sublimation. Using a sublimation and molecular diffusion model, and assuming that diffusion occurred through solid ice, the sublimation rate needed to fit the measured δ¹⁸O profile is 0.2 ? 10^? 3 mm yr^? 1, a value that is more similar to net ice removal rates derived from ³He data from cobbles in Beacon Valley till (7.0 ? 10^? 3 mm yr^? 1) than sublimation rates computed based on current climate (0.1–0.2 mm yr^?1). We suggest that the climate-based sublimation rates are offset due to potential ice recharge mechanisms or to missing parameters, particularly the nature and thermo-physical properties of the overlying sediments (i.e., temperature, humidity, pore structure and ice content, grain size).     

Regional geoid modeling in the area of subglacial Lake Vostok,Antarctica

We present a geoid model for the area of Lake Vostok, Antarctica, from a combination of local airborne gravity, ice-surface and ice-thickness data and a lake bathymetry model. The topography data are used for residual terrain modeling (RTM) in a remove–restore approach together with GOCE satellite data. The quasigeoid is predicted by least-squares collocation (LSC) and subsequently converted to geoid heights. Special aspects of that method in presence of an ice sheet are discussed.It is well known that a body freely floating in water is in a state of hydrostatic equilibrium (HE). This usually applies, e.g., to ice shelves or sea ice. However, it has been shown that this is valid also for the ice sheet covering the subglacial Lake Vostok. Thus, we demonstrate the use of such a refined regional geoid model for glaciological and geophysical applications by means of the HE surface of that lake. The mean quadratic residual geoid signal (0.56 m) w.r.t. the GOCE background model exceeds the residual variations of the estimated apparent lake level (ALL) (0.26 m) within the central part of the lake. An approach considering the actual geopotential at the ALL has been derived and subsequently applied. In this context, downward continuation of the potential field within the ice sheet as well as the latitudinal tilt of off-geoid equipotential surfaces are discussed. In view of the accuracy of the ice-thickness measurements that dominate the total error budget of the estimated ALL these effects are negligible. Thus, the HE surface of subglacial lakes may safely be described by a constant height bias in small-scale regional applications. However, field continuation is significant with respect to the formal uncertainty of the quasigeoid, which is at the level of 5 cm given that accurate airborne gravity data (±2 mGal) are available.     

Spatial and temporal variations of albedo on nine glaciers in western China from 2000 to 2011

This research demonstrates the spatiotemporal variations of albedo on nine glaciers in western China during 2000–2011, by the albedo derived from two types of datasets: Landsat TM/ETM + images and MOD10A1 product. Then, the influence factors of glacier albedo and its relationship with glacier mass balance are also analyzed by the correlation approach, which is frequently used in geostatistics. The paper finds that there are different spatiotemporal variations over the glaciers in western China: (1) For a single glacier, the albedo varies gently with altitude on its tongue and increases fast in the middle part, while in the accumulation zones, the albedo value appears in the form of fluctuation. This could provide a quantitative method to retrieve the snowline by determining the threshold albedo value of snowpack and bare ice. (2) For the glaciers in western China, the albedo decreases with distance to the center of Tibetan Plateau (TP). This may relate to the elevation of glacier, for the speed of glacier retreat highly depends on air temperature. (3) In the summer period, albedo on most glaciers declines over the last 12 years, and it decreases much faster in southeastern TP than other regions, for which air temperature overwhelms the black carbon concentration. In addition, the trend of glacier albedo in summer is greatly correlated with that of measured glacier mass balance, which implies that the long‐term albedo datasets by remote sensing technology could be used to monitor and predict the change of glacier mass balance in the future. Copyright © 2013 John Wiley & Sons, Ltd.     

Volume and mass changes of the Greenland ice sheet inferred from ICESat and GRACE

This study examines the recent evolution of the Greenland ice sheet and its six major drainage basins. Based on laser altimetry data acquired by the Ice, Cloud and Land Elevation Satellite (ICESat), covering the period September–November 2003 to February–March 2008, ice surface height changes and their temporal variations were inferred. Our refined repeat track analysis is solely based on ICESat data and is independent of external elevation models, since it accounts for both ice height changes and the local topography. From the high resolution ice height change pattern we infer an overall mean surface height trend of −0.12 ± 0.006 m yr⁻¹. Furthermore, the largest changes could be identified at coastal margins of the ice sheet, exhibiting rates of more than −2 m yr⁻¹. The total ice volume change of the entire ice sheet amounts to −205.4 ± 10.6 km³ yr⁻¹. In addition, we assessed mass changes from 78 monthly Gravity Recovery and Climate Experiment (GRACE) solutions. The Release-04 gravity field solutions of GeoForschungsZentrum Potsdam cover the period between August 2002 and June 2009. We applied an adjusted regional integration approach in order to minimize the leakage effects. Attention was paid to an optimized filtering which reduces error effects from different sources. The overall error assessment accounts for GRACE errors as well as for errors due to imperfect model reductions. In particular, errors caused by uncertainties in the glacial isostatic adjustment models could be identified as the largest source of errors. Finally, we determined both seasonal and long-term mass change rates. The latter amounts to an overall ice mass change of −191.2 ± 20.9 Gt yr⁻¹ corresponding to 0.53 ± 0.06 mm yr⁻¹ equivalent eustatic sea level rise. From the combination of the volume and mass change estimates we determined a mean density of the lost mass to be 930 ± 11 kg m⁻³. This value supports our applied density assumption 900 ± 30 kg m⁻³ which was used to perform the volume–mass-conversion of our ICESat results. Hence, mass change estimates from two independent observation techniques were inferred and are generally in good agreement.     

Two extraterrestrial dust horizons found in the Dome Fuji ice core,East Antarctica

Two silicate-rich dust layers were found in the Dome Fuji ice core in East Antarctica, at Marine Isotope Stages 12 and 13. Morphologies, textures, and chemical compositions of constituent particles reveal that they are high-temperature melting products and are of extraterrestrial origin. Because similar layers were found ～ 2000 km east of Dome Fuji, at EPICA (European Project for Ice Coring in Antarctica)-Dome C, particles must have rained down over a wide area 434 and 481 ka. The strewn fields occurred over an area of at least 3 × 10⁶ km². Chemical compositions of constituent phases and oxygen isotopic composition of olivines suggest that the upper dust layer was produced by a high-temperature interaction between silicate-rich melt and water vapor due to an impact explosion or an aerial burst of a chondritic meteoroid on the inland East Antarctic ice sheet. An estimated total mass of the impactor, on the basis of particle flux and distribution area, is at least 3 × 10⁹ kg. A possible parent material of the lower dust layer is a fragment of friable primitive asteroid or comet. A hypervelocity impact of asteroidal/cometary material on the upper atmosphere and an explosion might have produced aggregates of sub-μm to μm-sized spherules. Total mass of the parent material of the lower layer must exceed 1 × 10⁹ kg. The two extraterrestrial horizons, each a few millimeters in thickness, represent regional or global meteoritic events not identified previously in the Southern Hemisphere.     

Crustal velocity structure in the southern edge of the Central Alborz (Iran)

Inversion of local earthquake travel times and joint inversion of receiver functions and Rayleigh wave group velocity measurements were used to derive a simple model for the velocity crustal structure beneath the southern edge of the Central Alborz (Iran), including the seismically active area around the megacity of Tehran. The P and S travel times from 115 well-located earthquakes recorded by a dense local seismic network, operated from June to November 2006, were inverted to determine a 1D velocity model of the upper crust. The limited range of earthquake depths (between 2 km and 26 km) prevents us determining any velocity interfaces deeper than 25 km. The velocity of the lower crust and the depth of the Moho were found by joint inversion of receiver functions and Rayleigh wave group velocity data. The resulting P-wave velocity model comprises an upper crust with 3 km and 4 km thick sedimentary layers with P wave velocities (V_p) of ～5.4 and ～5.8 km s^?1, respectively, above 9 km and 8 km thick layers of upper crystalline crust (V_p ～6.1 and ～6.25 km s^?1 respectively). The lower crystalline crust is ～34 km thick (V_p ～ 6.40 km s^?1). The total crustal thickness beneath this part of the Central Alborz is 58 ± 2 km.     

The active fault system of SW Alps

Historical and active seismicity in the south-western Alps (France and Italy) shows the recurrence of relatively high-magnitude earthquakes (M ≥ 5.8), like the one that recently affected the Italian Apennine range (M = 6.3 on the 30th March 2009). However, up-to-date detailed mapping of the active fault network has been poorly established. The evaluation of seismological hazard in particular in the highly populated French and Italian coastal region cannot be done without this. Here, we present a detailed study of the main active fault system, based on geological observations along the south-western flank of the Alpine arc. This N140° right-lateral strike-slip active fault system runs along the edge of the Argentera-Mercantour range and can be followed down to the Mediterranean Sea. It is evidenced by (1) Holocene offsets of glacial geomorphology witnessing ongoing fault activity since 10 ka, (2) widespread recent (10–20 Ma) pseudotachylytes featuring long term activity of the faults, (3) active landslides along the main fault zone, (4) geothermal anomalies (hot springs) emerging in the active faults, (5) ongoing low-magnitude seismic activity and (6) localization of the main historical events. In the light of our investigations, we propose a new tectonic pattern for the active fault system in the south-western Alps.     

Resolving sea level contributions by identifying fingerprints in time-variable gravity and altimetry

We have studied the ability of the GRACE gravimetry mission and Jason-1 altimetry to resolve ice and glacier induced contributions to sea level rise, by means of a fingerprint method. Here, the signals from ice sheet and land glacier changes, steric changes, glacial isostatic adjustment and terrestrial hydrology are assumed to have fixed spatial patterns. In a joint inversion using GRACE and Jason-1 data the unknown temporal components can then be estimated by least-squares. In total, we estimate temporal components for up to ∼ 80 individual patterns. From a propagation of the full error-covariance from GRACE and a diagonal error-covariance from Jason-1 altimetry we find that: (1) GRACE almost entirely explains the mass related parameters in the joint inversion, (2) an inversion using only Jason-1 data has a marginal ability to estimate the mass related parameters, while the steric parameters have much better formal accuracy. In terms of mean sea level rise the steric patterns have a maximum formal accuracy of 0.01 mm for an 11 week running mean. In general, strong negative error correlations (ρ < − 0.9) exists between the high and low elevation parts of the ice sheet drainage basins, when those are estimated independently. The largest formal errors found are in the order of 40 Gton for small high elevation subbasins in the southern Greenland ice sheet, which are difficult to separate. In a simplified joint inversion, merging high and low elevation basins, we have investigated the ability of the GRACE and Jason-1 data to separate the geocenter motion into a present-day contribution and a contribution from glacial isostatic adjustment (GIA). We find that the GIA related signal is larger than the present-day component with a maximum of −0.71 mm/year in the Z direction. Total geocenter motion rates are found to be −0.28, 0.43, −1.08 mm/year for the X, Y and Z components, respectively. The inversion results have been propagated to the Jason-1 along-track measurements. Over the time period considered, we see that a large part of the variability in the Pacific, Atlantic and Indian ocean can be explained by our inversion results. The applied inversion method therefore seems a feasible way to separate steric from mass induced sea level changes. At the same time, the joint inversion would benefit from more advanced parameterizations, which may aid in fitting remaining signal from altimetry.     

A reassessment of current volcanic emissions from the Central American arc with specific examples from Nicaragua

The Central American volcanic arc supplies a significant proportion of the persistent annual global sulphur dioxide emissions from volcanoes. In November/December 2003, we completed a survey of the arc section from Mombacho to San Cristóbal in Nicaragua recording individual mean fluxes of 800, 530 and 220 Mg day^− 1 in the plumes from San Cristóbal, Telica and Masaya, respectively. An assessment of fluxes published since 1997 along the entire Central America arc yields a mean total arc flux of SO₂ of 4360 Mg day^− 1 or 8–16% of the annual estimated global volcanic SO₂ flux to the troposphere. New field data shows that Masaya volcano continues to show stable HCl/SO₂ and HF/SO₂ ratios, suggesting a sustained flux of these components of ∼ 220 and 30 Mg day^− 1, respectively (1997 to 2004). Masaya's plume composition also appears to have been stable, between 2001 and 2003, with respect to all the particulate species measured, with significant fluxes of SO₄²⁻ (4 Mg day^− 1), Na⁺ (0.9–1.3 Mg day^− 1) and K⁺ (0.7 Mg day^− 1). Extrapolating the Masaya plume species ratios to the entire Central American arc gives mean HCl and HF fluxes of 1300 and 170 Mg day^− 1 and a particulate sulphate flux of 40 Mg day^− 1 for 1997 to 2004, although without further understanding of the degassing processes and sources at depth of these different volatiles, these arc-scale estimates should be treated with caution. Combining our arc scale mean SO₂ flux with published measurements of volcanic gas compositions with respect to CO₂ and H₂O allows us to estimate mean CO₂ fluxes of 4400–9600 Mg day^− 1 and H₂O fluxes of 70,000–78,000 Mg day^− 1 for the arc. Preliminary comparisons of these estimates of outgassing rates with published volatile input fluxes into the Central American subduction zone, suggest that Cl is more efficiently recycled through the subduction zone than CO₂. The results for H₂O are inconclusive.     

Young eclogite from the Greater Himalayan Sequence,Arun Valley,eastern Nepal: P–T–t path and tectonic implications

Garnet geochronology was used to provide the first direct measurement of the timing of eclogitization in the central Himalaya. Lu–Hf dates from garnet separates in one relict eclogite from the Arun River Valley in eastern Nepal indicate an age of 20.7 ± 0.4 Ma, significantly younger than ultra-high pressure eclogites from the western Himalaya, reflecting either different origins or substantial time lags in tectonics along strike. Four proximal garnet amphibolites from structurally lower horizons are 14–15 Ma, similar to post-eclogitization ages published for rocks along strike in southern Tibet. P–T calculations indicate three metamorphic episodes for the eclogite: i) eclogite-facies metamorphism at ～ 670 °C and ≥ 15 kbar at 23–16 Ma; ii) a peak-T granulite event at ～ 780 °C and 12 kbar; and iii) late-stage amphibolite-facies metamorphism at ～ 675 °C and 6 kbar at ～ 14 Ma. The garnet amphibolites were metamorphosed at ～ 660 °C. Three models are considered to explain the observed P–T–t evolution. The first assumes that the Main Himalayan Thrust (basal thrust of the Himalayan thrust system) cuts deeper at Arun than elsewhere. While conceptually the simplest, this model has difficulty explaining both the granulite-facies overprint and the pulse of exhumation between 25 and 14 Ma. A second model assumes that (aborted) subduction, slab breakoff, and ascent of India's leading edge occurred diachronously: ～ 50 Ma in the western Himalaya, ～ 25 Ma in the central Himalaya of Nepal, and presumably later in the eastern Himalaya. This model explains the P–T–t path, particularly heating during initial exhumation, but implies significant along-strike diachroneity, which is generally lacking in other features of the Himalaya. A third model assumes repeated loss of mantle lithosphere, first by slab breakoff at ～ 50 Ma, and again by delamination at ～ 25 Ma; this model explains the P–T–t path, but requires geographically restricted tectonic behavior at Arun. The P–T–t history of the Arun eclogites may imply a change in the physical state of the Himalayan metamorphic wedge at 16–25 Ma, ultimately giving rise to the Main Central Thrust by 15–16 Ma.     

Distribution of dissolved labile and particulate iron and copper in Terra Nova Bay polynya (Ross Sea,Antarctica) surface waters in relation to nutrients and phytoplankton growth

The distribution of the dissolved labile and of the particulate Fe and Cu together with dissolved oxygen, nutrients, chlorophyll a and total particulate matter was investigated in the surface waters of Terra Nova Bay polynya in mid-January 2003. The measurements were conducted within the framework of the Italian Climatic Long-term Interactions of the Mass balance in Antarctica (CLIMA) Project activities. The labile dissolved fraction was operationally defined by employing the chelating resin Chelex-100, which retains free and loosely bound trace metal species. The dissolved labile Fe ranges from below the detection limit (0.15 nM) to 3.71 nM, while the dissolved labile Cu from below the detection limit (0.10 nM) to 0.90 nM. The lowest concentrations for both metals were observed at 20 m depth (the shallowest depth for which metals were measured). The concentration of the particulate Fe was about 5 times higher than the dissolved Fe concentration, ranging from 0.56 to 24.83 nM with an average of 6.45 nM. The concentration of the particulate Cu ranged from 0.01 to 0.71 nM with an average of 0.17 nM. The values are in agreement with the previous data collected in the same area. We evaluated the role of the Fe and Cu as biolimiting metals. The N:dissolved labile Fe ratios (18,900–130,666) would or would not allow a complete nitrate removal, on the basis of the N:Fe requirement ratios that we calculated considering the N:P and the C:P ratios estimated for diatoms. This finding partially agrees with the Si:N ratio that we found (2.29). Moreover we considered a possible influence of the dissolved labile Cu on the Fe uptake process.     

Factors controlling sediment yield in a major South American drainage basin: the Magdalena River,Colombia

The Magdalena River, a major fluvial system draining most of the Colombian Andes, has the highest sediment yield of any medium-sized or large river in South America. We examined sediment yield and its response to control variables in the Magdalena drainage basin based on a multi-year dataset of sediment loads from 32 tributary catchments. Various morphometric, hydrologic, and climatic variables were estimated in order to understand and predict the variation in sediment yield. Sediment yield varies from 128 to 2200 t km⁻² yr⁻¹ for catchments ranging from 320 to 59,600 km². The mean sediment yield for 32 sub-basins within the Magdalena basin is ∼690 t km⁻² yr⁻¹. Mean annual runoff is the dominant control and explains 51% of the observed variance in sediment yield. A multiple regression model, including two control variables, runoff and maximum water discharge, explains 58% of the variance. This model is efficient (ME=0.89) and is a valuable tool for predicting total sediment yield from tributary catchments in the Magdalena basin. Multiple correlations for those basins corresponding to the upper Magdalena, middle basin, Eastern Cordillera, and catchment areas greater than 2000 km², explain 75, 77, 89, and 78% of the variance in sediment yield, respectively. Although more variance is explained when dataset are grouped into categories, the models are less efficient (ME<0.72). Within the spatially distributed models, six catchment variables predict sediment yield, including runoff, precipitation, precipitation peakedness, mean elevation, mean water discharge, and relief. These estimators are related to the relative importance of climate and weathering, hillslope erosion, and fluvial transport processes. Time series analysis indicates that significant increases in sediment load have occurred over 68% of the catchment area, while 31% have experienced a decreasing trend in sediment load and thus yield. Land use analysis and increasing sediment load trends indicate that erosion within the catchment has increased over the last 10–20 years.     

40Ar / 39Ar geochronology constraints on late miocene weathering rates in Minas Gerais,Brazil

⁴⁰Ar / ³⁹Ar incremental heating ages for twenty one grains of cryptomelane, collected at 0, 42, 45, and 60 m depths in the Cachoeira Mine weathering profile, Minas Gerais, permit calculating long-term (10 Ma time scale) weathering rate (saprolitization rate) in SE Brazil. Pure well-crystallized cryptomelane grains with high K contents (3–5 wt.%) yield reliable geochronological results. The ⁴⁰Ar / ³⁹Ar plateau ages obtained decrease from the top to the bottom of the profile (12.7 ± 0.1 to 7.6 ± 0.1 Ma at surface; 7.6  ± 0.2 to 6.1 ± 0.2 Ma at 42 m; and 7.1 ± 0.2 to 5.9 ± 0.1 Ma at 45 m; 6.6 ± 0.1 to 5.2 ± 0.1 Ma at 60 m), yielding a weathering front propagation rate of 8.9 ± 1.1 m/m.y. From the geochronological results and the mineral transformations implicit by the current mineralogy in the weathering profiles, it is possible to calculate the saprolitization rate for the Cachoeira Mine lithologies and for adjacent weathering profiles developed on granodiorites and schists. The measured weathering front propagation rate yields a saprolitization rate of 24.9 ± 3.1 t/km²/yr. This average long-term (> 10 Ma) saprolitization rate is consistent with mass balance calculations results for present saprolitization rates in weathering watersheds. These results are also consistent with long-term saprolitization rates estimated by combining cosmogenic isotope denudation rates with mass balance calculations.     

The formation of supraglacial debris covers by primary dispersal from transverse englacial debris bands

Glaciological controls on debris cover formation are investigated from the perspective of primary dispersal of supraglacial debris across a melting ice surface. This involves the migration of angled debris septa outcrops across a melting, thinning glacier ablation zone. Three measures of a glacier's ability to evacuate supraglacial debris are outlined: (1) a concentration factor describing the focusing of englacial debris into specific supraglacial mass loads; (2) the rate of migration of a septum outcrop relative to the local ice surface; and (3) a downstream velocity differential between a slower septum outcrop and the faster ice surface velocity. Measures (1) and (2) are inversely related, while measure (3) increases down‐glacier to explain why slow‐moving, thinning ice rapidly becomes debris covered. Data from Glacier d'Estelette (Italian Alps) are used to illustrate these processes, and to explore the potential for debris cover formation and growth in different glaciological environments. The transition from a ‘clean’, transport‐dominated to a debris‐covered ablation‐dominated glacier is explained by the melting out of more closely‐spaced debris septa, in combination with the geometric interactions of angled septa and ice surface in a field of reducing flow and increasing ablation. The growth and shrinkage of debris covers are most sensitive to glaciological changes at glaciers with gently‐dipping debris‐bearing foliation, but less sensitive at high‐compression glaciers whose termini are constrained by moraine dams and other forms of obstruction. These findings show that a variety of debris‐covered glacier types will show a spectrum of response characteristics to negative mass balance. Copyright © 2013 John Wiley & Sons, Ltd.     

U–Pb dating of speleothems from Spannagel Cave,Austrian Alps: A high resolution comparison with U–series ages

U–Pb dating is increasingly used to date speleothems that are too old for precise U–Th disequilibrium dating; however there is little data that can independently validate its application to such material. This study presents U–Pb ages for speleothems from the Spannagel Cave in the Austrian Alps including a detailed comparison with U–Th ages from an unusually U–rich sample that yields precise ages by both methods. Sample SPA4 is a flowstone with three growth phases separated by distinct hiatuses. For the youngest growth phase the U–Pb and U–Th ages are 267 ± 1 ka and 267 ± 5 ka respectively; the middle growth phase is 291 ± 1 versus 295 ± 11 ka while for the oldest growth phase a single sub-sample, assuming the same initial Pb composition as for the younger phases, yields an age of 340 ± 2 ka compared to 353 ± 9 ka by U–Th. Correlation of these ages with the marine isotope stages confirms that these speleothems grew during glacial stages as suggested by previous work on the same sample. Sample SPA 15 has U–Th isotopic compositions indistinguishable from secular equilibrium; the U–Pb data on the main growth phase of this sample give an age of 551 ± 10 ka, whereas a single analysis from the oldest phase suggests it may be on the order of 40 ka older. This detailed comparison of U–Pb and U–Th ages provides important support for the potential validity of the U–Pb method in older samples beyond the range of U–Th.