Evolution and degradation of flat‐top mesas in the hyper‐arid Negev,Israel revealed from 10Be cosmogenic nuclides

Mesas are ubiquitous landforms in arid and semiarid regions and are often characterized by horizontal stratified erodible rocks capped by more resistant strata. The accepted conceptual model for mesa evolution and degradation considers reduction in the width of the mesa flat‐top plateau due to cliff retreat but ignores possible denudation of the mesa flat‐top and the rates and mechanism of erosion. In this study we examine mesas in the northeastern hyperarid Negev Desert where they appear in various sizes and morphologies and represent different stages of mesa evolution. The variety of mesas within a single climatic zone allows examination of the process of mesa evolution through time. Two of the four sites examined are characterized by a relatively wide (200–230 m) flat‐top and a thick caprock whereas the other two are characterized by a much narrower remnant flat‐top (several meters) and thinner caprock. We use the concentration of the cosmogenic nuclide ¹⁰Be for: (a) determining the chronology of the various geomorphic features associated with the mesa; and (b) understanding geomorphic processes forming the mesa. The ¹⁰Be data, combined with field observations, suggest a correlation between the width of flat‐top mesa and the denudation and cliff retreat rates. Our results demonstrate that: (a) cliff retreat rates decrease with decreasing width of the flat‐top mesa; (b) vertical denudation rates increase with decreasing width of the flat‐top mesa below a critical value (~60 m, for the Negev Desert); (c) the reduction in the width of the flat‐top mesa is driven mainly by cliff retreat accompanied by extremely slow vertical denudation rate which can persist for a very long time (>10⁶ Ma); and (d) when the width of the mesa decreases below a certain threshold, its rate of denudation increases dramatically and mesa degradation is completed in a short time. Copyright © 2014 John Wiley & Sons, Ltd.     

Impact of transient freshwater releases in the Southern Ocean on the AMOC and climate

The bipolar ocean seesaw is a process that explains the competition between deep waters formed in the North Atlantic (NA) and in the Southern Ocean (SO). In this picture, an increase in the rate of formation of one of these water masses is made at the expense of the other. However, recent studies have questioned the effectiveness of this process. Namely, they show that adding freshwater in the SO can reduce deep water formation in the SO as well as in the NA. In this study, we explore the mechanisms and time scales excited by such a SO freshwater release by performing sensitivity experiments where a freshwater input is added abruptly in the ocean, south of 60°S, with different rates and durations. For this purpose, we evaluate the separate effects of wind, temperature and salinity changes, and we put the emphasis on the time evolution of the system. We find three main processes that respond to these freshwater inputs and affect the NA Deep Water (NADW) production: (i) the deep water adjustment, which enhances the NADW cell, (ii) the salinity anomaly spread from the SO, which weakens the NADW cell, and (iii) the increase in the Southern Hemisphere wind stress, which enhances the NADW cell. We show that process (i) affects the Atlantic in a few years, due to an adjustment of the pycnocline depth through oceanic waves in response to the buoyancy perturbation in the SO. The salinity anomalies responsible for the NADW production decrease [process (ii)] invades the NA in around 30 years, while the wind stress from process (iii) increases in around 20 years after the beginning of the freshwater perturbation. Finally, by testing the response of the ocean to a large range of freshwater release fluxes, we show that for fluxes larger than 0.2 Sv, process (ii) dominates over the others and limits NADW production after a few centuries, while for fluxes lower than 0.2 Sv, process (ii) hardly affects the NADW production. On the opposite, the NADW export is increased by processes (i) and (iii) even for fluxes smaller than 0.1 Sv. The climatic impact of the freshwater release in the SO is mainly a cooling of the Southern Hemisphere, of up to 10°C regionally, which increases with freshwater release fluxes for a large range of values.     

A systematic study of the impact of freshwater pulses with respect to different geographical locations

The first comparative and systematic climate model study of the sensitivity of the climate response under Last Glacial Maximum (LGM) conditions to freshwater perturbations at various locations that are known to have received significant amounts of freshwater during the LGM (21 kyr BP) climate conditions is presented. A series of ten regions representative of those receiving most of the meltwater from decaying ice-sheets during the deglaciation is defined, comprising the border of LGM ice-sheets, outlets of rivers draining part of the melting ice-sheets and iceberg melt zones. The effect of several given freshwater fluxes applied separately in each of these regions on regional and global climate is subsequently tested. The climate response is then analysed both for the atmosphere and oceans. Amongst the regions defined, it is found that the area close by and dynamically upstream to the main deep water formation zone in the North Atlantic are most sensitive to freshwater pulses, as is expected. However, some important differences between Arctic freshwater forcing and Nordic Seas forcing are found, the former having a longer term response linked to sea-ice formation and advection whereas the latter exhibits more direct influence of direct freshening of the deep water formation sites. Combining the common surface temperature response for each respective zone, we fingerprint the particular surface temperature response obtained by adding freshwater in a particular location. This is done to examine if a surface climate response can be used to determine the origin of a meltwater flux, which is relevant for the interpretation of proxy data. We show that it is indeed possible to generally classify the fingerprints by their origin in terms of sea-ice modification and modification of deep-water formation. Whilst the latter is not an unambiguous characterization of each zone, it nonetheless provides important clues on the physical mechanisms at work. In particular, it is shown that in order to obtain a consistent see-saw temperature pattern, addition of freshwater in the Northern Hemisphere at sites dynamically close to the deep water formation zones is needed. Finally a preliminary data—model comparison for the time of the Heinrich event 1 suggests that those sites are indeed the most favourable to explain the pattern of climate variability recorded in proxy data for this period. More importantly, this model—data comparison enables us to clearly reject a substantial fraction of the zones tested as potential source for large freshwater entering the ocean at that time.     

New experimental determination of Li and B partition coefficients during upper mantle partial melting

Despite the growing interest for Li and B as geochemical tracers, especially for material transfer from subducting slabs to overlying peridotites, little is known about the behaviour of these two elements during partial melting of mantle sources. In particular, mineral/melt partition coefficients for B and to a lesser extent Li are still a matter of debate. In this work, we re-equilibrated a synthetic basalt doped with ~10 ppm B and ~6 ppm Li with an olivine powder from a spinel lherzolite xenolith at 1 GPa–1,330°C, and we analyzed Li and B in the run products by secondary ion mass spectrometry (SIMS). In our experiment, B behaved as a highly incompatible element with mineral/melt partition coefficients of the order of 10⁻² (D ^ol/melt = 0.008 (0.004–0.013); D ^opx/melt = 0.024 (0.015–0.033); D ^cpx/melt = 0.041 (0.021–0.061)), and Li as a moderately incompatible element (D ^ol/melt = 0.427 (0.418–0.436); D ^opx/melt = 0.211 (0.167–0.256); D ^cpx/melt = 0.246 (0.229–0.264)). Our partition coefficients for Li are in good agreement with previous determinations. In the case of B, our partition coefficients are equal within error to those reported by Brenan et al. (1998) for all the mineral phases analyzed, but are lower than other coefficients from literature for some of the phases (up to 5 times for cpx). Our measurements complement the data set of Ds for modelling partial melting of the upper mantle and basalt generation, and confirm that, in this context, B is more incompatible than previously anticipated.     

Long lasting interactions between tectonic loading, unroofing, post-rift thermal subsidence and sedimentary transfers along the western margin of the Gulf of Mexico: Some insights from integrated quantitative studies

After Jurassic rifting, numerous carbonate platforms (i.e., the Orizaba, Cordoba and Golden Lane-Tuxpan platforms) developed during the Lower and Middle Cretaceous episode of thermal subsidence along the western passive margin of the Gulf of Mexico, with intervening basinal domains (i.e., the Tampico-Misantla, Zongolica, Veracruz and Deep Gulf of Mexico - DGM - basins).During the Late Cretaceous-Paleocene, the east-verging Sierra Madre Oriental thrust belt developed, resulting in tectonic uplift and unroofing of the allochthonous units (i.e. tectonic units made up of former Orizaba and Cordoba platforms and Zongolica Basin series). This new topography provided also an important source of clastics to feed the adjacent foredeep, where coeval tectonic loading accounted for the bending of the foreland lithosphere. However, shallow water facies or even emersion persisted until the Eocene in the forebulge area (at the present location of the Golden Lane), preventing locally the clastics to reach the DGM. This topographic barrier was ultimately bypassed by the clastics only during the Oligocene and Neogene, once (1) the prograding clastic wedge had exceeded accommodation, and (2) the long lasting thermal subsidence of the passive margin could overpass the effect of the bending and force the former bulge to sink.Numerous paleo-thermo-meters (Tmax, Ro), paleo-thermo-barometers (fluid inclusions), PVT and coupled forward kinematic and thermal modeling have been used to calibrate and date the progressive unroofing of the thrust belt. Coupled tectonic and sedimentologic modeling was applied in the foreland to predict the distribution of sand versus shale ratios in the Oligocene to Plio-Quaternary clastic sedimentary wedge of the passive margin, where gravitational gliding of post-Eocene series occurred during the Neogene along major listric faults.Mantle dynamics are advocated as the main process accounting for post-orogenic uplift and regional tilting of the basement, which initiated a massive transfer of sediments from the Cordillera towards the Gulf of Mexico, from Oligocene onward, resulting in a destabilization and gravitational collapse of the western slope of the Gulf of Mexico in Neogene times.     

An archaeo-seismological study of the Nîmes Roman aqueduct,France: indirect evidence for an <Emphasis Type="Italic">M</Emphasis> > 6 seismic event?

This paper presents a synthesis of a multidisciplinary study carried out along the Nîmes Roman aqueduct, located in the southeast part of France. The study was motivated by archaeologists attempting to explain the partial destructions of only one aerial bridge of the aqueduct (Pont de la Lône). Given its close proximity to the Nîmes seismically active fault, a possible seismic origin for the destructions was invoked. Seismologists and structural engineers thus carried out a variety of field and numerical investigations to test the seismic hypothesis. Supporting field evidence was found first along the aerial bridge section of the aqueduct: broken stalactites, arch warping, cracks and destruction of the bridge just above a breast wall shortly after its construction. Secondly, the underground part of the canal was analysed: irregularity in the thickness of calcite deposits of the canal walls, presence of numerous cracks, horizontal shift of the otherwise linear structure of the canal and presence of calcite twins in the deposits, found only where the canal crosses the Nîmes fault system. Numerical modelling and experimental results show that (1) only the Pont de la Lône would have suffered serious damage under seismic solicitation (assuming an M6 earthquake at 10 km distance). The second aerial bridge, the three levels arches Pont du Gard, an historical monument still standing today, would have suffered less damage, due to its very different fundamental frequency of around 0.4 Hz, far from the amplified seismic near-field spectrum; (2) the numerical models also show that a fall of the canal wall would require a higher magnitude event (M > 6); (3) the presence of calcite twinning requires a differential static stress of 4 MPa or greater, which can only be achieved very close to the fault rupture of a M > 6 earthquake; (4) finally, local surface rupturing of such a fault would also corroborate the hypothesis that the observed offset of the canal may be partly seismically induced. Although other possible origins for each individual evidence may not be excluded, the observed spatiotemporal concentration of architectural/geological anomalies together with the numerical results allow us to support a possible co-seismic origin for these disorders, indirectly attesting to the potential seismic activity (M > 6) of the nearby Nîmes fault. Furthermore, following the conclusions of archaeological studies, the disorders occurred between 250 and 350 year AD, thus providing timing for this possible seismic event, an essential parameter in seismic hazard analysis. These results validate the archaeo-seismological approach as a tool that may help improve the knowledge of major infrequent earthquakes in areas of moderate seismic activity.     

Spectral Simulation of Isotropic Gaussian Random Fields on a Sphere

A spectral algorithm is proposed to simulate an isotropic Gaussian random field on a sphere equipped with a geodesic metric. This algorithm supposes that the angular power spectrum of the covariance function is explicitly known. Direct analytic calculations are performed for exponential and linear covariance functions. In addition, three families of covariance functions are presented where the calculation of the angular power spectrum is simplified (shot-noise random fields, Yadrenko covariance functions and solutions of certain stochastic partial differential equations). Numerous illustrative examples are given.