THE RELATIVE CONTRIBUTION OF SOIL TILLAGE AND OVERLAND FLOW EROSION TO SOIL REDISTRIBUTION ON AGRICULTURAL LAND

This study uses evidence for the long-term (35 years) pattern of soil redistribution within two agricultural fields in the UK to identify the relative importance of tillage and overland flow erosion. Spatially distributed long-term total soil redistribution data for the fields (Dalicott Farm and Rufford Forest Farm) were obtained using the caesium-137 (¹³⁷Cs) technique. These data were compared with predicted patterns of soil redistribution. Recent studies have demonstrated that the redistribution of soil by tillage may be described as a diffusive process. A two-component model was, therefore, developed which accounts for soil redistribution by both overland flow and diffusive processes. Comparison of the predicted patterns of overland flow erosion alone with the observed (¹³⁷Cs-derived) data indicated a poor agreement (r² = 0.17 and 0.11). In contrast, a good agreement exists between the predicted pattern of diffusive redistribution and the observed data (r² = 0.43 and 0.41). These results give a clear indication that diffusive processes are dominant in soil redistribution within these fields. Possible diffusive processes include splash erosion, soil creep and tillage. However, the magnitude of the diffusion coefficients for the optimum predicted pattern (c. 350–400 kg m⁻¹ a⁻¹) demonstrates that tillage is the only process capable of explaining the very significant soil redistribution which is indicated by the ¹³⁷Cs data. Consideration is given to the implications of these results for both soil erosion prediction and landscape interpretation.     

Quantifying preferential flow in soils: A review of different techniques

Preferential flow (PF) in soil has both environmental and human health implications since it favours contaminant transport to groundwater without interaction with the chemically and biologically reactive upper layer of soil. PF is, however, difficult to measure and quantify. This paper reviews laboratory and field techniques, such as breakthrough curves, dye tracing, and scanning techniques, for evaluating PF in soil at different scales. Advanced technologies, such as scanning techniques, have increased our capability to quantify transport processes within the soil with minimal soil disturbance. Important issues with respect to quantifying PF concern large-scale studies, frozen soil conditions, tracing techniques for particles and gases, a lack of simple mathematical tools for interpreting field data, and the lack of a systematic approach for comparing PF data resulting from different measurement techniques. Also, more research is required to quantify the relative importance of the various PF processes that occur in soil rather than the integrated result of all PF processes in soils.     

Acigöl rhyolite field,central Anatolia (part II): geochemical and isotopic (Sr–Nd–Pb, δ<Superscript>18</Superscript>O) constraints on volcanism involving two high-silica rhyolite suites

The Acigöl rhyolite field erupted the most recent high-silica rhyolites within the Cappadocian Volcanic Province of central Anatolia, Turkey. It comprises two sequences of domes and pyroclastic rocks with eruption ages of ~150–200 ka (eastern group) and ~20–25 ka (western group). Compositionally, the eastern rhyolite group lavas are less evolved (SiO₂ = 74–76 wt%), whereas the western group has higher silica abundance (SiO₂ = ~77 wt%) with extremely depleted feldspar-compatible trace elements. Within each group, compositional variability is small and ¹⁴³Nd/¹⁴⁴Nd (0.51257–0.51265) and Pb isotope compositions (²⁰⁶Pb/²⁰⁴Pb = 18.87–18.88, ²⁰⁷Pb/²⁰⁴Pb = 15.65–15.67 and ²⁰⁸Pb/²⁰⁴Pb = 38.94–38.98) are homogeneous. The western group rhyolites have δ¹⁸O(zircon) overlapping mantle values (5.7 ± 0.2‰), whereas eastern group rhyolites are enriched in δ¹⁸O by ~0.5‰, consistent with a tendency to lower εNd values. By contrast, western group rhyolites have markedly more radiogenic ⁸⁷Sr/⁸⁶Sr ratios (0.7065–0.7091) compared to those of the eastern group (0.7059–0.7065). The presence of angular granitic xenoliths and a correlation between hydration (based on loss on ignition data) and ⁸⁷Sr/⁸⁶Sr in the western lavas, however, indicates that Sr was added during the eruption or post-eruption alteration. Isotope constraints preclude the possibility that the rhyolite magmas formed by partial melting of any known regional crystalline basement rocks. Basalts and andesites erupted in the periphery of the Acigöl field are characterised by ⁸⁷Sr/⁸⁶Sr ratios between 0.7040 and 0.7053, ¹⁴³Nd/¹⁴⁴Nd = 0.51259–0.51300, ²⁰⁶Pb/²⁰⁴Pb = 18.85–18.87, ²⁰⁷Pb/²⁰⁴Pb = 15.646–15.655, ²⁰⁸Pb/²⁰⁴Pb = 38.90–38.97. The isotopic and trace element data favour an origin of the rhyolites by mixing of basaltic/andesitic magmas with minor amounts of crustal melts and followed by extensive fractional crystallization.     

GROUND-WATER DISPERSION CALCULATIONS WITH A PROGRAMMABLE CALCULATOR

Abstract. A simple and convenient method for calculating two-dimensional dispersion on a TI 58 or 59 programmable calculator utilizing the normal distribution program available in the Solid State Software Master Library module as a subroutine is presented with an example.     

Divergence and flow direction as indicators of subsurface heterogeneity and stage‐dependent storage in alluvial floodplains

Assuming homogeneity in alluvial aquifers is convenient, but limits our ability to accurately predict stream‐aquifer interactions. Research is needed on (i) identifying the presence of focused, as opposed to diffuse, groundwater discharge/recharge to streams and (ii) the magnitude and role of large‐scale bank and transient storage in alluvial floodplains relative to changes in stream stage. The objective of this research was to document and quantify the effect of stage‐dependent aquifer heterogeneity and bank storage relative to changes in stream stage using groundwater flow divergence and direction. Monitoring was performed in alluvial floodplains adjacent to the Barren Fork Creek and Honey Creek in northeastern Oklahoma. Based on results from subsurface electrical resistivity mapping, observation wells were installed in high and low electrical resistivity subsoils. Water levels in the wells were recorded real time using pressure transducers (August to October 2009). Divergence was used to quantify heterogeneity (i.e. variation in hydraulic conductivity, porosity, and/or aquifer thickness), and flow direction was used to assess the potential for large‐scale (100 m) bank or transient storage. Areas of localized heterogeneity appeared to act as divergence zones allowing stream water to quickly enter the groundwater system, or as flow convergence zones draining a large groundwater area. Maximum divergence or convergence occurred with maximum rates of change in flow rates or stream stage. Flow directions in the groundwater changed considerably between base and high flows, suggesting that the floodplains acted as large‐scale bank storage zones, rapidly storing and releasing water during passage of a storm hydrograph. During storm events at both sites, the average groundwater direction changed by at least 90° from the average groundwater direction during baseflow. Aquifer heterogeneity in floodplains yields hyporheic flows that are more responsive and spatially and temporally complex than would be expected compared to more common assumptions of homogeneity. Copyright © 2012 John Wiley & Sons, Ltd.     

Circulation in the southwestern part of the Kara Sea in September 2007

During cruise 54 of the R/V Akademik Mstislav Keldysh to the southwestern Kara Sea (September 6 to October 7, 2007), a large amount of hydrophysical data with unique spatial resolution was obtained on the basis of measurements using different instruments. The analysis of the data gave us the possibility to study the dynamics and hydrological structure of the southwestern Kara Sea basin. The main elements of the general circulation are the following: the Yamal Current, the Eastern Novaya Zemlya Current, and the St. Anna Trough Current. All these currents are topographically controlled; they flow over the bottom slopes along the isobaths. The Yamal Current begins at the Kara Gates Strait and turns to the east as part of the cyclonic circulation. Then, it turns to the north and propagates along the Yamal coast over the 100-m isobath. The Eastern Novaya Zemlya Current (its core is located over the eastern slope of the Novaya Zemlya Trough) flows to the northeast. Near the northern edge of Novaya Zemlya, it encounters the St. Anna Trough Current, separates from the coast, and flows practically to the east merging with the continuation of the Yamal Current. A strong frontal zone is formed in the region where the two currents merge above the threshold that separates the St. Anna Trough from the Novaya Zemlya Trough and divides the warm and saline Arctic waters from the cooler and fresher waters of the southwestern part of the Kara Sea. This threshold, whose depth does not exceed 100–150 m, is a barrier that prevents the spreading of the Barents Sea and Arctic waters to the southwestern part of the Kara Sea basin through the St. Anna Trough.     

The connection of geodynamics with seismicity patterns

When a study is to be made of seismic risks, the present-day geodynamic conditions are of fundamental importance: Earthquakes do not happen by themselves, they do have a cause. The cause of earthquakes is that the tectonic stresses exceed a critical limit. The build-up of these stresses is conditioned by the geodynamic processes occurring in the region in question. A knowledge of the geodynamics characteristic of a region is therefore fundamental for seismic risk studies. The general methodology for making such a geodynamic study is based on the entire set of manifestations of the plate-tectonic conditions of that region: these include the mechanism of earthquakes, the stresses observed in mines, the orientation of surface joints and even the direction of river valleys. Examples of geodynamic studies and their bearing on seismic risks are shown from various areas of the world, notably from China, India, and Colombia.     

Silica diagenesis in Eocene shallow-water platform carbonates, southern Pyrenees

Spectacularly developed lower Eocene chert in the Corones platform carbonates of the Spanish Pyrenees is concentrated within a restricted, brackish-water, laminated ostracod-rich facies, which also contains abundant sponge spicules. The chert occurs as nodular, bedded and mottled varieties, and four petrographic types of quartz are developed: microquartz; length-fast (LF) chalcedony; megaquartz; and microspheres. δ¹⁸O values of chert range from 29·6‰ to 30·9‰ (SMOW), which correspond to a broad isotope rank common for biogenic and diagenetic replacement cherts. Calcian dolomite crystals with high Fe and Na are disseminated within the microquartz and LF-chalcedony, but are absent from the megaquartz and host carbonate. The chert is closely associated with desiccation cracks and with interstratal dewatering structures. Load casts are silicified, and laminae rich in sponge spicules are convoluted. Early cracks related to dewatering are filled by microquartz and quartz cements. Ostracod shells within chert are locally fractured; those in the host carbonate are commonly flattened. Late fractures are filled by LF-chalcedony and megaquartz. There is much evidence for the dissolution of sponge spicules and their calcitization in the carbonate host rock. Silica for the Corones cherts was derived from sponges during early diagenesis and shallow burial. Early mechanical compaction and sediment dewatering played a major role in sponge spicule dissolution, migration of silica-rich fluids and the consequent precipitation of chert. Quartz cements continued to be precipitated into the burial environment.     

A Least-squares Window Curves Method for Interpretation of Magnetic Anomalies Caused by Dipping Dikes

We have developed a least-squares method to determine simultaneously the depth and the width of a buried thick dipping dike from residualized magnetic data using filters of successive window lengths. The method involves using a relationship between the depth and the half-width of the source and a combination of windowed observations. The relationship represents a family of curves (window curves). For a fixed window length, the depth is determined for each half-width value by solving one nonlinear equation of the form f (z) = 0 using the least-squares method. The computed depths are plotted against the width values representing a continuous curve. The solution for the depth and the width of the buried dike is read at the common intersection of the window curves. The method involves using a dike model convolved with the same moving average filter as applied to the observed data. As a result, this method can be applied to residuals as well as to measured magnetic data. Procedures are also formulated to estimate the amplitude coefficient and the index parameter. The method is applied to theoretical data with and without random errors. The validity of the method is tested on airborne magnetic data from Canada and on a vertical component magnetic anomaly from Turkey. In all cases examined, the model parameters obtained are in good agreement with the actual ones and with those given in the published literature.