Sedimentology of the Vietnam shelf

On the Vietnam Shelf more than 1000 miles of shallow high-resolution seismics were analyzed to unravel post-glacial evolution in a tropical, siliciclastic environment together with 25 sediment cores from water depths between 21 and 169 m to determine stratigraphy, distribution and style of sedimentation. Fourty-seven samples were dated with the AMS-14C technique. The shelf was grouped into three regions: a southern part, a central part, and a northern part. On the broad Southern Shelf, sedimentation is influenced by the Mekong River, which drains into the SCS in this area. Here, incised valley fills are abundant that were cut into the late Pleistocene land surface by the Paleo-Mekong River during times of sea level lowstand. Those valleys are filled with transgressive deposits. The Holocene sedimentation rate in this low gradient accommodation-dominated depositional system is in the range of 5-10 and 25-40 cm/ky at locations sheltered from currents. The Central Shelf is narrow and the sedimentary strata are conformable. Here, numerous small mountainous rivers reach the SCS and transport large amounts of detrital sediment onto the shelf. Therefore, the Holocene sedimentation rate is high with values of 50-100 cm/ky in this supply-dominated depositional system. The broad Northern Shelf in the vicinity of the Red River Delta shows, as on the Southern Shelf, incised valleys cut into the Pleistocene land surface by paleo river channels. In this accommodation-dominated shelf area, the sedimentation rate is low with values of 5-10 cm/ky. Where applicable, we assigned the sampled deposits to different paleo-facies. The latter are related to certain intervals of water depths at their time of deposition. Comparison with the sea-level curve of (Hanebuth et al., 2000, doi:10.1126/science.288.5468.1033) indicates subsidence on the Central Shelf, which is in agreement with the high sedimentation rates in this area. In contrast, data from the Northern Shelf suggest tectonic uplift that might be related to recent tectonic movements along the Ailao Shan-Red River Fault zone. Data from the Southern Shelf are generally in agreement with the sea-level curve mentioned above.

Establishment, expansion and regeneration of aspen (Populus tremula) in the Torneträsk area, subarctic Sweden

In subarctic Sweden, recent decadal colonization and expansion of aspen (Populus tremula L.) were recorded. Over the past 100 years, aspen became c. 16 times more abundant, mainly as a result of increased sexual regeneration. Moreover, aspen now reach tree-size (>2 m) at the alpine treeline, an ecotone that has been dominated by mountain birch (Betula pubescens ssp. czerepanovii) for at least the past 4000 years. We found that sexual regeneration in aspen probably occurred seven times or more within the last century. Whereas sexual regeneration occurred during moist years following a year with an exceptionally high June-July temperature, asexual regeneration was favored by warm and dry summers. Disturbance to the birch forest by cyclic moth population outbreaks was critical in aspen establishment in the subalpine area. At the treeline, aspen colonization was less determined by these moth outbreaks, and was mainly restricted by summer temperature. If summer warming persists, aspen spread may continue in subarctic Sweden, particularly at the treeline. However, changing disturbance regimes, future herbivore population dynamics and the responses of aspen's competitors birch and pine to a changing climate may result in different outcomes.

Global 2009 1km MODIS (MOD35/MOD09) cloud frequency and MOD35 processing path

Identifying cloud interference in satellite-derived data is a critical step toward developing useful remotely sensed products. Most MODIS land products use a combination of the MODIS (MOD35) cloud mask and the 'internal' cloud mask of the surface reflectance product (MOD09) to mask clouds, but there has been little discussion of how these masks differ globally. We calculated global mean cloud frequency for both products, for 2009, and found that inflated proportions of observations were flagged as cloudy in the Collection 5 MOD35 product. These erroneously categorized areas were spatially and environmentally non-random and usually occurred over high-albedo land-cover types (such as grassland and savanna) in several regions around the world. Additionally, we found that spatial variability in the processing path applied in the Collection 5 MOD35 algorithm affects the likelihood of a cloudy observation by up to 20% in some areas. These factors result in abrupt transitions in recorded cloud frequency across landcover and processing-path boundaries impeding their use for fine-scale spatially contiguous modeling applications. We show that together, these artifacts have resulted in significantly decreased and spatially biased data availability for Collection 5 MOD35-derived composite MODIS land products such as land surface temperature (MOD11) and net primary productivity (MOD17). Finally, we compare our results to mean cloud frequency in the new Collection 6 MOD35 product, and find that landcover artifacts have been reduced but not eliminated. Collection 6 thus increases data availability for some regions and land cover types in MOD35-derived products but practitioners need to consider how the remaining artifacts might affect their analysis.

(Tables 1-3) Water chemistry of cloud forest streams at baseflow conditions, Rio San Francisco, Ecuador

We investigated controls on the water chemistry of a South Ecuadorian cloud forest catchment which is partly pristine, and partly converted to extensive pasture. From April 2007 to May 2008 water samples were taken weekly to biweekly at nine different subcatchments, and were screened for differences in electric conductivity, pH, anion, as well as element composition. A principal component analysis was conducted to reduce dimensionality of the data set and define major factors explaining variation in the data. Three main factors were isolated by a subset of 10 elements (Ca2+, Ce, Gd, K+, Mg2+, Na+, Nd, Rb, Sr, Y), explaining around 90% of the data variation. Land-use was the major factor controlling and changing water chemistry of the subcatchments. A second factor was associated with the concentration of rare earth elements in water, presumably highlighting other anthropogenic influences such as gravel excavation or road construction. Around 12% of the variation was explained by the third component, which was defined by the occurrence of Rb and K and represents the influence of vegetation dynamics on element accumulation and wash-out. Comparison of base- and fast flow concentrations led to the assumption that a significant portion of soil water from around 30 cm depth contributes to storm flow, as revealed by increased rare earth element concentrations in fast flow samples. Our findings demonstrate the utility of multi-tracer principal component analysis to study tropical headwater streams, and emphasize the need for effective land management in cloud forest catchments.

Half-degree gridded nitrogen and phosphorus fertilizer use for global agriculture production during 1900-2013

In addition to enhance agricultural productivity, synthetic nitrogen (N) and phosphorous (P) fertilizer application in croplands dramatically altered global nutrient budget, water quality, greenhouse gas balance, and their feedbacks to the climate system. However, due to the lack of geospatial fertilizer input data, current Earth system/land surface modeling studies have to ignore or use over-simplified data (e.g., static, spatially uniform fertilizer use) to characterize agricultural N and P input over decadal or century-long period. We therefore develop a global time-series gridded data of annual synthetic N and P fertilizer use rate in croplands, matched with HYDE 3,2 historical land use maps, at a resolution of 0.5º latitude by longitude during 1900-2013. Our data indicate N and P fertilizer use rates increased by approximately 8 times and 3 times, respectively, since the year 1961, when IFA (International Fertilizer Industry Association) and FAO (Food and Agricultural Organization) survey of country-level fertilizer input were available. Considering cropland expansion, increase of total fertilizer consumption amount is even larger. Hotspots of agricultural N fertilizer use shifted from the U.S. and Western Europe in the 1960s to East Asia in the early 21st century. P fertilizer input show the similar pattern with additional hotspot in Brazil. We find a global increase of fertilizer N/P ratio by 0.8 g N/g P per decade (p< 0.05) during 1961-2013, which may have important global implication of human impacts on agroecosystem functions in the long run. Our data can serve as one of critical input drivers for regional and global assessment on agricultural productivity, crop yield, agriculture-derived greenhouse gas balance, global nutrient budget, land-to-aquatic nutrient loss, and ecosystem feedback to the climate system.

222Radon flux map for Europe in netCDF format

A high-resolution 222Radon (222Rn) flux map for Europe was developed, based on a parameterization of 222Rn production and transport in the soil. The 222Rn exhalation rate is parameterized based on soil properties, uranium content, and modelled soil moisture from two different land-surface reanalysis data sets. Spatial variations in exhalation rates are primarily determined by the uranium content of the soil, but also influenced by soil texture and local water table depth. Temporal variations are related to soil moisture variations as the molecular diffusion in the unsaturated soil zone depends on available air-filled pore space. Monthly 222Rn exhalation rates from European soils were calculated with a nominal spatial resolution of 0.083° x 0.083°. The two realizations of the 222Rn flux map, based on the different soil moisture data sets, both realistically reproduce the observed seasonality in the fluxes but yield considerable differences for absolute flux values. The mean 222Rn flux from soils in Europe is estimated to be 10 mBq/m**2/s (ERA-Interim/Land soil moisture) or 15 mBq/m**2/s (GLDAS-Noah soil moisture) for the period 2006-2010. The 222Rn flux maps for Europe are available for the application in atmospheric transport studies, e.g to evaluate the performance of atmospheric transport models.

Stable carbon isotope ratios and accumulation rates of organic carbon and nutrients of DSDP Hole 79-545 (Appendix A)

Future warming is predicted to shift the Earth system into a mode with progressive increase and vigour of extreme climate events possibly stimulating other mechanisms that invigorate global warming. This study provides new data and modelling investigating climatic consequences and biogeochemical feedbacks that happened in a warmer world ~112 Myr ago. Our study focuses on the Cretaceous Oceanic Anoxic Event (OAE) 1b and explores how the Earth system responded to a moderate ~25,000 yr lasting climate perturbation that is modelled to be less than 1 °C in global average temperature. Using a new chronological model for OAE 1b we present high-resolution elemental and bulk carbon isotope records from DSDP Site 545 from Mazagan Plateau off NW Africa and combine this information with a coupled atmosphere-land-ocean model. The simulations suggest that a perturbation at the onset of OAE 1b caused almost instantaneous warming of the atmosphere on the order of 0.3 °C followed by a longer (~45,000 yr) period of ~0.8 °C cooling. The marine records from DSDP Site 545 support that these moderate swings in global climate had immediate consequences for African continental supply of mineral matter and nutrients (phosphorous), subsequent oxygen availability, and organic carbon burial in the eastern subtropical Atlantic, however, without turning the ocean anoxic. The match between modelling results and stratigraphic isotopic data support previous studies [summarized in Jenkyns 2003, doi:10.1098/rsta.2003.1240] in that methane emission from marine hydrates, albeit moderate in dimension, may have been the trigger for OAE 1b, though we can not finally rule out alternative mechanisms. Following the hydrate mechanism a total of 1.15 * 10**18 g methane carbon (delta13C=-60 ?), equivalent to about 10% to the total modern gas hydrate inventory, generated the delta13Ccarb profile recorded in the section. Modelling suggests a combination of moderate-scale methane pulses supplemented by continuous methane emission at elevated levels over ~25,000 yr. The proposed mechanism, though difficult to finally confirm in the geological past, is arguably more likely to occur in a warmer world and apparently perturbs global climate and ocean chemistry almost instantaneously. This study shows that, once set-off, this mechanism can maintain Earth's climate in a perturbed mode over geological time leading to pronounced changes in regional climate.

Model-based changes in global annual mean surface temperature change (Delta T_g) and radiative forcing due to land ice albedo changes (Delta R_[LI]) over the last 5 Myr, supplementary material

It is still an open question how equilibrium warming in response to increasing radiative forcing - the specific equilibrium climate sensitivity S - depends on background climate. We here present palaeodata-based evidence on the state dependency of S, by using CO2 proxy data together with a 3-D ice-sheet-model-based reconstruction of land ice albedo over the last 5 million years (Myr). We find that the land ice albedo forcing depends non-linearly on the background climate, while any non-linearity of CO2 radiative forcing depends on the CO2 data set used. This non-linearity has not, so far, been accounted for in similar approaches due to previously more simplistic approximations, in which land ice albedo radiative forcing was a linear function of sea level change. The latitudinal dependency of ice-sheet area changes is important for the non-linearity between land ice albedo and sea level. In our set-up, in which the radiative forcing of CO2 and of the land ice albedo (LI) is combined, we find a state dependence in the calculated specific equilibrium climate sensitivity, S[CO2,LI], for most of the Pleistocene (last 2.1 Myr). During Pleistocene intermediate glaciated climates and interglacial periods, S[CO2,LI] is on average ~ 45 % larger than during Pleistocene full glacial conditions. In the Pliocene part of our analysis (2.6-5 Myr BP) the CO2 data uncertainties prevent a well-supported calculation for S[CO2,LI], but our analysis suggests that during times without a large land ice area in the Northern Hemisphere (e.g. before 2.82 Myr BP), the specific equilibrium climate sensitivity, S[CO2,LI], was smaller than during interglacials of the Pleistocene. We thus find support for a previously proposed state change in the climate system with the widespread appearance of northern hemispheric ice sheets. This study points for the first time to a so far overlooked non-linearity in the land ice albedo radiative forcing, which is important for similar palaeodata-based approaches to calculate climate sensitivity. However, the implications of this study for a suggested warming under CO2 doubling are not yet entirely clear since the details of necessary corrections for other slow feedbacks are not fully known and the uncertainties that exist in the ice-sheet simulations and global temperature reconstructions are large.

Static GPS measurements in the vicinity of EDML and an elevation model for Dronning Maud Land

The knowledge of ice sheet surface topography and the location of the ice divides are essential for ice dynamic modeling. An improved digital elevation model (DEM) of Dronning Maud Land (DML), Antarctica, is presented in this paper. It is based on ground-based kinematic GPS profiles, airborne radar altimetry, and data of the airborne radio-echo sounding system, as well as spaceborne laser altimetry from NASA's Ice, Cloud and land Elevation Satellite (ICESat). The accuracy of ICESat ice sheet altimetry data in the area of investigation is discussed. The location of the ice divides is derived from aspect calculation of the topography and is verified with several velocity data derived from repeated static GPS measurements.

(Table 1) Surface water chemistry of northern Victoria Land lakes

Concentrations of major ions, silicate and nutrients (total N and P) were measured in samples of surface water from 28 lakes in ice-free areas of northern Victoria Land (East Antarctica). Sixteen lakes were sampled during austral summers 2001/02, 2003/04, 2004/05 and 2005/06 to assess temporal variation in water chemistry. Although samples showed a wide range in ion concentrations, their composition mainly reflected that of seawater. In general, as the distance from the sea increased, the input of elements from the marine environment (through aerosols and seabirds) decreased and there was an increase in nitrate and sulfate concentrations. Antarctic lakes lack outflows and during the austral summer the melting and/or ablation of ice cover, water evaporation and leaching processes in dry soils determine a progressive increase in water ion concentrations. During the five-year monitoring survey, no statistically significant variation in the water chemistry were detected, except for a slight (hardly significant) increase in TN concentrations. However, Canonical Correspondence Analysis (CCA) indicated that other factors besides distance from the sea, the presence of nesting seabirds, the sampling time and percentage of ice cover affect the composition of water in Antarctic cold desert environments.