Current states of stress in the northern Andes as indicated by focal mechanisms of earthquakes

Systematic inversion of double couple focal mechanisms of shallow earthquakes in the northern Andes reveals relatively homogeneous patterns of crustal stress in three main regions. The first region, presently under the influence of the Caribbean plate, includes the northern segment of the Eastern Cordillera of Colombia and the western flank of the Central Cordillera (north of 4°N). It is characterized by WNW–ESE compression of dominantly reverse type that deflects to NW–SE in the Merida Andes of Venezuela, where it becomes mainly strike–slip in type. A major bend of the Eastern thrust front of the Eastern Cordillera, near its junction with the Merida Andes, coincides with a local deflection of the stress regime (SW–NE compression), suggesting local accommodation of the thrust belt to a rigid indenter in this area. The second region includes the SW Pacific coast of Colombia and Ecuador, currently under the influence of the Nazca plate. In this area, approximately E–W compression is mainly reverse in type. It deflects to WSW–ENE in the northern Andes south of 4°N, where it is accommodated by right-lateral displacement of the Romeral fault complex and the Eastern front of the northern Andes. The third, and most complex, region is the area of the triple junction between the South American, Nazca and Caribbean plates. It reveals two major stress regimes, both mainly strike–slip in type. The first regime involves SW–NE compression related to the interaction between the Nazca and Caribbean plates and the Panama micro-plate, typically accommodated in an E–W left-lateral shear zone. The second regime involves NW–SE compression, mainly related to the interaction between the Caribbean plate and the North Andes block which induces left-lateral displacement on the Uramita and Romeral faults north of 4°N.Deep seismicity (about 150–170 km) concentrates in the Bucaramanga nest and Cauca Valley areas. The inversion reveals a rather homogeneous attitude of the minimum stress axis, which dips towards the E. This extension is consistent with the present plunge of the Nazca and Caribbean slabs, suggesting that a broken slab may be torn under gravitational stresses in the Bucaramanga nest. This model is compatible with current blocking of the subduction in the western northern Andes, inhibiting the eastward displacement of slabs, which are forced to break and sink in to the asthenosphere under their own weight.     

The Okavango giant mafic dyke swarm (NE Botswana): its structural significance within the Karoo Large Igneous Province

The structural organization of a giant mafic dyke swarm, the Okavango complex, in the northern Karoo Large Igneous Province (LIP) of NE Botswana is detailed. This N110°E-oriented dyke swarm extends for 1500 km with a maximum width of 100 km through Archaean basement terranes and Permo-Jurassic sedimentary sequences. The cornerstone of the study is the quantitative analysis of N>170 (exposed) and N>420 (detected by ground magnetics) dykes evidenced on a ca. 80-km-long section lying in crystalline host-rocks, at high-angle to the densest zone of the swarm (Shashe area). Individual dykes are generally sub-vertical and parallel to the entire swarm. Statistical analysis of width data indicates anomalous dyke frequency (few data <5.0 m) and mean dyke thickness (high value of 17 m) with respect to values classically obtained from other giant swarms. Variations of mean dyke thicknesses from 17 (N110°E swarm) to 27 m (adjoining and coeval N70°E giant swarm) are assigned to the conditions hosting fracture networks dilated as either shear or pure extensional structures, respectively, in response to an inferred NNW–SSE extension. Both fracture patterns are regarded as inherited brittle basement fabrics associated with a previous (Proterozoic) dyking event. The Okavango N110°E dyke swarm is thus a polyphase intrusive system in which total dilation caused by Karoo dykes (estimated frequency of 87%) is 12.2% (6315 m of cumulative dyke width) throughout the 52-km-long projected Shashe section. Assuming that Karoo mafic dyke swarms in NE Botswana follow inherited Proterozoic fractures, as similarly applied for most of the nearly synchronous giant dyke complexes converging towards the Nuanetsi area, leads us to consider that the resulting triple junction-like dyke/fracture pattern is not a definitive proof for a deep mantle plume in the Karoo LIP.     

Quantum-thermodynamic treatment of intrinsic anharmonicity; Wallace’s theorem revisited

Wallace (in Thermodynamics of crystals, 1972) developed a theorem, rooted in rigid lattice dynamics, which incorporates intrinsic anharmonic effects in solids. The practical application of this theorem in mineral physics is computationally involved and this is the main reason for the theorem not getting the attention it deserves. Because intrinsic anharmonicity is an important issue at the extreme conditions in planetary mantles, we derived a method which removes the computational obstacles in applying this theorem. We extended the theorem to incorporate details of the phonon spectrum and tested our algorithm on forsterite (Mg₂SiO₄). Using a least squares inversion technique applied to all available experimental data, we show that it results in an accurate representation of thermodynamic properties and sound wave velocities of Mg₂SiO₄ in its complete pressure–temperature stability range. We also show that the accuracy of our results is not significantly affected by the use of a different equation of state.     

Low-pressure adsorption of Ar, Kr, and Xe on carbonaceous materials (kerogen and carbon blacks), ferrihydrite, and montmorillonite: Implications for the trapping of noble gases onto meteoritic matter

Noble gases trapped in meteorites are tightly bound in a carbonaceous carrier labeled “phase Q.” Mechanisms having led to their retention in this phase or in its precursors are poorly understood. To test physical adsorption as a way of retaining noble gases into precursors of meteoritic materials, we have performed adsorption experiments for Ar, Kr, and Xe at low pressures (10⁻⁴ mbar to 500 mbar) encompassing pressures proposed for the evolving solar nebula. Low-pressure adsorption isotherms were obtained for ferrihydrite and montmorillonite, both phases being present in Orgueil (CI), for terrestrial type III kerogen, the best chemical analog of phase Q studied so far, and for carbon blacks, which are present in phase Q and can be considered as possible precursors.Based on adsorption data obtained at low pressures relevant to the protosolar nebula, we propose that the amount of noble gases that can be adsorbed onto primitive materials is much higher than previously inferred from experiments carried out at higher pressures. The adsorption capacity increases from kerogen, carbon blacks, montmorillonite to ferrihydrite. Because of its low specific surface area, kerogen can hardly account for the noble gas inventory of Q. Carbon blacks in the temperature range 75 K-100 K can adsorb up to two orders of magnitude more noble gases than those found in Q. Irreversible trapping of a few percent of noble gases adsorbed on such materials could represent a viable process for incorporating noble gases in phase Q precursors. This temperature range cannot be ruled out for the zone of accretion of the meteorite precursors according to recent astrophysical models and observations, although it is near the lower end of the temperatures proposed for the evolving solar nebula.     

Development of a soil moisture‐based distributed hydrologic model for determining hydrologically based critical source areas

A simple grid cell‐based distributed hydrologic model was developed to provide spatial information on hydrologic components for determining hydrologically based critical source areas. The model represents the critical process (soil moisture variation) to run‐off generation accounting for both local and global water balance. In this way, it simulates both infiltration excess run‐off and saturation excess run‐off. The model was tested by multisite and multivariable evaluation on the 50‐km² Little River Experimental Watershed I in Georgia, U.S. and 2 smaller nested subwatersheds. Water balance, hydrograph, and soil moisture were simulated and compared to observed data. For streamflow calibration, the daily Nash‐Sutcliffe coefficient was 0.78 at the watershed outlet and 0.56 and 0.75 at the 2 nested subwatersheds. For the validation period, the Nash‐Sutcliffe coefficients were 0.79 at the watershed outlet and 0.85 and 0.83 at the 2 subwatersheds. The per cent bias was less than 15% for all sites. For soil moisture, the model also predicted the rising and declining trends at 4 of the 5 measurement sites. The spatial distribution of surface run‐off simulated by the model was mainly controlled by local characteristics (precipitation, soil properties, and land cover) on dry days and by global watershed characteristics (relative position within the watershed and hydrologic connectivity) on wet days when saturation excess run‐off was simulated. The spatial details of run‐off generation and travel time along flow paths provided by the model are helpful for watershed managers to further identify critical source areas of non‐point source pollution and develop best management practices.     

Segmentation and kinematics of the North America‐Caribbean plate boundary offshore Hispaniola

We explored the submarine portions of the Enriquillo–Plantain Garden Fault zone (EPGFZ) and the Septentrional–Oriente Fault zone (SOFZ) along the Northern Caribbean plate boundary using high‐resolution multibeam echo‐sounding and shallow seismic reflection. The bathymetric data shed light on poorly documented or previously unknown submarine fault zones running over 200 km between Haiti and Jamaica (EPGFZ) and 300 km between the Dominican Republic and Cuba (SOFZ). The primary plate‐boundary structures are a series of strike‐slip fault segments associated with pressure ridges, restraining bends, step overs and dogleg offsets indicating very active tectonics. Several distinct segments 50–100 km long cut across pre‐existing structures inherited from former tectonic regimes or bypass recent morphologies formed under the current strike‐slip regime. Along the most recent trace of the SOFZ, we measured a strike‐slip offset of 16.5 km, which indicates steady activity for the past ~1.8 Ma if its current GPS‐derived motion of 9.8 ± 2 mm a^?1 has remained stable during the entire Quaternary.     

Dating megafaunal extinction on the Pleistocene Darling Downs,eastern Australia: the promise and pitfalls of dating as a test of extinction hypotheses

A key to understanding Late Pleistocene megafaunal extinction dynamics is knowledge of megafaunal ecological response(s) to long-term environmental perturbations. Strategically, that requires targeting fossil deposits that accumulated during glacial and interglacial intervals both before and after human arrival, with subsequent palaeoecological models underpinned by robust and reliable chronologies. Late Pleistocene vertebrate fossil localities from the Darling Downs, eastern Australia, provide stratigraphically-intact, abundant megafaunal sequences, which allows for testing of anthropogenic versus climate change megafauna extinction hypotheses. Each stratigraphic unit at site QML796, Kings Creek Catchment, was previously shown to have had similar sampling potential, and the basal units contain both small-sized taxa (e.g., land snails, frogs, bandicoots, rodents) and megafauna. Importantly, sequential faunal horizons show stepwise decrease in taxonomic diversity with the loss of some, but not all, megafauna in the geographically-small palaeocatchment. The purpose of this paper is to present the results of our intensive, multidisciplinary dating study of the deposits (>40 dates). Dating by means of accelerator mass spectrometry (AMS) ¹⁴C (targeting bone, freshwater molluscs, and charcoal) and thermal ionisation mass spectrometry U/Th (targeting teeth and freshwater molluscs) do not agree with each other and, in the case of AMS ¹⁴C dating, lack internal consistency. Scanning electron microscopy and rare earth element analyses demonstrate that the dated molluscs are diagenetically altered and contain aragonite cements that incorporated secondary young C, suggesting that such dates should be regarded as minimum ages. AMS ¹⁴C dated charcoals provide ages that occur out of stratigraphic order, and cluster in the upper chronological limits of the technique (～40–48 ka). Again, we suggest that such results should be regarded as suspicious and only minimum ages. Subsequent OSL and U/Th (teeth) dating provide complimentary results and demonstrate that the faunal sequences actually span ～120–83 ka, thus occurring beyond the AMS ¹⁴C dating window. Importantly, the dates suggest that the local decline in biological diversity was initiated ～75,000 years before the colonisation of humans on the continent. Collectively, the data are most parsimoniously consistent with a pre-human climate change model for local habitat change and megafauna extinction, but not with a nearly simultaneous extinction of megafauna as required by the human-induced blitzkrieg extinction hypothesis. This study demonstrates the problems inherent in dating deposits that lie near the chronological limits of the radiocarbon dating technique, and highlights the need to cross-check previously-dated archaeological and megafauna deposits within the timeframe of earliest human colonisation and latest megafaunal survival.     

Analysis of the interannual variability of annual daily extreme water levels in the St Lawrence River and Lake Ontario from 1918 to 2010

We compared the interannual variability of annual daily maximum and minimum extreme water levels in Lake Ontario and the St Lawrence River (Sorel station) from 1918 to 2010, using several statistical tests. The interannual variability of annual daily maximum extreme water levels in Lake Ontario is characterized by a positive long‐term trend showing two shifts in mean (1929–1930 and 1942–1943) and a single shift in variance (in 1958–1959). In contrast, for the St Lawrence River, this interannual variability is characterized by a negative long‐term trend with a single shift in mean, which occurred in 1955–1956. As for annual daily minimum extreme water levels, their interannual variability shows no significant long‐term change in trend. However, for Lake Ontario, the interannual variability of these water levels shows two shifts in mean, which are synchronous with those for maximum water levels, and a single shift in variance, which occurred in 1965–1966. These changes in trend and stationarity (mean and variance) are thought to be due to factors both climatic (the Great Drought of the 1930s) and human (digging of the Seaway and construction of several dams and locks during the 1950s). Despite this change in means and variance, the four series are clearly described by the generalized extreme value distribution. Finally, annual daily maximum and minimum extreme water levels in the St Lawrence and Lake Ontario are negatively correlated with Atlantic multidecadal oscillation over the period from 1918 to 2010. Copyright © 2013 John Wiley & Sons, Ltd.