• <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| 又黄又爽又免费观看的视频| 欧美中文日本在线观看视频| 最近中文字幕高清免费大全6| 综合色av麻豆| 超碰av人人做人人爽久久| 日韩欧美免费精品| 亚洲精品久久国产高清桃花| 一级毛片我不卡| 日本欧美国产在线视频| 亚洲欧美中文字幕日韩二区| 97热精品久久久久久| 久久韩国三级中文字幕| 一个人观看的视频www高清免费观看| 一夜夜www| 菩萨蛮人人尽说江南好唐韦庄 | 热99在线观看视频| 色在线成人网| av国产免费在线观看| 成人三级黄色视频| 亚洲美女黄片视频| 婷婷精品国产亚洲av在线| 日韩一本色道免费dvd| 日韩欧美精品v在线| 欧美一区二区精品小视频在线| 人人妻人人澡人人爽人人夜夜 | 亚洲婷婷狠狠爱综合网| 欧美成人a在线观看| 久久亚洲精品不卡| 亚洲美女视频黄频| 草草在线视频免费看| 搡老岳熟女国产| 亚洲美女黄片视频| 成人av一区二区三区在线看| 男女啪啪激烈高潮av片| 神马国产精品三级电影在线观看| 欧美丝袜亚洲另类| 国产精品日韩av在线免费观看| 国产国拍精品亚洲av在线观看| www日本黄色视频网| 日本五十路高清| 精品一区二区三区视频在线观看免费| 欧美最新免费一区二区三区| 成人亚洲欧美一区二区av| 欧美人与善性xxx| 亚洲人成网站在线播放欧美日韩| 免费av不卡在线播放| 在线观看免费视频日本深夜| 国产成人精品久久久久久| 麻豆久久精品国产亚洲av| 99久久久亚洲精品蜜臀av| 九九在线视频观看精品| 一边摸一边抽搐一进一小说| 18禁在线无遮挡免费观看视频 | 久久久久久久久中文| 日日摸夜夜添夜夜添小说| 成人美女网站在线观看视频| 国产伦精品一区二区三区四那| 午夜日韩欧美国产| 搞女人的毛片| 亚洲精品一卡2卡三卡4卡5卡| h日本视频在线播放| 熟女电影av网| 三级国产精品欧美在线观看| 国内少妇人妻偷人精品xxx网站| 亚洲精品亚洲一区二区| 久久人人精品亚洲av| 国产精品av视频在线免费观看| 热99在线观看视频| 久久精品久久久久久噜噜老黄 | 一本精品99久久精品77| 啦啦啦观看免费观看视频高清| 国产色婷婷99| 国产亚洲精品综合一区在线观看| 免费av观看视频| 欧美绝顶高潮抽搐喷水| 熟女电影av网| 悠悠久久av| 日本免费一区二区三区高清不卡| 此物有八面人人有两片| 欧美日本亚洲视频在线播放| avwww免费| 国产精品久久久久久精品电影| 久久精品夜夜夜夜夜久久蜜豆| www日本黄色视频网| 男女下面进入的视频免费午夜| 五月伊人婷婷丁香| 国产精品久久久久久精品电影| 51国产日韩欧美| 人妻少妇偷人精品九色| 色视频www国产| 中文字幕精品亚洲无线码一区| 在现免费观看毛片| 国产一区二区激情短视频| 免费看美女性在线毛片视频| 69av精品久久久久久| 日韩av在线大香蕉| 给我免费播放毛片高清在线观看| 99在线视频只有这里精品首页| 亚洲成人久久爱视频| 亚洲成a人片在线一区二区| 日韩高清综合在线| 最后的刺客免费高清国语| 香蕉av资源在线| 一区福利在线观看| 久久99热这里只有精品18| 国产国拍精品亚洲av在线观看| 国产一区二区三区在线臀色熟女| 一级毛片久久久久久久久女| 成人午夜高清在线视频| 少妇猛男粗大的猛烈进出视频 | 伦理电影大哥的女人| 99久久久亚洲精品蜜臀av| 亚洲精品一区av在线观看| 草草在线视频免费看| 香蕉av资源在线| 亚洲熟妇熟女久久| 欧美性感艳星| 深爱激情五月婷婷| 18禁裸乳无遮挡免费网站照片| 国产欧美日韩一区二区精品| 国产精品久久久久久久久免| 日韩欧美国产在线观看| 国产精品野战在线观看| 国产免费男女视频| 99视频精品全部免费 在线| 白带黄色成豆腐渣| 国产精品一区二区三区四区免费观看 | 日韩中字成人| 亚洲自拍偷在线| 男人舔奶头视频| 亚洲成人中文字幕在线播放| eeuss影院久久| АⅤ资源中文在线天堂| 国产精品精品国产色婷婷| 男女视频在线观看网站免费| 亚洲va在线va天堂va国产| 最近2019中文字幕mv第一页| 蜜桃亚洲精品一区二区三区| 国产精品一二三区在线看| 亚洲久久久久久中文字幕| 久久久久精品国产欧美久久久| 午夜视频国产福利| 又黄又爽又免费观看的视频| 欧美区成人在线视频| 欧美日韩在线观看h| 色吧在线观看| 色哟哟·www| 国产高清有码在线观看视频| 1024手机看黄色片| 蜜桃亚洲精品一区二区三区| 天堂网av新在线| 尾随美女入室| 国产高清视频在线播放一区| 久久人人精品亚洲av| 在线观看一区二区三区| 天堂√8在线中文| 亚洲av第一区精品v没综合| 亚洲欧美日韩高清专用| 免费看av在线观看网站| 亚洲在线自拍视频| 久久久久免费精品人妻一区二区| 亚洲电影在线观看av| 神马国产精品三级电影在线观看| 午夜a级毛片| 免费高清视频大片| 中国国产av一级| 少妇猛男粗大的猛烈进出视频 | 黄色欧美视频在线观看| 国产伦一二天堂av在线观看| 亚洲av五月六月丁香网| 日本-黄色视频高清免费观看| 免费看a级黄色片| 亚州av有码| 俄罗斯特黄特色一大片| or卡值多少钱| 国模一区二区三区四区视频| 国产中年淑女户外野战色| 免费看日本二区| 亚洲丝袜综合中文字幕| 亚洲国产精品合色在线| 亚洲精品亚洲一区二区| 精品乱码久久久久久99久播| 偷拍熟女少妇极品色| 禁无遮挡网站| 亚洲成人久久爱视频| 欧美另类亚洲清纯唯美| 看黄色毛片网站| 成年版毛片免费区| 日本免费一区二区三区高清不卡| 婷婷六月久久综合丁香| 国产黄a三级三级三级人| 欧美成人一区二区免费高清观看| 在线免费观看不下载黄p国产| 国产综合懂色| 小说图片视频综合网站| 亚洲美女黄片视频| 精品午夜福利在线看| 亚洲欧美日韩卡通动漫| 国产人妻一区二区三区在| av黄色大香蕉| 精品少妇黑人巨大在线播放 | 亚州av有码| 国产成人a∨麻豆精品| 少妇人妻精品综合一区二区 | 亚洲国产精品成人综合色| 在线免费观看不下载黄p国产| 亚洲天堂国产精品一区在线| 美女被艹到高潮喷水动态| 大香蕉久久网| 日韩人妻高清精品专区| 免费不卡的大黄色大毛片视频在线观看 | 高清午夜精品一区二区三区 | 大型黄色视频在线免费观看| 亚洲国产欧洲综合997久久,| 夜夜爽天天搞| 精品久久久久久成人av| 看免费成人av毛片| 国产三级在线视频| 淫秽高清视频在线观看| 国产麻豆成人av免费视频| 真实男女啪啪啪动态图| 国语自产精品视频在线第100页| 日韩人妻高清精品专区| 亚州av有码| 日韩高清综合在线| 老师上课跳d突然被开到最大视频| 久久精品国产亚洲av涩爱 | 精品国内亚洲2022精品成人| 看十八女毛片水多多多| 欧美bdsm另类| 免费观看人在逋| 乱系列少妇在线播放| 观看免费一级毛片| 国产女主播在线喷水免费视频网站 | 色5月婷婷丁香| www日本黄色视频网| 欧美一区二区国产精品久久精品| 高清毛片免费观看视频网站| 非洲黑人性xxxx精品又粗又长| 精品久久久噜噜| 变态另类成人亚洲欧美熟女| 搡女人真爽免费视频火全软件 | 一卡2卡三卡四卡精品乱码亚洲| 麻豆一二三区av精品| 国产一级毛片七仙女欲春2| 久久精品国产亚洲av香蕉五月| 中文字幕熟女人妻在线| 亚洲欧美中文字幕日韩二区| 亚洲天堂国产精品一区在线| 婷婷精品国产亚洲av| www日本黄色视频网| 九九热线精品视视频播放| 欧美色视频一区免费| 欧美另类亚洲清纯唯美| 成人永久免费在线观看视频| 22中文网久久字幕| 日本精品一区二区三区蜜桃| а√天堂www在线а√下载| 在线看三级毛片| 黄色欧美视频在线观看| 国产精品亚洲一级av第二区| 精华霜和精华液先用哪个| 综合色av麻豆| 麻豆国产av国片精品| 国产大屁股一区二区在线视频| 大香蕉久久网| 美女免费视频网站| 舔av片在线| 国产黄色视频一区二区在线观看 | 亚洲自偷自拍三级| av视频在线观看入口| 麻豆av噜噜一区二区三区| 波多野结衣巨乳人妻| 1000部很黄的大片| 亚洲精品色激情综合| 午夜福利在线观看吧| 成人亚洲精品av一区二区| 免费av观看视频| 中文字幕av成人在线电影| 亚洲av免费在线观看| 中国美女看黄片| 国产精品久久电影中文字幕| 看黄色毛片网站| 国产乱人视频| 亚洲av成人av| 国内精品宾馆在线| 男女那种视频在线观看| 亚洲成人av在线免费| 亚洲精品色激情综合| 国产麻豆成人av免费视频| 亚洲av二区三区四区| a级毛片免费高清观看在线播放| 成人鲁丝片一二三区免费| 午夜福利在线观看免费完整高清在 | 99久久精品一区二区三区| 欧美日本亚洲视频在线播放| 俄罗斯特黄特色一大片| 天堂影院成人在线观看| 我的老师免费观看完整版| 有码 亚洲区| 国产视频内射| 久久久久久久久中文| 亚洲国产精品sss在线观看| 毛片女人毛片| 日韩,欧美,国产一区二区三区 | 亚洲美女搞黄在线观看 | 亚洲无线在线观看| 精品人妻熟女av久视频| 国内精品一区二区在线观看| 18禁黄网站禁片免费观看直播| 国产亚洲av嫩草精品影院| 青春草视频在线免费观看| 免费看a级黄色片| 少妇裸体淫交视频免费看高清| 国产精品久久久久久久久免| 亚洲成人av在线免费| 精品一区二区免费观看| 亚洲精品国产av成人精品 | 99九九线精品视频在线观看视频| 精品久久久久久成人av| 99视频精品全部免费 在线| 久久久久久久久大av| 日韩高清综合在线| 精品久久久久久久人妻蜜臀av| 精华霜和精华液先用哪个| 久99久视频精品免费| 看十八女毛片水多多多| 激情 狠狠 欧美| 啦啦啦观看免费观看视频高清| 高清午夜精品一区二区三区 | 别揉我奶头~嗯~啊~动态视频| 午夜免费男女啪啪视频观看 | av在线老鸭窝| 18禁黄网站禁片免费观看直播| 国产欧美日韩一区二区精品| 欧美日本视频| 国产高清视频在线播放一区| 麻豆国产97在线/欧美| 中文字幕精品亚洲无线码一区| 在线国产一区二区在线| av中文乱码字幕在线| 99九九线精品视频在线观看视频| 在线看三级毛片| 91狼人影院| 91av网一区二区| 欧美中文日本在线观看视频| 一级毛片aaaaaa免费看小| 成人毛片a级毛片在线播放| 亚洲精品粉嫩美女一区| 免费观看的影片在线观看| 国产aⅴ精品一区二区三区波| 小说图片视频综合网站| 久久亚洲精品不卡| 美女cb高潮喷水在线观看| 午夜免费男女啪啪视频观看 | 亚洲国产精品成人久久小说 | 少妇熟女aⅴ在线视频| 亚州av有码| 国产美女午夜福利| 国产精品一区二区三区四区免费观看 | 婷婷六月久久综合丁香| 成人漫画全彩无遮挡| 一本精品99久久精品77| 国产精品福利在线免费观看| 97热精品久久久久久| 午夜视频国产福利| 国产精品av视频在线免费观看| 毛片女人毛片| 国内精品久久久久精免费| 免费看a级黄色片| 一个人观看的视频www高清免费观看| 老师上课跳d突然被开到最大视频| 波多野结衣高清无吗| 日本爱情动作片www.在线观看 | 免费电影在线观看免费观看| 男人狂女人下面高潮的视频| 日韩精品青青久久久久久| 一级毛片我不卡| 亚洲三级黄色毛片| 男人舔奶头视频| 天美传媒精品一区二区| 久久久a久久爽久久v久久| 国产三级在线视频| 91在线精品国自产拍蜜月| 国内久久婷婷六月综合欲色啪| 欧美日韩一区二区视频在线观看视频在线 | 日韩一本色道免费dvd| 一区二区三区四区激情视频 | 亚洲在线自拍视频| 亚洲美女黄片视频| 老司机影院成人| 人人妻人人看人人澡| 国产成年人精品一区二区| 久久精品国产亚洲av天美| 人人妻人人澡欧美一区二区| 亚洲aⅴ乱码一区二区在线播放| 精品久久久久久久久久免费视频| 亚洲精品亚洲一区二区| 69av精品久久久久久| 十八禁国产超污无遮挡网站| 亚洲内射少妇av| 精品无人区乱码1区二区| 一个人看视频在线观看www免费| 成人国产麻豆网| 免费观看的影片在线观看| 真人做人爱边吃奶动态| 午夜福利18| 免费搜索国产男女视频| 国产欧美日韩精品一区二区| 久久久午夜欧美精品| 岛国在线免费视频观看| 久久久久精品国产欧美久久久| 国产一级毛片七仙女欲春2| 九九热线精品视视频播放| 国产精品伦人一区二区| 免费无遮挡裸体视频| 九色成人免费人妻av| 国产单亲对白刺激| 在线观看免费视频日本深夜| 我要看日韩黄色一级片| 亚洲精品乱码久久久v下载方式| 人人妻,人人澡人人爽秒播| 国产亚洲精品久久久com| 久久久久免费精品人妻一区二区| 国产高清激情床上av| 精品人妻熟女av久视频| 波野结衣二区三区在线| 精品久久久久久久久av| 久久草成人影院| 夜夜看夜夜爽夜夜摸| 亚洲最大成人av| 精品人妻一区二区三区麻豆 | 亚洲国产欧洲综合997久久,| 欧美另类亚洲清纯唯美| 嫩草影视91久久| 18禁在线播放成人免费| 丝袜喷水一区| 亚洲真实伦在线观看| 国产精品久久视频播放| 大又大粗又爽又黄少妇毛片口| 99热精品在线国产| 国产精品永久免费网站| 国产成人精品久久久久久| 国产精品女同一区二区软件| 久久久a久久爽久久v久久| 欧美日本亚洲视频在线播放| 给我免费播放毛片高清在线观看| 91在线观看av| 国产私拍福利视频在线观看| 夜夜夜夜夜久久久久| 国产老妇女一区| av在线天堂中文字幕| 老师上课跳d突然被开到最大视频| 此物有八面人人有两片| 美女内射精品一级片tv| 免费一级毛片在线播放高清视频| 欧美bdsm另类| 国产成人aa在线观看| 网址你懂的国产日韩在线| 一级毛片久久久久久久久女| 精品久久久久久久久久免费视频| 日产精品乱码卡一卡2卡三| 最近手机中文字幕大全| 91狼人影院| 波多野结衣高清作品| 国产av一区在线观看免费| 国产不卡一卡二| 三级毛片av免费| 九九在线视频观看精品| 97热精品久久久久久| 国产精品一区二区三区四区免费观看 | 91狼人影院| 日本精品一区二区三区蜜桃| 内地一区二区视频在线| 无遮挡黄片免费观看| 美女 人体艺术 gogo| 国产高清不卡午夜福利| 黄色视频,在线免费观看| 国产aⅴ精品一区二区三区波| 少妇人妻一区二区三区视频| 在线免费观看不下载黄p国产| 人妻夜夜爽99麻豆av| 嫩草影院入口| 特级一级黄色大片| 国产黄a三级三级三级人| 99久久成人亚洲精品观看| av在线观看视频网站免费| 日韩亚洲欧美综合| .国产精品久久| av卡一久久| 国内精品一区二区在线观看| 免费大片18禁| 国产高清三级在线| 国产伦精品一区二区三区视频9| 成年免费大片在线观看| 国产极品精品免费视频能看的| 国模一区二区三区四区视频| 看免费成人av毛片| 校园人妻丝袜中文字幕| 亚洲性久久影院| 十八禁网站免费在线| 男女那种视频在线观看| av免费在线看不卡| 中文字幕免费在线视频6| 国产一区二区三区av在线 | videossex国产| 三级经典国产精品| 免费大片18禁| 亚洲中文日韩欧美视频| or卡值多少钱| 免费人成在线观看视频色| 菩萨蛮人人尽说江南好唐韦庄 | 日韩国内少妇激情av| a级毛片a级免费在线| 免费一级毛片在线播放高清视频| 我的女老师完整版在线观看| 亚洲国产精品国产精品| 久久99热这里只有精品18| 免费黄网站久久成人精品| 婷婷色综合大香蕉| 国产中年淑女户外野战色| 日韩人妻高清精品专区| 51国产日韩欧美| 亚洲电影在线观看av| 国产亚洲av嫩草精品影院| 老司机影院成人| 亚洲成人av在线免费| 精品久久久久久久末码| 久久久久国产网址| 欧美日韩乱码在线| 日韩成人伦理影院| 国产精品三级大全| 亚洲av成人精品一区久久| 色综合亚洲欧美另类图片| 中出人妻视频一区二区| 久久人人精品亚洲av| 午夜福利视频1000在线观看| 婷婷精品国产亚洲av| 卡戴珊不雅视频在线播放| 免费无遮挡裸体视频| 成人无遮挡网站| 午夜福利在线观看免费完整高清在 | 日韩成人伦理影院| 乱人视频在线观看| 国产精品不卡视频一区二区| 99在线人妻在线中文字幕| 国产色爽女视频免费观看| 99久久精品热视频| 亚洲av电影不卡..在线观看| 日本五十路高清| 人妻少妇偷人精品九色| 成人av一区二区三区在线看| 99riav亚洲国产免费| 亚洲成人av在线免费| 亚洲高清免费不卡视频| 九九在线视频观看精品| 久久九九热精品免费| 国内精品美女久久久久久| 1024手机看黄色片| 12—13女人毛片做爰片一| 国产亚洲av嫩草精品影院| 国产精品久久久久久久电影| 男插女下体视频免费在线播放| 亚洲精品一卡2卡三卡4卡5卡| 亚洲五月天丁香| 婷婷精品国产亚洲av| 日韩欧美精品免费久久| 精品福利观看| 毛片一级片免费看久久久久| 小说图片视频综合网站| 男人的好看免费观看在线视频| 日本黄大片高清| 99热全是精品| 午夜福利18| 久久精品国产99精品国产亚洲性色| 91av网一区二区| 国产单亲对白刺激| 看非洲黑人一级黄片| 欧洲精品卡2卡3卡4卡5卡区| 久久这里只有精品中国| 色尼玛亚洲综合影院| 菩萨蛮人人尽说江南好唐韦庄 | 成人高潮视频无遮挡免费网站| 综合色丁香网| 久久精品国产亚洲av涩爱 | 久久综合国产亚洲精品| 亚洲欧美成人精品一区二区| 深爱激情五月婷婷| 中文资源天堂在线| 欧美性猛交╳xxx乱大交人| 91在线观看av| 超碰av人人做人人爽久久| 精品少妇黑人巨大在线播放 | 熟妇人妻久久中文字幕3abv| 日本免费a在线| 又黄又爽又免费观看的视频| 99国产精品一区二区蜜桃av| 少妇高潮的动态图| 亚洲婷婷狠狠爱综合网| 99视频精品全部免费 在线| 高清毛片免费观看视频网站| 欧美最新免费一区二区三区| 99riav亚洲国产免费| 日本黄大片高清| 精品福利观看| 亚洲欧美清纯卡通| 精品一区二区三区人妻视频| av女优亚洲男人天堂| 日韩精品青青久久久久久| 夜夜夜夜夜久久久久| 五月伊人婷婷丁香|