Xiongwei Niu, Pingchuan Tan, Weiwei Ding, Wei Wang, Yao Wei, Xiaodong Wei, Aiguo Ruan, Jie Zhang, Chunyang Wang, Yong Tang, Jiabiao Li. Oceanic crustal structure and tectonic origin of the southern Kyushu-Palau Ridge in the Philippine Sea[J]. Acta Oceanologica Sinica, 2022, 41(1): 39-49. doi: 10.1007/s13131-021-1978-9
Citation: Wei Qinsheng, Wang Baodong, Yao Qingzhen, Yu Zhigang, Fu Mingzhu, Sun Junchuan, Xu Bochao, Xie Linping, Xin Ming. Physical-biogeochemical interactions and potential effects on phytoplankton and Ulva prolifera in the coastal waters off Qingdao (Yellow Sea, China)[J]. Acta Oceanologica Sinica, 2019, 38(2): 11-23. doi: 10.1007/s13131-019-1344-3

Physical-biogeochemical interactions and potential effects on phytoplankton and Ulva prolifera in the coastal waters off Qingdao (Yellow Sea, China)

doi: 10.1007/s13131-019-1344-3
  • Received Date: 2017-07-02
  • In recent years, the spectacular massive green tide of Ulva prolifera has become a recurrent phenomenon appearing every summer in the coastal waters off Qingdao (Yellow Sea, China), attracting the attention of scientists and local government. Based on multidisciplinary data collected during summer and winter, this study focuses on the hydrological characteristics and regional biogeochemical processes in coastal waters off Qingdao. The results show that the boundary of the Yellow Sea Cold Water Mass (YSCWM) can reach the Qingdao coastal region in summer and is locally raised to the upper layers to form coastal upwelling beyond tidal mixing and favorable wind. The regional summer upwelling off the Qingdao coast effectively enriches the nutrient concentrations in the upper water column and thus promotes growth of phytoplankton but reduces the dissolved oxygen (DO) concentration and pH value in the bottom. The regional summer upwelling off Qingdao coast may facilitate the growth and regional blooming of the U. prolifera that migrate to this region with the southerly wind. Additionally, the effects of the front on the aggregation of U. prolifera may be significant. In winter, the Yellow Sea Warm Current (YSWC) extends and spreads along the offshore region off the Subei Shoal towards the Qingdao coastal sea. This tongue-shaped warm water meets the cold coastal water off Qingdao, which leads to the formation of a physical front. As a consequence, remarkable fronts of nutrient and chlorophyll a (Chl a) also form between the shoreward warm water and the cold coastal water. This study increases the understanding of the interactions between the regional physical, chemical, and biological processes off the Qingdao coast.
  • The Izu-Bonin-Marina arc (IBM) subduction system is a typical representative of the large-scale convergence boundary in the world, having strong lithospheric deformation and magmatic activity in the West Pacific. It forms a typical “trench-arc-basin” system with an eastward transition of the back-arc spreading due to the westward subducting Pacific Plate (Fig. 1). It is widely considered as an ideal area to study the “subduction factory” as well as the seismic activities of the subduction zone (Ishizuka et al., 2011a, b). In the western part, the IBM subduction system has a remnant arc, namely the Kyushu-Palau Ridge (KPR); whereas the modern IBM volcanic island arc belt (such as the Izu-Ogasawara Ridges and Mali Yana Island Arc) lies in the east. Several Cenozoic oceanic basins developed in-between, including the Shikoku Basin, Parece Vela Basin (PVB), Ogasawara Trough, and Mariana Trough (Fig. 1).

    Figure  1.  Bathymetry (a, Tozer et al., 2019) and free air gravity (1 Gal=1 cm/s2) (b, Sandwell et al., 2014) of the Mariana arc-back-arc region with locations of the wide-angle seismic refraction profiles, where the study area is indicated by a red dashed rectangle. The black solid lines show the profiles from Nishizawa et al. (2016), while the red solid line shows the profile of this study.

    The KPR extends from the Nankai Trough in the north to the northern part of the Palau Archipelago in the south. It is a bathymetric high extending through the Philippine Sea Plate with a total length of 2 600 km (Fig. 1). It was once connected with the modern IBM island arc belt before the splitting of the proto-IBM arc and the back-arc spreading with the continuous retreat of the Pacific Plate subduction at about 30 Ma BP (Stern et al., 2004). Therefore, the KPR should be firstly built by the initial subduction between the proto-Pacific Plate under the West Philippine Basin (WPB), and then modified by the magmatic activity after the breakup of the proto-IBM island arc. However, how magmatism controlled the growth of the KPR is still debatable.

    Seismic velocity models can be used to provide a crisis information for the mechanisms of the growth of the island arc crust. The previous deep reflection/refraction seismic explorations carried out on the KPR were completed mainly by Japanese scientists (Murauchi et al., 1968; Arisaka et al., 2003; Nishizawa et al., 2007, 2016). Especially during 2004–2008, the Japan Agency for Marine-Earth Science and Technology conducted a series of seismic experiments on the KPR and presented 27 wide-angle seismic profiles to understand its lithospheric structural changes from north (31°N) to south (13°N) (Nishizawa et al., 2016). Their results showed that the KPR has a large variation in crustal thickness (8–23 km), which is generally larger than those oceanic crust of the surrounding back-arc basin, such as the PVB in the east and the WPB in the west. The growth of the thick crust is accommodated mainly by accretion of the lower crust (6.8–7.0 km/s), whereas for crustal thickness greater than 20 km, it also has a relatively thick (>5 km) middle crust with a P-wave velocity ($V_{\rm{p}} $) of ~6 km/s. Some high velocity lower-crust layers are also observed under the KPR with a relatively high $V_{\rm{p}} $ of 7.2 km/s (Nishizawa et al., 2016). Two wide-angle seismic profiles cross the southern part of the KPR (south of 15°N) (Figs 1 and 2), but show a different crustal structure (Nishizawa et al., 2016). One profile (ID: Kpr40, near 14°N) indicate that the KPR crust has a crustal “root” with an average lower crustal velocity of 7.2 km/s. The crust is 15 km thick, and its thickening is attributed mainly to the accretion of the middle and lower crust. The other profile crossing the KPR ridge (ID: Kpr41, near 13 °N) lacks the high velocity at lower crustal layer and crustal “root” with a very thin (<2 km) middle crustal layer. The crustal thickness is 10 km (Nishizawa et al., 2016).

    Figure  2.  Seismic refraction profiles shown on bathymetry (a) and free-air gravity (1 Gal=1 cm/s2) (b) of the southern part of the KPR. The red line is the ocean bottom seismograph (OBS) profile KPR2020-2 in this study, whose circles show OBS positions, and yellow fill indicate useful data. The numbers show the OBS positions represented in Figs 3−7. Black solid lines (Kpr40, Kpr41) represent the seismic refraction profiles from Nishizawa et al. (2016).

    According to the above, seismic experiments focused mainly on the northern part of the KPR (north of 15°N) (Fig. 1), and the crustal structure of the southern part (south of 15°N) has limited data and is less understood. To better understand the crustal structure of the southern KPR, this study used a 217 km long wide-angle seismic profile (ID: KPR2020-2) to estimate the structure of the crust and the uppermost mantle across the KPR near 12°N. Comparing it with the crustal structure of the above two profiles will provide key information on the growth of the KPR in this area.

    Based on petrological and geomagnetic data, the tectonic evolution history of the Philippine Sea Plate has been well understood (Stern and Bloomer, 1992; Bloomer et al., 1995; Kobayashi et al., 1995; Okino et al., 1999; Hall, 2002; Stern et al., 2004). The first stage is the formation of the WPB, which is a back-arc basin developed with the seafloor spreading starting at around 55 Ma BP (Deschamps and Lallemand, 2002). Then, during 50–52 Ma BP, the Paleo-Pacific Plate began to subduct beneath the WPB, forming initial lithosphere subsidence along the eastern edge of the WPB (Ishizuka et al., 2011a). At 43 Ma BP, when the abrupt change of the Pacific Plate movement from NNW to NWW occurred, the subsidence of the lithosphere evolved into subduction (Richards and Lithgow-Bertelloni, 1996), indicating that the initial subduction lasted about 7 Ma (Ishizuka et al., 2011a). The International Ocean Discovery Program (IODP) has successively carried out two expeditions in this area (IODP Expeditions 351 and 352). Their results indicated that the fore-arc seafloor spreading occurred before the initiation of the subduction, then with the continuous subduction of the proto-Pacific Plate, the proto-IBM arc was formed within the fore-arc oceanic basin (Arculus et al., 2015; Reagan et al., 2017; Ishizuka et al., 2018). Based on the age and geochemical characteristics of the dredged samples along the entire KPR, Ishizuka et al. (2011b) concluded that the activity timing of KPR is 25–48 Ma BP, but most ages of the dredged samples fall into a narrow range of 25–28 Ma BP. The proto-IBM arc is dominated by island arc porphyry and calc-alkaline basalt (Ishizuka et al., 2018). After the formation of the proto-IBM arc, rifting and seafloor spreading occurred in the proto-IBM arc at about 30 Ma BP, forming the Shikoku Sea Basin (30–15 Ma BP; Okino et al., 1994) and the PVB (29–12 Ma BP; Okino et al., 1999) between the residual KPR and the present IBM arc (Fig. 1). These two spreading systems connected at 20 Ma BP and ceased spreading around 15 Ma BP. The spreading center jumped eastward with the retreat of the Pacific Plate, which opened the Mariana Trough at 3–4 Ma BP (Stern et al., 2004).

    The seismic source was a 4×32.7 L air-gun array shot at an average shooting interval of 200 m with a pressure of 10.79 MPa every 78–97 s. Here, this study focuses on the 217 km seismic profile KPR2020-2 with 1 083 shots (red line in Fig. 2). The OBS is composed of 4.5 Hz four-component geophones. Ten OBSs recorded good data sets, whereas OBS1 and OBS2 in the west and OBS13 in the east failed (Fig. 2). This study determined P wave velocities using the hydrophone component. For OBS5, the vertical component was used, as poor signals were recorded in its hydrophone component. The processing of OBS data included correction for OBS clock drift, relocation of OBSs and shots using direct arrivals, and bandpass filtering at 4–20 Hz. A final straight-line approximation of the profile was determined by least squares fit to all the shots on the profile. The depth of each OBS was first roughly determined from multi-beam bathymetric data, and then refined by fitting direct arrivals.

    In general, the seismic data quality is good. Phases are identified using initial travel time modeling. This study identified the direct water wave labeled Pw, the refracted wave from oceanic Layer 2 labeled P2, the refracted wave from oceanic Layer 3 labeled P3, the Moho reflection labeled PmP, and the refracted wave from the upper mantle labeled Pn.

    Clear Pw, P2, P3, PmP and Pn phases were observed at most OBSs, especially OBSs5–11 deployed above the KPR (Figs 37). A maximum offset of the Pn phase was observed up to 85 km, which made a good overlap control for the model. Travel-time uncertainties were dominated by picking uncertainties, which were estimated as 50 ms, 60 ms, 80 ms, 100 ms and 120 ms for Pw, P2, P3, PmP and Pn arrivals, respectively.

    Figure  3.  An original seismic record section of the vertical component of OBS7 (a), the record section with picked and calculated travel time overlain (b), and a simulation of ray-tracing (c). T-D represents vertical component data of OBS. In these diagrams, the reduction velocity is 8.0 km/s. The names of phases are explained in the text. In b, black dots represent the predicted travel time, and the colored vertical bars represent the observed travel time in the same color of rays in c. The size of the vertical bars indicates twice the uncertainty (Zelt and Smith, 1992). In c, the colored lines represent the ray paths of different phases. The black dashed lines represent the seabed, the interface between oceanic Layer 2 and Layer 3, and the Moho discontinuity respectively, from top to bottom.
    Figure  4.  An original seismic record section of the vertical component of OBS8 (a), the record section with picked and calculated travel time overlain (b), and a simulation of ray-tracing (c). T-D represents vertical component data of OBS. In these diagrams, the reduction velocity is 8.0 km/s. The names of phases are explained in the text. In b, black dots represent the predicted travel time, and the colored vertical bars represent the observed travel time in the same color of rays in c. The size of the vertical bars indicates twice the uncertainty (Zelt and Smith, 1992). In c, the colored lines represent the ray paths of different phases. The black dashed lines represent the seabed, the interface between oceanic Layer 2 and Layer 3, and the Moho discontinuity respectively, from top to bottom.
    Figure  5.  An original seismic record section of the vertical component of OBS9 (a), the record section with picked and calculated travel time overlain (b), and a simulation of ray-tracing (c). T-D represents vertical component data of OBS. In these diagrams, the reduction velocity is 8.0 km/s. The names of phases are explained in the text. In b, black dots represent the predicted travel time, and the colored vertical bars represent the observed travel time in the same color of rays in c. The size of the vertical bars indicates twice the uncertainty (Zelt and Smith, 1992). In c, the colored lines represent the ray paths of different phases. The black dashed lines represent the seabed, the interface between oceanic Layer 2 and Layer 3, and the Moho discontinuity respectively, from top to bottom.
    Figure  6.  An original seismic record section of the vertical component of OBS10 (a), the record section with picked and calculated travel time overlain (b), and a simulation of ray-tracing (c). T-D represents vertical component data of OBS. In these diagrams, the reduction velocity is 8.0 km/s. The names of phases are explained in the text. In b, black dots represent the predicted travel time, and the colored vertical bars represent the observed travel time in the same color of rays in c. The size of the vertical bars indicates twice the uncertainty (Zelt and Smith, 1992). In c, the colored lines represent the ray paths of different phases. The black dashed lines represent the seabed, the interface between oceanic Layer 2 and Layer 3, and the Moho discontinuity respectively, from top to bottom.
    Figure  7.  An original seismic record section of the vertical component of OBS11 (a), the record section with picked and calculated travel time overlain (b), and a simulation of ray-tracing (c). T-D represents vertical component data of OBS. In these diagrams, the reduction velocity is 8.0 km/s. The names of phases are explained in the text. In b, black dots represent the predicted travel time, and the colored vertical bars represent the observed travel time in the same color of rays in c. The size of the vertical bars indicates twice the uncertainty (Zelt and Smith, 1992). In c, the colored lines represent the ray paths of different phases. The black dashed lines represent the seabed, the interface between oceanic Layer 2 and Layer 3, and the Moho discontinuity respectively, from top to bottom.

    OBS7 was deployed at the western part of the KPR. Clearly, P2, P3, PmP and Pn phases were observed at both sides of OBS7, PmP phase was found at offsets about 20 km and 40 km of its western and eastern sides, respectively (Fig. 3). OBS8 was deployed 15 km east of OBS7. Clearly, P2, P3, PmP and Pn phases were found at its both sides, and PmP phase was found at offsets starting at 20 km and 33 km of its western and eastern sides, respectively (Fig. 4). OBS9 was deployed at the center of the KPR, 15 km east of OBS8, and P2, P3, and Pn phases were also found at both its sides, but PmP phase was found only at offsets starting at 35 km of its eastern side (Fig. 5). OBS10 was deployed at 15 km east of OBS9. Clearly, PmP phase was observed at offsets of 30 km and 15 km at its western and eastern sides, respectively (Fig. 6). Although PmP phase was found only at the western side of OBS11 with an offset starting at about 30 km (Fig. 7), P2, P3, and Pn phases were observed clearly at OBS10 and OBS11. A few arrivals may not be fitted well, which may be due to the rugged topography across the KPR.

    Travel time modelling and inversion approach of Zelt and Smith (1992) were used to construct velocity models. The initial model consisted of two crustal layers, representing oceanic layers 2 and 3, with linear gradients within these layers and no velocity discontinuity between them. Oceanic Layer 2 was 2 km thick and had velocities of 1.8 km/s at the top and 6.4 km/s at the bottom, whereas oceanic Layer 3 was 4 km thick and had velocities of 6.4 km/s at the top and 7.0 km/s at the bottom. The upper mantle had a velocity of 8.0 km/s at the top. The horizontal node spacings within oceanic Layers 2 and 3 were 5 km. In the upper mantle, the node spacing was 50 km. This model was based on global averages of the oceanic crustal structure (Kennett, 1982; White et al., 1992). As described by Zelt and Smith (1992), velocities and boundaries were initially adjusted manually by trial and error. After this forward modeling, a damped least squares inversion was conducted to optimize the velocity within each layer, while keeping the layer boundaries fixed.

    Along the profile KPR2020-2, the average thicknesses of oceanic crust Layer 2 and oceanic crust Layer 3 were ~3.5 km and 4.7 km, respectively, and the average velocities in Layer 2 and Layer 3 were 4.2 km/s and 6.7 km/s, respectively. As a result, the vertical velocity gradient within Layer 2 was high (average is 1.2 s−1), while within Layer 3 was low (average is 0.2 s−1). The velocity at the top of the upper mantle was 7.8 km/s. According to the tomography features and internal structures, this study divided the seismic profile into three parts: WPB (0–70 km), KPR (70–170 km), PVB (170–217 km).

    Beneath the WPB, Layers 2 and 3 had roughly constant thickness (2.7 km and 3.8 km, respectively), and different vertical velocity gradients (~1.6 s−1 and ~0.2 s−1, respectively). In general, lower crustal velocity of 7.0 km/s was observed at the eastern margin (60–70 km) of the WPB with few high-velocity lower crustal layers (7.3 km/s) (Fig. 8a). The Moho lies at ~12.2 km below sea level, resulting in a crustal thickness of 6.5 km. Part of the Moho reflector was controlled well by PmP (Fig. 8a).

    Figure  8.  The final velocity model along the profile (a) and the distribution of the ray density along the profile (b), the size of the statistic network is 0.5 km×0.2 km. Numbered red triangles represent the OBS stations. Thin black lines represent the velocity contours every 0.4 km/s. Thick black lines represent the seabed and the Moho discontinuity. On the Moho discontinuity, the red colored portions indicate sections constrained the PmP reflections. The thin dotted black line represents the interface between oceanic Layers 2 and 3. The thin dashed line represents P wave velocity contour of 4.0 km/s. The white lines with notes mark the crust thicknesses.

    Beneath the major part of the KPR, both crustal thickness and velocity show lateral heterogeneous. Thickness of Layer 2 changed little (between 3.7 km and 4.2 km), although the seafloor depth changed a lot (from 2.5 km to 5.0 km). The thickness of Layer 3 varied from 4.0 km at the edges of the KPR to 6.7 km in the center. Velocity varied at the top of Layer 2 (seafloor, from 1.6 km/s to 3.2 km/s) and at the bottom of Layer 3 (6.7 km/s to 7.2 km/s), and the velocity between Layers 2 and 3 showed no discontinuities at 6.1 km/s, which is lower than that beneath the WPB and PVB (6.4 km/s). The depth of the Moho reflector was well controlled by PmP (Fig. 8a), which varied from 12.7 km to 14.0 km below sea level, sketching a relatively thicker crust ranging from to 8.2 km to 11.3 km.

    The crustal structure beneath the PVB was found to be similar to that of the WPB, both Layers 2 and 3 have roughly constant thickness (3.0 km and 4.4 km, respectively), but different vertical velocity gradients (~1.5 s−1 and ~0.2 s−1, respectively). The average $V_{\rm{p}} $ of the lower crust along the western margin of the PVB was less than 7 km/s. The main difference between the PVB and the WPB was the thickness of Layer 3, the Layer 3 was ~0.6 km thicker in the PVB.

    Pick uncertainties and the travel-time misfit of the best model are provided in Table 1. The overall root-mean-square (RMS) misfit was 113 ms, and the misfit increased slightly with depth and offset. The overall χ2 value was 1.443, with 95.0% of picks fitted, and the normalized travel time misfit suggested that the model was suitably parameterized, with a travel-time misfit a slightly bigger than the pick uncertainty. The numbers of rays through each cell were generally larger than 10 and reached over 100 (Fig. 8b), indicating that the model was well constrained. By using a perturbation test (Zelt and Smith, 1992; Muller et al., 1997), the uncertainty of the Moho depth was found to be ±0.3 km.

    Table  1.  Statistics of travel-time analysis
    PhaseTotal picksInverted picksFit ratio/%RMS/ms$\ \chi^2$
    Pw45644798.0360.521
    P260454389.91012.813
    P32 1001 92991.9981.528
    PmP1 0691 05398.5900.913
    Pn2 1262 12499.91461.476
    All6 4166 09495.01131.443
    Note: RMS, root-mean-square.
     | Show Table
    DownLoad: CSV

    Two-dimension (2D) gravity modeling and isostatic analysis can give further constraints on the crustal model. The forward gravity response of the 2D gravity modeling is based on the 2D algorithm by Talwani et al. (1959). The free-air gravity anomaly is derived from Sandwell et al. (2014), which is a global 1′ × 1′ gravity grid. The crustal and upper mantle density model was made from the P wave velocity model in this study, where the densities for each polygon were derived from RAYINVR using its internal velocity to density conversion (Zelt and Smith, 1992). The result showed a reasonable agreement between the calculated and observed gravity anomalies along the WPB and PVB, whereas the calculated gravity under the KPR was higher than the observed (20–40 mGal, 1 Gal=1 cm/s2) (Fig. 9). This difference might be due to the low resolution of the long-wavelength satellite derived gravity anomalies. However, the overall trend of the calculated gravity filed fits well with the observed trend, which makes good constraints on the crustal structure of the seismic refraction profile.

    Figure  9.  Results of 2D gravity modeling of the profile KPR2020-2 from 40 km to 217 km. a. The pressure is calculated at the bottom of the model at 20 km depth; b. the results of the calculated gravity (red solid line) and observed gravity from Sandwell et al. (2014) (black dots); c. the crustal density model derived from the crustal velocity modeling. The ray density along the profile (Fig. 8b) is added as background picture in c.

    For isostatic analysis, this study calculated the pressure based on the crustal and upper mantle density model down to depth of 20 km (Fig. 9). The pressures along the WPB and PVB were 0.485 GPa and 0.495 GPa, respectively. There was no regional trend along both WPB and PVB. The difference between these two is probably due to the different oceanic lithospheric densities caused by the different ages of oceanic crust there. The pressure under the ridge was 0.51–0.52 GPa, 5.1% higher, suggesting that the KPR grew on oceanic lithosphere with some flexural strength.

    Early studies have shown that the velocity structure of the mature intra-oceanic island arc crust is quite different from that of submarine volcanoes that grow on the oceanic crust. The main difference is that the island arc has a thicker middle crust larger than 5 km (6.0–6.5 km/s) (Holbrook et al., 1999; Takahashi et al., 2008). The growth of the island arc results mainly from the accretion of both the middle and lower crust. The first stage includes the successive basaltic underplating caused by the partial melting of the underneath subduction plate. The second stage involves the formation of the middle crust, which comes from the initial basaltic crust differentiating to generate a crust with intermediate to felsic components with higher SiO2 content (Takahashi et al., 2007, 2008).

    Previous studies show that the development of the southern KPR has two stages (Stern et al., 2004). The first stage includes early initial island arc formation caused by the subduction of the paleo-Pacific Plate under the WPB starting at 45 Ma BP (Ishizuka et al., 2011a), forming the proto-IBM arc. The second stage involves the splitting of the proto-IBM arc and back-arc seafloor spreading of the PVB between 30–28 Ma BP (Stern et al., 2004; Okino et al., 1994; Ishizuka et al., 2011a; Grevemeyer et al., 2021), producing the relic KPR and the present IBM island arc. If the KPR is composed of the initial island arc, the crustal structure should resemble a velocity structure of the mature arc island. However, if the growth of the KPR is dominated by the magmatism related to the later stage back-arc spreading, the original velocity structure of the KPR could be modified, which could be similar to the structure of typical oceanic submarine volcanoes.

    Based on the seismic velocity model, the KPR crust of this study can be divided into the upper crust (<6.4 m/s), and lower crust (6.4–7.2 km/s). The crustal thickness of the KPR is up to 12 km, but it gradually decreases to 5–6 km on both sides adjoining the WPB and PVB. The crustal structures of the WPB and PVB are consistent with results obtained from earlier seismic studies in this area (Arisaka et al., 2003; Nishizawa et al., 2007, 2016). The crustal velocity structure under the KPR is comparable to the profile Kpr41, but shows remarkable difference with that of the profile Kpr40. The latter has a thick (>5 km) middle crustal layer with a thick (~5 km) high-velocity lower crustal layer of 7.2 km/s. (Nishizawa et al., 2016) (Fig. 2). Compared to the Atlantic Ocean 0–127 Ma oceanic crustal structures, the velocity structure under the KPR is similar to that of the oceanic crust (White et al., 1992) (Fig. 10). The thick middle crust (>5 km) that observed in the mature intra-oceanic island arc crust was not observed in the velocity model of this study, which indicates that the KPR in the study area has immature island arc nature.

    Figure  10.  A comparison of the 1D velocity curves beneath each OBS station, continental crust and Atlantic Ocean 0–127 Ma oceanic crust. In a, bold green solid line represents extended continental crust, shadow zone represents average continental crust (Christensen and Mooney, 1995). In b, shadow zone represents Atlantic Ocean 0–127 Ma oceanic crust (White et al., 1992).

    Back-arc spreading usually accompanies arc development (Martinez and Taylor, 2002; Dunn and Martinez, 2011). Consequently, magmatism during the arc splitting process could be controlled by pressure-release partial melting of fertile mantle. The crustal structure formed by this process is similar to the velocity structure of normal oceanic crust.

    The subducting slabs could also generate arc magmatism affecting KPR crustal accretion during the breakup of the arc island. In this case, it will form a thick high-velocity lower crustal layer (7.2–7.6 km/s) along the western margin of the PVB, which mainly has mafic-to-ultramafic crustal composition underplated the lower crustal layer (Eason and Dunn, 2015). However, the seismic velocity structure shows that the average $V_{\rm{p}} $ of the lower crust along the western margin of the PVB is less than 7 km/s. In addition, the crustal thickness changes gradually from the KPR to the PVB, indicating that the rifting of the KPR is accompanied by a large amount of magmatism. Thus, this study suggest that the absence of the middle crust and thick high-velocity lower crustal layers in the southern KPR might due to the limited domination of arc magmatism, and/or the modification by the later stage magmatism related with arc rifting.

    The main formation stage of the KPR can be estimated by the plate strength at the time of the KPR loading. The IODP results show that the initial stage of the KPR formation is built by the arc magmatism (~45 Ma BP), which grew on the oceanic basin created by the fore-arc seafloor spreading starting at 50–52 Ma BP (Arculus et al., 2015; Reagan et al., 2017; Ishizuka et al., 2018). If most of the crust accretion of the KPR comes from this process, then the KPR should be formed on an oceanic plate (oceanic age 5–7 Ma) with low lithospheric plate strength. In this case, the crust should have a mountain root with nearly local isostasy compensation. At around 30 Ma BP, the rejuvenate of the KPR grew on an oceanic plate with age older than 20 Ma. If that so, the crustal structure should have some regional compensation as the ridge grew on an oceanic plate with some flexural strength. The crustal model of this study showed that the pressure under the KPR was 0.51–0.52 GPa, which is 5.1% higher than the adjacent PVB and WPB. This implies that there is some regional compensation during the loading of the KPR. Thus, this study exclude the effects of initial arc magmatism during the formation of the KPR.

    In conclusion, the seismic crustal velocity structure shows that the velocity structure of KPR is more similar to an oceanic crust rather than a mature island arc. The estimation of the plate strength during the KPR loading suggests that the KPR is built mainly by the magmatism during the rifting of the proto-IBM arc between 30–28 Ma BP. The absence of the middle crust layer, larger than 5 km (6.0–6.5 km/s) and high velocity lower crustal layers (7.2–7.6 km/s) indicate that the arc magmatism plays a less important role in the KPR formation.

    Travel-time modeling of active-source wide-angle seismic P wave data from OBSs deployed crossing the WPB, KPR, and PVB led to the following conclusions:

    (1) This study show a high resolution crustal $V_{\rm{p}} $ structures beneath the WPB, KPR, and PVB. Beneath the KPR, $V_{\rm{p}} $ in crustal Layers 1–2, and crustal Layer 3 are 1.6–6.1 km/s and 6.1–7.2 km/s, respectively. The crustal thickness of the KPR is up to 12 km, but gradually decreases to 5–6 km on both sides adjoining the WPB and PVB. Similar crustal structures are also observed beneath the adjacent WPB in the west and PVB in the east. Only a little lateral velocity heterogeneity of less than 0.3 km/s is seen in the crust. 2D gravity modeling shows that the overall trend of the calculated gravity filed is reasonably well fit with the observed, which confirms the crustal structure of the seismic refraction profile.

    (2) 1 D velocity tests show that the seismic velocity structure of the KPR is comparable to the Atlantic oceanic crust (0–127 Ma), which implies a typical oceanic crust beneath the KPR. Isostatic analysis shows that there is some regional compensation during the loading of the KPR, indicating that the initial arc-magmatism is not the dominant process of the KPR formation. The absence of the thick middle crust (6.0–6.5 km/s) and thick high-velocity lower crust layer (7.2–7.6 km/s) suggest that arc-magmatism plays a less important role in the formation of the KPR.

    We thank the crew of the R/V Dayanghao, and the technical staff and researchers involved in the West Philippine Sea 2020 Survey.

  • Bao Min, Guan Weibing, Yang Yang, et al. 2015. Drifting trajectories of green algae in the western Yellow Sea during the spring and summer of 2012. Estuarine, Coastal and Shelf Science, 163:9-16, doi: 10.1016/j.ecss.2015.02.009
    Bao Shaowu, Li Xiaofeng, Shen Dongliang, et al. 2017. Ocean upwelling along the Yellow Sea coast of China revealed by satellite observations and numerical simulation. IEEE Transactions on Geoscience and Remote Sensing, 55(1):526-536, doi: 10.1109/TGRS.2016.2610761
    Castelao R M, Wang Yuntao. 2014. Wind-driven variability in sea surface temperature front distribution in the California Current System. Journal of Geophysical Research:Oceans, 119(3):1861-1875, doi: 10.1002/2013JC009531
    Chen C T A. 2009. Chemical and physical fronts in the Bohai, Yellow and East China seas. Journal of Marine Systems, 78(3):394-410, doi: 10.1016/j.jmarsys.2008.11.016
    Ekman V W. 1905. On the influence of the earth's rotation on ocean-currents. Arkiv för Matematik, Astronomi Och Fysik, 2(11):1-53
    Fu Mingzhu, Wang Zongling, Li Yan, et al. 2009. Phytoplankton biomass size structure and its regulation in the Southern Yellow Sea (China):Seasonal variability. Continental Shelf Research, 29(18):2178-2194, doi: 10.1016/j.csr.2009.08.010
    Gao Shengquan, Lin Yi'an, Jin Mingming, et al. 2003. Distribution of nutrient and its relationship with anchovy spawning ground in the southern waters of Shandong Peninsula. Haiyang Xuebao (in Chinese), 25(Suppl 2):157-166
    Ge Renfeng, Qiao Fangli, Yu Fei, et al. 2003. A method for calculating thermocline characteristic elements in shelf sea area-Quasi-step function approximation method. Advances in Marine Science (in Chinese), 21(4):393-400
    Guo Binghuo, Hu Xiaomin, Xiong Xuejun, et al. 2003. Study on interaction between the coastal water, shelf water and Kuroshio water in the Huanghai Sea and East China Sea. Acta Oceanologica Sinica, 22(3):351-367
    He Chongben, Wang Yuanxiang, Lei Zongyou, et al. 1959. A preliminary study of the formation of Yellow Sea Cold Water Mass and its properties. Oceanologia et Limnologia Sinica (in Chinese), 2(1):11-15
    Hu Chuanmin, Li Daqiu, Chen Changsheng, et al. 2010. On the recurrent Ulva prolifera blooms in the Yellow Sea and East China Sea. Journal of Geophysical Research, 115:C05017
    Huang Daji, Fan Xiaopeng, Xu Dongfeng, et al. 2005. Westward shift of the Yellow Sea warm salty tongue. Geophysical Research Letters, 32:L24613, doi: 10.1029/2005GL024749
    Hu Jianyu, Wang Xiaohua. 2016. Progress on upwelling studies in the China seas. Reviews of Geophysics, 54:653-673, doi: 10.1002/2015RG000505
    Hwang J H, Van Sy P, Choi B J, et al. 2014. The physical processes in the Yellow Sea. Ocean & Coastal Management, 102:449-457
    Lee J H, Pang I C, Moon I J, et al. 2011. On physical factors that controlled the massive green tide occurrence along the southern coast of the Shandong Peninsula in 2008:A numerical study using a particle-tracking experiment. Journal of Geophysical Research, 116:C12036, doi: 10.1029/2011JC007512
    Leliaert F, Zhang Xiaowen, Ye Naihao, et al. 2009. Research note:Identity of the Qingdao algal bloom. Phycological Research, 57(2):147-151, doi: 10.1111/pre.2009.57.issue-2
    Lie H J, Cho C H. 2016. Seasonal circulation patterns of the Yellow and East China Seas derived from satellite-tracked drifter trajectories and hydrographic observations. Progress in Oceanography, 146:121-141, doi: 10.1016/j.pocean.2016.06.004
    Lie H J, Cho C H, Lee S. 2009. Tongue-shaped frontal structure and warm water intrusion in the southern Yellow Sea in winter. Journal of Geophysical Research, 114:C01003
    Lin Xiaopei, Yang Jiayan, Guo Jingsong, et al. 2011. An asymmetric upwind flow, Yellow Sea Warm Current:1. New observations in the western Yellow Sea. Journal of Geophysical Research, 116:C04026
    Liu Qing. 2015. The interactions study between bloom-forming Ulva prolifera and phytoplankton in the Yellow Sea (in Chinese)[dissertation]. Qingdao:University of Chinese Academy of Sciences
    Liu Xin, Chiang K P, Liu Sumei, et al. 2015. Influence of the Yellow Sea Warm Current on phytoplankton community in the central Yellow Sea. Deep Sea Research Part I:Oceanographic Research Papers, 106:17-29, doi: 10.1016/j.dsr.2015.09.008
    Liu Dongyan, Keesing J K, Dong Zhijun, et al. 2010. Recurrence of the world's largest green-tide in 2009 in Yellow Sea, China:Porphyra yezoensis aquaculture rafts confirmed as nursery for macroalgal blooms. Marine Pollution Bulletin, 60(9):1423-1432, doi: 10.1016/j.marpolbul.2010.05.015
    Liu Dongyan, Keesing J K, He Peimin, et al. 2013. The world's largest macroalgal bloom in the Yellow Sea, China:formation and implications. Estuarine, Coastal and Shelf Science, 129:2-10, doi: 10.1016/j.ecss.2013.05.021
    Liu Dongyan, Keesing J K, Xing Qianguo, et al. 2009. World's largest macroalgal bloom caused by expansion of seaweed aquaculture in China. Marine Pollution Bulletin, 58(6):888-895, doi: 10.1016/j.marpolbul.2009.01.013
    Lü Xin'gang, Qiao Fangli. 2008. Distribution of sunken macroalgae against the background of tidal circulation in the coastal waters of Qingdao, China, in summer 2008. Geophysical Research Letters, 35:L23614, doi: 10.1029/2008GL036084
    Lü Xin'gang, Qiao Fangli, Xia Changshui, et al. 2010. Upwelling and surface cold patches in the Yellow Sea in summer:Effects of tidal mixing on the vertical circulation. Continental Shelf Research, 30(6):620-632, doi: 10.1016/j.csr.2009.09.002
    Lü Lian'gang, Wang Xiao, Wang Huiwu, et al. 2013. The variations of zooplankton biomass and their migration associated with the Yellow Sea Warm Current. Continental Shelf Research, 64:10-19, doi: 10.1016/j.csr.2013.05.007
    Parsons T R, Maita Y, Lalli C M. 1984. A Manual of Chemical and Biological Methods for Seawater Analysis. Oxford:Pergamon Press
    Qi Lin, Hu Chuanmin, Xing Qianguo, et al. 2016. Long-term trend of Ulva prolifera blooms in the western Yellow Sea. Harmful Algae, 58:35-44, doi: 10.1016/j.hal.2016.07.004
    Qi Jianhua, Su Yusong. 1998. Numerical simulation of the tide-induced continental front in the Yellow Sea. Oceanologia et Limnologia Sinica (in Chinese), 29(3):247-254
    Qiao Fangli, Ma Deyi, Zhu Mingyuan, et al. 2008. The green marcoalgal bloom in the Yellow Sea in 2008 and the scientific countmeasures. Advances in Marine Science (in Chinese), 26(3):409-410
    Qiao Fangli, Wang Guansuo, Lü Xin'gang, et al. 2011. Drift characteristics of green macroalgae in the Yellow Sea in 2008 and 2010. Chinese Science Bulletin, 56(21):2236-2242, doi: 10.1007/s11434-011-4551-7
    Ren Shihe, Xie Jiping, Zhu Jiang. 2014. The roles of different mechanisms related to the tide-induced fronts in the Yellow Sea in summer. Advances in Atmospheric Sciences, 31(5):1079-1089, doi: 10.1007/s00376-014-3236-y
    Su Yusong. 1986. A survey of geographical environment, circulation systems and the central fishing grounds in the Huanghai Sea and East China Sea. Journal of Shandong College of Oceanology (in Chinese), 16(1):12-27
    Sun Xiao, Wu Mengquan, Xing Qianguo, et al. 2018. Spatio-temporal patterns of Ulva prolifera blooms and the corresponding influence on chlorophyll-a concentration in the Southern Yellow Sea, China. Science of the Total Environment, 640-641:807-820
    Sun Yao, Yu Hong, Yang Qinfang, et al. 1990. Analysis and evaluation of nutritional condition and chemical indexes in Dingzi Bay waters. Journal of Fisheries of China (in Chinese), 14(1):33-39
    Tang Qisheng, Su Jilan. 2000. Study on Ecosystem Dynamics in Chinese Coastal Ocean I:Key Scientific Question and Study Stratagem (in Chinese). Beijing:Science Press
    Teague W J, Jacobs G A. 2000. Current observations on the development of the Yellow Sea Warm Current. Journal of Geophysical Research, 105(C2):3401-3411, doi: 10.1029/1999JC900301
    Wan Ruijing, Wei Hao, Sun Shan, et al. 2008. Spawning ecology of the anchovy Engraulis japonicus in the spawning ground of the Southern Shandong Peninsula I. Abundance and distribution characters of anchovy eggs and larvae. Acta Zoologica Sinica (in Chinese), 54(5):785-797
    Wang Yuntao, Castelao R M, Yuan Yeping. 2015. Seasonal variability of alongshore winds and sea surface temperature fronts in Eastern Boundary Current Systems. Journal of Geophysical Research:Oceans, 120(3):2385-2400, doi: 10.1002/2014JC010379
    Wang Rong, Gao Shangwu, Wang Ke, et al. 2003. Zooplankton indication of the Yellow Sea Warm Current in winter. Journal of Fisheries of China (in Chinese), 27(Suppl):39-48
    Wang Bin, Hirose N, Kang B, et al. 2014. Seasonal migration of the Yellow Sea Bottom Cold Water. Journal of Geophysical Research:Oceans, 119(7):4430-4443, doi: 10.1002/2014JC009873
    Wang Rong, Zuo Tao. 2004. The Yellow Sea Warm Current and the Yellow Sea Cold Bottom Water:their impact on the distribution of zooplankton in the southern Yellow Sea. Journal of the Korean Society of Oceanography, 39:1-13
    Wei Qinsheng. 2016. Characteristics and mechanisms of chemical hydrology and ecological responses in the southern Yellow Sea and off the Changjiang Estuary (in Chinese)[dissertation]. Qingdao:Ocean University of China
    Wei Qinsheng, Li Xiansen, Wang Baodong, et al. 2016a. Seasonally chemical hydrology and ecological responses in frontal zone of the central southern Yellow Sea. Journal of Sea Research, 112:1-12, doi: 10.1016/j.seares.2016.02.004
    Wei Qinsheng, Liu Lu, Zhan Run, et al. 2010a. Distribution features of the chemical parameters in the Southern Yellow Sea in summer. Periodical of Ocean University of China (in Chinese), 40(1):82-88
    Wei Qinsheng, Wang Baodong, Yao Qingzhen, et al. 2018. Hydro-biogeochemical processes and their implications for Ulva prolifera blooms and expansion in the world's largest green tide occurrence region (Yellow Sea, China). Science of the Total Environment, 645:257-266, doi: 10.1016/j.scitotenv.2018.07.067
    Wei Qinsheng, Yu Zhigang, Wang Baodong, et al. 2016b. Coupling of the spatial-temporal distributions of nutrients and physical conditions in the southern Yellow Sea. Journal of Marine Systems, 156:30-45, doi: 10.1016/j.jmarsys.2015.12.001
    Wei Qinsheng, Zhou Ming, Wei Xiuhua, et al. 2010b. Distribution features and variation tendency of chemical elements in the southern Yellow Sea in winter. Advances in Marine Science (in Chinese), 28(3):353-363
    Xia Zongwan, Guo Binghuo. 1983. The cold water and upwelling in the tip areas of Shandong peninsula and Liaodong peninsula. Journal of Oceanography of Huanghai & Bohai Seas (in Chinese), 1(1):13-19
    Xin Fuyan, Qu Keming, Cui Yi, et al. 2001. Distribution and variation of dissolved nutrient inorganic nitrogen and phosphorous in Aoshan Bay. Journal of Fishery Sciences of China (in Chinese), 8(4):79-82
    Xu Qing, Zhang Hongyuan, Ju Lian, et al. 2014. Interannual variability of Ulva prolifera blooms in the Yellow Sea. International Journal of Remote Sensing, 35(11-12):4099-4113
    Yu Fei, Zhang Zhixin, Diao Xinyuan, et al. 2006. Analysis of evolution of the Huanghai Sea Cold Water Mass and its relationship with adjacent water masses. Haiyang Xuebao (in Chinese), 28(5):26-34
    Yu Fei, Zhang Zhixin, Diao Xinyuan, et al. 2010. Observational evidence of the Yellow Sea Warm Current. Chinese Journal of Oceanology and Limnology, 28(3):677-683, doi: 10.1007/s00343-010-0006-2
    Yuan Dongliang, Li Yao, Wang Bin, et al. 2017a. Coastal circulation in the southwestern Yellow Sea in the summers of 2008 and 2009. Continental Shelf Research, 143:101-117, doi: 10.1016/j.csr.2017.01.022
    Yuan Yongquan, Yu Zhiming, Song Xiuxian, et al. 2017b. Temporal and spatial characteristics of harmful algal blooms in Qingdao Waters, China. Chinese Journal of Oceanology and Limnology, 35(2):400-414, doi: 10.1007/s00343-016-5279-7
    Yuan Dongliang, Zhu Jianrong, Li Chunyan, et al. 2008. Cross-shelf circulation in the Yellow and East China Seas indicated by MODIS satellite observations. Journal of Marine Systems, 70(1-2):134-149
    Zhang Shuwen, Wang Qingye, Lü Yan, et al. 2008. Observation of the seasonal evolution of the Yellow Sea Cold Water Mass in 1996-1998. Continental Shelf Research, 28(3):442-457, doi: 10.1016/j.csr.2007.10.002
    Zhao Baoren. 1985. The fronts of the Huanghai Sea Cold Water Mass induced by tidal mixing. Oceanologia et Limnologia Sinica (in Chinese), 16(6):451-460
    Zhao Baoren. 1987. A preliminary study of continental shelf fronts in the western part of southern Huanghai Sea and circulation structure in the front region of the Huanghai Cold Water Mass (HCWM). Oceanologia et Limnologia Sinica (in Chinese), 18(3):217-226
    Zhao Sheng, Yu Fei, Diao Xinyuan, et al. 2011. The path and mechanism of the Yellow Sea Warm Current. Marine Science (in Chinese), 35(11):73-80
    Zhou Mingjiang, Liu Dongyan, Anderson D M, et al. 2015. Introduction to the special issue on green tides in the Yellow Sea. Estuarine, Coastal and Shelf Science, 163:3-8, doi: 10.1016/j.ecss.2015.06.023
  • Relative Articles

  • Cited by

    Periodical cited type(3)

    1. Panfeng Li, Xuwen Qin, Yong Zhang, et al. Along-Strike Morphologic Variations of the Kyushu-Palau Ridge at 13°–17° N and Their Tectonic Implications. Russian Journal of Pacific Geology, 2024, 18(2): 220. doi:10.1134/S1819714024020052
    2. Xiaodong Wei, Weiwei Ding, Aiguo Ruan, et al. Crustal structure of the central and southern Kyushu-Palau Ridge: Implications for intra-oceanic arc evolution. Journal of Asian Earth Sciences, 2024, 267: 106151. doi:10.1016/j.jseaes.2024.106151
    3. Jie Zhang, Jiabiao Li, Weiwei Ding, et al. Crustal structure and magmatism of the southern Kyushu-Palau ridge. Tectonophysics, 2023, 858: 229862. doi:10.1016/j.tecto.2023.229862

    Other cited types(0)

  • 加载中

Catalog

    通讯作者: 陈斌, bchen63@163.com
    • 1. 

      沈阳化工大学材料科学与工程学院 沈阳 110142

    1. 本站搜索
    2. 百度学术搜索
    3. 万方数据库搜索
    4. CNKI搜索

    Article Metrics

    Article views (874) PDF downloads(346) Cited by(3)
    Proportional views
    Related

    /

    DownLoad:  Full-Size Img  PowerPoint
    Return
    Return