Methane emissions from floodplains in the Amazon Basin: challenges in developing a process-based model for global applications

Tropical wetlands are estimated to represent about 50% of the natural wetland methane (CH4) emissions and explain a large fraction of the observed CH4 variability on timescales ranging from glacial–interglacial cycles to the currently observed year-to-year variability. Despite their importance, however, tropical wetlands are poorly represented in global models aiming to predict global CH4 emissions. This publication documents a first step in the development of a process-based model of CH4 emissions from tropical floodplains for global applications. For this purpose, the LPX-Bern Dynamic Global Vegetation Model (LPX hereafter) was slightly modified to represent floodplain hydrology, vegetation and associated CH4 emissions. The extent of tropical floodplains was prescribed using output from the spatially explicit hydrology model PCR-GLOBWB. We introduced new plant functional types (PFTs) that explicitly represent floodplain vegetation. The PFT parameterizations were evaluated against available remote-sensing data sets (GLC2000 land cover and MODIS Net Primary Productivity). Simulated CH4 flux densities were evaluated against field observations and regional flux inventories. Simulated CH4 emissions at Amazon Basin scale were compared to model simulations performed in the WETCHIMP intercomparison project. We found that LPX reproduces the average magnitude of observed net CH4 flux densities for the Amazon Basin. However, the model does not reproduce the variability between sites or between years within a site. Unfortunately, site information is too limited to attest or disprove some model features. At the Amazon Basin scale, our results underline the large uncertainty in the magnitude of wetland CH4 emissions. Sensitivity analyses gave insights into the main drivers of floodplain CH4 emission and their associated uncertainties. In particular, uncertainties in floodplain extent (i.e., difference between GLC2000 and PCR-GLOBWB output) modulate the simulated emissions by a factor of about 2. Our best estimates, using PCR-GLOBWB in combination with GLC2000, lead to simulated Amazon-integrated emissions of 44.4 ± 4.8 Tg yr−1. Additionally, the LPX emissions are highly sensitive to vegetation distribution. Two simulations with the same mean PFT cover, but different spatial distributions of grasslands within the basin, modulated emissions by about 20%. Correcting the LPX-simulated NPP using MODIS reduces the Amazon emissions by 11.3%. Finally, due to an intrinsic limitation of LPX to account for seasonality in floodplain extent, the model failed to reproduce the full dynamics in CH4 emissions but we proposed solutions to this issue. The interannual variability (IAV) of the emissions increases by 90% if the IAV in floodplain extent is accounted for, but still remains lower than in most of the WETCHIMP models. While our model includes more mechanisms specific to tropical floodplains, we were unable to reduce the uncertainty in the magnitude of wetland CH4 emissions of the Amazon Basin. Our results helped identify and prioritize directions towards more accurate estimates of tropical CH4 emissions, and they stress the need for more research to constrain floodplain CH4 emissions and their temporal variability, even before including other fundamental mechanisms such as floating macrophytes or lateral water fluxes.

Reconstruction of a flash flood with large wood transport and its influence on hazard patterns in an ungauged mountain basin

The reconstruction of past flash floods in ungauged basins leads to a high level of uncertainty, which increases if other processes are involved such as the transport of large wood material. An important flash flood occurred in 1997 in Venero Claro (Central Spain), causing significant economic losses. The wood material clogged bridge sections, raising the water level upstream. The aim of this study was to reconstruct this event, analysing the influence of woody debris transport on the flood hazard pattern. Because the reach in question was affected by backwater effects due to bridge clogging, using only high water mark or palaeostage indicators may overestimate discharges, and so other methods are required to estimate peak flows. Therefore, the peak discharge was estimated (123 ± 18 m3 s–1) using indirect methods, but one-dimensional hydraulic simulation was also used to validate these indirect estimates through an iterative process (127 ± 33 m3 s–1) and reconstruct the bridge obstruction to obtain the blockage ratio during the 1997 event (~48%) and the bridge clogging curves. Rainfall–Runoff modelling with stochastic simulation of different rainfall field configurations also helped to confirm that a peak discharge greater than 150 m3 s–1 is very unlikely to occur and that the estimated discharge range is consistent with the estimated rainfall amount (233 ± 27 mm). It was observed that the backwater effect due to the obstruction (water level ~7 m) made the 1997 flood (~35-year return period) equivalent to the 50-year flood. This allowed the equivalent return period to be defined as the recurrence interval of an event of specified magnitude, which, where large woody debris is present, is equivalent in water depth and extent of flooded area to a more extreme event of greater magnitude. These results highlight the need to include obstruction phenomena in flood hazard analysis.

A review of dendrogeomorphological research applied to flood risk analysis in Spain

Over the last forty years, applying dendrogeomorphology to palaeoflood analysis has improved estimates of the frequency and magnitude of past floods worldwide. This paper reviews the main results obtained by applying dendrogeomorphology to flood research in several case studies in Central Spain. These dendrogeomorphological studies focused on the following topics: (1) anatomical analysis to understand the physiological response of trees to flood damage and improve sampling efficiency; (2) compiling robust flood chronologies in ungauged mountain streams, (3) determining flow depth and estimating flood discharge using two-dimensional hydraulic modelling, and comparing them with other palaeostage indicators; (4) calibrating hydraulic model parameters (i.e. Manning roughness); and (5) implementing stochastic-based, cost–benefit analysis to select optimal mitigation measures. The progress made in these areas is presented with suggestions for further research to improve the applicability of dendrogeochronology to palaeoflood studies. Further developments will include new methods for better identification of the causes of specific types of flood damage to trees (e.g. tilted trees) or stable isotope analysis of tree rings to identify the climatic conditions associated with periods of increasing flood magnitude or frequency.

A new methodological protocol for the use of dendrogeomorphological data in flood risk analysis

Dendrogeomorphology uses information sources recorded in the roots, trunks and branches of trees and bushes located in the fluvial system to complement (or sometimes even replace) systematic and palaeohydrological records of past floods. The application of dendrogeomorphic data sources and methods to palaeoflood analysis over nearly 40 years has allowed improvements to be made in frequency and magnitude estimations of past floods. Nevertheless, research carried out so far has shown that the dendrogeomorphic indicators traditionally used (mainly scar evidence), and their use to infer frequency and magnitude, have been restricted to a small, limited set of applications. New possibilities with enormous potential remain unexplored. New insights in future research of palaeoflood frequency and magnitude using dendrogeomorphic data sources should: (1) test the application of isotopic indicators (16O/18O ratio) to discover the meteorological origin of past floods; (2) use different dendrogeomorphic indicators to estimate peak flows with 2D (and 3D) hydraulic models and study how they relate to other palaeostage indicators; (3) investigate improved calibration of 2D hydraulic model parameters (roughness); and (4) apply statistics-based cost–benefit analysis to select optimal mitigation measures. This paper presents an overview of these innovative methodologies, with a focus on their capabilities and limitations in the reconstruction of recent floods and palaeofloods.

Counting scars on tree stems to assess rockfall hazards: A low effort approach, but how reliable?

Rockfall is a widespread and hazardous process in mountain environments, but data on past events are only rarely available. Growth-ring series from trees impacted by rockfall were successfully used in the past to overcome the lack of archival records. Dendrogeomorphic techniques have been demonstrated to allow very accurate dating and reconstruction of spatial and temporal rockfall activity, but the approach has been cited to be labor intensive and time consuming. In this study, we present a simplified method to quantify rockfall processes on forested slopes requiring less time and efforts. The approach is based on a counting of visible scars on the stem surface of Common beech (Fagus sylvatica L.). Data are presented from a site in the Inn valley (Austria), where rocks are frequently detached from an ~ 200-m-high, south-facing limestone cliff. We compare results obtained from (i) the “classical” analysis of growth disturbances in the tree-ring series of 33 Norway spruces (Picea abies (L.) Karst.) and (ii) data obtained with a scar count on the stem surface of 50 F. sylvatica trees. A total of 277 rockfall events since A.D. 1819 could be reconstructed from tree-ring records of P. abies, whereas 1140 scars were observed on the stem surface of F. sylvatica. Absolute numbers of rockfalls (and hence return intervals) vary significantly between the approaches, and the mean number of rockfalls observed on the stem surface of F. sylvatica exceeds that of P. abies by a factor of 2.7. On the other hand, both methods yield comparable data on the spatial distribution of relative rockfall activity. Differences may be explained by a great portion of masked scars in P. abies and the conservation of signs of impacts on the stem of F. sylvatica. Besides, data indicate that several scars on the bark of F. sylvatica may stem from the same impact and thus lead to an overestimation of rockfall activity.

Adaptation options under climate change for multifunctional agriculture: a simulation study for western Switzerland

Besides its primary role in producing food and fiber, agriculture also has relevant effects on several other functions, such as management of renewable natural resources. Climate change (CC) may lead to new trade-offs between agricultural functions or aggravate existing ones, but suitable agricultural management may maintain or even improve the ability of agroecosystems to supply these functions. Hence, it is necessary to identify relevant drivers (e.g., cropping practices, local conditions) and their interactions, and how they affect agricultural functions in a changing climate. The goal of this study was to use a modeling framework to analyze the sensitivity of indicators of three important agricultural functions, namely crop yield (food and fiber production function), soil erosion (soil conservation function), and nutrient leaching (clean water provision function), to a wide range of agricultural practices for current and future climate conditions. In a two-step approach, cropping practices that explain high proportions of variance of the different indicators were first identified by an analysis of variance-based sensitivity analysis. Then, most suitable combinations of practices to achieve best performance with respect to each indicator were extracted, and trade-offs were analyzed. The procedure was applied to a region in western Switzerland, considering two different soil types to test the importance of local environmental constraints. Results show that the sensitivity of crop yield and soil erosion due to management is high, while nutrient leaching mostly depends on soil type. We found that the influence of most agricultural practices does not change significantly with CC; only irrigation becomes more relevant as a consequence of decreasing summer rainfall. Trade-offs were identified when focusing on best performances of each indicator separately, and these were amplified under CC. For adaptation to CC in the selected study region, conservation soil management and the use of cropped grasslands appear to be the most suitable options to avoid trade-offs.

Chronostratigraphy of the 600,000 year old continental record of Lake Van

Lake Van sediment cores from the Ahlat Ridge and Northern Basin drill sites of the ICDP project PALEOVAN contain a wealth of information about past environmental processes. The sedimentary sequence was dated using climatostratigraphic alignment, varve chronology, tephrostratigraphy, argon-argon single-crystal dating, radiocarbon dating, magnetostratigraphy, and cosmogenic nuclides. Based on the lithostratigraphic framework, the different age constraints are compiled and a robust and precise chronology of the 600,000 year-old Lake Van record is constructed. Proxy records of total organic carbon content and sediment color, together with the calcium/potassium-ratios and arboreal pollen percentages of the 174-meter-long Ahlat Ridge record, mimic the Greenland isotope stratotype (NGRIP). Therefore, the proxy records are systematically aligned to the onsets of interstadials reflected in the NGRIP or synthesized Greenland ice-core stratigraphy. The chronology is constructed using 27 age control points derived from visual synchronization with the GICC05 timescale, an absolutely-dated speleothem record (e.g., Hulu, Sanbao, Linzhu cave) and the Epica Dome C timescale. In addition, the uppermost part of the sequence is complemented with four ages from Holocene varve chronology and two calibrated radiocarbon ages. Furthermore, nine argon-argon ages and a comparison of the relative paleointensity record of the magnetic field with reference curve PISO-1500 confirm the accuracy of the age model. Also the identification of the Laschamp event via measurements of 10Be in the sediment confirms the presented age model. The chronology of the Ahlat Ridge record is transferred to the 79-meter-long event-corrected composite record from the Northern Basin and supplemented by additional radiocarbon dating on organic marco-remains. The basal age of the Northern Basin record is estimated at ~90 ka. The variations of the time series of total organic carbon content, the Ca/K ratio, and the arboreal pollen percentages illustrate that the presented chronology and paleoclimate data are suited for reconstructions and modeling of the Quaternary and Pleistocene climate evolution in the Near East at millennial timescales. Furthermore, the chronology of the last 250 kyr can be used to test other dating techniques.

Galilee Blooming: First Palynological and Archaeological Data from an Early Byzantine Cistern at Horvat Kur

Preliminary archaeological and palynological results are presented from an early Byzantine cistern of the village Horvat Kur in eastern Lower Galilee/Israel. The rural site was settled from the Hellenistic until the Early Arab period, its synagogue was constructed shortly after 425 AD and renovated sometimes during the 2nd half of the 6th century AD. It was abandoned probably as a consequence of the earthquake of 749 AD. The intact and properly sealed cistern contained complete or fully restorable pottery. Two cooking pots from the early 5th century AD comprised sediments which was sampled for palynological purposes. Both samples, as well as a sample from the soil beneath one of the pots and a modern surface sample from the site, revealed well preserved palynomorphs in comparably high concentration showing a great potential of the cistern as a pollen archive. The pollen content points to an open, grassy semiarid landscape with an apparent scarcity of cultivars and trees in the vicinity of the site and an abundance of herbs, especially Asteraceae, which are still commonly found in modern regional vegetation.

Antarctic Volcanic Flux Ratios from Law Dome Ice Cores

Explosive volcanic eruptions can inject large quantities of sulphur dioxide into the stratosphere. The aerosols that result from oxidation of the sulphur dioxide can produce significant cooling of the troposphere by reflecting or absorbing solar radiation. It is possible to obtain an estimate of the relative stratospheric sulphur aerosol concentration produced by different volcanoes by comparing sulphuric acid fluxes determined by analysis of polar ice cores. Here, we use a non-sea-salt sulphate time series derived from three well-dated Law Dome ice cores to investigate sulphuric acid flux ratios for major eruptions over the period AD 1301-1995. We use additional data from other cores to investigate systematic spatial variability in the ratios. Only for the Kuwae eruption (Law Dome ice date AD 1459.5) was the H2SO4 flux larger than that deposited by Tambora (Law Dome ice date AD 1816.7).

Late Pleistocene Interactions of East and West Antarctic Ice-Flow Regimes: Evidence from the McMurdo Ice Shelf

We present new interpretations of deglaciation in McMurdo Sound and the western Ross Sea, with observationally based reconstructions of interactions between East and West Antarctic ice at the last glacial maximum (LGM), 16 000, 12 000, 8000 and 4000 sp. At the LGM? East Antarctic ice from Mulock Glacier split, one branch turned westward south of Ross Island but the other branch rounded Ross Island before flowing southwest into McMurdo Sound. This flow regime, constrained by an ice saddle north of Ross Island, is consistent with the reconstruction of Stuiver and others (1981a). After the LGM, grounding-line retreat was most rapid in areas with greatest water depth, especially along the Victoria Land coast. By 12 000 sp, the ice-now regime in McMurdo Sound changed to through-flowing Mulock Glacier ice, with lesser contributions from Koettlitz, Blue and Ferrar Glaciers, because the former ice saddle north of Ross Island was replaced by a dome. The modern flew regime was established similar to 4000 BP. Ice derived from high elevations on the Polar Plateau but now stranded on the McMurdo Ice Shelf, and the pattern of the Transantarctic Mountains erratics support our reconstructions of Mulock Glacier ice rounding Minna Bluff but with all ice from Skelton Glacier ablating south of the bluff. They are inconsistent with Drewry's (1979) LGM reconstruction that includes Skelton Glacier ice in the McMurdo-Sound through-flow. Drewry's (1979) model closely approximates our results for 12 000-4000 BP. Ice-sheet modeling holds promise for determining whether deglaciation proceeded by grounding-line retreat of an ice sheet that was largely stagnant, because it never approached equilibrium flowline profiles after the Ross Ice Shelf, grounded, or of a dynamic ice sheet with flowline profiles kept low by active ice streams that extended northward from present-day outlet glaciers after the Ross Ice Shelf grounded.

Holtite and Dumortierite from the Szklary Pegmatite, Lower Silesia, Poland

The Szklary holtite is represented by three compositional varieties: (I) Ta-bearing (up to 14.66 wt.% Ta(2)O(5)), which forms homogeneous crystals and cores within zoned crystals; (2) Ti-bearing (up to 3.82 wt.% TiO(2)), found as small domains within the core; and (3) Nb-bearing (up to 5.30 wt.% Nb(2)O(5),) forming the rims of zoned crystals. All three varieties show variable Sb+As content, reaching 19.18 wt.% Sb(2)O(3) (0.87 Sb a.p.f.u.) and 3.30 wt.% As(2)O(3) (0.22 As a.p.f.u.) in zoned Ta-bearing holtite, which constitutes the largest Sb+As content reported for the mineral. The zoning in holtite is a result of Ta-Nb fractionation in the parental pegmatite-forming melt together with contamination of the relatively thin Szklary dyke by Fe, Mg and Ti. Holtite and the As- and Sb-bearing dumortierite, which in places overgrows the youngest Nb-bearing zone, suggest the following crystallization sequence: Ta-bearing holtite -> Ti-bearing holtite -> Nb-bearing holtite -> As- and Sb-bearing, (Ta,Nb,Ti)-poor dumortierite -> As- and Sb-dominant, (Ta,Nb,Ti)-free dumortierite-like mineral (16.81 wt.% As(2)O(3) and 10.23 wt.% Sb(2)O(3)) with (As+Sb) > Si. The last phase is potentially a new mineral species, Al(6)rectangle B(Sb,As)(3)O(15). or Al(5)rectangle(2)B(Sb,As)(3)O(12)(OH)(3), belonging to the dumortierite group. The Szklary holtite shows no evidence of clustering of compositions around 'holtite I' and 'holtite II'. Instead, the substitutions of Si(4+) by Sb(3+)+As(3+) at the Si/Sb sites and of Ta(5+) by Nb(5+) or Ti(4+) at the Al(l) site suggest possible solid solutions between: (1) (Sb,As)-poor and (Sb,As)-rich holtite; (2) dumortierite and the unnamed (As+Sb)-dominant dumortierite-like mineral; and (3) Ti-bearing dumortierite and holtite, i.e. our data provide further evidence for miscibility between holtite and dumortierite, but leave open the question of defining the distinction between them. The Szklary holtite crystallized from the melt along with other primary Ta-Nb-(Ti) minerals such as columbite-(Mn), tantalite-(Mn), stibiotantalite and stibiocolumbite as the availability of Ta decreased. The origin of the parental melt can be related to anatexis in the adjacent Sowie Mountains complex, leading to widespread migmatization and metamorphic segregation in pelitic-psammitic sediments metamorphosed at similar to 390-380 Ma.

Stable-Isotope and Trace Element Time Series from Fedchenko Glacier (Pamirs) Snow/Firn Cores

In summer 2005, two pilot snow/firn cores were obtained at 5365 and 5206 m a.s.l. on Fedchenko glacier, Pamirs, Tajikistan, the world's longest and deepest alpine glacier. The well-defined seasonal layering appearing in stable-isotope and trace element distribution identified the physical links controlling the climate and aerosol concentration signals. Air temperature and humidity/precipitation were the primary determinants of stable-isotope ratios. Most precipitation over the Pamirs originated in the Atlantic. In summer, water vapor was re-evaporated from semi-arid regions in central Eurasia. The semi-arid regions contribute to non-soluble aerosol loading in snow accumulated on Fedchenko glacier. In the Pamir core, concentrations of rare earth elements, major and other elements were less than those in the Tien Shan but greater than those in Antarctica, Greenland, the Alps and the Altai. The content of heavy metals in the Fedchenko cores is 2-14 times lower than in the Altai glaciers. Loess from Afghan-Tajik deposits is the predominant lithogenic material transported to the Pamirs. Trace elements generally showed that aerosol concentration tended to increase on the windward slopes during dust storms but tended to decrease with altitude under clear conditions. The trace element profile documented one of the most severe droughts in the 20th century.

CO(2) Diffusion in Polar Ice: Observations from Naturally Formed CO(2) Spikes in the Siple Dome (Antarctica) Ice Core

One common assumption in interpreting ice-core CO(2) records is that diffusion in the ice does not affect the concentration profile. However, this assumption remains untested because the extremely small CO(2) diffusion coefficient in ice has not been accurately determined in the laboratory. In this study we take advantage of high levels of CO(2) associated with refrozen layers in an ice core from Siple Dome, Antarctica, to study CO(2) diffusion rates. We use noble gases (Xe/Ar and Kr/Ar), electrical conductivity and Ca(2+) ion concentrations to show that substantial CO(2) diffusion may occur in ice on timescales of thousands of years. We estimate the permeation coefficient for CO(2) in ice is similar to 4 x 10(-21) mol m(-1) s(-1) Pa(-1) at -23 degrees C in the top 287 m (corresponding to 2.74 kyr). Smoothing of the CO(2) record by diffusion at this depth/age is one or two orders of magnitude smaller than the smoothing in the firn. However, simulations for depths of similar to 930-950m (similar to 60-70 kyr) indicate that smoothing of the CO(2) record by diffusion in deep ice is comparable to smoothing in the firn. Other types of diffusion (e.g. via liquid in ice grain boundaries or veins) may also be important but their influence has not been quantified.

A 12,000 Year Record of Explosive Volcanism in the Siple Dome Ice Core, West Antarctica

Air mass trajectories in the Southern Hemisphere provide a mechanism for transport to and deposition of volcanic products on the Antarctic ice sheet from local volcanoes and from tropical and subtropical volcanic centers. This study extends the detailed record of Antarctic, South American, and equatorial volcanism over the last 12,000 years using continuous glaciochemical series developed from the Siple Dome A (SDMA) ice core, West Antarctica. The largest volcanic sulfate spike ( 280 mu g/L) occurs at 5881 B. C. E. Other large signals with unknown sources are observed around 325 B. C. E. ( 270 mu g/L) and 2818 B. C. E. ( 191 mu g/L). Ages of several large equatorial or Southern Hemisphere volcanic eruptions are synchronous with many sulfate peaks detected in the SDMA volcanic ice chemistry record. The microprobe "fingerprinting'' of glass shards in the SDMA core points to the following Antarctic volcanic centers as sources of tephra found in the SDMA core: Balenny Island, Pleiades, Mount Berlin, Mount Takahe, and Mount Melbourne as well as Mount Hudson and possibly Mount Burney volcanoes of South America. Identified volcanic sources provide an insight into the poorly resolved transport history of volcanic products from source volcanoes to the West Antarctic ice sheet.