Feifei Zhang, Dingding Wang, Xiaolin Ji, Fanghui Hou, Yuan Yang, Wanyin Wang. Structural features in the mid-southern section of the Kyushu–Palau Ridge based on satellite altimetry gravity anomaly[J]. Acta Oceanologica Sinica, 2024, 43(4): 50-60. doi: 10.1007/s13131-024-2341-8
Citation: Feifei Zhang, Dingding Wang, Xiaolin Ji, Fanghui Hou, Yuan Yang, Wanyin Wang. Structural features in the mid-southern section of the Kyushu–Palau Ridge based on satellite altimetry gravity anomaly[J]. Acta Oceanologica Sinica, 2024, 43(4): 50-60. doi: 10.1007/s13131-024-2341-8

Structural features in the mid-southern section of the Kyushu–Palau Ridge based on satellite altimetry gravity anomaly

doi: 10.1007/s13131-024-2341-8
Funds:  ‘Research on Deep Structural Differences between Potential Oil-rich Depressions in Offshore basins of China Sea’ from the scientific and technological project of CNOOC Research Institute Co., Ltd., under contract No.CCL2021RCPS0167KQN; ‘Resource Potential, Accumulation Mechanism and Breakthrough Direction of Potential Oil-rich Depressions in Offshore China Sea’, under contract No. 220226220101; the Project of China Geological Survey under contract No. DD20191003; the National Natural Science Foundation of Shandong Province of China under contract No. ZR2022MD047.
More Information
  • The Kyushu–Palau Ridge (KPR), an anti-S-shaped submarine highland at the center of the Philippine Sea Plate (PSP), is considered the residual arc of the Izu–Bonin–Mariana Island Arc, which retains key information about the cessation of the Western Philippine Basin (WPB) expansion and the Parece Vela Basin (PVB) breakup. Herein, using the new generation of satellite altimetry gravity data, high-precision seafloor topography data, and newly acquired ship-borne gravity data, the topographic and gravity characteristics of the KPR mid-southern section and adjacent region are depicted. The distribution characteristics of the faults were delineated using the normalized vertical derivative–total horizontal derivative method(NVDR-THDR) and the minimum curvature potential field separation method. The Moho depth and crustal thickness were inverted using the rapid inversion method for a double-interface model with depth constraints. Based on these results, the crust structure features in the KPR mid-southern section, and the “triangular” structure geological significance where the KPR and Central Basin Rift (CBR) of the WPB intersect are interpreted. The KPR crustal thickness is approximately 6–16 km, with a distinct discontinuity that is slightly thicker than the normal oceanic crust. The KPR mid-southern section crust structure was divided into four segments (S1–S4) from north to south, formed by the CBR eastward extension joint action and clockwise rotation of the PVB expansion axis and the Mindanao fault zone blocking effect.
  • The stable isotope and elemental compositions of planktonic foraminiferal shells are important paleoceanographic proxies (Lea et al., 2000; Mohtadi et al., 2014; Stott et al., 2002). Based on the use of such proxies, a series of important research findings, including significant hydrological changes in the low latitude and the discovery of the super El Niño/Southern Oscillation (ENSO) or ENSO-like cycles have been obtained (Koutavas and Joanides, 2012; Mohtadi et al., 2010; Stott et al., 2002). Globigerinoides ruber (white) and Trilobatus sacculifer are two dominant species in tropical and subtropical oceans, serving as commonly used paleoceanographic proxies (Chowdhury et al., 2003; Jayan et al., 2021; Rippert et al., 2016; Sijinkumar et al., 2010; Steinke et al., 2005). In this respect, a major portion of existing knowledge on past surface seawater characteristics is based on these two surface-dwelling planktonic foraminifera species, G. ruber (white) and T. sacculifer (Gussone et al., 2004; Mohtadi et al., 2010; Steinke et al., 2005). However, the two species both have morphotypes, and the difference of geochemical proxies between these morphological shell types has not been thoroughly examined.

    Globigerinoides ruber (white) has two main morphotypes, namely, the sensu stricto (s.s.) and sensu lato (s.l., or “Elongatus morphotypes”) morphotypes, which in most tropical sediments, are usually differentiated using taxonomic criteria (Jayan et al., 2021; Wang, 2000). Specifically, G. ruber s.l. is characterized by a more compact and higher degree of trochospiral coiling than G. ruber s.s. (Wang, 2000; Figs 1a and b). Their different isotopic compositions of these morphotype shells from surface sediments indicate that they inhabit at different water depths, G. ruber s.s. at a shallower depth than that of G. ruber s.l. (Antonarakou et al., 2015; Kawahata, 2005; Steinke et al., 2005; Wang, 2000). Thus, in the last two decades, G. ruber Mg/Ca ratios have been widely used as a proxy for reconstructing sea surface temperature (SST) in tropical and subtropical oceans (Dang et al., 2012; Jia et al., 2018; Lea et al., 2000; Medina-Elizalde and Lea, 2005; Stott et al., 2002). However, some studies have indicated that the Mg/Ca ratios of G. ruber s.s. and G. ruber s.l. in the South China Sea and the eastern Indian Ocean are significantly different (Steinke et al., 2005, 2008). Conversely, the results of some other studies have revealed that the two G. ruber morphotypes have relatively similar Mg/Ca ratios in some areas in the eastern Indian Ocean (Mohtadi et al., 2009, 2010) and western equatorial Pacific Ocean (Tachikawa et al., 2014). It appears that it is still debatable whether or not difference of shell Mg/Ca ratios exists between different morphotypes, hence SST reconstruction. Thereby, whether there are significant differences between the Mg/Ca ratios of G. ruber s.s. and G. ruber s.l. might be important for reconstructing SST changes and upper-ocean structure.

    Figure  1.  Scanning electron micrographs of G. ruber (white) and T. sacculifer morphotypes: G. ruber sensu stricto (s.s.) (a), G. ruber sensu lato (s.l.) (b), T. sacculifer (with a sac-like final chamber) (c), and T. sacculifer (without a sac-like final chamber) (d). The foraminferal shells are from core-top sample of Core MD06-3047B. Scale bars (500 μm) are indicated on each image.

    Additionally, modern T. sacculifer also has two basic morphotypes, i.e., with (w) and without (w/o) a sac-like final chamber) (André et al., 2013; Bijma and Hemleben, 1994; Spezzaferri et al., 2015). Based on the analysis of sediment samples from the Atlantic Ocean, Elderfield et al. (2002) and Anand et al. (2003) observed that there are slight differences in oxygen-carbon isotope composition and Mg/Ca ratios between T. sacculifer (w) and T. sacculifer (w/o). Nevertheless, great care has been taken to distinguish the two morphotypes of T. sacculifer when this species is used for SST reconstruction in most previous studies (Lynch-Stieglitz et al., 2015; Setiawan et al., 2015). However, few studies have referred whether there is significant difference through geological time. This is important to verify if it is necessary to strictly select the same morphotype of T. sacculifer for downcore Mg/Ca analytical testing.

    However, the spatial and temporal differences in Mg/Ca between different morphotypes of G. ruber and T. saccluifer and its influencing factors are still not clear enough. How much the uncertainty of Mg/Ca records of different morphotypes exists is important for deciding whether mixed morphotypes of the two species can be used for SST reconstruction, especially in study areas where there are insufficient single morphotype shells for Mg/Ca testing. In addition, what environmental parameters are responsible for the Mg/Ca differences between different morphotypes can help figure out the potential of the inter-morphotype Mg/Ca differences as indicators for upper-water structure in the western tropical Pacific. The western Philippine Sea (WPS) has become an important area for paleoceanography studies of the western tropical Pacific. Exploring the differences in Mg/Ca ratios between the different morphotypes of G. ruber (white) and T. sacculifer in this area is significant for using their Mg/Ca ratios to reflect changes in SST or thermal structure. In this study, we investigate the differences in Mg/Ca ratios and SST estimations between the different morphotypes of G. ruber (white) and T. sacculifer in the WPS. On this basis, by collecting and comparing relevant data from other sites over the western tropical Pacific, we further discuss the spatial and temporal variations of inter-morphotype G. ruber Mg/Ca between different morphotypes and its potential as indicators of upper-ocean structure.

    Core MD06-3047B (17°00.44'N, 124°47.93'E; 2510 m water depth) from the WPS was recovered using a Calypso Giant Piston Corer on board the R/V Marion‒Dufresne vessel during the Marco Polo 2 cruise of the International Marine Global Change Studies Program (IMAGES) in June 2006. The coring site, located on the Benham Rise formed at ~50–45 Ma, is approximately 240 km east off Luzon (Hilde and Lee, 1984). The core sediments primarily consist of yellowish-brown silty clay with no visible evidence of sediment deformation (Jia et al., 2018), and for this study, the sediments were sampled at 4 cm intervals from the upper 60 cm of the core. The age model was established by correlating the benthic foraminiferal δ18O record of this core to LR04 stack (Jia et al., 2018). The 0–60 cm portion of the core provided a continuous record covering the last ~47 ka, spanning from MIS 3 to MIS 1, with a mean temporal resolution of approximately 3 ka (Jia et al., 2018).

    Approximately 30–50 tests of G. ruber s.s., G. ruber s.l., T. sacculifer (w), and T. sacculifer (w/o) (Fig. 1) respectively were hand-picked from the 250–300 μm size fraction to determine Mg/Ca ratios. Samples corresponding to the intervals of 38–39 cm and 42–43 cm were mixed for G. ruber s.s. Mg/Ca testing, as the amounts of G. ruber s.s. specimens in individual samples are insufficient for analysis. These samples for Mg/Ca analysis were cleaned in accordance with previously reported treatment procedures without a reductive step (Barker et al., 2003). Next, the Mg/Ca ratios were analyzed by inductively coupled plasma-optical emission spectrometry (ICP-OES, iCAP6300 radial, Thermo Fisher, Waltham, MA, USA) at the Institute of Oceanology, Chinese Academy of Sciences. The long-term internal reproducibility (SD) is 0.07 mmol/mol based on a consistency standard (Mg/Ca = 3.33 mmol/mol). Additionally, the elemental ratios of Fe/Ca, Al/Ca, and Mn/Ca were analyzed at the same time to monitor any contamination by the oxides of iron and manganese, which could bias the Mg/Ca measurements. Notably, there was no significant positive correlation of Fe/Ca, Al/Ca, and Mg/Ca values, indicating negligible contamination of detrital materials and/or metal-oxides on the Mg/Ca ratios (Fig. 2).

    Figure  2.  Fe/Ca (left), Al/Ca (middle) and Mn/Ca (right) versus Mg/Ca of the G. ruber (white) and T. sacculifer. There is no obvious correlation between these records, indicating negligible influence of contamination on Mg/Ca ratios.

    First, the Mg/Ca (mmol/mol) ratios were converted to temperature (T, ℃) using a multispecies regression equation (Anand et al., 2003):

    $$ T = {\mathrm{ln}} [{\mathrm{Mg}}/{\mathrm{Ca}} /0.38]/0.09. $$ (1)

    The error of the temperature reconstruction was estimated by propagating the errors introduced by the Mg/Ca measurements and the Mg/Ca-temperature calibration, as introduced by Mohtadi et al. (2014). The resulting errors were on average of ± 1.0℃, similar to those associated with Mg/Ca-temperature calibration reported in other studies (Anand et al., 2003; Dekens et al., 2002; Sagawa et al., 2012).

    For comparison, we also use species-specific equations. including

    $$ T = {\mathrm{ln}}[{\mathrm{Mg}}/{\mathrm{Ca}}/0.31]/0.097 \;{\mathrm{for}}\;G.\;ruber\; {\mathrm{s}}.{\mathrm{s}} .$$ (2)
    $$ T ={\mathrm{ln}} [{\mathrm{Mg}}/{\mathrm{Ca}} /0.32]/0.097 \;{\mathrm{for}}\; G.\; ruber\; {\mathrm{s}}.{\mathrm{l}}. $$ (3)
    $$ T= {\mathrm{ln}} [{\mathrm{Mg}}/{\mathrm{Ca}} /1.06]/0.048 \;{\mathrm{for}} \;T. \;sacculifer\; ({\mathrm{w}}/{\mathrm{o}}). $$ (4)
    $$ T = {\mathrm{ln}} [{\mathrm{Mg}}/{\mathrm{Ca}} /0.67]/0.069 \;{\mathrm{for}} \;T.\; sacculifer \;({\mathrm{w}}.).$$ (5)

    from Anand et al. (2003). The errors estimated by error propagation are on average of ±2.1℃ and ±2.0℃ for temperature estimates of G. ruber s.s. and G. ruber s.l., ±7.4℃ and ±8.0℃ for that of T. sacculifer (w/o) and T. sacculifer (w), respectively. By comparison of temperature estimates and corresponding differences based on multispecies and species-specific equations, consistent variation trends were observed (Fig. S1). Considering the larger errors from temperature estimate based on species-specific equations, and to facilitate the comparison of temperature estimates from Mg/Ca ratios of different morphotypes of G. ruber (white) and T. sacculifer, we choose to use the estimates using the multispecies equation in the following analysis.

    To compare results from our study site with those from other regions across the western tropical Pacific, we collect published Mg/Ca ratios of G. ruber s.s. and G. ruber s.l. of the core-top sediment samples from the South China Sea, Offshore Papua New Guinea, Gulf of Papua, and the WPS (Hollstein et al., 2017; Regoli et al., 2015; Steinke et al., 2005). Then the difference between the Mg/Ca ratios of two morphotypes of G. ruber [Δ(Mg/Ca)G.ruber s.s.−s.l.] are estimated see Table S2. Data of mixed-layer depth (MLD) of annual average, average from June to August, and average from December to February, respectively, are retrieved from Monterey and Levitus (1997).

    Welch’s test is performed to analyze the Mg/Ca ratios of the different morphotypes of G. ruber and T. sacculifer. First, for comparison, the mean values of the two sets of data are assumed to be equivalent if the Mg/Ca ratios of G. ruber s.s. and G. ruber s.l., or T. sacculifer (w.) and T. sacculifer (w/o) are similar. Thus, the null hypothesis is accepted (H = H0), indicating no significant difference in Mg/Ca ratios between the different morphotypes of G. ruber and T. sacculifer. If significant mean differences were observed, the null hypothesis would be rejected (H = Ha), indicating significant difference in Mg/Ca ratios between different morphological types of these species.

    The Mg/Ca ratios and temperature estimates of G. ruber s.s., G. ruber s.l., T. sacculifer (w), and T. sacculifer (w/o) generally show similar variation trends since MIS 3 (Fig. 3; Table S1). The core-top Mg/Ca values of G. ruber s.s., G. ruber s.l., T. sacculifer (w), and T. sacculifer (w/o) yields temperature estimates of 27.2℃, 25.6℃, 24. 7℃, 24.5℃, respectively, matching well with previous conclusions that the calcification depth of G. ruber s.s. is shallower than that of G. ruber s.l., and that the calcification depth of G. ruber is shallower than that of T. sacculifer (Hollstein et al., 2017). The overall Mg/Ca ratios of G. ruber s.s. was higher than that of G. ruber s.l., and the Δ(Mg/Ca)G.ruber s.s.−s.l. varied in the range of ~0.05–0.59 mmol/mol, with a mean value of ~0.32 mmol/mol. The differences in temperature estimates based on Mg/Ca ratios of G. ruber s.s. and G. ruber s.l. varied between ~0.2℃ and 1.6℃, with a mean value of ~1.0℃. The difference between the Mg/Ca ratios of the two morphotypes of T. sacculifer [Δ(Mg/Ca)T.sacculifer w/o−w] varies in the range of –0.29 mmol/mol to 0.35 mmol/mol, with a mean value ~0.03 mmol/mol, and an average value of absolute Δ(Mg/Ca)T.sacculifer w/o−w of ~0.13 mmol/mol. The corresponding difference in temperature estimates varied in the range of –1.1–1.3℃, showing a mean difference of ~0.1℃. The variation in Δ(Mg/C)aT.sacculifer w/o−w could be divided into three stages: (1) since ~12.5 ka, Mg/Ca ratios of T. sacculifer (w) and T. sacculifer (w/o) show no obvious differences (the mean difference was only ~0.08 mmol/mol), (2) between ~21.7 ka and 12.5 ka, Mg/Ca of T. sacculifer (w/o) are slightly higher than those of T. sacculifer (w), with a mean difference of ~0.16 mmol/mol, (3) between ~40 ka and 21.7 ka, Mg/Ca ratios of T. sacculifer (w/o) are lower than those of T. sacculifer (w), with a mean difference of ~0.18 mmol/mol. Additionally, during the 47–40 ka period, Mg/Ca of T. sacculifer (w) and T. sacculifer (w/o) show no obvious differences, with a mean difference of ~0.07 mmol/mol.

    Figure  3.  Mg/Ca ratios (a, b), and Mg/Ca-based temperature (c, d) estimates (based on the multispecies equation from Anand et al. (2003)) of G. ruber s.s., G. ruber s.l., T. sacculifer (w), and T. sacculifer (w/o) of Core MD06-3047B , and the differences between these records for each species.

    By comparing the Δ(Mg/Ca)G.ruber s.s.−s.l. records obtained from our study and those from core-top samples in other areas across the western tropical Pacific, we found that the Δ(Mg/Ca)G.ruber s.s.−s.l. vary widely from one area to another (Antonarakou et al., 2015; Hollstein et al., 2017; Mohtadi et al., 2010). Based on the collected data, Δ(Mg/Ca)G.ruber s.s.−s.l. offshore Papua New Guinea and in the Gulf of Papua are relatively small, with an average difference of approximately 0.09 mmol/mol. The average Δ(Mg/Ca)G.ruber s.s.−s.l. from the WPS and the South China Sea are larger, with an average of ~0.47 mmol/mol. Generally, Δ(Mg/Ca)G.ruber s.s.−s.l. values from core-top samples collected around the center of western Pacific warm pool (WPWP) are relatively small, while those observed from the northern margin of the WPWP were relatively large (Fig. 4).

    Figure  4.  Comparison of Δ(Mg/Ca)G.ruber s.s.−s.l. from core-top sediment samples and the mixed-layer depth (MLD) in the western tropical Pacific. a. Location map showing the sites of the core-top Δ(Mg/Ca)G.ruber s.s.−s.l.. The red dots represent positive Δ(Mg/Ca)G.ruber s.s.−s.l. values, while the white dots indicate negative Δ(Mg/Ca)G.ruber s.s.−s.l. values. The sizes of the symbols denote the magnitude of the absolute values of the Δ(Mg/Ca)G.ruber s.s.−s.l. records. b–e. Plots of Δ(Mg/Ca)G.ruber s.s.−s.l. vs. average annual MLD (b), difference between winter and summer MLD (c), summer MLD (d), and winter MLD (e). Data of the South China Sea and WPS are obtained from Steinke et al. (2005), those in the Gulf of Papua and offshore Papua New Guinea from Hollstein et al. (2017) and Tachikawa et al. (2014). Detailed information regarding the data sites is provided in Table S2. Data of the MLD are from Monterey and Levitus (1997).

    The average value of Δ(Mg/Ca)G.ruber s.s.−s.l. is approximately ~0.32 mmol/mol, with a maximum value of ~0.59 mmol/mol, much larger than the reproducibility of Mg/Ca measurements (±0.07 mmol/mol). The difference between the Mg/Ca-based temperature estimates of these two morphotypes is ~1.0℃ on average, with a maximum difference of 1.6℃, also larger than the error of the temperature reconstructions (±1℃). Based on the results of the Welch’s test performed on the Mg/Ca ratios of G. ruber s.s. and G. ruber s.l. and the estimated water temperatures, we reject the null hypothesis (H = Ha), thus, the G. ruber Mg/Ca ratios and estimated temperatures of different morphotypes show significant difference in the WPS (Table 1). Therefore, when using the Mg/Ca ratios G. ruber for SST reconstruction, a monomorphic type of is better to be selected, to avoid the effect of the Mg/Ca differences between morphotypes on the reconstruction of SSTs. If mixed morphotypes must be used, the difference value (~1.0℃) between different morphotypes should be considered.

    Table  1.  Mean values of the Mg/Ca ratios (mmol·mol−1) of G. ruber s.s., G. ruber s.l., T. sacculifer (w), and T. sacculifer (w/o) based on geochemical measurements using Core MD06-3047B with outcomes obtained by performing Welch’s t test at a p < 0.05 significance level.
    Mg/Ca T/℃ Mg/Ca T/℃
    G. ruber s.s. 3.72 25.29 T. sacculifer (w/o) 3.17 23.64
    G. ruber s.l. 3.40 24.28 T. sacculifer (w) 3.20 23.53
    H Ha Ha H H0 H0
    Note: H = H0 implies the null hypothesis cannot be rejected; H = Ha implies the null hypothesis can be rejected.
     | Show Table
    DownLoad: CSV

    Additionally, the Mg/Ca ratios of T. sacculifer (w) and T. sacculifer (w/o) at MD06-3047B and the estimated water temperatures showed relatively small differences. The Δ(Mg/Ca)T.sacculifer w/o−w value is only ~0.03 mmol/mol on average, less than the reproducibility of Mg/Ca measurements (±0.07 mmol/mol). Further, Welch’s test also showed that the Mg/Ca ratios of T. sacculifer (w) and T. sacculifer (w/o) and the estimated temperatures are not significantly different (H = H0) (Table 1). However, the Δ(Mg/Ca)T.sacculifer w/o−w show staged variation, and the Δ(Mg/Ca)T.sacculifer w/o−w values are also greater than the reproducibility of the Mg/Ca measurements (±0.07 mmol/mol) within different periods. This possibly indicates that there are significant differences in the Mg/Ca ratios of T. sacculifer (w) and T. sacculifer (w/o) during different geological periods. Thus, when reconstructing water temperature using the Mg/Ca ratios of T. sacculifer in the WPS, it is better to use data based on a monomorphic type of T. sacculifer if the sedimentary records span over glacial-interglacial stages.

    A question is, “what causes the staged variation of Δ(Mg/Ca)T.sacculifer w/o-w.?”. In a previous study, it was observed that the abundance of the T. sacculifer (w) morphotype increases with increasing temperature compared with that of T. sacculifer (w/o) (André et al., 2013), possibly indicating a higher calcification temperature of T. sacculifer (w) than that of T. sacculifer (w/o). However, this feature cannot explain the step-wise negative/positive variations in the Δ(Mg/Ca)T.sacculifer w/o−w of MD06-3047B. Modern research indicated that there are little or no genetic differences between T. sacculifer morphotypes, and the development of the sac-like final chamber is just considered as a phenotype (André et al., 2013; Darling and Wade, 2008). However, based on fossil records, the two morphotypes have different ecological preferences (Spezzaferri et al., 2015). Thus, it seemed more complicated when the habitats of the different T. sacculifer morphotypes through geological time were considered. According to δ18O and Mg/Ca records of the two morphotypes from Anand et al. (2003), the calcification depth of T. sacculifer (w/o) seems shallower (deeper) than T. sacculifer (w) during winter (summer), which may be due to the different sensitivity of the two morphotypes to seasonal variation. Thus, we speculate that the staged negative/positive variations of ΔTT.sacculifer w/o−w since ~47 ka might be associated with alternate stronger winter/summer monsoons (Fig. 5; Cheng et al., 2016; Guo et al., 2009). That is higher (lower) ΔTT.sacculifer w/o−w corresponds to stronger (weaker) East Asian winter monsoon (EAWM) and weaker (stronger) East Asian summer monsoon (EASM) (Figs 5a, c and d). This seems conflict with the modern observation that deeper mixed layer depth and smaller vertical thermal gradient are developed during winter than summer around our study area (Locarnini et al., 2013). Thus, the relationship between ΔTT.sacculifer w/o−w and upper-ocean structure might be more complicated, that needs to be fully examined by further studies on core-top sediment, sediment trap and/or culture experiments.

    Figure  5.  Estimated water temperatures based on the Mg/Ca ratios of G. ruber s.s., G. ruber s.l., T. sacculifer (w), and T. sacculifer (w/o) (a). Differences in estimated temperatures based on the Mg/Ca ratios of G. ruber s.s. and G. ruber s.l. (b). Smectite crystallinity values of Core MD06-3047 (Xu et al., 2012) and the speleothem δ18O records from caves in China (Cheng et al., 2016) (c). Differences in estimated temperatures based on the Mg/Ca ratios of T. sacculifer (w/o) and T. sacculifer (w), and loess grain-size changes (content of the >32 μm fraction) from Xifeng Section, Loess Plateau of China (Guo et al., 2009). EASM and EAWM means East Asian summer and winter monsoons, respectively."+" means stronger monsoon, "-" means weaker monsoon.

    Considering that G. ruber predominantly inhabits in the mixed layer, we assumed that the regional distribution of core-top Δ(Mg/Ca)G.ruber s.s.−s.l. could be attributed to variability in the MLD of the upper-ocean. Accordingly, the core-top Δ(Mg/Ca)G.ruber s.s.−s.l. from the western tropical Pacific are compared with the annual mean MLD and the seasonal difference of MLD between winter and summer at these sites (Fig. 4b and c). As shown in Fig. 4b, no correlation is observed between the core-top Δ(Mg/Ca)G.ruber s.s.−s.l. and the annual mean MLD (R2 = 0.00004), indicating little influence of annual mean MLD on the Δ(Mg/Ca)G.ruber s.s.−s.l.. We also observed a negative relationship between the core-top Δ(Mg/Ca)G.ruber s.s.−s.l. and the difference between MLD in winter and summer (R2 = 0.269), possibly suggesting a certain role played by the seasonal fluctuation of MLD on the Δ(Mg/Ca)G.ruber s.s.−s.l. values. For further analysis, we plotted the core-top Δ(Mg/Ca)G.ruber s.s.−s.l. against the winter and summer MLD (Figs 4d and e), from which it can be observed that the Δ(Mg/Ca)G.ruber s.s.−s.l. are more significantly correlated with summer MLD (R2 = 0.523) than winter MLD (R2 = 0.001). This observation implies that changes in summer MLD may be the major factor for the Δ(Mg/Ca)G.ruber s.s.−s.l. over the western tropical Pacific.

    The negative relationship between the core-top Δ(Mg/Ca)G.ruber s.s.−s.l. and summer MLD indicates that larger Δ(Mg/Ca)G.ruber s.s.−s.l. might be related to a shallower summer MLD. Although there is evidence that G. ruber can generally reflect annual surface water conditions (Dekens et al., 2002; Lea et al., 2000), a plankton tow study indicates high abundance of G. ruber (white) in the warm summer surface waters of the western tropical Pacific (Troelstra and Kroon, 1989). Additionally, the two morphotypes of G. ruber (G. ruber s.s. and G. ruber s.l.) show different vertical distributions (Kuroyanagi and Kawahata, 2004), as G. ruber s.s. is consistently calcified in warmer waters to a greater degree than G. ruber s.l., due to the shallower habitat depth or a more summer-weighted seasonal distribution of G. ruber s.s (Antonarakou et al., 2015). It has also been reported that in the South China Sea, G. ruber s.s. is predominant in summer and in shallower water layers (0–40 m), while G. ruber s.l. predominant in deeper water layers (40–60 m) (Zhang et al., 2020). Thus, we could speculate that the summer-weighted seasonal distribution and different vertical distributions of G. ruber morphotypes can lead to relatively strong negative correlation between Δ(Mg/Ca)G.ruber s.s.−s.l. and summer MLD in the western tropical Pacific. This relationship lends strong support to propose that Δ(Mg/Ca)G.ruber s.s.−s.l. and the corresponding differences in Mg/Ca-based temperatures between G. ruber s.s. and G. ruber s.l. (ΔTG.ruber s.s.−s.l.) are indicators of changes in summer mixed-layer depth.

    To verify whether Δ(Mg/Ca)G.ruber s.s.−s.l. and ΔTG.ruber s.s.−s.l. records can reflect summer MLD in geological history, we compared the ΔTG.ruber s.s.−s.l. records from Core MD06-3047B since 47 ka with records reflecting changes in EASM, including the smectite crystallinity values from Site MD06-3047 (Xu et al., 2012) and the speleothem δ18O records from caves in China (Cheng et al., 2016). As shown in Fig. 5, the ΔTG.ruber s.s.−s.l. records show overall covariance with EASM, with higher values corresponding to a stronger EASM. In the WPS, a stronger EASM can lead to a shallower summer MLD (Hao et al., 2012; Qu, 2003; Tozuka et al., 2002), as in the boreal summer, the prevailing westly/southwesterly wind may produce a positive wind stress curl that results in Ekman divergence and upwelling in the study area (Xiong et al., 2022). Thus, in the WPS, the significant difference between estimated temperatures based on downcore Mg/Ca ratios of G. ruber s.s. and G. ruber s.l. might be affected by EASM-driven summer MLD changes. This result provides further evidence for the potential of ΔTG.ruber s.s.−s.l. records to serve as indicators of changes in summer MLD at orbital timescale. The higher Δ(Mg/Ca)G.ruber s.s.−s.l. indicates shallower MLD and stronger EASM. More records of the MLD and summer monsoon in this area are needed for further validation.

    In this study, we present the Mg/Ca records of G. ruber s.s., G. ruber s.l., T. sacculifer (w), and T. sacculifer (w/o) in the upper 60 cm of Core MD06-3047B obtained from the WPS. Based on the results obtained, we arrived at the following conclusions:

    (1) Since MIS 3, significant consistent differences exist in the Mg/Ca ratios and the estimated temperatures between the morphotypes of G. ruber from Site MD06-3047B; while staged changes in the differences of Mg/Ca ratios are observed between the morphotypes of T. sacculifer. We suggest that when reconstructing SST changes using the Mg/Ca ratios of G. ruber and T. sacculifer in the WPS, it should be better to select one single morphotype.

    (2) Systematic patterns are observed in the core-top Δ(Mg/Ca)G.ruber s.s.−s.l. across the western tropical Pacific, which show negative correlation with summer MLD. Combining with the covariance between downcore ΔTG.ruber s.s.−s.l. records of Core MD06-3047B and summer monsoon changes, we propose that Δ(Mg/Ca)G.ruber s.s.−s.l. may serve as a potential indicator of summer MLD changes.

    Acknowledgements: We thank the chief scientist Carlo Laj, the captain and the crew of R/V Marion Dufresne, and the French Polar Institute for their efforts and support during the IMAGES XIV-MD155-Marco Polo 2 cruise.
  • Chao Dingbo, Yao Yunsheng, Li Jiancheng, et al. 2002. Interpretaion on the tectonics and characteristics of altimeter-derived gravity anomalies in China South Sea. Geomatics and Information Science of Wuhan University (in Chinese), 27(4): 343–347
    Calvert A J. 2011. The seismic structure of island arc crust. In: Brown D, Ryan P D, eds. Arc-Continent Collision Berlin, Heidelberg: Springer, 87–119
    Deschamps A, Lallemand S. 2002. The West Philippine Basin: an Eocene to Early Oligocene back arc basin opened between two opposed subduction zones. Journal of Geophysical Research: Solid Earth, 107(B12): 2322
    Deschamps A, Okino K, Fujioka K. 2002. Late amagmatic extension along the central and eastern segments of the West Philippine Basin fossil spreading axis. Earth and Planetary Science Letters, 203(1): 277–293, doi: 10.1016/S0012-821X(02)00855-5
    Graw J H, Wood W T, Phrampus B J. 2021. Predicting global marine sediment density using the random forest regressor machine learning algorithm. Journal of Geophysical Research: Solid Earth, 126(1): e2020JB020135, doi: 10.1029/2020JB020135
    Hall R, Ali J R, Anderson C D, et al. 1995. Origin and motion history of the Philippine Sea Plate. Tectonophysics, 251(1–4): 229–250, doi: 10.1016/0040-1951(95)00038-0
    Hilde T W C, Lee C S. 1984. Origin and evolution of the west Philippine Basin: a new interpretation. Tectonophysics, 102(1–4): 85–104, doi: 10.1016/0040-1951(84)90009-X
    Hou Fanghui, Qin Ke, Lu Kai, et al. 2022. Tectono-sedimentary characteristics and subduction initiation in the middle Kyushu-Palau Ridge and adjacent basins: a comprehensive study of multichannel seismic reflection profiles. Marine Geology & Quaternary Geology (in Chinese), 42(5): 187–198
    Hu Zhengan, Zhang Chunguan, Ren Zusong, et al. 2014. Calculation of the Satellite Altimetric bouguer gravity anomaly based on regression analysis. Journal of Yangtze University (Natural Science Edition) (in Chinese), 11(1): 21–23
    Ishihara T, Koda K. 2007. Variation of crustal thickness in the Philippine Sea deduced from three-dimensional gravity modeling. Island Arc, 16(3): 322–337, doi: 10.1111/j.1440-1738.2007.00593.x
    Ishizuka O, Taylor R N, Yuasa M, et al. 2011. Making and breaking an island arc: a new perspective from the Oligocene Kyushu-Palau arc, Philippine Sea. Geochemistry, Geophysics, Geosystems, 12(5): Q05005
    Ji Xiaolin, Wang Wanyin, Du Xiangdong, et al. 2019a. Tectonic division by gravity and magnetic anomaly data of salt-bearing basins, south-central section of West Africa. Chinese Journal of Geophysics (in Chinese), 62(4): 1502–1514
    Ji Xiaolin, Wang Wanyin, Qiu Zhiyun. 2015. The research to the minimum curvature technique for potential field data separation. Chinese Journal of Geophysics (in Chinese), 58(3): 1042–1058
    Ji Xiaolin, Wang Wanyin, Qiu Zhiyun. 2019b. Parameter choose experimental research to the minimum curvature technique potential field data separation method. Progress in Geophysics (in Chinese), 34(4): 1441–1452
    Jin Xianglong, Gao Jinyao. 2001. The satellite altimetry gravity field and the geodynamic feature in the West Pacific. Marine Geology & Quaternary Geology (in Chinese), 21(1): 1–6
    Kobayashi K. 2004. Origin of the Palau and Yap trench-arc systems. Geophysical Journal International, 157(3): 1303–1315, doi: 10.1111/j.1365-246X.2003.02244.x
    Lallemand S, Arcay D. 2021. Subduction initiation from the earliest stages to self-sustained subduction: insights from the analysis of 70 Cenozoic sites. Earth-Science Reviews, 221: 103779, doi: 10.1016/j.earscirev.2021.103779
    Lei Shoumin. 1984. Calculation of generalized topographic and isostatic gravity corrections. Marine Geology & Quaternary Geology (in Chinese), 4(1): 101–111
    Li Changzhen, Li Naisheng, Lin Meihua. 2000. Terrain features of the Philippine sea. Marine Sciences (in Chinese), 24(6): 47–51
    Luo Zhicai, Zhong Bo, Zhou Hao, et al. 2022. Progress in determining the Earth’s gravity field model by satellite gravimetry. Geomatics and Information Science of Wuhan University (in Chinese), 47(10): 1713–1727
    Ning Jinsheng, Wang Zhengtao, Chao Nengfang. 2016. Research status and progress in international next-generation satellite gravity measurement missions. Geomatics and Information Science of Wuhan University (in Chinese), 41(1): 1–8
    Nishizawa A, Kaneda K, Katagiri Y, et al. 2007. Variation in crustal structure along the Kyushu-Palau Ridge at 15–21°N on the Philippine Sea plate based on seismic refraction profiles. Earth, Planets and Space, 59(6): e17–e20
    Nishizawa A, Kaneda K, Oikawa M. 2016. Crust and uppermost mantle structure of the Kyushu-Palau Ridge, remnant arc on the Philippine Sea Plate. Earth, Planets and Space, 68(1): 30
    Okino K, Fujioka K. 2003. The Central Basin Spreading Center in the Philippine Sea: structure of an extinct spreading center and implications for marginal basin formation. Journal of Geophysical Research: Solid Earth, 108(B1): 2040
    Park J O, Hori T, Kaneda Y. 2009. Seismotectonic implications of the Kyushu - Palau Ridge subducting beneath the westernmost Nankai forearc. Earth, Planets and Space, 61(8): 1013–1018
    Pasyanos M E, Masters T G, Laske G, et al. 2014. LITHO1.0: an updated crust and lithospheric model of the Earth. Journal of Geophysical Research: Solid Earth, 119(3): 2153–2173, doi: 10.1002/2013JB010626
    Sandwell D, Garcia E, Soofi K, et al. 2013. Toward 1-mGal accuracy in global marine gravity from CryoSat-2, Envisat, and Jason-1. The Leading Edge, 32(8): 892–899, doi: 10.1190/tle32080892.1
    Sandwell D T, Harper H, Tozer B, et al. 2021. Gravity field recovery from geodetic altimeter missions. Advances in Space Research, 68(2): 1059–1072, doi: 10.1016/j.asr.2019.09.011
    Sandwell D T, Müller R D, Smith W H F, et al. 2014. New global marine gravity model from CryoSat-2 and Jason-1 reveals buried tectonic structure. Science, 346(6205): 65–67, doi: 10.1126/science.1258213
    Shang Luning, Li Panfeng, Du Runlin, et al. 2021. Structural characteristics of the KPR-CBR triple-junction inferred from gravity and magnetic interpretations, Philippine Sea Plate. China Geology, 4(4): 541–552, doi: 10.31035/cg2021089
    Smith W H F, Sandwell D T. 1997. Global sea floor topography from satellite altimetry and ship depth soundings. Science, 277(5334): 1956–1962, doi: 10.1126/science.277.5334.1956
    Song Weiyu, Liu Yanan, Hu Bangqi, et al. 2021. Landform classification for the Philippine Sea Based on DEM data. Marine Geology & Quaternary Geology (in Chinese), 41(1): 192–198
    Stern R J. 2004. Subduction initiation: spontaneous and induced. Earth and Planetary Science Letters, 226(3–4): 275–292, doi: 10.1016/S0012-821X(04)00498-4
    Stern R J, Gerya T. 2018. Subduction initiation in nature and models: a review. Tectonophysics, 746: 173–198, doi: 10.1016/j.tecto.2017.10.014
    Straume E O, Gaina C, Medvedev S, et al. 2019. GlobSed: updated total sediment thickness in the world’s oceans. Geochemistry, Geophysics, Geosystems, 20(4): 1756–1772
    Taylor B, Goodliffe A M. 2004. The West Philippine Basin and the initiation of subduction, revisited. Geophysical Research Letters, 31(12): L12602
    Wang Wanyin, Pan Zuoshu. 1993. Fast solution of forward and inverse problems for gravity field in a dual interface model. Geophysical Prospecting for Petroleum (in Chinese), 32(2): 81–87, 123
    Wang Wanyin, Pan Yu, Qiu Zhiyun. 2009. A new edge recognition technology based on the normalized vertical derivative of the total horizontal derivative for potential field data. Applied Geophysics, 6(3): 226–233, doi: 10.1007/s11770-009-0026-x
    Wang Dingding, Wang Wanyin, Zhu Yingjie, et al. 2021. Extraction methods and application of feature points of edge recognition for potential field. Chinese Journal of Geophysics (in Chinese), 64(4): 1401–1411
    Wei Xiaodong, Ding Weiwei, Ruan Aiguo, et al. 2022. Crustal structure and variation along the southern part of the Kyushu-Palau Ridge. Acta Oceanologica Sinica, 41(1): 50–57, doi: 10.1007/s13131-021-1979-8
    Wu Shiguo, Fan Jianke, Dong Dongdong. 2013. Discussion on the tectonic division of the Philippine Sea Plate. Chinese Journal of Geology (in Chinese), 48(3): 677–692
    Wu J, Suppe J, Lu Renqi, et al. 2016. Philippine Sea and East Asian plate tectonics since 52 Ma constrained by new subducted slab reconstruction methods. Journal of Geophysical Research: Solid Earth, 121(6): 4670–4741, doi: 10.1002/2016JB012923
    Yamashita M, Tsuru T, Takahashi N, et al. 2007. Fault configuration produced by initial arc rifting in the Parece Vela Basin as deduced from seismic reflection data. Island Arc, 16(3): 338–347, doi: 10.1111/j.1440-1738.2007.00594.x
    Yang Huiliang, Wei Jia, Li Panfeng, et al. 2021. Characteristics of the stratigraphic architectures of the shallow sections in deep sea basin on both sides of Kyushu - Palau ridge. Marine Geology & Quaternary Geology (in Chinese), 41(1): 14–21
    Zhang Jie, Li Jiabiao, Ding Weiwei. 2012. Reviews of the study on crustal structure and evolution of the Kyushu-Palau ridge. Advances in Marine Science (in Chinese), 30(4): 595–607
    Zhang Feifei, Wang Hao, Zhang Yimi, et al. 2023. Accuracy analysis of satellite altimetry gravity data in the Western Pacific Area. Geomatics and Information Science of Wuhan University (in Chinese), ). https://doi.org/10.13203/j.whugis20220429 [2023/-05/-26]
    Zhang Minghua, Zhang Jiaqiang. 2005. Resolution of modern satellite altimetric gravity anomaly and its application to marine geological survey. Geophysical & Geochemical Exploration (in Chinese), 29(4): 295–298,303
    Zhang Chunguan, Zhang Minghua, An Yulin. 2007. The negative correlation correction between gravity anomaly and topography of Tibetan Plateau. Geophysical & Geochemical Exploration (in Chinese), 31(3): 236–238
  • Relative Articles

  • 加载中

Catalog

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

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

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

    Figures(11)  / Tables(1)

    Article Metrics

    Article views (296) PDF downloads(31) Cited by()
    Proportional views
    Related

    /

    DownLoad:  Full-Size Img  PowerPoint
    Return
    Return