Research on crustal structure and lithosphere property in the Okinawa Trough area

A mass of geological, geophysical and geochemical data and information from the Okinawa Trough area are collected for comprehensive research in the study area from East China to Okinawa Trough and then to Ryukyu Island Are region. According to the seismic tomography result (P and S wave) and the processing result of free-air and Bouguer gravity anomaly and magnetic anomaly data in the study area, the comprehensive interpretation is carried out. The Moho depth distribution of the study area is obtained by the inversion calculation based on gravity data using the Harmonious Series method. The crust properties are analyzed. Meantime, some Cenozoic basalt data from Kuandian (NE China), Hannuoba (North China), Minxi (South China), Penghu Islands (Taiwan Strait), Okinawa Trough and Japan Island Arc regions are chosen to make the comparison research on element- isotopes. The result indicates that the lithosphere thickness in the Okinawa Trough area has obviously decreased, where a Low -velocity layer of upper-mantle has reached the Moho interface and the metasometized asthenosphere has formed. The research result on element- isotopes shows that the characteristic of the crust in the Okinawa Trough area is different from that in East China area and the Ryukyu Island Arc area. It is considered that the crust in the Okinawa Trough area belongs to the transition type, which is quite similar to the feature of the oceanic crust.

Analysis on deep crustal structure along the onshore-offshore seismic profile across the Binhai (littoral) fault zone in northeastern South China Sea

The onshore-offshore deep seismic experiment was carried out for the first time and filled the blankness of the seismic surveys in the transition area between South China and northeastern South China Sea. The seismic data were analyzed and processed. The different seismic phases were identified and their travel time arrivals were modeled by ray-tracing to study the P-wave velocity crustal structure of this area. The crustal structure of this area is the continental crust. The crust thickness is gradually decreasing southward along the on-shore-offshore seismic line. The low-velocity layer (5.5 similar to 5.9 km (.) s(-1)) exists generally in the middle crust (about 10.0 similar to 18.0km)with about 2.5 similar to 4.0 km thickness, which is also thinning seaward. No obvious high-velocity layer appears in the lower crust. The Binhai (littoral) fault zone is a low velocity zone, which is located about 35km southeast to the Nan'ao station and corresponding to the gradient belt of gravity & magnetism anomalies. The depth of the fault zone is close to the Moho discontinuity. The littoral fault zone is a boundary between the normal continental crust of South China and the thinned continental crust of the sea area.

Damage Localization, Sensitivity of Energy Release and the Catastrophe Transition

Large earthquakes can be viewed as catastrophic ruptures in the earth’s crust. There are two common features prior to the catastrophe transition in heterogeneous media. One is damage localization and the other is critical sensitivity; both of which are related to a cascade of damage coalescence. In this paper, in an attempt to reveal the physics underlying the catastrophe transition, analytic analysis based on mean-field approximation of a heterogeneous medium as well as numerical simulations using a network model are presented. Both the emergence of damage localization and the sensitivity of energy release are examined to explore the inherent statistical precursors prior to the eventual catastrophic rupture. Emergence of damage localization, as predicted by the mean-field analysis, is consistent with observations of the evolution of damage patterns. It is confirmed that precursors can be extracted from the time-series of energy release according to its sensitivity to increasing crustal stress. As a major result, present research indicates that the catastrophe transition and the critical point hypothesis (CPH) of earthquakes are interrelated. The results suggest there may be two cross-checking precursors of large earthquakes: damage localization and critical sensitivity.

Long-range automaton models of earthquakes: Power-law accelerations, correlation evolution, and mode-switching

We introduce a conceptual model for the in-plane physics of an earthquake fault. The model employs cellular automaton techniques to simulate tectonic loading, earthquake rupture, and strain redistribution. The impact of a hypothetical crustal elastodynamic Green's function is approximated by a long-range strain redistribution law with a r(-p) dependance. We investigate the influence of the effective elastodynamic interaction range upon the dynamical behaviour of the model by conducting experiments with different values of the exponent (p). The results indicate that this model has two distinct, stable modes of behaviour. The first mode produces a characteristic earthquake distribution with moderate to large events preceeded by an interval of time in which the rate of energy release accelerates. A correlation function analysis reveals that accelerating sequences are associated with a systematic, global evolution of strain energy correlations within the system. The second stable mode produces Gutenberg-Richter statistics, with near-linear energy release and no significant global correlation evolution. A model with effectively short-range interactions preferentially displays Gutenberg-Richter behaviour. However, models with long-range interactions appear to switch between the characteristic and GR modes. As the range of elastodynamic interactions is increased, characteristic behaviour begins to dominate GR behaviour. These models demonstrate that evolution of strain energy correlations may occur within systems with a fixed elastodynamic interaction range. Supposing that similar mode-switching dynamical behaviour occurs within earthquake faults then intermediate-term forecasting of large earthquakes may be feasible for some earthquakes but not for others, in alignment with certain empirical seismological observations. Further numerical investigation of dynamical models of this type may lead to advances in earthquake forecasting research and theoretical seismology.

Tectonic Evolution and Dynamics of Deepwater Area of Pearl River Mouth Basin, Northern South China Sea

Quantitative studies on the evolution and dynamics of the deepwater area of Pearl River Mouth basin (PRMB) were carried out based on the latest geological and seismic data. The study area is generally in an extensional state during the Cenozoic. The major extension happened in the earlier syn-rift stages before 23 Ma and the extension after 23 Ma is negligible. Two rapid subsidence periods, 32-23 Ma and 5.3-2.6 Ma, are identified, which are related to the abrupt heat decay during margin breakup and the collision between the Philippine Sea plate and the Eurasian plate, respectively. The strongest crustal thinning in the Baiyun (sic) sag may trigger the syn-rift volcanism along the weak faulted belt around the sag. The Cenozoic tectonic evolution of the study area could be divided into five stages: rifting (similar to 50-40 Ma), rift-drift transition (similar to 40-32 Ma), early post-breakup (similar to 32-23 Ma), thermal subsidence (similar to 23-5.3 Ma) and neotectonic movement (similar to 5.3-0 Ma).

Geochemical properties of the variolitic andesite from Susong in the Dabie orogenic belt, China

The variolitic andesite from the Susong County in the Dabie Mountains implies that it was erupted in water. The mineralogy of the varioles is primarily radiate plagioclase (albite sind oligoclase), with little pyroxene, hornblende and quartz (derived from alteration). The pyroxene, hornblende and quartz are in the interstices between plagiocalse. The matrix consists of glass, hornblende, chlorite, epidote and zoisite. It is clearly subjected an extensive alteration. The andesite has an uncommon chemical composition. The SiO2 content is about 56.8%, TiO2 = 0.9%, MgO = 6.4%, Fe2O3 (tot) = 6.7%similar to 7.6%, 100Mg/(Mg + Fe) = 64.1 similar to 66.2. Mg-# is significantly high. The andesite has high abundances of large-lithophile trace elements (e.g. K, Ba. Sr, LREE), e.g. La/Nb = 5.56 similar to 6.07, low abundances of high-strength-field elements (HFSE e.g. Ta, Nb, P, Ti), particularly Ta and Nb strongly depleted. These are consistent with the characteristics of subduction-related magmas. In the spider diagram of trace elements, from Ce to right hand, the abundances of elements decrease quickly, showing a character of the continental margins. There has a strong punishment of light-rare-earth elements, with a significant diffraction of REEs (the mean value of (La/Yb)(N) is 32.84). No Eu anomaly, but there are anomaly high (La/Yb)(N) = 28.63 similar to 36.74, (La/Y)(N) = 70.33 similar to 82.4. The elements Y and Yb are depleted greatly, Y<20 g/g, Y-N = 2.74 similar to 2.84, Yb-N = 2.18 similar to 2.35. From the La-(La/Sm) diagram, the andesite is derived from partial melting. But the epsilone value of Nd is - 18.7 similar to -19.2, so that the material source may be the mantle materials affected by the crustal materials. The Nd model age is 1.9 Ga indicating that the variolitic basaltic andesite was resulted from the mantle wedge of North China block, which had the Nd model age of 2.5Ga, when the Yangze block which had the Nd model age of 1.7Ga subducted beneath it. So the variolitic andesite has characteristics of the island-are volconic rocks oil a continental basement in the vicinity of the destructive continental margin.

Genesis of Th-230 excess in basalts from mid-ocean ridges and ocean islands: Constraints from the global U-series isotope database and major and rare earth element geochemistry

Based on Th-230-U-238 disequilibrium and major element data from mid-ocean ridge basalts (MORBs) and ocean island basalts (OIBs), this study calculates mantle melting parameters, and thereby investigates the origin of Th-230 excess. (Th-230/U-238) in global MORBs shows a positive correlation with Fe-8, P (o), Na-8, and F-melt (Fe-8 and Na-8 are FeO and Na2O contents respectively after correction for crustal fractionation relative to MgO = 8 wt%, P (o)=pressure of initial melting and F (melt)=degree of melt), while Th-230 excess in OIBs has no obvious correlation with either initial mantle melting depth or the average degree of mantle melting. Furthermore, compared with the MORBs, higher (Th-230/U-238) in OIBs actually corresponds to a lower melting degree. This suggests that the Th-230 excess in MORBs is controlled by mantle melting conditions, while the Th-230 excess in OIBs is more likely related to the deep garnet control. The vast majority of calculated initial melting pressures of MORBs with excess Th-230 are between 1.0 and 2.5 GPa, which is consistent with the conclusion from experiments in recent years that D (U)> D (Th) for Al-clinopyroxene at pressures of > 1.0 GPa. The initial melting pressure of OIBs is 2.2-3.5 GPa (around the spinel-garnet transition zone), with their low excess Ra-226 compared to MORBs also suggesting a deeper mantle source. Accordingly, excess Th-230 in MORBs and OIBs may be formed respectively in the spinel and garnet stability field. In addition, there is no obvious correlation of K2O/TiO2 with (Th-230/U-238) and initial melting pressure (P (o)) of MORBs, so it is proposed that the melting depth producing excess Th-230 does not tap the spinel-garnet transition zone. OIBs and MORBs in both (Th-230/U-238) vs. K2O/TiO2 and (Th-230/U-238) vs. P (o) plots fall in two distinct areas, indicating that the mineral phases which dominate their excess Th-230 are different. Ce/Yb-Ce curves of fast and slow ridge MORBs are similar, while, in comparison, the Ce/Yb-Ce curve for OIBs shows more influence from garnet. The mechanisms generating excess Th-230 in MORBs and OIBs are significantly different, with formation of excess Th-230 in the garnet zone only being suitable for OIBs.

Geochemical and isotopic characteristics of volcanic rocks from the northern East China Sea shelf margin and the Okinawa Trough

Volcanic rocks both from the northern East China Sea (NECS) shelf margin and the northern Okinawa Trough are subalkaline less aluminous, and lower in High Field Strength Elements (HFSE). These rocks are higher in Large Ion Lithophile Elements (LILE), thorium and uranium contents, positive lead anomalies, negative Nb-Ta anomalies, and enrichment in Light Rare Earth Elements (LREE). Basalts from the NECS shelf margin are akin to Indian Ocean Mid-Ocean Ridge Basalt (MORB), and rhyolites from the northern Okinawa Trough have the highest Pb-207/Pb-208 and Pb-208/Pb-204 ratios. The NECS shelf margin basalts have lower Sr-87/Sr-86 ratios, epsilon(Nd) and sigma O-18 than the northern Okinawa Trough silicic rocks. According to K-40-Ar-40 isotopic ages of basalts from the NECS shelf margin, rifting of the Okinawa Trough may have been active since at least 3.65-3.86 Ma. The origin of the NECS shelf margin basalt can be explained by the interaction of melt derived from Indian Ocean MORB-like mantle with enriched subcontinental lithosphere. The basalts from both sides of the Okinawa Trough may have a similar origin during the initial rifting of the Okinawa Trough, and the formation of basaltic magmas closely relates to the thinning of continental crust. The source of the formation of the northern Okinawa Trough silicic rocks was different from that of the middle Okinawa Trough, which could have been generated by the interaction of basaltic melt with an enriched crustal component. From the Ryukyu island arc to East China, the Cenozoic basalts have apparently increasing trends of MgO contents and ratios of LREE to Heavy Rare Earth Elements (HREE), suggesting that the trace element variabilities of basalts may have been influenced by the subduction of the Philippine Sea plate, and that the effects of subduction of the Philippine Sea plate on the chemical composition of basaltic melts have had a decreasing effect from the Ryukyu island arc to East China.