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

    美女被艹到高潮喷水动态| 黄片小视频在线播放| 久久久久久国产a免费观看| 欧美成人一区二区免费高清观看| www.熟女人妻精品国产| 国产免费男女视频| 美女大奶头视频| 床上黄色一级片| 免费观看的影片在线观看| 亚洲男人的天堂狠狠| 亚洲,欧美,日韩| 97碰自拍视频| 露出奶头的视频| 99在线视频只有这里精品首页| 亚洲不卡免费看| 99国产精品一区二区三区| 欧美一区二区精品小视频在线| 国产精品免费一区二区三区在线| 久久性视频一级片| 日本一本二区三区精品| 亚洲中文字幕日韩| 欧美黑人欧美精品刺激| 麻豆国产av国片精品| 天堂影院成人在线观看| 亚洲成人免费电影在线观看| 亚洲,欧美精品.| 中亚洲国语对白在线视频| 亚洲欧美精品综合久久99| 亚洲无线观看免费| 18禁黄网站禁片免费观看直播| 舔av片在线| bbb黄色大片| 国产探花在线观看一区二区| 久久久久久九九精品二区国产| 国产成年人精品一区二区| 国模一区二区三区四区视频| 国产大屁股一区二区在线视频| 国产午夜精品论理片| 久久香蕉精品热| 婷婷六月久久综合丁香| 欧美xxxx黑人xx丫x性爽| 欧美极品一区二区三区四区| 国产黄片美女视频| 国产高清激情床上av| 亚洲欧美日韩卡通动漫| 欧美极品一区二区三区四区| 亚洲人成网站在线播放欧美日韩| 精品乱码久久久久久99久播| 最近在线观看免费完整版| 日本一本二区三区精品| 亚洲 国产 在线| 国产精品爽爽va在线观看网站| 欧美日韩福利视频一区二区| 欧美最黄视频在线播放免费| 国产乱人视频| a级毛片a级免费在线| 日韩有码中文字幕| 日日摸夜夜添夜夜添av毛片 | 又粗又爽又猛毛片免费看| 一二三四社区在线视频社区8| 97热精品久久久久久| 亚洲不卡免费看| 亚洲专区国产一区二区| 美女高潮喷水抽搐中文字幕| 亚洲欧美日韩无卡精品| 九色成人免费人妻av| 欧美色视频一区免费| 天堂网av新在线| 动漫黄色视频在线观看| 亚洲激情在线av| 亚洲熟妇中文字幕五十中出| 亚洲av成人av| 久久精品国产亚洲av香蕉五月| 国内揄拍国产精品人妻在线| 国产91精品成人一区二区三区| 长腿黑丝高跟| 18禁黄网站禁片午夜丰满| 国产精品野战在线观看| 亚洲成人中文字幕在线播放| 成人无遮挡网站| av黄色大香蕉| 97碰自拍视频| 亚洲国产精品成人综合色| 午夜福利成人在线免费观看| 亚洲av五月六月丁香网| 亚洲av免费在线观看| 国产亚洲精品久久久com| 亚洲av成人精品一区久久| 动漫黄色视频在线观看| 黄色视频,在线免费观看| 最新在线观看一区二区三区| 网址你懂的国产日韩在线| 国产欧美日韩一区二区三| 午夜老司机福利剧场| 精品久久久久久久久久免费视频| 日韩免费av在线播放| 国产亚洲欧美在线一区二区| xxxwww97欧美| h日本视频在线播放| 久久久成人免费电影| 一边摸一边抽搐一进一小说| 可以在线观看毛片的网站| 亚洲专区中文字幕在线| 69人妻影院| 99热这里只有精品一区| 不卡一级毛片| 欧美性感艳星| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 90打野战视频偷拍视频| 99久国产av精品| 国内精品一区二区在线观看| 成人av一区二区三区在线看| 国产69精品久久久久777片| bbb黄色大片| 天堂影院成人在线观看| 婷婷亚洲欧美| 亚洲真实伦在线观看| 琪琪午夜伦伦电影理论片6080| 午夜精品久久久久久毛片777| 亚洲在线观看片| 日韩精品中文字幕看吧| 久久国产乱子免费精品| 亚洲av日韩精品久久久久久密| 男女视频在线观看网站免费| 久久久久性生活片| 国产亚洲精品av在线| av天堂在线播放| 亚洲无线在线观看| 午夜福利在线观看吧| 中出人妻视频一区二区| 亚洲精品乱码久久久v下载方式| 老女人水多毛片| 高潮久久久久久久久久久不卡| av专区在线播放| 亚洲色图av天堂| av中文乱码字幕在线| 精品久久久久久久久亚洲 | avwww免费| 天堂√8在线中文| 一个人免费在线观看电影| 久久久久国内视频| 99热这里只有是精品在线观看 | 免费电影在线观看免费观看| 国产精品免费一区二区三区在线| 91在线观看av| 欧美激情国产日韩精品一区| 91久久精品电影网| 欧美在线黄色| 国产主播在线观看一区二区| 国产精品免费一区二区三区在线| 免费电影在线观看免费观看| 乱人视频在线观看| 90打野战视频偷拍视频| 亚洲五月婷婷丁香| bbb黄色大片| 国产色爽女视频免费观看| 国产综合懂色| 搞女人的毛片| 在线天堂最新版资源| 欧美日本亚洲视频在线播放| 一a级毛片在线观看| 欧美精品国产亚洲| 天天躁日日操中文字幕| 三级毛片av免费| 真实男女啪啪啪动态图| av福利片在线观看| 国产精品日韩av在线免费观看| 成年女人毛片免费观看观看9| 亚洲经典国产精华液单 | 亚洲美女搞黄在线观看 | 亚洲经典国产精华液单 | 亚洲最大成人中文| 欧美精品国产亚洲| 亚洲第一欧美日韩一区二区三区| 精品一区二区三区视频在线观看免费| 日本免费a在线| 别揉我奶头 嗯啊视频| 欧美一区二区精品小视频在线| 亚洲七黄色美女视频| 久久久久久久亚洲中文字幕 | 99久久成人亚洲精品观看| 三级国产精品欧美在线观看| 一区二区三区四区激情视频 | 在线播放国产精品三级| 两个人的视频大全免费| 啪啪无遮挡十八禁网站| 久久久久久久久中文| 少妇丰满av| 男人和女人高潮做爰伦理| 国产精品人妻久久久久久| 欧美bdsm另类| 国产真实乱freesex| 亚洲va日本ⅴa欧美va伊人久久| 久久精品国产清高在天天线| 欧美黄色淫秽网站| 国产精品综合久久久久久久免费| 欧美黑人巨大hd| 免费人成在线观看视频色| 男人舔奶头视频| 欧美成人a在线观看| 免费在线观看日本一区| 色综合亚洲欧美另类图片| 久99久视频精品免费| 亚洲va日本ⅴa欧美va伊人久久| 五月玫瑰六月丁香| 亚洲avbb在线观看| 午夜福利在线在线| 亚洲成av人片免费观看| 成人国产综合亚洲| 美女xxoo啪啪120秒动态图 | or卡值多少钱| 国产高清三级在线| 男女做爰动态图高潮gif福利片| 国产麻豆成人av免费视频| 757午夜福利合集在线观看| 精品午夜福利在线看| 婷婷丁香在线五月| 午夜免费成人在线视频| 亚洲成av人片免费观看| av天堂中文字幕网| 听说在线观看完整版免费高清| 国产av一区在线观看免费| 91午夜精品亚洲一区二区三区 | 国产精品永久免费网站| 91av网一区二区| 麻豆成人午夜福利视频| 又紧又爽又黄一区二区| 观看美女的网站| 日本一本二区三区精品| 美女大奶头视频| 国产欧美日韩一区二区精品| 在线观看66精品国产| 淫秽高清视频在线观看| 观看免费一级毛片| 最新中文字幕久久久久| 夜夜夜夜夜久久久久| 国产精品亚洲一级av第二区| 51午夜福利影视在线观看| 美女被艹到高潮喷水动态| 草草在线视频免费看| 无人区码免费观看不卡| 在线观看免费视频日本深夜| 国产高清有码在线观看视频| 免费av观看视频| 欧美午夜高清在线| 人妻久久中文字幕网| 91字幕亚洲| 中亚洲国语对白在线视频| 亚洲国产日韩欧美精品在线观看| 国产国拍精品亚洲av在线观看| 精品日产1卡2卡| 在线观看av片永久免费下载| 亚洲最大成人手机在线| 波野结衣二区三区在线| 69人妻影院| 赤兔流量卡办理| 国产精品久久久久久精品电影| 免费在线观看成人毛片| 国产高潮美女av| 久久亚洲真实| 成人国产综合亚洲| 嫩草影院新地址| 美女被艹到高潮喷水动态| 丰满的人妻完整版| 色av中文字幕| 亚洲成人中文字幕在线播放| 国产69精品久久久久777片| 1000部很黄的大片| 亚洲av免费高清在线观看| 看十八女毛片水多多多| 亚洲欧美日韩卡通动漫| 日韩国内少妇激情av| 麻豆国产97在线/欧美| 丰满人妻熟妇乱又伦精品不卡| 亚洲欧美日韩无卡精品| 美女 人体艺术 gogo| 国产视频一区二区在线看| 亚洲av第一区精品v没综合| 成人精品一区二区免费| 99在线人妻在线中文字幕| 别揉我奶头 嗯啊视频| 亚洲自拍偷在线| 国产精品美女特级片免费视频播放器| 国产一级毛片七仙女欲春2| 亚洲aⅴ乱码一区二区在线播放| 99久久99久久久精品蜜桃| 美女高潮喷水抽搐中文字幕| 久久人妻av系列| 国内揄拍国产精品人妻在线| 中文字幕免费在线视频6| 波多野结衣高清无吗| 1024手机看黄色片| 一个人看的www免费观看视频| www.色视频.com| 国产成人影院久久av| 性欧美人与动物交配| 国产免费男女视频| 日本 av在线| 亚洲国产色片| 国产黄片美女视频| 久久久久久大精品| 男人狂女人下面高潮的视频| 国内精品一区二区在线观看| av天堂在线播放| 国产成人福利小说| 国产成人av教育| 中文字幕久久专区| 国产午夜福利久久久久久| 国产私拍福利视频在线观看| 一级毛片久久久久久久久女| 桃色一区二区三区在线观看| 日韩精品中文字幕看吧| 久久精品国产自在天天线| 成人午夜高清在线视频| 成人亚洲精品av一区二区| 日本黄色视频三级网站网址| av视频在线观看入口| 在线播放无遮挡| 少妇的逼好多水| 午夜久久久久精精品| 亚洲片人在线观看| 久久伊人香网站| 少妇丰满av| 如何舔出高潮| ponron亚洲| 国产欧美日韩一区二区精品| 99久久久亚洲精品蜜臀av| 青草久久国产| 亚洲av二区三区四区| 成人午夜高清在线视频| 99国产精品一区二区三区| 欧美一区二区亚洲| 国产在线精品亚洲第一网站| 在线国产一区二区在线| 中文字幕精品亚洲无线码一区| 久久草成人影院| 在线播放无遮挡| 亚洲国产高清在线一区二区三| 午夜福利免费观看在线| 免费看日本二区| 国产视频一区二区在线看| 一二三四社区在线视频社区8| 久久九九热精品免费| 毛片女人毛片| 久久草成人影院| 国产伦在线观看视频一区| 精品久久久久久久久av| 高清在线国产一区| 一本综合久久免费| 九九久久精品国产亚洲av麻豆| 亚洲五月婷婷丁香| 99久久99久久久精品蜜桃| 亚洲成人精品中文字幕电影| 欧美日韩乱码在线| 久久精品国产亚洲av天美| av国产免费在线观看| 国产真实乱freesex| 日韩大尺度精品在线看网址| 成人av一区二区三区在线看| 欧美成人性av电影在线观看| 可以在线观看毛片的网站| 国产蜜桃级精品一区二区三区| 好男人在线观看高清免费视频| 看黄色毛片网站| 一区福利在线观看| 久久99热6这里只有精品| 色哟哟哟哟哟哟| 亚洲内射少妇av| 精品久久久久久成人av| 成人av一区二区三区在线看| 亚洲aⅴ乱码一区二区在线播放| 99久久精品一区二区三区| 亚洲激情在线av| 丰满人妻熟妇乱又伦精品不卡| 国产男靠女视频免费网站| 免费人成视频x8x8入口观看| 欧美最黄视频在线播放免费| 人人妻人人澡欧美一区二区| 最近中文字幕高清免费大全6 | 久久天躁狠狠躁夜夜2o2o| 国产亚洲精品av在线| www日本黄色视频网| 免费观看的影片在线观看| 国产精品国产高清国产av| 亚洲18禁久久av| 国产精品永久免费网站| 免费看a级黄色片| 91在线精品国自产拍蜜月| 久久人妻av系列| av中文乱码字幕在线| 国产高清有码在线观看视频| 免费人成在线观看视频色| 国产高清激情床上av| 美女高潮的动态| 全区人妻精品视频| 欧美性猛交╳xxx乱大交人| 观看免费一级毛片| 男人舔女人下体高潮全视频| 18美女黄网站色大片免费观看| 小蜜桃在线观看免费完整版高清| 亚洲av成人不卡在线观看播放网| 91字幕亚洲| 极品教师在线免费播放| 欧美黑人欧美精品刺激| 国产亚洲欧美在线一区二区| 国产精品久久视频播放| 亚洲av成人av| 久久天躁狠狠躁夜夜2o2o| 麻豆av噜噜一区二区三区| 久久天躁狠狠躁夜夜2o2o| 亚洲人成网站在线播| 韩国av一区二区三区四区| 亚洲 欧美 日韩 在线 免费| 久久这里只有精品中国| 日韩欧美免费精品| 精品无人区乱码1区二区| 国产精品电影一区二区三区| 一级黄色大片毛片| 香蕉av资源在线| 久久久久国内视频| 乱码一卡2卡4卡精品| 国产三级在线视频| 亚洲五月婷婷丁香| 国产精品亚洲美女久久久| 婷婷精品国产亚洲av在线| 一进一出好大好爽视频| 国产成人影院久久av| 久久久精品大字幕| 欧美日韩瑟瑟在线播放| 少妇被粗大猛烈的视频| 免费大片18禁| 女人被狂操c到高潮| 国产伦在线观看视频一区| 神马国产精品三级电影在线观看| 国产精品三级大全| 久久久久国内视频| av黄色大香蕉| 亚洲国产精品合色在线| 日韩欧美精品v在线| а√天堂www在线а√下载| 婷婷丁香在线五月| 成年版毛片免费区| 久久亚洲精品不卡| 成人无遮挡网站| 日本三级黄在线观看| 在线观看免费视频日本深夜| 丰满人妻一区二区三区视频av| 亚洲最大成人手机在线| 国产探花在线观看一区二区| 日韩精品青青久久久久久| 极品教师在线免费播放| 99久国产av精品| 亚洲欧美日韩东京热| av国产免费在线观看| 久久亚洲真实| 久久国产精品人妻蜜桃| ponron亚洲| 18禁在线播放成人免费| 精品久久久久久成人av| 99视频精品全部免费 在线| 国产国拍精品亚洲av在线观看| bbb黄色大片| 国产午夜精品论理片| eeuss影院久久| a在线观看视频网站| 啦啦啦韩国在线观看视频| 别揉我奶头 嗯啊视频| 伦理电影大哥的女人| 国产精品久久久久久精品电影| 午夜福利在线观看吧| 午夜久久久久精精品| 女同久久另类99精品国产91| 一二三四社区在线视频社区8| 国产成人aa在线观看| av专区在线播放| 成人av一区二区三区在线看| 有码 亚洲区| 两人在一起打扑克的视频| 国产精品国产高清国产av| 国产伦精品一区二区三区视频9| 男插女下体视频免费在线播放| 欧美潮喷喷水| 亚洲国产精品999在线| 两个人视频免费观看高清| 桃色一区二区三区在线观看| 成年版毛片免费区| 国产高清激情床上av| 麻豆国产av国片精品| 精品久久久久久久久久免费视频| 成人国产一区最新在线观看| 搡老妇女老女人老熟妇| 国产久久久一区二区三区| 久久久久久九九精品二区国产| 俄罗斯特黄特色一大片| 直男gayav资源| 国产精品综合久久久久久久免费| 偷拍熟女少妇极品色| www.熟女人妻精品国产| 老熟妇仑乱视频hdxx| 午夜亚洲福利在线播放| 黄色丝袜av网址大全| 久久久久久九九精品二区国产| 日日干狠狠操夜夜爽| 淫妇啪啪啪对白视频| 少妇人妻一区二区三区视频| www.www免费av| 人人妻人人看人人澡| 国产免费av片在线观看野外av| 久99久视频精品免费| 一个人看的www免费观看视频| 精品久久国产蜜桃| 一卡2卡三卡四卡精品乱码亚洲| 国产午夜福利久久久久久| av在线天堂中文字幕| 一个人看视频在线观看www免费| 色综合亚洲欧美另类图片| 亚洲久久久久久中文字幕| 又黄又爽又免费观看的视频| 国产高清激情床上av| 非洲黑人性xxxx精品又粗又长| 黄色日韩在线| 欧美一区二区国产精品久久精品| 内地一区二区视频在线| 日韩免费av在线播放| 精品久久久久久久末码| 亚洲成av人片在线播放无| 国产精品,欧美在线| 欧美3d第一页| 国产高清视频在线播放一区| 亚洲五月天丁香| 在线观看舔阴道视频| 亚洲av免费在线观看| 欧美成人免费av一区二区三区| 国内少妇人妻偷人精品xxx网站| 18美女黄网站色大片免费观看| 国产熟女xx| 欧美精品国产亚洲| 久久久久精品国产欧美久久久| 亚洲乱码一区二区免费版| 国产精品嫩草影院av在线观看 | 国产综合懂色| 国产不卡一卡二| 欧美另类亚洲清纯唯美| 日韩欧美国产一区二区入口| 久久午夜亚洲精品久久| 亚洲色图av天堂| 搡老妇女老女人老熟妇| 成人精品一区二区免费| 久久中文看片网| 欧美高清性xxxxhd video| av在线观看视频网站免费| 好男人在线观看高清免费视频| 三级男女做爰猛烈吃奶摸视频| 免费观看人在逋| 免费av观看视频| 十八禁网站免费在线| 亚洲天堂国产精品一区在线| 午夜a级毛片| 两人在一起打扑克的视频| 91字幕亚洲| 亚洲国产色片| 全区人妻精品视频| 亚洲黑人精品在线| 欧美日韩国产亚洲二区| 国产成人福利小说| 成人美女网站在线观看视频| 久久精品国产亚洲av香蕉五月| 首页视频小说图片口味搜索| 国产三级中文精品| 国产欧美日韩一区二区精品| 国产视频一区二区在线看| 啪啪无遮挡十八禁网站| 久久久久久久久久成人| a级毛片免费高清观看在线播放| 精品国内亚洲2022精品成人| 黄片小视频在线播放| 级片在线观看| 精品国产三级普通话版| 男人狂女人下面高潮的视频| 国产精品久久久久久亚洲av鲁大| 精品一区二区三区视频在线| 三级国产精品欧美在线观看| 俺也久久电影网| 免费av不卡在线播放| 亚洲av熟女| 欧美绝顶高潮抽搐喷水| 听说在线观看完整版免费高清| 精品久久久久久久人妻蜜臀av| 国产欧美日韩精品亚洲av| 亚洲在线自拍视频| 国产精品一及| 尤物成人国产欧美一区二区三区| 搡女人真爽免费视频火全软件 | 久久久久久久久中文| 久久久久久久午夜电影| 可以在线观看毛片的网站| 国产高清视频在线观看网站| 天堂影院成人在线观看| 国产伦精品一区二区三区视频9| 欧美高清性xxxxhd video| 午夜a级毛片| 亚洲精品乱码久久久v下载方式| 午夜福利欧美成人| 中出人妻视频一区二区| 日韩欧美精品v在线| 老司机午夜福利在线观看视频| 国产真实伦视频高清在线观看 | 精品一区二区三区视频在线| 两个人的视频大全免费| 真人一进一出gif抽搐免费| av在线蜜桃| 国产欧美日韩一区二区精品| 每晚都被弄得嗷嗷叫到高潮| 精品熟女少妇八av免费久了| 91午夜精品亚洲一区二区三区 |