Temporal Variations of Seismic Wave Velocity Associated with 1998 M6.1 Shizukuishi Earthquake

—?We attempt to detect temporal variations of seismic wave velocity before and after 1998 M6.1 Shizukuishi, northeastern Japan, earthquake by using waveform data from explosions and earthquake doublets spanning the period immediately before and after the earthquake. Direct P waves of the second explosion are delayed by ～20 ms at observation stations with epicentral distances less than 15 km. This tendency does not change if the analysis frequency band is changed. Our result suggests that the P-wave velocity decreased by at least 1% in the extremely shallow region of the hanging wall of the M6.1 thrust event after its occurrence. On the other hand, there was the frequency dependence of the coda wave delays for both artificial sources and for natural events. At 5–10 Hz, immediate sharp increases by more than 20 ms in time delays and lower coherency were observed at several stations. We estimated the region in which P-wave velocity might have decreased after the M6.1 earthquake. Maximum depth of the region is 13 km. The region includes the aftershock area of the M6.1 earthquake, but is eccentric to the earthquake in the west. Considering the frequency band analyzed (5–10 Hz), the scale of the spatial inhomogeneity which led to the coda wave delay is several hundreds meters. We investigated candidates for the cause of the direct P-wave and coda wave delay. Observed direct P-wave delay can be partly explained by the stress changes caused by coseismic fault slip. However, the coda wave delay cannot be explained by the stress changes that are limited to the superficial area. Crustal heterogeneity should have changed at coseismic time in the deeper area where aftershocks of the M6.1 earthquake occurred.     

Geochemical and stable isotopic evolution of the Guarani Aquifer System in the state of São Paulo, Brazil

The purpose of this report is to explain geochemical and stable isotopes trends in the Brazilian unit of the Guarani Aquifer System (Botucatu and Piramboia aquifers) in S?o Paulo State, Brazil. Trends of dissolved species concentrations and geochemical modeling indicated a significant role of cation exchange and dissolution of carbonates in downgradient evolution of groundwater chemistry. Loss of calcium by the exchange for sodium drives dissolution of carbonates and results in Na–HCO₃ type of groundwater. The cation-exchange front moves downgradient at probably much slower rate compared to the velocity of groundwater flow and at present is located near to the cities of Sert?ozinho and águas de Santa Barbara (wells PZ-34 and PZ-148, respectively) in a shallow confined area, 50–70 km from the recharge zone. Part of the sodium probably enters the Guarani Aquifer System. together with chloride and sulfate from the underlying Piramboia Formation by diffusion related to the dissolution of evaporates like halite and gypsum. High concentrations of fluorine (up to 13.3 mg/L) can be explained by dissolution of mineral fluoride also driven by cation exchange. However, it is unclear if the dissolution takes place directly in the Guarani Aquifer System or in the overlying basaltic Serra Geral Formation. There is depletion in δ ²H and δ ¹⁸O values in groundwater downgradient. Values of δ ¹³C(DIC) are enriched downgradient, indicating dissolution of calcite under closed system conditions. Values of δ ¹³C(DIC) in deep geothermal wells are very high (>–6.0‰) and probably indicate isotopic exchange with carbonates with δ ¹³C about –3.0‰. Future work should be based on evaluation of vertical fluxes and potential for penetration of contamination to the Guarani Aquifer System. Electronic Publication     

U–Pb zircon geochronology of the North Pole Dome adamellite in the eastern Pilbara Craton

Supracrustal rocks around the North Pole Dome area, Western Australia, provide valuable information regarding early records of the evolution of crustal processes, surface environments, and biosphere. Owing to the occurrence of the oldest known microfossils, the successions at the North Pole Dome area have attracted interest from many researchers. The Paleoarchean successions (Warrawoona Group) mainly comprise mafic‐ultramafic greenstones with intercalated cherts and felsic lavas. Age constraints on the sediments have been mainly based on zircon U–Pb geochronology. However, many zircon grains have suffered from metamictization and contain anomalously high contents of common Pb, which makes interpretation of the U–Pb data complicated. In order to provide more convincing chronological constraints, an U–Pb Concordia age is widely accepted as the best estimate. Most zircons separated from two adamellites also suffered from severe metamictization. In our analyses, less metamictized domains were selected using a pre‐ablation technique in conjunction with elemental mapping, and then their U–Pb isotopic compositions were determined with a laser ablation inductively coupled plasma mass spectrometry. Most analyzed domains contained certain amounts of common Pb (²⁰⁴Pb/²⁰⁶Pb > 0.000 1), whereas three and five U–Pb data points with less common Pb (²⁰⁴Pb/²⁰⁶Pb < 0.000 1) were obtained. These U–Pb datasets yielded U–Pb Concordia ages of ca 3 445 Ma and 3 454 Ma, respectively. These ages represent the timing of the adamellite intrusion, and constrain the minimum depositional age of the Warrawoona Group. In addition, a single xenocrystic zircon grain showed a ²⁰⁷Pb/²⁰⁶Pb age of ca 3 545 Ma, supporting the idea that the sialic basement of the Pilbara Craton existed prior to 3 500 Ma. The in situ U–Pb zircon dating combined with the pre‐ablation technique has the potentials to identify non‐metamictized parts and to yield precise and accurate geochronological data even from partially metamictized zircons.     

Effect of canopy openness and meteorological factors on spatial variability of throughfall isotopic composition in a Japanese cypress plantation

Spatial variability of throughfall (TF) isotopic composition, used as tracer input, influences isotope hydrological applications in forested watersheds. Notwithstanding, identification of the dominant canopy factors and processes that affect the patterns of TF isotopic variability remains ambiguous. Here, we examined the spatio‐temporal variability of TF isotopic composition in a Japanese cypress plantation, in which intensive strip thinning was performed and investigated whether canopy structure at a fine resolution of canopy effect analysis is related to TF isotopic composition and how this is affected by meteorological factors. Canopy openness, as an index of canopy structure, was calculated from hemispherical photographs at different zenith angles. TF samples were collected in a 10 × 10 m experimental plot in both pre‐thinning (from July to November 2010) and post‐thinning (from May 2012 to March 2013) periods. Our results show that thinning resulted in a smaller alteration of input δ¹⁸O of gross precipitation, whereas the changes in deuterium excess varied in both directions. Despite the temporal stability of spatial patterns in TF amount, the spatial variability of TF isotopic composition was not temporally stable in both pre‐ and post‐thinning periods. Additionally, after thinning, the isotopic composition of TF was best related to canopy openness calculated at the zenith angle of 7°, exhibiting three different relationships, that is, significantly negative, significantly positive, and nonsignificant. Changes in meteorological factors (wind speed, rainfall intensity, and temperature) were found to affect the relationships between TF δ¹⁸O and canopy openness. The observed shifts in the relationships reveal different dominant factors (partial evaporation and the selection), and canopy water flowpaths control such differences. This study provides useful insights into the spatial variability of TF isotopic composition and improves our understanding of the physical processes of interception through canopy passage.     

Geochemical characteristics of zircons in the Ashizuri A‐type granitoids: An additional granite topology tool for detrital zircon studies

Zircon is resistant to alteration over a wide range of geological environments, and isotopic ratios within the mineral provide constraints on ages and their parental magmas. Trace element compositions in zircon are also expected to reflect those of their parent magmas, and have a potential as essential indicators for their host rocks. Because most detrital zircons that accumulate at river mouths are derived primarily from granitoids, the classification of zircon within granitoids is potentially meaningful. This study employs the conventional classification scheme of granites (I‐, S‐, M‐, and A‐types). To clarify geochemical characteristics of zircons in A‐type granites, trace element compositions of zircons extracted from the A‐type Ashizuri granitoids were examined. Zircons from the Ashizuri granitoids commonly show enrichments of heavy rare earth elements and positive Ce anomalies, indicating that these zircons were igneous in origin. In addition, zircons in these A‐type granites are characterized by enrichments of Nb, Y, Ta, Th, and U and strong negative Eu anomalies, which exhibit good positive correlations with those in their whole rocks. This fact indicates that these signatures in zircons reflect well those in their parental bodies and are useful in identifying zircons derived from A‐type granite. Based on compilations of available data, zircons from A‐type granites can be clearly discriminated from other‐types of granites within Nb/Sr–Eu anomaly, U/Sr–Eu anomaly, Nb/Sr–U/Sr, and Nb/Sr–Ta/Sr cross‐plots. All indices used in these diagrams were selected based on the geochemical features of both zircon and whole rock of A‐type granites. Application of these discrimination diagrams to detrital zircons will likely provide further insights. For example, some Hadean detrital zircons plot in similar fields to A‐type granites, implying the existence of A‐type magmatism in the Earth's earliest history.     

Detrital zircon multi‐chronology,provenance, and low‐grade metamorphism of the Cretaceous Shimanto accretionary complex,eastern Shikoku,Southwest Japan: Tectonic evolution in response to igneous activity within a subduction zone

Detrital zircon multi‐chronology combined with provenance and low‐grade metamorphism analyses enables the reinterpretation of the tectonic evolution of the Cretaceous Shimanto accretionary complex in Southwest Japan. Detrital zircon U–Pb ages and provenance analysis defines the depositional age of trench‐fill turbidites associated with igneous activity in provenance. Periods of low igneous activity are recorded by youngest single grain zircon U–Pb ages (YSG) that approximate or are older than the depositional ages obtained from radiolarian fossil‐bearing mudstone. Periods of intensive igneous activity recorded by youngest cluster U–Pb ages (YC1σ) that correspond to the younger limits of radiolarian ages. The YC1σ U–Pb ages obtained from sandstones within mélange units provide more accurate younger depositional ages than radiolarian ages derived from mudstone. Determining true depositional ages requires a combination of fossil data, detrital zircon ages, and provenance information. Fission‐track ages using zircons estimated YC1σ U–Pb ages are useful for assessing depositional and annealing ages for the low‐grade metamorphosed accretionary complex. These new dating presented here indicates the following tectonic history of the accretionary wedge. Evolution of the Shimanto accretionary complex from the Albian to the Turonian was caused by the subduction of the Izanagi plate, a process that supplied sediments via the erosion of Permian and Triassic to Early Jurassic granitic rocks and the eruption of minor amounts of Early Cretaceous intermediate volcanic rocks. The complex subsequently underwent intensive igneous activity from the Coniacian to the early Paleocene as a result of the subduction of a hot and young oceanic slab, such as the Kula–Pacific plate. Finally, the major out‐of‐sequence thrusts of the Fukase Fault and the Aki Tectonic Line formed after the middle Eocene, and this reactivation of the Shimanto accretionary complex as a result of the subduction of the Pacific plate.     

Multifractal analysis of spatial distribution of microearthquakes in the Kanto region

b
The spatial distribution of earthquakes is a fractal, which is characterized by a fractal dimension. However, if a spatial distribution has a heterogeneous fractal structure, a single value of fractal dimension [e.g. D_o (capacity dimension) or D ₂ (correlation dimension)] is not enough to characterize it. From a multifractal viewpoint, we analysed the spatial distribution of microearthquakes in the Kanto region by using a local density function. Generalized dimensions, D_q , of the spatial distribution were calculated from the slopes of generalized correlation integrals, C_q(r) versus distance r , on a log-log plot, examining the self-similarity of the spatial distribution of microearthquakes. Self-similar structures are held well at scales from 1.26 to 12.6 km. Our results suggest that the spatial distribution of microearthquakes in the Kanto region is not a homogeneous fractal structure but a heterogeneous one with generalized dimensions D₂ = 2.2 ≤ D ₃≤…≤ D _∞= 1.7. The value of D _∞, the lower limit of fractal dimension, is the fractal dimension of the most intensive clustering in the heterogeneous fractal set. The fractal dimension of the most intensive clustering of microearthquakes in the Kanto region is 1.7.     

Thermal evolution of the Tseel terrane,SW Mongolia and its relation to granitoid intrusions in the Central Asian Orogenic Belt

The timing and thermal effects of granitoid intrusions into accreted sedimentary rocks are important for understanding the growth process of continental crust. In this study, the petrology and geochronology of pelitic gneisses in the Tseel area of the Tseel terrane, SW Mongolia, are examined to understand the relationship between igneous activity and metamorphism during crustal evolution in the Central Asian Orogenic Belt (CAOB). Four mineral zones are recognized on the basis of progressive changes in the mineral assemblages in the pelitic gneisses, namely: the garnet, staurolite, sillimanite and cordierite zones. The gneisses with high metamorphic grades (i.e. sillimanite and cordierite zones) occur in the central part of the Tseel area, where granitoids are abundant. To the north and south of these granitoids, the metamorphic grade shows a gradual decrease. The composition of garnet in the pelitic gneisses varies systematically across the mineral zones, from grossular‐rich garnet in the garnet zone to zoned garnet with grossular‐rich cores and pyrope‐rich rims in the staurolite zone, and pyrope‐rich garnet in the sillimanite and cordierite zones. Thermobarometric analyses of individual garnet crystals reveal two main stages of metamorphism: (i) a high‐P and low‐T stage (as recorded by garnet in the garnet zone and garnet cores in the staurolite zone) at 520–580 °C and 4.5–7 kbar in the kyanite stability field and (ii) a low‐P and high‐T stage (garnet rims in the staurolite zone and garnet in the sillimanite and cordierite zones) at 570–680 °C and 3.0–6.0 kbar in the sillimanite stability field. The earlier high‐P metamorphism resulted in the growth of kyanite in quartz veins within the staurolite and sillimanite zones. The U–Pb zircon ages of pelitic gneisses and granitoids reveal that (i) the protolith (igneous) age of the pelitic gneisses is c. 510 Ma; (ii) the low‐P and high‐T metamorphism occurred at 377 ± 30 Ma; and (iii) this metamorphic stage was coeval with granitoid intrusion at 385 ± 7 Ma. The age of the earlier low‐T and high‐P metamorphism is not clearly recorded in the zircon, but probably corresponds to small age peaks at 450–400 Ma. The low‐P and high‐T metamorphism continued for c. 100 Ma, which is longer than the active period of a single granitoid body. These findings indicate that an elevation of geotherm and a transition from high‐P and low‐T to low‐P and high‐T metamorphism occurred, associated with continuous emplacement of several granitoids, during the crustal evolution in the Devonian CAOB.