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

    Thermo-hydro-poro-mechanical responses of a reservoir-induced landslide tracked by high-resolution fiber optic sensing nerves

    2024-03-25 11:06:36XioYeHongHuZhuGngChengHuFuPeiBinShiLuShentoAlessnroPsuto

    Xio Ye,Hong-Hu Zhu,*,Gng Cheng,b,Hu-Fu Pei,Bin Shi,Lu Shento,Alessnro Psuto

    a School of Earth Sciences and Engineering, Nanjing University, Nanjing, 210023, China

    b School of Computer Science, North China Institute of Science and Technology (National Safety Training Center of Coal Mines), Beijing,101601, China

    c State Key Lab of Coastal and Offshore Engineering, Department of Geotechnical Engineering, Dalian University of Technology, Dalian,116024, China

    d National Research Council-Research Institute for Geo-Hydrological Protection (CNR-IRPI), Padova, 35127, Italy

    e Department of Information Engineering, University of Padova, Padova, 35131, Italy

    Keywords: Reservoir landslide Thermo-hydro-poro-mechanical response Ultra-weak fiber bragg grating (UWFBG) subsurface evolution Engineering geological interface Geotechnical monitoring

    ABSTRACT Thermo-poro-mechanical responses along sliding zone/surface have been extensively studied.However,it has not been recognized that the potential contribution of other crucial engineering geological interfaces beyond the slip surface to progressive failure.Here,we aim to investigate the subsurface multiphysics of reservoir landslides under two extreme hydrologic conditions (i.e.wet and dry),particularly within sliding masses.Based on ultra-weak fiber Bragg grating(UWFBG)technology,we employ specialpurpose fiber optic sensing cables that can be implanted into boreholes as“nerves of the Earth”to collect data on soil temperature,water content,pore water pressure,and strain.The Xinpu landslide in the middle reach of the Three Gorges Reservoir Area in China was selected as a case study to establish a paradigm for in situ thermo-hydro-poro-mechanical monitoring.These UWFBG-based sensing cables were vertically buried in a 31 m-deep borehole at the foot of the landslide,with a resolution of 1 m except for the pressure sensor.We reported field measurements covering the period 2021 and 2022 and produced the spatiotemporal profiles throughout the borehole.Results show that wet years are more likely to motivate landslide motions than dry years.The annual thermally active layer of the landslide has a critical depth of roughly 9 m and might move downward in warmer years.The dynamic groundwater table is located at depths of 9-15 m,where the peaked strain undergoes a periodical response of leap and withdrawal to annual hydrometeorological cycles.These interface behaviors may support the interpretation of the contribution of reservoir regulation to slope stability,allowing us to correlate them to local damage events and potential global destabilization.This paper also offers a natural framework for interpreting thermo-hydro-poro-mechanical signatures from creeping reservoir bank slopes,which may form the basis for a landslide monitoring and early warning system.

    1.Introduction

    The construction and operation of dams inevitably create reservoirs that can affect the stability of slopes distributed along riverbanks (Schuster,1979;Paronuzzi et al.,2013;Chen et al.,2022).Large-scale catastrophic reservoir landslides are a potential threat to local communities and infrastructures.The hazards are mainly related to the large-volume rock and soil debris (i.e.value of over 1 ×106m3) and potential high-speed motion.Moreover,it can be amplified by the occurrence of extreme rainfall events that are showing an increasing frequency in recent years due to climatic changes (Alonso et al.,2016;Jones et al.,2021;Ye et al.,2022).Unfortunately,these geohazard events might not be expected to lessen over time,for instance,the first half of 2022 was the deadliest on record for landslides (Westra et al.,2014;Cendrero et al.,2020;Ozturk et al.,2022;Ma and Yuan,2023).

    The Vajont landslide,which occurred on October 9,1963,was an iconic reservoir landslide event that provoked a worldwide sensation and served as a wake-up call for hazard monitoring and mitigation in the reservoir area(Müller,1964;Veveakis et al.,2007;Paronuzzi et al.,2013).For this,a large number of researchers have made significant efforts and contributions to this issue,such as the Qiangjianping landslide in Zigui,China,2003 (Wang et al.,2004),San Juan de Grijalva slide in Mexico,2007 (Alcántara-Ayala and Domínguez-Morales,2008),and Outang landslide in Fengjie,China (Yin et al.,2016).The kinematics of these events are quite different,but they have been widely classified into flow-like earth motion and sliding motion (Hungr et al.,2014).The flow-like rheology is mainly related to the hydrodynamic and rock/soil properties,which is supported by power-law viscoplastic flow theory(Iverson et al.,1997;Cascini et al.,2010;Di Maio et al.,2020;Hu et al.,2020).Instead,for creeping landslides,it is usually focused on the location and evolutionary characteristics of potential sliding surfaces (Veveakis et al.,2007;Ye et al.,2022).The impacts of the cyclic filling and drawdown operation of the reservoir are twofold:triggering a seepage process,and/or inducing a stabilizing load at the submerged slide toe (Costa and Schuster,1988;Iqbal et al.,2017).Slope stability can thus be time-dependently enhanced or worsened.Hydromechanical responses of the sliding surface (i.e.shear zone)to predisposing factors often help in understanding the landslide deformation mechanism (Finnegan et al.,2021).The temperature changes in the shear band can be linked to the confined aquifer and landslide deformation (Veveakis et al.,2007;Seguí and Veveakis,2021),which can be further used to evaluate the stability of deep-seated landslide and predict landslide displacement.Nevertheless,there is still a gap in understanding the thermo-hydro-poro-mechanical behaviors inside sliding masses of an activated slow-moving reservoir landslide.

    One of the major challenges in landslide hydrogeology is to identify the location and characterize the multi-physics of crucial underground interfaces that can evolve into a sliding surface(Zhu,2023).Fiber optics have been widely employed to obtain accurate temperature and strain data subjected to external forces or temperature variations (Soga and Luo,2018;Shi et al.,2021;Xu et al.,2022a).However,measuring hydraulically relevant variables is challenging but significant,as water is the main mover of earth surface process (Palmer,2017).More recently,fiber-optic multiphysical sensing has reached a new level in groundwater monitoring,which permits us to infer subsurface hydraulic flow paths by utilizing the thermophysical effects of the material (Selker et al.,2006;Zhang and Xue,2019;Sun et al.,2022;Ye et al.,2022).The feasibility to obtain accurate soil-water content measurements using the active-heated fiber-optic (AHFO) temperature sensing method through laboratory tests is demonstrated (Sayde et al.,2010).Encouragingly,the AHFO technique has successfully applied to filed-scale monitoring of hydrology and heat transport in porous media (Read et al.,2013;Sun et al.,2021).Fiber optics has demonstrated their potential in monitoring subsurface multiphysical behavior at high resolutions (Hasler,2022;Sidder,2022).Nonetheless,there are few reports on field landslide monitoring of thermo-hydro-mechanical responses of sliding masses.Actually,this particular feature is crucial in understanding the seepage processes and stress states that affected the slope stability during the filling-drawdown cycles of the reservoir(Paronuzzi et al.,2013).

    In this paper,we present the results of thermo-hydro-poromechanical responses tracked by ultra-weak fiber Bragg grating(UWFBG) array deployed at the toe of the Xinpu landslide,in the middle reach of the Three Gorges Reservoir Area.These specialpurpose UWFBG-based sensing cables,acting as a “fiber optic nerve system(FONS)”of the ground(Ye et al.,2022),were vertically buried in a 31 m-deep borehole near the reservoir shoreline to collect data on soil temperature,water content,pore water pressure,and strain.The field-scale multi-physical measurements almost covering the years 2021-2022 display the spatiotemporal profiles in wet and dry hydrologic extremes.Thermo-hydro-poromechanical representations with an unprecedented level of detail will contribute to providing new insights of other crucial engineering geological interfaces besides the slip surface to infer flow fluid migration and stress relief within the sliding masses.

    2.Studied landslide

    2.1.Site description

    The study area is situated at Xinpu Village,Fengjie County,Chongqing,in the middle reach of the Three Gorges Reservoir Area in China (Fig.1).The Xinpu landslide,almost encompassing the whole village,is located on the right bank of the mainstream of the Yangtze River.The Shengjibao landslide was monitored 15 years ago,which is the predecessor of the current Xinpu landslide(i.e.the lower part of the landslide).In October 2006,the activation of the prior Shengjibao landslide was observed when the Three Gorges Reservoir was impounded to 156 m,resulting in extensive deformation caused by the drawdown of the reservoir water level(RWL)and seasonal rainfall (Jiang et al.,2020).In response,16 on-site surface displacement-monitoring stations were installed over the hillslope in March 2007.Since then,localized accelerateddeformation or even acute collapse occurred almost every year within and beyond the original Shengjibao area,leading to several relocations of the residents.Consequently,the boundary was redefined and renamed as the Xinpu landslide.Regrettably,the possibility of additional resettlements cannot be dismissed,as localized large deformations and devastating failure within the landslide have become increasingly severe in recent years than ever.

    Fig.1.Location of the study area:(a)Digital elevation model(DEM)(data from https://srtm.csi.cgiar.org/srtmdata/)map of the Three Gorges Reservoir area,China;and(b)Location of the Xinpu landslide bounded by yellow curve,and the sliding towards nearly north-direction.

    Fig.3.Overview of the Xinpu landslide:(a)Aerial photo of the studied landslide that is zoned four blocks;and(b)Partially enlarged study site and instrumented borehole FOS3,the monitoring box is left on the surface for storing acquisition,communication and control equipment introduced residential electricity supply,rather than a solar panel.

    As shown in Figs.2 and 3,the Xinpu landslide is a large-scale complex landslide that ranges from 90 m above sea level (a.s.l.)(toe) to 705 m a.s.l.(crown).It covers an area of 1.94 km2,with a volume of 5.4×107m3and an average slope gradient of about 18°-28°.The main sliding orientation of the slide is almost 347°.The unstable slope has been divided into four blocks in accordance with the landslide style of activity,with the lower part(Blocks I and IV)being the most dominant in terms of landslide deformation.For this,the monitoring campaign has been focused on this zone.The upper and middle parts of the slide are relatively stable (Blocks II and III).Note that the slide toe is perennially submerged under the operating RWL (i.e.145-175 m a.s.l.),making it a typical wading landslide.

    Previous drilling exploration data reveal the geological structure of the Xinpu landslide.The slide mass,from top to bottom,mainly consists of gravelly silty clay(i.e.with a gravel content of 20%-40%),gravel with silty clay (i.e.,with a gravel content of 60%-80%),and crushed siltstone(Fig.7).The thicknesses of these layers vary from 3 m to 23.8 m,1.2 m-39.4 m,and 1 m-14 m,respectively.The closer to the slope foot,the thicker the sliding accumulation,reaching 80 m at certain locations.

    2.2.Downhole fiber-optic cable configurations

    To investigate the multi-physical evolution within the landslide masses,a bunch of specific-purpose fiber optic cables was installed in a 31 m-deep borehole FOS3 at the toe of the Xinpu landslide to measure soil temperature,water content,pore water pressure,and strain,respectively (Figs.2 and 3).Among them,temperature,water content,and strain sensors were encapsulated into cables at 1 m sampling intervals,forming separate sensing arrays (Fig.4).Whereas three piezometric pressure sensors were installed in a 2 m-deep sand filter at depths of 10 m,20 m,and 30 m down the borehole.The interrogator was in situ left in a protective box at the surface near the borehole,and the data were acquired remotely and in real time,typically with a sampling frequency of 1-60 min(Fig.4f and g).The initial reference measurements started on February 11,2021,which was approximately three weeks after backfilling allowing possible effects of compressional consolidation of the filling materials in the borehole to relax and therefore eliminating their effects on the results.

    Fig.4.Field installation of UWFBG sensing arrays:(a)Temperature sensing cable;(b)Moisture sensing cable;(c)Strain sensing cable;(d)Pore water pressure cell;(e)Instrumented borehole to be downwards placed and backfilled;(f) Monitoring box at the top of the borehole FOS3;and (g) In situ interrogator.

    3.Data and methods

    3.1.Principle of UWFBG

    Fiber optic Bragg grating (FBG) is one of the most mature fiber optic sensing technology that has been widely employed for geoengineering health monitoring (Kelam et al.,2016,2022).The principle of FBG sensing technology is to create a wavelength selective mirror,by periodically changing the refractive index of the fiber core(i.e.the grating)(Ho et al.,2006).External forces,and the consequent longitudinal strain or temperature variations applied to the grating determine a reproducible shift of the central reflected wavelength of the grating (Soga and Luo,2018).In general,the amount of wavelength shift of the central reflect wavelength encodes simultaneously the overall effects of strain and temperature.

    One of the main advantages of this technology is the multiplexing capability,with tens of sensors localized in the same fiber to implement an array of daisy-chained sensors.If each FBG of such an array has a distinct central wavelength,the entire array can be simultaneously interrogated by the so-called wavelength division multiplexing (WDM) technique,which allow addressing each FBG by its own central wavelength.However,it fails to provide a truly continuous description of strain or temperature changes because the number of sensing points(i.e.gratings)along an optical fiber is limited to some tens,depending on the bandwidth requirements of each FBG(Moyo et al.,2005).The limited number of FBGs allocable in a single fiber determines,in turn,a constraint on the spatial resolution,especially if the area to be monitored is wide.Although FBG sensing,in some cases,can provide useful and cost-effective information about global behavior,in some applications it can give only a coarse or partial representation of a specific area of the phenomena.

    The UWFBG refers to as a particular FBG technology implemented by gratings with reflectivity of less than 0.01% (-40 dB).Askins et al.(1992) developed a weak-reflection FBG with reflectivity less than 2%and linewidth of about 0.1 nm,and proposed its use as a unit to form a serialized,distributed and addressable sensor array.UWFBG technology has developed over decades to make it possible to apply in various engineering scenarios(Lindner et al.,2014;Yang et al.,2016).In one of the most effective implementations,a large-scale and dense UWFBG sensor array is made up of hundreds or thousands of identical-wavelength FBGs.Given the spectral overlapping of the gratings,the time-division multiplexing (TDM) technique is utilized to enable the simultaneous measurement of numerous gratings along a single fiber,as depicted in Fig.5.This greatly enhances the sensing capacity and offers a spatial resolution of less than 1.0 m that is well-suited to capture hydro-and geo-mechanical behavior of subsurface thin-layer interfaces (Liu et al.,2021;Ye et al.,2022).With the TDM interrogation technique,each grating of the fiber can be addressed from the time difference among the arrival of reflected light passing the entire grating sequence back.Given the vertical deployment of the fiber in the borehole,it allows to probe a specific depth in the borehole where the sensor is located.As specified above,four separated UWFBG arrays were installed,one for each of the measured parameters (i.e.temperature,water content,strain,and pore water pressure).

    Fig.5.Sensing principle of UWFBG array.

    By means of temperature compensation or hollow configuration between the optical fiber and sheath,strain and temperature are readily separated from each other(Zhu et al.,2017).It is noted that both the sheaths of temperature and moisture cable are highly thermally conductive,ensuring highly-sensitive temperature measurement (Fig.4a and b).The water content measurement is performed by an actively-heated UWFBG (AH-UWFBG) moisture sensing cable,which allows for estimation of the water content from the temperature dynamics of the heated cable according to a specific temperature model (Read et al.,2013;Sun et al.,2021;Ye et al.,2022).The strain cable jacket is designed with anchor sections at an interval of 1 m to improve coupling with the soils,and the inclusion of wire ropes enhances its suitability in harsh field conditions (Fig.4c).

    3.2.Data preparation

    3.2.1.Temperature,water content,pore water pressure,and strain measurements

    To obtain accurate measurements,these sensors installed in the borehole need to be calibrated in advance by means of laboratory tests.The temperature or strain coefficient is determined by taking the average of the results of repeated multiple calibrations (see Fig.S1).The calibration methods can be referred to as the previous works (e.g.Liu et al.,2021;Zhang et al.,2021a;Ye et al.,2022).

    To calibrate the temperature-related fiber optic cables (i.e.temperature and moisture),the cables were coiled and placed in a thermostat for gradual heating to obtain different levels of wavelength and temperature curves,which were then used to calculate temperature sensitivity coefficients.Further,a remaining section of the moisture sensing cable installed in the field borehole was buried in a laboratory-scale soil column model,and the temperature increments at different water contents were determined as temperature characteristic values (i.e.temperature increments).The heating setup used in the laboratory was consistent with the field setup.Finally,the water content was quantitatively estimated by fitting the relationship between the temperature characteristic value and the volumetric soil water content.

    To calibrate the pore pressure,the FBG pore pressure cell was placed inside a sealed water tank and gradually pressurized to determine the wavelength-pressure correlation.Similarly,for strain calibration,a section of the strain sensing cable was subjected to a tensile test while maintaining a constant temperature,and the strain sensitivity coefficient was obtained.According to the principle of FBG sensing,temperature,water content,pore water pressure,and strain were respectively inverted from the measured wavelengths.

    3.2.2.Surface displacement

    Surface displacements were measured by 41 global navigation satellite system (GNSS) benchmarks and 2 reference points distributed over the slope surface.The reference point is the iron tube placed in a stabilized area.The GNSS observations collected during the same time period(i.e.from January 1,2021,to November 10,2022) as our fiber optic sensing reveal the spatiotemporal pattern of the landslide surface deformation.It provides a valuable comparison for the subsurface monitoring results (Fig.2).

    3.2.3.Reservoir water level(RWL)

    The TGR was first impounded to~135 m a.s.l.in June 2003,then to~172 m a.s.l.in September 2008,and 175 m a.s.l.in October 2010.Since then,the RWL has been regulated to maintain a flood control level between 145 m a.s.l.and 175 m a.s.l.per year.More specifically,the reservoir level is routinely adjusted to a low level during the rainy season from May to September to mitigate the risk of flooding,and rises back to a high level during other months when precipitation is less.

    3.2.4.Meteorological data

    Meteorological data used in this study were collected from a weather station installed on the western part of the landslide.Daily records of precipitation and air temperature since 2021 are available,along with other climate-related variables such as wind velocity,atmospheric pressure,and ultraviolet intensity.The station can record these data every 30 min.

    4.Results and analyses

    In this section,we analyzed the monitoring data (February 11,2021-November 10,2022) to investigate the subsurface thermohydro-poro-mechanical behavior within the sliding masses of the lower Xinpu landslide.Note that the data were missing from June 13-28,2021 due to the anthropogenic interruption of the system power supply.

    4.1.Subsurface multi-physical evolution of sliding masses

    4.1.1.Soil temperature

    The full annual subsurface temperature cycle(ignoring the data gap in gray) since February 11,2021 is reported in Fig.6.The temperature profile,linearly interpolated,shows a depth of around 9 m of the annual thermally active layer,where the soil temperature varies in the range of 16-23°C.Surface temperature measured by the uppermost sensor corresponded well with air temperature records,exhibiting distinct hot and cold seasons.The near-surface temperature varied in the range of 13-38°C (Freifeld et al.,2008).In contrast,in the thermally inactive layer,variations in soil temperature were relatively sluggish,with only a magnitude of about 3°C.Additionally,it seems to be a temperature skew at depths of 9-14 m from early April to mid-October in 2021,but the reason is,at the time being,unknown.

    Fig.7 compares the monthly mean temperature distribution in the borehole FOS3 throughout the monitoring period of 2021 and 2022,along with a stratigraphic column illustrating the material composition from the borehole core samples.It is found that the temperature fluctuation in 2022 was less than that in 2021.This might be partly explained by the progressive consolidation of the borehole backfilling which contributes to more robust temperature measurement.This is also the reason that the thermal active layer interface has not been observed with an obvious downward shift in 2022,despite higher near-surface temperatures than those in 2021(Luo et al.,2022).In addition,the ideal linear temperature distribution (i.e.geothermal gradient) has not been found because this borehole is not deep enough (Rybach and Pfister,1994).

    Fig.7.Monthly mean subsurface temperature profiles in research borehole FOS3,along with the temperature range indicated in gray: (a) 2021 and (b) 2022.The left column profiles stratigraphic distribution according to the drilling core samples: A -C represent silty clay,crushed siltstone with gravel,and gravelly silty clay,respectively.

    4.1.2.Water content

    Fig.8a outlines the linearly interpolated water content distribution and its evolution in the considered time-span.We observed an interface of an annual active thermo-hydraulic layer(marked in red dashed dot-dot line),which was also mostly supported by the aforementioned temperature profile.Within these depths,the atmosphere strongly interacted with shallow soils,maintaining a highly active thermo-hydraulic energy conversion,varying from 0.13 m3/m3to 0.32 m3/m3of the water content.More specifically,it was featuring significant rainfall infiltration and soil evaporation during the rainy season (Blight,1997;Sun et al.,2021,2022).In contrast,water was broadly buffered in shallow soils due to weaker evaporation during the non-rainy season.We noted the evident leaps in water content at shallow depths indicating infiltrated rainfall,such as on July 7 and August 26,2021.Moreover,we argued that the water content within depths of 7 m was overestimated during the approximate initial three months of the monitoring campaign.This could be explained by a measurement artifact caused by the borehole construction,which would be naturally removed over time (Hugentobler et al.,2020).

    Fig.8.Annual hydraulic responses to the hydrometeorological records: (a) Water content profile,and (b) RWL and their fluctuations.

    Strikingly,the evolved water content profile revealed a dynamic groundwater table that fluctuated ranging from 9 m to 15 m(marked in yellow dashed dot-dot line).The groundwater level perfectly matched the RWL regulation,lagging approximately 40 d behind the filling-drawdown of the reservoir water (Prokeˇsová et al.,2013).The groundwater levels inferred from the water content measurements were quite reasonable,since the drilling,for installation of the fiber optic cables in mid-January 2021,had recorded a dive level at a depth of 13 m.The increased water content at the end of March 2021 was confirmed to be caused by agricultural irrigation lasting for several days.These findings,in turn,demonstrate the potential of the UWFBG-based water content evaluation for identifying hydrologic,meteorological,and anthropogenic activities (Froude and Petley,2018;Xu et al.,2022b;Ye et al.,2022).

    4.1.3.Pore water pressure

    As shown in Fig.9,pore water pressures measured by piezometers installed in the borehole within the landslide toe provided very different pressure values.Overall,changes in pore water pressure measured by the uppermost sensor at a depth of 10 m during the monitoring period show a good temporal correlation with the reservoir level,with slight time lagging.More specifically,the increase-decrease of pore pressure was positively correlated with the fill-drawdown of the RWL.It,however,seems that pore pressure changes did not correlate with fluctuations of reservoir water.In addition,it was not observed that some rainfall events should have caused a rapid increase in pore water pressures.This is possibly related to the fact that they were mostly short-duration and high-intensity rainfalls that almost played a role of surface runoff,and partially attributed to flat terrain and borehole top covered by concrete pier.Changes in pore pressures at deeper levels were much smaller than those at the shallower ones,likely due to the lag of pore pressure transfer in deeper soils,as well as the complexity of the stratum(Schulz et al.,2018).

    Fig.9.Annual variations in pore water pressures at depths of interest: (a) Pore pressure versus time;and (b) Daily rainfall,RWL,and their fluctuations.

    4.1.4.Strain

    Fig.10a depicts the annual full-depth strain profile measured since February 11,2021.We observed a prominent peak strain zone at depths of 11-15 m(i.e.close to the conceptualized groundwater table) with a peak strain of 2.1 × 10-3that occurred on July 25,lasting for 45 d.It is noted that the strain peak,after the rapid rise in the reservoir level to 170 m a.s.l.on September 10,2021,gradually decreased to 1.05 × 10-3within about two weeks and remained largely unchanged,with fluctuations of±1×10-4.It suggested that the stress concentration in this zone primarily resulted from the agricultural irrigation event at the end of March,as well as the continued drawdown of reservoir water.The presence of smectite clay at a depth of 12 m facilitated the shrinking-swelling of weathered siltstone,which in turn impacted the cohesion and friction of the soil(Baroˇn et al.,2004).The constant rise of reservoir water since late August increased buoyancy at the landslide foot,which contributed to enhancing the landslide stability and exhibited a strain relaxation(Ferrari et al.,2011;Tang et al.,2019).

    Fig.10.Annual variations in strain: (a) Strain profile;and (b) Records of rainfall and reservoir water levels.

    We further compared strain profiles characterized by monthly mean values throughout the monitoring period to clarify the spatial strain distribution and inter-annual evolution(Fig.11a).There was nearly 50% strain reduction at a depth of 13 m from August to December 2021,which was completely different from that within depths of 2.5-7.5 m (i.e.slightly increased compression strain).It suggests that the high-level reservoir water caused reversible deformation at the landslide toe.In contrast,the strain variability in 2022 is much smaller than that in 2021.One possible reason is that the overall landslide deformation in 2022 is much smaller(Fig.2a).Furthermore,the Yangtze River Basin underwent an unprecedented extreme drought in 2022(Ma and Yuan,2023),which led to a delay in the regulation of RWL compared to previous years (e.g.170 m a.s.l.on September 10,2021,compared to 148 m a.s.l.on September 10,2022).This is the reason that the strain contraction was not clearly observed in 2022.As mentioned earlier,such reservoirinduced reversible strain (i.e.implying reversible deformation)may be transformed into irreversible displacements during cyclic reservoir regulation,which macroscopically manifests as an increasingly flattened surface uplift at the landslide toe.

    Fig.11.Monthly mean subsurface strain profiles in borehole FOS3 in (a) 2021 and (b) 2022.The annual strain ranges are marked in gray.

    4.2.Slope kinematics and its potential controls

    According to the data provided by the monitoring system,the recorded surface deformation was much larger in 2021 than that in 2022 (Fig.12a),and this deformation was mostly concentrated in the area below the “anti-slide” piles row,thus confirming the effectiveness of such countermeasure (Zhang et al.,2023).In particular,the landslide kinematic behavior of the upper part of the unstable slope exhibits almost a constant slow deformation trend(e.g.G16,G20,G21),whereas the area below the piles shows a typical step-like deformation pattern in response to the drawdown of reservoir level and/or concentrated rainfall (e.g.G02,G03,G28,G15).It assumes that these stepwise surface displacements were primarily governed by the constant decrease of RWL with a time lag of about 1-1.5 months;a large amount of rainfall further accelerated the deformation.The period of rapid increase of surface displacement in 2021 was prolonged due to a few high-intensity rainfall events,triggering shallow landslides especially in the lower part of the slope (Fan et al.,2020).In contrast,the slight increase of surface displacements in 2022 might be induced by the drawdown of reservoir levels.The entire landslide exhibited an insignificant displacement trend,even within the most active lower part of the landslide (Fig.2a).The landslide movement behavior characterized by surface displacement at the slide toe and borehole strain is consistent (see Figs.10-12).In addition,the extreme drought also forced a change in water level regulation for the 2022 flood season,representing an extended low water period(i.e.150 m a.s.l.) of at least 1 month.Preliminary results suggest that the deformation mechanism of the Xinpu landslide driven by the drawdown of RWL and rainfall is similar to the Shuping landslide,Taping landslide,and Majiagou landslide (Song et al.,2018;Zhang et al.,2021b,2021c).

    Fig.12.Surface and subsurface kinematics linked to hydrometeorological factors: (a) Lateral surface displacements;(b) Strain time series at depths of interest;and (c) hydrometeorological records.The flood season in 2021 and 2022 are marked by red and blue shaded areas,respectively.

    5.Discussion

    5.1.Thermo-hydro-poro-mechanical behavior within sliding masses

    Temperature,water content,pore water pressure,and strain sensing arrays co-located in a single borehole reveal the effects of annual hydrometeorological cycles on the transient subsurface thermo-hydro-poro-mechanical regime.The annual thermally active layer of the landslide is roughly 9 m thick,which may appropriately shift downward with a warming climate.The water content measurement opens new perspectives for groundwater monitoring,showing an excellent correlation of peak water content evolution with reservoir level regulation and documenting that hydraulic activity within the sliding mass lags behind reservoir levels significantly.The data so far gathered indicate that seasonal pore pressure changes are primarily regulated by reservoir levels rather than rainfall and snowmelt events.The strain sensing arrays allow us to characterize the cumulative deformation reflecting the whole landslide activity,especially at its toe.The data presented confirm that the wet years are more likely to result in substantial deformation (or damage) than dry years (Zscheischler et al.,2018;De Luca et al.,2020).These potential drivers such as rapid drawdown in reservoir levels,intense storms,or extreme temperatures,often act simultaneously(Tichavsky et al.,2019;Kafle et al.,2022).The measured strain during these events usually involves irreversible and reversible components that cannot be ignored.The former is permanent damage at crucial engineering geological interfaces within the landslide driven by hydrometeorological factors,whereas the latter may be related to annual reservoir level fluctuations (Hugentobler et al.,2020).

    To date,it has been challenging to distinguish the irreversible portion from the deformation associated with progressive rock damage in a relatively short monitoring time series.Nevertheless,the aim of this paper is to provide a natural framework for interpreting the thermo-hydro-poro-mechanical responses of the creeping reservoir slide,with a specific focus on the sliding mass.Here we report a preliminary study based on the monitoring results and analyses spanning almost two years.We anticipate that continued monitoring over multiple annual cycles will allow us to develop an integrated framework for multi-physical monitoring and multi-level early warning of landslide hazards(Ye et al.,2022).

    5.2.Implications for understanding landslide kinematics

    As illustrated in Fig.13,in a landslide,there are often various engineering geological interfaces that involve migration of threephase media,transformation of state and progressive changes in physical and mechanical properties,such as base-cover interface,local weak surface,interface between fractured rock and clay,diving surface and infiltration front (Zhu,2023).These interfaces sometimes play a controlling role in the process of gestation,development and occurrence of geological hazards.In particular,the first three are often recognized as potential sliding surfaces,where remarkable thermo-hydro-poro-mechanical changes could occur prior to landslide initiation,and macroscopically manifest as large deformations (Hu et al.,2020;Seguí and Veveakis,2021;Scaringi and Loche,2022).The sliding surfaces and other crucial engineering geological interfaces between bedrock and earth surface may provide inspiring information related to landslide evolution processes (Zhu,2023).For instance,due to the morphological changes related to landslide activity,the slip surface near the slide toe is usually quite deep,but it provides an opportunity to investigate multi-physics of the interfaces between the layers of accumulated materials.Here we report inter-annual monitoring results of thermo-hydro-mechanical responses on the dynamic groundwater interface,as one of the crucial engineering geological interfaces,which reveals the evolution of the time-varying and timedelay water content,as well as the strain accumulation and withdrawal.These interface behaviors can support the interpretation of the contribution mechanism of the filling-drawdown of reservoir water on slope stability,allowing us to link them to local failure events and landslide kinematics(Paronuzzi et al.,2013;Tang et al.,2019).

    Fig.13.Illustration of crucial engineering geological interfaces in a reservoir landslide subjected to different hydrometeorological conditions.The schematic cross-section conceptually cuts through our borehole location (FOS3),revealing annual thermal active layer,base-cover interface,as well as dynamic groundwater interface that could evolve into a potential sliding surface.

    The strain measurements allow some possible mechanisms to be better understood and it seems that the rate of strain accumulation may exacerbate during wet spells after a dry spell(Tichavsky et al.,2019).The combined effect of extreme temperatures and rainfall has been correlated to landslide activity due to the shrinking-swelling effects of clay-rich materials,although the entire landslide kinematics tends to be quasi-steady state during dry spells (Lehmann and Or,2012;Bennett et al.,2016).The“forcing”of the hydrologic cycle,that is the increase of frequency of precipitation extremes following extreme high temperatures,is likely to amplify the magnitude of sliding events (Cohen-Waeber et al.,2018;De Luca et al.,2020).We hypothesize that the transition from a quasi-stable slope to a failing landslide is a long-term and continuous process.Therefore,the preliminary findings of this paper will be critically reviewed in future studies based on our continuously updated monitoring dataset.

    5.3.Applicability and limitations of UWFBG sensing nerves

    Compared to the mature FBG and DFOS technologies,UWFBG allows near real-time acquisition of large amounts of multiphysical data along the fiber optic nerves with high spatial resolution such as temperature,water content,and strain(see Table S1).Specifically,the sensing elements on a single fiber optic nerve can reach over a thousand at an interval of 1 m,which is far less than FBG with comparable measurement accuracy,despite that both of them can collect data remotely and in real time.UWFBG also significantly outperforms DFOS in terms of automatization and integration of monitoring system(Zeni et al.,2015;Shi et al.,2021).Therefore,UWFBG technique is suitable for deep-seated landslides or rockslides,especially for discrimination of critical engineering geological interface (Zhu,2023).

    A major limitation so far is that only a few pore water pressure cells can be placed in a borehole with a diameter of 110 mm,given that each pressure cell contains a pie-shaped permeable stones with a diameter of 50 mm.In this aspect,unlike temperature,moisture and strain,we can only acquire depth-specific pore water pressures.Therefore,the determination of the installation depth of the pore pressure cell is challengeable and requires prior knowledge of the formation conditions.In other words,deploying such multi-physical boreholes in critical deformation areas at the lower,middle,and upper parts of the landslide could provide a possibility of understanding the subsurface thermo-hydro-poro-mechanical evolution processes of global large-scale landslides.

    6.Conclusions

    The present paper describes a multi-physical dataset collected from a novel borehole-based high-resolution fiber optic nerve sensing system using UWFBG technology.The system was installed to continuously record changes in below-ground soil temperature,water content,pore water pressure,and strain with a resolution of 1 m at the toe of the Xinpu landslide in the Three Gorges Reservoir Area,China,to investigate thermo-hydro-poro-mechanical responses of the unstable processes subject to periodical reservoir regulation and seasonal weather events.Here we report a preliminary study that involves monitoring data covering two recent flood periods,and the following conclusions can be drawn:

    (1) The annual thermally active layer of the landslide had a critical depth of about 9 m and might appropriately move downward in warmer years.The spatio-temporal evolution of the soil water content revealed a dynamic groundwater table fluctuating from 9 m to 15 m with a time lag of about 40 d with respect of the RWL variations.

    (2) Unlike the sluggish water pressure variation in the deep layers,the pore water pressure in the shallow layer at a depth of 10 m showed a close connection with the fillingdrawdown of the RWL.

    (3) The peak strain along the borehole profile occurred roughly at a depth of 13 m and experienced a periodical tension and relaxation to annual hydrometeorological cycles.The magnitude of such strain variations was regulated by wet and dry hydrological extremes.As expected,wet years tend to cause much larger strain accumulation,triggering more significant landslide displacements,than dry years.

    Declaration of competing interest

    The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

    Acknowledgments

    We acknowledge the funding support from the National Science Fund for Distinguished Young Scholars of National Natural Science Foundation of China(Grant No.42225702)and the National Natural Science Foundation of China (Grant No.42077235).Special thanks go to the Institute of Exploration Technology,China Geological Survey,Chengdu,China,for providing geological information of the study area.We also appreciate Xing Zheng,Jia Wang and Bing Wu in Nanjing University for their help in data processing and analyses.The first author particularly thanks Dr.Yu-Zhou Hou in 3rd Construction Co.,Ltd.of China Construction 5th Engineering Bureau for his insightful comments to result analyses.

    Appendix A.Supplementary data

    Supplementary data to this article can be found online at https://doi.org/10.1016/j.jrmge.2023.04.004.

    欧美精品国产亚洲| 国产一级毛片在线| 午夜免费观看性视频| 少妇的逼水好多| 午夜免费男女啪啪视频观看| 色网站视频免费| 男的添女的下面高潮视频| 搡女人真爽免费视频火全软件| 人人妻人人添人人爽欧美一区卜| 国产乱来视频区| 成年动漫av网址| 国产精品一区二区在线观看99| 亚洲一码二码三码区别大吗| 成人手机av| 国产熟女午夜一区二区三区| 欧美日韩视频精品一区| 少妇人妻 视频| 少妇人妻久久综合中文| 国产成人精品婷婷| 少妇的逼好多水| 丰满迷人的少妇在线观看| 亚洲国产精品一区三区| 日韩成人av中文字幕在线观看| 中国国产av一级| 91精品伊人久久大香线蕉| 成人无遮挡网站| 又黄又粗又硬又大视频| 五月天丁香电影| 五月天丁香电影| 国产精品99久久99久久久不卡 | 久久99精品国语久久久| 26uuu在线亚洲综合色| 日本wwww免费看| av片东京热男人的天堂| 一本—道久久a久久精品蜜桃钙片| 久久99蜜桃精品久久| 久久久久精品久久久久真实原创| 亚洲成人手机| 在线观看人妻少妇| kizo精华| 久久久亚洲精品成人影院| 波多野结衣一区麻豆| 菩萨蛮人人尽说江南好唐韦庄| 久久久国产一区二区| 啦啦啦中文免费视频观看日本| 国产亚洲一区二区精品| 久久鲁丝午夜福利片| 亚洲天堂av无毛| 人人澡人人妻人| 麻豆乱淫一区二区| 精品人妻熟女毛片av久久网站| 国产一区二区三区av在线| av天堂久久9| 秋霞在线观看毛片| 国产探花极品一区二区| 搡女人真爽免费视频火全软件| 在线观看免费高清a一片| 国产成人欧美| 国产男女超爽视频在线观看| 精品久久蜜臀av无| 51国产日韩欧美| 久久精品aⅴ一区二区三区四区 | 99热国产这里只有精品6| 啦啦啦中文免费视频观看日本| 蜜桃国产av成人99| 国产高清不卡午夜福利| 又粗又硬又长又爽又黄的视频| 久久国内精品自在自线图片| 午夜免费观看性视频| 91午夜精品亚洲一区二区三区| 午夜久久久在线观看| 亚洲欧美中文字幕日韩二区| 巨乳人妻的诱惑在线观看| 最近最新中文字幕大全免费视频 | 91精品三级在线观看| 亚洲美女黄色视频免费看| 久久热在线av| 国产在线一区二区三区精| 人人妻人人澡人人爽人人夜夜| 国产一区二区在线观看日韩| 性色avwww在线观看| 精品亚洲成a人片在线观看| 成人免费观看视频高清| 十分钟在线观看高清视频www| 99久久精品国产国产毛片| 国产色婷婷99| 国产精品无大码| 欧美日韩视频高清一区二区三区二| 性色av一级| 亚洲精华国产精华液的使用体验| 人人妻人人澡人人爽人人夜夜| 永久网站在线| freevideosex欧美| 久久久久久伊人网av| 少妇的丰满在线观看| 男女国产视频网站| 国产视频首页在线观看| 成人漫画全彩无遮挡| 久久韩国三级中文字幕| 女性被躁到高潮视频| 久久久久国产精品人妻一区二区| 韩国高清视频一区二区三区| 老熟女久久久| 在线观看免费视频网站a站| 深夜精品福利| 国产av码专区亚洲av| 91aial.com中文字幕在线观看| 人人妻人人添人人爽欧美一区卜| 国产亚洲精品第一综合不卡 | 最近的中文字幕免费完整| 久久久欧美国产精品| 亚洲精品美女久久av网站| 最近最新中文字幕大全免费视频 | 亚洲美女黄色视频免费看| 精品久久久久久电影网| 制服人妻中文乱码| xxxhd国产人妻xxx| 国内精品宾馆在线| 天天躁夜夜躁狠狠躁躁| 久久婷婷青草| 人妻一区二区av| a级片在线免费高清观看视频| 欧美日韩一区二区视频在线观看视频在线| av在线观看视频网站免费| 亚洲国产av新网站| 国产精品欧美亚洲77777| 美女内射精品一级片tv| 亚洲欧美中文字幕日韩二区| 天天躁夜夜躁狠狠躁躁| 久久精品国产自在天天线| 日本wwww免费看| 97在线视频观看| 卡戴珊不雅视频在线播放| 国产精品嫩草影院av在线观看| 欧美成人午夜精品| 久久婷婷青草| 中国美白少妇内射xxxbb| 黄色一级大片看看| 高清黄色对白视频在线免费看| 国产精品无大码| 久久久久久伊人网av| 免费大片18禁| 高清欧美精品videossex| 国产精品久久久久久久久免| 宅男免费午夜| 伊人亚洲综合成人网| xxx大片免费视频| 精品久久久精品久久久| 欧美日韩亚洲高清精品| 成年人午夜在线观看视频| 国产免费一区二区三区四区乱码| 国产69精品久久久久777片| 80岁老熟妇乱子伦牲交| 一级,二级,三级黄色视频| 欧美精品国产亚洲| 人妻少妇偷人精品九色| av视频免费观看在线观看| 丝袜人妻中文字幕| 免费av中文字幕在线| 狂野欧美激情性bbbbbb| 久久久欧美国产精品| 超色免费av| 在线观看国产h片| 久久久久国产精品人妻一区二区| 少妇被粗大的猛进出69影院 | 少妇高潮的动态图| 精品99又大又爽又粗少妇毛片| 人妻 亚洲 视频| 国产精品久久久久久久电影| 日韩制服骚丝袜av| 久久影院123| 9热在线视频观看99| 国产在线一区二区三区精| 亚洲国产看品久久| 成人亚洲精品一区在线观看| 日韩av不卡免费在线播放| 精品国产一区二区久久| 久久影院123| 国产精品熟女久久久久浪| 天天躁夜夜躁狠狠躁躁| 精品视频人人做人人爽| 波野结衣二区三区在线| 一区二区三区精品91| 日韩三级伦理在线观看| 美女内射精品一级片tv| 国产精品蜜桃在线观看| 亚洲精华国产精华液的使用体验| 1024视频免费在线观看| xxxhd国产人妻xxx| 亚洲精品久久成人aⅴ小说| 亚洲熟女精品中文字幕| 97人妻天天添夜夜摸| 十八禁网站网址无遮挡| 亚洲欧洲日产国产| 日日啪夜夜爽| 男人操女人黄网站| 久久久久精品性色| 九九在线视频观看精品| 日韩成人av中文字幕在线观看| 国产亚洲精品第一综合不卡 | 热99国产精品久久久久久7| 国产毛片在线视频| 草草在线视频免费看| 国产成人av激情在线播放| 亚洲伊人色综图| 久久鲁丝午夜福利片| 成人国产麻豆网| 国产精品久久久久久精品电影小说| 高清毛片免费看| 插逼视频在线观看| 免费观看性生交大片5| 亚洲欧美色中文字幕在线| 91国产中文字幕| 亚洲综合精品二区| 亚洲精品美女久久av网站| 久久人人爽人人爽人人片va| 香蕉精品网在线| 男女免费视频国产| 国产成人a∨麻豆精品| 亚洲丝袜综合中文字幕| 日日撸夜夜添| 在线观看www视频免费| 精品国产露脸久久av麻豆| 精品第一国产精品| 日韩电影二区| 亚洲av.av天堂| 人妻系列 视频| 亚洲成人av在线免费| 亚洲av电影在线观看一区二区三区| 精品亚洲乱码少妇综合久久| 国产精品一国产av| 美女国产高潮福利片在线看| 国产高清不卡午夜福利| 久久精品久久精品一区二区三区| 熟女av电影| 成人国产av品久久久| 美女内射精品一级片tv| 欧美人与善性xxx| 少妇的丰满在线观看| 国产成人欧美| 蜜臀久久99精品久久宅男| 最后的刺客免费高清国语| 深夜精品福利| 久久久国产一区二区| 亚洲美女搞黄在线观看| 91aial.com中文字幕在线观看| 色婷婷久久久亚洲欧美| 香蕉丝袜av| 综合色丁香网| a级片在线免费高清观看视频| 18禁裸乳无遮挡动漫免费视频| 国产69精品久久久久777片| 亚洲av日韩在线播放| 欧美日韩视频高清一区二区三区二| 日本91视频免费播放| av免费观看日本| 久久久久久人人人人人| freevideosex欧美| 97人妻天天添夜夜摸| 久久久久久久久久人人人人人人| 久久 成人 亚洲| 久久久欧美国产精品| 久久久久视频综合| 十八禁高潮呻吟视频| 两性夫妻黄色片 | 99精国产麻豆久久婷婷| 大码成人一级视频| 日韩一区二区三区影片| 天天躁夜夜躁狠狠久久av| 国产成人a∨麻豆精品| 国产精品不卡视频一区二区| 自拍欧美九色日韩亚洲蝌蚪91| 日本欧美国产在线视频| 亚洲精品美女久久av网站| 97在线视频观看| 黄网站色视频无遮挡免费观看| 久久鲁丝午夜福利片| 男女国产视频网站| 成人午夜精彩视频在线观看| 成人国语在线视频| 极品少妇高潮喷水抽搐| 免费播放大片免费观看视频在线观看| 亚洲国产欧美日韩在线播放| 日韩人妻精品一区2区三区| 9191精品国产免费久久| 国产一区二区激情短视频 | a级毛色黄片| 狠狠婷婷综合久久久久久88av| 久久人人爽av亚洲精品天堂| 赤兔流量卡办理| 亚洲一级一片aⅴ在线观看| 国产欧美日韩一区二区三区在线| 波多野结衣一区麻豆| 51国产日韩欧美| 欧美日韩一区二区视频在线观看视频在线| 成人无遮挡网站| 久久鲁丝午夜福利片| 成人影院久久| 赤兔流量卡办理| 日韩中文字幕视频在线看片| 国产亚洲av片在线观看秒播厂| 久久精品久久久久久噜噜老黄| 国产亚洲一区二区精品| 欧美日韩成人在线一区二区| 韩国av在线不卡| 天天操日日干夜夜撸| av电影中文网址| 看十八女毛片水多多多| 卡戴珊不雅视频在线播放| 亚洲性久久影院| 国国产精品蜜臀av免费| 日本欧美视频一区| 国产男人的电影天堂91| 欧美97在线视频| 日本vs欧美在线观看视频| 纯流量卡能插随身wifi吗| 日本与韩国留学比较| 国产精品99久久99久久久不卡 | 亚洲经典国产精华液单| 人人妻人人添人人爽欧美一区卜| 欧美+日韩+精品| www.熟女人妻精品国产 | 女人精品久久久久毛片| 国产日韩欧美亚洲二区| 色婷婷av一区二区三区视频| 国产爽快片一区二区三区| 最新的欧美精品一区二区| 人妻少妇偷人精品九色| 精品人妻偷拍中文字幕| 精品人妻偷拍中文字幕| 晚上一个人看的免费电影| 日韩一本色道免费dvd| 久久久久网色| 国产成人精品无人区| 午夜福利影视在线免费观看| 日本免费在线观看一区| 国产亚洲精品久久久com| 色网站视频免费| 丝袜美足系列| 色网站视频免费| 嫩草影院入口| 最近中文字幕2019免费版| 毛片一级片免费看久久久久| 免费av不卡在线播放| 亚洲欧美日韩卡通动漫| 久久久久久久久久久免费av| 久久国产亚洲av麻豆专区| 欧美日本中文国产一区发布| 日本爱情动作片www.在线观看| 国产伦理片在线播放av一区| 观看美女的网站| 日韩成人伦理影院| 成人漫画全彩无遮挡| 亚洲精品一二三| √禁漫天堂资源中文www| 巨乳人妻的诱惑在线观看| 各种免费的搞黄视频| 老女人水多毛片| 国产一区二区在线观看av| 日韩制服丝袜自拍偷拍| 看非洲黑人一级黄片| 晚上一个人看的免费电影| 波多野结衣一区麻豆| 国产无遮挡羞羞视频在线观看| 一区二区三区精品91| 日韩一区二区三区影片| 丝袜人妻中文字幕| 热re99久久精品国产66热6| 18在线观看网站| 一级黄片播放器| 日韩av不卡免费在线播放| 熟女人妻精品中文字幕| 这个男人来自地球电影免费观看 | 一本久久精品| 制服丝袜香蕉在线| 国产 精品1| 国产精品 国内视频| 日韩 亚洲 欧美在线| 日韩人妻精品一区2区三区| 国产精品99久久99久久久不卡 | av.在线天堂| 欧美人与性动交α欧美软件 | 久久人人爽人人片av| 亚洲av国产av综合av卡| 国产精品国产三级国产专区5o| 汤姆久久久久久久影院中文字幕| 啦啦啦啦在线视频资源| 成人黄色视频免费在线看| 亚洲成人一二三区av| 中文欧美无线码| 999精品在线视频| 日本色播在线视频| 天天躁夜夜躁狠狠躁躁| 国产成人免费观看mmmm| 日韩视频在线欧美| 亚洲欧洲日产国产| a级片在线免费高清观看视频| 精品第一国产精品| 在线天堂中文资源库| 在线观看免费高清a一片| 久久女婷五月综合色啪小说| 又黄又粗又硬又大视频| 亚洲av欧美aⅴ国产| 最近的中文字幕免费完整| 最新的欧美精品一区二区| 国产亚洲一区二区精品| 精品国产乱码久久久久久小说| 超碰97精品在线观看| 国产一区二区三区av在线| 热re99久久精品国产66热6| 男人操女人黄网站| 亚洲图色成人| 一本色道久久久久久精品综合| 免费黄频网站在线观看国产| 日本猛色少妇xxxxx猛交久久| 免费人成在线观看视频色| 欧美国产精品一级二级三级| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 午夜免费鲁丝| 少妇高潮的动态图| 亚洲精品国产av成人精品| 久久99一区二区三区| 日韩 亚洲 欧美在线| 国产精品国产三级专区第一集| 韩国高清视频一区二区三区| 国产欧美日韩综合在线一区二区| 乱码一卡2卡4卡精品| xxx大片免费视频| 日韩欧美一区视频在线观看| 99热这里只有是精品在线观看| av播播在线观看一区| 国产精品蜜桃在线观看| 女人精品久久久久毛片| 成年动漫av网址| 久久午夜综合久久蜜桃| 五月开心婷婷网| 久久精品国产综合久久久 | 日韩制服丝袜自拍偷拍| 各种免费的搞黄视频| 精品一品国产午夜福利视频| 精品第一国产精品| 久久久国产精品麻豆| 亚洲欧美清纯卡通| 天天躁夜夜躁狠狠久久av| 天堂俺去俺来也www色官网| 69精品国产乱码久久久| 精品亚洲乱码少妇综合久久| 日日撸夜夜添| 精品人妻在线不人妻| 两个人看的免费小视频| 大香蕉久久网| 99久国产av精品国产电影| 91午夜精品亚洲一区二区三区| 啦啦啦啦在线视频资源| 精品一区二区免费观看| 日韩精品免费视频一区二区三区 | 午夜影院在线不卡| 欧美精品av麻豆av| 大香蕉久久网| 免费日韩欧美在线观看| freevideosex欧美| 免费黄网站久久成人精品| 一级毛片 在线播放| 成年女人在线观看亚洲视频| 天天躁夜夜躁狠狠躁躁| 久久精品国产亚洲av香蕉五月 | 欧美不卡视频在线免费观看 | 精品亚洲成a人片在线观看| 久久香蕉激情| 亚洲欧美日韩另类电影网站| 亚洲七黄色美女视频| netflix在线观看网站| 一区二区日韩欧美中文字幕| 久久久久久亚洲精品国产蜜桃av| 飞空精品影院首页| 一边摸一边抽搐一进一小说 | 国产精品一区二区在线观看99| 成人18禁高潮啪啪吃奶动态图| 国产不卡一卡二| 99久久国产精品久久久| 国产高清视频在线播放一区| 嫁个100分男人电影在线观看| 久热爱精品视频在线9| 亚洲七黄色美女视频| 午夜福利影视在线免费观看| 五月开心婷婷网| 久久 成人 亚洲| 亚洲av美国av| 亚洲精品国产色婷婷电影| 欧美av亚洲av综合av国产av| 日本wwww免费看| 黑人巨大精品欧美一区二区mp4| 国产精品自产拍在线观看55亚洲 | 欧美激情 高清一区二区三区| 国产欧美日韩综合在线一区二区| 精品国产亚洲在线| 欧美日韩亚洲综合一区二区三区_| 高潮久久久久久久久久久不卡| 欧美精品高潮呻吟av久久| 国产一区二区三区综合在线观看| 黑人猛操日本美女一级片| 久久久国产欧美日韩av| 久99久视频精品免费| 中文亚洲av片在线观看爽 | 下体分泌物呈黄色| 久久人妻福利社区极品人妻图片| 窝窝影院91人妻| 国产精品二区激情视频| 免费在线观看日本一区| 1024香蕉在线观看| av一本久久久久| 久热这里只有精品99| 色综合欧美亚洲国产小说| 乱人伦中国视频| 欧美日韩一级在线毛片| 12—13女人毛片做爰片一| 91国产中文字幕| 国产精品亚洲一级av第二区| 亚洲精华国产精华精| 水蜜桃什么品种好| 丝袜美腿诱惑在线| 99国产精品99久久久久| 99在线人妻在线中文字幕 | 久久亚洲精品不卡| 国产亚洲欧美98| 韩国精品一区二区三区| 国产熟女午夜一区二区三区| 精品少妇一区二区三区视频日本电影| 高清黄色对白视频在线免费看| 一本大道久久a久久精品| 女性生殖器流出的白浆| 国产激情久久老熟女| 久久精品国产a三级三级三级| 女警被强在线播放| 精品国产亚洲在线| 国产精品成人在线| 黑人巨大精品欧美一区二区mp4| 天天操日日干夜夜撸| 欧美日韩精品网址| 久久久久久久精品吃奶| 成年人午夜在线观看视频| 国产亚洲精品久久久久5区| 欧美日韩瑟瑟在线播放| 欧美日韩视频精品一区| 精品一区二区三卡| 久久青草综合色| 99国产精品一区二区蜜桃av | 欧美一级毛片孕妇| 在线天堂中文资源库| 手机成人av网站| 黄色片一级片一级黄色片| 久久久精品免费免费高清| 色综合婷婷激情| 最新的欧美精品一区二区| 国产亚洲精品第一综合不卡| 91老司机精品| 男女下面插进去视频免费观看| 午夜精品国产一区二区电影| 亚洲国产中文字幕在线视频| 女性生殖器流出的白浆| 亚洲精品成人av观看孕妇| 国产亚洲欧美98| 高清视频免费观看一区二区| 757午夜福利合集在线观看| 久99久视频精品免费| 午夜福利一区二区在线看| 91在线观看av| 成年人黄色毛片网站| 99re在线观看精品视频| 一区二区三区国产精品乱码| 18在线观看网站| 波多野结衣av一区二区av| 99久久综合精品五月天人人| 人人妻人人爽人人添夜夜欢视频| 超碰成人久久| 久99久视频精品免费| 在线永久观看黄色视频| netflix在线观看网站| 久久性视频一级片| 精品一区二区三卡| 色综合婷婷激情| 国产精华一区二区三区| 免费少妇av软件| 国产精品电影一区二区三区 | 午夜福利乱码中文字幕| netflix在线观看网站| 青草久久国产| 精品视频人人做人人爽| 99久久99久久久精品蜜桃| 国产精品偷伦视频观看了| 国产真人三级小视频在线观看| 真人做人爱边吃奶动态| 嫁个100分男人电影在线观看| 国产又爽黄色视频| 男男h啪啪无遮挡| 国产精品偷伦视频观看了| 精品国产一区二区三区四区第35| 美女福利国产在线| 18禁裸乳无遮挡免费网站照片 | 精品无人区乱码1区二区| 欧美黄色片欧美黄色片| 亚洲一区中文字幕在线| 性少妇av在线| 日本黄色日本黄色录像| 在线观看66精品国产| 可以免费在线观看a视频的电影网站| 19禁男女啪啪无遮挡网站| 在线观看免费视频日本深夜| 国产一区二区三区视频了| 国产成人啪精品午夜网站| 叶爱在线成人免费视频播放| 飞空精品影院首页| 搡老乐熟女国产| 欧美精品高潮呻吟av久久| 丁香六月欧美|