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

    Short- and Long-Term Response of Phytoplankton to ENSO in Prydz Bay,Antarctica:Evidences from Field Measurements,Remote Sensing Data and Stratigraphic Biomarker Records

    2014-04-17 02:42:46ZHAOJunHansUlrichPeterZHANGHaishengHANZhengbingHUChuanyuYUPeisongLUBingandThomasBianchi
    Journal of Ocean University of China 2014年3期

    ZHAO Jun,Hans-Ulrich Peter,ZHANG Haisheng,,HAN Zhengbing,HU Chuanyu,YU Peisong,LU Bing,and Thomas S.Bianchi

    1)Second Institute of Oceanography,State Oceanic Administration(SOA),Hangzhou 310012,P. R.China

    2) Key Laboratory of Marine Ecosystem and Biogeochemistry, SOA,Hangzhou 310012,P. R.China

    3) Institute of Ecology,Friedrich Schiller University, Jena D-07743, Germany

    4) Department of Geological Sciences,University of Florida,Gainesville FL 32611-2120, USA

    ? Ocean University of China,Science Press and Springer-Verlag Berlin Heidelberg 2014

    1 Introduction

    ENSO(El Ni?o-Southern Oscillation)with warm(El Ni?o)and cold(La Ni?a)phases is rooted in the tropical Pacific,involving large-scale ocean-atmosphere interactions,and is the strongest natural climate variation at interannual scales(Changet al.,2007; McPhadenet al.,2006; Philander and Fedorov,2003; Trenberth and Hoar,1997).The impact of El Ni?o and La Ni?a on marine phytoplankton is evolving as a spotlight with great interests in understanding the causes and consequences of ENSO(Chavezet al.,2011; Doneyet al.,2012).Past work has shown that marine phytoplankton growth rate and community structure change are correlated significantly with global warming and El Ni?o/La Ni?a events(Behrenfeldet al.,2006; Rittenouret al.,2000).More recent studies similarly have corroborated that variation of marine phytoplankton primary production is strongly correlated with(extended)Multivariate ENSO Index(MEI or MEI.ext)(Wolter and Timlin,1983,2011); an annual response to El Ni?o/La Ni?a was especially visible(Hernández-de-la-Torreet al.,2003; Mackeyet al.,2010).

    The Southern Ocean is a typical High Nutrient-Low Chlorophyll(HNLC)region with a simple marine food chain and remains one of the least studied of all the oceans.We do know that marine phytoplankton primary production and community structure have an important effect on this important global marine ecosystem and appear to be extremely sensitive to climate change(Arrigoet al.,2008; Smetacek and Nicol,2005).The relationship between ENSO events and annual variation of SST,as well as phytoplankton biomass(chlorophyll-a,Chl-a)around the Antarctic Peninsula region and the Ross Sea,has shown more subtle consequences of ENSO variability on marine ecosystem responses(Arrigo and van Dijken,2004; Loebet al.,2009; Reisset al.,2009).These studies highlight that high latitude oceans are sensitive to global impact of ENSO events.Although there is no unified understanding on the biological responses and feedback mechanisms in low and high latitude regions to ENSO,we do recognize the need for a better understanding in the polar oceans,which have been impacted by changes in both oceanic and atmospheric circulations(Sarmientoet al.,2004).While the Arctic has received considerably more attention in recent years(Stein and Macdonald,2004),there is perhaps greater need to better understand marine ecosystem variability associated with ENSO in the Southern Ocean forcing by examining internally consistent short- and long-term multidisciplinary data.

    Diagnostic lipid biomarkers in stratigraphic records have been successfully employed to reconstruct historical changes of phytoplankton relative abundance and community structure(Sicreet al.,2000; Bianchi and Canuel,2011).For example,brassicasterol,dinosterol,and alkenones(C37and C38)are indicators of diatoms,dinoflagellates and haptophytes(mainlyEmiliania huxleyi),respectively(Volkman,2006).Relative proportions of biomarkers can be used to reveal phytoplankton community variation,while their bulk concentrations are indicators of phytoplankton production since diatoms,dinoflagellates and haptophytes are dominant algae in the Antarctic phytoplankton community(Wright and van den Enden,2000; Wrightet al.,2010).

    In this study,we present 1)annual field measurements of nutrients,Chl-a,phytoplankton abundance and community in austral summer during 1990 to 2002; 2)monthly remote sensing data of Chl-a,particulate organic carbon(POC)and sea surface temperature(SST)in austral summer during 2007 to 2011; and 3)long-term stratigraphic lipid biomarker records during 1900 to 2005 in Prydz Bay,East Antarctica.Then we compare these data with ENSO variation indicated by MEI or MEI.ext.The primary goal of this work is to explore short- and long-term response of phytoplankton to ENSO variability using algal biomarkers in the sediment record.

    2 Materials and Methods

    2.1 Site Description

    Prydz Bay is the largest embayment and a typical marginal ice zone in the Indian sector of the Southern Ocean,East Antarctic Margin,between approximately 69?E and 80?E,70?S and 67?S(Fig.1).It is bounded to the northeast by the Four Ladies Bank,to the southeast by the Ingrid Christensen Coast-Larsemann Hills Oasis,to the south by the Amery Ice Shelf-Lambert Glacier ice drainage system,to the west by Mac.Robertson Land,to the northwest by the Fram Bank,and to the north by the continental shelf edge(Pu and Dong,2003).Prydz Bay is trenched by the Prydz Channel to the north and the Svenner Channel to the east,and the Four Ladies Bank and Fram Bank form partial barriers to deep-water exchange between the bay and the deep ocean.Usually,the bay is covered by sea ice in winter and spring,but several polynyas occur far earlier than is usual for such high latitudes(Smithet al.,1984).The Chinese Zhongshan Station(69?22′24′S,76?22′40′′E)is located on Prydz Bay,and has been in operation for the past 20 years.

    Fig.1 Study area located in Prydz Bay,East Antarctica.Black dot shows where sediment core sample III-12(67.49?S,73.00?E)was collected by an undisturbed multicorer during CHINARE-21,2005.

    2.2 Sampling and Data Collection

    Sediment core(III-12,67.49?S,73.00?E)was collected using an undisturbed multi-corer onboard theR/V Xuelongfrom Prydz Bay on February 2,2005 during the 21stChinese National Antarctic Research Expedition(CHINARE-21)(Fig.1).The sample was then stored frozen onboard(?20℃)and transferred to a ?20℃ freezer in the laboratory for storage,prior to biomarkers and radionuclides analyses.Annual ship-based field data of sea water temperature,dissolved oxygen(DO),Chl-a,nutrients and information of phytoplankton species in Prydz Bay from 1991 to 2002(during CHINARE-6 to 18)were cited from references.Monthly ocean color remote sensing concentrations of Chl-a,POC and SST were from satellite remote sensing database of MODIS.

    2.3 Sediment Chronolo gy

    Core III-12 was 28 cm in length.It was subsampled at 0.5 cm intervals on the top 1 cm,at 1 cm from 1 to 10 cm,and at 2 cm below 10 cm.The subsamples were then freeze-dried with a Freezone-6 freeze-dryer(Labconco,USA).The210Pb activities(half life = 22.3 years)were determined by alpha particle spectrometry methods according to Chenet al.(2001),and whole down-core excess210Pb(210Pbxs)vales were utilized for calculating sediment accumulation rates based on the calculation methods of Nittrouer and Sternberg(1981)since no surface mixed-layer was observed.Briefly,a constant initial concentration(CIC)model was employed,and the function is

    wheresis the sedimentation rate(cm year?1),λis the ra-dioactive decay constant for210Pb(0.0311 year?1),kis the slope(linear coefficient)between natural logarithms of210Pbxsradioactivity(x-axis)and depth(y-axis).

    2.4 Biomarker Extraction and GC Analysis

    Sedimentary biomarkers were extracted according to Zhanget al.(2012).Briefly,approximately 10 g homogenized freeze-dried sediment was mixed with methylene chloride/methanol(v/v = 3/1),and 20μL(16.19 mg)r-C30alkane(squalane)and 20μL(16.74 mg)n-C19alcohol as internal standards.The mixtures were then extracted for 48 h by Soxhlet extraction with water?bath temperature ca.62℃ and return velocity 4?5 times h1.The extracts were desulfurized by active copper overnight and concentrated by rotary evaporation.After that,10 mL 1 mol L?1KOH/ methanol(v/v = 5/95)solvent were added to alkaline hydrolyse the extracts and to remove acidic components.Then the extracts were centrifuged and re-extracted with hexane for four times.The supernatant was decanted into a rotary evaporator so as to be concentrated,blown to dryness with nitrogen and re-dissolved with hexane.Silica column chromatography method was employed to separate different components.Sterols and ketenes were separated by methylene chloride/methanol(v/v = 95/5),blown to dryness with nitrogen and re- dissolved by hexane to 1 mL vials.They were then transformed into trimethylsilyl derivative(TMS-ether)by BS-TFA [Bis(trimethyl)-trifluoroa-cetamide] for GC analysis.

    Fig.2(A)Structures of sedimentary biomarkers:brassicasterol,dinosterol and alkenones(C37 and C38),and(B)their GC chromatograms.

    The extract was injected into an HP5890 GC coupled with an elastic quartz capillary column(DB-5,30 m × 0.25 mm i.d.,0.25μm thickness of membrane).The temperature was 260℃ at sample entrance,and temperature-programming was:80℃ at the very beginning,followed by being ramped to 200℃ with a speed of 25℃min?1,to 250℃ with a speed of 4℃ min?1,then to 300℃with a speed of 1.8℃ min?1,and finally to 310℃ with a speed of 5℃min?1to hold isocratically for 2 min before next injection.Hydrogen was employed as carrier gas,with a flow velocity 1.2 mL min?1.Identification of biomarkers was performed according to the retention time of chromatograms,and their quantifications were based on the integration areas of chromatograms using response factor of internal standards(Volkman,2006).The biomarker structures and sample chromatograms are shown in Fig.2.

    2.5 Field Measurements of Chl-a

    Chl-awas determined by extraction fluorometric method.Briefly,water sample(100 ? 500 mL)was filtered by GF/F filter.Chl-awas extracted by 90% acetone for 24 h under ?20℃,and then determined using Turner Designs Fluorometer.

    3 Results and Discussion

    3.1 Annual Field Measured Phytoplankton and Environmental Factors(1990-2002),and Their Relationship with ENSO Events

    Annual ship-based field data of water temperature,salinity,dissolved oxygen,nutrients,Chl-a,phytoplankton standing stock and community of 7 years within the period from 1990 to 2002 in Prydz Bay,East Antarctica(Table 1)have been employed to study the responding of marine environment and ecosystem change to El Ni?o and La Ni?a events.The interannual variability of the above parameters was significant in the bay according to the data in Table 1 and references therein.During El Ni?o years(March 1991? June 1992 and March 1993? September 1993),the proportion of bacillariophyta increased(84.8% and 78.6%)relative to normal years(77.9% and 73.0%),while those of other algae(pyrrhophyta,chrysophyta,etc.)decreased;Nitzschia barkleyiandNitzschia cylindrusbloomed rapidly and dominated the phytoplankton community(Ning,1998).During La Ni?a years(December 1998 ? July 2000),the proportion of bacillariophyta decreased(61.1% and 70.5%)relative to normal years,while those of other algae increased,withNitzschia curtaas the dominant phytoplankton species(Caiet al.,2003; Liuet al.,2001).Although bacillariophyta dominated the algal community,its proportion varied significantly during El Ni?o or La Ni?a events.Phytoplankton species decreased during both El Ni?o(78 and 95)and La Ni?a(36 and 65)years,compared with normal years(145 and 122)(Caiet al.,2003; Liuet al.,2001; Ning,1998; Zhuet al.,1994).These results suggest that changes in marine phytoplankton composition are linked to both warm and cold climate anomalies(El Ni?o and La Ni?a events)in the Southern Ocean.

    Table 1 Field measurements of environmental factors,phytoplankton community and relative abundance in Prydz Bay,1990-2002,and their relationship with ENSO events

    In January and February,2002,the normal year between La Ni?a(December 1998 ? July 2000)and El Ni?o(May 2002 ? March 2003)events,the proportion of bacillariophyta was lower(68.6%)than in other normal years(January ? March 1990 and December 1990 ? February 1991),andFragilariopsis kerguelensiswas the dominant species(Liuet al.,2004).Field measurements of phytoplankton species(86)during this period increased sig-nificantly compared to those in previous La Ni?a events,but were still less than those in other normal years(Caiet al.,2003; Liuet al.,2004; Liuet al.,2001; Ning,1998; Zhuet al.,1994).This indicates that the recovery time from these climate fluctuations for these phytoplankton may be quite long in the Antarctic,based on the connections with El Ni?o during May 2002 and March 2003.These results have significant implications for future impacts of climate on Southern Ocean phytoplankton communities.However,further research is necessary to better understand and document the linkages between physical forcings and biotic response times.

    Physicochemical factors,such as nutrients,dissolved oxygen and salinity,were also associated with the changes in phytoplankton composition and varied among different years(normal,El Ni?o and La Ni?a).Although we only had a very limited set of data for these parameters,a distinct seasonal pattern was not apparent(Table 1 and references therein).However,salinity was lower in El Ni?o years due to accelerated melting of Antarctic ice shelves and sea ice than in La Ni?a or normal years.This is similar to the situation in the central and eastern tropical Pacific where salinity decreases due to increased precipitation during El Ni?o years(McPhadenet al.,2006).

    3.2 Monthly Remote Sensing Data of Chl-a,POC and SST(2007-2011),and Their Response to ENSO Events

    Monthly remote sensing observations of Chl-a,POC and SST during 2007 – 2011(austral summer,November– March)were employed to estimate the relationship between phytoplankton production and ENSO events(MEI index,http://www.esrl.noaa.gov/psd/enso/mei/)in Prydz Bay(72? ? 78?E,67 ? 69?S).Chl-aand the other environmental factors have been shown to be sensitive to climate change in Antarctica(Loebet al.,2009).In our study phytoplankton bloom areas and Chl-aconcentrations were different due to ENSO-related SST variation,sea ice variability,and nutrients supply changes between years.This was particularly obvious in polynya regions,which also supports previous work(Arrigo and van Dijken,2003; Tremblay and Smith,2007).During La Ni?a years(shown in cyan color).Phytoplankton blooms occurred in January(both Chl-aand POC were higher than in December and were similar to those in February),and lagged behind the elevated SST.During the El Ni?o year(shown in light red color),phytoplankton blooms occurred in December(both Chl-aand POC were higher than in January and February)and preceded the elevated SST(Fig.3).These results indicated that the relationships between variability of Chl-a,POC and SST probably arose from changes in temperature,stratification and nutrients(e.g.,Fe)in upper Southern Ocean which were mainly ENSO-driven(Arrigoet al.,2008).The mechanism for this kind of relationships still remains unknown,with further research needed to better explore the relationship between ENSO,marine phytoplankton composition and abundance,and SST.

    Fig.3 Relationship between remote sensing monthly average concentrations of Chl-a(mg m?3),POC(mg m?3)as well as SST(℃)and Multivariate ENSO Index(MEI,monthly)(http://www.esrl.noaa.gov/psd/enso/mei/)in Prydz Bay(72?–78?E,67?–69?S)during austral summer(Nov.– Mar.),2007 – 2011.El Ni?o event(Jun.2009 – Apr.2010)is in light red color and La Ni?a events(Sept.2007 – Apr.2008,Oct.2008 – Mar.2009 and Jul.2010 – Apr.2011)are in cyan color.

    Both Chl-aand SST in Prydz Bay were significantly higher in austral summer during El Ni?o(shown in light red color)compared to La Ni?a years(shown in cyan color).For example,the highest Chl-aconcentration in an El Ni?o year(that in December)was 1.7 to 2.0 folds higher than that in La Ni?a years(January).However,POC did not vary significantly between El Ni?o and La Ni?a years(Fig.3).It indicated that phytoplankton biomass in Prydz Bay was impacted significantly by different ENSO events,and that the timing of phytoplankton blooms and SST was quite different between warm and cold phases of climate change pattern(Loebet al.,2009; Philander and Fedorov,2003; Yuan,2004).

    3.3 Long-Term Marine Lipid Biomarker Stratigraphic Records(1900–2005)of Phytoplankton,and Their Response to ENSO Events

    Fig.4 Excess 210Pb concentrations(dpm g?1)in sediment core III-12 from Prydz Bay,East Antarctica,collected in Feb.2005 during CHINARE-21.The sediment accumulation rate is 0.17 cm yr?1(R2 = 0.90).

    Based on the down-core distribution of210Pbxsactivity in sediment core III-12(Fig.4),the sediment accumulation rate was 0.17 cm yr?1.Thus,each centimeter of the core represented approximately 5.9 years since there was no surface mixed-layer.Accordingly,the top of core III-12 should represent the year 2005,and a depth of 18 cm the year 1899.We have shown here that the long-term variation of marine phytoplankton community structure,with climate change in Antarctica,can be effectively examined by employing historical trends of total concentration and relative proportion of diagnostic lipid biomarker stratigraphic records of core III-12 and MEI.ext index(Fig.5).Over the past century,marine phytoplankton composition and abundance,as indicated by total concentration of sedimentary brassicasterol(diatoms),dinosterol(dinoflagellates)and C37alkenones(haptophytes),showed increasing trends.The relative proportion of diatoms was the highest,followed by dinoflagellates,and haptophytes being the lowest.Previous work has shown that Antarctic phytoplankton community structure varied due to global climate change(Arrigoet al.,1999).Diatoms(low ratio of C to nutrient drawdown)are sensitive to SST variation,and therefore their overall abundance increased significantly with SST increase and was higher than that ofPhaeocystis antarctica(high ratio of C to nutrient drawdown)(Arrigoet al.,1999).Warm or cold climate changes have been shown to impact the utilization of C and atmosphericpCO2(Marinovet al.,2006).The shift in dominant phytoplankton fromP.antarcticato diatoms inferred that the Antarctic marine ecosystem in this bay is showing a response to climate changes.These changes at the level of primary production will likely have significant impacts at higher trophic levels.In fact,recent research on lipid biomarkers in penguin droppings from the Ardley Island,West Antarctica,showed that ENSO impacted significantly on lower plants and soil microorganism as well as penguin community in the Antarctic terrestrial environment(Zhanget al.,2012).

    Fig.5 Bulk biomarker concentrations and relative concentrations of individual biomarker in each layer of core III-12,and their relationship with ENSO(indicated by monthly MEI.ext,Wolter and Timlin(2011)).Strong El Ni?o events(May 1997– May 1998,Sept.1982 – Sept.1983,Jun.1972 – Mar.1973 and Apr.1957 – Aug.1958)are in light red color and strong La Ni?a events(Dec.1998 – Feb.2000,Jun.1988 – Jun.1989,May 1975 – Mar.1976 and Oct.1955 – Jul.1956)are in cyan color.

    Variations of brassicasterol,dinosterol percentages,brassicasterol:dinosterol ratio,and MEI.ext index showed that phytoplankton community structure changes(e.g.,percentage of diatoms in phytoplankton)were affected by ENSO.When examining downcore changes in such lipids it is also important to consider the effects of remineralization on biomarkers during early diagenesis(Sicreet al.,2000; Volkman,2006).As shown in Fig.5,the relatively high resolution stratigraphic record of lipid biomarker in the top 10 cm of core III-12 revealed similar trends with MEI.ext index.This supports the notion that these biomarkers were stable throughout this section of the core.Based on these lipids,we see that strong El Ni?o(shown in light red color)and La Ni?a(shown in cyan color) events have been recorded in sediments by these biomarkers by high(86.5%,85.1%,83.1% and 78.1%)and low(70.6%,72.5%,81.2% and 62.0%)percentages of sedimentary brassicasterol,respectively.It is important to note that the sedimentary biomarker concentration in each layer was an integration of ca.6 years.Further corroboration of the veracity of the biomarker data is that the annual field data counts agreed with the percentages of diatoms based on biomarker abundance in sediments(Table 1).This change in diatoms likely resulted from differences in ENSO-driven marine environment variability,e.g.,surface water stratification(SST-derived).For example,surface water stratification was more intensified during El Ni?o than during La Ni?a events,since SST was higher(see Fig.3).As stated earlier,the higher abundance of diatoms during El Ni?o(Arrigoet al.,1999)was confirmed by a higher concentration of brassicasterol in the sediment record(Figs.5B,E).It should be noted that there was clearly a lower resolution of sedimentary records at bottom of the core,making the ecological signals of El Ni?o and La Ni?a indistinguishable.Thus,there is a need for higher-resolution sediment sampling in future studies in this region.In general,our results not only infer that diatoms blooms occur during El Ni?o rather than La Ni?a events,but also show that the blooms had an increasing trend with time during the past century.We can only speculate at this time that these changes have occurred from a combination of ENSO and global climate warming(e.g.,intensified surface ocean stratification)(Arrigoet al.,1999; McPhadenet al.,2006).

    4 Conclusions

    Field measured and remote sensing data as well as sedimentary biomarker records provided comprehensive evidences of both short- and long-term impacts of ENSO on phytoplankton production and community structure in Prydz Bay,East Antarctica:

    1)Annual field measured phytoplankton community showed significant ENSO-related succession during 1990 to 2002.Algae species decreased significantly during El Ni?o/La Ni?a years compared to normal years.

    2)Monthly variation of remote sensing Chl-aand SST indicated that ENSO impacted the timing of phytoplankton blooms during 2007 to 2011.Phytoplankton blooms(indicated by Chl-a)preceded increases in SST in El Ni?o years,while they lagged behind increases in SST in La Ni?a years.

    3)Long-term stratigraphic lipid biomarker records inferred that the proportion of different algae changed significantly between El Ni?o and La Ni?a events.The relative proportion of diatoms increased with dinoflagellates being decreased during El Ni?o,but the reverse occurred during La Ni?a.

    Acknowledgements

    The authors wish to thank team members of CHINAREs,captains and the crew ofR/V Xuelongfor helping collect samples.This study was financially supported by the National Natural Science Foundation of China(NSFC)(40876104,41306202,41376193,41076134 and 41006118),the scientific research fund of Second Institute of Oceanography,SOA(JT1208 and JG1218),Chinese Arctic and Antarctic Administration Foundation(20110208),and the special fund for polar environment comprehensive investigation and assessment(CHINARE 2014-04-04,2014-01-04 and 2014-04-01).

    Arrigo,K.R.,and van Dijken,G.L.,2003.Phytoplankton dynamics within 37 Antarctic coastal polynya systems.Journal of Geophysical Research,108(C8):3271,DOI:10.1029/2002JC001739.

    Arrigo,K.R.,and van Dijken,G.L.,2004.Annual changes in sea-ice,chlorophylla,and primary production in the Ross Sea,Antarctica.Deep Sea Research II,51(1-3):117-138.

    Arrigo,K.R.,Robinson,D.H.,Worthen,D.L.,Dunbar,R.B.,DiTullio,G.R.,van Woert,M.,and Lizotte,M.P.,1999.Phytoplankton community structure and the drawdown of nutrients and CO2in the Southern Ocean.Science,283(5400):365-367.

    Arrigo,K.R.,van Dijken,G.L.,and Bushinsky,S.,2008.Primary production in the Southern Ocean,1997?2006.Journal of Geophysical Research,113(C8):C08004,DOI:10.01029/02007jc004551.

    Behrenfeld,M.J.,O’Malley,R.T.,Siegel,D.A.,McClain,C.R.,Sarmiento,J.L.,Feldman,G.C.,Milligan,A.J.,Falkowski,P.G.,Letelier,R.M.,and Boss,E.S.,2006.Climate-driven trends in contemporary ocean productivity.Nature,444(7120):752-755.

    Bianchi,T.S.,and Canuel,E.A.,2011.Chemical Biomarkers in Aquatic Ecosystems.Princeton University Press,Princeton,396pp.

    Cai,Y.,Ning,X.,Zhu,G.,and Shi,J.,2003.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 1998/1999.Acta Oceanologica Sinica,22(4):651-670.

    Chang,P.,Zhang,L.,Saravanan,R.,Vimont,D.J.,Chiang,J.C.H.,Ji,L.,Seidel,H.,and Tippett,M.K.,2007.Pacific meridional mode and El Ni?o-Southern Oscillation.Geophysical Research Letters,34(16):L16608,DOI:10.11029/12007GL030302.

    Chavez,F.P.,Messié,M.,and Pennington,J.T.,2011.Marine primary production in relation to climate variability and change.Annual Review of Marine Science,3(1):227-260.

    Chen,N.,Bianchi,T.S.,McKee,B.A.,and Bland,J.M.,2001.Historical trends of hypoxia on the Louisiana shelf:Application of pigments as biomarkers.Organic Geochemistry,32(4):543-561.

    Doney,S.C.,Ruckelshaus,M.,Emmett Duffy,J.,Barry,J.P.,Chan,F.,English,C.A.,Galindo,H.M.,Grebmeier,J.M.,Hollowed,A.B.,Knowlton,N.,Polovina,J.,Rabalais,N.N.,Sydeman,W.J.,and Talley,L.D.,2012.Climate change impacts on marine ecosystems.Annual Review of Marine Science,4(1):11-37.

    Hernández-de-la-Torre,B.,Gaxiola-Castro,G.,Alvarez-Borrego,S.,Gómez-Valdés,J.,and Nájera-Mart??nez,S.,2003.Interannual variability of new production in the southern region of the California Current.Deep Sea Research II,50(14-16): 2423-2430.

    Liu,C.,Ning,X.,Sun,J.,Cai,Y.,and Cao,J.,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 Oceanologica Sinica,26(6):107-117(in Chinese with English abstract).

    Liu,Z.,Cai,Y.,Ning,X.,Liu,C.,Zhu,G.,and Wang,X.,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.Chinese Journal of Polar Research,13(1):1-12(in Chinese with English abstract).

    Loeb,V.J.,Hofmann,E.E.,Klinck,J.M.,Holm-Hansen,O.,and White,W.B.,2009.ENSO and variability of the Antarctic Peninsula pelagic marine ecosystem.Antarctic Science,21(2):135-148.

    Mackey,K.R.M.,van Dijken,G.L.,Mazloom,S.,Erhardt,A.M.,Ryan,J.,Arrigo,K.R.,and Paytan,A.,2010.Influence of atmospheric nutrients on primary productivity in a coastal upwelling region.Global Biogeochemical Cycles,24(4):GB4027,DOI:10.1029/2009GB003737.

    Marinov,I.,Gnanadesikan,A.,Toggweiler,J.R.,and Sarmiento,J.L.,2006.The Southern Ocean biogeochemical divide.Nature,441(7096):964-967.

    McPhaden,M.J.,Zebiak,S.E.,and Glantz,M.H.,2006.ENSO as an integrating concept in Earth science.Science,314(5806):1740-1745.

    Ning,X.,1998.Phytoplankton Standing Stock and Primary Productivity in Specific Areas of the Southern Ocean.Ocean Press,Beijing,98pp(in Chinese).

    Nittrouer,C.A.,and Sternberg,R.W.,1981.The formation of sedimentary strata in an allochthonous shelf environment:The Washington continental shelf.Marine Geology,42(1-4):201-232.

    Philander,S.G.,and Fedorov,A.,2003.Is El Ni?o Sporadic or Cyclic?Annual Review of Earth and Planetary Sciences,31(1):579-594.

    Pu,S.,and Dong,Z.,2003.Progress in physical oceanographic studies of Prydz Bay and its adjacent oceanic area.Chinese Journal of Polar Research,15(1):53-64(in Chinese with English abstract).

    Reiss,C.S.,Hewes,C.D.,and Holm-Hansen,O.,2009.Influence of atmospheric teleconnections and Upper Circumpolar Deep Water on phytoplankton biomass around Elephant Island,Antarctica.Marine Ecology Progress Series,377:51-62.

    Rittenour,T.M.,Brigham-Grette,J.,and Mann,M.E.,2000.El Ni?o-like climate teleconnections in New England during the late Pleistocene.Science,288(5468):1039-1042.

    Sarmiento,J.L.,Gruber,N.,Brzezinski,M.A.,and Dunne,J.P.,2004.High-latitude controls of thermocline nutrients and low latitude biological productivity.Nature,427(6969):56-60.

    Sicre,M.A.,Ternois,Y.,Paterne,M.,Boireau,A.,Beaufort,L.,Martinez,P.,and Bertrand,P.,2000,Biomarker stratigraphic records over the last 150 kyears off the NW African coast at 25?N.Organic Geochemistry,31(6):577-588.

    Smetacek,V.,and Nicol,S.,2005.Polar ocean ecosystems in a changing world.Nature,437(7057):362-368.

    Smith,N.R.,Dong,Z.,Kerry,K.R.,and Wright,S.,1984.Water masses and circulation in the region of Prydz Bay,Antarctica.Deep Sea Research I,31(9):1121-1147.

    Stein,R.,and Macdonald,R.W.,2004.The Organic Carbon Cycle in the Arctic Ocean.Springer,Berlin,363pp.

    Tremblay,J.E.,and Smith,W.O.,2007.Primary production and nutrient dynamics in polynyas.In:Polynyas:Windows to the World.Smith,W.O.,and Barber,D.G.,eds.,Elsevier,Amsterdam,239-269.

    Trenberth,K.E.,and Hoar,T.J.,1997.El Ni?o and climate change.Geophysical Research Letters,24(23):3057-3060,DOI:10.1029/3097GL03092.

    Volkman,J.K.,2006.Lipid markers for marine organic matter.In:Marine Organic Matter:Biomarkers,Isotopes and DNA.Volkman,J.K.,ed.,Springer,Berlin/Heidelberg,27-70.

    Wolter,K.,and Timlin,M.S.,1983.Monitoring ENSO in COADS with a seasonally adjusted principal component index.Proceeding of the 17th Climate Diagnostics Workshop,52-57.

    Wolter,K.,and Timlin,M.S.,2011.El Ni?o/Southern Oscillation behaviour since 1871 as diagnosed in an extended multivariate ENSO index(MEI.ext).International Journal of Climatology,31(7):1074-1087.

    Wright,S.W.,and van den Enden,R.L.,2000.Phytoplankton community structure and stocks in the East Antarctic marginal ice zone(BROKE survey,January?March 1996)determined by CHEMTAX analysis of HPLC pigment signatures.Deep Sea Research II,47(12-13):2363-2400.

    Wright,S.W.,van den Enden,R.L.,Pearce,I.,Davidson,A.T.,Scott,F.J.,and Westwood,K.J.,2010.Phytoplankton community structure and stocks in the Southern Ocean(30?-80?E)determined by CHEMTAX analysis of HPLC pigment signatures.Deep Sea Research II,57(9-10):758-778.

    Yuan,X.,2004.ENSO-related impacts on Antarctic sea ice:A synthesis of phenomenon and mechanisms.Antarctic Science,16(4):415-425.

    Zhang,H.,Lu,D.,Yu,P.,Zhang,W.,Lu,B.,Peter,H.U.,and Vetter,W.,2012.Biomarker records in penguin droppings and observed changes in penguin communities and their response to the ENSO in the Western Antarctic.Science China:Earth Sciences,55(8):1238-1247.

    Zhu,G.,Liu,Z.,Ning,X.,and Wang,M.,1994.Distribution characteristics of planktonic nano-and microalgae in the Prydz Bay and its adjacent Southern Indian Ocean,during austral summer.Antarctic Research,5(2):33-44.

    Zhu,G.,Ning,X.,Liu,Z.,Cai,Y.,Liu,C.,and Pan,J.,2006.A study on phytoplankton in the Prydz Bay and its adjacent sea area,Antarctica during the austral summer of 2000.Acta Oceanologica Sinica,28(1):118-126(in Chinese with English abstract).

    国产真实伦视频高清在线观看 | 午夜激情欧美在线| 一进一出好大好爽视频| 国产亚洲精品久久久久久毛片| 黄色丝袜av网址大全| 无人区码免费观看不卡| 99riav亚洲国产免费| 噜噜噜噜噜久久久久久91| 88av欧美| 亚洲国产精品sss在线观看| 两个人的视频大全免费| 白带黄色成豆腐渣| 色尼玛亚洲综合影院| 中文字幕精品亚洲无线码一区| 一个人免费在线观看电影| 欧美bdsm另类| 国产精品影院久久| 三级男女做爰猛烈吃奶摸视频| 欧美最新免费一区二区三区 | 欧美不卡视频在线免费观看| 精品久久久久久久久久久久久| 在线a可以看的网站| 午夜福利成人在线免费观看| 嫁个100分男人电影在线观看| 嫩草影院精品99| 天堂√8在线中文| 亚洲av成人不卡在线观看播放网| 欧美日韩综合久久久久久 | 成年人黄色毛片网站| 毛片一级片免费看久久久久 | 日韩国内少妇激情av| 国产中年淑女户外野战色| 性插视频无遮挡在线免费观看| 综合色av麻豆| 精品午夜福利在线看| 久久这里只有精品中国| 日本五十路高清| 国产极品精品免费视频能看的| 国产黄片美女视频| 99国产精品一区二区三区| 欧美最黄视频在线播放免费| 日本免费a在线| ponron亚洲| 亚洲人成网站在线播放欧美日韩| av黄色大香蕉| 十八禁人妻一区二区| 国产欧美日韩一区二区三| 在线天堂最新版资源| 日日摸夜夜添夜夜添小说| 特大巨黑吊av在线直播| 国产精品爽爽va在线观看网站| 国产欧美日韩一区二区三| 国产精品乱码一区二三区的特点| 99久国产av精品| 最近最新免费中文字幕在线| 黄色女人牲交| 欧美日韩黄片免| 日日夜夜操网爽| 国产成人影院久久av| 最近中文字幕高清免费大全6 | 狠狠狠狠99中文字幕| 亚洲自拍偷在线| 一进一出抽搐gif免费好疼| 色综合婷婷激情| 老女人水多毛片| 真实男女啪啪啪动态图| 国产av一区在线观看免费| 国产成人aa在线观看| 桃红色精品国产亚洲av| 国产高清激情床上av| 国产精品综合久久久久久久免费| 黄色配什么色好看| 日本黄色视频三级网站网址| 亚洲精品456在线播放app | 全区人妻精品视频| 久久久久精品国产欧美久久久| 国产视频内射| 亚洲av第一区精品v没综合| 久久中文看片网| 深夜精品福利| 村上凉子中文字幕在线| 简卡轻食公司| 色综合亚洲欧美另类图片| 国产亚洲精品久久久久久毛片| 18+在线观看网站| 特大巨黑吊av在线直播| 亚洲精华国产精华精| 亚洲最大成人av| 色吧在线观看| 欧美在线黄色| 国产精品综合久久久久久久免费| 婷婷亚洲欧美| av专区在线播放| 天堂影院成人在线观看| 精品日产1卡2卡| 欧美中文日本在线观看视频| 日本黄色片子视频| 久久久久性生活片| 亚洲欧美精品综合久久99| 亚洲欧美激情综合另类| 中国美女看黄片| 欧美性猛交╳xxx乱大交人| 国产黄a三级三级三级人| 91麻豆精品激情在线观看国产| 一a级毛片在线观看| 男人舔奶头视频| av天堂中文字幕网| 亚洲一区二区三区色噜噜| 国产欧美日韩一区二区三| 两个人的视频大全免费| 一进一出抽搐gif免费好疼| 免费av观看视频| 可以在线观看毛片的网站| 国产久久久一区二区三区| 国产麻豆成人av免费视频| 国产成人av教育| 老司机午夜十八禁免费视频| 国产亚洲精品久久久com| 国产欧美日韩精品亚洲av| avwww免费| 国内揄拍国产精品人妻在线| 免费高清视频大片| 精品久久久久久久久亚洲 | 久久久精品欧美日韩精品| 91字幕亚洲| 亚洲精品一卡2卡三卡4卡5卡| 真人做人爱边吃奶动态| 亚洲久久久久久中文字幕| 亚洲精品粉嫩美女一区| 亚洲国产高清在线一区二区三| 国产毛片a区久久久久| 免费一级毛片在线播放高清视频| 国产精品一区二区免费欧美| 亚洲第一电影网av| 色av中文字幕| 国产淫片久久久久久久久 | 丁香欧美五月| 三级男女做爰猛烈吃奶摸视频| 午夜日韩欧美国产| 高清在线国产一区| 人妻久久中文字幕网| 日本撒尿小便嘘嘘汇集6| 国产精品爽爽va在线观看网站| 搡女人真爽免费视频火全软件 | 真人做人爱边吃奶动态| 日本一本二区三区精品| 女人被狂操c到高潮| 国产伦人伦偷精品视频| 噜噜噜噜噜久久久久久91| 久久久久久久久大av| 热99在线观看视频| 亚洲国产精品999在线| 精品人妻1区二区| 国内精品美女久久久久久| 久久香蕉精品热| 一区二区三区激情视频| av在线老鸭窝| 久久99热6这里只有精品| 真人一进一出gif抽搐免费| 亚洲熟妇中文字幕五十中出| 国产精品,欧美在线| 色综合亚洲欧美另类图片| 丰满人妻一区二区三区视频av| 久久中文看片网| 日韩欧美一区二区三区在线观看| 久久6这里有精品| 亚洲国产欧美人成| 午夜福利免费观看在线| a级毛片a级免费在线| 国产三级黄色录像| 五月玫瑰六月丁香| 国产高清激情床上av| 国产精品av视频在线免费观看| 午夜福利欧美成人| 午夜影院日韩av| 精品久久久久久,| 久久人人精品亚洲av| 成人av一区二区三区在线看| 成年版毛片免费区| 久久精品国产亚洲av涩爱 | 国产蜜桃级精品一区二区三区| 大型黄色视频在线免费观看| 嫩草影视91久久| 国产三级在线视频| 午夜影院日韩av| 久久久成人免费电影| 此物有八面人人有两片| 精品人妻偷拍中文字幕| 我的老师免费观看完整版| 99久久久亚洲精品蜜臀av| 99久久久亚洲精品蜜臀av| 国产高潮美女av| 嫩草影视91久久| 夜夜看夜夜爽夜夜摸| 久久精品人妻少妇| 久久国产乱子伦精品免费另类| 国产三级在线视频| 亚洲成人中文字幕在线播放| 最新中文字幕久久久久| 国产蜜桃级精品一区二区三区| 午夜福利18| 热99在线观看视频| 欧美精品国产亚洲| 欧美性猛交╳xxx乱大交人| 免费电影在线观看免费观看| 久久亚洲精品不卡| 日本 av在线| 88av欧美| 国产精品自产拍在线观看55亚洲| 香蕉av资源在线| 亚洲av美国av| www.999成人在线观看| 久久精品国产亚洲av香蕉五月| 国产伦精品一区二区三区视频9| 久久久成人免费电影| 亚洲精品色激情综合| 少妇的逼好多水| 久久久久久久久久黄片| 欧美日韩中文字幕国产精品一区二区三区| 国产黄a三级三级三级人| 九九热线精品视视频播放| 最近中文字幕高清免费大全6 | 亚洲精品日韩av片在线观看| 精品人妻1区二区| 国产精品av视频在线免费观看| 十八禁网站免费在线| 中文字幕熟女人妻在线| 性色av乱码一区二区三区2| 免费av毛片视频| 天堂动漫精品| 51午夜福利影视在线观看| 午夜福利在线观看吧| 中亚洲国语对白在线视频| 午夜福利视频1000在线观看| 麻豆成人午夜福利视频| 亚洲欧美激情综合另类| 最近在线观看免费完整版| 毛片女人毛片| 51国产日韩欧美| 级片在线观看| 麻豆一二三区av精品| 日日摸夜夜添夜夜添av毛片 | 性欧美人与动物交配| 亚洲在线观看片| 亚洲av中文字字幕乱码综合| 日日摸夜夜添夜夜添小说| 精品久久久久久久久久免费视频| netflix在线观看网站| 在线a可以看的网站| 波野结衣二区三区在线| 亚洲第一区二区三区不卡| 91九色精品人成在线观看| 国内精品一区二区在线观看| 成年女人看的毛片在线观看| 91午夜精品亚洲一区二区三区 | 国产熟女xx| 婷婷精品国产亚洲av| 人人妻人人澡欧美一区二区| 免费av毛片视频| av视频在线观看入口| 亚洲黑人精品在线| 欧美乱妇无乱码| 美女高潮的动态| 久久伊人香网站| 国产精品人妻久久久久久| 亚洲专区国产一区二区| av在线观看视频网站免费| 国产精品,欧美在线| 超碰av人人做人人爽久久| 久久午夜福利片| 在线看三级毛片| 一本精品99久久精品77| 欧美日韩国产亚洲二区| 亚洲成人久久性| 欧美在线黄色| 国产三级黄色录像| 精品久久久久久久末码| 一本一本综合久久| 在线免费观看不下载黄p国产 | 日日干狠狠操夜夜爽| 亚洲国产色片| 精品人妻熟女av久视频| 色综合欧美亚洲国产小说| 日韩欧美国产在线观看| 天堂av国产一区二区熟女人妻| 亚洲欧美日韩东京热| 黄色一级大片看看| 午夜福利在线观看免费完整高清在 | 欧美激情国产日韩精品一区| 757午夜福利合集在线观看| 亚洲一区二区三区色噜噜| 最新中文字幕久久久久| 淫秽高清视频在线观看| 午夜免费激情av| 欧美黑人巨大hd| 很黄的视频免费| 亚洲av成人精品一区久久| 欧美日韩国产亚洲二区| 国产精品一区二区性色av| 美女免费视频网站| 欧美一区二区国产精品久久精品| 女人十人毛片免费观看3o分钟| 国产真实乱freesex| 国产成人啪精品午夜网站| 国产欧美日韩精品一区二区| 色哟哟哟哟哟哟| 欧美在线黄色| 国产一区二区亚洲精品在线观看| 久久久久国内视频| 亚洲男人的天堂狠狠| 精品人妻一区二区三区麻豆 | 色5月婷婷丁香| 午夜视频国产福利| 国产精品98久久久久久宅男小说| 国产高清视频在线播放一区| 亚洲美女视频黄频| 极品教师在线免费播放| 69人妻影院| 国产视频内射| 成人一区二区视频在线观看| 国产乱人伦免费视频| 夜夜看夜夜爽夜夜摸| 亚洲乱码一区二区免费版| 亚洲熟妇中文字幕五十中出| 成人午夜高清在线视频| 美女 人体艺术 gogo| 亚洲av电影不卡..在线观看| 99久久99久久久精品蜜桃| 身体一侧抽搐| 九九久久精品国产亚洲av麻豆| 亚洲无线在线观看| 亚洲五月婷婷丁香| 欧美中文日本在线观看视频| 99久久精品一区二区三区| 国产精品一区二区性色av| 宅男免费午夜| 日韩av在线大香蕉| 校园春色视频在线观看| 亚洲,欧美,日韩| 国产精品影院久久| 亚州av有码| 免费人成视频x8x8入口观看| 在线看三级毛片| 人妻丰满熟妇av一区二区三区| 看十八女毛片水多多多| 国产不卡一卡二| 亚洲av第一区精品v没综合| aaaaa片日本免费| 亚洲av中文字字幕乱码综合| 国产极品精品免费视频能看的| 99国产极品粉嫩在线观看| 色哟哟·www| 真人一进一出gif抽搐免费| 精品久久久久久成人av| 免费一级毛片在线播放高清视频| 色哟哟哟哟哟哟| 久久这里只有精品中国| 岛国在线免费视频观看| 国产视频内射| 波多野结衣巨乳人妻| 国产高清视频在线播放一区| 国产精品,欧美在线| 男插女下体视频免费在线播放| 中文字幕人妻熟人妻熟丝袜美| 搞女人的毛片| 亚洲国产精品sss在线观看| 亚洲av.av天堂| 五月玫瑰六月丁香| 亚洲自偷自拍三级| 亚洲精品粉嫩美女一区| 久久久久国内视频| 亚洲av中文字字幕乱码综合| 亚洲真实伦在线观看| 亚洲aⅴ乱码一区二区在线播放| 国产成人啪精品午夜网站| 亚洲人成网站在线播放欧美日韩| 国产精品一区二区三区四区久久| 精品熟女少妇八av免费久了| av专区在线播放| 此物有八面人人有两片| 国产精品精品国产色婷婷| 丝袜美腿在线中文| 久久久久性生活片| 中文字幕久久专区| a级一级毛片免费在线观看| 久久热精品热| 欧美精品国产亚洲| 午夜精品在线福利| 直男gayav资源| 日本熟妇午夜| av在线观看视频网站免费| 一个人免费在线观看的高清视频| 18禁黄网站禁片午夜丰满| 国产 一区 欧美 日韩| 婷婷精品国产亚洲av| 黄色女人牲交| 欧美成人a在线观看| 亚洲不卡免费看| 亚洲va日本ⅴa欧美va伊人久久| 亚洲国产精品999在线| 国产黄色小视频在线观看| 男人狂女人下面高潮的视频| 免费av观看视频| 一区二区三区高清视频在线| 国产精品三级大全| 非洲黑人性xxxx精品又粗又长| 欧美一级a爱片免费观看看| 在线天堂最新版资源| 草草在线视频免费看| a在线观看视频网站| 国产大屁股一区二区在线视频| 老熟妇乱子伦视频在线观看| 少妇被粗大猛烈的视频| 淫秽高清视频在线观看| 男人和女人高潮做爰伦理| 啦啦啦观看免费观看视频高清| 中文字幕av成人在线电影| 国产精品爽爽va在线观看网站| 久99久视频精品免费| 看黄色毛片网站| 午夜福利在线在线| 亚洲性夜色夜夜综合| 成年免费大片在线观看| 中文字幕人成人乱码亚洲影| 亚洲熟妇熟女久久| 真人一进一出gif抽搐免费| 国产在视频线在精品| 欧美潮喷喷水| 国产精品,欧美在线| 亚洲成人久久爱视频| 国产精品乱码一区二三区的特点| 少妇丰满av| 午夜老司机福利剧场| 亚洲国产日韩欧美精品在线观看| 精品国产三级普通话版| bbb黄色大片| av在线蜜桃| 桃红色精品国产亚洲av| 丝袜美腿在线中文| 中文字幕av在线有码专区| 宅男免费午夜| 直男gayav资源| 欧美在线黄色| 欧美3d第一页| 老熟妇仑乱视频hdxx| 最后的刺客免费高清国语| 亚洲色图av天堂| 一级毛片久久久久久久久女| 小蜜桃在线观看免费完整版高清| 十八禁网站免费在线| 欧洲精品卡2卡3卡4卡5卡区| 一a级毛片在线观看| 99久国产av精品| 国产真实伦视频高清在线观看 | 亚洲男人的天堂狠狠| 日韩成人在线观看一区二区三区| 两人在一起打扑克的视频| 神马国产精品三级电影在线观看| 又爽又黄a免费视频| 精品人妻视频免费看| 国产精品久久电影中文字幕| 丰满人妻一区二区三区视频av| 啦啦啦观看免费观看视频高清| 一二三四社区在线视频社区8| 天堂av国产一区二区熟女人妻| 亚洲性夜色夜夜综合| 少妇熟女aⅴ在线视频| 午夜老司机福利剧场| 国产高清视频在线播放一区| 夜夜夜夜夜久久久久| 男人舔女人下体高潮全视频| 午夜福利成人在线免费观看| 97超视频在线观看视频| 中文字幕av在线有码专区| 91在线精品国自产拍蜜月| 久久国产精品影院| 成人av在线播放网站| av在线观看视频网站免费| 成人亚洲精品av一区二区| 国产欧美日韩一区二区三| 亚洲三级黄色毛片| av在线天堂中文字幕| 白带黄色成豆腐渣| 在现免费观看毛片| 亚洲精品一区av在线观看| 精品99又大又爽又粗少妇毛片 | 久久亚洲真实| 免费在线观看影片大全网站| 首页视频小说图片口味搜索| 色综合亚洲欧美另类图片| 国产av麻豆久久久久久久| 免费观看人在逋| 变态另类成人亚洲欧美熟女| 日韩有码中文字幕| 精品久久久久久久久亚洲 | 久久香蕉精品热| 99精品久久久久人妻精品| 国产野战对白在线观看| 精品久久久久久久久久久久久| 长腿黑丝高跟| 一个人看视频在线观看www免费| 亚洲精华国产精华精| 一区二区三区高清视频在线| 一级av片app| 一进一出好大好爽视频| 又黄又爽又刺激的免费视频.| 嫩草影院精品99| 网址你懂的国产日韩在线| 欧美色视频一区免费| 色在线成人网| 午夜久久久久精精品| 丰满乱子伦码专区| 看十八女毛片水多多多| 成人性生交大片免费视频hd| 亚洲成人中文字幕在线播放| 日本一本二区三区精品| 黄色视频,在线免费观看| 亚洲av二区三区四区| 久久性视频一级片| 亚洲乱码一区二区免费版| 欧美国产日韩亚洲一区| 亚洲国产欧美人成| 国内精品一区二区在线观看| 午夜精品在线福利| 亚洲欧美日韩无卡精品| 最近最新中文字幕大全电影3| 亚洲无线在线观看| 国产白丝娇喘喷水9色精品| av女优亚洲男人天堂| 久久国产乱子免费精品| 亚洲在线自拍视频| 少妇裸体淫交视频免费看高清| 久久香蕉精品热| 久久99热6这里只有精品| 一进一出好大好爽视频| 中文资源天堂在线| 在线观看午夜福利视频| 久久欧美精品欧美久久欧美| 特大巨黑吊av在线直播| 亚洲av成人av| 99国产精品一区二区三区| 欧美乱妇无乱码| 日韩成人在线观看一区二区三区| 欧美一级a爱片免费观看看| 超碰av人人做人人爽久久| 夜夜看夜夜爽夜夜摸| 国产视频内射| 国产精品,欧美在线| 亚洲一区二区三区不卡视频| 十八禁人妻一区二区| 午夜免费成人在线视频| 真实男女啪啪啪动态图| 我要搜黄色片| 色哟哟哟哟哟哟| 欧美日韩亚洲国产一区二区在线观看| АⅤ资源中文在线天堂| 成人国产综合亚洲| 国产成人aa在线观看| avwww免费| 日日摸夜夜添夜夜添小说| 久久亚洲真实| 99国产极品粉嫩在线观看| 亚洲第一电影网av| 乱人视频在线观看| 午夜两性在线视频| 麻豆av噜噜一区二区三区| 美女黄网站色视频| 免费看光身美女| av女优亚洲男人天堂| 亚洲美女视频黄频| 亚洲午夜理论影院| 亚洲五月婷婷丁香| 亚洲成人中文字幕在线播放| 国产av麻豆久久久久久久| 午夜福利在线观看吧| av在线观看视频网站免费| 成人特级av手机在线观看| av福利片在线观看| 两人在一起打扑克的视频| 亚洲18禁久久av| 757午夜福利合集在线观看| 日韩 亚洲 欧美在线| 99久久精品热视频| 日日摸夜夜添夜夜添av毛片 | 国产色爽女视频免费观看| 天美传媒精品一区二区| 三级国产精品欧美在线观看| 人妻久久中文字幕网| 国产亚洲欧美在线一区二区| 波野结衣二区三区在线| 久久99热6这里只有精品| 欧美精品国产亚洲| 成人国产综合亚洲| 亚洲欧美激情综合另类| 免费av毛片视频| 欧美高清成人免费视频www| 青草久久国产| av国产免费在线观看| 国产欧美日韩精品一区二区| 欧美xxxx黑人xx丫x性爽| 欧美日本亚洲视频在线播放| 午夜精品一区二区三区免费看| 亚洲国产精品成人综合色| av在线观看视频网站免费| 亚洲在线自拍视频| 91在线精品国自产拍蜜月| 日韩精品青青久久久久久| 乱码一卡2卡4卡精品| 日本免费a在线| 每晚都被弄得嗷嗷叫到高潮| 日本在线视频免费播放| 观看美女的网站| 日本免费一区二区三区高清不卡| 丰满人妻熟妇乱又伦精品不卡| 我要搜黄色片|