Understanding the Marginal Seas of Northeast Asia for Tectonics and Submarine Geohazards

This study focuses on revisiting the tectostratigraphic framework of the Ulleung Basin and conceptualizing neotectonics around the western East Sea margin. Based on the analysis of 2D and 3D multi-channel seismic reflection data and offshore drill wells, we divided the entire sedimentary successions of the Ulleung Basin into four tectostratigraphic sequences, named TS1 (c. 23–16 Ma), TS2 (c. 16–9 Ma), TS3 (c. 9–4 Ma), and TS4 (c. 4 Ma–present), in ascending order. The results show that each sequence has been deformed once or multiple times in different periods by juxtaposing two major compressional structures named the Dolgorae Thrust-Fold Belt and the Gorae Anticline. Interpretation of the stratal deformation and termination patterns of the syn- and post-deformational sequences of each structures suggests that the thrusting and folding of the Dolgorae Thrust-Fold Belt was active from c. 16 Ma to c. 9 Ma under the NNW–SSE compressional stress regime (Stage-2), whereas the Gorae Anticline was active from 4 Ma to the present under the ENE–WSW compressional stress regime (Stage-4). Between these two compressional events, there was an intervening period of regional slow subsidence driven by thermal contraction of the back-arc lithosphere and isostatic sedimentary loading (Stage-3). Based on the stratigraphic and structural reconstruction, we propose a 4-stage tectonic model: Stage-1) back-arc opening stage associated with the southward drift of the Japanese islands (c. 23–16 Ma), Stage-2) tectonic-inversion stage in association with the reorganization of the Pacific and Philippine Sea plates and clockwise rotation of SW Japan (c. 16–9 Ma), Stage-3) post-inversion stage with regional thermal and isostatic subsidence (c. 9–4 Ma), and Stage-4) neotectonic stage in which embryonic subduction is nucleating on the East Sea margins under the E–W compressional stress regime (c. 4 Ma–present).

Scholte-wave dispersion analysis is effective at imaging the relatively low shear-wave velocity of shallow marine sediments in marginal seas. The combination of a four-component ocean-bottom-seismometer (OBS) and a towed air-gun source can economically and effectively acquire the marine dispersive seismic data. Extracting higher-order dispersive Scholte wave modes is the most critical problem in the dispersion analysis method. The extremely low shear-wave velocity and severe attenuation in the top hundreds of meters of marginal sea sediment provide an uneven dispersive energy distribution for the four components of the Scholte wave data. The fundamental mode dispersive energy dominates in the vertical component and higher-order modes dominate in the horizontal component. We developed the method of the four-component OBS Scholte velocity-spectra stacking, which can effectively, rapidly, and robustly extract higher-order modes. We imaged the shear-wave velocity structure of complicated shallow marine sediment in the North Yellow Sea using an active OBS seismic profile with a large-volume air-gun array. The fourth higher-order Scholte wave mode can be imaged with the four-component velocity-spectra stacking method with a lower frequency range of 1.0–7.0 Hz. Only the second-order mode can be recognized from the dispersion energy image of the single vertical component. The joint inversion of multimode dispersion curves can provide more accuracy and deeper constraints for the inverted model; thus, the constraint depth with five modes increases by a factor of 1.9 compared with single fundamental mode inversion. The inverted profile suggests a low shear-wave velocity of 123–670 m/s and strong lateral variations within 350 m. The main regional geological structures are shown by the inverted shear-wave velocity structure.

As an important segment of the North China Craton, the Trans-North China Orogen (TNCO) has experienced strong tectonic deformation and magmatic activities since the Cenozoic and is characterized by significant seismicity. To understand the mechanism of the crustal deformation and seismic hazards, we determined the crustal thickness (H), Vp/Vs ratio (κ) and crustal anisotropy (the fast polarization direction φ and splitting time τ) beneath the TNCO and its adjacent areas by analyzing receiver function data recorded by a dense seismic array. The (H, κ) and (φ, τ) at a total of 309 stations were measured, respectively. The Moho depth varies from ∼30 km beneath the western margin of the Bohai bay basin to the maximum value of ∼48 km beneath the northern Lüliang Mountain, which shows the positive and negative correlations with the elevation and the Bouguer anomaly. The average φ is roughly parallel to the strikes of the faults, grabens and Mountains in this study area, whereas a rotating distribution is shown around the Datong-Hannuoba volcanic regions. Based on the φ measured from the Moho Ps and SKS/SKKS phases, we propose that the crustal deformation and seismic hazards beneath the TNCO could be due to the counterclockwise rotation of the Ordos block driven by the far-field effects of the India-Eurasian collision.