RegCM3 nested in HadAM3 scenarios A2 and B2: projected changes in extratropical cyclogenesis, temperature and precipitation over the South Atlantic Ocean

The RegCM3 (Regional Climate Model-version 3) was nested in HadAM3 model to simulate present (1975-1989, referred hereafter as RegHad) and two future climate scenarios (A2 and B2 from 2071 to 2085, referred as RegA2 and RegB2) over the South America (SA) and South Atlantic Ocean (SAO). Projected changes in the air temperature, precipitation, low level circulation and cyclogenesis climatology were investigated. The cyclogenesis were identified using an automatic scheme for tracking based on the minimum of relative vorticity (zeta) from 10-m height wind. During summer, a general decrease (increase) in the precipitation is projected by RegA2 and RegB2 over the northeastern SA (center-west and south Brazil, north Argentina and Uruguay). For winter, an anomalous low level anticyclonic circulation is associated with the reduction in the rainfall over the central part of southern Brazil in RegA2 and RegB2 scenarios. Similar to HadAM3, RegCM3 projects larger warming in A2 scenario. For the present climate, when compared to HadAM3, RegHad defines better both the location of the main cyclogenetic areas and its annual cycle near southwestern SAO. The projections indicate a reduction in the total number of cyclones of -7.2% and -4.7% for RegA2 and RegB2, respectively, while HadAM3 reduction is -4.5% for both scenarios. The decrease is larger for initially intense cyclones (zeta <=-<= 2.5 x 10(-5) s-(1)): -20.9% (RegA2) and -11.3% (RegB2). For the lifetime, distance traveled and mean velocity of the cyclones, the A2 and B2 scenarios present mean values close to the present climate ( 3 days, 1900 km, and 9 m s(-1), respectively). Regarding the initial mean vorticity of the systems, RegB2 simulates values similar to the present climate, but they are initially weaker in RegA2. In general, RegA2 and RegB2 show a large decrease in the number of cyclones over the southern SAO due to an anticyclonic anomaly covering SAO between 30-55A degrees S. The reduction is larger in the scenario with higher concentrations of greenhouse gases (RegA2).

Estimation of vertical groundwater fluxes into a streambed through continuous temperature profile monitoring and the relationship of groundwater fluxes to coaster brook trout spawning habitat

We hypothesized that the spatial distribution of groundwater inflows through river bottom sediments is a critical factor associated with the selection of coaster brook trout (a life history variant of Salvelinus fontinalis,) spawning sites. An 80-m reach of the Salmon Trout River, in the Huron Mountains of the upper peninsula of Michigan, was selected to test the hypothesis based on long-term documentation of coaster brook trout spawning at this site. Throughout this site, the river is relatively similar along its length with regard to stream channel and substrate features. A monitoring well system consisting of an array of 27 wells was installed to measure subsurface temperatures underneath the riverbed over a 13-month period. The monitoring well locations were separated into areas where spawning has and has not been observed. Over 200,000 total temperature measurements were collected from 5 depths within each of the 27 monitoring wells. Temperatures within the substrate at the spawning area were generally cooler and less variable than river temperatures. Substrate temperatures in the non-spawning area were generally warmer, more variable, and closely tracked temporal variations in river temperatures. Temperature data were inverted to obtain subsurface groundwater velocities using a numerical approximation of the heat transfer equation. Approximately 45,000 estimates of groundwater velocities were obtained. Estimated velocities in the spawning and non-spawning areas confirmed that groundwater velocities in the spawning area were primarily in the upward direction, and were generally greater in magnitude than velocities in the non-spawning area. In the non-spawning area there was a greater occurrence of velocities in the downward direction, and velocity estimates were generally lesser in magnitude than in the spawning area. Both the temperature and velocity results confirm the hypothesis that spawning sites correspond to areas of significant groundwater influx to the river bed.

Metamorphic and geochronogical study of the Triassic El Oro metamorphic complex, Ecuador: Implications for high-temperature metamorphism in a forearc zone

In the forearc of the Andean active margin in southwest Ecuador, the El Oro metamorphic complex exhibits a well exposed tilted forearc section partially migmatized. We used Raman spectroscopy on carbonaceous matter (RSCM) thermometry and pseudosections coupled with mineralogical and textural studies to constrain the pressure–temperature (P–T) evolution of the El Oro metamorphic complex during Triassic times. Our results show that anatexis of the continental crust occurred by white-mica and biotite dehydration melting along a 10 km thick crustal domain (from 4.5 to 8 kbar) with increasing temperature from 650 to 700 °C. In the biotite dehydration melting zone, temperature was buffered at 750–820 °C in a 5 km thick layer. The estimated average thermal gradient during peak metamorphism is of 30 °C/km within the migmatitic domain can be partitioned into two apparent gradients parts. The upper part from surface to 7 km depth records a 40–45 °C/km gradient. The lower part records a quasi-adiabatic geotherm with a 10 °C/km gradient consistent with an isothermal melting zone. Migmatites U–Th–Pb geochronology yielded zircon and monazite ages of 229.3 ± 2.1 Ma and 224.5 ± 2.3 Ma, respectively. This thermal event generated S-type magmatism (the Marcabeli granitoid) and was immediately followed by underplating of the high-pressure low-temperature (HP-LT) Arenillas–Panupalí unit at 225.8 ± 1.8 Ma. The association of high-temperature low-pressure (HT-LP) migmatites with HP-LT unit constitutes a new example of a paired metamorphic belt along the South American margin. We propose that in addition to crustal thinning, underplating of the Piedras gabbroic unit before 230 Ma provided the heat source necessary to foster crustal anatexis. Furthermore, its MORB signature shows that the asthenosphere was involved as the source of the heat anomaly. S-type felsic magmatism is widespread during this time and suggests that a large-scale thermal anomaly affected a large part of the South American margin during the late Triassic. We propose that crustal anatexis is related to an anomaly that arose during subduction of the Panthalassa ocean under the South American margin. Slab verticalization or slab break-off can be invoked as the origin of the upwelling of the asthenosphere.

Large-scale temperature response to external forcing in simulations and reconstructions of the last millennium

Understanding natural climate variability and its driving factors is crucial to assessing future climate change. Therefore, comparing proxy-based climate reconstructions with forcing factors as well as comparing these with paleoclimate model simulations is key to gaining insights into the relative roles of internal versus forced variability. A review of the state of modelling of the climate of the last millennium prior to the CMIP5–PMIP3 (Coupled Model Intercomparison Project Phase 5–Paleoclimate Modelling Intercomparison Project Phase 3) coordinated effort is presented and compared to the available temperature reconstructions. Simulations and reconstructions broadly agree on reproducing the major temperature changes and suggest an overall linear response to external forcing on multidecadal or longer timescales. Internal variability is found to have an important influence at hemispheric and global scales. The spatial distribution of simulated temperature changes during the transition from the Medieval Climate Anomaly to the Little Ice Age disagrees with that found in the reconstructions. Thus, either internal variability is a possible major player in shaping temperature changes through the millennium or the model simulations have problems realistically representing the response pattern to external forcing. A last millennium transient climate response (LMTCR) is defined to provide a quantitative framework for analysing the consistency between simulated and reconstructed climate. Beyond an overall agreement between simulated and reconstructed LMTCR ranges, this analysis is able to single out specific discrepancies between some reconstructions and the ensemble of simulations. The disagreement is found in the cases where the reconstructions show reduced covariability with external forcings or when they present high rates of temperature change.

Dating groundwater in the Bohemian Cretaceous Basin: Understanding tracer variations in the subsurface

The northern section of the Bohemian Cretaceous Basin has been the site of intensive U exploitation with harmful impacts on groundwater quality. The understanding of groundwater flow and age distribution is crucial for the prediction of the future dispersion and impact of the contamination. State of the art tracer methods (3H, 3He, 4He, 85Kr, 39Ar and 14C) were, therefore, used to obtain insights to ageing and mixing processes of groundwater along a north–south flow line in the centre of the two most important aquifers of Cenomanian and middle Turonian age. Dating of groundwater is particularly complex in this area as: (i) groundwater in the Cenomanian aquifer is locally affected by fluxes of geogenic and biogenic gases (e.g. CO2, CH4, He) and by fossil brines in basement rocks rich in Cl and SO4; (ii) a thick unsaturated zone overlays the Turonian aquifer; (iii) a periglacial climate and permafrost conditions prevailed during the Last Glacial Maximum (LGM), and iv) the wells are mostly screened over large depth intervals. Large disagreements in 85Kr and 3H/3He ages indicate that processes other than ageing have affected the tracer data in the Turonian aquifer. Mixing with older waters (>50 a) was confirmed by 39Ar activities. An inverse modelling approach, which included time lags for tracer transport throughout the unsaturated zone and degassing of 3He, was used to estimate the age of groundwater. Best fits between model and field results were obtained for mean residence times varying from modern up to a few hundred years. The presence of modern water in this aquifer is correlated with the occurrence of elevated pollution (e.g. nitrates). An increase of reactive geochemical indicators (e.g. Na) and radiogenic 4He, and a decrease in 14C along the flow direction confirmed groundwater ageing in the deeper confined Cenomanian aquifer. Radiocarbon ages varied from a few hundred years to more than 20 ka. Initial 14C activity for radiocarbon dating was calibrated by means of 39Ar measurements. The 14C age of a sample recharged during the LGM was further confirmed by depleted stable isotope signatures and near freezing point noble gas temperature. Radiogenic 4He accumulated in groundwater with concentrations increasing linearly with 14C ages. This enabled the use of 4He to validate the dating range of 14C and extend it to other parts of this aquifer. In the proximity of faults, 39Ar in excess of modern concentrations and 14C dead CO2 sources, elevated 3He/4He ratios and volcanic activity in Oligocene to Quaternary demonstrate the influence of gas of deeper origin and impeded the application of 4He, 39Ar and 14C for groundwater dating.

Seasonal Climatology of Hydrographic Conditions in the Upwelling Region Off Northern Chile

Over 30 years of hydrographic data from the northern Chile (18 degreesS-24 degreesS) upwelling region are used to calculate the surface and subsurface seasonal climatology extending 400 km offshore. The data are interpolated to a grid with sufficient spatial resolution to preserve cross-shelf gradients and then presented as means within four seasons: austral winter (July-September), spring (October-December), summer (January-March), and fall (April-June). Climatological monthly wind forcing, surface temperature, and sea level from three coastal stations indicate equatorward (upwelling favorable) winds throughout the year, weakest in the north. Seasonal maximum alongshore wind stress is in late spring and summer (December-March). Major water masses of the region are identified in climatological T-S plots and their sources and implied circulation discussed. Surface fields and vertical transects of temperature and salinity confirm that upwelling occurs year-round, strongest in summer and weakest in winter, bringing relatively fresh water to the surface nearshore. Surface geostrophic flow nearshore is equatorward throughout the year. During summer, an anticyclonic circulation feature in the north which extends to at least 200 m depth is evident in geopotential anomaly and in both temperature and geopotential variance fields. Subsurface fields indicate generally poleward flow throughout the year, strongest in an undercurrent near the coast. This undercurrent is strongest in summer and most persistent and organized in the south (south of 21 degreesS), A subsurface oxygen minimum, centered at similar to 250 m, is strongest at lower latitudes. Low-salinity subsurface water intrudes into the study area near 100 m, predominantly in offshore regions, strongest during summer and fall and in the southernmost portion of the region. The climatological fields are compared to features off Baja within the somewhat analogous California Current and to measurements from higher latitudes within the Chile-Peru Current system.

Reconstructing the subsurface ocean decadal variability using surface nudging in a perfect model framework

Initialising the ocean internal variability for decadal predictability studies is a new area of research and a variety of ad hoc methods are currently proposed. In this study, we explore how nudging with sea surface temperature (SST) and salinity (SSS) can reconstruct the threedimensional variability of the ocean in a perfect model framework. This approach builds on the hypothesis that oceanic processes themselves will transport the surface information into the ocean interior as seen in ocean-only simulations. Five nudged simulations are designed to reconstruct a 150 years ‘‘target’’ simulation, defined as a portion of a long control simulation. The nudged simulations differ by the variables restored to, SST or SST + SSS, and by the area where the nudging is applied. The strength of the heat flux feedback is diagnosed from observations and the restoring coefficients for SSS use the same time-scale. We observed that this choice prevents spurious convection at high latitudes and near sea-ice border when nudging both SST and SSS. In the tropics, nudging the SST is enough to reconstruct the tropical atmosphere circulation and the associated dynamical and thermodynamical impacts on the underlying ocean. In the tropical Pacific Ocean, the profiles for temperature show a significant correlation from the surface down to 2,000 m, due to dynamical adjustment of the isopycnals. At mid-tohigh latitudes, SSS nudging is required to reconstruct both the temperature and the salinity below the seasonal thermocline. This is particularly true in the North Atlantic where adding SSS nudging enables to reconstruct the deep convection regions of the target. By initiating a previously documented 20-year cycle of the model, the SST + SSS nudging is also able to reproduce most of the AMOC variations, a key source of decadal predictability. Reconstruction at depth does not significantly improve with amount of time spent nudging and the efficiency of the surface nudging rather depends on the period/events considered. The joint SST + SSS nudging applied verywhere is the most efficient approach. It ensures that the right water masses are formed at the right surface density, the subsequent circulation, subduction and deep convection further transporting them at depth. The results of this study underline the potential key role of SSS for decadal predictability and further make the case for sustained largescale observations of this field.

Sea surface temperature and sea ice variability in the subpolar North Atlantic from explosive volcanism of the late thirteenth century

In this study, we use IP and alkenone biomarker proxies to document the subdecadal variations of sea ice and sea surface temperature in the subpolar North Atlantic induced by the decadally paced explosive tropical volcanic eruptions of the second half of the thirteenth century. The short- and long-term evolutions of both variables were investigated by cross analysis with a simulation of the IPSL-CM5A LR model. Our results show short-term ocean cooling and sea ice expansion in response to each volcanic eruption. They also highlight that the long response time of the ocean leads to cumulative surface cooling and subsurface heat buildup due to sea ice capping. As volcanic forcing relaxes, the surface ocean rapidly warms, likely amplified by subsurface heat, and remains almost ice free for several decades

Effect of surface restoring on subsurface variability in a climate model during 1949–2005

Initializing the ocean for decadal predictability studies is a challenge, as it requires reconstructing the little observed subsurface trajectory of ocean variability. In this study we explore to what extent surface nudging using well-observed sea surface temperature (SST) can reconstruct the deeper ocean variations for the 1949–2005 period. An ensemble made with a nudged version of the IPSLCM5A model and compared to ocean reanalyses and reconstructed datasets. The SST is restored to observations using a physically-based relaxation coefficient, in contrast to earlier studies, which use a much larger value. The assessment is restricted to the regions where the ocean reanalyses agree, i.e. in the upper 500 m of the ocean, although this can be latitude and basin dependent. Significant reconstruction of the subsurface is achieved in specific regions, namely region of subduction in the subtropical Atlantic, below the thermocline in the equatorial Pacific and, in some cases, in the North Atlantic deep convection regions. Beyond the mean correlations, ocean integrals are used to explore the time evolution of the correlation over 20-year windows. Classical fixed depth heat content diagnostics do not exhibit any significant reconstruction between the different existing bservation-based references and can therefore not be used to assess global average time-varying correlations in the nudged simulations. Using the physically based average temperature above an isotherm (14°C) alleviates this issue in the tropics and subtropics and shows significant reconstruction of these quantities in the nudged simulations for several decades. This skill is attributed to the wind stress reconstruction in the tropics, as already demonstrated in a perfect model study using the same model. Thus, we also show here the robustness of this result in an historical and observational context.