A 2300-year record of environmental change from SW Anatolia, Lake Burdur, Turkey

This work deals with changes recorded by lacustrine sediments from SW Anatolia in Turkey in the context of increasing stress on the Mediterranean environment in relation to human-climate-environmental interactions. Paleolimnological investigations were carried out on Lake Burdur (Lake District geographical subregion of Turkey), which has been subject to rapid changes in its hydrological system that caused a ~10-m water-level drop in the last 30 years. Study of a 5-m-long sediment core, taken from the recently dried out part of the lake, shows significant variation in granulometry, clay mineralogy, nitrogen and organic carbon content and its isotopic composition (?¹³C_org) throughout the sedimentary sequence, which represents the last 2,300 years. Chronology is based on radiocarbon dates obtained from plant fragments. The results provide a record of environmental changes, including biological productivity and erosion intensity, in relation to changes in water level and humidity of the area. The lowest previous water level, which is the same as observed today, occurred at approximately 300 BC, after which time a water-level increase led to a maximum stage between AD 200 and 1200. From AD 1200 to the present, the climate became drier accompanied by an overall lowering of water level, with periods of water-level fluctuations superimposed. The lake never dried out as is observed today at the coring site. The present drop in the lake water level is attributed to human activity rather than to climate change. This record of climate and environmental change in the Lake Burdur area appears to be stratigraphically complete. This study provides rare data for the late Holocene in Anatolia as well as information about rapid climate changes during specific periods. The influence of both the Northern hemisphere atmospheric circulation and Indian monsoon on the east Mediterranean remains unclear. Nevertheless, the North Sea—Caspian atmospheric teleconnection could be an important factor that explains local differences in climatic evolution of Anatolia during the period considered.     

Daily streamflow modelling and assessment based on the curve‐number technique

A cell‐based long‐term hydrological model (CELTHYM) that can be integrated with a geographical information system (GIS) was developed to predict continuous stream flow from small agricultural watersheds. The CELTHYM uses a cell‐by‐cell soil moisture balance approach. For surface runoff estimation, the curve number technique considering soil moisture on a daily basis was used, and release rate was used to estimate baseflow. Evapotranspiration was computed using the FAO modified Penman equation that considered land‐use‐based crop coefficients, soil moisture and the influence of topography on radiation. A rice paddy field water budget model was also adapted for the specific application of the model to East Asia. Model sensitivity analysis was conducted to obtain operational information about the model calibration parameters. The CELTHYM was calibrated and verified with measured runoff data from the WS#1 and WS#3 watersheds of the Seoul National University, Department of Agricultural Engineering, in Hwaseong County, Kyounggi Province, South Korea. The WS#1 watershed is comprised of about 35·4% rice paddy fields and 42·3% forest, whereas the WS#3 watershed is about 85·0% forest and 11·5% rice paddy fields. The CELTHYM was calibrated for the parameter release rate, K, and soil moisture storage coefficient, STC, and results were compared with the measured runoff data for 1986. The validation results for WS#1 considering all daily stream flow were poor with R², E² and RMSE having values of 0·40, ?6·63 and 9·69 (mm), respectively, but validation results for days without rainfall were statistically significant (R² = 0·66). Results for WS#3 showed good agreement with observed data for all days, and R², E² and RMSE were 0·92, 0·91 and 2·23 (mm), respectively, suggesting potential for CELTHYM application to other watersheds. The direct runoff and water balance components for watershed WS#1 with significant areas of paddy fields did not perform well, suggesting that additional study of these components is needed. Copyright © 2002 John Wiley & Sons, Ltd.     

4. Building of a Habitable Planet

Except the old Jack Hills zircon crystals, it does not exit direct record of the first 500 Ma of the Earth history. Consequently, the succession of events that took place during this period is only indirectly known through geochemistry, comparison with other telluric planets, and numerical modelling. Just after planetary accretion several episodes were necessary in order to make life apparition and development possible and to make the Earth surface habitable. Among these stages are: the core differentiation, the formation of a magma ocean, the apparition of the first atmosphere, oceans and continents as well as the development of magnetic field and of plate tectonics. In the same time, Earth has been subject to extraterrestrial events such as the Late Heavy Bombardment (LHB) between 3.95 and 3.8 Ga. Since 4.4–4.3 Ga, the conditions for pre-biotic chemistry and appearance of life were already met (liquid water, continental crust, no strong meteoritic bombardment, etc...). This does not mean that life existed as early, but this demonstrates that all necessary conditions assumed for life development were already present on Earth.     

Noble gases in the Martian meteorite Northwest Africa 2737: A new chassignite signature

Abstract— We report noble gas data for the second chassignite, Northwest Africa (NWA) 2737, which was recently found in the Moroccan desert. The cosmic ray exposure (CRE) age based on cosmogenic ³He, ²¹Ne, and ³⁸Ar around 10–11 Ma is comparable to the CRE ages of Chassigny and the nakhlites and indicates ejection of meteorites belonging to these two families during a discrete event, or a suite of discrete events having occurred in a restricted interval of time. In contrast, U‐Th/He and K/Ar ages <0.5 Ga are in the range of radiometric ages of shergottites, despite a Sm‐Nd signature comparable to that of Chassigny and the nakhlites (Misawa et al. 2005). Overall, the noble gas signature of NWA 2737 resembles that of shergottites rather than that of Chassigny and the nakhlites: NWA 2737 does not contain, in detectable amount, the solar‐like xenon found in Chassigny and thought to characterize the Martian mantle nor apparently fission xenon from ²⁴⁴Pu, which is abundant in Chassigny and some of the nakhlites. In contrast, NWA 2737 contains Martian atmospheric noble gases trapped in amounts comparable to those found in shergottite impact glasses. The loss of Martian mantle noble gases, together with the trapping of Martian atmospheric gases, could have occurred during assimilation of Martian surface components, or more likely during shock metamorphism, which is recorded in the petrology of this meteorite.     

Multibeam bathymetry and sidescan imaging of the Rivera Transform–Moctezuma Spreading Segment junction, northern East Pacific Rise: New constraints on Rivera–Pacific relative plate motion

To better understand the recent motion of the Pacific plate relative to the Rivera plate and to better define the limitations of the existing Rivera–Pacific plate motion models for accurately predicting this motion, total-field magnetic data, multibeam bathymetric data and sidescan sonar images were collected during the BART and FAMEX campaigns of the N/O L'Atalante conducted in April and May 2002 in the area surrounding the Moctezuma Spreading Segment of the East Pacific Rise, located offshore of Manzanillo, Mexico, at 106°16′W, between 17.8°N and 18.5°N. Among the main results are: (1) the principle transform displacement zone of the Rivera Transform is narrow and well defined east of 107o15′W and these azimuths should be used preferentially when deriving new plate motion models, and (2) spreading rates along the Moctezuma Spreading Segment should not be used in plate motion studies as either seafloor spreading has been accommodated at more than one location since the initiation of seafloor spreading in the area of the Moctezuma Spreading Segment, or this spreading center is not a Rivera–Pacific plate boundary as has been previously assumed. Comparison of observed transform azimuths with those predicted by the best-fit poles of six previous models of Rivera–Pacific relative motion indicate that, in the study area, a significant systematic bias is present in the predictions of Rivera–Pacific motion. Although the exact source of this bias remains unclear, this bias indicates the need to derive a new Rivera–Pacific relative plate motion model.     

The Grader layered intrusion (Havre-Saint-Pierre Anorthosite, Quebec) and genesis of nelsonite and other Fe–Ti–P ores

The Grader layered intrusion is part of the Havre-Saint-Pierre anorthosite in the Grenville Province (Quebec, Canada). This intrusion has a basin-like morphology and contains significant resources of Fe–Ti–P in ilmenite and apatite. Outcropping lithologies are massive oxide alternating with anorthosite layers, banded ilmenite–apatite–plagioclase rocks and layered oxide apatite (gabbro-)norites. Drill cores provide evidence for stratigraphic variations of mineral and whole rock compositions controlled by fractional crystallization with the successive appearance of liquidus phases: plagioclase and ilmenite followed by apatite, then orthopyroxene together with magnetite, and finally clinopyroxene. This atypical sequence of crystallization resulted in the formation of plagioclase–ilmenite–apatite cumulates or “nelsonites” in plagioclase-free layers. Fine-grained ferrodiorites that cross-cut the cumulates are shown to be in equilibrium with the noritic rocks. The high TiO₂ and P₂O₅ contents of these assumed liquids explains the early saturation of ilmenite and apatite before Fe–Mg silicates, thus the nelsonites represent cumulates rather than crystallized Fe–Ti–P-rich immiscible melts. The location of the most evolved mineral and whole rock compositions several tens of meters below the top of the intrusion, forming a sandwich horizon, is consistent with crystallization both from the base and top of the intrusion. The concentrations of V and Cr in ilmenite display a single fractionation path for the different cumulus assemblages and define the cotectic proportion of ilmenite to 21 wt.%. This corresponds to bulk cotectic cumulates with ca. 8 wt.% TiO₂, which is significantly lower than what is commonly observed in the explored portion of the Grader intrusion. The proposed mechanism of ilmenite-enrichment is the lateral removal of plagioclase due to its relative buoyancy in the dense ferrodiorite melt. This plagioclase has probably accumulated in other portions of the intrusion or has not been distinguished from the host anorthosite.     

Effects of heterogeneous land use/cover types on river channel morphology in the Solo River catchment,Eastern Uganda

In the tropics, unmonitored land use/cover types cause significant effects on the narrowing and widening of river channels which affects the integrity of water resources. River channel planform extent was characterized using Landsat images, while water and bedload samples were collected and analysed for a period of one year. The results revealed that in 1986, the channel planform covered 3.7 sq km in length than in 2013 where it increased to 4.2 sq km. Wetland (537.1mgl^?1) and bushland (186.3mgl^?1) cover types had the highest concentration of suspended sediments. Fine sand (0.25 mm), silty sand (1 mm) and silty clay (0.125 mm) bedload particle types dominated the riverbed along the channel from the sampled land use/cover types. The high concentration of sediments, bedload materials, bank instability, and streamflow were significant contributors to the narrowing and widening of the channel (p < 0.05). Agricultural land use was the major contributor of channel aggradation (0.8 m) and degradation (0.25 m) compared to tree plantations, bushlands, forest and wetland cover types.     

Evidence of quaternary transtensional tectonics in the Nekor basin (NE Morocco)

The geodynamic processes in the western Mediterranean are driven by both deep (mantle) processes such as slab‐rollback or delamination, oblique plate convergence and inherited structures. The present‐day deformation of the Alboran Sea and in particular the Nekor basin area is linked to these coeval effects. The seismically active Nekor basin is an extensional basin formed in a convergent setting at the eastern part of the Rif Chain whose boundaries extend both onshore and offshore Morocco. We propose a new structural model of the Nekor basin based on high‐resolution offshore data compiled from recent seismic reflection profiles, swath bathymetry acquisitions and industrial seismic reflection profiles. The new data set shows that the northern limit of the basin is oriented N49° with right‐stepping faults from the Bousekkour–Aghbal fault to the sinistral Bokkoya fault zone. This pattern indicates the presence of an inherited left‐lateral basement fault parallel to the major inherited Nekor fault. This fault has been interpreted as a Quaternary active left‐lateral transfer fault localized on weak structural discontinuities inherited from the orogenic period. Onshore and offshore active faults enclose a rhombohedral tectonic Nekor Basin. Normal faults oriented N155° offset the most recent Quaternary deposits in the Nekor basin, and indicate the transtensional behaviour of this basin. The geometry of these faults suggests a likely rollover structure and the presence at depth of a crustal detachment. Inactive Plio‐Quaternary normal faults to the east of the Ras Tarf promontory and geometries of depocentres seem to indicate the migration of deformation from east to west. The local orientations of horizontal stress directions deduced from normal fault orientations are compatible with the extrusion of the Rifian units and coherent with the westward rollback of the Tethyan slab and the localization of the present‐day slab detachment or delamination.     

Computation and visualisation of the accuracy of old maps using differential distortion analysis

The accuracy of old maps can hold interesting historical information, and is therefore studied using distortion analysis methods. These methods start from a set of ground control points that are identified both on the old map and on a modern reference map or globe, and conclude with techniques that compute and visualise distortion. Such techniques have advanced over the years, but leave room for improvement, as the current ones result in approximate values and a coarse spatial resolution. We propose a more elegant and more accurate way to compute distortion of old maps by translating the technique of differential distortion analysis, used in map projection theory, to the setting where an old map and a reference map are directly compared. This enables the application of various useful distortion metrics to the study of old maps, such as the area scale factor, the maximum angular distortion and the Tissot indicatrices. As such a technique is always embedded in a full distortion analysis method we start by putting forward an optimal analysis method for a general-purpose study, which then serves as the foundation for the development of our technique. Thereto, we discuss the structure of distortion analysis methods and the various options available for every step of the process, including the different settings in which the old map can be compared to its modern counterpart, the techniques that can be used to interpolate between both, and the techniques available to compute and visualise the distortion. We conclude by applying our general-purpose method, including the differential distortion analysis technique, to an example map also used in other literature.     

A stochastic inversion workflow for monitoring the distribution of CO<Subscript>2</Subscript> injected into deep saline aquifers

We present a two-step stochastic inversion approach for monitoring the distribution of CO₂ injected into deep saline aquifers for the typical scenario of one single injection well and a database comprising a common suite of well logs as well as time-lapse vertical seismic profiling (VSP) data. In the first step, we compute several sets of stochastic models of the elastic properties using conventional sequential Gaussian co-simulations (SGCS) representing the considered reservoir before CO₂ injection. All realizations within a set of models are then iteratively combined using a modified gradual deformation algorithm aiming at reducing the mismatch between the observed and simulated VSP data. In the second step, these optimal static models then serve as input for a history matching approach using the same modified gradual deformation algorithm for minimizing the mismatch between the observed and simulated VSP data following the injection of CO₂. At each gradual deformation step, the injection and migration of CO₂ is simulated and the corresponding seismic traces are computed and compared with the observed ones. The proposed stochastic inversion approach has been tested for a realistic, and arguably particularly challenging, synthetic case study mimicking the geological environment of a potential CO₂ injection site in the Cambrian-Ordivician sedimentary sequence of the St. Lawrence platform in Southern Québec. The results demonstrate that the proposed two-step reservoir characterization approach is capable of adequately resolving and monitoring the distribution of the injected CO₂. This finds its expression in optimized models of P- and S-wave velocities, density, and porosity, which, compared to conventional stochastic reservoir models, exhibit a significantly improved structural similarity with regard to the corresponding reference models. The proposed approach is therefore expected to allow for an optimal injection forecast by using a quantitative assimilation of all available data from the appraisal stage of a CO₂ injection site.