Reply to comment of Legates et al.

In the previous comment, Legates et al. express concern about the statistical reliability of the positive runoff–temperature relationship presented by Labat et al. We are grateful for this opportunity to respond to these concerns. As Legates et al. correctly points out, the effect of temperature on runoff is a complex relationship, which involves precipitation, evaporation, anthropomorphic affects, among others. As such, the effect of increased temperature on runoff is strongly dependent on the identity of the watershed of interest. For example, a watershed located in a glaciated region, such as Iceland, exhibits a strong positive correlation between runoff and temperature, whereas a watershed located in a arid climate, such as the Sahara desert, exhibits a negative correlation; often there is no run off at all during the summer months in such watersheds

Critical loads for lead in France: First results on forest soils

Within the framework of the United Nation Convention on Long Range Transboundary Air Pollution, France is part of the Working Group on Effect which aims at evaluating the impact of atmospheric deposition on ecosystems by calculating critical loads. The critical loads are the highest deposition of compounds that will not cause chemical changes in soil leading to long-term harmful effects on ecosystem structure and function. A guidance manual for calculation of critical loads for heavy metals (lead and cadmium) has been proposed by the Coordination Center for Effects (executive body of the WGE). French National Focal Center (CNRS and ADEME) aims in this study at evaluating the accuracy of the european methodology for calculation of critical loads for french forest soils. It appears that critical load approach is adapted for France but need to be calibrated at least for calculation of weathering fluxes and determination of critical limits. Stand-still on the contrary is not adequate because of inherent contradictions in the method and too much uncertainties in the transfer functions

Evidence for global runoff increase related to climate warming

Ongoing global climatic change initiated by the anthropogenic release of carbon dioxide is a matter of intense debate. We focus both on the impact of these climatic changes on the global hydrological cycle and on the amplitude of the increase of global and continental runoff over the last century, in relation to measured temperature increases. In this contribution, we propose an original statistical wavelet-based method for the reconstruction of the monthly discharges of worldwide largest rivers. This method provides a data-based approximation of the evolution of the annual continental and global runoffs over the last century. A consistent correlation is highlighted between global annual temperature and runoff, suggesting a 4% global runoff increase by 1 C global temperature rise. However, this global trend should be qualified at the regional scale where both increasing and decreasing trends are identified. North America runoffs appear to be the most sensitive to the recent climatic changes. Finally, this contribution provides the first experimental data-based evidence demonstrating the link between the global warming and the intensification of the global hydrological cycle. This corresponds to more intense evaporation over oceans coupled to continental precipitation increase or continental evaporation decrease. This process finally leads to an increase of the global continental runoff

A new approach for modeling Cenozoic oceanic lithium isotope paleo-variations: the key role of climate

International audienceThe marine record of ocean lithium isotope composition may provide important information constraining the factors that control continental weathering and how they have varied in the past. However, the equations establishing links between the continental flux of Li to the ocean, the continental Li isotope composition and the ocean Li isotope composition are under-constrained, and their resolution are related to significant uncertainties. In order to partially reduce this uncertainty, we propose a new approach that couples the C and Li cycles, such that our proposed reconstruction of the Cenozoic Li cycle is compatible with the required stability of the exospheric carbon cycle on geological timescales. The results of this exercise show, contrary to expectations, that the Cenozoic evolution of the Li isotope composition of rivers did not necessarily mimic the oceanic δ7Li rise. In contrast, variations in the continental flux of Li to the ocean are demonstrated to play a major role in setting the ocean δ7Li. We also provide evidence that Li storage in secondary phases is an important element of the global Li cycle that cannot be neglected, in particular during the early Cenozoic. Our modeling of the published foraminifera record highlights a close link between soil formation rate and indexes recording the climate evolution during the Cenozoic, such as foraminifera δ18O and pCO2 reconstructions. This leads us to conclude that the Li isotope record does not provide persuasive, unique evidence for erosional forcing of Cenozoic change because it could alternatively be consistent with a climatic control on soil production rates

Designing a Suite of Models to Explore Critical Zone Function

Critical Zone; weathering; hydrology; ecology; watershedsThe Critical Zone (CZ) incorporates all aspects of the earth's environment from the vegetation canopy to the bottom of groundwater. CZ researchers target processes that cross timescales from that of water fluxes (milliseconds to decades) to that of the evolution of landforms (thousands to tens of millions of years). Conceptual and numerical models are used to investigate the important fluxes: water, energy, solutes, carbon, nitrogen, and sediments. Depending upon the questions addressed, these models must calculate the distribution of landforms, regolith structure and chemistry, biota, and the chemistry of water, solutes, sediments, and soil atmospheres. No single model can accomplish all these objectives. We are designing a group of models or model capabilities to explore the CZ and testing them at the Susquehanna Shale Hills CZ Observatory. To examine processes over different timescales, we establish the core hydrologic fluxes using the Penn State Integrated Hydrologic Model (PIHM) – and then augment PIHM with simulation modules. For example, most land-atmosphere models currently do not incorporate an accurate representation of the geologic subsurface. We are exploring what aspects of subsurface structure must be accurately modelled to simulate water, carbon, energy, and sediment fluxes accurately. Only with a suite of modeling tools will we learn to forecast – earthcast -- the future CZ

Designing a Suite of Models to Explore Critical Zone Function

Critical Zone; weathering; hydrology; ecology; watershedsThe Critical Zone (CZ) incorporates all aspects of the earth's environment from the vegetation canopy to the bottom of groundwater. CZ researchers target processes that cross timescales from that of water fluxes (milliseconds to decades) to that of the evolution of landforms (thousands to tens of millions of years). Conceptual and numerical models are used to investigate the important fluxes: water, energy, solutes, carbon, nitrogen, and sediments. Depending upon the questions addressed, these models must calculate the distribution of landforms, regolith structure and chemistry, biota, and the chemistry of water, solutes, sediments, and soil atmospheres. No single model can accomplish all these objectives. We are designing a group of models or model capabilities to explore the CZ and testing them at the Susquehanna Shale Hills CZ Observatory. To examine processes over different timescales, we establish the core hydrologic fluxes using the Penn State Integrated Hydrologic Model (PIHM) – and then augment PIHM with simulation modules. For example, most land-atmosphere models currently do not incorporate an accurate representation of the geologic subsurface. We are exploring what aspects of subsurface structure must be accurately modelled to simulate water, carbon, energy, and sediment fluxes accurately. Only with a suite of modeling tools will we learn to forecast – earthcast -- the future CZ

Assessing Volcanic Controls on Miocene Climate Change

The Miocene period saw substantially warmer Earth surface temperatures than today, particularly during a period of global warming called the Mid Miocene Climatic Optimum (MMCO; ∼17–15 Ma). However, the long-term drivers of Miocene climate remain poorly understood. By using a new continuous climate-biogeochemical model (SCION), we can investigate the interaction between volcanism, climate and biogeochemical cycles through the Miocene. We identify high tectonic CO2 degassing rates and further emissions associated with the emplacement of the Columbia River Basalt Group as the primary driver of the background warmth and the MMCO respectively. We also find that enhanced weathering of the basaltic terrane and input of explosive volcanic ash to the oceans are not sufficient to drive the immediate cooling following the MMCO and suggest that another mechanism, perhaps the change in ocean chemistry due to massive evaporite deposition, was responsible

The role of sea-level change and marine anoxia in the Frasnian-Famennian (Late Devonian) mass extinction

Johnson et al. (Johnson, J.G., Klapper, G., Sandberg, C.A., 1985. Devonian eustatic fluctuations in Euramerica. Geological Society of America Bulletin 96, 567–587) proposed one of the first explicit links between marine anoxia, transgression and mass extinction for the Frasnian–Famennian (F–F, Late Devonian) mass extinction. This cause-and-effect nexus has been accepted by many but others prefer sea-level fall and cooling as an extinction mechanism. New facies analysis of sections in the USA and Europe (France, Germany, Poland), and comparison with sections known from the literature in Canada, Australia and China reveal several high-frequency relative sea-level changes in the late Frasnian to earliest Famennian extinction interval. A clear signal of major transgression is seen within the Early rhenana Zone (e.g. drowning of the carbonate platform in the western United States). This is the base of transgressive–regressive Cycle IId of the Johnson et al. (Johnson, J.G., Klapper, G., Sandberg, C.A., 1985. Devonian eustatic fluctuations in Euramerica. Geological Society of America Bulletin 96, 567–587) eustatic curve. This was curtailed by regression and sequence boundary generation within the early linguiformis Zone, recorded by hardground and karstification surfaces in sections from Canada to Australia. This major eustatic fall probably terminated platform carbonate deposition over wide areas, especially in western North America. The subsequent transgression in the later linguiformis Zone, recorded by the widespread development of organic-rich shale facies, is also significant because it is associated with the expansion of anoxic deposition, known as the Upper Kellwasser Event. Johnson et al.'s (Johnson, J.G., Klapper, G., Sandberg, C.A., 1985. Devonian eustatic fluctuations in Euramerica. Geological Society of America Bulletin 96, 567–587) original transgression-anoxia–extinction link is thus supported, although some extinction losses of platform carbonate biota during the preceeding regression cannot be ruled out. Conodont faunas suffered major losses during the Upper Kellwasser Event, with deep-water taxa notably affected. This renders unreliable any eustatic analyses utilising changes in conodont biofacies. Claims for a latest Frasnian regression are not supported, and probably reflect poor biostratigraphic dating of the early linguiformis Zone sequence boundary

Atmospheric turbulence triggers pronounced diel pattern in karst carbonate geochemistry

CO2 exchange between terrestrial ecosystems and the atmosphere is key to understanding the feedbacks between climate change and the land surface. In regions with carbonaceous parent material, CO2 exchange patterns occur that cannot be explained by biological processes, such as disproportionate outgassing during the daytime or nighttime CO2 uptake during periods when all vegetation is senescent. Neither of these phenomena can be attributed to carbonate weathering reactions, since their CO2 exchange rates are too small. Soil ventilation induced by high atmospheric turbulence is found to explain atypical CO2 exchange between carbonaceous systems and the atmosphere. However, by strongly altering subsurface CO2 concentrations, ventilation can be expected to influence carbonate weathering rates. By imposing ventilation-driven CO2 outgassing in a carbonate weathering model, we show here that carbonate geochemistry is accelerated and does play a surprisingly large role in the observed CO2 exchange pattern of a semi-arid ecosystem. We found that by rapidly depleting soil CO2 during the daytime, ventilation disturbs soil carbonate equilibria and therefore strongly magnifies daytime carbonate precipitation and associated CO2 production. At night, ventilation ceases and the depleted CO2 concentrations increase steadily. Dissolution of carbonate is now enhanced, which consumes CO2 and largely compensates for the enhanced daytime carbonate precipitation. This is why only a relatively small effect on global carbonate weathering rates is to be expected. On the short term, however, ventilation has a drastic effect on synoptic carbonate weathering rates, resulting in a pronounced diel pattern that exacerbates the non-biological behavior of soil–atmosphere CO2 exchanges in dry regions \mbox{with carbonate soils}.M. Roland was granted by the Institute for Promotion of Innovation through Science and Technology in Flanders (IWT-Vlaanderen). I. A. Janssens and R. Van Grieken acknowledge the Research Foundation – Flanders (FWO). P. Serrano-Ortiz is funded by a postdoctoral fellowship from the Spanish Ministry of Science and Innovation. S. Cuezva was funded by a postdoctoral fellowship from the Spanish Ministry of Science and Innovation, research programme Juan de la Cierva