Neoproterozoic alkaline magmatism at Sivamalai, southern India

The Sivamalai alkaline complex lies at the southern margin of the Cauvery Shear System that separates the Archaean and Proterozoic domains of the Southern Granulite Terrain in India. U–Pb TIMS dating of zircon from a pegmatitic syenite sample in the complex yields a concordant age of 590.2 ± 1.3 (2σ) Ma which is interpreted to date the intrusion of the alkaline rocks. A lower concordia intercept at 168 ± 210 Ma defined by two grains with high common lead may indicate post-magmatic disturbances due to recrystallisation which is also evident in the CL images of the zircons. EPMA dating of monazite from a post-kinematic pegmatite which intrudes the crystalline basement hosting the alkaline rocks yields an age of 478 ± 29 (2σ) Ma and provides a lower bracket for the main phase of tectonism in this part of the Southern Granulite Terrain. The Pan-African high-grade metamorphism and ductile deformation has thus most likely affected the alkaline rocks. This is supported by the presence of a metamorphic foliation and extensive recrystallisation textures seen in the rocks. The major and trace element concentrations measured on selected samples reveals the presence of both enriched and depleted rock types. The enriched group includes ferrosyenite and nepheline syenite while the depleted group has only nepheline syenites. The trace element depletion of some nepheline syenites is interpreted to be a result of fractional crystallization involving the removal of accessory phases like zircon, titanite, apatite and allanite.     

Aragonite cycles: diagenesis caught in the act

Organic biomarker and nannofossil assemblages are used in combination with sedimentary petrology to identify the role of diagenesis for the formation of a rhythmic calcareous succession. A record from the Early Pliocene of the Maldives (Ocean Drilling Program, ODP Leg 115) is presented containing precession and eccentricity cycles expressed as variable aragonite content. Additional sub‐Milankovitch cycles are caused by rhythmic precipitation of calcite cement in the lower part of the interval. Comparison with palaeo‐productivity indicators (nannofossils, chlorin, total organic carbon) suggests that cementation occurs preferentially in intervals characterized by increasing or decreasing productivity. The coupled variability in productivity and carbonate diagenesis is attributed to the effect of organic matter degradation in the sediment. The observed combination of primary and diagenetic factors hampers the interpretation of the cyclicity on Milankovitch and sub‐Milankovitch scales. Diagenetically stable proxies for palaeo‐productivity, such as nannofossil assemblage data, were used to distinguish between palaeoclimate and diagenetic influences.     

Three-dimensional model and late stage warping of the Plattengneis Shear Zone in the Eastern Alps

The Plattengneis shear zone is a 250–600 m thick, flat lying, Cretaceous, eclogite facies, mylonitic shear zone, with north-over-south transport direction, that is exposed over almost 1000 km² in the Koralpe region along the eastern margin of the Alps. Although the shear zone is one of the largest in the Alps, its role in the Eoalpine metamorphic evolution and the subsequent exhumation of the region, remain enigmatic and its large-scale geometry is not well understood. The outcrop pattern suggests that the shear zone is made up of a single sheet that is folded into a series of open syn- and antiforms with wavelengths of about 10 km. Eclogite bodies occur above, within and below the shear zone and there is no metamorphic grade change across the shear zone. In the south, the fold axes strike east–west and plunge shallowly to the east. In the north, the fold axes are oriented in north–south direction and form a dome shaped structure of the shear zone. Total shortening during this late stage warping event was of the order of 5%. Indirect evidence constrains this folding event to have occurred between 80 and 50 Ma and the fold geometry implies that the final exhumation in the Koralpe occurred somewhat later than further north. Interestingly, the shear zone appears to strike out of the topography in the south and dip into the topography in the north, so that north of the shear zone only hanging-wall rocks are exposed and south of it only foot-wall rocks. Possibilities for the geometric relationship of the Plattengneis shear zone with the surrounding south dipping detachments are discussed.     

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.     

Non-volcanic deep low-frequency tremors accompanying slow slips in the southwest Japan subduction zone

Non-volcanic deep low-frequency tremors in southwest Japan exhibit a strong temporal and spatial correlation with slow slip detected by the dense seismic network. The tremor signal is characterized by a low-frequency vibration with a predominant frequency of 0.5–5 Hz without distinct P- or S-wave onset. The tremors are located using the coherent pattern of envelopes over many stations, and are estimated to occur near the transition zone on the plate boundary on the forearc side along the strike of the descending Philippine Sea plate. The belt-like distribution of tremors consists of many clusters. In western Shikoku, the major tremor activity has a recurrence interval of approximately six months, with each episode lasting over a week. The tremor source area migrates during each episode along the strike of the subducting plate with a migration velocity of about 10 km/day. Slow slip events occur contemporaneously with this tremor activity, with a coincident estimated source area that also migrates during each episode. The coupling of tremor and slow slip in western Shikoku is very similar to the episodic tremor and slip phenomenon reported for the Cascadia margin in northwest North America. The duration and recurrence interval of these episodes varies between tremor clusters even on the same subduction zone, attributable to regional difference in the frictional properties of the plate interface.     

Deep structure of the Nojima Fault, southwest Japan, estimated from borehole observations of fault-zone trapped waves

To estimate the deep structure of the southern part of the Nojima Fault, southwest Japan without the influence of near-surface structures, we analyzed the Love-wave-type fault-zone trapped waves (LTWs) recorded by a borehole seismometer at 1800 m depth. We examined the polarization, dispersion, and dominant frequency of the wavetrain following the direct S-wave in each seismogram to identify the LTW. We selected eight candidates for typical LTWs from 462 records. Because the duration of the LTW increases with hypocentral distance, we infer that the low velocity fault-zone of the Nojima Fault continues towards the seismogenic depth. In addition, since the duration of the LTW increases nonlinearly with hypocentral distance, we infer that the S-wave velocity of the fault-zone increases with depth. The location of events showing the LTW indicates that the fault-zone dips to the southeast at 75° and continues to a depth of approximately 10 km. We assumed a uniform low-velocity waveguide to estimate the average structure of the fault-zone. We estimated the average width, S-wave velocity, and Q_s of the fault-zone by comparing an analytical solution of the LTW with measured data. The average width, S-wave velocity, and Q_s of the fault-zone are 150 to 290 m, 2.5 to 3.2 km/s, and 40 to 90, respectively. Hence the fault-zone structure with a larger width and smaller velocity reduction than the fault-zone model estimated by previous surface observation is more suitable to represent the average fault-zone structure of the Nojima fault. The present study also indicated that the shallow layers and/or a shallow fault-zone structure drastically changes the characteristics of the LTW recorded at the surface, and therefore cause a discrepancy in the fault-zone model between the borehole observation and surface observation.     

Geophysical methods applied to fault characterization and earthquake potential assessment in the Lower Tagus Valley, Portugal

The study region is located in the Lower Tagus Valley, central Portugal, and includes a large portion of the densely populated area of Lisbon. It is characterized by a moderate seismicity with a diffuse pattern, with historical earthquakes causing many casualties, serious damage and economic losses. Occurrence of earthquakes in the area indicates the presence of seismogenic structures at depth that are deficiently known due to a thick Cenozoic sedimentary cover. The hidden character of many of the faults in the Lower Tagus Valley requires the use of indirect methodologies for their study. This paper focuses on the application of high-resolution seismic reflection method for the detection of near-surface faulting on two major tectonic structures that are hidden under the recent alluvial cover of the Tagus Valley, and that have been recognized on deep oil-industry seismic reflection profiles and/or inferred from the surface geology. These are a WNW–ESE-trending fault zone located within the Lower Tagus Cenozoic basin, across the Tagus River estuary (Porto Alto fault), and a NNE–SSW-trending reverse fault zone that borders the Cenozoic Basin at the W (Vila Franca de Xira–Lisbon fault). Vertical electrical soundings were also acquired over the seismic profiles and the refraction interpretation of the reflection data was carried out. According to the interpretation of the collected data, a complex fault pattern disrupts the near surface (first 400 m) at Porto Alto, affecting the Upper Neogene and (at least for one fault) the Quaternary, with a normal offset component. The consistency with the previous oil-industry profiles interpretation supports the location and geometry of this fault zone. Concerning the second structure, two major faults were detected north of Vila Franca de Xira, supporting the extension of the Vila Franca de Xira–Lisbon fault zone northwards. One of these faults presents a reverse geometry apparently displacing Holocene alluvium. Vertical offsets of the Holocene sediments detected in the studied geophysical data of Porto Alto and Vila Franca de Xira–Lisbon faults imply minimum slip rates of 0.15–0.30 mm/year, three times larger than previously inferred for active faults in the Lower Tagus Valley and maximum estimates of average return periods of 2000–5000 years for M 6.5–7 co-seismic ruptures.     

A way to synchronize models with seismic faults for earthquake forecasting: Insights from a simple stochastic model

Numerical models are starting to be used for determining the future behaviour of seismic faults and fault networks. Their final goal would be to forecast future large earthquakes. In order to use them for this task, it is necessary to synchronize each model with the current status of the actual fault or fault network it simulates (just as, for example, meteorologists synchronize their models with the atmosphere by incorporating current atmospheric data in them). However, lithospheric dynamics is largely unobservable: important parameters cannot (or can rarely) be measured in Nature. Earthquakes, though, provide indirect but measurable clues of the stress and strain status in the lithosphere, which should be helpful for the synchronization of the models.The rupture area is one of the measurable parameters of earthquakes. Here we explore how it can be used to at least synchronize fault models between themselves and forecast synthetic earthquakes. Our purpose here is to forecast synthetic earthquakes in a simple but stochastic (random) fault model. By imposing the rupture area of the synthetic earthquakes of this model on other models, the latter become partially synchronized with the first one. We use these partially synchronized models to successfully forecast most of the largest earthquakes generated by the first model. This forecasting strategy outperforms others that only take into account the earthquake series. Our results suggest that probably a good way to synchronize more detailed models with real faults is to force them to reproduce the sequence of previous earthquake ruptures on the faults. This hypothesis could be tested in the future with more detailed models and actual seismic data.     

Great earthquakes of variable magnitude at the Cascadia subduction zone

Comparison of histories of great earthquakes and accompanying tsunamis at eight coastal sites suggests plate-boundary ruptures of varying length, implying great earthquakes of variable magnitude at the Cascadia subduction zone. Inference of rupture length relies on degree of overlap on radiocarbon age ranges for earthquakes and tsunamis, and relative amounts of coseismic subsidence and heights of tsunamis. Written records of a tsunami in Japan provide the most conclusive evidence for rupture of much of the plate boundary during the earthquake of 26 January 1700. Cascadia stratigraphic evidence dating from about 1600 cal yr B.P., similar to that for the 1700 earthquake, implies a similarly long rupture with substantial subsidence and a high tsunami. Correlations are consistent with other long ruptures about 1350 cal yr B.P., 2500 cal yr B.P., 3400 cal yr B.P., 3800 cal yr B.P., 4400 cal yr B.P., and 4900 cal yr B.P. A rupture about 700-1100 cal yr B.P. was limited to the northern and central parts of the subduction zone, and a northern rupture about 2900 cal yr B.P. may have been similarly limited. Times of probable short ruptures in southern Cascadia include about 1100 cal yr B.P., 1700 cal yr B.P., 3200 cal yr B.P., 4200 cal yr B.P., 4600 cal yr B.P., and 4700 cal yr B.P. Rupture patterns suggest that the plate boundary in northern Cascadia usually breaks in long ruptures during the greatest earthquakes. Ruptures in southernmost Cascadia vary in length and recurrence intervals more than ruptures in northern Cascadia.     

三江平原旱田耕作对湿地土壤理化性质的累积影响初探

在中国科学院三江平原湿地生态试验站综合实验场,对不同耕作年限的湿地土壤(0~20 cm)进行环境累积效应分析。结果表明,随着开垦时间的增加,土壤的理化性质发生渐变,物理性质方面,土壤容重和比重逐渐增大,而孔隙度和田间持水量逐年减少;化学性质方面,土壤pH值随开垦时间的增加而增加,有机质和其他养分则随开垦时间的增加而逐年降低。弃耕后土壤性质有所恢复。土壤性质在开垦初期变化较明显,而后逐渐变缓。