Soil discrimination using diffuse reflectance Vis–NIR spectroscopy in a local toposequence

Vis–NIR spectroscopy is nowadays presented as a possible routine method for soil sample analysis. However, there is still no consensus on which is the best multivariate statistical method to use. We propose to use principal component analysis to complete the spectral data treatment. The soil samples came from a pedological cover made up of red–yellow Latosols: 88 samples of 11 soil profiles on four toposequences were collected; clay, organic matter, silica, iron, aluminum and titanium total contents were determined; the contents of goethite, hematite, gibbsite, and kaolinite were calculated. Diffuse reflectance Vis–NIR spectroscopy at wavelengths from 400 to 2400 nm combined with principal component analysis (PCA) was sufficiently sensitive to discriminate different Latosols. Wavelengths of 700 nm and 2200 to 2300 nm were influenced by content ratios of organic matter and iron oxides (700 nm), and kaolinite and gibbsite absorption (2200 and 2300 nm). The spectral responses were affected not only by the content of these constituents, but also by the composition of the minerals, so that the same class of Latosol may have different or similar spectral responses. The role of microaggregation is discussed.     

Kinematics of the Malatya–Ovacik Fault Zone

Abstract

We investigate the left-lateral slip on the 240-km- long, NE-SW-trending, Malatya-Ovacik fault zone in eastern Turkey. This fault zone splays southwestward from the North Anatolian fault zone near Erzincan, then follows the WSW-trending Ovacik valley between the Munzur and Yilan mountain ranges. It bends back to a SW orientation near Arapkir, from where we trace its main strand SSW beneath the Plio-Quaternary sediment of the Malatya basin. We propose that this fault zone was active during ~5–3 Ma, when it took up 29 km of relative motion between the Turkish and Arabian plates; it ceased to be active when the East Anatolian fault zone formed at ~3 Ma. The geometry of the former Erzincan triple junction, which differs from the modem Karliova triple junction, where the North and East Anatolian fault zones intersect, suggests a possible explanation for why slip on the Malatya- Ovacik fault zone was unable to continue. We interpret the SW- and SSW-trending segments of the Malatya-Ovacik fault zone as transform faults, which define an Euler pole ~1 400 km to the southeast. Its central part along the Ovacik valley, which is ~30° oblique to the adjoining transform faults, is interpreted as the internal fault of a stepover. The adjoining mountain ranges, which now rise up to ~3 300 m, ~2 000 m above the surrounding land surface, are largely the result of the surface uplift which accompanied the components of shortening and thickening of the upper crustal brittle layer that occurred around this stepover while the left-lateral faulting was active. © 2001 Éditions scientifiques et médicales Elsevier SAS     

Small valley glaciers and the effectiveness of the glacial buzzsaw in the northern Basin and Range, USA

The glacial buzzsaw hypothesis suggests that efficient erosion limits topographic elevations in extensively glaciated orogens. Studies to date have largely focussed on regions where large glaciers (tens of kilometres long) have been active. In light of recent studies emphasising the importance of lateral glacial erosion in lowering peaks and ridgelines, we examine the effectiveness of small glaciers in limiting topography under both relatively slow and rapid rock uplift conditions. Four ranges in the northern Basin and Range, Idaho, Montana, and Wyoming, USA, were chosen for this analysis. Estimates of maximum Pleistocene slip rates along normal faults bounding the Beaverhead–Bitterroot Mountains (~ 0.14 mm y^− 1), Lemhi Range (~ 0.3 mm y^− 1) and Lost River Range (~ 0.3 mm y^− 1) are an order of magnitude lower than those on the Teton Fault (~ 2 mm y^− 1). We compare the distribution of glacial erosion (estimated from cirque floor elevations and last glacial maximum (LGM) equilibrium line altitude (ELA) reconstructions) and fault slip rate with three metrics of topography in each range: the along-strike maximum elevation swath profile, hypsometry, and slope-elevation profiles. In the slowly uplifting Beaverhead–Bitterroot Mountains, and Lemhi and Lost River Ranges, trends in maximum elevation parallel ELAs, independent of variations in fault slip rate. Maximum elevations are offset ~ 500 m from LGM ELAs in the Lost River Range, Lemhi Range, and northern Beaverhead–Bitterroot Mountains, and by ~ 350 m in the southern Beaverhead–Bitterroot Mountains, where glacial extents were less. The offset between maximum topography and mean Quaternary ELAs, inferred from cirque floor elevations, is ~ 350 m in the Lost River and Lemhi Ranges, and 200–250 m in the Beaverhead–Bitterroot Mountains. Additionally, slope-elevation profiles are flattened and hypsometry profiles show a peak in surface areas close to the ELA in the Lemhi Range and Beaverhead–Bitterroot Mountains, suggesting that small glaciers efficiently limit topography. The situation in the Lost River Range is less clear as a glacial signature is not apparent in either slope-elevation profiles or the hypsometry. In the rapidly uplifting Teton Range, the distribution of ELAs appears superficially to correspond to maximum topography, hypsometry, and slope-elevations profiles, with regression lines on maximum elevations offset by ~ 700 and ~ 350 m from the LGM and mean Quaternary ELA respectively. However, Grand Teton and Mt. Moran represent high-elevation “Teflon Peaks” that appear impervious to glacial erosion, formed in the hard crystalline bedrock at the core of the range. Glacier size and drainage density, rock uplift rate, and bedrock lithology are all important considerations when assessing the ability of glaciers to limit mountain range topography. In the northern Basin and Range, it is only under exceptional circumstances in the Teton Range that small glaciers appear to be incapable of imposing a fully efficient glacial buzzsaw, emphasising that high peaks represent an important caveat to the glacial buzzsaw hypothesis.     

Calibration of decomposition of serine to alanine in Bithynia opercula as a quantitative dating technique for Middle and Late Pleistocene sites in Britain

The concentration ratio of the amino acids alanine and serine in calcite opercula of freshwater gastropods of genus Bithynia has been assessed as an age proxy for British Middle and Late Pleistocene sites. Due to decomposition of peptide-bound serine within intra-crystalline protein molecules, alanine/serine ratios increase with sample age, enabling resolution of the marine oxygen-isotope stage of almost all sites and substage resolution for some sites. However, slight overlaps are evident in alanine/serine ratios between some sites from late in one temperate stage and early in the next; this effect results from the low decomposition rates during cold stages, combined with an evident slight dependence of decomposition rates on site conditions. The technique, which can be extended to the Early Pleistocene without technical refinements, is amenable to quantitative calibration using a first-order rate equation, making it feasible in future to present results as numerical ages, subject to an assumed temperature history. Most Pleistocene sites investigated are shown using this technique to have the same ages as previously deduced using biostratigraphy. However, the serine–alanine ages obtained preclude glaciation of southeast England during MIS 16 or of central England during MIS 10.     

Channel Representation in Physically Based Models Coupling Groundwater and Surface Water: Pitfalls and How to Avoid Them

Recent models that couple three‐dimensional subsurface flow with two‐dimensional overland flow are valuable tools for quantifying complex groundwater/stream interactions and for evaluating their influence on watershed processes. For the modeler who is used to defining streams as a boundary condition, the representation of channels in integrated models raises a number of conceptual and technical issues. These models are far more sensitive to channel topography than conventional groundwater models. On all spatial scales, both the topography of a channel and its connection with the floodplain are important. For example, the geometry of river banks influences bank storage and overbank flooding; the slope of the river is a primary control on the behavior of a catchment; and at the finer scale bedform characteristics affect hyporheic exchange. Accurate data on streambed topography, however, are seldom available, and the spatial resolution of digital elevation models is typically too coarse in river environments, resulting in unrealistic or undulating streambeds. Modelers therefore perform some kind of manual yet often cumbersome correction to the available topography. In this context, the paper identifies some common pitfalls, and provides guidance to overcome these. Both aspects of topographic representation and mesh discretization are addressed. Additionally, two tutorials are provided to illustrate: (1) the interpolation of channel cross‐sectional data and (2) the refinement of a mesh along a stream in areas of high topographic variability.     

Forecast of Natural Aquifer Discharge Using a Data‐Driven,Statistical Approach

In the Western United States, demand for water is often out of balance with limited water supplies. This has led to extensive water rights conflict and litigation. A tool that can reliably forecast natural aquifer discharge months ahead of peak water demand could help water practitioners and managers by providing advanced knowledge of potential water‐right mitigation requirements. The timing and magnitude of natural aquifer discharge from the Eastern Snake Plain Aquifer (ESPA) in southern Idaho is accurately forecast 4 months ahead of the peak water demand, which occurs annually in July. An ARIMA time‐series model with exogenous predictors (ARIMAX model) was used to develop the forecast. The ARIMAX model fit to a set of training data was assessed using Akaike's information criterion to select the optimal model that forecasts aquifer discharge, given the previous year's discharge and values of the predictor variables. Model performance was assessed by application of the model to a validation subset of data. The Nash‐Sutcliffe efficiency for model predictions made on the validation set was 0.57. The predictor variables used in our forecast represent the major recharge and discharge components of the ESPA water budget, including variables that reflect overall water supply and important aspects of water administration and management. Coefficients of variation on the regression coefficients for streamflow and irrigation diversions were all much less than 0.5, indicating that these variables are strong predictors. The model with the highest AIC weight included streamflow, two irrigation diversion variables, and storage.     

A skeleton design theory for spatial data infrastructure

In this paper, we look into the theory of designing geoservice systems, i.e., SDI networks and their constituent SDI nodes. As the field of SDI is strongly about bridging between geoservice systems, interoperability and harmonisation, it is not surprising that standardisation efforts are of crucial importance in it. These efforts have historically addressed abstract and concrete content models for data and metadata exchange, as well as abstract and concrete behavioural models for computational processes. The list of standards that are in use in the SDI field continues to expand, and reaches out to neighbouring fields such as sensor nets. We argue that given these trends, the resulting levels of standardisation in actual systems, and the complexity of geoservice systems in general, it appears only natural to look into the possibility to define a standardised design theory for SDI and its nodes, which addresses the function base and the communication base. Specifically, we provide an overview of those components that need to be designed, and what are their relationships. We do so in an abstract way, focussing on the concern of information content in this paper, and only hinting at an appropriate theory of realisation based on our skeleton theory.     

Late Cenozoic surface uplift,basaltic volcanism,and incision by the River Tigris around Diyarbakır,SE Turkey

We document the staircase of terraces of the River Tigris in the Diyarbakır area of SE Turkey, in the northern Arabian Platform, and improve control on the ages of these terrace deposits by dating of overlying basalt flows using the unspiked K–Ar technique. These fluvial terraces are formed of polymict gravel, including clasts derived from the Anatolian metamorphic terrane farther north as well as of local basalt. At least 9 Tigris terraces have been recognised so far, the highest of which, ∼200 m above present river level, marks the local transition from stacked deposition to fluvial incision, the timing of which is bounded between the mid Late Miocene and the Middle Pliocene. Our K–Ar dating indicates a hiatus in fluvial incision in the late Early Pleistocene, as basalts dated to 1.22 ± 0.02 and 1.07 ± 0.03 Ma overlie Tigris gravels at very similar levels, ∼60–70 m above the present river. The lower terraces record the subsequent entrenchment of the modern Tigris valley following an increase in incision rates in the early Middle Pleistocene, evident from the disposition of younger basalt, dated to 0.43 ± 0.02 Ma, capping fluvial gravel only ∼21–22 m above the present river level. Numerical modelling can account for the observed uplift history, as the response to coupling between surface processes and induced flow in the lower crust, with the mobile lower-crust thin (∼5–7 km thick), consistent with the known presence of a thick layer of mafic underplating at the base of the crust beneath the Arabian Platform. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Regional isostatic response to Messinian Salinity Crisis events

The salinity crisis of the Mediterranean during Messinian time was one of the most dramatic episodes of oceanic change of the past 20 or so million years, resulting in the deposition of kilometer thick evaporitic sequences. A large and rapid drawdown of the Mediterranean water level caused erosion and deposition of non-marine sediments in a large ‘Lago Mare’ basin. Both the surface loading by the Lower Messinian evaporites, and the removal of the water load resulted in isostatic/flexural rebound that significantly affected river canyons and topographic slopes. We use flexure models to quantitatively predict possible signatures of these events, and verify these expectations at well-studied margins. The highly irregular shape of the reconstructed basin calls for a three-dimensional model. Near basin margins, plate-bending effects are most pronounced which is why flexure is particularly important for a relatively narrow basin like the Mediterranean. We focus on one specific sea level scenario for the Messinian Salinity Crisis, where most of the evaporite load was deposited during a sea level highstand, followed by a rapid desiccation. Evaporite loading at current sea level is expected to cause subsidence of the deep basins by hundreds of meters and simultaneous uplift of continental parts of the margins. Differential uplift may lead to significant slope angle changes and thus gravity flows. The relative scarcity of Lower Evaporite sequences along the margins may be a result of these phenomena. Normal faulting of Lower Evaporite and older sediments and rocks is expected on the margins. Desiccation enhances erosion of the freshly exposed continental shelf and slope. Subsidence and riverbed sedimentation occurs on the continental margins, and significant uplift towards the basin center. Reverse faulting is predicted at the margins. Finally, regional isostatic uplift following Zanclean flooding is predicted to destabilize margin slope deposits, and to cause marginal uplift, river down-cutting, and normal faulting.