• <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

    亚洲av电影不卡..在线观看| av福利片在线观看| 日韩欧美国产一区二区入口| 天天躁狠狠躁夜夜躁狠狠躁| 哪里可以看免费的av片| 成人午夜高清在线视频| 亚洲第一欧美日韩一区二区三区| 91在线观看av| 久久久久亚洲av毛片大全| 男女做爰动态图高潮gif福利片| 黄色视频不卡| 欧美色视频一区免费| 国内精品久久久久久久电影| 日日摸夜夜添夜夜添小说| 一级黄色大片毛片| 亚洲国产看品久久| 一个人观看的视频www高清免费观看 | 亚洲国产精品合色在线| 男女午夜视频在线观看| 成在线人永久免费视频| 亚洲人成伊人成综合网2020| 精品久久久久久久久久免费视频| 99精品在免费线老司机午夜| 国产片内射在线| 三级毛片av免费| 无人区码免费观看不卡| 久久国产精品人妻蜜桃| 老司机在亚洲福利影院| 人人妻人人澡欧美一区二区| 亚洲18禁久久av| 亚洲五月天丁香| 久久这里只有精品19| 高清在线国产一区| 精品少妇一区二区三区视频日本电影| 午夜a级毛片| 久久中文字幕一级| 久9热在线精品视频| 午夜精品久久久久久毛片777| 亚洲avbb在线观看| 国产精品国产高清国产av| 亚洲成人久久性| 免费av毛片视频| 最近最新中文字幕大全免费视频| 免费在线观看亚洲国产| 欧美高清成人免费视频www| 日本黄色视频三级网站网址| 国产又色又爽无遮挡免费看| 亚洲一码二码三码区别大吗| 美女扒开内裤让男人捅视频| 亚洲欧美日韩高清在线视频| 日本熟妇午夜| 久久人妻av系列| 久久久久久久久中文| 欧美精品啪啪一区二区三区| 国产精品一区二区三区四区免费观看 | 美女扒开内裤让男人捅视频| 一级作爱视频免费观看| 欧美成人一区二区免费高清观看 | 色精品久久人妻99蜜桃| 精品欧美国产一区二区三| 精品久久蜜臀av无| 天堂影院成人在线观看| 久久久国产成人精品二区| 国产区一区二久久| 校园春色视频在线观看| 可以在线观看毛片的网站| 日韩国内少妇激情av| 在线观看舔阴道视频| 好看av亚洲va欧美ⅴa在| 成人欧美大片| 国产aⅴ精品一区二区三区波| 日韩欧美在线二视频| 男人的好看免费观看在线视频 | 亚洲男人的天堂狠狠| 亚洲欧美精品综合久久99| 一卡2卡三卡四卡精品乱码亚洲| 久久午夜综合久久蜜桃| 制服人妻中文乱码| 国产探花在线观看一区二区| 91九色精品人成在线观看| 国内精品一区二区在线观看| 婷婷精品国产亚洲av| 久久精品91蜜桃| 无遮挡黄片免费观看| 一区福利在线观看| av在线播放免费不卡| 色噜噜av男人的天堂激情| 久久精品影院6| 日本免费a在线| 亚洲国产精品成人综合色| 最近最新中文字幕大全免费视频| 国产欧美日韩一区二区三| www日本黄色视频网| 人人妻人人澡欧美一区二区| 亚洲成人中文字幕在线播放| 精品国产超薄肉色丝袜足j| 成人国产综合亚洲| 婷婷亚洲欧美| 欧美精品亚洲一区二区| 色老头精品视频在线观看| 久久久久久亚洲精品国产蜜桃av| 中国美女看黄片| aaaaa片日本免费| 国产激情欧美一区二区| 我的老师免费观看完整版| 欧美黄色片欧美黄色片| 欧美色视频一区免费| 可以在线观看毛片的网站| 日韩免费av在线播放| 欧美一区二区国产精品久久精品 | 亚洲av成人一区二区三| 男人舔奶头视频| 极品教师在线免费播放| 女人被狂操c到高潮| 成人欧美大片| www日本黄色视频网| 69av精品久久久久久| 91麻豆av在线| 动漫黄色视频在线观看| 亚洲熟妇熟女久久| 精品无人区乱码1区二区| 国产亚洲精品综合一区在线观看 | 香蕉国产在线看| 级片在线观看| 成人三级做爰电影| 久久精品夜夜夜夜夜久久蜜豆 | 亚洲片人在线观看| 一个人观看的视频www高清免费观看 | 日韩有码中文字幕| 欧美激情久久久久久爽电影| 美女黄网站色视频| 三级毛片av免费| 欧美 亚洲 国产 日韩一| 村上凉子中文字幕在线| a级毛片在线看网站| 亚洲第一电影网av| 久久精品人妻少妇| 亚洲国产欧美一区二区综合| 日韩大尺度精品在线看网址| 亚洲精品一区av在线观看| 国产精品爽爽va在线观看网站| 99热6这里只有精品| 国产熟女xx| 国产真人三级小视频在线观看| 99国产精品一区二区蜜桃av| 国产精品一区二区免费欧美| 久久久久亚洲av毛片大全| 免费观看人在逋| 国产精品久久电影中文字幕| 不卡一级毛片| 亚洲七黄色美女视频| 亚洲男人的天堂狠狠| 欧美成人午夜精品| 久久精品综合一区二区三区| 成人三级做爰电影| 男人舔女人的私密视频| 最近最新中文字幕大全免费视频| 午夜精品一区二区三区免费看| 成人av一区二区三区在线看| 国产精品精品国产色婷婷| 亚洲成人久久爱视频| 香蕉久久夜色| 午夜a级毛片| 1024手机看黄色片| 欧美激情久久久久久爽电影| 九色成人免费人妻av| 国产精品美女特级片免费视频播放器 | 亚洲国产日韩欧美精品在线观看 | 日本黄大片高清| 国产亚洲欧美在线一区二区| 两人在一起打扑克的视频| 欧美绝顶高潮抽搐喷水| 色综合亚洲欧美另类图片| tocl精华| 国产精品一区二区三区四区免费观看 | 国产精品久久久久久久电影 | 男女做爰动态图高潮gif福利片| 一本久久中文字幕| 亚洲自偷自拍图片 自拍| 一区福利在线观看| 热99re8久久精品国产| 国产伦一二天堂av在线观看| 美女午夜性视频免费| 麻豆av在线久日| 日日夜夜操网爽| 亚洲国产精品sss在线观看| 欧美黑人巨大hd| 99热这里只有是精品50| 午夜福利18| 国模一区二区三区四区视频 | 国产一区二区三区视频了| 国产精品爽爽va在线观看网站| 亚洲第一欧美日韩一区二区三区| 黄色 视频免费看| 日韩欧美免费精品| 免费观看精品视频网站| 国产精品综合久久久久久久免费| 岛国在线免费视频观看| 日本一本二区三区精品| 亚洲精品在线观看二区| 一本久久中文字幕| 曰老女人黄片| 91成年电影在线观看| www国产在线视频色| 国语自产精品视频在线第100页| 免费在线观看完整版高清| 中文字幕精品亚洲无线码一区| 久久久久免费精品人妻一区二区| 午夜免费观看网址| 国产成人欧美在线观看| 亚洲精品中文字幕在线视频| av在线天堂中文字幕| 特大巨黑吊av在线直播| 亚洲成人久久性| 亚洲aⅴ乱码一区二区在线播放 | 欧美性猛交╳xxx乱大交人| 中出人妻视频一区二区| 亚洲一区二区三区色噜噜| 成年人黄色毛片网站| 中文字幕最新亚洲高清| 好男人电影高清在线观看| 国产伦人伦偷精品视频| 欧美国产日韩亚洲一区| 男人舔奶头视频| 日韩大码丰满熟妇| 手机成人av网站| 女生性感内裤真人,穿戴方法视频| 老司机午夜福利在线观看视频| 婷婷精品国产亚洲av| 18禁国产床啪视频网站| 性欧美人与动物交配| 亚洲欧美精品综合一区二区三区| 亚洲av日韩精品久久久久久密| 可以免费在线观看a视频的电影网站| 小说图片视频综合网站| 男女之事视频高清在线观看| 国产探花在线观看一区二区| 国产真人三级小视频在线观看| 男人的好看免费观看在线视频 | 制服丝袜大香蕉在线| www.精华液| 好男人电影高清在线观看| 啪啪无遮挡十八禁网站| 亚洲中文日韩欧美视频| 午夜影院日韩av| 久久性视频一级片| 三级男女做爰猛烈吃奶摸视频| АⅤ资源中文在线天堂| 欧美日韩亚洲国产一区二区在线观看| 日韩欧美 国产精品| 男女之事视频高清在线观看| 99热这里只有是精品50| 国产男靠女视频免费网站| 国产精品久久视频播放| 99国产精品一区二区三区| 99热只有精品国产| 日本三级黄在线观看| 日韩免费av在线播放| 久久亚洲真实| 黄色毛片三级朝国网站| 国产精品1区2区在线观看.| 悠悠久久av| 在线观看一区二区三区| 黄色视频不卡| 正在播放国产对白刺激| 亚洲国产精品sss在线观看| 好男人在线观看高清免费视频| 欧美乱色亚洲激情| 免费高清视频大片| 精品久久久久久久久久免费视频| 欧美黑人精品巨大| 色av中文字幕| 国产成人系列免费观看| 18禁黄网站禁片午夜丰满| 黄色成人免费大全| 成人午夜高清在线视频| 好看av亚洲va欧美ⅴa在| 美女高潮喷水抽搐中文字幕| 国产精品美女特级片免费视频播放器 | 又黄又爽又免费观看的视频| av有码第一页| 法律面前人人平等表现在哪些方面| 国产精品综合久久久久久久免费| 在线观看午夜福利视频| 亚洲av成人一区二区三| 亚洲激情在线av| 午夜福利成人在线免费观看| 亚洲美女黄片视频| 日本撒尿小便嘘嘘汇集6| 成人特级黄色片久久久久久久| 制服丝袜大香蕉在线| 国产探花在线观看一区二区| 亚洲美女视频黄频| 久久草成人影院| 精品电影一区二区在线| 久久天躁狠狠躁夜夜2o2o| 香蕉久久夜色| 99国产综合亚洲精品| 50天的宝宝边吃奶边哭怎么回事| 久久精品aⅴ一区二区三区四区| x7x7x7水蜜桃| 91字幕亚洲| 久久性视频一级片| 精品一区二区三区视频在线观看免费| 亚洲精品久久国产高清桃花| 超碰成人久久| 制服丝袜大香蕉在线| 丁香六月欧美| av有码第一页| 脱女人内裤的视频| 日韩三级视频一区二区三区| 91老司机精品| 美女高潮喷水抽搐中文字幕| 午夜久久久久精精品| 亚洲真实伦在线观看| 久久性视频一级片| 亚洲中文字幕一区二区三区有码在线看 | 国产精品 国内视频| 人人妻人人看人人澡| 亚洲18禁久久av| 最近最新中文字幕大全电影3| 国产成人系列免费观看| 成年免费大片在线观看| 黄色 视频免费看| 91国产中文字幕| 一二三四在线观看免费中文在| 一个人免费在线观看的高清视频| 国产在线精品亚洲第一网站| 久久婷婷人人爽人人干人人爱| 非洲黑人性xxxx精品又粗又长| 夜夜看夜夜爽夜夜摸| а√天堂www在线а√下载| 波多野结衣巨乳人妻| 欧美黑人欧美精品刺激| 丰满人妻熟妇乱又伦精品不卡| 欧美日本亚洲视频在线播放| 国产成人av教育| 动漫黄色视频在线观看| 男人的好看免费观看在线视频 | 欧美成人一区二区免费高清观看 | 免费在线观看影片大全网站| 十八禁人妻一区二区| 欧美乱色亚洲激情| 亚洲精品中文字幕一二三四区| 亚洲午夜理论影院| 久久香蕉国产精品| 精品国内亚洲2022精品成人| 久久久久久久久久黄片| 精品久久久久久久久久久久久| 丁香欧美五月| 老司机午夜福利在线观看视频| 大型黄色视频在线免费观看| 亚洲狠狠婷婷综合久久图片| 夜夜躁狠狠躁天天躁| 久久精品国产清高在天天线| 国产精品 欧美亚洲| 美女午夜性视频免费| 欧美丝袜亚洲另类 | 亚洲18禁久久av| 观看免费一级毛片| 亚洲精品色激情综合| 99热这里只有精品一区 | 人人妻人人澡欧美一区二区| 欧美色欧美亚洲另类二区| 在线观看www视频免费| 在线观看美女被高潮喷水网站 | av天堂在线播放| 18禁黄网站禁片午夜丰满| 色哟哟哟哟哟哟| 两性夫妻黄色片| 亚洲电影在线观看av| 麻豆av在线久日| 国产精品1区2区在线观看.| 亚洲欧美一区二区三区黑人| 一本大道久久a久久精品| 色综合婷婷激情| 天堂动漫精品| 精品久久久久久久久久久久久| 五月伊人婷婷丁香| 蜜桃久久精品国产亚洲av| 99国产极品粉嫩在线观看| 男女视频在线观看网站免费 | 免费人成视频x8x8入口观看| 欧美性长视频在线观看| 两个人视频免费观看高清| 欧美乱色亚洲激情| 国产亚洲精品久久久久5区| 国内精品久久久久精免费| 国产91精品成人一区二区三区| 日韩大尺度精品在线看网址| 欧美又色又爽又黄视频| 日本一区二区免费在线视频| 亚洲av成人不卡在线观看播放网| 黄色片一级片一级黄色片| 老司机靠b影院| 国产精品美女特级片免费视频播放器 | 国内少妇人妻偷人精品xxx网站 | 一级毛片女人18水好多| 久久婷婷人人爽人人干人人爱| 蜜桃久久精品国产亚洲av| 男女下面进入的视频免费午夜| 精品日产1卡2卡| 国产又色又爽无遮挡免费看| 午夜成年电影在线免费观看| 国产精品av视频在线免费观看| 欧美成人免费av一区二区三区| 搡老岳熟女国产| 狂野欧美激情性xxxx| 在线视频色国产色| 美女扒开内裤让男人捅视频| 国产真人三级小视频在线观看| 中出人妻视频一区二区| 男人舔女人下体高潮全视频| 成人三级做爰电影| 国产91精品成人一区二区三区| 欧美大码av| 色综合亚洲欧美另类图片| 欧洲精品卡2卡3卡4卡5卡区| 1024香蕉在线观看| 男女那种视频在线观看| 日韩欧美国产一区二区入口| 大型av网站在线播放| 色综合欧美亚洲国产小说| 日韩精品中文字幕看吧| 99热只有精品国产| 亚洲第一欧美日韩一区二区三区| 日本黄色视频三级网站网址| 日本a在线网址| 精品久久久久久久人妻蜜臀av| 黄色女人牲交| 欧美3d第一页| 一区二区三区高清视频在线| av超薄肉色丝袜交足视频| 看片在线看免费视频| 色综合亚洲欧美另类图片| 国产成人aa在线观看| 免费无遮挡裸体视频| 黄色视频不卡| 日本在线视频免费播放| 大型黄色视频在线免费观看| 欧美日韩福利视频一区二区| av免费在线观看网站| 听说在线观看完整版免费高清| 国产成人精品久久二区二区免费| 国内毛片毛片毛片毛片毛片| 日本熟妇午夜| 一本久久中文字幕| 日韩精品青青久久久久久| 欧美色欧美亚洲另类二区| 精品久久久久久,| 法律面前人人平等表现在哪些方面| 久久亚洲精品不卡| 亚洲人成网站高清观看| 九色成人免费人妻av| 99国产精品99久久久久| 欧美一区二区国产精品久久精品 | 亚洲国产看品久久| 欧美大码av| 天堂动漫精品| 欧美一级a爱片免费观看看 | 日韩成人在线观看一区二区三区| 老司机靠b影院| 男人舔女人的私密视频| 久久香蕉激情| 91九色精品人成在线观看| 一级a爱片免费观看的视频| 久久国产乱子伦精品免费另类| 成人午夜高清在线视频| 成人午夜高清在线视频| 日日摸夜夜添夜夜添小说| 国产精品久久久人人做人人爽| 精品无人区乱码1区二区| 我的老师免费观看完整版| 亚洲七黄色美女视频| 好看av亚洲va欧美ⅴa在| 午夜免费观看网址| 大型av网站在线播放| 国产伦在线观看视频一区| 免费av毛片视频| 婷婷丁香在线五月| 999精品在线视频| 99riav亚洲国产免费| 麻豆成人av在线观看| www.www免费av| 精品一区二区三区av网在线观看| 午夜福利18| 狂野欧美激情性xxxx| 少妇粗大呻吟视频| 色精品久久人妻99蜜桃| 亚洲国产精品合色在线| 国产成人精品无人区| 欧美日韩瑟瑟在线播放| or卡值多少钱| 久久久精品大字幕| 久久精品夜夜夜夜夜久久蜜豆 | 国产1区2区3区精品| 十八禁网站免费在线| 制服人妻中文乱码| 草草在线视频免费看| 久久久久国产一级毛片高清牌| 亚洲专区中文字幕在线| 欧美色视频一区免费| 老司机靠b影院| 好男人在线观看高清免费视频| 欧美三级亚洲精品| 这个男人来自地球电影免费观看| 国产精品98久久久久久宅男小说| 亚洲黑人精品在线| 国产黄片美女视频| 国产精品久久视频播放| 精品电影一区二区在线| 法律面前人人平等表现在哪些方面| 亚洲熟妇熟女久久| 精品国产乱子伦一区二区三区| 99国产综合亚洲精品| 国产精品亚洲一级av第二区| 别揉我奶头~嗯~啊~动态视频| 久久国产乱子伦精品免费另类| 成人18禁高潮啪啪吃奶动态图| 成在线人永久免费视频| 午夜两性在线视频| 欧美黄色淫秽网站| 18禁国产床啪视频网站| 男男h啪啪无遮挡| 人成视频在线观看免费观看| 日本免费一区二区三区高清不卡| 亚洲成a人片在线一区二区| 精品久久久久久,| 99久久综合精品五月天人人| av在线播放免费不卡| 国产高清videossex| 精品乱码久久久久久99久播| 亚洲人成网站高清观看| 亚洲精品国产精品久久久不卡| 在线国产一区二区在线| 免费看美女性在线毛片视频| 亚洲va日本ⅴa欧美va伊人久久| 成人午夜高清在线视频| 婷婷亚洲欧美| 亚洲黑人精品在线| 亚洲av五月六月丁香网| 一个人观看的视频www高清免费观看 | 亚洲精品久久成人aⅴ小说| 高清在线国产一区| 亚洲av美国av| 50天的宝宝边吃奶边哭怎么回事| 亚洲熟妇熟女久久| 亚洲国产欧美网| 国产伦人伦偷精品视频| 亚洲自偷自拍图片 自拍| 久久精品国产亚洲av高清一级| 在线观看免费视频日本深夜| e午夜精品久久久久久久| 久久天堂一区二区三区四区| 91在线观看av| 青草久久国产| 日韩 欧美 亚洲 中文字幕| 脱女人内裤的视频| 夜夜爽天天搞| 精品人妻1区二区| 欧美久久黑人一区二区| 亚洲欧美一区二区三区黑人| 少妇粗大呻吟视频| 欧美国产日韩亚洲一区| 国产精品自产拍在线观看55亚洲| 亚洲电影在线观看av| 精品久久蜜臀av无| 久久香蕉国产精品| 麻豆成人av在线观看| 欧美黄色淫秽网站| 99在线人妻在线中文字幕| 午夜激情av网站| aaaaa片日本免费| 国产精品一及| 国产精品亚洲美女久久久| 青草久久国产| 在线观看一区二区三区| 欧美高清成人免费视频www| 日本一二三区视频观看| 巨乳人妻的诱惑在线观看| 天堂av国产一区二区熟女人妻 | 欧美 亚洲 国产 日韩一| 精品欧美国产一区二区三| 国产精品香港三级国产av潘金莲| 色综合婷婷激情| 青草久久国产| avwww免费| 中文字幕高清在线视频| 又大又爽又粗| av中文乱码字幕在线| 高潮久久久久久久久久久不卡| 啪啪无遮挡十八禁网站| 男女午夜视频在线观看| 免费在线观看日本一区| 久久精品91无色码中文字幕| 日韩欧美精品v在线| 看免费av毛片| 免费看十八禁软件| av天堂在线播放| 中文字幕久久专区| 麻豆av在线久日| 国产精品久久久av美女十八| 黄色毛片三级朝国网站| 18禁裸乳无遮挡免费网站照片| 精华霜和精华液先用哪个| 啦啦啦观看免费观看视频高清| 露出奶头的视频| 激情在线观看视频在线高清| 亚洲欧美一区二区三区黑人| 久久久精品国产亚洲av高清涩受| 999精品在线视频| 欧美色欧美亚洲另类二区| 美女午夜性视频免费| 黑人欧美特级aaaaaa片| 19禁男女啪啪无遮挡网站|