Oxygen‐18 dynamics in precipitation and streamflow in a semi‐arid agricultural watershed,Eastern Washington,USA

Understanding flow pathways and mechanisms that generate streamflow is important to understanding agrochemical contamination in surface waters in agricultural watersheds. Two environmental tracers, δ¹⁸O and electrical conductivity (EC), were monitored in tile drainage (draining 12 ha) and stream water (draining nested catchments of 6‐5700 ha) from 2000 to 2008 in the semi‐arid agricultural Missouri Flat Creek (MFC) watershed, near Pullman Washington, USA. Tile drainage and streamflow generated in the watershed were found to have baseline δ¹⁸O value of ?14·7‰ (VSMOW) year round. Winter precipitation accounted for 67% of total annual precipitation and was found to dominate streamflow, tile drainage, and groundwater recharge. ‘Old’ and ‘new’ water partitioning in streamflow were not identifiable using δ¹⁸O, but seasonal shifts of nitrate‐corrected EC suggest that deep soil pathways primarily generated summer streamflow (mean EC 250 µS/cm) while shallow soil pathways dominated streamflow generation during winter (EC declining as low as 100 µS/cm). Using summer isotopic and EC excursions from tile drainage in larger catchment (4700‐5700 ha) stream waters, summer in‐stream evaporation fractions were estimated to be from 20% to 40%, with the greatest evaporation occurring from August to October. Seasonal watershed and environmental tracer dynamics in the MFC watershed appeared to be similar to those at larger watershed scales in the Palouse River basin. A 0·9‰ enrichment, in shallow groundwater drained to streams (tile drainage and soil seepage), of δ¹⁸O values from 2000 to 2008 may be evidence of altered precipitation conditions due to the Pacific Decadal Oscillation (PDO) in the Inland Northwest. Copyright © 2009 John Wiley & Sons, Ltd.     

Weichselian interstadials at Riipiharju,northern Sweden – interpretation of vegetation and climate from fossil and modern pollen records

Hättestrand, M. & Robertsson, A.‐M. 2010: Weichselian interstadials at Riipiharju, northern Sweden – interpretation of vegetation and climate from fossil and modern pollen records. Boreas, 10.1111/j.1502‐3885.2009.00129.x. ISSN 0300‐9483. The most complete records of Weichselian ice‐free conditions in northern Sweden have been retrieved from kettleholes in the Riipiharju esker. In an earlier study, the Riipiharju I core was described as containing two Weichselian interstadials and Riipiharju was chosen as type site for the second Weichselian interstadial in northern Sweden. Here, we present a palynological investigation of two new sediment cores (Riipiharju II and III) retrieved from Riipiharju. Together, the new cores comprise a late cold part of the first Weichselian interstadial recorded in northeastern Sweden (Tärendö I, earlier correlated with Peräpohjola in Finland) as well as a long sequence of the second Weichselian interstadial (Tärendö II, earlier named Tärendö). The results indicate that the climate during deposition of the Tärendö II sequence was more variable than earlier suggested. According to the present interpretation it was relatively warm in the early part of Tärendö II; thereafter a long cold phase persisted, and finally the climate was warmer again in the late part of Tärendö II. The warm phases are characterized by Betula‐dominant pollen assemblages, while the cold phase is characterized by high percentages of Artemisia and Gramineae pollen. Since there is still no firm chronology established of the interstadials in northeastern Sweden, two possible correlations are discussed; either Tärendö I and II are correlated with Brörup (MIS 5c) and Odderade (MIS 5a), or, perhaps more likely, they are correlated with Odderade and early Middle Weichselian (MIS 3) time.     

The data‐driven approach as an operational real‐time flood forecasting model

Accurate water level forecasts are essential for flood warning. This study adopts a data‐driven approach based on the adaptive network–based fuzzy inference system (ANFIS) to forecast the daily water levels of the Lower Mekong River at Pakse, Lao People's Democratic Republic. ANFIS is a hybrid system combining fuzzy inference system and artificial neural networks. Five ANFIS models were developed to provide water level forecasts from 1 to 5 days ahead, respectively. The results show that although ANFIS forecasts of water levels up to three lead days satisfied the benchmark, four‐ and five‐lead‐day forecasts were only slightly better in performance compared with the currently adopted operational model. This limitation is imposed by the auto‐ and cross‐correlations of the water level time series. Output updating procedures based on the autoregressive (AR) and recursive AR (RAR) models were used to enhance ANFIS model outputs. The RAR model performed better than the AR model. In addition, a partial recursive procedure that reduced the number of recursive steps when applying the AR or the RAR model for multi‐step‐ahead error prediction was superior to the fully recursive procedure. The RAR‐based partial recursive updating procedure significantly improved three‐, four‐ and five‐lead‐day forecasts. Our study further shows that for long lead times, ANFIS model errors are dominated by lag time errors. Although the ANFIS model with the RAR‐based partial recursive updating procedure provided the best results, this method was able to reduce the lag time errors significantly for the falling limbs only. Improvements for the rising limbs were modest. Copyright © 2011 John Wiley & Sons, Ltd.     

Influence of heterogeneity on unsaturated hydraulic properties: (1) local heterogeneity and scale effect

Spatial heterogeneity is ubiquitous in nature, which may significantly affect the soil hydraulic property curves. The models of a closed‐form functional relationship of soil hydraulic property curves (e.g. VG model or exponential model) are valid at point or local scale based on a point‐scale hydrological process, but how do scale effects of heterogeneity have an influence on the parameters of these models when the models are used in a larger scale process? This paper uses a two‐dimensional variably saturated flow and solute transport finite element model (VSAFT2) to simulate variations of pressure and moisture content in the soil flume under a constant head boundary condition. By changing different numerical simulation block sizes, a quantitative evaluation of parameter variations in the VG model, resulting from the scale effects, is presented. Results show that the parameters of soil hydraulic properties are independent of scale in homogeneous media. Parameters of α and n in homogeneous media, which are estimated by using the unsaturated hydraulic conductivity curve (UHC) or the soil water retention curve (WRC), are identical. Variations of local heterogeneities strongly affect the soil hydraulic properties, and the scale affects the results of the parameter estimations when numerical experiments are conducted. Furthermore, the discrepancy of each curve becomes considerable when moisture content becomes closer to a dry situation. Parameters estimated by UHC are totally different from the ones estimated by WRC. Copyright © 2012 John Wiley & Sons, Ltd.     

Geoarchaeological Investigations in Mesoamerica Move into the 21st Century: A Review

Over the past thirty years, geoarchaeology has moved from the fringe to mainstream status within Mesoamerican archaeological investigations. This review focuses on works published since the year 2000. Five themes are identified as central to recent studies: (1) the correlation of environmental change and cultural history; (2) anthropogenic environmental impacts; (3) ancient land cover, land use, and diet; (4) archaeological prospection; and (5) provenance studies. These themes are often interwoven in the application of complex systems approaches that allow scientists to more accurately model the intricacies of ancient human–environment interactions.     

Early solar system aqueous activity: K isotope evidence from Allende

The alkali element K is moderately volatile and fluid mobile; thus, it can be influenced by both primary processes (evaporation and recondensation) in the solar nebula and secondary processes (thermal and aqueous alteration) in the parent body. Since these primary and secondary processes would induce different isotopic fractionations, K isotopes could become a potential tracer to distinguish them. Using recently developed methods with improved precision (0.05‰, 95% confidence interval), we systematically measured the K isotopic compositions and major/trace elemental compositions of chondritic components (18 chondrules, 3 CAIs, 2 matrices, and 5 bulks) in the carbonaceous chondrite fall Allende. Among all the components analyzed in this study, CAIs, which formed initially under high‐temperature conditions in the solar nebula and were dominated by nominally K‐free refractory minerals, have the highest K₂O content (average 0.53 wt%) and have K isotope compositions most enriched in heavy isotopes (δ⁴¹K: ?0.30 to ?0.25‰). Such an observation is consistent with previous petrologic studies that show CAIs in Allende have undergone alkali enrichment during metasomatism. In contrast, chondrules contain lower K₂O content (0.003–0.17 wt%) and generally lighter K isotope compositions (δ⁴¹K: ?0.87‰ to ?0.24‰). The matrix and bulks are nearly identical in K₂O content and K isotope compositions (0.02–0.05 wt%; δ⁴¹K: ?0.62 to ? 0.46‰), which are, as expected, right in the middle of CAIs and chondrules. This strongly indicates that most of the chondritic components of Allende suffered aqueous alteration and their K isotopic compositions are the ramification of Allende parent‐body processing instead of primary nebular signatures. Nevertheless, we propose the small K isotope fractionations observed (< 1‰) among Allende components are likely similar to the overall range of K isotopic fractionation that occurred in nebular environment. Furthermore, the K isotope compositions seen in the components of Allende in this study are consistent with MC‐ICP‐MS analyses of the components in ordinary chondrites, which also show an absence of large (10‰) isotope fractionations. This is not expected as evaporation experiments in nebular conditions suggest there should be large K isotopic fractionations. Nevertheless, possible nebular processes such as chondrules back exchanging with ambient gas when they formed could explain this lack of large K isotopic variation.     

Unique achondrite Northwest Africa 11042: Exploring the melting and breakup of the L chondrite parent body

Northwest Africa (NWA) 11042 is a heavily shocked achondrite with medium‐grained cumulate textures. Its olivine and pyroxene compositions, oxygen isotopic composition, and chromium isotopic composition are consistent with L chondrites. Sm‐Nd dating of its primary phases shows a crystallization age of 4100 ± 160 Ma. Ar‐Ar dating of its shocked mineral maskelynite reveals an age of 484.0 ± 1.5 Ma. This age coincides roughly with the breakup event of the L chondrite parent body evident in the shock ages of many L chondrites and the terrestrial record of fossil L chondritic chromite. NWA 11042 shows large depletions in siderophile elements (<0.01×CI) suggestive of a complex igneous history involving extraction of a Fe‐Ni‐S liquid on the L chondrite parent body. Due to its relatively young crystallization age, the heat source for such an igneous process is most likely impact. Because its mineralogy, petrology, and O isotopes are similar to the ungrouped achondrite NWA 4284 (this work), the two meteorites are likely paired and derived from the same parent body.     

Strain analysis of a seismically imaged mass‐transport complex,offshore Uruguay

Strain style, magnitude and distribution within mass‐transport complexes (MTCs) are important for understanding the process evolution of submarine mass flows and for estimating their runout distances. Structural restoration and quantification of strain in gravitationally driven passive margins have been shown to approximately balance between updip extensional and downdip contractional domains; such an exercise has not yet been attempted for MTCs. We here interpret and structurally restore a shallowly buried (c. 1,500 mbsf) and well‐imaged MTC, offshore Uruguay using a high‐resolution (12.5 m vertical and 15 × 12.5 m horizontal resolution) three‐dimensional seismic‐reflection survey. This allows us to characterise and quantify vertical and lateral strain distribution within the deposit. Detailed seismic mapping and attribute analysis shows that the MTC is characterised by a complicated array of kinematic indicators, which vary spatially in style and concentration. Seismic‐attribute extractions reveal several previously undocumented fabrics preserved in the MTC, including internal shearing in the form of sub‐orthogonal shear zones, and fold‐thrust systems within the basal shear zone beneath rafted‐blocks. These features suggest multiple transport directions and phases of flow during emplacement. The MTC is characterised by a broadly tripartite strain distribution, with extensional (e.g. normal faults), translational and contractional (e.g. folds and thrusts) domains, along with a radial frontally emergent zone. We also show how strain is preferentially concentrated around intra‐MTC rafted‐blocks due to their kinematic interactions with the underlying basal shear zone. Overall, and even when volume loss within the frontally emergent zone is included, a strain difference between extension (1.6–1.9 km) and contraction (6.7–7.3 km) is calculated. We attribute this to a combination of distributed, sub‐seismic, ‘cryptic’ strain, likely related to de‐watering, grain‐scale deformation and related changes in bulk sediment volume. This work has implications for assessing MTCs strain distribution and provides a practical approach for evaluating structural interpretations within such deposits.     

Northwest Africa 11024—A heated and dehydrated unique carbonaceous (CM) chondrite

Based on the high abundance of fine‐grained material and its dark appearance, NWA 11024 was recognized as a CM chondrite, which is also confirmed by oxygen isotope measurements. But contrary to known CM chondrites, the typical phases indicating aqueous alteration (e.g., phyllosilicates, carbonates) are missing. Using multiple analytical techniques, this study reveals the differences and similarities to known CM chondrites and will discuss the possibility that NWA 11024 is the first type 3 CM chondrite. During the investigation, two texturally apparent tochilinite–cronstedtite intergrowths were identified within two thin sections. However, the former phyllosilicates were recrystallized to Fe‐rich olivine during a heating event without changing the textural appearance. A peak temperature of 400–600 °C is estimated, which is not high enough to destroy or recrystallize calcite grains. Thus, calcites were never constituents of the mineral paragenesis. Another remarkable feature of NWA 11024 is the occurrence of unknown clot‐like inclusions (UCLIs) within fine‐grained rims, which are unique in this clarity. Their density and S concentration are significantly higher than of the surrounding fine‐grained rim and UCLIs can be seen as primary objects that were not formed by secondary alteration processes inside the rims. Similarities to chondritic and cometary interplanetary dust particles suggest an ice‐rich first‐generation planetesimal for their origin. In the earliest evolution, NWA 11024 experienced the lowest degree of aqueous alteration of all known CM chondrites and subsequently, a heating event dehydrated the sample. We suggest to classify the meteorite NWA 11024 as the first type 3 CM chondrite similar to the classification of CV3 chondrites (like Allende) that could also have lost their matrix phyllosilicates by thermal dehydration.     

Geological and geophysical studies of the Agoudal impact structure (Central High Atlas,Morocco): New evidence for crater size and age

Since the discovery of shatter cones (SCs) near the village of Agoudal (Morocco, Central High Atlas Mountains) in 2013, the absence of one or several associated circular structures led to speculation about the age of the impact event, the number, and the size of the impact crater or craters. Additional constraints on the crater size, age, and erosion rates are obtained here from geological, structural, and geophysical mapping and from cosmogenic nuclide data. Our geological maps of the Agoudal impact site at the scales of 1:30,000 (6 km²) and 1:15,000 (2.25 km²) include all known occurrences of SCs in target rocks, breccias, and vertical to overturned strata. Considering that strata surrounding the impact site are subhorizontal, we argue that disturbed strata are related to the impact event. Three types of breccias have been observed. Two of them (br1‐2 and br2) could be produced by erosion–sedimentation–consolidation processes, with no evidence for impact breccias, while breccia (br1) might be impact related. The most probable center of the structure is estimated at 31°59′13.73?N, 5°30′55.14?W using the concentric deviation method applied to the orientation of strata over the disturbed area. Despite the absence of a morphological expression, the ground magnetic and electromagnetic surveys reveal anomalies spatially associated with disturbed strata and SC occurrences. The geophysical data, the structural observations, and the area of occurrence of SCs in target rocks are all consistent with an original size of 1.4–4.2 km in diameter. Cosmogenic nuclide data (³⁶Cl) constrain the local erosion rates between 220 ± 22 m Ma^?1 and 430 ± 43 m Ma^?1. These erosion rates may remove the topographic expression of such a crater and its ejecta in a time period of about 0.3–1.9 Ma. This age is older than the Agoudal iron meteorite age (105 ± 40 kyr). This new age constraint excludes the possibility of a genetic relationship between the Agoudal iron meteorite fall and the formation of the Agoudal impact site. A chronolgy chart including the Atlas orogeny, the alternation of sedimentation and erosion periods, and the meteoritic impacts is presented based on all obtained and combined data.