| Citation: | Jun Dai, Huizan Wang, Weimin Zhang, Pinqiang Wang, Tengling Luo. Three-dimensional structure of an observed cyclonic mesoscale eddy in the Northwest Pacific and its assimilation experiment[J]. Acta Oceanologica Sinica, 2021, 40(5): 1-19. doi: 10.1007/s13131-021-1810-6
|
Magmatism at mid-ocean ridges is one of our planet’s most important geological processes as it forms the oceanic crust, which covers nearly two-thirds of the Earth’s surface. Unlike fast spreading ridges, it is widely accepted that ultraslow–spreading ridges have relatively low extents of melting and magmatism, thinner crust, and experience limited crustal processes (Dick, 1989). Although ultraslow–spreading ridges are characterized by low magma supplies, there is also evidence of substantial magmatic processes in magma reservoirs, rather than the stable magma chambers present under some segments (Dick et al., 2003; Coogan et al., 2001; Jian et al., 2017). Until now, how magmas behave in the crust and how the thermal and dynamic regimes of the magma reservoirs and conduit systems effect the eruption mechanism of magmas from ultraslow mid-ocean ridges (MORs) has received comparatively little attention. As the product of erupted magma, mid-ocean ridge basalts (MORBs) can provide information about the magmatism under ridges (Yang et al., 2013). Therefore, unraveling the chemical effects of magmatic processes reflected by the MORBs is key to determining the details of the magma plumbing system beneath ultraslow oceanic spreading centers.
Previous studies have mainly focused on the whole–rock geochemical and isotopic compositions of MORBs from ultraslow–spreading ridges (Dick et al., 2003; Yang et al., 2013, 2017). However, whole–rock geochemistry cannot provide information about the magmatic processes occurring within the crust, which we are attempting to understand, because the bulk–rock compositions of MORBs reflect the variable overprinting of the low–pressure differentiation of mantle–derived primary liquids (Stolper, 1980; Yang et al., 2013). Unlike whole–rock compositions, the major and trace element contents of plagioclase phenocrysts are easily preserved and are sensitive to the physicochemical conditions of the melt from which they crystallized (Ginibre et al., 2004; Mollo et al., 2011; Mutch et al., 2019; Bennett et al., 2019). Their geochemical compositions and morphology are commonly used to obtain information about the conditions of the magmatic system from which they crystallized. For example, Mutch et al. (2019) established an element diffusion model for plagioclase to constrain the timescales of magmatic processes. Bennett et al. (2019) demonstrated that various plagioclase textures in mid–ocean ridge basalts can be used as indicators of various magmatic processes. In addition, the composition of plagioclase can provide information about the conditions of the magmatic system, such as the water content, temperature, and melt composition (Hellevang and Pedersen, 2008; Lange et al., 2013; Coote and Shane, 2016). These conditions are related to magma mixing, fractional crystallization, and assimilation processes. Therefore, plagioclase phenocrysts have the potential to record a magmatic history that might be obscured at the whole–rock geochemical scale. However, compared with continental and other oceanic settings, plagioclase phenocrysts are still underutilized in the study of MORBs from ultraslow spreading ridges.
In recent years, the China Ocean Mineral Resources R&D Association (COMRA) has provided support for scientific expeditions to the Southwest Indian Ridge (SWIR), during which a large number of new samples were collected (Tao et al., 2012). Among the various types of MORBs, plagioclase ultraphyric basalts (PUBs), defined by Cullen et al. (1989), are valuable due to their high plagioclase contents (10%–54%). These plagioclase phenocrysts contain unique information about crustal processes and the characteristics of the magmatic array present in the lower oceanic crust (Hellevang and Pedersen, 2008; Lange et al., 2013; Bennett et al., 2019). In this study, we examined the compositions of plagioclase phenocrysts from PUBs erupted on the SWIR (51°E), which is a typical ultraslow spreading ridge. Microanalysis of these phenocrysts provides an opportunity to investigate the geochemical changes that occurred, and thus, provides insights into the details of the magma reservoirs beneath the ultraslow oceanic spreading centers.
The SWIR is a typical ultraslow spreading ridge, with a half–spreading rate of around 7–9 mm/a (Dick et al., 2003). It separates the African Plate from the Antarctic Plate. The SWIR extends for 8 000 km from the Rodrigues Triple Junction (RTJ, 70°E) at its eastern end to the Bouvet Triple Junction at its western end (BTJ, 0°) (Fig. 1a). The SWIR is characterized by strong segmentation and discontinuous magmatism. Bathymetric data has revealed a shallow central region between the Prince Edward Transform Fracture Zone (35.5°E) and the Gallieni Fracture Zone (52.3°E), with an average depth of ~3 200 m compared with the deeper western (~4 000 m) and the eastern (~4 500 m) sections of the ridge (Cannat et al., 2008). As a result of the Marion hotspot to the southwest, this region has strong negative residual mantle Bouguer gravity anomalies, indicating relatively active crust–mantle exchange, deep magmatism, moderate levels of melting, and a moderate heat supply (Georgen et al., 2001; Sauter et al., 2009).
The study area and the sample sites are located between the Indomed (46°E) and Gallieni fracture zones (IFZ–GFZ) on the shallow central region of the SWIR (Fig. 1b). Previous geophysical and geochemical studies have been conducted on the ridge segment between the IFZ and GFZ. The center of this segment has anomalously thick crust (up to 10 km) (Niu et al., 2015) compared with that of the neighboring ridge sections. This thick oceanic crust indicates a robust magma supply in this area, which has been inferred to be associate with the Crozet hotspot (Sauter et al., 2009; Zhang et al., 2013) or with tectonic processes (Jian et al., 2017). Segments with robust magma supplies are also promising areas for hydrothermal activity (Tao et al., 2012).
The PUB sample examined in this study was collected by television–guided grabs (TVGs) during the R/V Dayang Yihao Cruise DY115–21 to the 51°E magmatic segment in 2010. The sample was collected at a water depth of about 1 655 m. Optical microscopy analysis indicates that lavas from the SWIR are porphyritic and contain ~15% phenocrysts. The phenocrysts are mostly euhedral to subhedral plagioclase with polysynthetic twinning. The other main phase is subhedral to anhedral olivine (<1%). The plagioclase crystals contain abundant melt inclusions. The groundmass is primarily composed of suhedral, lath–like plagioclase and subordinate, anhedral olivine (Fig. 2).
Backscattered electron (BSE) imaging of the plagioclase in polished thin sections was used to characterize the textures of the crystals. Mineral analysis of the plagioclase was conducted using the JEOL JXA 8100 electron microprobe at the Key Laboratory of Submarine Geosciences (KLSG), Ministry of Natural Resources (MNR). The analytical conditions were as follows: a 15 kV accelerating voltage, a 20 nA specimen current, and a 1 μm focused beam. The peak counting times were 90 s for Fe and Mg, and 20 s for all of the other major elements. The detection limits for most of the elements, except for Ti, were lower than 400×10–6, depending on the abundance of the elements. The detection limit of Ti was 600×10–6 due to its lower content. The following natural and synthetic standards were used for the specified elements: Olivine (Si, Mg), Apatite (Ca, P), Hematite (Fe), Albite (Na, Al), Orthoclase (K), Rhodonite (Mn), Rutile (Ti), and Tugtupite (Cl). The raw data was corrected using the ZAF correction. The chemical formulas of the plagioclase phenocrysts were calculated from the mineral analysis results based on 24 anions.
The plagioclase phenocrysts (typically >0.5 mm) from the SWIR are predominantly euhedral to subhedral in shape with tabular habits. They exhibit three textural types (Figs 3a–c). Type 1 plagioclase crystals are characterized by oscillatory zoning, surrounded by a thin rim (<50 μm). Type 2 plagioclase crystals are also characterized by oscillatory zoning, but have wider rims (50–150 μm) than Type 1 plagioclase (<50 μm). Type 3 plagioclase crystals contain numerous circular melt inclusions and do not exhibit oscillatory zoning. The plagioclase in the groundmass varies in size. The relatively large microphenocrysts typically have sizes of 0.01–0.50 mm and are primarily unzoned, whereas the relatively small groundmass microlites are <0.01 mm (Fig. 3d).
A total of 25 microprobe analyses were performed on the studied sample. Representative chemical data for the plagioclases are presented in Table 1. These phenocrysts have An contents of 58 to 82. No distinct compositional differences exist among the three types of plagioclase phenocrysts. All of the phenocrysts have calcic cores (An74–82) and sodic rim growth (~An67–71). The variation from the core to the rim is 10–20 mol% An. The plagioclase crystals found within the microphenocrysts and microlites are sodic (An58–63), similar to the rims of the phenocrysts. On the ternary classification diagram, the plagioclase in the studied sample display a continuous range from bytownite to labradorite, with An decreasing from 82 to 58. The cores of the plagioclase phenocrysts are bytownite, while the rims range from bytownite to labradorite. All of the plagioclase crystals in the groundmass are labradorite with lower An contents (Fig. 4).
| Element | Phenocrystal | Groundmass1) | ||||||||||
| Type I | Type II | Type III | Pl-m | Pl-g | ||||||||
| Core | Rim | Core | Rim | Core | Rim | |||||||
| SiO2 | 48.25 | 51.42 | 50.73 | 52.24 | 48.80 | 51.74 | 53.70 | 54.69 | ||||
| TiO2 | 0.06 | 0.10 | 0.05 | 0.00 | 0.00 | 0.11 | 0.08 | 0.17 | ||||
| Al2O3 | 31.49 | 30.03 | 30.28 | 29.61 | 31.70 | 29.81 | 28.41 | 26.72 | ||||
| FeO | 0.29 | 0.40 | 0.50 | 0.54 | 0.38 | 0.62 | 0.58 | 1.41 | ||||
| MnO | 0.04 | 0.00 | 0.04 | 0.00 | 0.00 | 0.00 | 0.00 | 0.00 | ||||
| MgO | 0.15 | 0.19 | 0.24 | 0.18 | 0.20 | 0.20 | 0.18 | 0.37 | ||||
| CaO | 16.80 | 14.52 | 14.94 | 14.01 | 16.29 | 14.13 | 12.77 | 11.92 | ||||
| Na2O | 1.99 | 3.25 | 2.92 | 3.42 | 2.21 | 3.22 | 4.15 | 4.64 | ||||
| K2O | 0.01 | 0.02 | 0.01 | 0.03 | 0.01 | 0.02 | 0.04 | 0.06 | ||||
| P2O5 | 0.01 | 0.03 | 0.00 | 0.02 | 0.01 | 0.00 | 0.01 | 0.05 | ||||
| Total | 99.08 | 99.95 | 99.71 | 100.04 | 99.59 | 99.86 | 99.92 | 100.02 | ||||
| Calculated atoms based on 24 oxygens | ||||||||||||
| Si | 6.704 | 7.031 | 6.962 | 7.127 | 6.736 | 7.072 | 7.311 | 7.457 | ||||
| Ti | 0.006 | 0.010 | 0.005 | 0.000 | 0.000 | 0.011 | 0.008 | 0.017 | ||||
| Al | 5.157 | 4.840 | 4.899 | 4.761 | 5.157 | 4.803 | 4.559 | 4.294 | ||||
| Fe | 0.034 | 0.046 | 0.057 | 0.062 | 0.043 | 0.071 | 0.066 | 0.161 | ||||
| Mn | 0.005 | 0.000 | 0.004 | 0.000 | 0.000 | 0.000 | 0.000 | 0.000 | ||||
| Mg | 0.031 | 0.039 | 0.049 | 0.036 | 0.041 | 0.041 | 0.037 | 0.074 | ||||
| Ca | 2.501 | 2.127 | 2.198 | 2.047 | 2.409 | 2.070 | 1.863 | 1.742 | ||||
| Na | 0.536 | 0.862 | 0.776 | 0.905 | 0.591 | 0.852 | 1.096 | 1.225 | ||||
| K | 0.002 | 0.003 | 0.002 | 0.005 | 0.001 | 0.003 | 0.007 | 0.011 | ||||
| P | 0.002 | 0.003 | 0.000 | 0.002 | 0.001 | 0.000 | 0.001 | 0.005 | ||||
| Total | 14.977 | 14.961 | 14.953 | 14.945 | 14.980 | 14.924 | 14.947 | 14.987 | ||||
| An | 0.82 | 0.71 | 0.74 | 0.69 | 0.80 | 0.71 | 0.63 | 0.58 | ||||
| Ab | 0.18 | 0.29 | 0.26 | 0.31 | 0.20 | 0.29 | 0.37 | 0.41 | ||||
| Or | 0.00 | 0.00 | 0.00 | 0.00 | 0.00 | 0.00 | 0.00 | 0.01 | ||||
| FeO/MgO | 1.91 | 2.12 | 2.08 | 3.04 | 1.90 | 3.06 | 3.16 | 3.87 | ||||
| Ca/Na | 4.67 | 2.47 | 2.83 | 2.26 | 4.07 | 2.43 | 1.70 | 1.42 | ||||
| T/°C2) | 1296 | 1108 | 1237 | 1098 | 1280 | 1129 | 1063 | 1042 | ||||
| Note: 1) Pl-m and Pl-g represent plagioclase microphenocrysts and microlites in groundmass, respectively; 2) crystallization temperature is calculated according to Kudo and Weill (1983). | ||||||||||||
DownLoad:
CSV
In terms of the major elements, the FeO and MgO concentrations of the plagioclase do not vary significantly. The FeO and MgO concentrations are 0.29%–1.41% (wt) and 0.15%–0.37% (wt), respectively. The FeO and MgO concentrations of the plagioclase cores are relatively depleted compared with those of the plagioclase rims and groundmass.
Plagioclase compositions are a useful indicator of crystallization temperature (Kudo and Weill, 1970; Mollo et al., 2011). In this study, the plagioclase–melt geothermometry method proposed by Kudo and Weill (1970) was used to estimate the crystallization temperature of the plagioclase. Before applying the geothermometry method, the pressures must be determined. According to Chen et al. (2002), the pressures of the plagioclase rims and cores are approximately 0.5×108 and 1.0×108 Pa, respectively. Similarly, the plagioclase in the groundmass is estimated to have crystallized at shallower depths, within the upper crust or on the seabed. The pressure of the groundmass is also assumed to be 0.5×108 Pa according to Chen et al. (2002). Based on the assumptions stated above, Eqs (1) and (2) were used to calculate the crystallization temperatures of the plagioclase crystals:
| $$ \frac{{{\rm{ln}}\lambda }}{{{\sigma }}} + 1.29 \times \frac{{{{10}^4}\varphi }}{T} = 9.87 \times {10^{ - 3}}T - 15.21,\left( {{P_{{{\rm{H}}_{\rm{2}}}{\rm{O}}}} = 0.5\times 10^8\;{\rm {Pa}}} \right), $$ | (1) |
| $$ \frac{{{\rm{ln}}\lambda }}{\sigma } + 1.29 \times \frac{{{{10}^4}\varphi }}{T} = 9.60 \times {10^{ - 3}}T - 15.76,\left( {{P_{{{\rm{H}}_{\rm{2}}}{\rm{O}}}} = 1.0\times 10^8\;{\rm {Pa}} } \right), $$ | (2) |
where
According to the plagioclase geothermometer described above, the crystallization temperatures of the phenocryst cores and rims are (1 273±18)°C and (1 099±10)°C, respectively. From the core to the rim of the phenocryst, the crystallization temperature decreases by about 200°C. The crystallization temperatures of the microphenocrysts and microlites in the groundmass range from 1 063°C to 1 087°C (average = 1 072°C), which is similar to the crystallization temperature of the phenocryst rims (Table 1).
In previous studies, the crystallization temperatures of high–An (An≥70%) plagioclase phenocrysts, which were estimated from the entrapment temperature of melt inclusions in samples from ultraslow spreading ridges, were found to range from 1 230°C to 1 260°C (Nielsen et al., 1995; Drignon et al., 2019). Whereas the minimum estimation of the crystallization temperature of low An (An<70%) plagioclase is 1 100°C (Yang et al., 2019). The crystallization temperatures in our study are consistent with these results, which suggests that our calculated results are reasonable and the plagioclases thermometer by Kudo and Weill (1970) can be used to calculate the crystallization temperatures of plagioclase phenocrysts from MORBs beneath ultraslow ridges.
In plagioclase, the diffusion rates of major elements, such as the NaSi–CaAl exchange, are extremely slow (Grove et al., 1984). The geochemical zoning of the crystals likely reflects the magmatic conditions, and thus, it can be used to investigate the likelihood of plagioclase–melt equilibrium. Previous experiments have demonstrated that the partition coefficient of Ca/Na (KCa/Na) between plagioclase and melt mainly positively depends on the water content of the magma (Sisson and Grove, 1993; Martel et al., 2006). For basaltic magmas from mid-ocean ridges, when the magmatic water content is 3%, KCa/Na is typically ~1 for mid to upper lithospheric pressures (<10×108 Pa) (Sisson and Grove, 1993). The magmatic water contents of the SWIR lavas are close to the global average value for the upper mantle (0.3%–0.4%) (Robinson et al., 2001), which indicates that their KCa/Na value is less than 1.
In addition to the KCa/Na value, a representative melt composition is required to assess the equilibrium composition of the plagioclase phenocrysts. The groundmass of these rocks is microcrystalline and lacks glass, so it represents the final melt. Thus, as was discussed above, the average whole–rock composition of the aphyric basalts from the same segment is representative of the groundmass composition. The reported composition of the basalts and glass from this ridge segment do not vary significantly, with molar Ca/Na values ranging from 2.11 to 2.68 (Yang et al., 2014; Bézos and Humler, 2005). The molar Ca/Na values of the plagioclase phenocrysts in equilibrium with the melt are consistently≤2.68 (Fig. 5). The molar Ca/Na values of the plagioclase rims and groundmass range from 1.42 to 2.47, which is in equilibrium with the melt (Fig. 5). However, the compositions of the plagioclase cores (Ca/Na = 3 to 4) are higher than the upper limit of the equilibrium melt (Ca/Na = 2.68) (Fig. 5). In addition, the crystallization temperatures of the plagioclase cores (average of 1 273°C) are close to the experimental melting point of basaltic magma (~1 300°C) and the estimated upper mantle potential temperature of the SWIR lavas (~1 280°C) (Kamenetsky et al., 2000; Robinson et al., 2001). Therefore, the plagioclase cores are unlikely to have formed in the host magma. Instead, they are most likely xenocrysts, which crystallized from a more calcic melt.
It is worth noting that compared with the plagioclase rims, the plagioclase phenocrysts in the groundmass have lower Ca/Na ratios (<2). Their more sodic compositions are most likely due to decreasing magmatic water contents caused by the fact that the KCa/Na between the plagioclase and the melt decreases as of the water content of the magma increases (Martel and Schmidt, 2003). Water is lost during the late-stage of magma ascent through the conduit due to ascent–driven decompression of the water-saturated magma. Thus, the compositional variations can be explained by the variations in the magmatic water content.
The MgO zoning patterns of the plagioclase phenocrysts have the potential to record the magmatic composition because the Mg contents of the plagioclase phenocrysts reflect the composition of the host melt (Ginibre et al., 2002). Although it is difficult to accurately estimate the Mg partition coefficient (KMg) between plagioclase and melt, empirical studies have suggested an Arrhenius-like relationship of decreasing KMg with increasing XAn (Bindeman et al., 1998). To investigate the equilibrium relationship between the SWIR plagioclase and the host magma, we estimated the KMg based on this empirical relationship, which is often used to determine whether plagioclase phenocrysts are in equilibrium with the melt in magmatic systems (Bindeman et al., 1998; Coote and Shane, 2016). The majority of the plagioclase phenocrysts analyzed in this study are more enriched in MgO than the modeled melt compositions for a crystallization temperature of ~1099°C, which is based on the plagioclase geothermometer described above (Fig. 6a). This implies that the plagioclase in the PUBs analyzed in this study could not have crystallized from the basaltic host melt (MgO of 6.38%–8.87% (wt), average of 7.77% (wt)) and would require a more mafic melt (MgO of up to 14% (wt)). This agrees with the major element modeling (Ca/Na) of the plagioclase cores, but does not agree with the equilibrium between the rims and the melt suggested by the Ca/Na ratios (Fig. 5).
The rimward increase in Mg and decrease in An appear to be similar to the trends described in previous diffusion studies (Costa et al., 2003; Moore et al., 2014). When minerals crystallize due to large degrees of undercooling and the diffusion in the melt cannot keep pace with the crystal formation, it is possible for late–stage rapid disequilibrium crystallization to produce a boundary layer melt enriched in incompatible (Fe, Mg) elements (Ginibre et al., 2002; Coote et al., 2018). This could explain the elevated MgO contents of the plagioclase rims and the groundmass plagioclase relative to the equilibrium values (Fig. 6a). The internal disequilibrium in the MgO contents can be explained by the mixing of primitive (MgO=14% (wt)) and more evolved magmas (MgO=7.77%(wt)).
Additionally, in plagioclase, the post–crystallization diffusion of Fe is slower than that of Mg (Costa et al., 2003). Thus, the highest FeO/MgO ratios occur in the outermost parts of the rims and in the groundmass plagioclase (Fig. 6b) due to the fact that FeO diffuses into the plagioclase more slowly than MgO, and thus, more FeO accumulates in the boundary layer during a short residence time. However, it is worth noting that since the FeO content of the plagioclase depends on both the melt composition and the oxygen fugacity, it is difficult to assess the plagioclase–melt equilibrium using only the FeO content (Coote et al., 2018). Therefore, disequilibrium diffusion within a short time period can result in elevated Mg contents compared with the equilibrium values and the highest FeO/MgO ratios occurring in the plagioclase rims (Coote and Shane, 2016; Moore et al., 2014). Whereas the discrepancies in the plagioclase–melt equilibrium inferred from the major elements (Ca/Na) and the Mg contents are likely an artifact of disequilibrium diffusion (Coote and Shane, 2016).
In the SWIR, the crystallization temperature and composition of the plagioclase cores are distinctly different from those of the rims (Table 1), which indicates that the plagioclase cores and rims have different thermal histories.
The plagioclase cores are usually uniform and exhibit oscillatory zoning. The oscillatory zoning results in small-scale compositional variations, which suggests a regime of near–constant intensive parameters (pressure, temperature) (Landi et al., 2004; Shcherbakov et al., 2011). These characteristics suggest the cores crystallized from a stable environment and do not have complex crystallization histories. Besides, the presence of cores in plagioclase that are too primitive to be in equilibrium with the host magma indicates that they crystallized from a more primitive region in the plumbing system and were picked up by a more evolved melt later. Thus, we propose that the plagioclase cores grew in a stable mush zone where the temperature was high and constant, and were later entrained into a more evolved melt.
The thin plagioclase rims (normally < 150 μm) have lower crystallization temperatures and lower An values than the cores, and their major elements (Ca/Na) are in equilibrium with the host magma. These characteristics demonstrate that they crystallized from the host melt, which is more evolved than the magma from which the cores crystalized. Similarly, the plagioclase microphenocrysts and microlites are also in equilibrium (Ca/Na) values with the host magma, which suggests that they crystallized from the host magma as well. Compared with fast to intermediate spreading ridges, there are generally no stable magma chambers and the volume of melt may be low under ultraslow mid–ocean ridges (Dick et al., 2003). Thus, the transport of magma into cooler regions of the reservoir would be expected to result in abrupt, strong undercooling of the magma. Due to these large degrees of rapid cooling, the plagioclase phenocrysts in SWIR PUBs generally have thin rims. The plagioclase rims of the phenocrysts exhibit major element (Fe, Mg) enrichment because these elements cannot reach equilibrium when diffusing from the host magma into the plagioclase phenocrysts. The MgO contents of the outermost rims of the plagioclase can be used to calculate the maximum time between the incorporation of the plagioclase into the host melt and quenching on the seafloor (Costa et al., 2003; Moore et al., 2014). We used Eq. (8) in Costa et al. (2003) to determine the Mg diffusion coefficient of the plagioclase. The maximum residence time of the plagioclases in the host melt was nearly 5–8 d. The close temporal relationship between the evolved magma replenishment and the eruption suggests that replenishment plays an important role in driving the eruption, which has also been suggested for other MORB eruptions on slow and intermediate ridges (Costa et al., 2010). Overall, our favored model is that replenishment by an evolved melt under the SWIR ridges (51°E) drives the eruption over a short period of time.
(1) The plagioclase cores with high An values have higher crystallization temperatures (1 273±18)°C than the rims (1 099±10)°C. The range of crystallization temperatures for the microphenocrysts and microlites in the groundmass is similar to that of the phenocryst rims.
(2) The compositions of the plagioclase cores from the SWIR indicate that they did not form in the host magma, but xenocrysts are crystallized from a more mafic melt composition. Whereas the plagioclase rims and the microphenocrysts and microlites in the groundmass are in equilibrium with the host basaltic melts.
(3) The disequilibrium MgO contents and the higher FeO/MgO ratios of the rims of the plagioclase phenocrysts reflect shorter magmatic residence time periods than would be resulting from equilibrium diffusion.
(4) An evolved melt replenished the magma under the SWIR ridges (51°E), driving the eruption over a short period of time.
We are grateful to Yin-Jia Jin and two anonymous reviewers for their careful editing and constructive comments, which improved the manuscript. We also thank the crew and scientists involved in the R/V Dayang Yihao Cruise DY115–21.
| [1] |
Barth A, Beckers J M, Troupin C, et al. 2014. Divand-1.0: n-dimensional variational data analysis for ocean observations. Geoscientific Model Development, 7(1): 225–241. doi: 10.5194/gmd-7-225-2014
|
| [2] |
Beismann J O, Käse R H, Lutjeharms J R E. 1999. On the influence of submarine ridges on translation and stability of Agulhas rings. Journal of Geophysical Research: Oceans, 104(C4): 7897–7906. doi: 10.1029/1998JC900127
|
| [3] |
Chaigneau A, Gizolme A, Grados C. 2008. Mesoscale eddies off Peru in altimeter records: Identification algorithms and eddy spatio-temporal patterns. Progress in Oceanography, 79(2–4): 106–119
|
| [4] |
Chaigneau A, Le Texier M, Eldin G, et al. 2011. Vertical structure of mesoscale eddies in the eastern South Pacific Ocean: A composite analysis from altimetry and Argo profiling floats. Journal of Geophysical Research: Oceans, 116(C11): C11025. doi: 10.1029/2011JC007134
|
| [5] |
Chelton D B, Schlax M G, Samelson R M. 2011. Global observations of nonlinear mesoscale Eddies. Progress in Oceanography, 91(2): 167–216. doi: 10.1016/j.pocean.2011.01.002
|
| [6] |
Chen Gengxin, Hou Yijun, Chu Xiaoqing. 2011. Mesoscale eddies in the South China Sea: Mean properties, spatiotemporal variability, and impact on thermohaline structure. Journal of Geophysical Research: Oceans, 116(C6): C06018
|
| [7] |
Dai Jun, Wang Huizan, Zhang Weimin, et al. 2020. Observed spatiotemporal variation of three-dimensional structure and heat/salt transport of anticyclonic mesoscale eddy in Northwest Pacific. Journal of Oceanology and Limnology, 38(6): 1654–1675. doi: 10.1007/s00343-019-9148-z
|
| [8] |
Dong Di, Brandt P, Chang Ping, et al. 2017. Mesoscale eddies in the northwestern Pacific Ocean: Three-dimensional eddy structures and heat/salt transports. Journal of Geophysical Research: Oceans, 122(12): 9795–9813. doi: 10.1002/2017JC013303
|
| [9] |
Early J J, Samelson R M, Chelton D B. 2011. The evolution and propagation of quasigeostrophic ocean eddies. Journal of Physical Oceanography, 41(8): 1535–1555. doi: 10.1175/2011JPO4601.1
|
| [10] |
Ferrari R, Wunsch C. 2009. Ocean circulation kinetic energy: Reservoirs, sources, and sinks. Annual Review of Fluid Mechanics, 41: 253–282. doi: 10.1146/annurev.fluid.40.111406.102139
|
| [11] |
Ferron B. 2011. A 4D-variational approach applied to an eddy-permitting North Atlantic configuration: Synthetic and real data assimilation of altimeter observations. Ocean Modelling, 39(3–4): 370–385
|
| [12] |
Gao Shan, Wang Fan, Li Mingkui, et al. 2008. Application of altimetry data assimilation on mesoscale eddies simulation. Science in China Series D: Earth Sciences, 51(1): 124–151
|
| [13] |
He Yinghui, Cai Shuqun, Wang Dongxiao, et al. 2015. A model study of Luzon cold eddies in the northern South China Sea. Deep Sea Research Part I: Oceanographic Research Papers, 97: 107–123. doi: 10.1016/j.dsr.2014.12.007
|
| [14] |
He Zhongjie, Xie Yuanfu, Li Wei, et al. 2008. Application of the sequential three-dimensional variational method to assimilating SST in a global ocean model. Journal of Atmospheric and Oceanic Technology, 25(6): 1018–1033. doi: 10.1175/2007JTECHO540.1
|
| [15] |
Hoskins B J, Draghici I, Davies H C. 1978. A new look at the ω-equation. Quarterly Journal of the Royal Meteorological Society, 104(439): 31–38. doi: 10.1002/qj.49710443903
|
| [16] |
Kamenkovich V M, Leonov Y P, Nechaev D A, et al. 1996. On the influence of bottom topography on the Agulhas eddy. Journal of Physical Oceanography, 26(6): 892–912. doi: 10.1175/1520-0485(1996)026<0892:OTIOBT>2.0.CO;2
|
| [17] |
Li Cheng, Zhang Zhiwei, Zhao Wei, et al. 2017. A statistical study on the subthermocline submesoscale eddies in the northwestern Pacific Ocean based on Argo data. Journal of Geophysical Research: Oceans, 122(5): 3586–3598. doi: 10.1002/2016JC012561
|
| [18] |
Lin Xiayan, Guan Yuping, Liu Yu. 2013. Three-dimensional structure and evolution process of Dongsha Cold Eddy during autumn 2000. Journal of Tropical Oceanography (in Chinese), 32(2): 55–65
|
| [19] |
Liu Danian, Wang Fan, Zhu Jiang, et al. 2020. Impact of assimilation of moored velocity data on low-frequency current estimation in Northwestern Tropical Pacific. Journal of Geophysical Research: Oceans, 125(9): e2019JC015829
|
| [20] |
Liu Danian, Zhu Jiang, Shu Yeqiang, et al. 2018a. Targeted observation analysis of a Northwestern Tropical Pacific Ocean mooring array using an ensemble-based method. Ocean Dynamics, 68(9): 1109–1119. doi: 10.1007/s10236-018-1188-y
|
| [21] |
Liu Danian, Zhu Jiang, Shu Yeqiang, et al. 2018b. Model-based assessment of a Northwestern Tropical Pacific moored array to monitor intraseasonal variability. Ocean Modelling, 126: 1–12. doi: 10.1016/j.ocemod.2018.04.001
|
| [22] |
Ma Xiaohui, Jing Zhao, Chang Ping, et al. 2016. Western boundary currents regulated by interaction between ocean eddies and the atmosphere. Nature, 535(7613): 533–537. doi: 10.1038/nature18640
|
| [23] |
Mahadevan A. 2016. The impact of submesoscale physics on primary productivity of plankton. Annual Review of Marine Science, 8: 161–184. doi: 10.1146/annurev-marine-010814-015912
|
| [24] |
Mahadevan A, Thomas L N, Tandon A. 2008. Comment on “Eddy/wind interactions stimulate extraordinary mid-ocean plankton blooms”. Science, 320(5875): 448
|
| [25] |
Martin A P, Richards K J. 2001. Mechanisms for vertical nutrient transport within a North Atlantic mesoscale eddy. Deep Sea Research Part II: Topical Studies in Oceanography, 48(4–5): 757–773
|
| [26] |
McGillicuddy D J Jr, Anderson L A, Bates N R, et al. 2007. Eddy/wind interactions stimulate extraordinary mid-ocean plankton blooms. Science, 316(5827): 1021–1026. doi: 10.1126/science.1136256
|
| [27] |
McWilliams J C, Graves L P, Montgomery M T. 2003. A formal theory for vortex rossby waves and vortex evolution. Geophysical & Astrophysical Fluid Dynamics, 97(4): 275–309
|
| [28] |
Moore A M, Arango H G, Broquet G, et al. 2011a. The Regional Ocean Modeling System (ROMS) 4-dimensional variational data assimilation systems: Part I—System overview and formulation. Progress in Oceanography, 91(1): 34–49. doi: 10.1016/j.pocean.2011.05.004
|
| [29] |
Moore A M, Arango H G, Broquet G, et al. 2011b. The Regional Ocean Modeling System (ROMS) 4-dimensional variational data assimilation systems: Part II—performance and application to the California Current System. Progress in Oceanography, 91(1): 50–73. doi: 10.1016/j.pocean.2011.05.003
|
| [30] |
Nardelli B B. 2013. Vortex waves and vertical motion in a mesoscale cyclonic eddy. Journal of Geophysical Research: Oceans, 118(10): 5609–5624. doi: 10.1002/jgrc.20345
|
| [31] |
Ni Qinbiao. 2014. Statistical characteristics and composite three-dimensional structures of mesoscale eddies near the Luzon Strait (in Chinese) [dissertation]. Xiamen: Xiamen University
|
| [32] |
Ni Qinbiao. 2019. Study on eddy movement in the ocean (in Chinese) [dissertation]. Xiamen: Xiamen University
|
| [33] |
Oka E, Kouketsu S, Toyama K, et al. 2011. Formation and subduction of central mode water based on profiling float data, 2003–08. Journal of Physical Oceanography, 41(1): 113–129. doi: 10.1175/2010JPO4419.1
|
| [34] |
Powell B S, Arango H G, Moore A M, et al. 2008. 4DVAR data assimilation in the intra-Americas Sea with the Regional Ocean Modeling System (ROMS). Ocean Modeling, 25(3–4): 173–188
|
| [35] |
Qiu Bo, Chen Shuiming, Klein P, et al. 2020. Reconstructing upper-ocean vertical velocity field from sea surface height in the presence of unbalanced motion. Journal of Physical Oceanography, 50(1): 55–79. doi: 10.1175/JPO-D-19-0172.1
|
| [36] |
Rubio A, Barnier B, Jordá G, et al. 2009. Origin and dynamics of mesoscale eddies in the Catalan Sea (NW Mediterranean): Insight from a numerical model study. Journal of Geophysical Research: Oceans, 114(C6): C06009
|
| [37] |
Sasaki H, Klein P, Qiu Bo, et al. 2014. Impact of oceanic-scale interactions on the seasonal modulation of ocean dynamics by the atmosphere. Nature Communications, 5: 5636. doi: 10.1038/ncomms6636
|
| [38] |
Shu Yeqiang, Chen Ju, Li Shuo, et al. 2019. Field-observation for an anticyclonic mesoscale eddy consisted of twelve gliders and sixty-two expendable probes in the northern South China Sea during summer 2017. Science China Earth Sciences, 62(2): 451–458. doi: 10.1007/s11430-018-9239-0
|
| [39] |
Shu Yeqiang, Wang Dongxiao, Zhu Jiang, et al. 2011. The 4-D structure of upwelling and Pearl River plume in the northern South China Sea during summer 2008 revealed by a data assimilation model. Ocean Modelling, 36(3–4): 228–241
|
| [40] |
Souza J M A C, De Boyer Montégut C, Le Traon P Y. 2011. Comparison between three implementations of automatic identification algorithms for the quantification and characterization of mesoscale eddies in the South Atlantic Ocean. Ocean Science, 7(3): 317–334. doi: 10.5194/os-7-317-2011
|
| [41] |
Thompson P D. 2010. Reduction of analysis error through constraints of dynamical consistency. Journal of Applied Meteorology, 8(5): 738–742
|
| [42] |
Troupin C, Barth A, Sirjacobs D, et al. 2012. Generation of analysis and consistent error fields using the Data Interpolating Variational Analysis (DIVA). Ocean Modelling, 52–53: 90–101
|
| [43] |
Wang Guihua, Su Jilan, Chu P C. 2003. Mesoscale eddies in the South China Sea observed with altimeter data. Geophysical Research Letters, 30(21): 2121. doi: 10.1029/2003GL018532
|
| [44] |
Wang Lei, Gan Jianping. 2014. Delving into three-dimensional structure of the West Luzon Eddy. Deep Sea Research Part I: Oceanographic Research Papers, 90: 48–61. doi: 10.1016/j.dsr.2014.04.011
|
| [45] |
Warner J C, Sherwood C R, Arango H G, et al. 2005. Performance of four turbulence closure models implemented using a generic length scale method. Ocean Modelling, 8(1–2): 81–113
|
| [46] |
Xu Dazhi. 2012. Research on predictability of mesoscale eddies in the northern South China Sea based on data assimilation (in Chinese) [dissertation]. Guangzhou: South China Sea Institute of Oceanology
|
| [47] |
Yang Guang, Wang Fan, Li Yuanlong, et al. 2013. Mesoscale eddies in the northwestern subtropical Pacific Ocean: Statistical characteristics and three-dimensional structures. Journal of Geophysical Research: Oceans, 118(4): 1906–1925. doi: 10.1002/jgrc.20164
|
| [48] |
Zhang Zhiwei, Li Peiliang, Xu Lixiao, et al. 2015. Subthermocline eddies observed by rapid-sampling Argo floats in the subtropical northwestern Pacific Ocean in Spring 2014. Geophysical Research Letters, 42(15): 6438–6445. doi: 10.1002/2015GL064601
|
| [49] |
Zhang Wenzhou, Ni Qinbiao, Xue Huijie. 2018. Composite eddy structures on both sides of the Luzon Strait and influence factors. Ocean Dynamics, 68(11): 1527–1541. doi: 10.1007/s10236-018-1207-z
|
| [50] |
Zhang Zhiwei, Tian Jiwei, Qiu Bo, et al. 2016. Observed 3D structure, generation, and dissipation of oceanic mesoscale eddies in the South China Sea. Scientific Reports, 6: 24349. doi: 10.1038/srep24349
|
| [51] |
Zhang Zhengguang, Wang Wei, Qiu Bo. 2014. Oceanic mass transport by mesoscale eddies. Science, 345(6194): 322–324. doi: 10.1126/science.1252418
|
| [52] |
Zhang Zhiwei, Zhang Xincheng, Qiu Bo, et al. 2020a. Submesoscale currents in the subtropical upper ocean observed by long-term high-resolution mooring arrays. Journal of Physical Oceanography, 51(1): 187–206. doi: 10.1175/JPO-D-20-0100.1
|
| [53] |
Zhang Zhiwei, Zhang Yuchen, Qiu Bo, et al. 2020b. Spatiotemporal characteristics and generation mechanisms of submesoscale currents in the northeastern South China Sea revealed by numerical simulations. Journal of Geophysical Research: Oceans, 125(2): e2019JC015404
|
| [54] |
Zhang Zhengguang, Zhang Yu, Wang Wei, et al. 2013. Universal structure of mesoscale eddies in the ocean. Geophysical Research Letters, 40(14): 3677–3681. doi: 10.1002/grl.50736
|
| [55] |
Zhao Fu, Zhang Yunfei, Zhu Xueming, et al. 2017. An assimilative numerical study of the paired cold and warm mesoscale eddies during winter in the southwest of Taiwan. Marine Forecasts (in Chinese), 34(5): 1–15
|
| [56] |
Zhong Yisen, Bracco A, Tian Jiwei, et al. 2017. Observed and simulated submesoscale vertical pump of an anticyclonic eddy in the South China Sea. Scientific Reports, 7: 44011. doi: 10.1038/srep44011
|
| 1. | Abhijit Shee, Sourav Sil, Rahul Deogharia. Three-dimensional characteristics of mesoscale eddies in the western boundary current region of the Bay of Bengal using ROMS-NPZD. Dynamics of Atmospheres and Oceans, 2024, 105: 101424. doi:10.1016/j.dynatmoce.2023.101424 | |
| 2. | Vishwanath Boopathi, Sachiko Mohanty. Mesoscale eddy variability over the Bay of Bengal in response to the contrasting phases of extreme Indian Ocean Dipole events. Regional Studies in Marine Science, 2024. doi:10.1016/j.rsma.2024.103604 | |
| 3. | Zheliang Zhang, Yunxia Zheng, Hao Li. Imprints of tropical cyclone on three-dimensional structural characteristics of mesoscale oceanic eddies. Frontiers in Earth Science, 2023, 10 doi:10.3389/feart.2022.1057798 | |
| 4. | Gongfu Zhou, Guijun Han, Wei Li, et al. High‐Resolution Gridded Temperature and Salinity Fields From Argo Floats Based on a Spatiotemporal Four‐Dimensional Multigrid Analysis Method. Journal of Geophysical Research: Oceans, 2023, 128(5) doi:10.1029/2022JC019386 | |
| 5. | Zhihui Chen, Pinqiang Wang, Senliang Bao, et al. Rapid reconstruction of temperature and salinity fields based on machine learning and the assimilation application. Frontiers in Marine Science, 2022, 9 doi:10.3389/fmars.2022.985048 |
Figures(21) / Tables(1)
Supported by:
Beijing Renhe Information Technology Co. Ltd
Jun Dai, Huizan Wang, Weimin Zhang, Pinqiang Wang, Tengling Luo. Three-dimensional structure of an observed cyclonic mesoscale eddy in the Northwest Pacific and its assimilation experiment[J]. Acta Oceanologica Sinica, 2021, 40(5): 1-19. doi: 10.1007/s13131-021-1810-6
| Element | Phenocrystal | Groundmass1) | ||||||||||
| Type I | Type II | Type III | Pl-m | Pl-g | ||||||||
| Core | Rim | Core | Rim | Core | Rim | |||||||
| SiO2 | 48.25 | 51.42 | 50.73 | 52.24 | 48.80 | 51.74 | 53.70 | 54.69 | ||||
| TiO2 | 0.06 | 0.10 | 0.05 | 0.00 | 0.00 | 0.11 | 0.08 | 0.17 | ||||
| Al2O3 | 31.49 | 30.03 | 30.28 | 29.61 | 31.70 | 29.81 | 28.41 | 26.72 | ||||
| FeO | 0.29 | 0.40 | 0.50 | 0.54 | 0.38 | 0.62 | 0.58 | 1.41 | ||||
| MnO | 0.04 | 0.00 | 0.04 | 0.00 | 0.00 | 0.00 | 0.00 | 0.00 | ||||
| MgO | 0.15 | 0.19 | 0.24 | 0.18 | 0.20 | 0.20 | 0.18 | 0.37 | ||||
| CaO | 16.80 | 14.52 | 14.94 | 14.01 | 16.29 | 14.13 | 12.77 | 11.92 | ||||
| Na2O | 1.99 | 3.25 | 2.92 | 3.42 | 2.21 | 3.22 | 4.15 | 4.64 | ||||
| K2O | 0.01 | 0.02 | 0.01 | 0.03 | 0.01 | 0.02 | 0.04 | 0.06 | ||||
| P2O5 | 0.01 | 0.03 | 0.00 | 0.02 | 0.01 | 0.00 | 0.01 | 0.05 | ||||
| Total | 99.08 | 99.95 | 99.71 | 100.04 | 99.59 | 99.86 | 99.92 | 100.02 | ||||
| Calculated atoms based on 24 oxygens | ||||||||||||
| Si | 6.704 | 7.031 | 6.962 | 7.127 | 6.736 | 7.072 | 7.311 | 7.457 | ||||
| Ti | 0.006 | 0.010 | 0.005 | 0.000 | 0.000 | 0.011 | 0.008 | 0.017 | ||||
| Al | 5.157 | 4.840 | 4.899 | 4.761 | 5.157 | 4.803 | 4.559 | 4.294 | ||||
| Fe | 0.034 | 0.046 | 0.057 | 0.062 | 0.043 | 0.071 | 0.066 | 0.161 | ||||
| Mn | 0.005 | 0.000 | 0.004 | 0.000 | 0.000 | 0.000 | 0.000 | 0.000 | ||||
| Mg | 0.031 | 0.039 | 0.049 | 0.036 | 0.041 | 0.041 | 0.037 | 0.074 | ||||
| Ca | 2.501 | 2.127 | 2.198 | 2.047 | 2.409 | 2.070 | 1.863 | 1.742 | ||||
| Na | 0.536 | 0.862 | 0.776 | 0.905 | 0.591 | 0.852 | 1.096 | 1.225 | ||||
| K | 0.002 | 0.003 | 0.002 | 0.005 | 0.001 | 0.003 | 0.007 | 0.011 | ||||
| P | 0.002 | 0.003 | 0.000 | 0.002 | 0.001 | 0.000 | 0.001 | 0.005 | ||||
| Total | 14.977 | 14.961 | 14.953 | 14.945 | 14.980 | 14.924 | 14.947 | 14.987 | ||||
| An | 0.82 | 0.71 | 0.74 | 0.69 | 0.80 | 0.71 | 0.63 | 0.58 | ||||
| Ab | 0.18 | 0.29 | 0.26 | 0.31 | 0.20 | 0.29 | 0.37 | 0.41 | ||||
| Or | 0.00 | 0.00 | 0.00 | 0.00 | 0.00 | 0.00 | 0.00 | 0.01 | ||||
| FeO/MgO | 1.91 | 2.12 | 2.08 | 3.04 | 1.90 | 3.06 | 3.16 | 3.87 | ||||
| Ca/Na | 4.67 | 2.47 | 2.83 | 2.26 | 4.07 | 2.43 | 1.70 | 1.42 | ||||
| T/°C2) | 1296 | 1108 | 1237 | 1098 | 1280 | 1129 | 1063 | 1042 | ||||
| Note: 1) Pl-m and Pl-g represent plagioclase microphenocrysts and microlites in groundmass, respectively; 2) crystallization temperature is calculated according to Kudo and Weill (1983). | ||||||||||||
DownLoad:
CSV
DownLoad:
