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

    The Congo Basin Zonal Overturning Circulation

    2016-11-14 05:56:13DepartmentofGeologicalSciencesJacksonSchoolofGeosciencesTheUniversityofTexasatAustinAustinTexas78712USA
    Advances in Atmospheric Sciences 2016年6期

    Department of Geological Sciences,Jackson School of Geosciences,The University of Texas at Austin, Austin,Texas 78712,USA

    (Received 28 August 2015;revised 19 November 2015;accepted 30 November 2015)

    The Congo Basin Zonal Overturning Circulation

    Naresh NEUPANE×

    Department of Geological Sciences,Jackson School of Geosciences,The University of Texas at Austin, Austin,Texas 78712,USA

    (Received 28 August 2015;revised 19 November 2015;accepted 30 November 2015)

    The Gulf of Guinea in the equatorial Atlantic is characterized by the presence of strong subsidence at certain times of the year.This subsidence appears in June and becomes well established from July to September.Since much of the West African monsoon flow originates over the Gulf,Guinean subsidence is important for determining moisture sources for the monsoon. Using reanalysis products,I contribute to a physical understanding of what causes this seasonal subsidence,and how it relates to precipitation distributions across West Africa.

    There is a seasonal zonal overturning circulation above the Congo basin and the Gulf of Guinea in the ERA-Interim, ERA-40,NCEP2,and MERRA reanalyses.The up-branch is located in the Congo basin around 20°E.Mid-tropospheric easterly flows constitute the returning-branch and sinking over the Gulf of Guinea forms the down-branch,which diverges at 2°W near the surface,with winds to the east flowing eastward to complete the circulation.This circulation is driven by surface temperature differences between the eastern Gulf and Congo basin.Land temperatures remain almost uniform, around 298 K,throughout a year,but the Guinean temperatures cool rapidly from 294 K in May to about 290 K in August. These temperature changes increase the ocean/land temperature contrast,up to 8 K,and drive the circulation.

    I hypothesize that when the overturning circulation is anomalously strong,the northward moisture transport and Sahelian precipitation are also strong.This hypothesis is supported by ERA-Interim and PERSIANN-CDR(Precipitation Estimation from Remotely Sensed Information using Artificial Neural Networks-Climate Data Record)data.

    climate,atmospheric dynamics,monsoon

    1.Introduction

    Past studies(Trenberth et al.,2000;Wang,2002a,2002b, 2005;Nicholson and Webster,2007;Pokam et al.,2014)indicate that atmospheric subsidence exists over the equatorial Atlantic,with the strongest sinking occurringover the central andeastern equatorialAtlantic(i.e.,overthe Gulf ofGuinea). However,the relationship between the Gulf of Guinea subsidence and local and regional climate variability over continental Africa is not yet well explored.For example,the Gulf of Guinea is known to be a primary moisture source for the West African monsoon system(Vizy and Cook,2001; Grist and Nicholson,2001;Fontaine et al.,2003;Cook and Vizy,2006),hencefluctuationsinsubsidencemayaffectrainfall patterns over the adjacent African continent by altering the low-level circulation and moisture convergence.Given Africa’s strong reliance on rainfall for agricultural needs,it is important to improve our understanding of any mechanism(s)that influence rainfall variations over equatorial and West Africa.The purpose of this study is to use availabledata to documentthe seasonal-scalesubsidence in the Gulf of Guinea,and improve our physical understanding of the relationship between this subsidence and boreal summer rainfall in West Africa.

    The low-level Guinean subsidence(Hastenrath,2001; Cook and Vizy,2006;Hastenrath and Polzin,2011)is often related to the rainfall over West Africa,including the onset of monsoon(Hagos and Cook,2007;Caniaux et al.,2011). Using reanalyses and satellite observations,Leduc-Leballeur et al.(2013)showed that increases in the Guineansubsidence are related to increases in convection over the northern Gulf of Guinea.In an observational study,Segele et al.(2009) suggested that intensification of the Guinean surface pressure is associated with abundant rainfall in northern Africa.Similarly,Nicholson and Webster(2007),using NCEP reanalysis,observed precipitation enhancement over the Sahel when the Guinean subsidence strengthened.Studies have also indicated that this subsidence inhibits deep vertical transport of moisture overthe Gulf and supportspolewardmoisture transport by the monsoon flow(Vizy and Cook,2001;Neupane and Cook,2013).

    The presence of a deep atmospheric meridional overturning circulation over the Atlantic,often referred to as“the At-lantic Hadley circulation”,has also been revealed,and the descending branch of this circulation,whose up-branchis located over the Sahel,has been shown to be related to the Guinean subsidence.Low-level flows from the equator toward the Sahel form the lower branch,while the equatorward flows above 600 hPa form the returning branch of the circulation.This circulation becomes most active during the boreal summer(Trenberth et al.,2000;Wang,2004;Wang, 2005).An example can be found in Wang(2002a),who used 1950–1999 NCEP reanalysis climatology and averaged tropospheric circulation profiles from 10°W to 10°E to show the subsiding branch of the Atlantic Hadley circulation.In addition,some other studies relate this subsidence to the descending branch of a shallow meridional overturning circulation that is embedded within the deep meridional overturning circulation.The shallow circulation consists of surface southerly inflows from the Guinean Gulf and the northerly returning outflows between 700 and 500 hPa from the Sahel as well(Thorncroftand Blackburn,1999;Zhang et al.,2006; Nolan et al.,2007;Hagos and Zhang,2010;Thorncroftet al., 2011).

    ?Institute of Atmospheric Physics/Chinese Academy of Sciences,and Science Press and Springer-Verlag Berlin Heidelberg 2016

    The deep and shallow meridional overturning circulations,as discussed above,are affected by Coriolis acceleration.As the onshore southerly winds flow northward onto the continent,theytendto deviateeastwardbecauseof theaction of zonal Coriolis acceleration.It is difficult for a purely meridional flow to be maintained.The Coriolis force appears as soon as an air parcel is deflected in a meridional direction and imparts zonal acceleration to the parcel.

    The zonal circulations along the equator,e.g.,Walker circulation,remainunaffectedbytheCoriolisacceleration.Vizy and Cook(2001)hypothesized that a Walker-type circulation exists above the Congo basin and the Gulf of Guinea. They used an AGCM forced with warm SST anomalies in the Gulf.Rainfall increased along the coast and decreased overthe Congo basin when the SSTs were anomalouslyhigh, and it was found to be associated with a weakening of the circulation.This circulation can also be seen in zonal and vertical wind streamlines along the equator in Lau and Yang (2003,Fig.2c).Theyusedthe 1949–1999NCEP climatology for July to show the Walker circulations around the globe.In addition,sketches indicating the presence of this circulation can be seen in Hastenrath(2006)and the references therein.

    A number of recent studies have indicated the existence of low-levelwesterly flows from the Gulf to the Congo basin. The circulation shows that the increased rainfall over the Congo basin is associated with anomalously stronger lowlevel westerly(LLW)flows(Pokam et al.,2012;Dezfuli and Nicholson,2013;Nicholson and Dezfuli,2013;Dezfuli et al.,2015);plus,the circulation is often associated with the northward flows and rainfall across West Africa,as discussed in Cook and Vizy(2015).These studies mostly focused on equatorial atmospheric circulation in the spring and fall.Pokam et al.(2014)investigated the seasonal variability and driver of the LLW flows using four reanalyses for the period 1989–2003.The LLW,driven by the land–ocean temperature contrast,forms the lower branch of the Walker-type circulation.This is weakest in spring and strongest in fall. This also suggests that the West African monsoonsystem and the atmospheric circulation over eastern equatorial Atlantic are closely connected.In this study,a zonal overturning circulation is identified over the Gulf of Guinea and the Congo basin.The variability of the circulation is examined based on an index obtained from the rising branch.This index is further used to understand the relationship between the circulation and rainfall and moisture distributions across West Africa.The focus is on the West African monsoon season, from July through September,when the overturning circulation is most pronounced.

    2.Data

    Since there is some degree of uncertainty in the observations,various observational and reanalysis datasets are analyzed and compared to each other to build confidence in the findings.Each dataset is discussed briefly below.

    Four reanalysis datasets are utilized to provide monthly mean information regardingthe circulation,temperature,and geopotential height fields:the 1.125°-resolution ERA-40 (Uppala et al.,2005);the 1.5°-resolution ERA-Interim(Dee et al.,2011);the 2.5°-resolution NCEP-2(Kanamitsu et al., 2002);and the 1.25°-resolution MERRA(Rienecker et al., 2011).ERA-Interim,MERRA and NCEP-2 cover the same period(1979–2013),while ERA-40 has data available for the period 1958–2002.The analysis will be conducted over the full available period.

    The number of vertical levels differs among the reanalyses,ranging from 17 levels in NCEP-2,23 levels in ERA-40, 37 levels in ERA-Interim,and42 levels in MERRA.Analysis will focus on comparing standard levels among the different reanalyses.Also,the 2.5°-resolution NOAA OLR data are used(Liebmann and Smith,1996).The reanalysis is available for the period 1974–2013.

    Precipitation estimates from two datasets are used:the 0.25°-resolution NASA TRMM satellite-derived monthly rainfall product for the period 1998–2013(TRMM 3B42V7; Kummerowetal.,1998),andversion1ofthe0.25°-resolution PERSIANN-CDR(Precipitation Estimation from Remotely Sensed Information using ArtificialNeural Networks-Climate Data Record)dataset(Sorooshian et al.,2000)for the period 1983–2013.The 0.5°-resolution CRU gridded monthly surface temperature dataset(CRUTS3.21;Mitchell and Jones,2005)is also used.This dataset provides landbased monthly rainfall and surface temperature estimates for the period 1901–2012.

    3.Results

    3.1.Identification of the Congo basin zonal overturning circulation

    Fig.1.ERA-Interim climatological monthly 900-hPa p-velocity(omega;×102Pa s-1)for the(a)JFM,(b)AMJ,(c) JAS and(d)OND means.Positive values indicate sinking motion.

    Figure 1 shows the ERA-Interim 1979–2013climatological vertical p-velocity[i.e.,omega,ω×10-2(Pa s-1)]at 900 hPa for the January–February–March(JFM,Fig.1a),April– May–June(AMJ,Fig.1b),July–August–September(JAS, Fig.1c),and October–November–December(OND,Fig.1d) means.In JFM,there is rising motion along the equatorial Atlantic(Fig.1a).In AMJ,this rising is replaced with sinking motion,up to 4×10-2Pa s-1,over the central equatorial Atlantic(7°–25°W)and the Gulf of Guinea near 8°E,clearly distinguishable from the central Atlantic maxima(Fig.1b). In JAS,both the central Atlantic and Gulf of Guinea subsidence maxima expand in spatial extent,with the latter extending from 4°E to 10°E between 5°S and 3°N(Fig.1c). The Guinean subsidence maximum reaches its peak strength of approximately 8×10-2Pa s-1,e.g.,near(1°S,8°E),and subsidence near 0°E extends northward up to 4°N.The maximum over the central equatorial Atlantic remains robust in size,but aloof from the Guinean maximum—for example,at (0°,8°E).This subsidence becomes sufficiently weak in winter,and is replaced with rising motion,up to-6×10-2Pa s-1(OND,Fig.1d).Over land,rising vertical motions occur year round.

    It is clear from Fig.1 that boreal summer is associated with subsidenceoverthe equatorialAtlantic between5°S and 3°N,with two centers of maximum intensity—one over the central Atlantic(Wang,2004)and the other over the Gulf of Guinea in the eastern Atlantic.Sinking motion maximizes over both regions in summer(Fig.1c).The Guinean subsidence exhibits a distinct seasonal cycle(not shown).The subsidence is primarily located to the south of the equator from January through May.In June,it appears over the Gulf of Guinea and starts its northward progression,maximizes and reaches its northernmost location of around 5°N from July to September,and then weakens and retreats south of the equator by October.ERA-40,NCEP-2 and MERRA exhibit a similar seasonal cycle(not shown).

    Figure 2 shows the longitude–height cross section of zonal and vertical p-velocity wind components[u(m s-1); ω×10-2(Pa s-1)]averaged between 5°S and 3°N from the ERA-Interimclimatology(1979–2013)for theJFM(Fig.2a), AMJ(Fig.2b),JAS(Fig.2c)and OND(Fig.2d)means. In JFM,there is rising motion,as indicated by the upwardpointing streamlines in the Gulf of Guinea(centered around 10°E;Fig.2a).Compared with the JFM mean,rising motion becomes weaker in AMJ(Fig.2b).In JAS(Fig.2c),there is a circulation near the surface between 2°W and 30°E,with rising motions from the surface to 650 hPa over the Congo basin from 12°E to 25°E,westward flow around 800 hPa from 15°E to 2°W,and subsidence along the equatorial West African coast at 2°W.Westerly onshoreflow from the Gulf of Guinea into the Congo basin is shallow,generally confinedto between the surface and 850 hPa,extending inland to about 12°E.The vertical extent of the circulation appears to be predominantlyconfinedtothe boundarylayer(i.e.,thesurfaceto 800 hPa),with evidenceof some mid-tropospheric(i.e.,700–400 hPa)rising(sinking)motion east(west)of 10°E.This circulation is identified as the Congo basin zonal overturning circulation.The circulation contracts in size and it almost disappears in OND(Fig.2d),although the rising branch remains strong over the Congo during this time(Fig.2d).This seasonality is also depicted in ERA-40,NCEP-2,MERRA, and JRA(1958–2014;not shown).

    Fig.2.ERA-Interim climatological monthly longitude–height cross section of streamlines[u(m s-1);-ω×102(Pa s-1)]averaged over 5°S–3°N for the(a)JFM,(b)AMJ,(c)JAS,and(d)OND means.

    A comparison of Figs.1 and 2 indicates that the seasonality of the Guinean subsidence east of 2°W and the Congo basin overturning circulation behave in a similar manner,especiallyduringtheborealsummermonthsofJunetoOctober, with the subsidence and Walker-type circulation’s strength being strongest from July to September.Thus,the rest of the analysis presented below will focus on physically understanding the maintenance of the Congo basin zonal overturning circulation and its relationship to subsidence over the Gulf of Guinea for the July–September period.Figures 3a–d show climatological July–September vertical cross sections of the zonal–omegawind components(vectors),and the meridionalwindcomponent(contours),averagedfrom5°S to 3°N forERA-Interim,ERA-40,NCEP-2andMERRA.All of the reanalyses show Walker-type circulation with rising motions over the Congo basin centered at 20°E and extending from near the surface to the upper troposphere,easterly flow around 650 hPa between 0°E and 20°E,and subsidence over the Gulf of Guinea that diverges at 2°W near the surface.Below 800 hPa from 2°W to 25°E,the flow is predominantly westerly,constitutingtheloweronshorebranchofthecirculation.Note that near the surface over this region there is also a southerly componentto the flow,but it is generally weaker in magnitude(e.g.,up to 3 m s-1in ERA-Interim).The circulation is embedded within a large-scale incomplete circulation that rises over the Congo basin,flows westward above 600 hPa,subsides over the Atlantic,but lacks the lower branch to complete the circulation.

    While the basic structure of the Congo basin zonal overturning circulation is captured by all of the reanalyses,there are still some differences among the datasets.For instance, the rising branch is stronger in ERA-Interim and ERA-40 than in NCEP-2 and MERRA.At 15°E and 700 hPa,the upward p-velocity is around 4×10-2Pa s-1in ERA-Interim, ERA-40 and NCEP-2,while it is replaced with a downward p-velocity of about-1×10-2Pa s-1in MERRA.

    Based upon Figs.2 and 3,the boundaries of the Congo basin zonal overturning circulation during July–September can be defined.The up-branch is located within(5°S–3°N, 15°–25°E)over the Congo basin,while the down-branch is within(5°S–3°N,2°W–8°E)over the Gulf of Guinea.Therefore,unless mentioned,the Gulf of Guinea and Congo basin regions in this paper are(5°S–3°N,2°W–8°E)and(5°S–3°N,15°–25°E),respectively.The 650 hPa easterly and the surface westerly constitute the up-branch and lower-branch (Figs.3a–d).

    The idea of there being a Walker-type circulation is not a new concept;previous studies have suggested the possibility of such a circulation(Vizy and Cook,2001;Lau and Yang, 2003;Hastenrath,2006).Pokam et al.(2014)reported the presence of this circulation.The lower branch is driven by the land–ocean heating contrast,and it develops well from September to November.As stated in the introduction,the focus here is on the West African monsoon season,July–September.

    Fig.3.Vertical cross section of zonal vertical p-velocity[u(m s-1);-ω×102(Pa s-1);vectors],vertical p-velocity [-ω×102(Pa s-1);shading]and meridional winds[v(m s-1);contours]averaged over 5°S–3°N for the JAS mean from(a)ERA-Interim,(b)ERA-40,(c)NCEP-2and(d)MERRA climatologies.Blue shading indicates upward motion. Values under the topography are masked white in(a),(c),and(d).

    The lower branch of the circulation is explored.Figures 4a–d display the July–September climatological 925-hPageopotentialheightsandwindsfromERA-Interim,ERA-40,NCEP-2 and MERRA,respectively.Heights are highover the equatorial Atlantic and low over the Congo basin, associated with the zonal height gradient between the Gulf of Guinea and the Congo basin(e.g.,heights approximately 20 gpm lower over the Congo basin than the Gulf of Guinea). Associated with this zonal gradient is predominantly westerly flow in the vicinity of the equator from about 10°E to 25°E,and the zonal flow is stronger than the meridional flow. For example,westerly wind reaches up to 6 m s-1,whereas southerly wind reaches up to 2 m s-1.Also,as shown in Fig. 3,the southerly flow is much weaker than the westerly flow in the region of interest from the surface to 800 hPa.These zonal flows remain undeflected by the Coriolis force.For these reasons,and also as discussed in the introduction,the focus here is on zonal circulation.

    To understand the maintenance of the zonal overturning circulation,the geopotential height gradient and its relationship with the low-level circulation,the horizontalmomentum equations are analyzed.The zonal(i.e.,u-)and meridional (i.e.,v-)components of the equation are as follows:

    where u is the zonal wind and v is the meridional wind.In Eq.(1)[Eq.(2)],the total zonal(meridional)acceleration is balanced on the right-hand side by the acceleration due to the zonal(meridional)pressure gradient-(?Φ/?x)[-(?Φ/?y)],where Φ is geopotential,the Coriolis acceleration component[for which f?2Ωsinφis the Coriolis parameter, where Ω is the angular speed of the rotation of Earth(= 7.3×10-5rad s-1)andφis latitude],and the zonal(meridional)frictional acceleration,represented by the terms Rx(Ry).Six-hourly values are used to calculate the daily and monthly values of the variables in the above equations.The acceleration term(Lagrangianacceleration)on the left is calculated as the sum of the Eulerian and convective accelerations.Note that Rxand Ryare calculated as residuals in the analysis,and therefore may contain errors due to the estimation of derivatives by finite differencing.

    Figure 5 shows the 925-hPaJuly–Septemberclimatological value of each term in the horizontal momentum equations [Eqs.(1)and(2)]from ERA-Interim,ERA-40,NCEP-2 and MERRA.Closetotheequator,thetotalacceleration(toprow) is small.Likewise,at the equator,the acceleration associated with the Coriolis force(third row)is small as f approaches zero.Thus,the primary balance occurs between the acceleration associated with the pressure gradient(second row)and the residual/friction term(fourth row).The pressure gradient force is primarily eastward over equatorial Africa,with the strongest magnitudes occurring between 12°E and 18°E in the Congo basin.In contrast,over the equatorial Atlantic, the pressuregradientaccelerationis muchweaker.Values are greater over coastal inland areas than over the ocean.For example,in ERA-Interim(Fig.5b),it is about 6×10-5m s-2over the equatorial Gulf of Guinea,while it is 18×10-5m s-2over the Congo basin.For the most part,the accelerationassociated with the residual/friction term opposes the wind direction(Fig.4),and has larger magnitude over land with complex topography,such as over western equatorial Africa. This shows that,associated with the friction due to topography,there is convergenceof winds over the Congo basin.

    Fig.5.Vectors showing combinations of terms in Eqs.(1)and(2),with contours displaying magnitudes(10-5m s-2)at 925 hPa for the JAS mean,as follows:total acceleration(first row);acceleration associated with the geopotential height gradient (second row);acceleration associated with the Coriolis force(third row);and residual terms(fourth row).

    In general,the 925-hPa geopotential heights are much lower in the Congo basin than in the Gulf of Guinea and,associated with this gradient,the flow is directed eastward and there is strongconvergenceoverthe Congo basin(Figs.4 and 5).Investigationof the underlyingsurface temperaturedistributions helps to understand the low-level geopotential height differences between the Congo basin and the Gulf.Figure 6 shows the climatological July–September surface temperatures for ERA-Interim(Fig.6a),ERA-40(Fig.6b),NCEP-2 (Fig.6c)and MERRA(Fig.6d),as well as the CRUTS.3.21 1901–2012 observational dataset(Fig.6e).Over the tropical eastern Atlantic,a cold tongue of relatively cooler SSTs forms,beginning in the late spring and reaching a minimum by late summer(Hastenrath and Lamb,2004).This feature is represented in all of the reanalyses shown in Fig.6 by the area of relatively cooler temperatures of 296 K between 5°S and the equator spanning from 10°W to the African coast. North of the equator,Atlantic SSTs are slightly warmer(i.e., up to 298–299 K).Generally,all of the reanalyses exhibit similar SST distributions over the Atlantic,which is not surprising since these datasets are typically forced with the observed SSTs.Over continental Africa,surface temperatures vary amongthe datasets,attributable to the differentland surface models included in predicting them.For example,compared with the observed CRUTS3.21 surface temperatures, the Congo basin temperatures are about 2 K cooler in ERAInterim,ERA-40 and NCEP-2,while around 4 K warmer in MERRA(Fig.6e).Despite these differences,the reanalyses generally indicate the existence of a surface temperature difference near the equatorbetween the relativelycooler Gulf of Guinea and the relatively warmer Congo basin.

    Fig.6.Climatological JAS surface temperatures(K)from(a)ERA-Interim,(b)ERA-40,(c)NCEP-2 and(d) MERRA,as well as from(e)CRUTS3.21 observations.Contour interval:1 K.

    Next,climatological monthly 925-hPa temperatures are area-averaged over the Gulf of Guinea(5°S–3°N,2°W–8°E) and the Congo basin(5°S–3°N,15°–25°E)and compared to one another for each reanalysis.Figure 7a shows the monthly Congo basin(long-dashed line)and the Gulf of Guinea(short-dashed line)surface temperatures,and the differences in the Congo basin and the Gulf of Guinea temperatures are shown in Fig.7b.Note that it is chosen to conduct this evaluation at 925 hPa to reduce the impact that differences in surface elevation may have on the comparison.

    The results in Fig.7a indicate that all of the reanalyses generally agree that the Congo basin 925-hPa temperatures are warmer than those of the Gulf of Guinea for all months. Furthermore,the seasonal variability over the annual cycle is much less over the Congo basin compared to the Gulf of Guinea,as temperatures generally fluctuate by 2 K or less for the former compared to around 5–6 K for the latter.

    Fig.7.(a)925-hPa climatological monthly mean temperatures(K)in the Congo basin(5°S–3°N,15°–25°E;long-dashed lines)and the Gulf of Guinea (5°S–3°N,2°W–8°E;short-dashed lines)from ERA-Interim(black),ERA-40(green),NCEP-2(red)and MERRA(blue).(b)925-hPa climatological monthly mean temperature differences(K)between the Congo basin and the Gulf of Guinea from ERA-Interim(black),ERA-40(green),NCEP-2(red)and MERRA(blue).

    In terms of the Congo basin and Gulf of Guinea temperature difference,it is positive year round,but relatively small (i.e.,3 K or less)from November to April for three of the fourreanalyses(Fig.7b).Theexception,MERRA,hasdifferences between 4 and 5.5 K during November–April.During the boreal summer months,the differences increase in May and June,peaking in July and August(at around5 K in ERAInterim and NCEP-2,7.3 K in ERA-40,and above 8 K in MERRA),and then after August the difference begins to decreaseasthecoldupwellingweakensintheGulfofGuinea.It is this low-level temperature difference that drives the overturning circulation.The circulation appears when the temperature difference starts increasing in June,becomes well developed from July to September when the temperature difference is largest,and,as the eastern equatorial Atlantic cold tongue weakens,the circulation starts weakening and disappears in the fall(Figs.2 and 7).

    These results suggest that it is important to have an accurate representation of the Congo basin/Gulf of Guinea lowlevel temperature gradient associated with the zonal circulation.Guinean temperatures are similar in the reanalyses; therefore,gradient differences are largely determined by the Congo basin temperatures.Many potential factors can influence the low-level temperatures over land,which can result in the inter-reanalysis spread shown in Fig.7.One of the important factors is how each reanalysis treats the prediction of land surface conditions;however,understandingthis is better left to the individual reanalysis groups,as they will have a better understanding of the intricacies of their respective reanalysis modeling algorithms.A more basic factor is to evaluate the precipitation in the different reanalyses,since temperature and precipitation are tightly related.For this reason, precipitation in the reanalyses are examined next.

    Climatological July–September precipitation rates are shown in Fig.8 for the TRMM(Fig.8a)and PERSIANN (Fig.8b)observations,as well as for ERA-Interim(Fig.8c), ERA-40(Fig.8d),NCEP-2(Fig.8e)and MERRA(Fig.8f). TRMM and PERSIANN indicate a zonally oriented band of precipitation across Africa during July–September,with maxima centered at(2°N,18°W)and(5°N,8°E)over the Cameroon highlands(Figs.8a and b).Rainfall rates over the northern Democratic Republic of Congo(DRC)are over 7 mm d-1and gradually decrease to about 3 mm d-1around 5°S.

    The zonally elongated band of rainfall with the two maxima is represented in each reanalysis(Figs.8c–f),albeit the magnitudes of the maxima vary among the differentresolutiondatasets.OvertheDRC,therearedifferencesinthe rainfall distributions among the reanalyses,with MERRA, ERA-40 and NCEP-2 demonstrating considerably drier results than TRMM.Of the four reanalyses,ERA-Interim captures the July–September distribution of rainfall over the DRC most realistically.

    Next,the 600-hPa winds and temperature are examined (figures not shown).The flow is primarily easterly from the Congo basin to the Gulf of Guinea near the equator,and this easterly flow constitutes the returning branch of the overturning circulation.Similar to the lower branch,this branch is driven by the temperature gradient.The Gulf of Guinea is warmer than the Congo basin.The difference is positive, up to 0.06 K in ERA-Interim,0.27 K in ERA-40,0.3 K inNCEP-2,and 0.35 K in MERRA.Exploring the relationship between the upper branch of this equatorial overturning circulationalong the equatorand the African easterly jet located further north in the Sahel is beyond the scope of this study.

    3.2.Relationship between the Congo basin zonal overturning circulation and West African atmospheric circulation

    There are differences in the representation of ocean–land temperature gradient in the reanalyses.For example,during the summer,temperature differences vary from 5 to 8 K among the reanalyses.These gradient differences are associated with variations in representing land temperature,as discussed above(Fig.7).Since land temperature is closely tied with precipitation,the precipitation values in the reanalyses are investigated.Compared to the other three reanalyses (MERRA,ERA-40 and NCEP-2),ERA-Interim accurately represents the precipitation values(Fig.8).This is consistent with the findings of Simmons et al.(2010),who also reported that the precipitation field is generally better represented in ERA-Interim compared with ERA-40.Also,in representing the Congo basin precipitation,large differencesexist among datasets and models,as shown by Washington et al.(2013).In total,including the spring,which was the focus of their analysis,ERA-Interim and ERA-40 yield wetter results than the NCEP reanalysis.Lower summertime rainfall rates(for example in MERRA)are likely associated with increased surface heating due to a lack of convection,warmerthan-observed surface temperatures predicted over the DRC(e.g.,F(xiàn)ig.6b),andhencealargerlow-levelCongobasin-Gulf of Guinea temperaturecontrast(Fig.7b).Since ERA-Interim provides the most realistic representation of climatological July-September surface temperature and precipitation fields,as reported in section 3.1(Figs.6 and 8),this reanalysis is used to examine the relationship.I hypothesize that,when the rising branch of the circulation is strong,the northward moisture transport and rainfall becomes strong across West Africa,including the Sahel.

    Monthly omega wind anomalies for July-September from 1979 to 2013 are calculated by subtracting the 1979-2013 climatological mean.Doing so removes the seasonal cycle from the time series to foster a better comparison over different months.From the above anomalies field,an 800-hPa area-average index is created.The area used is(5°S-3°N,15°-25°E).This area corresponds to the rising branch of the overturningcirculation,discussed earlier in section 3.1(Fig.2),while 800 hPa is chosen to be sufficiently above the topography of the region.The black box shown in Fig.9 denotes the averaging region.

    Figure 9 shows the correlation coefficients(multiplied by 100)between the Congo basin index and the July-September p-velocity field at 850 hPa for ERA-Interim.The 850-hPa level is selected because the Guinean subsidence is much larger at this level(see Fig.3a).Correlations are negative,up to 50%(significant at the 5%level),in the equatorial Atlantic(5°S-3°N,25°W-7°E),including the Gulf of Guinea. This implies that,when the rising motions over the Congo basin become stronger,subsidence over the Gulf of Guinea becomes strong,and vice versa.North of the equator,significant positive correlations are centered along 10°N over West Africa and southern Chad.This suggests that stronger rising motion over the Congo basin is also associated with stronger rising motion in the Sahel.The Guinean index(5°S-3°N,2°W-8°E)shows much stronger positive correlations(up to 100%)in the Gulf and negative correlations in the Congo basin,lending further support to the existence of a relationship between the Guinean subsidence and the Congo basin rising,as discussed above(figure not shown).

    Figure 10a shows the time series of the 800-hPa vertical p-velocityindexforJAS onlyfromERA-Interim.We investigate the physics of this variability.A negativevalue indicates a stronger rising motion.The index is mainly negative from 1979to1994.In1995,thereis anabruptshiftandtheindexis found to be primarily positive after 1994.For example,from 1979 to 1994,37 of the total 48 JAS periods are negative;while from 1995 to 2013,49 of the total 57 are positive.This discontinuity is likely associated with various adjustments madeto theERA-Interimalgorithm(e.g.,newbias correction schemes,new moisture analysis and model physics schemes havebeenusedtomaintainconsistencyamongobservingsystems)and data input(Simmons et al.,2007;Uppala et al.,2008;Dee and Uppala,2009;Kobayashi et al.,2009).

    Fig.9.Correlation coefficients(multiplied by 100)between the ERA-Interim 800-hPa Congo basin p-velocity index and the vertical p-velocity for JAS.Correlation coefficients exceeding the 95%confidence level are enclosed within red and greendashed lines.Region(5°S-3°N,15°-25°E)indicates the Congo basin.

    To remove this shift,the time series is divided into two parts(1979-1994 and 1995-2013),and each part of the time series is scaled with its own climatology.For example,the 1979-1994 July climatology is subtracted from the 16 months of July during 1979-1994 to obtain the deseasonalised value for July,with August and September created in a similar manner.Likewise,the 1995-2013 JAS periods are similarly produced,but by utilizing the 1995-2013 climatology for the appropriate month.

    Figure 10b shows the adjusted Congo basin index.The dashed lines correspond to±1 standard deviation.Positive and negative index values are now more evenly distributed over the time series and there is no longer any evidence of a discontinuity in 1995.The index shown in Fig.10b is used to detect individual months when the Congo basin vertical p-velocity index is strong and weak for JAS,and hence identify particular months when the circulation is strong or weak.A strong circulation month is defined as when the index value is less than-0.6×10-2Pa s-1,which corresponds to-1 standard deviation from the mean;whereas,a weak circulation month is defined as when the index is greater than +0.6×10-2Pa s-1,which corresponds to+1 standard deviation from the mean.Table 1 shows the months identified for each case.This index represents the zonal overturning circulation.

    Fig.10.(a)800-hPa average(5°S–3°N,15°–25°E)vertical p-velocity(×102Pa s-1)differences in JAS and their climatologies from ERA-Interim(1979–2013).(b)As in(a)but the vertical p-velocity values before 1994 are subtracted from the 1979–94 mean,and the values after 1995 are subtracted from the 1995–2013 mean.The dashed lines in(b)indicate±1 standard deviation.

    Figures 11a-c show the JAS vertical cross sections of the zonal and vertical components of the wind(streamlines)and vertical winds(shaded)averaged from 5°S to 3°N for the reanalysis climatology(Fig.11a),as well as the weak(Fig. 11b)and strong(Fig.11c)circulation composites.In the climatology,the circulation extends from the surface to 700 hPa,where the upward p-velocity is about 5×10-2Pa s-1(Fig.11a).Intheweakcomposite,therising branchweakens; for example,the rising motion at(15°E,800 hPa)is about 2×10-2Pa s-1weaker than in the climatology(Fig.11b). In contrast,upward vertical motions are stronger than in the climatology,by up to 2×10-2Pa s-1,in the strong circulation composite case(Fig.11c).The sinking in the Gulf and rising in the Congo basin in the strong composite is significantly different(at the 95%confidence level)from the weak composite(notshown),meaningthatastrongcirculationcase is different from a weak circulation case.

    The relationship between the overturning equatorial circulation and regional-scale circulation and rainfall in West Africa is explored by analyzing the low-level circulation and moisture fields in Figs.11d–f.Figure 11d shows the climatological800-hPageopotentialheightsandmoisturetransport vectorsfortheJASclimatology.The800-hPalevelisselected for consistency with the level used to create the two composite cases,as discussed earlier.Associated with the North Atlantic subtropical high,the largest heights are located off the coast of West Africa.The smallest height is located over the continent,east of 15°E.Moisture transport is predominantly westerly–southwesterly over the Guinean coast and coastal inland areas(0°–7°N,0°–15°E).Transport vectors are primarily easterly to the south of the equator.

    Fig.11.ERA-Interim longitude–height cross section of streamlines[u(m s-1);-ω×102(Pa s-1)]and vertical winds (×102Pa s-1;shaded),averaged over 5°S–3°N from the(a)JAS climatology,and the(b)weak and(c)strong Congo basin Walker circulation composites.Blue coloring indicates upward vertical motion.(d)ERA-Interim JAS climatological 800-hPa moisture transport(q. uuu;vectors)and geopotential heights(shaded).(e)Moisture transport and geopotential height differences in the weak composite and the climatology.(f)As in(e)but for the strong composite.Units for geopotential height and moisture transport are gpm and kg kg-1m s-1.

    Figure 11e displays the 800-hPa geopotential heights and moisture transport differences in the weak composite and the climatology,and the differences in the strong composite and climatology are shown in Fig.11f.Geopotential heights are larger,up to 2 gpm,in the weak composite than in the climatology in the Sahel(10°–20°N,25°W–20°E),and the heights are up to-1.5 gpm lower over the Gulf and the Congo basin (Fig.11e).Associated with these height differences,easterly and northeasterly moisture transport anomalies develop over the Congo basin and the Sahel.Also,with weakeningof the west African westerly jet(Pu and Cook,2010),westerly moisture transport from the West African coast to the central Sahel lowers,as shown by the easterly moisture transport anomalies around11°N,and the moisture level reduces in the central Sahel(not shown).In contrast,there is a general reduction,up to-2.5 gpm,of height over the Sahel and the Congo basin(?5°S,18°E)in the strong composite.This is associated with increases in moisture transport,as indicated by the westerly and southwesterly anomalies over the western coast of West Africa and the Gulf(Fig.11f).As the zonal circulationbecomesstronger,surfaceto 800-hPamoisture increases in the Sahel and coastal inland areas within(2°–7°N,0°–20°E)(not shown).As discussed in L′el′e et al.(2015), the low-level flows—especially below 850 hPa—constitute an important source of moisture supplying the West African monsoon system.This means that when the rising branch of the circulation is anomalously strong,the Congo basin and Sahelian moisture transports become strong.

    The JAS cloud cover in ERA-Interim matches fairly well with ISCCP(1983–2009)observations(not shown).Most of the cloud in ERA-Interim over the analysis region is low cloud(?70%,below 800 hPa).Changes in the low-cloud cover are examined.The JAS mean low-cloud cover differences displayed in Fig.12a are for the weak composite and climatology,while Fig.12b shows the results for the strong composite and climatology.The cloud over land,including over the Sahel and the Congo basin,and along the coastal upwelling region to the south of the equator,reduces,up to -6%,in the weak composite(Fig.12a).Meanwhile,the cloud cover increases,up to 4%,over the northern Guinean Gulf,resulting in lower amounts of solar flux to the surface. In contrast to the weak composite,the cloud amounts in the strong composite increase generally over land but decrease over the ocean(Fig.12b).An enhancement,up to 12%,occurs over the Congo basin and Sahel,permitting relatively smaller amounts of solar radiation.Over the Gulf of Guinea, cloud amounts lower,up to-8%,allowing more incoming solar radiation to heat up the surface(not shown).

    Figure 12c displays the differences in OLR in the weak compositeandclimatology(fromtheNOAA OLR),whilethe differences in the strong composite and climatology are displayed in Fig.12d.Changes in cloud cover in ERA-Interim are consistent with the NOAA OLR.Associated with the decreased cloud cover over the Sahel and the Congo basin,the OLR valueis,upto 10W m-2,greaterin the weakcomposite than in the climatology(Figs.12a and c).Also,associated with the increased cloud cover,reductions,up to-6 W m-2, occur in the strong composite(Figs.12b and d).

    The difference between the weak and strong composites is largely due to land surface heating.Associated with weak rising motion over the Congo basin,cloud cover reduces, allowing larger amounts of incoming solar radiation to increase the surface temperature in the weak composite.In contrast,cloud amounts increase generally over land,reducing the amount of incoming solar heat flux,and this cools the surface in the strong composite.

    Figure 12e displays the PERSIANN-CDR JAS precipitation differences in the weak composite and climatology, whilethe differencesin thestrongcompositeandclimatology are displayed in Fig.12f.Associated with the weakening of the Congo basin circulation,rainfall reduces,by up to-2.5 mm d-1,in the central Sahel,southern Chad.In contrast, as the circulation becomes strong,rainfall rates enhance,by up to 2.5 mm d-1in the central Sahel and the Congo basin. Also,rainfall rates increase,byup to 2mm d-1,overthewest coast of Africa in the strong composite,while reductions occur in the weak composite.Investigating the causes of these differences is beyond the scope of the present study.

    In summary,the Congo basin zonal overturning circulation is tightly connected with the West African monsoon system.A strong circulation is associated with increased subsidence over the Gulf of Guinea,and this supports greater northward moisture transport and increases of precipitation in the Sahel(Figs.11 and 12f).

    4.Conclusions

    This paper comprehensively documents the Gulf of Guinea subsidence and identifies a Congo basin zonal overturningcirculation,whosedownbranchis associated withthe presenceofthis subsidence.Theinterannualvariabilityofthe circulation and its relationship with moisture and precipitation distributions across West Africa are explored.

    The Guinean subsidence develops in the low levels at the beginning of June,attains a maximum intensity of 6×10-2Pa s-1,reaching as far as 5°N during boreal summer(JAS), and dissipates in October.This subsidence is associated with a complete zonal circulation,which has the same seasonality as that of the Guinean subsidence in ERA-Interim,ERA-40,NCEP-2 and MERRA.The circulation is identified as a Congo basin zonal overturning circulation.The circulation comprises of:(1)a low-level westerly flow from the Gulf to the Congo basin(lower branch);(2)rising winds over the basin around 20°E(up-branch);(3)easterly return flow from the basin to the Gulf of Guinea at 600 hPa(upper branch); and(4)subsiding winds over the Gulf centered around 2°W (down-branch)—with winds to the east flowing eastward to complete the circulation.This circulation reaches its maximum intensity in the summer from July to September.

    Analysis of the momentum budget shows that the lowlevel flow is strongly ageostrophic,with a small acceleration term and Coriolis term,which is zero at the equator, and large pressure gradient and friction terms.The eastwarddirected pressure gradient force drives the flow from the Gulf of Guinea toward the Congo basin,constituting the lower branchofthezonaloverturningcirculation.Thepressuregradient term is larger over the Congo basin than over the Gulf of Guinea.Over the basin,it is comparable to the residual term,which opposes the wind direction.The residual is interpreted as friction,and is expected to be much larger over land with complex topography.Associated with this friction, winds converge and rise over the Congo basin forming the up-branch of the circulation.Causes of the pressure gradient force are explored.

    Analysis ofsurface temperatureexplainsthe eastward decreasing pressure gradient force that drives the circulation and much of its seasonal variation.Temperature is(in ERAInterim,for example)up to-5 K lower in the Gulf of Guinea than in the Congo basin,and,associated with this gradient, the surface pressure is lower over land than over the ocean.

    The differences in the Gulf of Guinea–Congo basin temperatures are primarily caused by differences in the land temperatures,which are model-dependent,in the reanalyses. Land temperature and precipitation in ERA-40,NCEP-2 and MERRA are muchwarmer and drierthan in the observations. ERA-Interim is selected to understand the variability of thecirculation based on its realistic representation of temperature and precipitation.

    Fig.12.(a)Low-cloud cover(%)differences in the weak composite and the JAS mean from the ERA-Interim climatology.(b)As in(a)but for the strong composite.(c)OLR(W m-2)differences in the weak composite and the JAS mean from the NOAA climatology.(d)As in(c)but for the strong composite.(e)Precipitation(mm d-1)differences in the weak composite and the JAS mean from the PERSIANN-CDR climatology.(f)As in(e)but for the strong composite. Contour interval is every 2%in(a,b)and 2 W m-2in(c,d).

    Correlation between the Congo basin p-velocity index and the low-level vertical velocity field on interannual timescales reveals that the Congo basin rising motion is significantly correlated,up to-0.5,with the Guinean subsidence in ERA-Interim.Based upon the strength of the upbranch,two composites are created to explore the relationship between the equatorial overturning circulation and regional circulation and rainfall over the Sahel.Examination of the low-level flows shows that the Congo basin circulation is related to moisture distributions across West Africa and the Congo basin.Along with the weakening of the circulation,the Guinean subsidence becomes weak.This supports the development of convection over the Gulf of Guinea,and northward moisture transport and precipitation over the Sa-hel reduces.In contrast,when the circulation is stronger,the Guinean subsidence becomes stronger,and this inhibits convection over the Gulf of Guinea,supporting northward and eastward moisture transport associated with enhancementsof precipitation over the Sahel and the Congo basin.

    Having established the important role of the Congo basin zonal overturning circulation in moisture and precipitation distributions across West Africa and the Congo basin in observations,furtherinvestigationof the circulationin the stateof-the-art climate models will be addressed in a subsequent paper.

    Acknowledgements.This project is supported by the Jackson School of Geosciences/University of Texas at Austin.The ERAInterim and ERA-40 data are available at http://apps.ecmwf.int/ datasets/data/interim-full-daily/and http://apps.ecmwf.int/datasets/ data/era40-daily/,respectively;the NCEP-2 data are available from the ESRLwebsite,http://www.esrl.noaa.gov/psd/data/gridded/ data.ncep.reanalysis2.html;the MERRA data are available at http:// disc.sci.gsfc.nasa.gov/daac-bin/FTPSubset.pl?LOOKUPID List= MAIMCPASM;the NOAA OLR data are available at http:// www.esrl.noaa.gov/psd/data/gridded/data.interp OLR.html;and the PERSIANN-CDR data can be downloaded from ftp://data.ncdc. noaa.gov/cdr/persiann/files.Thanks to the two independent reviewers for their helpful comments.

    REFERENCES

    Caniaux,G.,H.Giordani,J.-L.Redelsperger,F.Guichard,E. Key,and M.Wade,2011:Coupling between the Atlantic cold tongue and the West African monsoon in boreal spring and summer.J.Geophys.Res.,116,C04003,doi:10.1029/2010 JC006570.

    Cook,K.H.,and E.K.Vizy,2006:Coupled model simulations of the West African monsoon system:Twentieth-and twentyfirst-century simulations.J.Climate,19,3681–3703.

    Cook,K.H.,and E.K.Vizy,2015:The Congo Basin walker circulation:Dynamics and connections to precipitation.Climate Dyn.,1–21,doi:10.1007/s00382-015-2864-y.

    Dee,D.P.,and S.Uppala,2009:Variational bias correction of satellite radiance data in the ERA-Interim reanalysis.Quart. J.Roy.Meteor.Soc.,135,1830–1841.

    Dee,D.P.,and Coauthors,2011:The ERA-Interim reanalysis: Configuration and performance of the data assimilation system.Quart.J.Roy.Meteor.Soc.,137,553–597.

    Dezfuli,A.K.,and S.E.Nicholson,2013:The relationship of rainfall variability in western equatorial Africa to the tropical Oceans and atmospheric circulation.Part II:The boreal autumn.J.Climate,26(1),66–84,doi:10.1175/JCLI-D-11-00686.1.

    Dezfuli,A.K.,B.F.Zaitchik,and A.Gnanadesikan,2015: Regional Atmospheric circulation and rainfall variability in south equatorial Africa.J.Climate,28(2),809–818,doi: 10.1175/JCLI-D-14-00333.1.

    Fontaine,B.,P.Roucou,and S.Trzaska,2003:Atmospheric water cycle and moisture fluxes in the West African monsoon: mean annual cycles and relationship using NCEP/NCAR reanalysis.Geophys.Res.Lett.,30,1117,doi:10.1029/2002GL 015834.

    Grist,J.P.,and S.E.Nicholson,2001:A study of the dynamic factors influencing the rainfall variability in the west African Sahel.J.Climate,14,1337–1359.

    Hagos,S.M.,and K.H.Cook,2007:Dynamics of the West African monsoon jump.J.Climate,20,5264–5284.

    Hagos,S.M.,and C.D.Zhang,2010:Diabatic heating,divergent circulation and moisture transport in the African monsoon system.Quart.J.Roy.Meteor.Soc.,136,411–425.

    Hastenrath,S.,2001:In search of zonal circulations in the equatorial Atlantic sector from the NCEP-NCAR reanalysis.Int.J. Climatol.,21,37–47.

    Hastenrath,S.,2006:Circulation and teleconnection mechanisms of northeast Brazil droughts.Prog.Oceanogr.,70,407–415.

    Hastenrath,S.,and P.J.Lamb,2004:Climate dynamics of atmosphere and ocean in the equatorial zone:A synthesis.Int.J Climatol.,24,1601–1612.

    Hastenrath,S.,and D.Polzin,2011:Long-term variations of circulation in the tropical Atlantic sector and Sahel rainfall.Int. J.Climatol,31,649–655.

    Kanamitsu M.,W.Ebisuzaki,J.Woollen,S.K.Yang,J.J.Hnilo, M.Fiorino,and G.L.Potter,2002:NCEP-DOE AMIP-II Reanalysis(R-2).Bull.Amer.Meteor.Soc.,83,1631–1643.

    Kobayashi,S.,M.Matricardi,D.Dee,and S.Uppala,2009:Toward a consistent reanalysis of the upper stratosphere based on radiance measurements from SSU and AMSU-A.Quart. J.Roy.Meteor.Soc.,135,2086–2099.

    Kummerow,C.,W.Barnes,T.Kozu,J.Shiue,and J.Simpson, 1998:The Tropical Rainfall Measuring Mission(TRMM) sensor package.J.Atmos.Oceanic Technol.,15,809–817.

    Lau,K.-M.,and S.Yang,2003:Walker circulation.Encyclopedia of Atmospheric Sciences,J.R.Holton et al.,Eds.,Academic Press,2505–2510.

    Leduc-Leballeur,M.,G.de Co¨etlogon,and L.Eymard,2013:Air–Sea interaction in the Gulf of Guinea at intraseasonal timescales:Wind bursts and coastal precipitation in boreal spring. Quart.J.Roy.Meteor.Soc.,139,387–400,doi:10.1002/qj. 1981.

    L′el′e,M.I.,L.M.Leslie,and P.J.Lamb,2015:Analysis of lowlevel atmospheric moisture transport associated with the West African monsoon.J.Climate,28,4414–4430,doi:10.1175/ JCLI-D-14-00746.1.

    Liebmann,B.,and C.A.Smith,1996:Description of a complete(interpolated)outgoing longwave radiation dataset.Bull. Amer.Meteor.Soc.,77,1725–1277.

    Mitchell,T.D.,and P.D.Jones,2005:An improved method of constructing a database of monthly climate observations and associated high-resolution grids.Int.J.Climatol.,25,693–712,doi:10.1002/joc.1181.

    Neupane,N.,and K.H.Cook,2013:A nonlinear response of Sahel rainfall to Atlantic warming.J.Climate,26,7080–7096, doi:10.1175/JCLI-D-12-00475.1.

    Nicholson,S.E.,and P.J.Webster,2007:A physical basis for the interannual variability of rainfall in the Sahel.Quart.J.Roy. Meteor.Soc.,133,2065–2084.

    Nicholson,S.E.,and A.K.Dezfuli,2013:The relationship of rainfall variability in western equatorial Africa to the tropical Oceans and atmospheric circulation.Part I:The boreal spring. J.Climate,26(1),45–65,doi:10.1175/JCLI-D-11-00653.1.

    Nolan,D.S.,C.D.Zhang,and S.H.Chen,2007:Dynamics of the shallow meridional circulation around intertropical convergence zones.J.Atmos.Sci.,64,2262–2285.

    Pokam,W.M.,L.A.T.Djiotang,and F.K.Mkankam,2012:At-mospheric water vapor transport and recycling in equatorial central Africa through NCEP/NCAR reanalysis data.Climate Dyn.,38(9–10),1715–1729,doi:10.1007/s00382-011-1242-7.

    Pokam,W.M.,C.L.Bain,R.S.Chadwick,R.Graham,D.J. Sonwa,and F.M.Kamga,2014:Identification of processes driving low-level westerlies in west equatorial Africa.J.Climate,27(11),4245–4262,doi:10.1175/JCLI-D-13-00490.1.

    Pu,B.,and K.H.Cook,2010:Dynamics of the West African westerly jet.J.Climate,23(23),6263–6276.

    Rienecker,M.M.,and Coauthors,2011:MERRA:NASA’s modern-era retrospective analysis for research and applications.J.Climate,24,3624–3648.

    Segele,Z.T.,P.J.Lamb,and L.M.Leslie,2009:Large-scale atmospheric circulation and global sea surface temperature associations with Horn of Africa June–September rainfall.Int. J.Climatol.,29(8),1075–1100.

    Simmons,A.,S.Uppala,D.Dee,and S.Kobayashi,2007: ERA-Interim:New ECMWF Reanalysis Products from 1989 Onwards.ECMWF Newsletter,No.110,ECMWF,Reading, United Kingdom,25–35.

    Simmons,A.J.,K.M.Willett,P.D.Jones,P.W.Thorne,and D.P.Dee,2010:Low-frequency variations in surface atmospheric humidity,temperature,and precipitation:Inferences from reanalyses and monthly gridded observational data sets. Journal of Geophysical Research:Atmospheres(1984–2012), 115,D01110,doi:10.1029/2009JD012442.

    Sorooshian,S.,K.-L.Hsu,X.G.Gao,H.V.Gupta,B.Imam, and D.Braithwaite,2000:Evaluation of PERSIANN system satellite-based estimates of tropical rainfall.Bull.Amer.Meteor.Soc.,81,2035–2046.

    Thorncroft,C.D.,and M.Blackburn,1999:Maintenance of the African easterly jet.Quart.J.Roy.Meteor.Soc.,125,763–786.

    Thorncroft,C.D.,H.Nguyen,C.D.Zhang,and P.Peyrill′e,2011: Annual cycle of the West African monsoon:Regional circulations and associated water vapour transport.Quart.J.Roy. Meteor.Soc.,137,129–147.

    Trenberth,K.E.,D.P.Stepaniak,and J.M.Caron,2000:The global monsoon as seen through the divergent atmospheric circulation.J.Climate,13,3969–3993.

    Uppala,S.M.,and Coauthors,2005:The ERA-40 re-analysis. Quart.J.Roy.Meteor.Soc.,131,2961–3012.

    Uppala,S.M.,D.P.Dee,S.Kobayashi,P.Berrisford,and A.J. Simmons,2008:Towards a climate data assimilation system: status update of ERA-Interim.ECMWF Newsletter,No.115, ECMWF,Reading,United Kingdom,12–18.

    Vizy,E.K.,and K.H.Cook,2001:Mechanisms by which Gulf of Guinea and eastern North Atlantic Sea surface temperature anomalies can influence African rainfall.J.Climate,14,795–821.

    Wang,C.,2005:ENSO,Atlantic climate variability,and the Walker and Hadley circulations.The Hadley Circulation: Present,Past,and Future,H.F.Diaz and R.S.Bradley,Eds., Kluwer Academic Publishers,173–202.

    Wang,C.Z.,2002a:Atlantic climate variability and its associated atmospheric circulation cells.J.Climate,15,1516–1536.

    Wang,C.Z.,2002b:Atmospheric circulation cells associated with the El Ni?o–southern oscillation.J.Climate,15,399–419.

    Wang,C.Z.,2004:ENSO,Atlantic climate variability,and the Walker and Hadley circulations.The Hadley Circulation: Present,Past,and Future,H.F.Diaz and R.S.Bradley, Eds.,Advances inGlobal Change Research,Vol.21,Springer, Netherlands,85–120.

    Washington,R.,R.James,H.Pearce,W.M.Pokam,and W. Moufouma-Okia,2013:Congo Basin rainfall climatology: can we believe the climate models?Philosophical Transactions of the Royal Society of London.Series B,Biological Sciences,368(1625),20120296,doi:10.1098/rstb.2012.0296.

    Zhang,C.D.,P.Woodworth,and G.J.Gu,2006:The seasonal cycle in the lower troposphere over West Africa from sounding observations.Quart.J.Roy.Meteor.Soc.,132,2559–2582.

    Neupane,N.,2016:The Congo basin zonal overturning circulation.Adv.Atmos.Sci.,33(6),767–782,

    10.1007/ s00376-015-5190-8.

    Naresh NEUPANE

    Email:nareshneupane@utexas.edu

    欧美日韩在线观看h| 女性被躁到高潮视频| 国产精品国产三级专区第一集| 天堂中文最新版在线下载| 亚洲真实伦在线观看| 人妻夜夜爽99麻豆av| 亚洲国产欧美在线一区| 欧美xxxx性猛交bbbb| 精品国产一区二区三区久久久樱花| 啦啦啦视频在线资源免费观看| 亚洲不卡免费看| 看非洲黑人一级黄片| 黄色视频在线播放观看不卡| 欧美日韩亚洲高清精品| 大又大粗又爽又黄少妇毛片口| 午夜福利,免费看| 欧美xxⅹ黑人| 一区二区三区精品91| 人人澡人人妻人| 观看美女的网站| 国产在线男女| 欧美日韩国产mv在线观看视频| 亚洲精品视频女| 亚洲欧美中文字幕日韩二区| 99国产精品免费福利视频| 免费观看无遮挡的男女| 免费av不卡在线播放| 精品午夜福利在线看| 性色av一级| 赤兔流量卡办理| 国产精品福利在线免费观看| 国产亚洲午夜精品一区二区久久| 亚洲av电影在线观看一区二区三区| 高清不卡的av网站| 亚洲人与动物交配视频| 在线观看人妻少妇| 国产在线视频一区二区| 中国三级夫妇交换| 色婷婷久久久亚洲欧美| 国产片特级美女逼逼视频| 你懂的网址亚洲精品在线观看| 久热久热在线精品观看| 又大又黄又爽视频免费| 男女无遮挡免费网站观看| 男女边吃奶边做爰视频| 五月开心婷婷网| 丝袜脚勾引网站| 色婷婷av一区二区三区视频| 99热这里只有是精品50| 一区二区三区免费毛片| 男女边摸边吃奶| 91久久精品国产一区二区三区| 丝袜在线中文字幕| 国产探花极品一区二区| 自拍欧美九色日韩亚洲蝌蚪91 | 成人午夜精彩视频在线观看| 99热网站在线观看| a级片在线免费高清观看视频| 亚洲一区二区三区欧美精品| 国产免费又黄又爽又色| 国产日韩一区二区三区精品不卡 | .国产精品久久| 欧美精品一区二区免费开放| 成人18禁高潮啪啪吃奶动态图 | 啦啦啦啦在线视频资源| 高清视频免费观看一区二区| 久久99蜜桃精品久久| 国产成人精品婷婷| 最近手机中文字幕大全| 久久久久久久久久久丰满| 伦精品一区二区三区| a级毛片在线看网站| 亚州av有码| 嫩草影院入口| 伊人久久国产一区二区| 欧美区成人在线视频| 国产欧美日韩综合在线一区二区 | 国产色爽女视频免费观看| 五月开心婷婷网| 少妇人妻一区二区三区视频| 黄色怎么调成土黄色| 少妇人妻精品综合一区二区| 久久毛片免费看一区二区三区| 成年美女黄网站色视频大全免费 | 一二三四中文在线观看免费高清| 超碰97精品在线观看| 两个人免费观看高清视频 | 午夜福利视频精品| 中文字幕免费在线视频6| 国产伦在线观看视频一区| 国产黄片美女视频| 日韩欧美 国产精品| 十分钟在线观看高清视频www | 亚洲性久久影院| 男人舔奶头视频| 黑人猛操日本美女一级片| 女性被躁到高潮视频| 中文字幕久久专区| 这个男人来自地球电影免费观看 | 亚洲精品色激情综合| 伦理电影大哥的女人| 高清欧美精品videossex| 亚洲av成人精品一二三区| 男的添女的下面高潮视频| 在线天堂最新版资源| 国产成人午夜福利电影在线观看| 一个人免费看片子| 中文字幕制服av| 免费观看的影片在线观看| 最近手机中文字幕大全| 亚洲,一卡二卡三卡| 欧美国产精品一级二级三级 | 国产美女午夜福利| 欧美 亚洲 国产 日韩一| 少妇裸体淫交视频免费看高清| 9色porny在线观看| 欧美国产精品一级二级三级 | 久久精品久久久久久噜噜老黄| 人人妻人人添人人爽欧美一区卜| 久久精品国产鲁丝片午夜精品| 成人免费观看视频高清| 精品一区二区免费观看| 色婷婷久久久亚洲欧美| 亚洲无线观看免费| 亚洲欧美成人综合另类久久久| av国产精品久久久久影院| 欧美成人精品欧美一级黄| 日本黄色日本黄色录像| 亚洲国产av新网站| 久久久久网色| 18+在线观看网站| 在线观看免费日韩欧美大片 | 美女视频免费永久观看网站| 久久女婷五月综合色啪小说| 啦啦啦在线观看免费高清www| 国产精品一区www在线观看| 国产色婷婷99| 亚洲精品,欧美精品| 综合色丁香网| 中文在线观看免费www的网站| 一级毛片黄色毛片免费观看视频| 2018国产大陆天天弄谢| 日本黄色片子视频| 精品国产一区二区三区久久久樱花| 日韩大片免费观看网站| 一级黄片播放器| 国产又色又爽无遮挡免| 婷婷色综合大香蕉| 国产精品熟女久久久久浪| 国产精品一区www在线观看| 丝袜在线中文字幕| 大香蕉久久网| 国产一区二区三区综合在线观看 | 看十八女毛片水多多多| 亚州av有码| a 毛片基地| 国产精品国产av在线观看| 日产精品乱码卡一卡2卡三| 色视频在线一区二区三区| 中文字幕制服av| 欧美日韩综合久久久久久| 女人精品久久久久毛片| 久久精品夜色国产| 亚洲成人av在线免费| 精品午夜福利在线看| 免费黄色在线免费观看| 夫妻性生交免费视频一级片| 久久久久人妻精品一区果冻| av黄色大香蕉| freevideosex欧美| 欧美日韩国产mv在线观看视频| 欧美精品高潮呻吟av久久| 男女边摸边吃奶| 亚洲欧美日韩卡通动漫| 亚洲国产精品国产精品| 最新的欧美精品一区二区| 国产成人免费无遮挡视频| 精品亚洲成a人片在线观看| 天天操日日干夜夜撸| 青春草亚洲视频在线观看| 自线自在国产av| 国产女主播在线喷水免费视频网站| 亚洲情色 制服丝袜| 蜜桃久久精品国产亚洲av| 日本黄大片高清| 久久影院123| 黑人巨大精品欧美一区二区蜜桃 | 国产高清不卡午夜福利| 在线观看国产h片| 晚上一个人看的免费电影| 欧美三级亚洲精品| 18禁在线播放成人免费| 国产高清三级在线| av免费观看日本| 我要看黄色一级片免费的| 少妇被粗大的猛进出69影院 | 一个人看视频在线观看www免费| 亚洲美女视频黄频| 国产亚洲一区二区精品| 亚洲国产精品一区二区三区在线| 久久午夜综合久久蜜桃| 夜夜爽夜夜爽视频| 日韩制服骚丝袜av| 国产淫片久久久久久久久| 久久久欧美国产精品| 又大又黄又爽视频免费| 久久99热6这里只有精品| 夫妻性生交免费视频一级片| 人人妻人人添人人爽欧美一区卜| 久久久久久久亚洲中文字幕| 最黄视频免费看| av女优亚洲男人天堂| 亚洲国产日韩一区二区| 国产国拍精品亚洲av在线观看| 热re99久久精品国产66热6| 你懂的网址亚洲精品在线观看| 久久午夜综合久久蜜桃| 青春草视频在线免费观看| 街头女战士在线观看网站| 插阴视频在线观看视频| 看免费成人av毛片| 中文精品一卡2卡3卡4更新| 在线观看www视频免费| 免费观看av网站的网址| 免费观看无遮挡的男女| 国产伦理片在线播放av一区| 国产成人精品婷婷| 国产91av在线免费观看| 免费高清在线观看视频在线观看| 黄色怎么调成土黄色| 如何舔出高潮| 日本黄色日本黄色录像| 如日韩欧美国产精品一区二区三区 | videossex国产| 美女中出高潮动态图| 国产黄频视频在线观看| 91在线精品国自产拍蜜月| 最近的中文字幕免费完整| 99热这里只有是精品50| 色网站视频免费| 日日啪夜夜撸| 欧美另类一区| 中文字幕精品免费在线观看视频 | 夫妻午夜视频| 日韩亚洲欧美综合| 91成人精品电影| 国产精品一区二区在线观看99| 少妇人妻久久综合中文| 2018国产大陆天天弄谢| 我的女老师完整版在线观看| 最近手机中文字幕大全| 精品视频人人做人人爽| 欧美亚洲 丝袜 人妻 在线| av国产精品久久久久影院| 高清毛片免费看| 久久精品国产亚洲av涩爱| 国产午夜精品久久久久久一区二区三区| 美女xxoo啪啪120秒动态图| 老熟女久久久| 男人舔奶头视频| 毛片一级片免费看久久久久| 十八禁高潮呻吟视频 | 久久久国产一区二区| 99热全是精品| 亚洲国产日韩一区二区| 一级毛片 在线播放| 一级毛片我不卡| 精品亚洲成a人片在线观看| 岛国毛片在线播放| 久久精品国产自在天天线| 国产精品三级大全| 人妻人人澡人人爽人人| 韩国高清视频一区二区三区| av视频免费观看在线观看| 天天操日日干夜夜撸| 最近最新中文字幕免费大全7| 国产av国产精品国产| 内地一区二区视频在线| 日本黄色日本黄色录像| 国产一区二区三区av在线| 一二三四中文在线观看免费高清| 久久久久久人妻| 波野结衣二区三区在线| 久久久久精品久久久久真实原创| 一级毛片久久久久久久久女| 欧美3d第一页| 青春草视频在线免费观看| 午夜影院在线不卡| 看十八女毛片水多多多| av一本久久久久| 亚洲图色成人| 日韩成人伦理影院| 免费观看a级毛片全部| 亚洲va在线va天堂va国产| 欧美精品高潮呻吟av久久| 你懂的网址亚洲精品在线观看| a级毛片免费高清观看在线播放| 国产69精品久久久久777片| 五月开心婷婷网| 国产精品三级大全| 韩国av在线不卡| 岛国毛片在线播放| 国产成人91sexporn| 一区二区三区乱码不卡18| videos熟女内射| 日韩一本色道免费dvd| 嫩草影院新地址| tube8黄色片| 日本黄大片高清| √禁漫天堂资源中文www| 在现免费观看毛片| 亚洲av成人精品一区久久| 全区人妻精品视频| 国产精品三级大全| 久久99热6这里只有精品| 日韩视频在线欧美| 一级,二级,三级黄色视频| 有码 亚洲区| 日韩一区二区三区影片| 九草在线视频观看| 色婷婷av一区二区三区视频| 美女福利国产在线| 国产日韩欧美亚洲二区| 日韩av不卡免费在线播放| 91久久精品国产一区二区成人| 亚洲精华国产精华液的使用体验| 狂野欧美激情性xxxx在线观看| 男人添女人高潮全过程视频| 久久久国产一区二区| 一个人免费看片子| 亚洲自偷自拍三级| 精品一区二区三卡| 日本猛色少妇xxxxx猛交久久| 亚洲国产毛片av蜜桃av| 国产精品成人在线| 2021少妇久久久久久久久久久| 伊人久久精品亚洲午夜| 五月玫瑰六月丁香| 狂野欧美激情性bbbbbb| 97精品久久久久久久久久精品| 午夜久久久在线观看| 妹子高潮喷水视频| 人体艺术视频欧美日本| 久久久精品94久久精品| 又爽又黄a免费视频| 亚洲av中文av极速乱| 午夜福利视频精品| 精品一区二区三卡| 深夜a级毛片| 中文字幕av电影在线播放| 老女人水多毛片| 男女免费视频国产| 搡老乐熟女国产| 国产精品一二三区在线看| 国产白丝娇喘喷水9色精品| 午夜久久久在线观看| 午夜日本视频在线| 国产精品国产三级国产av玫瑰| 大香蕉97超碰在线| 欧美+日韩+精品| 一区二区av电影网| 黄色视频在线播放观看不卡| 国产伦精品一区二区三区四那| 日日啪夜夜撸| 精品久久久精品久久久| 另类亚洲欧美激情| 亚洲伊人久久精品综合| 麻豆成人午夜福利视频| 大又大粗又爽又黄少妇毛片口| av在线观看视频网站免费| 精品少妇内射三级| 在线观看人妻少妇| 男女啪啪激烈高潮av片| 9色porny在线观看| 国产色爽女视频免费观看| 美女大奶头黄色视频| 久久久久人妻精品一区果冻| 特大巨黑吊av在线直播| 熟女电影av网| 久久午夜综合久久蜜桃| 涩涩av久久男人的天堂| 国产亚洲5aaaaa淫片| 中文字幕av电影在线播放| tube8黄色片| 国产一区二区三区av在线| 婷婷色综合www| 国产精品.久久久| 黑人猛操日本美女一级片| 日本91视频免费播放| av福利片在线| 午夜激情福利司机影院| www.av在线官网国产| 久久久国产一区二区| 国产一级毛片在线| 大话2 男鬼变身卡| 日日啪夜夜撸| 大片电影免费在线观看免费| 欧美丝袜亚洲另类| 人体艺术视频欧美日本| 人妻一区二区av| 亚洲国产精品一区二区三区在线| 一级片'在线观看视频| 晚上一个人看的免费电影| 寂寞人妻少妇视频99o| 免费av中文字幕在线| 国产精品国产三级国产专区5o| 国产精品99久久久久久久久| 汤姆久久久久久久影院中文字幕| 亚洲国产精品一区三区| 亚洲精品中文字幕在线视频 | 中文字幕免费在线视频6| av网站免费在线观看视频| 国产成人91sexporn| 天堂8中文在线网| 久久久久网色| 亚洲国产毛片av蜜桃av| 免费观看在线日韩| 中国三级夫妇交换| 男男h啪啪无遮挡| 中文字幕制服av| 亚洲中文av在线| 日韩不卡一区二区三区视频在线| 卡戴珊不雅视频在线播放| 亚洲av在线观看美女高潮| 黑丝袜美女国产一区| 日韩大片免费观看网站| 晚上一个人看的免费电影| 亚洲精品亚洲一区二区| 午夜福利影视在线免费观看| 伦精品一区二区三区| 免费久久久久久久精品成人欧美视频 | 免费在线观看成人毛片| 两个人免费观看高清视频 | 欧美xxxx性猛交bbbb| 人人妻人人澡人人爽人人夜夜| 国产精品国产三级国产av玫瑰| 亚洲国产欧美日韩在线播放 | 一本一本综合久久| 99久久精品一区二区三区| av播播在线观看一区| 国产精品久久久久成人av| 在线精品无人区一区二区三| 久久久久视频综合| 国产视频内射| 免费观看av网站的网址| 嫩草影院新地址| 少妇被粗大的猛进出69影院 | 国产中年淑女户外野战色| 人人妻人人看人人澡| 日韩制服骚丝袜av| 亚洲人与动物交配视频| 99精国产麻豆久久婷婷| 亚洲精品久久午夜乱码| 日日摸夜夜添夜夜添av毛片| 久久人人爽人人爽人人片va| 好男人视频免费观看在线| 51国产日韩欧美| 肉色欧美久久久久久久蜜桃| 国产免费一区二区三区四区乱码| 国产毛片在线视频| 亚洲精品国产成人久久av| 黄色日韩在线| 成年美女黄网站色视频大全免费 | 欧美激情国产日韩精品一区| 99热网站在线观看| 亚洲国产欧美日韩在线播放 | 香蕉精品网在线| 永久网站在线| 大片免费播放器 马上看| 午夜福利,免费看| 欧美另类一区| 免费少妇av软件| 桃花免费在线播放| 亚洲人成网站在线播| 免费观看a级毛片全部| 韩国高清视频一区二区三区| 男女边摸边吃奶| 国产精品久久久久久精品电影小说| 精品国产露脸久久av麻豆| 少妇猛男粗大的猛烈进出视频| 中文资源天堂在线| 亚洲精品国产色婷婷电影| h视频一区二区三区| 日韩三级伦理在线观看| 久久精品国产鲁丝片午夜精品| 久久精品国产自在天天线| 免费看日本二区| 国产精品偷伦视频观看了| 免费在线观看成人毛片| 日韩一本色道免费dvd| 简卡轻食公司| 久久久久久久久久久久大奶| 亚洲精品日本国产第一区| 国产免费福利视频在线观看| 在线 av 中文字幕| 国产熟女欧美一区二区| 成人无遮挡网站| 欧美区成人在线视频| 18禁在线无遮挡免费观看视频| 不卡视频在线观看欧美| 国产高清国产精品国产三级| 两个人免费观看高清视频 | 男人添女人高潮全过程视频| 黑人高潮一二区| 国产色爽女视频免费观看| 一级毛片 在线播放| av黄色大香蕉| 日日撸夜夜添| 观看免费一级毛片| 我的女老师完整版在线观看| 极品少妇高潮喷水抽搐| 国内揄拍国产精品人妻在线| 国产精品伦人一区二区| 国产极品天堂在线| 亚洲欧美日韩另类电影网站| 日韩欧美 国产精品| 热99国产精品久久久久久7| 肉色欧美久久久久久久蜜桃| 免费观看的影片在线观看| 欧美日韩一区二区视频在线观看视频在线| 久久午夜福利片| av福利片在线观看| 精品国产国语对白av| 久久女婷五月综合色啪小说| 国精品久久久久久国模美| 日韩熟女老妇一区二区性免费视频| 卡戴珊不雅视频在线播放| 国产欧美亚洲国产| 亚洲色图综合在线观看| 国产成人一区二区在线| 中国国产av一级| 午夜视频国产福利| 国产视频首页在线观看| 91精品一卡2卡3卡4卡| 国产精品一区二区在线不卡| 日韩在线高清观看一区二区三区| 少妇猛男粗大的猛烈进出视频| 欧美精品一区二区大全| 亚洲人成网站在线播| 精品久久国产蜜桃| 三级经典国产精品| 能在线免费看毛片的网站| 久久人妻熟女aⅴ| 曰老女人黄片| 欧美精品人与动牲交sv欧美| 简卡轻食公司| 国产亚洲午夜精品一区二区久久| tube8黄色片| 97超视频在线观看视频| 夜夜爽夜夜爽视频| √禁漫天堂资源中文www| 高清视频免费观看一区二区| 国产精品久久久久成人av| 国产一区二区三区av在线| 日本91视频免费播放| 国产男人的电影天堂91| 99久久精品一区二区三区| 国产视频首页在线观看| 91精品一卡2卡3卡4卡| 99九九线精品视频在线观看视频| 日韩成人伦理影院| 一二三四中文在线观看免费高清| 亚洲av二区三区四区| 美女cb高潮喷水在线观看| 在线看a的网站| 国产精品久久久久久精品电影小说| 日韩人妻高清精品专区| 青春草国产在线视频| 亚洲天堂av无毛| 久久久久久人妻| 精品一区二区三卡| 成人亚洲欧美一区二区av| 亚洲国产av新网站| tube8黄色片| 九色成人免费人妻av| 免费看av在线观看网站| www.色视频.com| 国产伦精品一区二区三区四那| a级毛片免费高清观看在线播放| 亚洲精品一二三| 久久久久久久精品精品| 色94色欧美一区二区| 一级毛片我不卡| 久久久久久久久久人人人人人人| 下体分泌物呈黄色| 美女主播在线视频| 91aial.com中文字幕在线观看| 嘟嘟电影网在线观看| 午夜日本视频在线| 亚洲不卡免费看| 亚洲内射少妇av| 久久久久国产网址| 99九九在线精品视频 | 男女边吃奶边做爰视频| 久久亚洲国产成人精品v| 啦啦啦视频在线资源免费观看| 免费av不卡在线播放| a 毛片基地| 在线观看免费日韩欧美大片 | 女的被弄到高潮叫床怎么办| 少妇的逼水好多| 自拍欧美九色日韩亚洲蝌蚪91 | 高清av免费在线| 亚洲无线观看免费| 在线观看免费高清a一片| 国产av码专区亚洲av| 国产美女午夜福利| 欧美日韩视频精品一区| 亚洲精品乱码久久久v下载方式| 五月玫瑰六月丁香| 在线播放无遮挡| 热re99久久国产66热| 亚洲国产精品一区二区三区在线| 麻豆乱淫一区二区| 搡女人真爽免费视频火全软件|