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

    Distinguishing the Regional Atmospheric Controls on Precipitation Isotopic Variability in the Central-Southeast Portion of Brazil

    2022-10-27 09:45:02VinciusdosSANTOSPeterMARSHALLFLEMINGLuHENRIQUEMANCINIStelaDALVASANTOSCOTAGrazieleBeatrizdeLIMARafaelaRODRIGUESGOMESRobertoEduardoKIRCHHEIMRicardoSANCHMURILLOandDidierGASTMANS
    Advances in Atmospheric Sciences 2022年10期

    Vinícius dos SANTOS, Peter MARSHALL FLEMING, Luís HENRIQUE MANCINI,Stela DALVA SANTOS COTA, Graziele Beatriz de LIMA,Rafaela RODRIGUES GOMES, Roberto Eduardo KIRCHHEIM,Ricardo SANCHéZ-MURILLO, and Didier GASTMANS*

    1S?o Paulo State University (UNESP), Environmental Studies Center. Av. 24ABased,1515 - Bela Vista, 13.506-900, Rio Claro (SP), Brazil

    2Center for the Development of the Nuclear Technology (CDTN - CNEN), Av. Presidente Ant?nio Carlos,6.627, Campus da UFMG - Pampulha, 31270-90, Belo Horizonte (MG), Brazil

    3Universidade de Brasília - UnB, Instituto de Geociências, Laboratório de Isótopos Estáveis.Campus "Darcy Ribeiro", Asa Norte, 70.910-900, Brasília (DF), Brasil

    4Hydrology and Territorial Management Directory (DHT), The Geological Survey of Brazil (CPRM-SGB), S?o Paulo (SP), Brazil

    5University of Texas at Arlington, Earth and Environmental Sciences Department,500 Yates Street, Arlington, Texas, 76019, USA

    6Stable Isotopes Research Group and Water Resources Management Laboratory,Universidad Nacional, Heredia, Costa Rica

    ABSTRACT

    Key words: precipitation, stable water isotopes, d-excess, moisture source and transport, southern Atlantic Ocean and Amazon

    1. Introduction

    The stable water isotopologues constitute an important tool to underpin the climatic controls involved in rainfall formation processes, considered one of the most vulnerable components of the hydrological cycle under a changing climate(Power et al., 2017; Wasko and Sharma, 2017). Different isotopologues are formed due to the isotopic fractionation during aqueous phase changes. Therefore, it is possible to trace and combine their isotopic compositions with physical and microphysical processes within distinct water reservoirs at the surface-atmosphere continuum (Bowen et al., 2019).

    Since the establishment of the Global Network of Isotopes in Precipitation in 1961 (known as GNIP; https://www.iaea.org/services/networks/gnip), multiple efforts have addressed the key climatic drivers governing the isotopic composition of precipitation (Rozanski et al., 1993; Clark and Fritz, 1997; Aggarwal et al., 2012; Terzer et al., 2013;Konecky et al., 2019). However, the isotopic drivers are still not well defined in tropical areas, whereby weak correlations with temperature and moderate negative correlations with the amount of precipitation (known as the “amount effect”)(Dansgaard, 1964; Kurita et al., 2009; Winnick et al., 2014;Sánchez-Murillo et al., 2016) are the most recurrent findings. Recently, convective and stratiform precipitation processes (Risi et al., 2008; Aggarwal et al., 2016; Tharammal et al., 2017; Lekshmy et al., 2018; Munksgaard et al., 2019;Sánchez-Murillo et al., 2019) have been considered to be the main factor that explains precipitation isotopic variability in the tropics and a large portion of extratropical areas worldwide (Aggarwal et al., 2016). However, the linear regression model used for rainfall correlation has already been debated(Konecky et al., 2019) since the humidity convergence processes and rain re-evaporation were not considered.

    The combination of more robust analytical instrumentation and techniques combined with the greater efficiency of reliable precipitation passive collectors (Michelsen et al.,2018, 2019) allow for the exploration of precipitation processes at different timescales (i.e., monthly, weekly, daily,and intra-events). Kurita et al. (2009) demonstrated that daily isotopic compositions reflect: 1) the regional atmospheric conditions (i.e., moisture source, rain-out history during transport) such as on a monthly timescale; 2) local effects (i.e., topography, surface meteorological conditions);and 3) post-condensation processes (i.e., the interaction between raindrops and the surrounding water vapor under the cloud and/or re-evaporation during the raindrops descent). It is necessary to emphasize that in daily or subdaily monitoring strategies, the correlation between δ18O and the rain intensity (recognized as the amount effect)tends to be weak (Kurita et al., 2009). This lack of correlation on daily basis (or amount effect) was directly related to regional factors, such as air mass origin, its trajectory, and the precipitation history or local factors (below-cloud processes, which depend on rain intensity, relative humidity,and drop sizes). During intense rainfall events, drops tend to be larger and hence less affected by isotopic exchanges with vapor (Muller et al., 2015; Graf et al., 2019; Han et al.,2020, 2021).

    Our study combines the isotopic compositions of monthly and daily precipitation data, collected over different timescales and obtained from meteorological stations located in the central-southeastern region of Brazil: Belo Horizonte-MG, Brasília-DF, Rio Claro-SP, and Rio de Janeiro-RJ, hereafter BH, BR, RC, and RJ, respectively (Fig. 1 and Table 1). This region is characterized by a high population density and high concentrations of large industrial and agroindustrial activities; therefore, sustaining water resources is one of the main challenges for managers and government agencies (Ferreira et al., 2018). Belo Horizonte, Brasília,and Rio de Janeiro are socio-climatic hot spots (Torres et al.,2012), while Rio Claro is located in the Parana Basin, close to important recharge areas of the Guarani Aquifer System(Batista et al., 2018; Kirchheim et al., 2019; Santarosa et al.,2021).

    The main objective of this study is to evaluate the multitemporal rainfall isotopic variability over the central-southeastern region of Brazil and identify how distinct isotope sampling timescales enhance our understanding of key drivers controlling rainfall generation. The main research questions are 1) How do different climatic regimes/conditions affect the isotopic composition at each location on monthly and daily sampling scales? 2) Do the atmospheric systems promote the same influence over isotopic variability at all stations?3) Does changing the sampling timescale, monthly or daily,affect the insights for isotopic interpretation in tropical areas?

    Our results offer an enhanced isotopic assessment of tropical rain, provide novel information and knowledge to improve the current regional and global interpretations of stable isotope variations in tropical rainfall under climate change scenarios. Furthermore, this study promotes an interpretation of new data using two timescales in a large regional area. From this perspective, the study of the isotopic composition of precipitation in Brazil represents a research area that has received limited attention.

    The remainder of this paper is organized as follows. Section 2 discusses the background settings, section 3 outlines the materials and methods with isotopes analysis and meteorological data, section 4 presents the results in two parts, from monthly to daily timescales, section 5 discusses the results,and section 6 presents the conclusions from this research.

    2. General settings

    The rainfall pattern in central-southeastern Brazil is characterized by an annual average of around 1500 mm and is divided into two distinct seasons, the dry season (April to September) and the wet season (October to March) (Reboita et al., 2010; Alvares et al., 2013). Recycled moisture (i.e.,evapotranspiration) from the Amazon rainforest is the main water vapor source during the wet season (Fig. 1). Amazon moisture is transported by the Low-Level Jet (LLJ), providing moisture flux from the east of the Andes range (i.e., Amazon region) to the southeast region of South America (Marengo et al., 2004; Montini et al., 2019; Ramos et al., 2019). The interactions between the incursions of Bolivian High (BoH),LLJ, and Cold Fronts (CF) on the southeastern coast contribute to the formation of the South American Convergence Zone (SACZ), an area of enhanced convective activity with a wide range of cloud coverage, humidity, and heavy rainfall during the Brazilian summer (Kodama, 1992; Lima et al.,2010) (Fig.1). These climatic features are all part of the seasonal dynamic of the South America Summer Monsoon system (Vera et al., 2006; Reboita et al., 2010; Marengo et al.,2012). During the dry season, there is decreased moisture flux into the Amazon, a displacement of the southern Atlantic Subtropical High (SASH) over the continent, and,consequently, a decrease in rainfall due to strong CF incursions (Reboita et al., 2010).

    3. Materials and methods

    3.1. Isotopic composition of monthly and daily precipitation

    Monthly isotopic compositions at BH, BR, and RJ were obtained from the GNIP dataset (reference code in Table 1,available at: https://www.iaea.org/services/networks/gnip).Both monthly and daily samples were collected in RC at the Environmental Studies Center - S?o Paulo State University(CEA-UNESP) and were analyzed using Cavity Ring-Down Laser Spectroscopy at the Hydrogeology and Hydrochemistry laboratory of UNESP’s Department of Applied Geology,according to dos Santos et al. (2019c), available at https://data.mendeley.com/datasets/72z5cyhg9g/1. The BR daily isotopic composition was collected at the Laboratory of the Geochronological, Geodynamics and Environmental Studies- Brasília University (UnB). The samples were analyzed at the Laboratory of the Stable Isotopes using a Picarro L2120-i Analyzer - Vaporization Module A0211 and Autosampler PAL HTC-xt. The BH daily samples were collected and analyzed at the Center for the Development of the Nuclear Technology (CDTN) using a Delta V IRMS Thermo Fischer system, model Delta V Advantage Isotope Ratio MS.

    Fig. 1. Circulation patterns and atmospheric systems: (a) ET: Evapotranspiration; ITCZ: Intertropical Convergence Zone; BoH: Bolivian High; LLJ: Low-Level Jet; SACZ: South America Convergence Zone;SASH: southern Atlantic Subtropical High; CF: Cold Front (CF1 = of extratropical origin; CF2 = associated with the formation of SACZ). Figure based on Marengo et al. (2004) and modified from Santos et al. (2019a);(b) Stable isotope sampling locations. Annual precipitation (color-coded in mm) was obtained from the dataset of the Climate Prediction Center Merged Analysis of Precipitation, available at: http://www.cpc.ncep.noaa.gov/products/global_precip/html/wpage.cmap.html.

    3.2. Meteorological data and El Ni?o-southern Oscillation (ENSO) analysis

    Precipitation, temperature, and relative humidity surface data from BH (OMM code: 83587, 19.93°S, 43.93°W), BR(OMM code: 83377, 15.78°S, 47.92°W), and RJ (OMM Code: 83743, 22.89°S, 43.18°W) were obtained from the Meteorological Dataset for Education and Research(BDMEP) at the National Institute of Meteorology(INMET), available at: http://www.inmet.gov.br/portal/index.php?r=bdmep/bdmep. Surface meteorological data for RC were obtained from a meteorological station located at the Centro de Análise e Planejamento Ambiental-S?o Paulo State University Rio Claro (22°23′S, 47°32′W). Historical precipitation data ranges are 1961 to 2018 (BH), 1965 to 2018 (BR), 1994 to 2018 (RC), and 1961 to 1983/2002 to 2016 (RJ).

    The influence of ENSO in the study area was assessed using the Oceanic Ni?o Index (ONI), which identifies different ENSO phases and classifies each intensity based on a continuous average of sea surface temperature anomalies (TSM)of three months in the Ni?o region 3.4. over the Pacific Ocean (5°N-5°S, 120°-170°W) (Huang et al., 2017). According to the National Oceanic and Atmospheric Administration(NOAA) (http://origin.cpc.ncep.noaa.gov/products/analysis_monitoring/ensostuff/ONI_v5.php), an ENSO event is defined as SST values exceeding a threshold of ±0.5°C for a minimum of five consecutive overlapping months. Positive and negative anomalies are defined as La Ni?a (cold phase)and El Ni?o (warm phase), respectively.

    Due to different periods covered by the isotopic data(Table 1), locations were separated to compare ENSO events during the same period. At stations BR and RJ, the following ENSO events were selected: 1) warm phase: March 1965-April 1966 and October 1968-December 1969; and 2)cold phase: May 1973-June 1974, October 1974-April 1976). At BH and RC, the following ENSO events were selected: 1) warm phase: November 2014-May 2016; and 2)cold phase: August-December 2016 and October 2017-March 2018. A regression analysis was carried out for these ENSO events, including the monthly isotopic compositions and the ONI-Index. Historical precipitation values were used to determine the standard precipitation anomaly during ENSO phases. At the RC station, the precipitation data (1950-2016) were obtained from dos Santos et al.(2019b) and were updated to include 2017 and 2018.

    3.3. HYSPLIT and synoptic analysis

    The origin of the air masses, their trajectories, and the subsequent precipitation for the daily isotopic composition at BH, BR, and RC were analyzed using the HYSPLIT model (Hybrid Single-Particle Lagrangian Integrated Trajectory) (Stein et al., 2015; Rolph et al., 2017). This model uses a 3D Lagrangian air mass vertical velocity algorithm to determine the position of the air masses and reports these values at hourly time resolution over the trajectory (Soderberg et al., 2013; Sánchez-Murillo et al., 2017). The trajectories of the air masses were estimated for 240 hours prior to precipitation onset, considering the estimated time of residence of the water vapor (Gimeno et al., 2010; van der Ent and Tuinenburg, 2017; Gimeno et al., 2021). The trajectories were computed using NOAA′s meteorological data (global data assimilation system, GDAS: 1 degree, global, 2006-present), with ending elevations of the trajectories at 1500 m above the surface, taking into account the climatological height of the LLJ, within 1000-2000 m (Marengo et al., 2004). The HYSPLIT model also calculates meteorological outputs along the trajectories, such as ambient temperature (K), relative humidity (RH%), rainfall intensity (mm hr-1), and mixed layer depth (mix depth m). The mixing depth is calculated based on the potential temperature profile (Draxler et al.,2020) Trajectory frequency maps were generated as the sum of the number of times the trajectories passed through one 3°grid cell using ArcGIS software.

    The synoptic atmosphere analysis was carried out based on the NCEP/NCAR (National Center for Environmental Prediction/National Center for Atmospheric Research)reanalysis of data from the National Oceanic and Atmospheric Administration - NOAA (Kalnay et al., 1996) (https://psl.noaa.gov/data/composites/day/). This dataset was used to generate plots of vertically integrated precipitable water (kg m-2)and the 500 hPa vertical velocity field or Omega (Pa s-1)with grid size 2.5° × 2.5°, and Outgoing Longwave Radiation(OLR W m-2), grid size 1.0° × 1.0°.

    The OLR is a thermal radiative flux emitted by the earth’s surface and was used as an indicator for the convective activity (Vuille and Werner, 2005), since values below 240 W m-2indicate organized deep convection, while values above 240 W m-2suggest no organized convection occurrence (Gadgil, 2003). The positive and negative Omega corresponded to downward and upward motion (convective activity), respectively.

    Table 1. General information of precipitation sampling locations.

    3.4. Convective and stratiform precipitation by ERAInterim

    The ERA-Interim is a global atmosphere reanalysis from the European Centre for Medium-Range Weather Forecasts (ECMRWF). This model and reanalysis system uses ECMWF’s Integrated Forecast System (IFS-31r2), configured for the following spatial resolution: 60 levels in the vertical, with the top level at 0.1 hPa, T255 spherical-harmonic representation for the basic dynamical fields, and a reduced Gaussian grid with approximately uniform 79 km spacing for the surface and other grid-point fields (Berrisford et al.,2011). The ERA-Interim provides daily convective and large-scale precipitation (stratiform), with 3-hr resolution and a spatial resolution grid of 0.75° × 0.75°.

    The convective and large-scale (stratiform) precipitation data were obtained online (https://apps.ecmwf.int/datasets/data/interim-full-daily/levtype=sfc/). Every grid pixel represents the amount of convective and large-scale precipitation[m of water and was converted to millimeters (mm)] based on the prognostic parameterization of cloud schemes by ECMWF (de Leeuw et al., 2015). Data were obtained for spatial grids corresponding to the BH, BR, and RC stations(Table 1), according to local precipitation intervals recorded by the meteorological stations, which correspond to cumulative precipitation samples from the previous day at 1200 UTC to the next day at 1200 UTC (equivalent to 24 houyrs).In the partitioning of ECMWF scheme by ERA-interim, the large-scale precipitation is considered stratiform precipitation(de Leeuw et al., 2015).

    3.5. Statistical Analysis

    A Kruskal-Wallis (Kruskal and Wallis, 1952) nonparametric test was applied to test statistical differences (p-value< 0.05) between the isotopic compositions (δ18O, δ2H and dexcess) for each locality. The D’Agostino-Pearson K2statistic(n ≥ 20 and k samples) (D’Agostino et al., 1990) was applied to verify that the data distribution was normal (parametric) or non-normal (non-parametric). A significant difference (p-value < 0.05) indicates a nonparametric distribution.A Spearman correlation test was used for nonparametric distribution data, whereas Pearson’s linear correlation test was applied for parametric data. These statistical tests (analysis of variance, normality, and correlation) were conducted between monthly and daily isotope ratios (δ18O, δ2H, and dexcess) and meteorological data (i.e., local station data, ONIIndex, HYSPLIT model, and rainfall type from the ERAInterim). All tests were performed with significance levels defined by a p-value (p) < 0.05, using the Corrplot R package(R Core Team, 2020).

    4. Results

    4.1. Climatic patterns and precipitation

    During the monitoring period (Table 1), the average monthly precipitation (AMPrecip) in all sites presented evident seasonal distributions, following the pattern observed in long-term monthly precipitation (Fig. 2). In terms of the AMPrecip values, December was the rainiest month at BH(275 mm) and BR (236 mm), while at RJ and RC, January featured greater values, at 207 and 204 mm, respectively. In contrast, July was the driest month at RC (18 mm), June in BR(19 mm), June and August in BR (30 mm), and August in RJ (41 mm) (Fig. 2). The wet season represents the largest contribution for annual precipitation, ranging from 83% in BH, 79% in BR, 71% in RC, and 68% in RJ.

    The influence of ENSO on precipitation was determined based on the long-term monthly precipitation (1960-2018)for all sites. Weak annual precipitation anomalies were found during both phases. During the warm phase, rainfall increases were reported in RC (12.6%), RJ (2.5%), and BR(1.9%), and a decrease was noted at BH (-5.6%). During the cold phase, the net rainfall anomalies were as smaller in all sites, with decreases at BH (-5.8%), RC (-1.4%), BR(-0.3%), and a small increase at RJ (0.1%) (Fig. S1 in the Electronic Supplementary Materials, ESM).

    Fig. 2. Monthly averages of meteorological data in Central-Southeast of Brazil. Long-term monthly precipitation(left bars) versus arithmetic mean monthly precipitation records during the isotope sampling period (right bars). The red line is the arithmetic mean temperature during the isotope monitoring period.

    Temperature also presented a seasonal pattern. The difference between the hottest monthly average, January for BH(21°C), RC (24°C), and RJ (26°C), and the coldest monthly average did not exceed 7°C in RC and BH, while at RJ and BR seasonal average temperature only differed by 5°C and 3°C, respectively (Fig. S2 in the ESM).

    Fig. 3. Monthly isotopic compositions of precipitation at each site: (a) δ18O, (b) δ2H, and (c) d-excess dispersion boxplots; (d) Dual δ18O and δ2H plot including the Global Meteoric Water Line (GMWL)for reference and Local Meteoric Water Lines (LMWLs) regressions; (e) monthly average precipitation in bars, and weighted mean isotopes given by lines.

    4.2. Monthly and seasonal isotopic variability

    Differences in the average monthly isotopic composition(δ18O and δ2H) are presented in Figs. 3a-b. Overall, δ18O mean values were: -4.91‰ (BH), -4.58‰ (RC), -4.08‰(BR), and -3.72‰ (RJ). The BR station presented higher isotopic variability, while a lower variation was observed at BH, RC, and RJ. As expected, δ2H values exhibited a similar distribution. Differences were also observed for d-excess(Fig. 3c). Greater d-excess values were observed at RC and RJ sites (d-excess > 12‰), while BR reported lower dexcess values. Amount-weighted δ18O averages were lower than the arithmetic δ18O means for all sites, ranging from-6.25‰ in BH, -5.07‰ in BR, -4.87‰ in RC, and RJ(-2.64‰). Isotope values were plotted along to the Global Meteoric Water Line (GMWL), with the slope of Local Meteoric Water Lines (LMWLs) close to 8 and intercepts greater than 11 for RJ, BH, and RC sites, except for BR that had an intercept near 10‰ (Fig. 3d).

    On a seasonal basis, the isotopic compositions exhibited depleted values during the wet season and enriched values during the dry season, with the most depleted δ18O values occurring during January, March, and December (~ -6‰ to-8‰) and enriched values observed from July to September(~ -4‰ to -1‰). Despite the similarity of the seasonal distributions, the Kruskal-Wallis test revealed statistically significant monthly differences between the sites for δ18O (p =0.03) and d-excess (p = < -0.0001) (Fig. 3e). For δ2H, no significant statistical differences were observed.

    Overall, significant negative correlations between the amount of precipitation and δ18O values were observed across all sites (Table 2). Relationships between the isotopic values and temperature (except in RJ; Table 2) exhibited weak correlations. Elevation and latitudinal effects were also not observed. The most depleted values were recorded at the BH station (857 m elevation). Despite the moderate and positive correlation between the isotopic composition and the ONI-Index during the cold phase of ENSO at BH and RC (Table 2), the p-value is not significant (> 0.05 for all localities). For δ2H and d-excess, the correlations with meteorological parameters are not strong, except for δ2H-temperature (Table 2).

    4.3. Daily isotopic variability

    The assessment of the daily isotopic composition at BH(September 2016-February 2017), BR, and RC (September 2016-March 2017) indicates relatively small differences as demonstrated by the amplitudes (Figs. 4a-c) and the average δ18O values: -4.03‰±3.72‰ (BH), -3.39‰±3.88‰ (BR)and -4.85‰±4.34‰ (RC). Isotopic composition at RC was the most depleted during this period (Figs. 4a and b), while greater d-excess values were observed at BR (Fig. 4c). The Kruskal-Wallis test indicated a statistically significant difference for d-excess among all sites, while for δ18O and δ2H,no significant differences were observed.

    Computed daily LMWLs presented similar slopes and intercepts: BH (δ2H = 7.79 ± 0.16δ18O + 10.74 ± 0.89); BR(δ2H = 7.42 ± 0.12δ18O + 12.68 ± 0.64); RC (δ2H = 7.96 ±0.12δ18O + 12.93 ± 0.80).

    The temporal distribution of the isotopic composition of precipitation is different between locations (Fig. 4). At BH (Fig. 4d) and BR (Fig. 4e), δ18O exhibited two depleted incursions (October-November and February-March).Strong depletions were observed at BR during February 2017. At RC, moderate depletions were observed until late December, with strong incursions between January and March (Fig. 4f). In general, low δ18O values were observed during months with larger precipitation amounts (Fig. 4).Most portions of rain events were characterized by a greater convective precipitation-ERA (Cv) fraction (Fig. S3 in the ESM); however, the Cv and Stratiform (St) precipitation of the ERA does not correlate with δ18O-δ2H values.

    Greater δ18O-δ2H correlations (r ≥ 0.50) were observed with regional meteorological data (Table 3), such as δ18OOmega for BH and BR; δ18O-Precipitable water for BR and RC; δ18O-OLR only for BR; δ18O-Mix depth along HYSPLIT trajectories for BR and RC; the accumulated rainfall along HYSPLIT trajectories presented greater δ18O-correlations for all localities (Table 3).

    The d-excess values were associated with large temporal variability. However, d-excess values can be divided into two periods (Fig. 4g). Until December, a large range of values was observed in RC (4.0‰-24.5‰), BR (1.15‰-20.96‰),and BH (2.09‰-18.12‰), while from January 2017 to March 2017, the observed variability was smaller, as observed at RC (4.2‰-17.1‰) and BH (5.68‰-15.71‰),while BR presented a similar variation (2.94‰-22.72‰).Only Cv (0.60) in BH and local precipitation (0.64) for BR exhibited strong correlations with d-excess (Table 3).

    Table 2. Correlations between precipitation isotopes and meteorological data.

    5. Discussion

    The isotopic composition of precipitation in centralsoutheastern Brazil can be associated with the different processes in the atmosphere operating at different timescales(monthly or daily). Evaluation of the monthly isotopic composition revealed that seasonality is a common rainfall denominator across the region. The depleted (enriched) values (Fig.3c) occurred during the rainy (dry) season (Fig. 2), in agreement with a previous study based on monthly timescales for RC (Gastmans et al., 2017; dos Santos et al., 2019a).

    Fig. 4. Daily isotopic composition of precipitation at each site: (a) δ18O, (b) δ2H, and (c) dexcess dispersion boxplots; Blue bars represent the daily precipitation, dots δ18O daily values for (d) BH, (e) BR, (f) RC, and (g) d-excess values for all sites. The average δ18O value is represented by a dotted line, and depleted values are plotted below this line and indicated with blue circles. In (g) arrows indicate the two periods described in the text.

    Table 3. Correlation between daily isotopic composition and meteorological data. Bolded values are significant (p-value < 0.05) and highlighted color in text are better correlations (r > 0.50).

    Monthly seasonality is not related to a temperature effect since there is a relatively small temperature seasonality between the warmer and colder months (highest value is 7°C) (Rozanski et al., 1993; Vuille et al., 2003; Jasechko,2019). The lack of a temperature effect in tropical regions,where the temperature ranges from 20°C to 30°C, is often attributed to the dominance of the monthly amount effect(Rozanski et al., 1993) and the seasonal changes of the locations of storm origin.

    The influence of ENSO did not result in significant differences in isotopic composition. Despite the existence of positive rainfall anomalies during the dry season (Fig. S1), there were no significant differences in historical rainfall rates or significant correlations between isotopic values and the ONI Index (Table 2). In fact, previous studies revealed that the central-southeastern region of Brazil is a transition region between the drier conditions to the north-northeast and wetter conditions over the southern portion of Brazil (Pezzi and Cavalcanti, 2001; Coelho et al., 2002; Barros et al., 2008; dos Santos et al., 2019b).

    Greater average precipitation amounts were observed at BH and BR compared to RC and RJ, which may be associated with the SACZ position during the wet season (Zilli et al.,2019). Likewise, RC and RJ are the most often influenced by cold-frontal incursions during the dry season, which results in higher precipitation amounts compared to BH and BR. This trend may explain the differences between monthly isotopic compositions, such as the relative depletion in BH and enrichment in RJ, as well as the range of values(Fig. 3). In addition, the small isotope variability in RJ is associated with the proximity to the Atlantic Ocean, whereas in BR, the large variability is more related to the continental effect.

    The distance to the Atlantic Ocean also affects moisture transport, highlighted by the LMWLs. Considering the monthly timescale, higher values of intercepts and slopes observed for BH and RC sites indicate that the isotopic composition of precipitation is influenced by the processes of moisture recirculation in the atmosphere (Wu et al., 2015;Tang et al., 2017; Putman et al., 2019).

    Despite the short period of sampling, which covered most of the wet season over the region across all the sites,the observed daily isotopic composition is similar to the observed monthly data, with more enriched precipitation at BR compared to sampling period, which covered most of the wet season across all sites, the observed daily isotopic composition is similar to the observed monthly data, with more enriched precipitation observed at BR than BH and RC. The intercept of LMWLs for BH and RC sites showed lower values relative to the values derived from monthly data. According to Kurita et al. (2009), this is an expected trend because, on the monthly scale, the sum of rainfall events inhibits the local processes (post-condensation processes) that tend to decrease intercept values.

    For this reason, the daily isotopic compositions (Fig. 4)are related to specific climatic features, acting over the region in different periods of the wet season. Contrary to those observed on a monthly timescale (Table 2), there was no amount effect observed on daily timescales (Table 3),and the most depleted values are not related to the larger precipitation events. During the wet season, both the radiation flux and temperatures are normally higher, which favor increased convective activity and the occurrence of largescale SACZ events (Reboita et al., 2010).

    Fig. 5. Time series between δ18O and NOAA meteorological data for all sites: (a) Outgoing Longwave Radiation (OLR, W m-2); (b) Precipitable Water (kg m-2); (c) Omega (Vertical velocity)at 500 hPa. Arrows indicate the limit of values associated with convective activity (< 240 W m-2 for OLR and negative for Omega).

    Thereby, the differing displacements of daily isotopic variations (Fig. 4) are explained by convective activity(OLR < 240 W m-2and negative values of Omega) and moisture contribution (HYSPLIT trajectories) (Fig. 6). The depleted δ18O values (Fig. 5) in BH and BR were coupled with the lower OLR values between November-December(avg: 198 and 195 W m-2, respectively), between February-March at BR (avg: 204 W m-2), and between January (avg:188 W m-2) and March at RC (avg: 204 W m-2) (Fig. 5a).Omega values (Fig. 5b) have similar distributions with OLR values, confirming the convective activity. This is evidenced by higher correlations between δ18O-OLR (r > 0.55; p <0.05) and δ18O-Omega (r > 0.61; p < 0 .05) (Table 3).Depleted δ18O values were observed during negative peaks of omega values, representing a rising motion (Fig. 5b) and for low OLR values that occurred during the same period(Fig. 5a), while positive omega values were observed with enriched δ18O values during September and early November(Fig. 5b), along with their associated higher OLR values (>240 W m-2) (Fig. 5a).

    Large amounts of moisture in the atmosphere [high precipitable water values (PW)] were also observed during the periods characterized by depleted δ18O values and high convective activity (Fig. 5c). The first incursions of more depleted values in RC (Fig. 5c) were not observed before December (PW average between September-December:32.18 kg m-2) when an increase of the convective activity and PW occurred (PW avg during δ18O depleted period in Fig. 5c: 41.59 kg m-2). On the other hand, the first depleted incursions at BR and BH were observed in early November(Fig. 5c), associated with an early increase in PW and convective activity (Fig. 5). The observed differences indicate that regional climatic features are modulated in different ways at each site.

    The temporal variability on moisture availability also impacts d-excess values, as indicated by the differences observed in the daily LMWL intercepts and the Kruskal-Wallis test. Thus, despite the constant increase in the availability of moisture in the atmosphere from September to December 2016, large variations in the PW values were observed for all locations: BH (19.23 kg m-2-41.49 kg m-2), BR(30.25 kg m-2-43.14 kg m-2), and RC (19.68 kg m-2-40.77 kg m-2). These large PW variations led to a greater variability for the observed d-excess values for this period(Fig. 4g). In contrast, from January to March, the PW variability was lower for BH (31.23 kg m-2-39.54 kg m-2) and RC(30.49 kg m-2-44.94 kg m-2), while BR was similar(30.49 kg m-2- 44.94 kg m-2), which led to lower observed variability in d-excess values (Fig. 4g).

    Thereby, the large variability observed for PW and consequent d-excess values can be associated with different moisture sources and transport conditions that interact with atmospheric dynamics and regional circulation features, such as the LLJ (moisture transport from Amazon to central-south of Brazil) and the SASH (moisture transport from the Atlantic Ocean to central-south of Brazil). The trajectory maps display this difference between the sites with four moisture sources: land, ocean, and two with pathway over Amazon Forest (influenced by reevaporated moisture), originating from Land (named Land-Amazon; Fig. 6) or Ocean (Ocean-Amazon). The trajectories influenced by Amazon are accounted in the same percentage. Thus,. the trajectory contributions to the BH are 17% from the land and 83% from the ocean (9% influenced by Amazon Basin), for BR, 32%from the land and 68% from the ocean (22% from the Amazon Basin), and for RC 40% from the land and 60% from the ocean (21% from the Amazon basin) (Fig. 6).

    Amazon moisture flux (Land-Amazon and Ocean-Amazon in Fig. 6a) traveled short distances to reach BR and,once associated with the SASH circulation (Ocean trajectories in Fig. 6a), they favor the moisture passage over BR (Fig.6a), resulting in higher variations of PW and elevated dexcess values (avg: +14.40‰). On the other hand, the moisture from the Amazon (Land-Amazon and Ocean-Amazon;Fig. 6c) that reaches the RC site does experience larger transit distances. In addition, the site’s association with coldfrontal activities (Ocean trajectories in the southern portion of map; Fig. 6c) produces a mixture of different moisture sources, resulting in lower d-excess values (avg: +13.09‰).Most of the moisture resulting in precipitation over BH originated from the Atlantic Ocean (Fig. 6b) and had a shorter pathway compared to the moisture arriving directly from the Atlantic Ocean to BR and RC; this is thought to have resulted in their lower d-excess values (avg: +11.55‰).

    In addition, the higher correlations (r ≥ 0.50) observed between δ18O at the daily timescale and the regional parameters (Rainfall and Mix depth along HYSPLIT trajectories,PW, Omega, and OLR) (Table 3 and Fig. 5) confirm the importance of regional factors on governing isotope variations. They also reveal that seasonal effect related to the amount effect at a monthly timescale (r ≥ 0.60 between δ18O-monthly precipitation) is explained by a combination of moisture sources, recirculation processes, and convective activity, better observed at a daily timescale.

    Despite the relationship between δ18O values and convective activity in RC, no statistical correlations were observed between δ18O and OLR, Omega, and rainfall type by the ERA-interim. Thus, PW and Rainfall-HYSPLIT exhibited the best correlations, indicating that moisture availability is the key climate driver for RC, characterized by distinct moisture contributions and atmospheric systems, confirming the results demonstrated in previous studies (dos Santos et al.,2019a, c).

    The effect of regional climatic characteristics was more evident in BR through correlations between δ18O and regional meteorological parameters (Omega, OLR, PW, Rainfall-HYSPLIT, and Mix depth-HYSPLIT). The daily analysis confirmed the importance of the continental effect on a monthly scale, represented by a group of climatic features characterized by the interaction with convective activity, moisture transport, moisture source, and atmospheric systems(mainly the SACZ); thus, it is not possible to determine a single climatic driver for BR.

    For BH, convective activity and moisture transport are important climatic features, confirmed by the strong correlation between Omega and Rainfall-HYSPLIT. The explanation for convective activity (movement of winds by atmospheric systems) is associated with SACZ and CF incursions, complemented by a weak correlation with δ18O and OLR. The importance of air movement (convection or subsidence) during the formation of precipitation systems can be explained by the strong correlation between d-excess and Cv-Era in BH. These findings are coherent with previous assessments done by Risi et al. (2010), in which d-excess values were associated with processes of subsidence in the transition zones of convective-stratiform precipitation. However,few discussions have involved d-excess values and rainfall type due to the inherent difficulty in associating meteorological parameters with d-excess variability.

    Fig. 6. Times series between δ18O and NOAA meteorological data for all sites: (a) Outgoing Longwave Radiation (OLR, W m-2); (b) Precipitable Water (kg m-2); (c) Omega (vertical velocity) at 500 hPa. Arrows indicate the limit of values to convective activity (< 240 W m-2 for OLR and negative for Omega).

    One possible explanation for weak correlations between daily δ18O-δ2H and Cv-St precipitation is the sampling frequency and monitoring design since the classifications by Aggarwal et al. (2016) were specific to monthly data. The monthly δ18O may be biased because the precipitation amount (P) also tends to increase with higher stratiform rain fraction, as mentioned by Konecky et al. (2019), such that the roles of moisture convergence and microphysical processes like rain re-evaporation (which also increases with P)remain hidden variables, mainly on monthly scales. The biased P was represented by similar distributions between daily δ18O-local precipitation (Fig. 4) and δ18O-Cv precipitation (Fig. S3), which also occurred with their correlations(Table 3). In addition, the daily data from RC was excluded from the convective/stratiform classification data using Global Precipitation Measurement (GPM) by Munksgaard et al. (2019) due to weak correlation.

    Moreover, SACZ and CF are responsible for organizing the convective activity over Brazil, generating precipitation systems formed by convective and stratiform precipitation(Machado and Rossow, 1993; Romatschke and Houze,2013). These precipitation systems are estimated to be 45%convective and 46% stratiform in South America, noting that rainfall type is also difficult to partition (Romatschke and Houze, 2013).

    Therefore, a combination of convective activity and moisture availability to the formation of precipitation systems agrees with other studies (Kurita et al., 2011; Kurita, 2013;Torri et al., 2017; Lacour et al., 2018) that demonstrate the major importance of climatic mechanisms, such as convection in relation to cloud microphysics (Aggarwal et al., 2016).

    Our results demonstrate a relevant application of the role of convective activity on the isotopic variability of daily precipitation as a response to the interaction between moisture-convection, advancing the understanding of isotopic depletion and its poor correlation with the amount effect. However, it also highlights the strong relationship with δ18O-OLR, as shown by Vuille and Werner (2005). Furthermore, it complements the assessment conducted by Kurita et al. (2011) and Lacour et al. (2018), whose focus was an assessment of the convection and isotopic composition of vapor and its response to precipitation. The spatial distribution of convection in locations over the continent differs from the regions over the tropical oceans (Torri et al., 2017)and those over a moisture convergence zone. The latter features intense convective activity, heavy rainfall, and large spatial distribution, such as the SACZ, which is generally used to explain the depleted isotopic composition of precipitation over the Americas (Gastmans et al., 2017; dos Santos et al.,2019c; Santos et al., 2021).

    6. Conclusion

    Sampling precipitation assessments of different timescales indicate how the different climatic features affect the isotopic composition of precipitation. While it was possible to visualize the isotopic variability associated with classic effects (seasonal, precipitation amount, and continentality)at the monthly scale, the daily isotopic composition revealed stronger relations to specific climatic features, such as convective activity and moisture transport.

    Therefore, the variability observed at the daily scale was very helpful to understand how the available moisture interacts with the convective activity resulting in depleted δ18O incursions, which are observed frequently during the wet season at different days at each site, and how the superposition of several events associated with these processes can result in monthly-scale depleted precipitation, leading to similar monthly seasonal variations for all sites, as well as for ENSO effects.

    The observed differences in d-excess values, even on a monthly and daily basis, confirm the importance of recirculation processes that are governed by different moisture sources (the Atlantic Ocean and Amazon Forest) and pathways (modulated by the dynamics of atmospheric systems),resulting in depleted δ18O values and large variability of dexcess values.

    Thus, to quantify and improve the knowledge of tropical climate dynamics in the central-southeastern region of Brazil, there is a need to investigate the climatic (mechanisms) related to convection and the processes responsible for rainfall type (convective vs. stratiform) formation with the most robust meteorological data, such as radar or satellite products with enhanced resolutions, and intra-event collection of isotopic composition precipitation to investigate events related to the SACZ and CF.

    Our results provide insights for new isotopic interpretation within the central-southeastern region of Brazil. Finally,our study may help to explain the hydrometeorological processes in a tropical region that is highly dependent on rainfall for hydroelectric power generation and agricultural production.

    Acknowledgements. This work was funded by grants from the S?o Paulo Research Foundation (FAPESP) under Process 2018/06666-4, and by the International Atomic Energy Agency grant BRA-17984 under the initiative CRP-F31004 “Stable isotopes in precipitation and paleoclimatic archives in tropical areas to improve regional hydrological and climatic impact models”, and BRA-23531 under the initiative CRP-F31006 “Isotope Variability of Rain for Assessing Climate Change Impacts”. VS thanks FAPESP for the scholarship provided under the Process 2019/03467-3. We would also like to thank Ling ZHOU and Xiayu YUAN , editors of the Advances in Atmospheric Sciences and two anonymous reviewers for their comments that have improved the initial manuscript.

    Electronic supplementary material: Supplementary material is available in the online version of this article at https//doi.org/10.1007/s00376-022-1367-0.

    久久久久久国产a免费观看| 身体一侧抽搐| 亚洲第一区二区三区不卡| 肉色欧美久久久久久久蜜桃 | 午夜免费观看性视频| 夜夜爽夜夜爽视频| 成人午夜精彩视频在线观看| a级一级毛片免费在线观看| 高清视频免费观看一区二区 | 免费看美女性在线毛片视频| 综合色丁香网| 免费观看性生交大片5| 少妇裸体淫交视频免费看高清| 成人午夜高清在线视频| 亚洲国产精品成人综合色| 深爱激情五月婷婷| 国产一区亚洲一区在线观看| 亚洲av男天堂| 男女下面进入的视频免费午夜| 亚洲欧美日韩无卡精品| 日韩伦理黄色片| 亚洲av中文av极速乱| 国产精品一区二区在线观看99 | 国产精品一区二区性色av| 免费看光身美女| 国产毛片a区久久久久| 国产免费一级a男人的天堂| 国产真实伦视频高清在线观看| 99热网站在线观看| 日本午夜av视频| 国产极品天堂在线| 美女主播在线视频| 婷婷色综合大香蕉| 久久久久久伊人网av| 黑人高潮一二区| 亚洲熟妇中文字幕五十中出| 国产欧美另类精品又又久久亚洲欧美| 国国产精品蜜臀av免费| 最近视频中文字幕2019在线8| 国产黄a三级三级三级人| 我要看日韩黄色一级片| 国产精品精品国产色婷婷| 啦啦啦中文免费视频观看日本| 一级毛片我不卡| 国产中年淑女户外野战色| 18禁动态无遮挡网站| 又爽又黄a免费视频| 波多野结衣巨乳人妻| 日韩,欧美,国产一区二区三区| 国产老妇女一区| 日韩av在线免费看完整版不卡| 日本午夜av视频| 精品午夜福利在线看| 午夜亚洲福利在线播放| 亚洲国产最新在线播放| 亚洲国产色片| 国产精品久久久久久久久免| 又黄又爽又刺激的免费视频.| 中文字幕av在线有码专区| 久久久久久久久久久免费av| 亚洲av日韩在线播放| 只有这里有精品99| 日韩制服骚丝袜av| 国产精品三级大全| 亚洲欧洲日产国产| 免费大片18禁| 插阴视频在线观看视频| 熟女人妻精品中文字幕| 欧美 日韩 精品 国产| 国产在线男女| 亚洲经典国产精华液单| 神马国产精品三级电影在线观看| 少妇被粗大猛烈的视频| 又粗又硬又长又爽又黄的视频| 综合色丁香网| 亚洲国产色片| 久久久久久久大尺度免费视频| 极品少妇高潮喷水抽搐| 观看免费一级毛片| 蜜桃久久精品国产亚洲av| 在线播放无遮挡| 一个人看视频在线观看www免费| 亚洲在久久综合| 国产一区二区三区综合在线观看 | eeuss影院久久| 只有这里有精品99| av国产免费在线观看| 欧美精品国产亚洲| 成人午夜高清在线视频| 国内精品宾馆在线| 成人性生交大片免费视频hd| 最近的中文字幕免费完整| 亚洲美女搞黄在线观看| 中国美白少妇内射xxxbb| 三级国产精品欧美在线观看| 日韩不卡一区二区三区视频在线| 美女xxoo啪啪120秒动态图| 男人爽女人下面视频在线观看| 麻豆乱淫一区二区| 天堂网av新在线| 成人欧美大片| 国产伦精品一区二区三区四那| 中文字幕av在线有码专区| 欧美不卡视频在线免费观看| 亚洲精品aⅴ在线观看| 国产精品av视频在线免费观看| 激情 狠狠 欧美| 男人和女人高潮做爰伦理| 女人被狂操c到高潮| 精品一区二区三卡| 国产白丝娇喘喷水9色精品| 日日啪夜夜撸| 国产欧美另类精品又又久久亚洲欧美| 国产精品久久视频播放| 亚洲精品456在线播放app| 亚洲国产日韩欧美精品在线观看| 男女啪啪激烈高潮av片| 精品久久国产蜜桃| 午夜福利视频精品| 韩国av在线不卡| 国产色婷婷99| 亚洲色图av天堂| 3wmmmm亚洲av在线观看| or卡值多少钱| 久久久久精品性色| 久久久久久久大尺度免费视频| 97超碰精品成人国产| 国产老妇女一区| 美女大奶头视频| 久久久精品94久久精品| 国产精品嫩草影院av在线观看| 狠狠精品人妻久久久久久综合| 亚洲在久久综合| 亚洲精品视频女| 久久久国产一区二区| 亚洲婷婷狠狠爱综合网| 人妻系列 视频| 国产日韩欧美在线精品| 国产精品一区二区在线观看99 | 免费播放大片免费观看视频在线观看| 国产亚洲最大av| 精品欧美国产一区二区三| 中文乱码字字幕精品一区二区三区 | 国产av国产精品国产| 午夜精品在线福利| 久久人人爽人人片av| 免费看日本二区| 最近最新中文字幕大全电影3| 国产精品一区二区在线观看99 | 一级毛片久久久久久久久女| 国产精品久久久久久精品电影小说 | 成人毛片a级毛片在线播放| 日本色播在线视频| 中国国产av一级| 国产精品美女特级片免费视频播放器| 日韩电影二区| 永久免费av网站大全| 美女国产视频在线观看| 久久99热这里只频精品6学生| 水蜜桃什么品种好| 丝袜喷水一区| 白带黄色成豆腐渣| 哪个播放器可以免费观看大片| 狂野欧美激情性xxxx在线观看| 亚洲精品日本国产第一区| 亚洲国产精品国产精品| 国产精品1区2区在线观看.| 少妇被粗大猛烈的视频| 搡老妇女老女人老熟妇| 麻豆av噜噜一区二区三区| 国产伦在线观看视频一区| av网站免费在线观看视频 | 国产精品久久久久久精品电影小说 | 日韩av在线免费看完整版不卡| 高清av免费在线| 一个人观看的视频www高清免费观看| 91精品国产九色| 日韩,欧美,国产一区二区三区| 亚洲精品,欧美精品| 99九九线精品视频在线观看视频| 秋霞在线观看毛片| 精品久久久久久久人妻蜜臀av| freevideosex欧美| 深爱激情五月婷婷| 晚上一个人看的免费电影| 少妇高潮的动态图| 少妇丰满av| 97热精品久久久久久| 久久久久久久久久久丰满| 欧美日韩综合久久久久久| 国产成人91sexporn| 一区二区三区乱码不卡18| 日韩视频在线欧美| 久久热精品热| 深爱激情五月婷婷| 久久99精品国语久久久| 汤姆久久久久久久影院中文字幕 | 日本av手机在线免费观看| 久久久久网色| 国产精品久久久久久久久免| 精品久久久噜噜| 欧美激情久久久久久爽电影| 美女大奶头视频| 国产av不卡久久| 少妇猛男粗大的猛烈进出视频 | 最后的刺客免费高清国语| 在现免费观看毛片| 精品久久久久久成人av| 搞女人的毛片| 男人舔女人下体高潮全视频| 毛片一级片免费看久久久久| 亚洲av在线观看美女高潮| 久久久久久久久久人人人人人人| 国内精品美女久久久久久| 国产不卡一卡二| 床上黄色一级片| 久久久久九九精品影院| av在线老鸭窝| 边亲边吃奶的免费视频| 中国国产av一级| 久久久久久久国产电影| 少妇熟女aⅴ在线视频| 美女国产视频在线观看| 中文乱码字字幕精品一区二区三区 | 日本免费在线观看一区| 欧美精品国产亚洲| 在线免费观看不下载黄p国产| 国产在视频线在精品| 欧美激情国产日韩精品一区| 哪个播放器可以免费观看大片| 99热网站在线观看| 午夜福利视频1000在线观看| 国产淫语在线视频| 久久99蜜桃精品久久| 69人妻影院| 日日啪夜夜爽| 国产免费一级a男人的天堂| 综合色av麻豆| 极品教师在线视频| 精品国产一区二区三区久久久樱花 | 熟女人妻精品中文字幕| 18+在线观看网站| 亚洲av中文av极速乱| 国产精品日韩av在线免费观看| 国产老妇伦熟女老妇高清| 精品国内亚洲2022精品成人| 国内精品美女久久久久久| 国产av不卡久久| 午夜爱爱视频在线播放| 午夜免费观看性视频| 国内精品美女久久久久久| 国产成人精品婷婷| 男女视频在线观看网站免费| 特大巨黑吊av在线直播| 欧美人与善性xxx| 国产精品福利在线免费观看| 欧美一区二区亚洲| 成人二区视频| 国内精品一区二区在线观看| 婷婷六月久久综合丁香| 欧美激情在线99| 亚洲av电影不卡..在线观看| 一个人看的www免费观看视频| 亚洲欧美成人综合另类久久久| 在线观看免费高清a一片| 亚洲av在线观看美女高潮| 国产老妇女一区| 欧美bdsm另类| 精品久久久精品久久久| 国产日韩欧美在线精品| 大话2 男鬼变身卡| 在线免费观看的www视频| 午夜免费激情av| 欧美 日韩 精品 国产| 国产精品99久久久久久久久| 亚洲av男天堂| 少妇熟女欧美另类| 日韩在线高清观看一区二区三区| 国产精品无大码| 国产黄片视频在线免费观看| 亚洲精品一区蜜桃| 国产成人精品久久久久久| 一级毛片电影观看| 99久国产av精品| 99视频精品全部免费 在线| 日韩三级伦理在线观看| 亚洲成人久久爱视频| 看黄色毛片网站| 亚洲精品,欧美精品| 嫩草影院入口| 联通29元200g的流量卡| 日本与韩国留学比较| 美女内射精品一级片tv| 日韩电影二区| 五月伊人婷婷丁香| 美女cb高潮喷水在线观看| 人妻一区二区av| 免费不卡的大黄色大毛片视频在线观看 | 国产日韩欧美在线精品| 午夜日本视频在线| 69人妻影院| 免费黄频网站在线观看国产| 国产精品一区二区性色av| 国产视频首页在线观看| 国产片特级美女逼逼视频| 久久久久免费精品人妻一区二区| 免费播放大片免费观看视频在线观看| 午夜免费激情av| 精品久久久久久久人妻蜜臀av| eeuss影院久久| 好男人视频免费观看在线| 中文字幕久久专区| 免费看a级黄色片| 美女内射精品一级片tv| 青春草视频在线免费观看| 五月伊人婷婷丁香| 国产淫片久久久久久久久| 尾随美女入室| 色综合亚洲欧美另类图片| 欧美xxxx性猛交bbbb| 天天一区二区日本电影三级| 亚洲18禁久久av| 在线 av 中文字幕| 免费看光身美女| 床上黄色一级片| 精品久久久精品久久久| 国产精品麻豆人妻色哟哟久久 | 好男人视频免费观看在线| 免费观看精品视频网站| eeuss影院久久| 人妻少妇偷人精品九色| 国产毛片a区久久久久| 亚洲伊人久久精品综合| av女优亚洲男人天堂| 大又大粗又爽又黄少妇毛片口| 国产男女超爽视频在线观看| 又粗又硬又长又爽又黄的视频| 亚洲精品成人av观看孕妇| 成人一区二区视频在线观看| 夜夜看夜夜爽夜夜摸| 精品久久久久久电影网| 观看美女的网站| av线在线观看网站| 亚洲色图av天堂| 国产一区二区三区av在线| 日日干狠狠操夜夜爽| 亚洲av在线观看美女高潮| 天天躁夜夜躁狠狠久久av| 成人漫画全彩无遮挡| 九九爱精品视频在线观看| 三级国产精品欧美在线观看| 国产一区二区亚洲精品在线观看| 中文乱码字字幕精品一区二区三区 | 国产午夜精品一二区理论片| 国产精品综合久久久久久久免费| 免费观看无遮挡的男女| 久久久精品免费免费高清| 日本午夜av视频| 国产精品一区www在线观看| 国产黄片美女视频| 激情 狠狠 欧美| 免费少妇av软件| 成人国产麻豆网| 免费黄色在线免费观看| 国产成人一区二区在线| 最近最新中文字幕免费大全7| 成人鲁丝片一二三区免费| 三级经典国产精品| 色哟哟·www| 欧美极品一区二区三区四区| 日韩人妻高清精品专区| 少妇猛男粗大的猛烈进出视频 | av在线亚洲专区| av又黄又爽大尺度在线免费看| 欧美 日韩 精品 国产| 国产黄色免费在线视频| 亚洲最大成人手机在线| 人人妻人人澡人人爽人人夜夜 | 亚洲国产精品国产精品| 日韩中字成人| 日韩一区二区三区影片| 麻豆久久精品国产亚洲av| 大片免费播放器 马上看| 久久精品国产亚洲av涩爱| 日本猛色少妇xxxxx猛交久久| videos熟女内射| 美女主播在线视频| 日本黄色片子视频| 国产爱豆传媒在线观看| 国模一区二区三区四区视频| 亚洲,欧美,日韩| av一本久久久久| 青春草视频在线免费观看| 亚洲综合色惰| 国语对白做爰xxxⅹ性视频网站| 一级毛片黄色毛片免费观看视频| 男女边摸边吃奶| 欧美最新免费一区二区三区| 亚洲第一区二区三区不卡| 久久久成人免费电影| 毛片一级片免费看久久久久| 精品久久久久久成人av| 精品久久久久久久人妻蜜臀av| 91久久精品电影网| 亚洲最大成人手机在线| 日本黄大片高清| 好男人在线观看高清免费视频| 身体一侧抽搐| 国产成人aa在线观看| 精品人妻偷拍中文字幕| 大香蕉97超碰在线| 日本-黄色视频高清免费观看| 18禁在线无遮挡免费观看视频| 伊人久久精品亚洲午夜| 在线观看免费高清a一片| 女人久久www免费人成看片| 少妇猛男粗大的猛烈进出视频 | 精品一区二区免费观看| 最近最新中文字幕大全电影3| 国产精品国产三级国产av玫瑰| 久久精品国产亚洲av天美| 九九久久精品国产亚洲av麻豆| 狠狠精品人妻久久久久久综合| 五月伊人婷婷丁香| 日本三级黄在线观看| 免费av毛片视频| 老女人水多毛片| 亚洲av免费在线观看| 欧美一级a爱片免费观看看| 亚洲av成人av| 中文字幕av在线有码专区| 嫩草影院精品99| 在线播放无遮挡| 在线观看美女被高潮喷水网站| 久久久久久久久久成人| 一个人观看的视频www高清免费观看| 成人亚洲精品一区在线观看 | 亚洲av男天堂| 国产免费又黄又爽又色| 久久久精品94久久精品| 欧美激情国产日韩精品一区| 国产免费视频播放在线视频 | 国产精品美女特级片免费视频播放器| 2022亚洲国产成人精品| 欧美3d第一页| 亚洲欧洲国产日韩| 午夜久久久久精精品| av又黄又爽大尺度在线免费看| 联通29元200g的流量卡| 我要看日韩黄色一级片| 久久久久国产网址| 18禁裸乳无遮挡免费网站照片| 肉色欧美久久久久久久蜜桃 | 国产欧美日韩精品一区二区| 最近视频中文字幕2019在线8| 久久精品国产鲁丝片午夜精品| 看非洲黑人一级黄片| 国产在线一区二区三区精| 777米奇影视久久| a级一级毛片免费在线观看| 国产精品日韩av在线免费观看| 国产亚洲精品av在线| 日本爱情动作片www.在线观看| 成年版毛片免费区| 日本一二三区视频观看| 中文字幕av成人在线电影| 99热这里只有是精品在线观看| av免费观看日本| av线在线观看网站| 国产精品一区二区三区四区久久| 中文字幕av在线有码专区| 国产黄片视频在线免费观看| 欧美日韩综合久久久久久| 久久久色成人| 蜜臀久久99精品久久宅男| 国产v大片淫在线免费观看| 亚洲怡红院男人天堂| 国产日韩欧美在线精品| 成人特级av手机在线观看| 麻豆成人av视频| 91精品伊人久久大香线蕉| 日本免费在线观看一区| 97精品久久久久久久久久精品| 亚洲成人精品中文字幕电影| 国产久久久一区二区三区| 亚洲av免费在线观看| 免费看不卡的av| 91在线精品国自产拍蜜月| 欧美日韩综合久久久久久| 水蜜桃什么品种好| 亚洲成人中文字幕在线播放| 老师上课跳d突然被开到最大视频| 高清av免费在线| 久久久久免费精品人妻一区二区| 毛片一级片免费看久久久久| 欧美精品一区二区大全| 一级av片app| 精品一区二区免费观看| 精品久久久久久久人妻蜜臀av| 日本爱情动作片www.在线观看| 国产av码专区亚洲av| 久99久视频精品免费| 日韩欧美精品v在线| av线在线观看网站| 一级毛片电影观看| 国产成人精品福利久久| 亚洲精品国产av成人精品| 男女啪啪激烈高潮av片| 成年av动漫网址| 国产 亚洲一区二区三区 | 晚上一个人看的免费电影| 国产精品福利在线免费观看| av线在线观看网站| 日本猛色少妇xxxxx猛交久久| 一级爰片在线观看| 日韩三级伦理在线观看| 免费看不卡的av| 六月丁香七月| 国内少妇人妻偷人精品xxx网站| 欧美 日韩 精品 国产| 日韩大片免费观看网站| 国产精品人妻久久久久久| 国产精品久久久久久av不卡| 纵有疾风起免费观看全集完整版 | 男的添女的下面高潮视频| 亚洲成人一二三区av| 日本wwww免费看| 成人一区二区视频在线观看| 日本爱情动作片www.在线观看| 亚洲欧美精品自产自拍| 麻豆av噜噜一区二区三区| 菩萨蛮人人尽说江南好唐韦庄| 成年女人在线观看亚洲视频 | 亚洲av电影不卡..在线观看| 亚洲国产欧美在线一区| 国产伦在线观看视频一区| 欧美成人一区二区免费高清观看| 日韩人妻高清精品专区| 极品教师在线视频| 亚洲18禁久久av| 人妻一区二区av| 18禁裸乳无遮挡免费网站照片| 久久99精品国语久久久| 别揉我奶头 嗯啊视频| 国产精品三级大全| 欧美激情久久久久久爽电影| 一本一本综合久久| 深夜a级毛片| 男女那种视频在线观看| 欧美一区二区亚洲| 一级毛片黄色毛片免费观看视频| 又爽又黄无遮挡网站| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 欧美日本视频| 精品一区在线观看国产| 一级毛片aaaaaa免费看小| 我的女老师完整版在线观看| 欧美成人午夜免费资源| 2021天堂中文幕一二区在线观| 毛片女人毛片| 精品少妇黑人巨大在线播放| xxx大片免费视频| 亚洲av在线观看美女高潮| 国产午夜精品久久久久久一区二区三区| 亚洲四区av| 欧美最新免费一区二区三区| 国产高清有码在线观看视频| 日本猛色少妇xxxxx猛交久久| 亚洲国产精品国产精品| 亚洲熟女精品中文字幕| 午夜激情福利司机影院| ponron亚洲| 国产亚洲精品久久久com| 精品国产露脸久久av麻豆 | 精品国产露脸久久av麻豆 | ponron亚洲| 99热这里只有精品一区| 国国产精品蜜臀av免费| 国产精品1区2区在线观看.| 亚洲精品第二区| 亚洲美女视频黄频| 午夜免费男女啪啪视频观看| 六月丁香七月| 日本与韩国留学比较| 日本一本二区三区精品| 女人十人毛片免费观看3o分钟| 大话2 男鬼变身卡| 欧美极品一区二区三区四区| 国产亚洲91精品色在线| 不卡视频在线观看欧美| 99久国产av精品国产电影| 中文欧美无线码| 大又大粗又爽又黄少妇毛片口| 嫩草影院新地址| 国产高潮美女av| 亚洲av电影不卡..在线观看| 亚洲精品国产成人久久av| 99视频精品全部免费 在线| 国产亚洲一区二区精品| 国产av不卡久久| 免费大片黄手机在线观看| 亚洲av免费高清在线观看| 日韩中字成人| 99久久精品一区二区三区| 久久午夜福利片| 亚洲国产精品sss在线观看| 亚洲三级黄色毛片| 在线观看免费高清a一片| 国产亚洲5aaaaa淫片| 久久久久久久亚洲中文字幕| 免费观看性生交大片5| 亚洲av中文av极速乱| 国产日韩欧美在线精品| 国产精品久久久久久av不卡| 日本与韩国留学比较|