Memories of Fire: Eduardo Galeano and the Geography of Guatemala

In Memory of Fire, a poetic narration of the history of the Americas from pre-Columbian times to the late 20th century, Eduardo Galeano furnishes readers with over 1200 of his trademark vignettes, some 35 of which pertain to Guatemala. Galeano evokes disparate aspects of the geography of Guatemala, past and present, in grounded miniatures of time, place, and episode. His sketches of the experiences of Maya peoples allow us to see them as survivors of three cycles of conquest: (1) conquest by imperial Spain; (2) conquest by local and international capitalism; and (3) conquest by state terror. Composed in the literary mode of creative non-fiction, Memory of Fire serves as an inspirational classroom text, exposing students not only to factual detail but also a powerful artistic imagination.     

Fracture analysis in the south-western Corinth rift (Greece) and implications on fault hydraulic behavior

This paper reviews the data concerning the fracture network and the hydraulic characteristics of faults in an active zone of the Gulf of Corinth. Pressure gap measured through fault planes shows that in this area the active normal faults (Aigion, Helike) act, at least temporarily and locally, as transversal seal. The analysis of the carbonate cements in the fractures on both the hangingwall and the footwall of the faults also suggests that they have acted as local seals during the whole fault zone evolution. However, the pressure and the characteristics of the water samples measured in the wells indicate that meteoric water circulates from the highest part of the relief to the coast, which means it goes through the fault zones. Field quantitative analysis and core studies from the AIG-10 well have been performed to define both regional and fault-related fracture networks. Then laboratory thin section observations have been done to recognize the different fault rocks characterizing the fault zone components. These two kinds of approach give information on the permeability characteristics of the fault zone. To synthesize the data, a schematic conceptual 3D fluid flow modeling has been performed taking into account fault zone permeability architecture, sedimentation, fluid flow, fault vertical offset and meteoric water influx, as well as compaction water flow. This modeling allows us to fit all the data with a model where the fault segments act as a seal whereas the relays between these segments allow for the regional flow from the Peloponnese topographic highs to the coast.     

Increasing critical sensitivity of the Load/Unload Response Ratio before large earthquakes with identified stress accumulation pattern

The Load/Unload Response Ratio (LURR) method is proposed for short-to-intermediate-term earthquake prediction [Yin, X.C., Chen, X.Z., Song, Z.P., Yin, C., 1995. A New Approach to Earthquake Prediction — The Load/Unload Response Ratio (LURR) Theory, Pure Appl. Geophys., 145, 701–715]. This method is based on measuring the ratio between Benioff strains released during the time periods of loading and unloading, corresponding to the Coulomb Failure Stress change induced by Earth tides on optimally oriented faults. According to the method, the LURR time series usually climb to an anomalously high peak prior to occurrence of a large earthquake. Previous studies have indicated that the size of critical seismogenic region selected for LURR measurements has great influence on the evaluation of LURR. In this study, we replace the circular region usually adopted in LURR practice with an area within which the tectonic stress change would mostly affect the Coulomb stress on a potential seismogenic fault of a future event. The Coulomb stress change before a hypothetical earthquake is calculated based on a simple back-slip dislocation model of the event. This new algorithm, by combining the LURR method with our choice of identified area with increased Coulomb stress, is devised to improve the sensitivity of LURR to measure criticality of stress accumulation before a large earthquake. Retrospective tests of this algorithm on four large earthquakes occurred in California over the last two decades show remarkable enhancement of the LURR precursory anomalies. For some strong events of lesser magnitudes occurred in the same neighborhoods and during the same time periods, significant anomalies are found if circular areas are used, and are not found if increased Coulomb stress areas are used for LURR data selection. The unique feature of this algorithm may provide stronger constraints on forecasts of the size and location of future large events.     

High-resolution sedimentary record of lead in the coastal mud of the East China Sea in the past one hundred and fifty years

In the East China Sea （ECS）, there are some mud areas, including the south coastal mud area, the north coastal mud area, and the mud area to the southwest of Cheju Island. X-ray fluorescence （XRF） techniques and Thermal Ionization Mass Spectrometry （TIMS） were used to study the high-resolution sedimentary record of Pb concentrations and Pb stable isotopic compositions in the past one hundred and fifty years in the coastal mud of the ECS. Pb concentrations of a ^210Pb dating S5 core in the study area have increased rapidly since 1980, and reached the maximal value with 65.08 μg/g in 2000, corresponding to the fast economic development of China since the implementation of the ＂Reform and Open Policy＂ in 1978; ^206Pb/^207Pb ratios generally had stabilized at 1.195 from 1860 to 1966, and decreased gradually from 1966 to 2000, indicating that the anthropogenic source Pb contribution to the ECS has increased gradually since 1966, especially since 1980. Pb concentrations decreased distinctly from 2000 to 2003 and ^206Pb/^207Pb ratios increased from 2001 to 2003, corresponding closely to the ban of lead gasoline from 2000 in China. From 1950 to 2003, there occurred four distinct decrease events of ^206Pb/^207Pb, possibly responding to the Changjiang River （Yangtze River） catastrophic floods in 1998, 1991, 1981 and 1954; from 1860 to 1966, there were two decrease periods of ^206Pb/^207Pb, which may respond to the catastrophic floods of Changjiang River in 1931 and 1935, and 1870. As a result of the erosion and drowning by the catastrophic floods, the anthropogenic lead accumulated in soil and water environments over a long period of time was brought into the Changjiang River, then part of them was finally transported into the ECS, which leads to changes in Pb stable isotopic composition.     

Mesozoic transtensional basin history of the Eastern Cordillera, Colombian Andes: Inferences from tectonic models

Backstripping analysis and forward modeling of 162 stratigraphic columns and wells of the Eastern Cordillera (EC), Llanos, and Magdalena Valley shows the Mesozoic Colombian Basin is marked by five lithosphere stretching pulses. Three stretching events are suggested during the Triassic–Jurassic, but additional biostratigraphical data are needed to identify them precisely. The spatial distribution of lithosphere stretching values suggests that small, narrow (<150 km), asymmetric graben basins were located on opposite sides of the paleo-Magdalena–La Salina fault system, which probably was active as a master transtensional or strike-slip fault system. Paleomagnetic data suggesting a significant (at least 10°) northward translation of terranes west of the Bucaramanga fault during the Early Jurassic, and the similarity between the early Mesozoic stratigraphy and tectonic setting of the Payandé terrane with the Late Permian transtensional rift of the Eastern Cordillera of Peru and Bolivia indicate that the areas were adjacent in early Mesozoic times. New geochronological, petrological, stratigraphic, and structural research is necessary to test this hypothesis, including additional paleomagnetic investigations to determine the paleolatitudinal position of the Central Cordillera and adjacent tectonic terranes during the Triassic–Jurassic. Two stretching events are suggested for the Cretaceous: Berriasian–Hauterivian (144–127 Ma) and Aptian–Albian (121–102 Ma). During the Early Cretaceous, marine facies accumulated on an extensional basin system. Shallow-marine sedimentation ended at the end of the Cretaceous due to the accretion of oceanic terranes of the Western Cordillera. In Berriasian–Hauterivian subsidence curves, isopach maps and paleomagnetic data imply a (>180 km) wide, asymmetrical, transtensional half-rift basin existed, divided by the Santander Floresta horst or high. The location of small mafic intrusions coincides with areas of thin crust (crustal stretching factors >1.4) and maximum stretching of the subcrustal lithosphere. During the Aptian–early Albian, the basin extended toward the south in the Upper Magdalena Valley. Differences between crustal and subcrustal stretching values suggest some lowermost crustal decoupling between the crust and subcrustal lithosphere or that increased thermal thinning affected the mantle lithosphere. Late Cretaceous subsidence was mainly driven by lithospheric cooling, water loading, and horizontal compressional stresses generated by collision of oceanic terranes in western Colombia. Triassic transtensional basins were narrow and increased in width during the Triassic and Jurassic. Cretaceous transtensional basins were wider than Triassic–Jurassic basins. During the Mesozoic, the strike-slip component gradually decreased at the expense of the increase of the extensional component, as suggested by paleomagnetic data and lithosphere stretching values. During the Berriasian–Hauterivian, the eastern side of the extensional basin may have developed by reactivation of an older Paleozoic rift system associated with the Guaicáramo fault system. The western side probably developed through reactivation of an earlier normal fault system developed during Triassic–Jurassic transtension. Alternatively, the eastern and western margins of the graben may have developed along older strike-slip faults, which were the boundaries of the accretion of terranes west of the Guaicáramo fault during the Late Triassic and Jurassic. The increasing width of the graben system likely was the result of progressive tensional reactivation of preexisting upper crustal weakness zones. Lateral changes in Mesozoic sediment thickness suggest the reverse or thrust faults that now define the eastern and western borders of the EC were originally normal faults with a strike-slip component that inverted during the Cenozoic Andean orogeny. Thus, the Guaicáramo, La Salina, Bitúima, Magdalena, and Boyacá originally were transtensional faults. Their oblique orientation relative to the Mesozoic magmatic arc of the Central Cordillera may be the result of oblique slip extension during the Cretaceous or inherited from the pre-Mesozoic structural grains. However, not all Mesozoic transtensional faults were inverted.     

Carbonates calibrated against oceanographic parameters along a latitudinal transect in the Gulf of California, Mexico

Trophic resources are an important control governing carbonate production. Though this importance has long been recognized, no calibration exists to quantitatively compare biogenic assemblages within trophic resource fields. This study presents a field calibration of carbonate producers in a range of settings against high‐resolution in situ measurements of nutrients, temperature and salinity. With its latitudinal extent from 30° to 23° N, the Gulf of California, Mexico, spans the warm‐temperate realm and encompasses nutrient regimes from oligo‐mesotrophic in the south to eutrophic in the north. Accordingly, from south to north carbonates are characterized by: (i) coral‐dominated shallow carbonate factories (5–20 m water depth) with average sea‐surface temperatures of 25 °C (min. 18 °C, max. 31 °C), average salinities of 35·06‰ and average chlorophyll a levels, which are a proxy for nutrients, of 0·25 mg Chl a m^?3 (max. 0·48, min. 0·1). (ii) Red algal‐dominated subtidal to inner‐shelf carbonate formation (10–25 m) in the central Gulf of California exhibiting average temperatures of 23 °C (min. 18 °C, max. 30 °C), average salinities of 35·25‰, and average Chl a levels of 0·71 Chl a m^?3 (max. 5·62, min. 0). (iii) Molluskan bryozoan‐rich inner to outer shelf factories in the northern Gulf of California (20–50 m) with average sea surface temperatures of only 20 °C (min. 13 °C, max 29 °C), average salinities of 35·01‰, and average contents of 2·2 mg Chl a m^?3 (max. 8·38, min. 0). By calibrating sedimentological data with in situ measured oceanographic information in different environments, the response of carbonate producers to environmental parameters was established and extrapolated to carbonates on a global scale. The results demonstrate the importance of recognizing and quantifying trophic resources as a dominant control determining the biogenic composition and facies character of both modern and fossil carbonates.     

Hydrology and hydrogeology of Sakaryabaşı Karstic springs, Çifteler, Turkey

The Sakarya River is one of the largest rivers in Turkey and is fed mainly from Sakaryabaşı springs. The Sakaryabaşı springs are located in the Central Anatolia and issue from confined/semi-confined karst having a thermal component and therefore, having quite different hydrogeological characteristics as compared to the Taurus Karst region, a typical example of the Mediterranean type of karst. The karstic carbonate rocks that form the groundwater reservoir are overlain by a thick semi-pervious overburden of mainly clastics of Neogene age. Tectonics is the major factor controlling the occurrence of the karst springs in the area where topography is rather flat. This study aimed at explaining the occurrence and movement of the karst groundwater within the system by use of hydrogeological, chemical, and isotopic tools. Isotopic composition of the waters revealed that all waters in the region are of meteoric origin and the thermal component is due to deep circulation. The catchment area of the hydrogeological system extends to the south and groundwater movement is towards the outlets, which are in a depression along a major fault. The movement of the groundwater, based on analysis of remotely sensed images, is controlled mainly by structural elements.     

Correlations in aftershock and seismicity patterns

Correlations in space and time play a fundamental role in earthquake processes. One direct manifestation of the effects of correlations is the occurrence of aftershocks due to the stress transfer in the vicinity of a main shock. Less obvious and more speculative changes in correlations may occur in the background seismicity before large earthquakes. Using statistical physics it is possible to introduce a measure of spatial correlations through a correlation length. This quantity characterizes how local fluctuations can influence the occurrence of earthquakes over distances comparable with the correlation length. In this work, the physical basis of spatial correlations of earthquakes is discussed in the context of critical phenomena and the percolation problem. The method of two-point correlation function is applied to the seismicity of California. Well defined variations in time of the correlation length are found for aftershock sequences and background seismicity. The scaling properties of our obtained distributions are analyzed with respect to changes in several scaling parameters such as lower magnitude cutoff of earthquakes, the maximum time interval between earthquakes, and the spatial size of the area considered. This scaling behavior can be described in a unified manner by utilizing the multifractal fit. Utilizing the percolation approach the time evolution of clusters of earthquakes is studied with the correlation length defined in terms of the radius of gyration of clusters. This method is applied to the seismicity of California.     

A laboratory experiment to monitor the contact state of a fault by transmission waves

We performed a series of laboratory experiments in which elastic waves were transmitted across a simulated fault. Two types of experiments were carried out: (1) Normal Stress Holding Test (NSHT): normal stress was kept constant for about 3 h without shear stress and transmission waves were observed. (2) Shear Stress Increasing Test (SSIT): shear stress was gradually increased until a stick-slip event occurred. Transmission waves were continuously observed throughout the process of stress accumulation. We focused on the change in transmission waves during the application of shear stress and especially during precursory slips.It was found in NSHT that the amplitude of transmission waves linearly increased with the logarithm of stationary contact time. The increase amounted to a few percent after about 3 h. Creep at asperity contacts is responsible for this phenomenon. From a theoretical consideration, it was concluded that the real contact area increased with the logarithm of stationary contact time.We observed in SSIT a significant increase in wave amplitude with shear stress application. This phenomenon cannot be attributed to the time effect observed in NSHT. Instead, it can be explained by the mechanism of “junction growth” proposed by Tabor. Junction growth yields an increase in real contact area. It is required for junction growth to occur that the material in contact is already plastic under a purely normal loading condition. A computer simulation confirmed that this requirement was satisfied in our experiments. We also found that the rate at which the amplitude increased was slightly reduced prior to a stick-slip event. The onset time of the reduction well coincides with the onset of precursory slip. The cause of the reduction is attributed to the reset of stationary contact time due to displacement. This interpretation is supported by the result of NSHT. Taking the time of stationary contact in SSIT into account, we may expect the change in wave amplitude to be, at most, only a few percent. The observed slight reduction in increasing rate is, in this sense, reasonable. The static stiffness of the fault also decreases with precursory slip. It was also found that low frequency waves are a better indicator of precursory slip than high frequency waves. This might suggest that low frequency waves with longer wavelength are a better indicator of average behavior of faults. The problem, however, merits a further investigation. The shifts in phase were also found to be a good indicator of the change in contact state of the fault. The changes in both amplitude and phase of transmission waves are unifyingly understood through the theory of transmission coefficient presented by Pyrak-Nolte et al. Rough surfaces have a tendency to give larger stick-slips than smooth surfaces. The amount of precursory slip is larger for rough surfaces than for smooth surfaces. Although it was confirmed by a computer simulation that rough surfaces have larger contact diameters than smooth surfaces, the rigorous relationship between the surface roughness (contact diameter) and the amount of precursory slips was not established.     

Effects of magnetic interactions in anisotropy of magnetic susceptibility: Models, experiments and implications for igneous rock fabrics quantification

Besides granites of the ilmenite series, in which the anisotropy of magnetic susceptibility (AMS) is mainly controlled by paramagnetic minerals, the AMS of igneous rocks is commonly interpreted as the result of the shape-preferred orientation of unequant ferromagnetic grains. In a few instances, the anisotropy due to the distribution of ferromagnetic grains, irrespective of their shape, has also been proposed as an important AMS source. Former analytical models that consider infinite geometry of identical and uniformly magnetized and coaxial particles confirm that shape fabric may be overcome by dipolar contributions if neighboring grains are close enough to each other to magnetically interact. On these bases we present and experimentally validate a two-grain macroscopic numerical model in which each grain carries its own magnetic anisotropy, volume, orientation and location in space. Compared with analytical predictions and available experiments, our results allow to list and quantify the factors that affect the effects of magnetic interactions. In particular, we discuss the effects of (i) the infinite geometry used in the analytical models, (ii) the intrinsic shape anisotropy of the grains, (iii) the relative orientation in space of the grains, and (iv) the spatial distribution of grains with a particular focus on the inter-grain distance distribution. Using documented case studies, these findings are summarized and discussed in the framework of the generalized total AMS tensor recently introduced by Cañon-Tapia (Cañon-Tapia, E., 2001. Factors affecting the relative importance of shape and distribution anisotropy in rocks: theory and experiments. Tectonophysics, 340, 117–131.). The most important result of our work is that analytical models far overestimate the role of magnetic interaction in rock fabric quantification. Considering natural rocks as an assemblage of interacting and non-interacting grains, and that the effects of interaction are reduced by (i) the finite geometry of the interacting clusters, (ii) the relative orientation between interacting grains, (iii) their heterogeneity in orientation, shape and bulk susceptibility, and (iv) their inter-distance distribution, we reconcile analytical models and experiments with real case studies that minimize the role of magnetic interaction onto the measured AMS. Limitations of our results are discussed and guidelines are provided for the use of AMS in geological interpretation of igneous rock fabrics where magnetic interactions are likely to occur.