Results of modelled Arctic river discharge trends since 7 ka BP

The Arctic hydrological cycle throughout the Holocene is analyzed based on the results of transient simulations with the coupled atmosphere-ocean circulation model ECHO-G. The results suggest a ~ 2 % increase of mid-Holocene to preindustrial Arctic river discharges for the Eurasian continent. However, rivers of the North America Arctic realm show a moderate runoff decline of approximately 4 to 5 % for the same period. The total river discharge into the Arctic Ocean has remained at an approximately constant preindustrial level since the mid Holocene. The positive discharge trend within Eurasia is caused by a more rapid decrease in local net evaporation compared to a smaller decline in advected moisture and hence precipitation. This effect is neither recognized within the North American Arctic domain nor in the far eastern part of the Eurasian Arctic realm. A detailed comparison of these model findings with a variety of proxy studies is conducted. The collected proxy records show trends of continental surface temperatures and precipitation rates that are consistent with the simulations. A continuation of the transient Holocene runs for the 19th and 20th century with increased greenhouse gases indicates an increase of the total river influx into the Arctic Ocean of up to 7.6 %. The Eurasian river discharges increase by 7.5 %, the North American discharges by up to 8.4 %. The most rapid increases have been detected since the beginning of the 20th century. These results are corroborated by the observed rising of Arctic river discharges during the last century which is attributed to anthropogenic warming. The acceleration of the Arctic hydrological cycle in the 20th century is without precedence in the Holocene.

Eocene sedimentary calcium carbonate contents and stable isotope composition of benthic foraminifera

'Hyperthermals' are intervals of rapid, pronounced global warming known from six episodes within the Palaeocene and Eocene epochs (~65-34 million years (Myr) ago) (Zachos et al., 2005, doi:10.1126/science.1109004; 2008, doi:10.1038/nature06588; Roehl et al., 2007, doi:10.1029/2007GC001784; Thomas et al., 2000; Cramer et al., 2003, doi:10.1029/2003PA000909; Lourens et al., 2005, doi:10.1038/nature03814; Petrizzo, 2005, doi:10.2973/odp.proc.sr.198.102.2005; Sexton et al., 2006, doi:10.1029/2005PA001253; Westerhold et al., 2007, doi:10.1029/2006PA001322; Edgar et al., 2007, doi:10.1038/nature06053; Nicolo et al., 2007, doi:10.1130/G23648A.1; Quillévéré et al., 2008, doi:10.1016/j.epsl.2007.10.040; Stap et al., 2010, doi:10.1130/G30777.1). The most extreme hyperthermal was the 170 thousand year (kyr) interval (Roehl et al., 2007) of 5-7 °C global warming (Zachos et al., 2008) during the Palaeocene-Eocene Thermal Maximum (PETM, 56 Myr ago). The PETM is widely attributed to massive release of greenhouse gases from buried sedimentary carbon reservoirs (Zachos et al., 2005; 2008; Lourenbs et al., 2005; Nicolo et al., 2007; Dickens et al., 1995, doi:10.1029/95PA02087; Dickens, 2000; 2003, doi:10.1016/S0012-821X(03)00325-X; Panchuk et al., 2008, doi:10.1130/G24474A.1) and other, comparatively modest, hyperthermals have also been linked to the release of sedimentary carbon (Zachos et al., 2008, Lourens et al., 2005; Nicolo et al., 2007; Dickens, 2003; Panchuk et al., 2003). Here we show, using new 2.4-Myr-long Eocene deep ocean records, that the comparatively modest hyperthermals are much more numerous than previously documented, paced by the eccentricity of Earth's orbit and have shorter durations (~40 kyr) and more rapid recovery phases than the PETM. These findings point to the operation of fundamentally different forcing and feedback mechanisms than for the PETM, involving redistribution of carbon among Earth's readily exchangeable surface reservoirs rather than carbon exhumation from, and subsequent burial back into, the sedimentary reservoir. Specifically, we interpret our records to indicate repeated, large-scale releases of dissolved organic carbon (at least 1,600 gigatonnes) from the ocean by ventilation (strengthened oxidation) of the ocean interior. The rapid recovery of the carbon cycle following each Eocene hyperthermal strongly suggests that carbon was resequestered by the ocean, rather than the much slower process of silicate rock weathering proposed for the PETM (Zachos et al., 2005; 2003). Our findings suggest that these pronounced climate warming events were driven not by repeated releases of carbon from buried sedimentary sources (Zachos et al., 2008, Lourens et al., 2005; Nicolo et al., 2007; Dickens, 2003; Panchuk et al., 2003) but, rather, by patterns of surficial carbon redistribution familiar from younger intervals of Earth history.

Modeling the sensitivity of precipitation oxygen isotopes to the Last Glacial Maximum boundary conditions: A study using Community Atmosphere Model (CAM3.0)

To understand the validity of d18O proxy records as indicators of past temperature change, a series of experiments was conducted using an atmospheric general circulation model fitted with water isotope tracers (Community Atmosphere Model version 3.0, IsoCAM). A pre-industrial simulation was performed as the control experiment, as well as a simulation with all the boundary conditions set to Last Glacial Maximum (LGM) values. Results from the pre-industrial and LGM simulations were compared to experiments in which the influence of individual boundary conditions (greenhouse gases, ice sheet albedo and topography, sea surface temperature (SST), and orbital parameters) were changed each at a time to assess their individual impact. The experiments were designed in order to analyze the spatial variations of the oxygen isotopic composition of precipitation (d18Oprecip) in response to individual climate factors. The change in topography (due to the change in land ice cover) played a significant role in reducing the surface temperature and d18Oprecip over North America. Exposed shelf areas and the ice sheet albedo reduced the Northern Hemisphere surface temperature and d18Oprecip further. A global mean cooling of 4.1 °C was simulated with combined LGM boundary conditions compared to the control simulation, which was in agreement with previous experiments using the fully coupled Community Climate System Model (CCSM3). Large reductions in d18Oprecip over the LGM ice sheets were strongly linked to the temperature decrease over them. The SST and ice sheet topography changes were responsible for most of the changes in the climate and hence the d18Oprecip distribution among the simulations.

Data from the EPICA Dome C ice core EDC

The Antarctic Vostok ice core provided compelling evidence of the nature of climate, and of climate feedbacks, over the past 420,000 years. Marine records suggest that the amplitude of climate variability was smaller before that time, but such records are often poorly resolved. Moreover, it is not possible to infer the abundance of greenhouse gases in the atmosphere from marine records. Here we report the recovery of a deep ice core from Dome C, Antarctica, that provides a climate record for the past 740,000 years. For the four most recent glacial cycles, the data agree well with the record from Vostok. The earlier period, between 740,000 and 430,000 years ago, was characterized by less pronounced warmth in interglacial periods in Antarctica, but a higher proportion of each cycle was spent in the warm mode. The transition from glacial to interglacial conditions about 430,000 years ago (Termination V) resembles the transition into the present interglacial period in terms of the magnitude of change in temperatures and greenhouse gases, but there are significant differences in the patterns of change. The interglacial stage following Termination V was exceptionally long - 28,000 years compared to, for example, the 12,000 years recorded so far in the present interglacial period. Given the similarities between this earlier warm period and today, our results may imply that without human intervention, a climate similar to the present one would extend well into the future.

(Supplemental Table 2) Plant characteristics, GHG fluxes, and soil chemical and microbial properties in experimental treatments in Alexandra Fiord

We evaluated above- and belowground ecosystem changes in a 16 year, combined fertilization and warming experiment in a High Arctic tundra deciduous shrub heath (Alexandra Fiord, Ellesmere Island, NU, Canada). Soil emissions of the three key greenhouse gases (GHGs) (carbon dioxide, methane, and nitrous oxide) were measured in mid-July 2009 using soil respiration chambers attached to a FTIR system. Soil chemical and biochemical properties including Q10 values for CO2, CH4, and N2O, Bacteria and Archaea assemblage composition, and the diversity and prevalence of key nitrogen cycling genes including bacterial amoA, crenarchaeal amoA, and nosZ were measured. Warming and fertilization caused strong increases in plant community cover and height but had limited effects on GHG fluxes and no substantial effect on soil chemistry or biochemistry. Similarly, there was a surprising lack of directional shifts in the soil microbial community as a whole or any change at all in microbial functional groups associated with CH4 consumption or N2O cycling in any treatment. Thus, it appears that while warming and increased nutrient availability have strongly affected the plant community over the last 16 years, the belowground ecosystem has not yet responded. This resistance of the soil ecosystem has resulted in limited changes in GHG fluxes in response to the experimental treatments.

Stable isotope ratios of foraminifera from various DSDP/ODP sites

The equator to high southern latitude sea surface and vertical temperature gradients are reconstructed from oxygen isotope values of planktonic and benthic foraminifers for the following five time intervals: late Paleocene, early Eocene, early middle Eocene, late Eocene, and early Oligocene. Paleotemperatures are calculated using standard oxygen isotope/temperature equations with adjustments to account for (1) variations in sea water delta18O related to changes in global ice volume over time and (2) latitudinal gradients in surface water delta18O. These reconstructions indicate that sea-surface temperatures (SST) of the Southern Oceans in the early Eocene were as high as 15°C, whereas temperatures during the late Paleocene and early middle Eocene reached maximum levels of 10°-12°C. By the late Eocene and early Oligocene high latitude SST had declined to 6 and 4°C, respectively. For most of the early Paleogene, low latitude sub-tropical temperatures remained constant and well within the range of Holocene temperatures (24°-25°C) but by the late Eocene and early Oligocene declined to values in the range of 18° to 22°C. The late Paleogene apparent decline in tropical temperatures, however, might be artificial because of dissolution of near-surface foraminifera tests which biased sediment assemblages toward deeper-dwelling foraminifera. Moreover, according to recent plate reconstructions, it appears that the majority of sites upon which the late Eocene and early Oligocene tropical temperatures were previously established were located either in or near regions likely to have been influenced by upwelling. Global deepwater temperature on average paralleled southern ocean SST for most of the Paleogene. We speculate based on the overall timing and character of marine sea surface temperature variation during the Paleogene that some combination of both higher levels of greenhouse gases and increased heat transport was responsible for the exceptional high-latitude warmth of the early Eocene.

B-isotope and Mg/Ca ratios of planktonic foraminifera from ODP Hole 108-668B and esimates of salinity, alkalinity and pCO2 (Table 1)

Knowledge of past atmospheric pCO2 is important for evaluating the role of greenhouse gases in climate forcing. Ice core records show the tight correlation between climate change and pCO2, but records are limited to the past ~900 kyr. We present surface ocean pH and pCO2 data, reconstructed from boron isotopes in planktonic foraminifera over two full glacial cycles (0-140 and 300-420 kyr). The data co-vary strongly with the Vostok pCO2-record and demonstrate that the coupling between surface ocean chemistry and the atmosphere is recorded in marine archives, allowing for quantitative estimation of atmospheric pCO2 beyond the reach of ice cores.

Ba/Sr ratio in marine barite from five DSDP and ODP holes

The Paleocene-Eocene Thermal Maximum (PETM), ca. 55 Ma, was a period of extreme global warming caused by rapid emission of greenhouse gases. It is unknown what ended this episode of greenhouse warming, but high oceanic export productivity over thousands of years (as indicated by high accumulation rates of barium, Ba) may have been a factor in ending this warm period by carbon sequestration. However, Ba has a short oceanic residence time (~10 k.y.), so a prolonged global increase in Ba accumulation rates requires an increase in input of Ba to the ocean, increasing barite saturation. We use a novel proxy for barite saturation (Sr/Ba in marine barite) to demonstrate that the seawater saturation state with respect to barite did not change across the PETM. The observations of increased barite burial, no change in saturation, and the short residence time can be reconciled if Ba burial decreased at continental margin and shelf sites due to widespread occurrence of suboxic conditions, leading to Ba release into the water column, combined with increased biological export production at some pelagic sites, resulting in Ba sink reorganization.

Stable carbon and oxygen isotope compositions of Eocene benthic foraminifera from ODP Site 207-1258

'Hyperthermals' are intervals of rapid, pronounced global warming known from six episodes within the Palaeocene and Eocene epochs (~65-34 million years (Myr) ago) (Zachos et al., 2005, doi:10.1126/science.1109004; 2008, doi:10.1038/nature06588; Roehl et al., 2007, doi:10.1029/2007GC001784; Thomas et al., 2000; Cramer et al., 2003, doi:10.1029/2003PA000909; Lourens et al., 2005, doi:10.1038/nature03814; Petrizzo, 2005, doi:10.2973/odp.proc.sr.198.102.2005; Sexton et al., 2006, doi:10.1029/2005PA001253; Westerhold et al., 2007, doi:10.1029/2006PA001322; Edgar et al., 2007, doi:10.1038/nature06053; Nicolo et al., 2007, doi:10.1130/G23648A.1; Quillévéré et al., 2008, doi:10.1016/j.epsl.2007.10.040; Stap et al., 2010, doi:10.1130/G30777.1). The most extreme hyperthermal was the 170 thousand year (kyr) interval (Roehl et al., 2007) of 5-7 °C global warming (Zachos et al., 2008) during the Palaeocene-Eocene Thermal Maximum (PETM, 56 Myr ago). The PETM is widely attributed to massive release of greenhouse gases from buried sedimentary carbon reservoirs (Zachos et al., 2005; 2008; Lourenbs et al., 2005; Nicolo et al., 2007; Dickens et al., 1995, doi:10.1029/95PA02087; Dickens, 2000; 2003, doi:10.1016/S0012-821X(03)00325-X; Panchuk et al., 2008, doi:10.1130/G24474A.1) and other, comparatively modest, hyperthermals have also been linked to the release of sedimentary carbon (Zachos et al., 2008, Lourens et al., 2005; Nicolo et al., 2007; Dickens, 2003; Panchuk et al., 2003). Here we show, using new 2.4-Myr-long Eocene deep ocean records, that the comparatively modest hyperthermals are much more numerous than previously documented, paced by the eccentricity of Earth's orbit and have shorter durations (~40 kyr) and more rapid recovery phases than the PETM. These findings point to the operation of fundamentally different forcing and feedback mechanisms than for the PETM, involving redistribution of carbon among Earth's readily exchangeable surface reservoirs rather than carbon exhumation from, and subsequent burial back into, the sedimentary reservoir. Specifically, we interpret our records to indicate repeated, large-scale releases of dissolved organic carbon (at least 1,600 gigatonnes) from the ocean by ventilation (strengthened oxidation) of the ocean interior. The rapid recovery of the carbon cycle following each Eocene hyperthermal strongly suggests that carbon was resequestered by the ocean, rather than the much slower process of silicate rock weathering proposed for the PETM (Zachos et al., 2005; 2003). Our findings suggest that these pronounced climate warming events were driven not by repeated releases of carbon from buried sedimentary sources (Zachos et al., 2008, Lourens et al., 2005; Nicolo et al., 2007; Dickens, 2003; Panchuk et al., 2003) but, rather, by patterns of surficial carbon redistribution familiar from younger intervals of Earth history.

Sea surface temperature reconstruction based on alkenones of sediment core M35003-4

Evidence for abrupt climate changes on millennial and shorter timescales is widespread in marine and terrestrial climate records (Dansgard et al., 1993, doi:10.1038/364218a0; Bond et al., 1993, doi:10.1038/365143a0; Charles et al., 1996, doi:10.1016/0012-821X(96)00083-0, Bard et al., 1997, doi:10.1038/385707a0). Rapid reorganization of ocean circulation is considered to exert some control over these changes (Broecker et al., 1985, doi:10.1038/315021a0), as are shifts in the concentrations of atmospheric greenhouse gases (Broecker, 1994, doi:10.1038/372421a0). The response of the climate system to these two influences is fundamentally different: slowing of thermohaline overturn in the North Atlantic Ocean is expected to decrease northward heat transport by the ocean and to induce warming of the tropical Atlantic (Crowley, 1992, doi:10.1029/92PA01058; Manabe and Stouffer, 1997, doi:10.1029/96PA03932), whereas atmospheric greenhouse forcing should cause roughly synchronous global temperature changes (Manabe et al., 1991, doi:10.1175/1520-0442(1991)004<0785:TROACO>2.0.CO;2). So these two mechanisms of climate change should be distinguishable by the timing of surface-water temperature variations relative to changes in deep-water circulation. Here we present a high-temporal-resolution record of sea surface temperatures from the western tropical North Atlantic Ocean which spans the past 29,000 years, derived from measurements of temperature-sensitive alkenone unsaturation in sedimentary organic matter. We find significant warming is documented for Heinrich event H1 (16,900-15,400 calendar years bp) and the Younger Dryas event (12,900-11,600 cal. yr bp), which were periods of intense cooling in the northern North Atlantic. Temperature changes in the tropical and high-latitude North Atlantic are out of phase, suggesting that the thermohaline circulation was the important trigger for these rapid climate changes.

In situ pore water oxygen microprofiles from the Polar Front, Southern South Atlantic Ocean

Without doubt, global climate change is directly linked to the anthropogenic release of greenhouse gases such as carbon dioxide (CO2) and methane (UN IPCC-Report 2007). Therefore, research efforts to comprehend the global carbon cycle have increased during the last years. In the context of the observed changes, it is of particular interest to decipher the role of the hydro-, bio- and atmospheres and how the different compartments of the earth system are affected by the increase of atmospheric CO2. Due to its huge carbon inventory, the marine carbon cycle represents the most important component in this respect. Numerous findings suggest that the Southern Ocean plays a key role in terms of oceanic CO2 uptake. However, an exact quantification of such fluxes of material is hard to achieve for large areas, not least on account of the inaccessibility of this remote region. In particular, there exist so far only few accurate data for benthic carbon fluxes. The latter can be derived from high resolution pore water oxygen profiles, as one possible method. However the ex situ flux determinations carried out on sediment cores, tend to suffer from temperature and pressure artefacts. Alternatively, oxygen microprofiles can be measured in situ, i.e. at the seafloor. Until now, no such data have been published for the Southern Ocean. During the Antarctic Expedition ANT-XXI/4, within the framework of this thesis, in situ and ex situ oxygen profiles were measured and used to derive benthic organic carbon fluxes. Having both types of measurements from the same locations, it was possible to establish a depth-related correction function which was applied subsequently to revise published and additional unpublished carbon fluxes to the seafloor. This resulted in a consistent data base of benthic carbon inputs covering many important sub-regions of the Southern Ocean including the Amundsen and Bellingshausen Seas (southern Pacific), Scotia and Weddell Seas (southern South Atlantic) as well as the Crozet Basin (southern Indian Ocean). Including additional locations on the Antarctic Shelf, there are now 134 new and revised measurement locations, covering almost 180° of the Southern Ocean, for which benthic organic carbon fluxes and sedimentary oxygen penetration depth values are available. Further, benthic carbon fluxes were empirically related to dominant diatom distributions in surface sediments as well as to long-term remotely sensed chlorophyll-a estimates. The comparison of these results with benthic carbon fluxes of the entire Atlantic Ocean reveals significantly higher export efficiencies for the Southern Ocean than have previously been assumed, especially for the area of the opal belt.

Peatland depths and radiocarbon dates, recent decades, North America

Peatland ecosystems store about 500-600 Pg of organic carbon, largely accumulated since the last glaciation. Whether they continue to sequester carbon or release it as greenhouse gases, perhaps in large amounts, is important in Earth's temperature dynamics. Given both ages and depths of numerous dated sample peatlands, their rate of carbon sequestration can be estimated throughout the Holocene. Here we use average values for carbon content per unit volume, the geographical extent of peatlands, and ecological models of peatland establishment and growth, to reconstruct the time-trajectory of peatland carbon sequestration in North America and project it into the future. Peatlands there contain ~163 Pg of carbon. Ignoring effects of climate change and other major anthropogenic disturbances, the rate of carbon accumulation is projected to decline slowly over millennia as reduced net carbon accumulation in existing peatlands is largely balanced by new peatland establishment. Peatlands are one of few long-term terrestrial carbon sinks, probably important for global carbon regulation in future generations. This study contributes to a better understanding of these ecosystems that will assist their inclusion in earth-system models, and therefore their management to maintain carbon storage during climate change.

The Lena River Delta - Land cover classification of tundra environments based on Landsat 7 ETM+ data and its application for upscaling of methane emissions

The Lena River Delta, situated in Northern Siberia (72.0 - 73.8° N, 122.0 - 129.5° E), is the largest Arctic delta and covers 29,000 km**2. Since natural deltas are characterised by complex geomorphological patterns and various types of ecosystems, high spatial resolution information on the distribution and extent of the delta environments is necessary for a spatial assessment and accurate quantification of biogeochemical processes as drivers for the emission of greenhouse gases from tundra soils. In this study, the first land cover classification for the entire Lena Delta based on Landsat 7 Enhanced Thematic Mapper (ETM+) images was conducted and used for the quantification of methane emissions from the delta ecosystems on the regional scale. The applied supervised minimum distance classification was very effective with the few ancillary data that were available for training site selection. Nine land cover classes of aquatic and terrestrial ecosystems in the wetland dominated (72%) Lena Delta could be defined by this classification approach. The mean daily methane emission of the entire Lena Delta was calculated with 10.35 mg CH4/m**2/d. Taking our multi-scale approach into account we find that the methane source strength of certain tundra wetland types is lower than calculated previously on coarser scales.

Time series of air chemistry measurements at Georg-von-Neumayer station, Dronning Maud Land, Antarctica in the years 1982 to 1991

The first Air Chemistry Observatory at the German Antarctic station Georg von Neumayer (GvN) was operated for 10 years from 1982 to 1991. The focus of the established observational programme was on characterizing the physical properties and chemical composition of the aerosol, as well as on monitoring the changing trace gas composition of the background atmosphere, especially concerning greenhouse gases. The observatory was designed by the Institut für Umweltphysik, University of Heidelberg (UHEIIUP). The experiments were installed inside the bivouac lodge, mounted on a sledge and put upon a snow hill to prevent snow accumulation during blizzards. All experiments were under daily control and daily performance protocols were documented. A ventilated stainless steel inlet stack (total height about 3-4 m above the snow surface) with a 50% aerodynamic cut-off diameter around 7-10 µm at wind velocities between 4-10 m/s supplied all experiments with ambient air. Contamination free sampling was realized by several means: (i) The Air Chemistry Observatory was situated in a clean air area about 1500 m south of GvN. Due to the fact that northern wind directions are very rare, contamination from the base can be excluded for most of the time. (ii) The power supply (20 kW) is provided by a cable from the main station, thus no fuel-driven generator is operated in the very vicinity. (iii) Contamination-free sampling is controlled by the permanently recorded wind velocity, wind direction and by condensation particle concentration. Contamination was indicated if one of the following criteria were given: Wind direction within a 330°-30° sector, wind velocity <2.2 m/s or >17.5 m/s, or condensation particle concentrations >2500/cm**3 during summer, >800/cm**3 during spring/autumn and >400/cm**3 during winter. If one or a definable combination of these criteria were given, high volume aerosol sampling and part of the trace gas sampling were interrupted. Starting at 1982 through 1991-01-14 surface ozone was measured with an electrochemical concentration cell (ECC). Surface ozone mixing ratio are given in ppbv = parts per 10**9 by volume. The averaging time corresponds to the given time intervals in the data sheet. The accuracy of the values are better than ±1 ppbv and the detection limit is around 1.0 ppbv. Aerosols were sampled on two Whatman 541 cellulose filters in series and analyzed by ion chromatography at the UHEI-IUP. Generally, the sampling period was seven days but could be up to two weeks on occasion. The air flow was around 100 m**3/h and typically 10000-20000 m**3 of ambient air was forced through the filters for one sample. Concentration values are given in nanogram (ng) per 1 m**3 air at standard pressure and temperature (1013 mbar, 273.16 K). Uncertainties of the values were approximately ±10% to ±15% for the main components MSA, chloride, nitrate, sulfate and sodium, and between ±20% and ±30% for the minor species bromide, ammonium, potassium, magnesium and calcium.

Table 1. pH, Temperature, conductivity and hydrochemistry of samples

Carbon dioxide deep geological storage, especially in deep saline aquifers, is one of the preferred technological options to mitigate the effects of greenhouse gases emissions. Thus, in the last decade, studies characterising the behaviour of potential CO2 deep geological storage sites along with thorough safety assessments have been considered essential in order to minimise the risks associated with these sites. The study of natural analogues represents the best source of reliable information about the expected hydrogeochemical processes involved in the CO2 storage in such deep saline aquifers. In this work, a comprehensive study of the hydrogeochemical features and processes taking place at the natural analogue of the Alicún de las Torres thermal system (Betic Cordillera) has been conducted. Thus, the main water/CO2/rock interaction processes occurring at the thermal system have been identified, quantified and modelled, and a principle conclusion is that the hydrogeochemical evolution of the thermal system is controlled by a global dedolomitization process triggered by gypsum dissolution. This geochemical process generates a different geochemical environment to that which would result from the exclusive dissolution of carbonates from the deep aquifer, which is generally considered as the direct result of CO2 injection in a deep carbonate aquifer. Therefore, discounting of the dedolomitization process in any CO2 deep geological storage may lead to erroneous conclusions. This process will also influence the porosity evolution of the CO2 storage formation, which is a very relevant parameter when evaluating a reservoir for CO2 storage. The geothermometric calculation performed in this work leads to estimate that the thermal water reservoir is located between 650 and 800 m depth, which is very close to the minimum required to inject CO2 in a deep geological storage. It is clear that the proper characterisation of the features and hydrogeochemical processes taking place at a natural system analogous to a man-made deep geological storage will provide useful conceptual, semi-quantitative and even quantitative information about the processes and consequences that may occur at the artificial storage system.