• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    Particle dynamics revealed by 210Po/210Pb disequilibria around Prydz Bay, the Southern Ocean in summer

    2022-07-20 01:34:04CHENMengyaCHENMinZHENGMinfangQIUYushengZHUJingQIANQiankun
    Advances in Polar Science 2022年1期

    CHEN Mengya, CHEN Min, ZHENG Minfang, QIU Yusheng, ZHU Jing & QIAN Qiankun

    Particle dynamics revealed by210Po/210Pb disequilibria around Prydz Bay, the Southern Ocean in summer

    CHEN Mengya, CHEN Min*, ZHENG Minfang, QIU Yusheng, ZHU Jing & QIAN Qiankun

    College of Ocean and Earth Sciences, Xiamen University, Xiamen 361102, China

    Seawater samples were collected around Prydz Bay in summer of 2014, dissolved and particulate210Po and210Pb were measured to reveal the disequilibrium characteristics and particle dynamics. Our results show that the distribution of210Po and210Po/210Pb activity ratio in the upper water is mainly affected by biological absorption or particle adsorption. An abnormal excess of210Po relative to210Pb was observed in the surface water at stations P1-2 and P2-2, which is likely to be the horizontal transport of water mass with high DPo/DPb)A.R.and TPo/TPb)A.R.. In this study, the removal of particulate210Po is mainly controlled by the scavenging of dissolved210Po and the two have a linear positive correlation with the salinity, a negative linear correlation with the content of dissolved oxygen and a reciprocal relationship with the content of POC. The export flux of POC at 100 m is estimated to be 1.8–4.4 mmol·m?2·d?1(avg. 2.9 mmol·m?2·d?1) based on210Po/210Pb disequilibria, with the highest value in the shelf, which is consistent with the distribution of biological productivity.

    210Po/210Pb disequilibria,210Po excess, particle dynamics, POC export, Prydz Bay

    1 Introduction

    210Po and210Pb belong to the238U decay series, of which210Pb is the direct parent of210Po. Since the half-life of210Po (138.4 d) is much shorter than210Pb (22.3 a), the two will reach radioactive equilibrium in about 2.5 a without biogeochemical fractionation. However,210Po is more easily absorbed by plankton or adsorbed by particulate organic matter, while210Pb has a stronger affinity for inorganic components (Bacon et al., 1976; Nozaki et al., 1998). The difference in the biogeochemical behavior of210Po and210Pb leads to the variation in spatial distribution and the210Po/210Pb disequilibria.210Po in the upper ocean is generally deficit with respect to210Pb. Two reasons are proposed to be responsible for this, one is that the atmospheric deposition flux of210Pb is much higher than that of210Po, and the other is that210Po is more easily absorbed by organisms or adsorbed by particles (Nozaki et al., 1998; Yang et al., 2006; Roca-Martí et al., 2018). However, excessive210Po was occasionally found in the upper waters of some sea areas(Kadko, 1993; Radakovitch et al., 1998), but the reason is not clear.Although extensive studies have been conducted on the distribution of210Po and210Pb in the ocean in recent years, due to sampling limitations, there are few reports in the Southern Ocean, especially around Prydz Bay. Yin et al. (2004) reported that the specific activities of total210Po and210Pb in the surface water around Prydz Bay ranged from 0.80 to 1.52 Bq·m?3and from 1.34 to 2.15 Bq·m?3, with averages of 1.14 and 1.66 Bq·m?3, respectively. The210Po/210Pb activity ratio in the surface water of Prydz Bay seems to be higher than the average value of other sea areas, which is attributed to the low atmospheric deposition flux of210Pb, the long residence time of210Po, or the horizontal transport of water masses with high210Po/210Pb activity ratio (Yin et al., 2004; Yang et al., 2009; Hu, 2016).

    Marine particles are carriers of biological pump. The photosynthesis of plankton in the euphotic zone drives the degradation and sedimentation of particulate organic matter in deep water. The export flux of particulate organic carbon (POC) is one of the important indicators for evaluating the efficiency of biological pump.210Po/210Pb disequilibria has been widely used in estimating the export flux of POC in recent years. These estimates of the POC flux exported from the euphotic zone are between 0.04 mmol·m?2·d?1and 12.6 mmol·m?2·d?1, showing a general decrease with the increase of latitude, which may be related to the changes in biological productivity caused by light and nutrients (Charette et al., 1999; Kim and Church, 2001; Murray et al., 2005; Stewart et al., 2007; Hu et al., 2014; Roca-Martí et al., 2016; Ma et al., 2017). Although the data around Prydz Bay is limited, the POC flux exported from the euphotic zone estimated from the210Po/210Pb disequilibria shows a pattern that the flux inside the bay is greater than that outside (Yang et al., 2009; Hu, 2016), which is generally consistent with the spatial variation of biological activities. Regarding the210Po and210Pb in Prydz Bay, some issues such as the reason why210Po is excess with respect to210Pb in surface water at some sites, and the regulation factors of scavenging and removal of210Po are still unclear and deserve more in-depth study.

    In this study, dissolved and particulate210Po and210Pb in the upper 100 m around Prydz Bay were determined to reveal the biogeochemical processes affecting the distribution of210Po and210Pb. In addition, the scavenging and removal rates of210Po are calculated based on the210Po/210Pb disequilibria. The export flux of POC was further estimated to reveal the spatial variation of biological pump. In the context of the scarcity of data on the210Po/210Pb disequilibria and POC export flux, our study is of great significance for understanding the particle dynamics and biological pump in the high-latitude Southern Ocean.

    2 Method

    2.1 Study area

    Prydz Bay is located in the eastern Antarctica and belongs to the Indian Ocean Sector of Southern Ocean. It is a semi-open embayment, shaped like a triangle. Frame Bank (near 68°S, 69°E) and Four Ladies Bank (near 70°S, 76°E) are located in the west and east of Prydz Bay, respectively. The Lambert Glacier extends from land to the bay and forms the Amery Ice Shelf in the southern boundary of Prydz Bay (Anderson, 1999). The interaction of atmosphere, ice shelf, and ocean affects the hydrological characteristics of Prydz Bay (Pu et al., 2001). The water in the upper 100 m around Prydz Bay includes two water masses, Antarctic Surface Water (AASW) and Winter Water (WW). Their temperature and salinity characteristics are as follows: AASW: 31.8 << 34.0, ?2°C << 1.5°C; WW: 34.0 << 34.4, ?2°C << ?0.5°C (Smith et al., 1984). The AASW is usually distributed at a depth less than 50 m, while the WW is located at a depth of 50–100 m below AASW(Jia, 2019). In addition, the Circumpolar Deep Water (CDW) with the characteristics of high temperature and high salinity at a depth of 100–2000 m rises from the Antarctic divergence zone, transporting nutrient-rich deep water to the surface (Su, 1987; Yabuki et al., 2006; Williams et al., 2016; Jia, 2019). Previous studies have shown that chlorophylland primary productivity in surface water in the inner bay and on the slope of Prydz Bay in summer are significantly higher than those in the open ocean (Liu et al., 2004). The concentration of POC also showed a decrease from the inner bay to the outer bay, and from the surface to the deep (Hu et al., 2001). In fact, it has been found that there is a good positive correlation between POC and chlorophyllaround Prydz Bay, indicating the biological source of POC(Hu et al., 2001). In addition, studies on the composition of phytoplankton communities around Prydz Bay have shown that diatoms are the main dominant species in summer (Cai et al., 2005). In summary, in addition to physical processes, biological activities such as POC production and remineralization also play an important role in biogeochemical cycle around Prydz Bay in summer.

    2.2 Sampling

    The seawater samples were collected from the 30th Chinese National Antarctic Research Expedition (CHINARE) onboard R/Vfrom February 26 to March 5, 2014. The sampling locations are located from 72.9°E to 76.5°E in longitude and 62.4°S to 69.3°S in latitude (Figure 1). A total of 34 samples with a depth of less than 100 m were collected at 7 stations around Prydz Bay. Among them, the station IS0 is located in the bay, and the others are outside the bay. Seawater samples with a volume of 4 to 10 dm3at different depths were collected by CTD-Rosette. The sample was immediately filtered through a mixed cellulose ester membrane with a diameter of 47 mm and a pore size of 0.4 μm to separate the dissolved and particulate phases. The filtrate was collected in a polyethylene bottle, and then by 10–20 cm3of 1:1 hydrochloric acid was added to adjust the pH to ~2. The filter membrane containing the particulate matter was placed in a sealed bag and stored frozen. Both the filtrate and the filter were brought back to the land laboratory for the determination of210Po and210Pb.

    2.3 Determination of 210Po and 210Pb

    The determination of210Po and210Pb refers to Fleer and Bacon (1984) and Yang et al. (2011). In brief, after the internal standards of209Po and stable Pb were added, the210Po and210Pb in the filtrate were enriched by Fe(OH)3co-precipitation. With the addition of shielding reagents and certain pH and temperature conditions, Po was self-deposited on a silver disc. The radioactivity of209Po and210Po was measured by an alpha spectrometer. The filter was digested with a mixed acid of nitric acid, perchloric acid and hydrofluoric acid at a high temperature. The subsequent treatment was the same as that of dissolved phase. The radioactivity of210Pb was calculated by measuring the210Po generated from210Pb over a period of time. That is, after the first self-deposition of Po, the sample solution was sealed and stored for more than one year to allow210Po growth. The radioactivity of the grown210Po was also measured with an alpha spectrometer. The reported activity concentrations of210Po and210Pb were corrected back to the sampling time according to Fleer and Bacon (1984). The reported error is ± 1σ counting error, and the propagation calculation has been performed. For convenience, the activity concentration of210Po in the dissolved and particle phases are denoted as DPo and PPo, respectively. Similarly, the activity concentration of210Pb in the dissolved and particle phases is represented by DPb and PPb, respectively. TPo and TPb are used to represent the total activity concentration of210Po and210Pb, respectively, where TPo = DPo + PPo and TPb = DPb + PPb.

    Figure 1 Sampling locations around Prydz Bay, Antarctica. The left panel shows the location of Prydz Bay. The red frame in the right panel is customized with IS0, P1-3, P1-2, P2-2, P2-A1, P2-A3, P2-A5 according to the order of station locations. The cross-sectional distribution of various parameters takes the offshore distance as horizontal axis.

    2.4 Calculation of kinetic parameters of 210Po

    Here, by ignoring the effects of advection and diffusion, a simplest one-dimensional steady-state irreversible model is used to describe the mass balance of210Po (Bacon et al., 1976). The equations are as follows:

    whereDPo,PPo,TPo,DPb,PPb, andTPbrepresent the activity concentration of DPo, PPo, TPo, DPb, PPb, and TPb (Bq·m?3), respectively. λPois the decay constant of210Po (1.828 a?1).Porepresents the atmospheric deposition rate of210Po (Bq·m?3·a?1).Porepresents the scavenging rate of dissolved210Po from dissolved to particulate phase (Bq·m?3·a?1).Porepresents the removal rate of particulate210Po due to particle sedimentation (Bq·m?3·a?1). The contribution of atmospheric deposition to210Po in seawater is very small, usually accounting for only about 2% of210Po generated from210Pb decay (Masqué et al., 2002). Therefore,Pois regarded as 0 in this study, as in many previous studies (Poet et al., 1972; Nozaki et al., 1997; Verdeny et al., 2009).

    whereDPo,PPoandTPorepresent the residence time of DPo, PPo and TPo in seawater (a), respectively.

    2.5 Other parameters

    The content of POC was determined in this study. Briefly, approximately 5 dm3of seawater samples were filtered through a pre-weighed and pre-burned (4 h at 450 °C) GF/F membrane to collect suspended particles. The sample was fumigated with hydrochloric acid for 48 h to remove inorganic carbon and washed to neutrality. After 24 h of drying, the POC was determined by an Element Analyzer (Delta V, Thermo Fisher Scientific). The standard material C6H6N2O was used to quantify POC content. The detection limit of POC was 0.1 μmol C, and the accuracy was better than 0.2%.

    Dissolved oxygen (DO) and silicate data were downloaded from the National Arctic and Antarctic Data Center of China. The DO was measured by iodometry. When the DO concentration is greater than 550 mmol·m?3, the detection limit is 5.3 mmol·m?3, and the accuracy is ±4.0 mmol·m?3.The silicate was determined by silicomolybdic blue spectrophotometry. When the silicate concentration is 4.5 mmol·m?3, the detection limit is 0.10 mmol·m?3, and the accuracy is ± 4.0%.

    3 Results

    3.1 Temperature and salinity

    The spatial variation of temperature and salinity at depths of above 200 m at our sites is shown in Figure 2. The temperature and salinity at the depth of 0–200 m are between ?1.78°C to 1.94°C and 33.18 to 34.64, with an average of ?0.27 °C and 4.10, respectively. The surface water at station IS0 located in the front of the Amery Ice Shelf has high temperature and low salinity, showing a characteristic of AASW. However, at depths greater than 100 m, the temperature decreases as the depth increases, while the salinity increases (Figure 2), which is related to the sinking and accumulation of high-density brines formed by freezing seawater in winter. A minimum temperature appears at a depth of 50–100 m at stations P1-3, P1-2, P2-2, and P2-A1, indicating the effect of cold water transporting northward. Previous studies have shown that a cold tongue with a temperature of ?1.8°C is observed between AASW and CDW at a depth of about 50 m north of 67 °S, resulting in a homogenous temperature layer with a thickness of about 50 m north of 66 °S (Pu et al., 2001). The temperature and salinity contours below 100 m at stations P2-A3 and P2-A5 protrude upward (Figure 2), indicating an upwelling of CDW.

    Figure 2 The sectional distribution of temperature, salinity, DO, silicate and POC in the upper 200 m water column.

    3.2 DO, silicate and POC

    The concentration of DO ranges from 176.01 mmol·m?3to 372.27 mmol·m?3, with an average of 305.33 mmol·m?3. The DO at station IS0 is relatively high, especially at a depth of 0–50 m (Figure 2), which may reflect the influence of primary productivity. The melting of a large number of icebergs and floating ice along the coast of Prydz Bay in summer promotes the water stratification, which is conducive to the growth of phytoplankton (Sun et al., 2012).The DO concentration at stations outside the bay is higher at the surface, and decreases with increasing depth (Figure 2). The DO at stations P2-A3 and P2-A5 decreases sharply with the increasing depth at a depth of 50–200 m, similar to changes in temperature and salinity, reflecting the effect of CDW upwelling (Figure 2).

    The concentration of silicate varies from 28.89 to 85.10 mmol·m?3, with an average of 56.40 mmol·m?3. The silicate in the bay is significantly lower than that outside the bay, which corresponds to the spatial variability of biological activities. The profiles of silicate show that the concentration generally increases with depth (Figure 2), reflecting the combined effect of absorption by phytoplankton and degradation of organic matter. The vertical changes of silicate at stations P2-A3 and P2-A5 show a significant increase at depths of 50–200 m (Figure 2), which also reflected the impact of CDW upwelling.

    The POC content varies from 0.12 mmol·m?3to 5.34 mmol·m?3, with an average of 0.78 mmol·m?3. The spatial variation of POC indicates that the POC content in the bay is significantly higher than that outside the bay (Figure 2), which is consistent with previous reports (Han et al., 2010). The profile of POC basically shows the characteristics of high surface and low depth layer, whether in the bay or outside the bay.The POC in water above 100 m varies widely, while below 100 m, the content is low and basically stable(Figure 2). The POC around Prydz Bay may mainly come from photosynthesis of plankton, which shows the differences between the inner and outer bay, and between the euphotic zone and the deep layer.

    3.3 210Po

    The activity concentration of TPo at a depth of 0–100 m around Prydz Bay ranges from 0.59 to 3.67 Bq·m?3, with an average of 1.26 Bq·m?3(Table S1). The dissolved210Po makes up a large proportion of TPo (86% ± 10%), resulting in similar spatial variability between TPo and DPo (Figure 3).

    Figure 3 The sectional distribution of activity concentration and activity ratio of dissolved, particulate and total210Po and210Pb in the upper 100 m water column.

    The maximum of DPo, PPo and TPo (1.43, 0.34, and 1.77 Bq·m?3, respectively) at station IS0 in the bay appears at 25 m (Figure 3), corresponding to the minimum of temperature (?0.06°C), salinity (33.54) and POC content (3.52 μmol·dm?3). Solar radiation, sea ice meltwater, and winter residual water together form the low temperature and low salinity in the subsurface, which weakens biological activity and particle scavenging, and retains a high210Po. In addition, the active plankton activities in surface water produce particles with high PPo, whose settling and retention forms the high PPo in the subsurface. The lowest DPo (0.75 Bq·m?3), PPo (0.01 Bq·m?3), and TPo (0.76 Bq·m?3) at a depth of 50 m correspond to the highest POC (3.92 μmol·dm?3), which may reflect the effect of enhanced biological activity on the absorption and adsorption of210Po. At stations P1-3 and P1-2, the activity concentrations of DPo and TPo in the surface layer are higher than those in the subsurface layer (Figure 3). In particular, abnormally high DPo and TPo appear in the waters above 40 m at station P2-2, which may be related to the extra input of DPo during the northward transport of the cold water mass on the shelf (Figure 3). In general, except for station P2-2, the activity concentration of DPo is the lowest in the surface layer, and increases with depth below 100 m, but the PPo at stations P2-A1, P2-A3 and P2-A5 are higher in the surface(Figure 3). The low DPo and high PPo in surface water are related to biological activity. The absorption by organisms and the adsorption by particles lead to a decrease in DPo and an increase in PPo. In the subsurface layer, the increase in DPo is accompanied by the decrease in PPo, reflecting the effect of degradation of particulate organic matter. In terms of spatial variability, the TPo and DPo in the shelf are lower than those in the slope and basin (Figure 3), indicating more active biological activities and stronger particle scavenging on the shelf.

    3.4 210Pb

    The activity concentrations of DPb and PPb vary from 1.11 to 4.00 Bq·m?3and from 0.06 to 0.67 Bq·m?3, with an average of 2.08 Bq·m?3and 0.24 Bq·m?3, respectively (Table S1). The distribution of DPb, PPb and TPb in the upper 100 m water column shows that DPb and TPb in surface water are lower, while PPb is higher (Figure 3), reflecting the effect of atmospheric deposition and particle scavenging, which is similar to those in the Weddell Sea and the Atlantic Ocean affected by the CDW (Somayajulu and Craig, 1976; Chung and Applequist, 1980). The atmospheric deposition of210Pb (Somayajulu and Craig, 1976; Baskaran, 2011) and the change in particle scavenging caused by the decrease in primary productivity with depth (Liu et al., 2004; Qiu et al., 2004; Sun et al., 2012) together form a mirror image between DPb (TPb) and PPb. Spatially, the activity concentrations of DPb and TPb in the bay are lower than those outside the bay, reflecting the stronger scavenging and removal of210Pb by biogenic particles in the bay. The high DPb and TPb and low PPb at depths of 50–100 m at stations P2-A3 and P2-A5 may be related to the upwelling of the CDW. Previous studies have found that the activity concentration of DPb in CDW is as high as 2–3 Bq·m?3(Chung and Applequist, 1980).

    4 Discussion

    4.1 Characteristics of 210Po/210Pb activity ratio

    The activity ratios of DPo/DPb (DPo/DPb)A.R.) and PPo/PPb (PPo/PPb)A.R.)are between 0.13 to 1.71 and between 0.01 to 3.23, and the average values are 0.59 and 0.82, respectively (Table S1). The activity ratio of TPo/TPb (TPo/TPb)A.R.) varies from 0.21 to 1.64, with an average of 0.58. Figure 3 shows the distribution of DPo/DPb)A.R., PPo/PPb)A.R.and TPo/TPb)A.R.in the upper 100 m water column at our stations. Although PPo/PPb)A.R.varies greatly, the average PPo/PPb)A.R.is greater than that of DPo/DPb)A.R., and the activity concentration of particulate210Po is always higher than that of particulate210Pb, indicating that biological absorption or particle adsorption preferentially remove210Po compared to210Pb.

    Among the stations in the bay, DPo/DPb)A.R., PPo/PPb)A.R.and TPo/TPb)A.R.have maximum values at a depth of 25 m, with an average of 1.15, 1.20, and 1.16, respectively, which also correspond to the maximum values of TPo. The excess of210Po at this depth indicates that organic matter degradation provides more210Po than those removed by bioabsorption and particle scavenging, thus resulting in a net input of210Po. Among the stations outside the bay, there are two cases of disequilibria between210Po and210Pb. The first case appears in the water above 20 m at station P1-2 and above 30 m at station P2-2, where DPo and TPo are excess with respect to DPb and TPb respectively, and PPo is deficient relative to PPb. The second case appears in the water above 50 m at other stations except P1-2 and P2-2, where DPo and TPo are deficient relative to DPb and TPb, respectively, indicating that210Po has a stronger affinity for particles.

    The spatial variability shows that DPo/DPb)A.R.and TPo/TPb)A.R.in the waters above 50 m at each station inside and outside the bay are both less than 1, indicating that the deficits of DPo and TPo outside the bay are greater than those in the bay. Interestingly, the PPo/PPb)A.R.in the bay is different from that outside the bay, where it is less than 1.0 in the bay and greater than 1.0 outside the bay. The difference in the spatial distribution of DPo/DPb)A.R.and PPo/PPb)A.R.may be related to the spatial variation of biological activities. In the upper water in the bay, biological activities are intense, resulting in abundant biogenic particles, which effectively remove210Po and210Pb during sedimentation process. However, biological activities outside the bay is relatively weak, resulting the particles’ removal of210Po prior to210Pb. In the 50–100 m depth interval at stations P2-A3 and P2-A5 affected by CDW upwelling, DPo/DPb)A.R.and TPo/TPb)A.R.are significantly less than 1.0, while DPo/DPb)A.Ris significantly greater than 1.0.

    4.2 Abnormal excess of 210Po in surface water

    In the surface water of stations P1-2 and P2-2, DPo/DPb)A.R.are 1.55 and 1.43, and TPo/TPb)A.R.are 1.31 and 1.28, respectively, showing that210Po significantly exceeds210Pb, which is rare in previous studies. The possible reasons for the abnormal excess of210Po include four aspects, namely, the effect of upwelling, the input of sea ice meltwater, the degradation of organic matter, and the horizontal transport of water masses with high DPo/DPb)A.R.and TPo/TPb)A.R.. Below we will discuss these four possibilities one by one.

    The upwelling of CDW promotes the growth of phytoplankton in surface water, leading to an increase in sedimentation of organic debris. These organic debris preferentially release210Po to the deep water through the degradation of organic matter, and then enter the surface layer through the upwelling, which leads to a continuous accumulation of210Po relative to210Pb in surface water (Kadko, 1993). Although CDW upwelling was observed near 65°S outside Prydz Bay during the 9th, 15th, and 18th CHINAREs, distribution of temperature, salinity and nutrients in this voyage showed that the upwelling of CDW at stations P1-2 and P2-2 did not reach the surface(Figure 2). In addition,210Po/210Pb)A.R.in the subsurface water at stations P1-2 and P2-2 is not higher than the surface. Therefore, the upwelling of CDW may not be a main reason for the excess of210Po in the surface water.

    The input of sea ice meltwater will lead to an increase of210Po/210Pb)A.R.in seawater, because210Pb deposited and accumulated on sea ices by atmospheric deposition undergoes radioactive decay to produce210Po. However, the reported210Po/210Pb)A.R.in sea ice meltwater falls between 0.4 and 1.0 (Masqué et al., 2007; Roca-Martí et al., 2016; RocaMartí et al., 2018), which cannot support the210Po excess (i.e.210Po/210Pb)A.R.>1) observed in this study. Moreover, sea ice meltwater often promotes the growth of plankton by enhancing stratification, leading to preferential scavenging and removal of210Po, which in turn reduces210Po/210Pb)A.R.(RocaMartí et al., 2018). Therefore, the input of sea ice meltwater may not be the reason for the obvious excess of210Po over210Pb in the surface water we observed here.

    The degradation of particulate organic matter preferentially releases210Po into seawater, resulting in an excess of210Po with respect to210Pb (Thomson and Turekian, 1976; Cochran et al., 1983; Shimmield et al., 1995; Nozaki et al., 1997; Sarin et al., 1999; Yang et al., 2009). Considering that the bioabsorption or particle scavenging in surface water is usually stronger than the degradation of organic matter, the obvious excess of210Po in surface water is also unlikely to be mainly formed by the degradation of organic matter.

    After excluding the above three possibilities, the remaining possible cause is the horizontal transport of water masses with high DPo/DPb)A.R.and TPo/TPb)A.R.. Note that station P2-2 is located at approximately the same latitude (65.48°S) east of station P1-2. If water masses with high210Po/210Pb)A.R.are transported horizontally from the east or west, it can explain the obvious excess of210Po in the surface water at stations P1-2 and P2-2. However, due to the limitation of spatial resolution of our data, we are unable to determine the source of the water masses with high210Po/210Pb)A.R., which needs further research in the future.

    4.3 Particle dynamics of 210Po

    Based on the one-dimensional steady-state model proposed by Bacon et al. (1976), the scavenging rate, removal rate and residence time of210Po around Prydz Bay were calculated, as shown in Figure 4. Note that at 25 m at station IS0, 0 m at station P1-2, and 0 m and 25 m at station P2-2,210Po is excess with respect to210Pb, which cannot be calculated by above model. The scavenging rate, removal rate and residence time of210Po in these layers are treated as zero in Figure 4. Our calculation shows that the scavenging rate and removal rate of210Po in the waters around Prydz Bay are between 0.40 to 5.94 Bq·m?3·a?1and 0.31 to 5.89 Bq·m?3·a?1, with an average of 2.20 and 2.35 Bq·m?3·a?1, respectively. The residence time of DPo and PPo varies from 0.09 a to 2.36 a and from 0 to 1.71 a, respectively, with an average of 0.71 a and 0.13 a, respectively. The residence time of TPo ranges from 0.14 a to 5.47 a, with an average of 0.76 a. There is a good linear positive correlation between the scavenging rate and the removal rate:Po= 0.89×Po+ 0.38 (2= 0.95,< 0.0001) (Figure 5), indicating that the removal of particulate210Po is mainly controlled by the scavenging of dissolved210Po, which is similar to that in the North Atlantic (Bacon et al., 1976). ThePoandPoare lower in the entire 100 m water column at station IS0 at the front of the ice shelf, despite the higher phytoplankton biomass and primary productivity (Figure 4). The reason may be the lagging response of210Po to particle removal under the influence of sea ice and residual water in winter. In terms of spatial variability,PoandPogradually increase from the inside to the outside of the bay, and increase with the increasing depth (Figure 4), indicating that the particles continuously scavenge and remove210Po during their settling process. The highPoandPocorrespond spatially to low DPo and TPo, and high DPb, TPb and PPo/PPb)A.R., reflecting the upwelling of CDW not only brings dissolved210Pb, but also strengthens the scavenging of210Po. The residence times of DPo and TPo are longer in the bay, but shorter in the areas affected by the CDW upwelling (Figure 4). Except for the surface water at station P2-A3, the residence time of PPo in the entire water column is relatively short. The longer residence time of PPo in surface water of station P2-A3 corresponds to a significant excess of PPo relative to PPb, which may be the result of high PPo/PPb)A.R.particle input.

    Figure 4 The sectional distribution of scavenging and removal rates of210Po, and residence times of DPo, PPo and TPo in the upper 100 m water column.

    Figure 5 The relationship between scavenging rate and removal rate of210Po around Prydz Bay.

    In order to reveal the factors affecting the particle dynamics of210Po around Prydz Bay, the relationship betweenPo,Poand temperature, salinity, dissolved oxygen (DO), silicate, and POC was explored, as shown in Figure 6 and Figure 7. In summary,PoandPoare not significantly correlated with temperature and silicate, but are significantly correlated with salinity, dissolved oxygen and POC. The positive correlation betweenPo,Poand salinity shows that as the salinity increases, the scavenging and removal rates of210Po increase. Since the salinity generally increases with the increasing depth, this relationship indicates that under the combined action of biological absorption and particle adsorption, the scavenging and removal rates of210Po increase with depth in the upper 100 m water column. There is a good negative linear relationship betweenPo,Poand dissolved oxygen. The fitting equations are as follows:Po= ?0.019 × DO + 8.18 (= 0.0003), andPo= ?0.016 × DO + 7.37 (= 0.001). In addition, there is a negative correlation betweenPo,Poand POC, and the fitting equations are as follows:Po= 0.63/POC + 0.94 (= 0.001), andPo= 0.52/POC + 1.32 (= 0.005). These relationships amongPo,Po, DO, and POC may reflect the effect of biological processes on the geochemical fractionation between210Po and210Pb. When biological photosynthesis is stronger than the degradation of organic matter, biogenic organic matter and DO increase, which may lead to a decrease in the fractionation between210Po and210Pb and a decrease inPoandPo. On the contrary, when the degradation of organic matter is stronger than biological photosynthesis, biogenic organic matter and DO decrease, which may lead to an increase in the fractionation between210Po and210Pb and an increase inPoandPo(Chen et al., 2021).

    Figure 6 The relationship between scavenging rate of210Po and (a) temperature, (b) salinity, (c) DO, (d) silicate, and (e) POC around Prydz Bay.

    4.4 Estimation of POC export flux

    The210Po/210Pb disequilibrium has been widely used to estimate the export flux of POC in the ocean (Kim and Church, 2001; Murray et al., 2005; Hu et al., 2014; Roca-Martí et al., 2016). Here, an empirical method given by Shimmield et al. (1995) is used to calculate the POC export flux around Prydz Bay based on the removal flux of210Po. The calculation formula is as follows:

    ,(9)

    ,(10)

    wherePOCrepresents POC export flux (mmol·m?2·d?1).POCandPPorepresent the inventory of POC (mmol·m?2) and PPo (Bq·m?2) in the water column, respectively.Po,removalrepresents the removal flux of particulate210Po from particle settling (Bq·m?2·d?1).TPoandTPbrepresent the inventory of TPo and TPb (Bq·m?2) in the water column, respectively. As for the selection of depth interface, considering that previous studies in the Bellingshausen Sea (Shimmield et al., 1995) and the Antarctic Circumpolar Current (Friedrich and van der Loeff, 2002) used 100 m as the interface, for the convenience of comparison, the POC export flux at 100 m is reported here. In addition, it has been pointed out in the above discussion that the210Po in the surface water at stations P1-2 and P2-2 may be affected by horizontal transport, which does not meet the conditions of the adopted model, so these two stations are not included in the calculation.

    Figure 7 The relationship between removal rate of210Po and (a) temperature, (b) salinity, (c) DO, (d) silicate, and (e) POC around Prydz Bay.

    Our calculations show that the export flux of POC at 100 m around Prydz Bay varies from 1.8 mmol·m?2·d?1to 4.4 mmol·m?2·d?1, with an average of 2.9 mmol·m?2·d?1. As shown in Figure 8, the highest POC export occurs at station IS0 at the front of the ice shelf, which is consistent with the highest phytoplankton biomass and primary productivity. The lowest POC export appears at station P2-A3 outside the bay. Note that the POC export flux in the continental slope (station P1-3) is lower than that in the open ocean area (stations P2-A1 and P2-A3), which indicates that primary productivity in the open ocean may be higher than that in the continental slope during the survey period. Previous studies have shown that although the phytoplankton biomass in the continental slope outside Prydz Bay is higher than that in the open ocean, the open ocean has a higher primary productivity due to a thicker euphotic zone (Liu et al., 2001). Therefore, the spatial variability of POC export flux may be closely related to primary production around Prydz Bay.

    Figure 8 The POC export flux at each station around Prydz Bay estimated from210Po/210Pb disequilibria.

    The export flux of POC around Prydz Bay has been estimated through210Po/210Pb disequilibria prior to this study. The comparison shows that the estimated POC export flux here is consistent with the report in the open ocean outside Prydz Bay (2.3 mmol·m?2·d?1, one site at 64.0 °S, 73.0 °E, Yang et al., 2009), but slightly lower than the reports containing more shelf stations (4.2–9.0 mmol·m?2·d?1, avg. 6.9 mmol·m?2·d?1, 68.0°S–64.0°S, 76.0°E–79.0°E) (Hu et al., 2021). Table S2 shows that the estimated POC export fluxes are lower in the high-latitude regions and increase with the decreasing latitude due to the changes in biological productivity induced by light or nutrients. It is worth noting that the export flux of POC based on210Po/210Pb disequilibria is significantly lower than that based on234Th/238U disequilibria (52.6– 185.6 mmol·m?2·d?1, avg. 104.7 mmol·m?2·d?1) (He et al., 2008), which may reflect the difference in time scale and geochemical behavior between210Po and234Th. Ceballos- Romero et al. (2016) compared the POC export fluxes obtained by the234Th/238U method,210Po/210Pb method and the sediment trap in the North Atlantic, and found that the results obtained by different methods differ by orders of magnitude except for the consistency in the late bloom, which is attributed to the different time scales of nuclides and the different growth stages of phytoplankton. Regarding the three methods for estimating POC export flux, sediment trap is a most direct tool to obtain POC export, which is susceptible to several factors such as hydrodynamics, biological predation, and particle dissolution (Buesseler et al., 2007). The234Th/238U disequilibria is more suitable for flux estimation under short time scale (seasonal) and high POC export. Studies have shown that in the case of low POC export flux, the results obtained by234Th/238U disequilibria are often higher than those by sediment trap (Stewart et al., 2010).210Po (1/2= 138.4 d) has a longer half-life than234Th (1/2= 24.1 d) and is used to track the POC export on a seasonal or interannual time scale. In addition,210Po is more biologically active than234Th, which means that210Po may be more suitable for tracking POC export (Le Moigne et al., 2013).

    5 Conclusions

    The spatial variation of210Po and210Pb and their disequilibria in the upper 100 m water column around Prydz Bay in the summer of 2014 were studied. We found the abnormal phenomenon in the surface water at some stations that210Po is excess with respect to210Pb, which probably reflects the impact of horizontal transport of water masses with highDPo/DPb)A.R.and high TPo/TPb)A.R.The scavenging and removal rates of210Po vary from 0.40 to 5.94 Bq·m?3·a?1(average 2.2 Bq·m?3·a?1) and from 0.31 to 5.89 Bq·m?3·a?1(average 2.35 Bq·m?3·a?1), respectively. There is a good linear positive correlation between the two rates, indicating that the removal of particulate210Po is mainly controlled by the scavenging process from the dissolved to particulate phase. The scavenging and removal rates of210Po is positively correlated with salinity and negatively correlated with dissolved oxygen and POC, indicating that biological activities have caused biogeochemical fractionation between210Po and210Pb. Based on the210Po/210Pb disequilibria, the POC export flux at 100 m around Prydz Bay is estimated to vary from 1.8 to 4.4 mmol·m?2·d?1, showing a feature of being higher in the bay than outside the bay.

    Acknowledgements We would like to thank the captains and crews of the R/Vfor their assistance in sampling. Data on temperature, salinity, dissolved oxygen and silicate are provided by National Arctic and Antarctic Data Center of China.This work was financially supported by National Polar Special Program “Impact and Response of Antarctic Seas to Climate Change” (Grant nos. IRASCC 01-01-02C, 02-01-01), was supported by National Natural Science Foundation of China (Grant no. 41721005), and the program funded by China Ocean Mineral Resources R & D Association (Grant no. DY135-13-E2-03). The authors also thank two anonymous reviewers, and Guest Editor Dr. Jianfeng He for their constructive comments.

    Anderson J B. 1999. Antarctic marine geology. Cambridge: Cambridge University Press, doi: 10.1016/S0025-3227(00)00098-0.

    Bacon M P, Spencer D W, Brewer P G. 1976.210Pb/226Ra and210Po/210Pb disequilibria in seawater and suspended particulate matter. Earth Planet Sci Lett, 32(2): 277-296, doi:10.1016/0012-821X(76)90068-6.

    Baskaran M. 2011. Po-210 and Pb-210 as atmospheric tracers and global atmospheric Pb-210 fallout: a review. J Environ Radioact, 102(5): 500-513, doi:10.1016/j.jenvrad.2010.10.007.

    Buesseler K O, Antia A N, Chen M, et al. 2007. An assessment of the use of sediment traps for estimating upper ocean particle fluxes. J Mar Res, 65(3): 345-416, doi:10.1357/002224007781567621.

    Cai Y M, Ning X R, Zhu G H, et al. 2005. Size fractionated biomass and productivity of phytoplankton and new production in the Prydz Bay and the adjacent Indian sector of the Southern Ocean during the austral summer of 1998/1999. Acta Oceanol Sin, 27(4): 135-147 (in Chinese with English abstract).

    Ceballos-Romero E, Le Moigne F A C, Henson S, et al. 2016. Influence of bloom dynamics on particle export efficiency in the North Atlantic: a comparative study of radioanalytical techniques and sediment traps. Mar Chem, 186: 198-210, doi:10.1016/j.marchem.2016.10.001.

    Charette M A, Moran S B, Bishop J K B. 1999.234Th as a tracer of particulate organic carbon export in the subarctic northeast Pacific Ocean. Deep Sea Res Part II Top Stud Oceanogr, 46(11-12): 2833-2861, doi:10.1016/S0967-0645(99)00085-5.

    Chen M Y, Chen M, Zheng M F, et al. 2021.210Po/210Pb disequilibria influenced by production and remineralization of particulate organic matter around Prydz Bay, Antarctica. Deep Sea Res Part II Top Stud Oceanogr, 191-192: 104961, doi:10.1016/j.dsr2.2021.104961.

    Chung Y, Applequist M D. 1980.226Ra and210Pb in the Weddell Sea. Earth Planet Sci Lett, 49(2): 401-410, doi:10.1016/0012-821X(80)90082-5.

    Cochran J K, Bacon M P, Krishnaswami S, et al. 1983.210Po and210Pb distributions in the central and eastern Indian Ocean. Earth Planet Sci Lett, 65(2): 433-452, doi:10.1016/0012-821X(83)90180-2.

    Fleer A P, Bacon M P. 1984. Determination of210Pb and210Po in seawater and marine particulate matter. Nucl Instrum Methods Phys Res, 223(2-3): 243-249, doi:10.1016/0167-5087(84)90655-0.

    Friedrich J, van der Loeff M M R. 2002. A two-tracer (210Po-234Th) approach to distinguish organic carbon and biogenic silica export flux in the Antarctic Circumpolar Current. Deep Sea Res Part I Oceanogr Res Pap, 49(1): 101-120, doi:10.1016/S0967-0637(01)00045-0.

    Han Z B, Hu C Y, Yu W, et al. 2010. Decmposition of organic carbon and inorganic carbon beneath euphotic zone in Prydz Bay, Antarctica. Chin J Polar Res, 22(3): 254-261 (in Chinese with English abstract).

    He J H, Ma H, Chen L Q, et al. 2008. The investigation on particulate organic carbon fluxes with disequilibria between thorium-234 and uranium-238 in the Prydz Bay, the Southern Ocean. Acta Oceanol Sin, 27(2): 21-29 (in Chinese with English abstract).

    Hu C Y, Zhang H S, Pan J M. 2001. The biogeochemistry of carbon cycle in summer of the Prydz Bay, Antarctica Ⅱ: characteristics of POC distribution. Chin J Polar Res, 13(3): 195-204 (in Chinese with English abstract).

    Hu H N. 2016. Distribution of210Po and210Pb in the Prydz Bay and its adjacent sea areas and their estimate of particulate organic matter export. M.S. thesis, Xiamen, China: Xiamen University (in Chinese).

    Hu H N, Liu X, Ren C Y, et al. 2021.210Po/210Pb disequilibria and its estimate of particulate organic carbon export around Prydz Bay, Antarctica. Front Mar Sci, 8: 701014, doi:10.3389/fmars.2021. 701014.

    Hu W J, Chen M, Yang W F, et al. 2014. Enhanced particle scavenging in deep water of the Aleutian Basin revealed by210Po-210Pb disequilibria. J Geophys Res Oceans, 119(6): 3235-3248, doi:10.1002/2014JC009 819.

    Jia R M. 2019. Spatial and temporal variations of freshwater components around the Prydz Bay, Antarctica and its implication for marine processes. Ph. D. thesis, Xiamen, China: Xiamen University (in Chinese).

    Kadko D. 1993. Excess210Po and nutrient recycling within the California coastal transition zone. J Geophys Res, 98(C1): 857-864, doi:10.1029/92JC01932.

    Kim G, Church T M. 2001. Seasonal biogeochemical fluxes of234Th and210Po in the upper Sargasso Sea: influence from atmospheric iron deposition. Global Biogeochem Cycles, 15(3): 651-661, doi:10.1029/ 2000GB001313.

    Le Moigne F A C, Villa-Alfageme M, Sanders R J, et al. 2013. Export of organic carbon and biominerals derived from234Th and210Po at the Porcupine Abyssal Plain. Deep Sea Res Part I Oceanogr Res Pap, 72: 88-101, doi: 10.1016/j.dsr.2012.10.010.

    Liu C G, Ning X R, Sun J, et al. 2004. Size structure of standing stock and productivity and new production of phytoplankton in the Prydz Bay and the adjacent Indian sector of the Southern Ocean during the austral summer of 2001/2002. Acta Oceanol Sin, 26(6): 107-117 (in Chinese with English abstract).

    Liu Z L, Cai Y M, Ning X R, et al. 2001. Primary productivity and standing stock of phytoplankton in the Prydz Bay and the adjacent northern sea area during the austral summer of 1999/2000. Chin J Polar Res, 13(1): 1-12 (in Chinese with English abstract).

    Ma H Y, Yang W F, Zhang L H, et al. 2017. Utilizing210Po deficit to constrain particle dynamics in mesopelagic water, western South China Sea. Geochem Geophys Geosyst, 18(4): 1594-1607, doi:10.1002/2017GC006899.

    Masqué P, Cochran J K, Hirschberg D J, et al. 2007. Radionuclides in Arctic sea ice: tracers of sources, fates and ice transit time scales. Deep Sea Res Part I Oceanogr Res Pap, 54(8): 1289-1310, doi:10.1016/j.dsr.2007.04.016.

    Masqué P, Sanchez-Cabeza J A, Bruach J M, et al. 2002. Balance and residence times of210Pb and210Po in surface waters of the northwestern Mediterranean Sea. Cont Shelf Res, 22(15): 2127-2146, doi:10.1016/S0278-4343(02)00074-2.

    Murray J W, Paul B, Dunne J P, et al. 2005.234Th,210Pb,210Po and stable Pb in the central equatorial Pacific: tracers for particle cycling. Deep Sea Res Part I Oceanogr Res Pap, 52(11): 2109-2139, doi:10.1016/j.dsr.2005.06.016.

    Nozaki Y, Dobashi F, Kato Y, et al. 1998. Distribution of Ra isotopes and the210Pb and210Po balance in surface seawaters of the mid Northern Hemisphere. Deep Sea Res Part I Oceanogr Res Pap, 45(8): 1263-1284, doi:10.1016/S0967-0637(98)00016-8.

    Nozaki Y, Zhang J, Takeda A. 1997.210Pb and210Po in the equatorial Pacific and the Bering Sea: the effects of biological productivity and boundary scavenging. Deep Sea Res Part II Top Stud Oceanogr, 44(9-10): 2203-2220, doi:10.1016/S0967-0645(97)00024-6.

    Poet S E, Moore H E, Martell E A. 1972. Lead 210, bismuth 210, and polonium 210 in the atmosphere: accurate ratio measurement and application to aerosol residence time determination. J Geophys Res, 77(33): 6515-6527, doi:10.1029/JC077i033p06515.

    Pu S Z, Dong Z Q, Hu X M, et al. 2001. Upper layer waters and their northward extension from Prydz Bay in summer. Chin J Polar Sci, 12(2): 89-98 (in Chinese with English abstract).

    Qiu Y S, Huang Y P, Liu G S, et al. 2004. Spatial and temporal variations of primary productivity in Prydz Bay and its adjacent sea area, Antarctica. J Xiamen Univ Nat Sci, 43(5): 676-681 (in Chinese with English abstract).

    Radakovitch O, Cherry R D, Heyraud M, et al. 1998. Unusual210Po/210Pb ratios in the surface water of the Gulf of Lions. Oceanol Acta, 21(3): 459-468, doi:10.1016/S0399-1784(98)80030-3.

    Roca-Martí M, Puigcorbé V, Friedrich J, et al. 2018. Distribution of210Pb and210Po in the Arctic water column during the 2007 sea-ice minimum: particle export in the ice-covered basins. Deep Sea Res Part I Oceanogr Res Pap, 142: 94-106, doi:10.1016/j.dsr.2018.09.011.

    Roca-Martí M, Puigcorbé V, van der Loeff M M R, et al. 2016. Carbon export fluxes and export efficiency in the central Arctic during the record sea-ice minimum in 2012: a joint234Th/238U and210Po/210Pb study. J Geophys Res Oceans, 121(7): 5030-5049, doi:10.1002/ 2016jc011816.

    Sarin M M, Kim G, Church T M. 1999.210Po and210Pb in the south-equatorial Atlantic: distribution and disequilibrium in the upper 500 m. Deep Sea Res Part II Top Stud Oceanogr, 46(5): 907-917, doi:10.1016/S0967-0645(99)00008-9.

    Shimmield G B, Ritchie G D, Fileman T W. 1995. The impact of marginal ice zone processes on the distribution of210Pb,210Po and234Th and implications for new production in the Bellingshausen Sea, Antarctica. Deep Sea Res Part II Top Stud Oceanogr, 42(4-5): 1313-1335, doi:10.1016/0967-0645(95)00071-W.

    Smith N R, Dong Z Q, Kerry K R, et al. 1984. Water masses and circulation in the region of Prydz Bay, Antarctica. Deep Sea Res A Oceanogr Res Pap, 31(9): 1121-1147, doi:10.1016/0198-0149(84) 90016-5.

    Somayajulu B L K, Craig H. 1976. Particulate and soluble210Pb activities in the deep sea. Earth Planet Sci Lett, 32(2): 268-276, doi:10. 1016/0012-821X(76)90067-4.

    Stewart G, Cochran J K, Xue J H, et al. 2007. Exploring the connection between210Po and organic matter in the northwestern Mediterranean. Deep Sea Res Part I Oceanogr Res Pap, 54(3): 415-427, doi:10.1016/j. dsr.2006.12.006.

    Stewart G M, Moran S B, Lomas M W. 2010. Seasonal POC fluxes at BATS estimated from210Po deficits. Deep Sea Res Part I Oceanogr Res Pap, 57(1): 113-124, doi:10.1016/j.dsr.2009.09.007.

    Su Y F. 1987. Upwelling of the deep water in the Prydz Bay and offshore in the Antarctic during the austral summer. Trans Oceanol Limnol, (2): 17-24, doi:10.13984/j.cnki.cn37-1141.1987.02.003 (in Chinese with English abstract).

    Sun W P, Hu C Y, Han Z B, et al. 2012. Distribution of nutrients and chlin Prydz Bay during the austral summer of 2011. Chin J Polar Res, 24(2): 178-186 (in Chinese with English abstract).

    Thomson J, Turekian K K. 1976.210Po and210Pb distributions in ocean water profiles from the Eastern South Pacific. Earth Planet Sci Lett, 32(2): 297-303, doi:10.1016/0012-821X(76)90069-8.

    van der Loeff M M R, Friedrich J, Bathmann U V. 1997. Carbon export during the Spring Bloom at the Antarctic Polar Front, determined with the natural tracer 234Th. Deep Sea Res Part II Top Stud Oceanogr, 44(1-2): 457-478, doi:10.1016/S0967-0645(96)00067-7.

    Verdeny E, Masqué P, Garcia-Orellana J, et al. 2009. POC export from ocean surface waters by means of234Th/238U and210Po/210Pb disequilibria: a review of the use of two radiotracer pairs. Deep Sea Res Part II Top Stud Oceanogr, 56(18): 1502-1518, doi:10.1016/j. dsr2.2008.12.018.

    Williams G, Herraiz-Borreguero L, Roquet F, et al. 2016. The suppression of Antarctic bottom water formation by melting ice shelves in Prydz Bay. Nat Commun, 7: 1-9, doi: 10.1038/ncomms12577.

    Yabuki T, Suga T, Hanawa K, et al. 2006. Possible source of the Antarctic bottom water in the Prydz Bay Region. J Oceanogr, 62(5): 649-655, doi:10.1007/s10872-006-0083-1.

    Yang W F, Huang Y P, Chen M, et al. 2006. Disequilibria between210Po and210Pb in surface waters of the southern South China Sea and their implications. Sci China Ser D, 49(1): 103-112, doi:10.1007/s11430- 004-5233-y.

    Yang W F, Huang Y P, Chen M, et al. 2009. Export and remineralization of POM in the Southern Ocean and the South China Sea estimated from210Po/210Pb disequilibria. Chin Sci Bull, 54(12): 2118-2123, doi:10.1007/s11434-009-0043-4.

    Yang W F, Huang Y P, Chen M, et al. 2011. Carbon and nitrogen cycling in the Zhubi coral reef lagoon of the South China Sea as revealed by210Po and210Pb. Mar Pollut Bull, 62(5): 905-911, doi:10.1016/j. marpolbul.2011.02.058.

    Yin M D, Zeng W Y, Wu S Y, et al. 2004. Distribution of uranium-series isotopes in the Prydz Bay, Antarctica. Chin J Polar Res, 16(1): 11-21 (in Chinese with English abstract).

    Table S1 The activity concentrations and activity ratio of dissolved, particulate and total210Po and210Pb around Prydz Bay

    StationsLayerDPo/(Bq·m?3)PPo/(Bq·m?3)TPo/(Bq·m?3)DPb/(Bq·m?3)PPb/(Bq·m?3)TPb/(Bq·m?3)DPo/DPb)A.R.PPo/PPb)A.R.TPo/TPb)A.R. IS000.97±0.090.11±0.031.08±0.091.28±0.100.25±0.031.53±0.100.75±0.090.45±0.150.70±0.08 251.43±0.090.34±0.051.77±0.101.24±0.090.28±0.041.53±0.101.15±0.111.20±0.241.16±0.10 500.75±0.050.01±0.040.76±0.071.11±0.080.54±0.071.65±0.110.68±0.070.02±0.080.46±0.05 750.94±0.060.12±0.031.06±0.071.16±0.080.20±0.031.36±0.090.81±0.080.60±0.160.78±0.07 1000.84±0.060.03±0.020.86±0.061.35±0.100.21±0.031.56±0.100.62±0.060.13±0.110.55±0.06 P1-301.95±0.090.21±0.062.17±0.112.47±0.190.67±0.093.14±0.210.79±0.070.32±0.090.69±0.06 500.88±0.050.17±0.021.06±0.051.94±0.140.13±0.022.07±0.140.46±0.041.32±0.260.51±0.04 750.79±0.050.43±0.051.21±0.072.35±0.180.13±0.022.48±0.180.33±0.033.23±0.590.49±0.04 1001.05±0.050.08±0.021.13±0.062.11±0.140.10±0.022.21±0.140.49±0.040.83±0.250.51±0.04 P1-201.94±0.070.25±0.062.19±0.091.25±0.100.43±0.061.68±0.111.55±0.130.59±0.151.31±0.10 250.70±0.070.25±0.060.94±0.091.23±0.090.55±0.071.78±0.120.56±0.070.45±0.120.53±0.06 500.79±0.050.23±0.041.01±0.062.17±0.160.26±0.042.42±0.170.36±0.040.88±0.190.42±0.04 751.28±0.060.01±0.021.29±0.072.19±0.170.26±0.032.46±0.170.58±0.050.05±0.090.52±0.05 1000.78±0.050.01±0.020.79±0.052.12±0.160.21±0.032.33±0.170.37±0.040.05±0.090.34±0.03 P2-202.30±0.110.07±0.022.37±0.111.61±0.130.25±0.031.86±0.131.43±0.130.27±0.091.28±0.11 253.58±0.170.09±0.023.67±0.182.09±0.150.15±0.022.24±0.151.71±0.150.62±0.171.64±0.14 500.88±0.060.12±0.021.00±0.072.22±0.160.16±0.022.38±0.160.40±0.040.74±0.170.42±0.04 750.57±0.040.09±0.020.65±0.052.06±0.160.06±0.012.13±0.160.28±0.031.34±0.360.31±0.03 1000.75±0.060.11±0.020.86±0.062.06±0.160.07±0.012.13±0.160.36±0.041.55±0.420.40±0.04 P2-A100.60±0.070.24±0.050.84±0.081.23±0.100.34±0.041.57±0.110.49±0.070.71±0.180.54±0.07 250.76±0.080.20±0.040.96±0.091.73±0.130.28±0.042.01±0.140.44±0.060.73±0.190.48±0.05 500.83±0.090.07±0.020.90±0.092.39±0.200.12±0.022.51±0.200.35±0.050.58±0.190.36±0.05 751.88±0.140.03±0.021.91±0.142.48±0.180.15±0.022.63±0.180.76±0.080.22±0.160.73±0.07 1000.56±0.080.03±0.020.59±0.082.49±0.180.12±0.022.61±0.180.23±0.030.22±0.130.22±0.03 P2-A301.16±0.080.53±0.051.68±0.091.53±0.110.32±0.041.85±0.120.75±0.081.67±0.260.91±0.08 251.18±0.080.22±0.041.40±0.091.67±0.130.39±0.052.06±0.140.71±0.070.56±0.120.68±0.06 501.00±0.090.06±0.021.05±0.092.75±0.210.19±0.032.94±0.210.36±0.040.30±0.100.36±0.04 750.52±0.060.15±0.020.67±0.072.41±0.180.13±0.022.54±0.180.22±0.031.19±0.240.27±0.03 1002.25±0.130.10±0.012.35±0.133.29±0.260.07±0.013.36±0.260.68±0.071.43±0.350.70±0.07 P2-A500.55±0.060.17±0.020.72±0.072.57±0.200.28±0.042.85±0.210.21±0.030.60±0.120.25±0.03 251.05±0.080.12±0.041.17±0.092.07±0.160.34±0.042.41±0.170.51±0.060.37±0.120.49±0.05 500.93±0.080.15±0.021.08±0.082.70±0.210.28±0.042.98±0.210.34±0.040.55±0.110.36±0.04 750.45±0.060.33±0.040.78±0.073.33±0.250.13±0.023.46±0.250.13±0.022.65±0.610.23±0.03 1000.75±0.080.10±0.020.85±0.084.00±0.300.07±0.014.07±0.300.19±0.021.42±0.420.21±0.03

    Table S2 The export flux of POC from the euphotic zone based on210Po/210Pb and234Th/238U disequilibria

    RegionlatitudeExport interface/m210Po/210Pb)A.R.*234Th–POC flux*/ (mmol·m?2·d?1)210Po–POC flux*/ (mmol·m?2·d?1)References Aleutian Basin64 oN–65oN1000.05–1.32(avg. 0.65)nd0.6–1.7(avg. 1.1)Hu et al. (2014) Mediterranean Sea40 oN–45oN200nd3.8–17.5(avg. 9.7)4.4–7.0(avg. 5.7)Stewart et al. (2007a, 2007b) Sargasso Sea30 oN–32oN1500.21–0.94(avg. 0.53)0.3–24.5(avg. 9.3)2.4–9.1(avg. 4.7)Kim and Church (2001) Antarctic Circumpolar Current47oS–57oS1000.47–0.99(avg. 0.66)3.9–38.4(avg. 16.7)4.8–17.1(avg. 10.8)van der Loeff et al. (1997); Friedrich and van der Loeff (2002) Bellingshausen Sea67oS1000.58–1.12(avg. 0.82)212.2Shimmield et al. (1995) Prydz Bay64oS100ndnd2.3Yang et al. (2009) Prydz Bay64oS–68oS100/2000.36–1.20(avg. 0.65)nd4.2–9.0(avg. 6.9)Hu et al. (2021) Prydz Bay62oS–69oS1000.21–1.16(avg. 0.53)nd1.8–4.4(avg. 2.9)This study

    Notes: * Value in parentheses indicate average value; nd represents no data

    5 September 2021;

    15 November 2021;

    30 March 2022

    10.13679/j.advps.2021.0045

    , ORCID: 0000-0003-0369-694X, E-mail:mchen@xmu.edu.cn

    : Chen M Y, Chen M, Zheng M F, et al.Particle dynamics revealed by210Po/210Pb disequilibria around Prydz Bay , the Southern Ocean in summer. Adv Polar Sci, 2022, 33(1):71-85,doi:10.13679/j.advps.2021.0045

    久久久精品免费免费高清| 欧美国产精品va在线观看不卡| 热re99久久精品国产66热6| 免费黄色在线免费观看| 丰满饥渴人妻一区二区三| 日韩 亚洲 欧美在线| av又黄又爽大尺度在线免费看| 性高湖久久久久久久久免费观看| 在线 av 中文字幕| 国产熟女午夜一区二区三区| 精品人妻熟女毛片av久久网站| 久久国内精品自在自线图片| 亚洲精品日本国产第一区| 男女免费视频国产| 亚洲欧美成人精品一区二区| 久久这里有精品视频免费| 在线观看人妻少妇| 亚洲伊人色综图| 欧美 日韩 精品 国产| 在线免费观看不下载黄p国产| 欧美日韩精品网址| 97人妻天天添夜夜摸| 下体分泌物呈黄色| 国产成人aa在线观看| 美女福利国产在线| 26uuu在线亚洲综合色| 国产探花极品一区二区| 免费黄网站久久成人精品| 成人影院久久| 国产成人av激情在线播放| 免费久久久久久久精品成人欧美视频| 可以免费在线观看a视频的电影网站 | 99香蕉大伊视频| 观看av在线不卡| av免费观看日本| 成人手机av| 欧美日韩视频精品一区| 美女主播在线视频| 欧美日韩视频精品一区| 国产一区二区激情短视频 | 大话2 男鬼变身卡| 亚洲精品久久成人aⅴ小说| av不卡在线播放| 精品久久久精品久久久| 精品久久久精品久久久| 国产免费一区二区三区四区乱码| 青青草视频在线视频观看| 久久午夜综合久久蜜桃| 男女高潮啪啪啪动态图| 高清av免费在线| 亚洲精品乱久久久久久| 国产成人91sexporn| 亚洲 欧美一区二区三区| 日韩中文字幕欧美一区二区 | 欧美日韩亚洲高清精品| 国产精品.久久久| 国产精品无大码| 制服人妻中文乱码| 亚洲精品中文字幕在线视频| 成年人免费黄色播放视频| 国产精品成人在线| 一区二区三区四区激情视频| 久久人人97超碰香蕉20202| 久久久久精品人妻al黑| 日韩不卡一区二区三区视频在线| 午夜老司机福利剧场| www.自偷自拍.com| 欧美xxⅹ黑人| 黑人猛操日本美女一级片| 欧美精品亚洲一区二区| 美女午夜性视频免费| 亚洲成人av在线免费| 日韩中字成人| 王馨瑶露胸无遮挡在线观看| 午夜福利视频在线观看免费| 亚洲av欧美aⅴ国产| xxx大片免费视频| 国产精品久久久久成人av| 老女人水多毛片| 久久国产精品大桥未久av| 国产麻豆69| 亚洲国产欧美网| 国产深夜福利视频在线观看| 亚洲欧美精品自产自拍| 精品一区在线观看国产| 18禁动态无遮挡网站| 国产欧美亚洲国产| 七月丁香在线播放| 亚洲欧美精品综合一区二区三区 | 99热网站在线观看| av女优亚洲男人天堂| 国产精品香港三级国产av潘金莲 | 老司机影院毛片| 女人精品久久久久毛片| 99久久精品国产国产毛片| 免费少妇av软件| 成人毛片60女人毛片免费| 国产男女内射视频| 麻豆精品久久久久久蜜桃| 深夜精品福利| 亚洲av.av天堂| 欧美国产精品va在线观看不卡| 国产精品香港三级国产av潘金莲 | 伦理电影大哥的女人| 极品人妻少妇av视频| av电影中文网址| 国产成人精品一,二区| 成人毛片60女人毛片免费| 国产精品 欧美亚洲| 最近的中文字幕免费完整| 亚洲,欧美,日韩| 最黄视频免费看| 免费观看无遮挡的男女| 亚洲av在线观看美女高潮| 丝袜美腿诱惑在线| 自拍欧美九色日韩亚洲蝌蚪91| 国产 精品1| 久久亚洲国产成人精品v| 久久精品亚洲av国产电影网| 午夜福利视频精品| 亚洲精品久久成人aⅴ小说| 亚洲av国产av综合av卡| 青青草视频在线视频观看| 黄色 视频免费看| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 久久人人97超碰香蕉20202| kizo精华| 性高湖久久久久久久久免费观看| 人成视频在线观看免费观看| 亚洲国产av影院在线观看| 久久人人爽av亚洲精品天堂| 色婷婷久久久亚洲欧美| 黑人巨大精品欧美一区二区蜜桃| 国产精品不卡视频一区二区| 成人毛片a级毛片在线播放| 日本wwww免费看| 在线观看美女被高潮喷水网站| 欧美日韩av久久| 涩涩av久久男人的天堂| 黄色一级大片看看| 大香蕉久久成人网| 中文字幕色久视频| 亚洲人成77777在线视频| 高清视频免费观看一区二区| 欧美黄色片欧美黄色片| 国产成人一区二区在线| 亚洲第一区二区三区不卡| 两性夫妻黄色片| 精品国产乱码久久久久久小说| 亚洲精品国产色婷婷电影| 亚洲一区中文字幕在线| 极品少妇高潮喷水抽搐| 性高湖久久久久久久久免费观看| 国产一区二区三区av在线| 男女高潮啪啪啪动态图| 亚洲欧美中文字幕日韩二区| 高清在线视频一区二区三区| 91精品三级在线观看| 日产精品乱码卡一卡2卡三| 精品久久蜜臀av无| 99九九在线精品视频| 晚上一个人看的免费电影| 欧美日韩综合久久久久久| 男女边吃奶边做爰视频| 婷婷成人精品国产| 寂寞人妻少妇视频99o| 成人漫画全彩无遮挡| 日韩制服丝袜自拍偷拍| 18在线观看网站| 啦啦啦视频在线资源免费观看| 一个人免费看片子| 欧美日韩视频精品一区| 人妻少妇偷人精品九色| 国产成人精品在线电影| 91久久精品国产一区二区三区| 国产欧美日韩一区二区三区在线| 电影成人av| 亚洲精品久久午夜乱码| 精品人妻在线不人妻| 制服诱惑二区| 亚洲一区中文字幕在线| 国产成人精品在线电影| 一级毛片电影观看| 精品人妻一区二区三区麻豆| 亚洲美女黄色视频免费看| 男男h啪啪无遮挡| 国产精品麻豆人妻色哟哟久久| 电影成人av| 青青草视频在线视频观看| 国产野战对白在线观看| 日本午夜av视频| 91精品伊人久久大香线蕉| 极品少妇高潮喷水抽搐| a级毛片在线看网站| 一区二区三区四区激情视频| 午夜免费男女啪啪视频观看| xxx大片免费视频| 亚洲精品中文字幕在线视频| 欧美xxⅹ黑人| 大片免费播放器 马上看| 久久精品熟女亚洲av麻豆精品| 黄片无遮挡物在线观看| 美女国产高潮福利片在线看| 大陆偷拍与自拍| av一本久久久久| 边亲边吃奶的免费视频| 久久 成人 亚洲| 久久人人97超碰香蕉20202| 十八禁网站网址无遮挡| 国产激情久久老熟女| av不卡在线播放| 精品国产一区二区三区四区第35| 99热国产这里只有精品6| 成年人午夜在线观看视频| 精品少妇黑人巨大在线播放| 中文精品一卡2卡3卡4更新| 亚洲成人一二三区av| 久久人妻熟女aⅴ| 国产成人av激情在线播放| 亚洲国产毛片av蜜桃av| 日韩精品有码人妻一区| 亚洲精品在线美女| 美女大奶头黄色视频| 少妇猛男粗大的猛烈进出视频| 成年女人在线观看亚洲视频| 成人国语在线视频| 国产欧美日韩一区二区三区在线| 国产精品偷伦视频观看了| 一本久久精品| 亚洲中文av在线| 欧美国产精品一级二级三级| 波多野结衣av一区二区av| 日本av免费视频播放| 国产老妇伦熟女老妇高清| 五月开心婷婷网| 一级毛片电影观看| 美女午夜性视频免费| 久久久久精品久久久久真实原创| 超碰97精品在线观看| 亚洲欧美一区二区三区黑人 | 国产成人午夜福利电影在线观看| xxx大片免费视频| 午夜精品国产一区二区电影| 日韩人妻精品一区2区三区| 少妇猛男粗大的猛烈进出视频| 国产午夜精品一二区理论片| 国产精品人妻久久久影院| 如何舔出高潮| 久久久久久久久久久久大奶| 亚洲四区av| 伦理电影大哥的女人| 中文字幕最新亚洲高清| 99久久人妻综合| 日韩av在线免费看完整版不卡| 日韩成人av中文字幕在线观看| 91成人精品电影| 91aial.com中文字幕在线观看| 亚洲精品久久午夜乱码| 免费在线观看视频国产中文字幕亚洲 | 午夜免费观看性视频| 日韩一本色道免费dvd| 亚洲一级一片aⅴ在线观看| 国产老妇伦熟女老妇高清| kizo精华| 日本免费在线观看一区| 视频在线观看一区二区三区| 国产97色在线日韩免费| 久久午夜福利片| 日本91视频免费播放| 边亲边吃奶的免费视频| 久久久久人妻精品一区果冻| 丝袜在线中文字幕| 夫妻性生交免费视频一级片| 日韩制服骚丝袜av| 国产精品熟女久久久久浪| 欧美变态另类bdsm刘玥| 国产精品一国产av| 啦啦啦视频在线资源免费观看| 男女午夜视频在线观看| 在线看a的网站| 高清av免费在线| 久久久精品区二区三区| 国产免费现黄频在线看| 亚洲情色 制服丝袜| 一区福利在线观看| 最近最新中文字幕免费大全7| 黄色 视频免费看| 午夜日本视频在线| 国产成人aa在线观看| av.在线天堂| 黄色毛片三级朝国网站| www.av在线官网国产| 午夜久久久在线观看| 精品亚洲成国产av| 99热网站在线观看| 高清在线视频一区二区三区| 精品久久久久久电影网| 国产在线一区二区三区精| av免费观看日本| 麻豆精品久久久久久蜜桃| 欧美精品一区二区大全| 五月天丁香电影| 国产精品一区二区在线观看99| 我的亚洲天堂| 亚洲久久久国产精品| 国产片内射在线| 国产淫语在线视频| 一边亲一边摸免费视频| 亚洲第一av免费看| 亚洲av福利一区| 久久精品国产鲁丝片午夜精品| 亚洲天堂av无毛| 成年人午夜在线观看视频| 国产白丝娇喘喷水9色精品| 国产成人av激情在线播放| 午夜免费观看性视频| 国产精品人妻久久久影院| 高清在线视频一区二区三区| 七月丁香在线播放| 亚洲av综合色区一区| 国产黄色免费在线视频| 久久久久国产精品人妻一区二区| a级毛片在线看网站| 在线 av 中文字幕| 欧美激情高清一区二区三区 | 国产精品久久久久久久久免| 国产熟女午夜一区二区三区| 亚洲av综合色区一区| 久久99精品国语久久久| 丰满饥渴人妻一区二区三| 亚洲一区二区三区欧美精品| 亚洲欧美一区二区三区久久| 最近中文字幕2019免费版| 热re99久久国产66热| 日韩制服丝袜自拍偷拍| av在线播放精品| 老汉色∧v一级毛片| 99国产精品免费福利视频| 国产成人91sexporn| 亚洲国产精品成人久久小说| 国产色婷婷99| 美女大奶头黄色视频| 综合色丁香网| 久久久久久久国产电影| 欧美激情极品国产一区二区三区| 欧美精品国产亚洲| 国产成人欧美| 男女免费视频国产| 亚洲精品久久午夜乱码| 巨乳人妻的诱惑在线观看| 观看av在线不卡| 男人爽女人下面视频在线观看| 亚洲,欧美精品.| 精品一区二区免费观看| 久久久久精品性色| 捣出白浆h1v1| 女人精品久久久久毛片| 国产日韩欧美视频二区| 一级毛片黄色毛片免费观看视频| 午夜久久久在线观看| 国产成人精品一,二区| 少妇的丰满在线观看| 国产野战对白在线观看| 亚洲国产欧美在线一区| 久久狼人影院| 黑人猛操日本美女一级片| 精品人妻在线不人妻| 欧美国产精品一级二级三级| 国产午夜精品一二区理论片| 久久综合国产亚洲精品| 精品人妻熟女毛片av久久网站| 天美传媒精品一区二区| 亚洲中文av在线| 91精品伊人久久大香线蕉| 日韩熟女老妇一区二区性免费视频| 啦啦啦视频在线资源免费观看| 岛国毛片在线播放| av卡一久久| 国产精品一区二区在线观看99| 亚洲成人手机| 人体艺术视频欧美日本| 国产精品嫩草影院av在线观看| 在线看a的网站| 国产精品熟女久久久久浪| 最近最新中文字幕大全免费视频 | 国产97色在线日韩免费| 亚洲国产毛片av蜜桃av| 在线观看三级黄色| 少妇被粗大的猛进出69影院| 成人18禁高潮啪啪吃奶动态图| 国产不卡av网站在线观看| 91久久精品国产一区二区三区| 欧美变态另类bdsm刘玥| 久久人人爽av亚洲精品天堂| 国产成人精品福利久久| 91精品伊人久久大香线蕉| 免费高清在线观看日韩| 高清欧美精品videossex| 免费女性裸体啪啪无遮挡网站| av不卡在线播放| 国产女主播在线喷水免费视频网站| 人妻少妇偷人精品九色| 国产黄频视频在线观看| 中文字幕制服av| 亚洲欧美成人综合另类久久久| 乱人伦中国视频| 日韩制服骚丝袜av| 欧美在线黄色| 搡老乐熟女国产| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 精品亚洲乱码少妇综合久久| 性少妇av在线| 国产成人精品一,二区| 亚洲国产欧美在线一区| av在线app专区| 亚洲av.av天堂| 巨乳人妻的诱惑在线观看| www.精华液| 熟女av电影| 欧美日韩视频高清一区二区三区二| 成年人免费黄色播放视频| 亚洲综合色网址| 成年女人在线观看亚洲视频| 两个人看的免费小视频| 在线天堂中文资源库| 中文字幕av电影在线播放| 人人妻人人爽人人添夜夜欢视频| 激情五月婷婷亚洲| 黄色视频在线播放观看不卡| 国产精品一二三区在线看| 欧美激情高清一区二区三区 | 制服人妻中文乱码| 黄色 视频免费看| 亚洲欧美一区二区三区黑人 | 免费看不卡的av| 国产成人免费观看mmmm| 黄色 视频免费看| 美女大奶头黄色视频| 中文字幕精品免费在线观看视频| 亚洲国产看品久久| 成人国产av品久久久| 大片电影免费在线观看免费| 国产野战对白在线观看| 亚洲欧美成人综合另类久久久| 蜜桃国产av成人99| 免费高清在线观看视频在线观看| 人妻少妇偷人精品九色| 亚洲精品中文字幕在线视频| 婷婷色综合大香蕉| 91成人精品电影| 免费久久久久久久精品成人欧美视频| 热99久久久久精品小说推荐| 黑人欧美特级aaaaaa片| 伦理电影大哥的女人| 黄片播放在线免费| 欧美变态另类bdsm刘玥| 天堂8中文在线网| 日本黄色日本黄色录像| 亚洲av电影在线观看一区二区三区| 久久婷婷青草| 在线观看美女被高潮喷水网站| 国产免费现黄频在线看| a级毛片黄视频| 最新的欧美精品一区二区| 国产精品蜜桃在线观看| 七月丁香在线播放| 赤兔流量卡办理| 黄色怎么调成土黄色| 久久久久久久久久人人人人人人| 老司机亚洲免费影院| 日韩av不卡免费在线播放| 激情视频va一区二区三区| 免费播放大片免费观看视频在线观看| 丝袜喷水一区| 青春草亚洲视频在线观看| 捣出白浆h1v1| 欧美 亚洲 国产 日韩一| 这个男人来自地球电影免费观看 | 亚洲,欧美,日韩| 王馨瑶露胸无遮挡在线观看| 久久久久久久久久久久大奶| 欧美 日韩 精品 国产| 最近最新中文字幕免费大全7| 久久精品aⅴ一区二区三区四区 | a级片在线免费高清观看视频| 校园人妻丝袜中文字幕| 狠狠婷婷综合久久久久久88av| 不卡av一区二区三区| 欧美日韩成人在线一区二区| 国产深夜福利视频在线观看| 色哟哟·www| 欧美人与善性xxx| 午夜老司机福利剧场| 日韩一本色道免费dvd| 晚上一个人看的免费电影| 一本色道久久久久久精品综合| 久久久精品免费免费高清| 免费观看性生交大片5| 日本免费在线观看一区| 亚洲五月色婷婷综合| 又大又黄又爽视频免费| 国产xxxxx性猛交| 青青草视频在线视频观看| 亚洲成色77777| 中文乱码字字幕精品一区二区三区| 两性夫妻黄色片| 亚洲伊人色综图| 夫妻午夜视频| 久久久国产精品麻豆| 国产老妇伦熟女老妇高清| 国产精品无大码| 日本91视频免费播放| 中国国产av一级| 波野结衣二区三区在线| 国产 精品1| 26uuu在线亚洲综合色| 男人添女人高潮全过程视频| 亚洲精品久久午夜乱码| 免费观看av网站的网址| 久久久久久免费高清国产稀缺| 一级毛片 在线播放| 考比视频在线观看| 国产av精品麻豆| 欧美黄色片欧美黄色片| 18禁裸乳无遮挡动漫免费视频| 亚洲国产精品国产精品| 性少妇av在线| 精品亚洲成a人片在线观看| 99热国产这里只有精品6| 免费看av在线观看网站| 超色免费av| 97人妻天天添夜夜摸| 亚洲av国产av综合av卡| 最近的中文字幕免费完整| 日产精品乱码卡一卡2卡三| 人妻系列 视频| 国产精品偷伦视频观看了| 婷婷成人精品国产| 欧美日韩一级在线毛片| 韩国精品一区二区三区| 亚洲国产av新网站| videos熟女内射| 777米奇影视久久| 考比视频在线观看| 成人漫画全彩无遮挡| 国产一区二区三区综合在线观看| 啦啦啦在线免费观看视频4| 亚洲成色77777| 又黄又粗又硬又大视频| 国产成人精品久久久久久| 久久久久精品性色| 2022亚洲国产成人精品| 热re99久久精品国产66热6| 欧美日韩视频精品一区| 精品99又大又爽又粗少妇毛片| 人人妻人人澡人人看| 午夜福利视频在线观看免费| 最近最新中文字幕免费大全7| 美女中出高潮动态图| 国产1区2区3区精品| 男女啪啪激烈高潮av片| 极品少妇高潮喷水抽搐| 国产成人av激情在线播放| 男女免费视频国产| 18禁观看日本| 国产精品欧美亚洲77777| 亚洲经典国产精华液单| 亚洲五月色婷婷综合| 十分钟在线观看高清视频www| 欧美+日韩+精品| 中文乱码字字幕精品一区二区三区| 欧美bdsm另类| 国产又色又爽无遮挡免| 一区二区三区乱码不卡18| 街头女战士在线观看网站| 国产免费又黄又爽又色| 91精品国产国语对白视频| 日韩人妻精品一区2区三区| 韩国精品一区二区三区| 久久热在线av| 国产福利在线免费观看视频| 精品99又大又爽又粗少妇毛片| 国产片内射在线| www.av在线官网国产| 一本大道久久a久久精品| 久久久亚洲精品成人影院| 亚洲国产精品一区三区| 飞空精品影院首页| 看非洲黑人一级黄片| 精品亚洲成a人片在线观看| 桃花免费在线播放| 日本欧美视频一区| 黑人猛操日本美女一级片| 十八禁网站网址无遮挡| 日韩欧美精品免费久久| 国产女主播在线喷水免费视频网站| 麻豆乱淫一区二区| 又大又黄又爽视频免费| 亚洲国产毛片av蜜桃av| 亚洲精品国产色婷婷电影| 不卡视频在线观看欧美| 一级毛片我不卡| a级毛片在线看网站| 久热久热在线精品观看| 日韩一卡2卡3卡4卡2021年| 男女无遮挡免费网站观看| 美女xxoo啪啪120秒动态图| 黄色一级大片看看| 日本av免费视频播放| 亚洲伊人色综图| 午夜久久久在线观看| 久久久久国产网址| 免费日韩欧美在线观看| 视频在线观看一区二区三区| 欧美日韩亚洲高清精品|