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

    Deforestation and fragmentation trends of seasonal dry tropical forest in Ecuador:impact on conservation

    2021-10-12 08:11:52CarlosRivasJosGuerreroCasadoandRafaelNavarroCerillo
    Forest Ecosystems 2021年3期

    Carlos A.Rivas,José Guerrero-Casado and Rafael M.Navarro-Cerillo

    Abstract Background:Fragmentation and deforestation are one of the greatest threats to forests,and these processes are of even more concern in the tropics,where the seasonal dry forest is possibly one of the most threatened ecosystems with the least remaining surface area.Methods: The deforestation and fragmentation patterns that had occurred in Ecuadorian seasonal dry forests between 1990 and 2018 were verified,while geographic information systems and land cover shapes provided by the Ecuadorian Ministry of the Environment were employed to classify and evaluate three types of seasonal dry forests:deciduous,semi-deciduous,and transition.The study area was tessellated into 10 km2 hexagons,in which six fragmentation parameters were measured:number of patches,mean patch size,median patch size,total edge,edge density and reticular fragmentation index (RFI).The RFI was also measured both outside and inside protected natural areas (unprotected,national protected areas and protected forest).Moreover,the areas with the best and worst conservation status,connectivity and risk of disappearance values were identified by means of a Getis-Ord Gi*statistical analysis.Results: The deforestation of seasonal dry forests affected 27.04% of the original surface area still remaining in 1990,with an annual deforestation rate of ?1.12% between 1990 and 2018.The RFI has increased by 11.61% as a result of the fact that small fragments of forest have tended to disappear,while the large fragments have been fragmented into smaller ones.The semi-deciduous forest had the highest levels of fragmentation in 2018.The three categories of protection had significantly different levels of fragmentation,with lower RFI values in national protected areas and greater values in protected forests.Conclusions:The seasonal dry forest is fragmenting,deforesting and disappearing in some areas.An increased protection and conservation of the Ecuadorian seasonal dry forest is,therefore,necessary owing to the fact that not all protection measures have been effective.

    Keywords:Deciduous forest,Semi-deciduous forest,Remnant forest,Patch isolation,Habitat loss,Protected areas

    Introduction

    The term ‘forest fragmentation’ refers to the spatial configuration and amount of treed-vegetation (Hermosilla et al.2018),a landscape-level process during which anthropogenic factors progressively subdivide forest tracts into (initially,but not necessarily ultimately) smaller,geometrically more complex and more isolated patches as a result of natural processes and land use activities(McGarigal and Marks 1995;Chakraborty et al.2017).This concept can refer to the entire process of forest loss and isolation or,more specifically,to changes in the spatial configuration of remnants of forest that are the result of deforestation (Fahrig 2003;Kupfer 2006).The fragmentation process involves changes in the composition,structure and function of the landscape,and occurs on a mosaic background of natural patches created by changing landforms and natural disturbances(McGarigal and Marks 1995;Asbjornsen et al.2004).

    At the landscape level,the most common effect of fragmentation is the formation of new edges or the modification of existing ones,which play a fundamental role in the structure and functioning of ecosystems (Forman and Godron 1989;Asbjornsen et al.2004).These changes can alter ecological functions related to biodiversity,the nutrient cycle and the hydrological cycle,and may even affect the microclimate of the area (Asbjornsen et al.2004;Taubert et al.2018).

    The increase in forest fragmentation is one of the main threats to natural tree populations in the tropics around the world (Trejo and Dirzo 2000;Fuchs et al.2003),where large areas of forests have been transformed into pastures and crops,thus creating a mosaic of agricultural areas and forests in which forests remain as small scattered patches (Asbjornsen et al.2004;Taubert et al.2018).The tropical forests in South America underwent a net loss of 2.6 million hectares in the 2010–2020 period,although the deforestation rate has decreased significantly when compared to 2000–2010(FAO 2020a).More specifically,Ecuador maintained the highest deforestation rates in South America during the periods 1990–2010,with annual rates of between ?1.5%to ?1.8% (FAO 2011) and with an overall deforestation of 21,340 km2between 1990 and 2020 (FAO 2020a).One consequence of this intensive fragmentation is that 47 ecosystems of mainland Ecuador have been classified as very-highly or highly fragmented,i.e.30% of the natural areas (Ministerio del ambiente de Ecuador 2015).Those most affected are located in the coastal region,in which there was an area of annual deforestation of 678.13 km2between 1990 and 2008 (Sierra 2013).The deforestation and degradation of the seasonal dry forests in this region have been particularly intense,thus making them the most threatened type of forest in the country,in addition to being less protected than the evergreen forests (Manchego et al.2018;Rivas et al.2020).Deforestation has,in fact,become the greatest threat to seasonal dry forest ecosystems in Ecuador,with an average change in area reduction of 1.4% per year between 2008 and 2014 (Tapia-Armijos et al.2015;Manchego et al.2018).Indeed,tropical dry forests are among the most threatened ecosystems in the world (Hoekstra et al.2005;Portillo-Quintero and Sánchez-Azofeifa 2010),and are the ecosystems of which the least amount of original surface remains (less than 25%) (Ferrer-Paris et al.2018).This deforestation has,according to the IUCN criteria,led the equatorial dry forest to be classified as in critical danger of extinction (Ferrer-Paris et al.2018),and approximately 70% of the remaining surface has very high levels of fragmentation (Rivas et al.2020).Intense deforestation is consequently considered to be the main threat to the biodiversity of the tropical seasonal dry forests of the Tumbension region,which are characterised by a high degree of endemism since they harbour 16 endemic mammals (Loaiza 2013) and 39 endemic bird species (Bird Life International 2019).The Ecuadorian seasonal dry forests contain high levels of floristic diversity,and approximately 80% of their components are regionally endemic as part of the Tumbesian Endemism Centre (Ministerio del Ambiente del Ecuador 2012).Seasonal dry forest areas also provide local communities with wood and food products,which results in the degradation of the structure,functionality and dynamics of the forest (Ministerio del Ambiente del Ecuador 2012).

    Despite the worrying state of conservation,tropical seasonal dry forests have traditionally been studied to a lesser degree than their neighbours,humid forests,with a ratio of approximately one study in dry forests to six in humid forests (Lessmann et al.2014).One issue that has not been addressed in any great depth is the fragmentation of the Ecuadorian tropical dry forest in the last few decades,and how this fragmentation has transformed the landscape.According to the framework of the Convention on Biological Diversity (Biodiversity Indicators Partnership 2011),international organisations worldwide,such as the EBONE“Europe Biodiversity Observation Nature”(Parr et al.2010) or the BIP“Biodiversity Indicators Partnership”,have recommended analysing ecosystems through the use of fragmentation indices.Class indices separately quantify the quantity and distribution of each type of patch in the landscape,and fragmentation indices can,therefore,be considered for each type of patch (McGarigal and Marks 1995).The objective of this study was consequently to assess the fragmentation of the Ecuadorian seasonal dry forest between 1990 and 2018.The specific objectives of this work were the following:i) to study the deforestation and fragmentation of Ecuadorian seasonal dry forests during five different periods (1990,2000,2008,2014,2016 and 2018);ii) to describe the spatial patterns of fragmentation during these study periods;iii) to analyse different parameters of fragmentation (e.g.edge density,number of patches,mean patch size) in the three types of dry forests (deciduous,semi-deciduous and transition forest) in the region between the years 1990–2018;iv) to analyse fragmentation in order to find patterns that indicate the most vulnerable areas;and v) to compare the fragmentation index in protected and unprotected areas.The intention of this was to provide useful information on the state of the Ecuadorian dry forest and the areas with the worst conditions and conservation,which would be useful as regards developing effective protection measures according to the present conservation status and future trends.

    Materials and methods

    Study area

    Our study area included the seasonal dry forest in the coastal region of Ecuador (Fig.1a),also known as Western Ecuador,located along the Pacific Ocean and the west slope of the Andes mountain range.The coastal region is characterised by three large structural elements that influence the distribution patterns of the biota:the Guayas River,the Esmeraldas River and the Coastal mountain range.This region has a total of 24 ecosystems,22 of which are divided into two biogeographic regions that are clearly distinguishable as regards their composition and floristic structure,in addition to their bioclimate:the predominantly humid region of Chocó and the region of the Equatorial Pacific,which are mostly dry (Ministerio del Ambiente del Ecuador 2013).Seasonal dry forests in Ecuador thrive in extreme climatic conditions,with an annual rainfall of 400–600 mm in a period of 3–4 months,generally in February,March and April;the average annual temperature is 24.9°C,and the potential evapotranspiration is 1783 mm·year?1(Ministerio del ambiente del Ecuador 2012).In the present study,we considered the seasonal dry forest of the Ecuadorian Pacific,which is divided into deciduous and semi-deciduous areas.In deciduous forests,75% of individuals of the arboreal or shrub species lose their leaves during the dry period,which lasts between 6 and 8 months,whereas in the semi-deciduous forest,between 75% and 25% of individuals of the arboreal or shrub species lose their leaves and are located in areas in which the dry periods last between 1 to 6 months a year(Prentice 1990;Ministerio del ambiente del Ecuador 2013;Rivas et al.2020).

    Fig.1 a Map of continental Ecuador showing its three main geographical regions;b Division of the dry forest by phenologies in hexagons of 10 km2;c Protected areas

    GIS sources

    In order to limit the potential extent of the seasonal dry forests,the layers of phenology and land use were obtained from the Ecuadorian Ministry of the Environment(available at http://ide.ambiente.gob.ec/mapainteractivo).Land uses have been obtained by the Ecuadorian Ministry of the Environment,using Landsat satellite images and Advanced Spaceborne Thermal Emission and Reflection Radiometer (ASTER),orthorectification and these have later been certified by experts and by means of fieldwork,with a pixel size of 30 m (Peralvo and Delgado 2010;Ministerio del Ambiente 2012;MAE and MAGAP 2015;Ministerio del Ambiente 2017).The Kappa index is approximately 0.7 (Ministerio del Ambiente 2012).We selected those land-uses classified as native forests in zones with a deciduous and semideciduous phenology,since seasonal dry forest predominates in these areas.Flooded areas (mangrove areas) and areas without vegetation cover or without woody vegetation were eliminated.Once the area of seasonal dry forest had been obtained,a transition zone was created between the deciduous and semi-deciduous forest by applying a 10 km buffer in the area that divided both ecosystems.We eventually obtained three analysis zones:deciduous,semi-deciduous and transition forest(Fig.1b).Land-uses are available for the years 1990,2000,2008,2014,2016 and 2018,and they were reclassified into two main land-uses:native forest and non-forest zones.Shrub and Herbaceous Vegetation,Agricultural Land,Body of Water,Anthropic Zone,Other Land,No Information and Forest Plantations were classified as nonforest areas,while the ‘native forest’ land use was classified as forest (Ministerio del Ambiente 2017).

    Deforestation and fragmentation analysis

    We calculated the changes in these two main land-uses throughout the temporal periods.The deforestation rate was calculated by employing the formula proposed by Puyravaud (2003) (Eq.1).

    where A1and A2are the forest cover at times t1and t2,respectively.

    The study area was divided into 10 km2tiles(Fig.1b) made of hexagonal polygons,since this is considered the most suitable geometry when studying interaction and connectivity (Birch et al.2007).The use of polygons improves the ability to assess landscape metrics in a more homogeneous manner.We selected 10 km2because 99.8% of the world’s forest fragments cover less than 10 km2(FAO and PNUMA 2020).One of the three types of forests (deciduous,semi-deciduous,and transition) was assigned to each tile on the basis of the predominant type within each segment.The Patch Analysis Tool (Rempel et al.2012) in ArcGIS was employed in order to calculate different landscape metrics for each tile (Table 1).These were:Number of patches (NumP),average patch size (MPS),median patch size (MedPS),total edge (TE),and edge density (ED).These parameters were then used to calculate the reticular fragmentation index (RFI) on the basis of the percentage without forest (PSB%) and the percentage of edge density(ED%) (Table 1),using the formula proposed by Leautaud Valenzuela (2014).A 1990 forest fragment of 0.2 ha was used as a reference value in order to determine 100% of the PSB% and ED% metrics.This size was used because smaller sizes distorted the calculation.The RFI was divided into five categories:very high (>80%),high (60%–80%),medium (40%–60%),low (40%–20%) and very low (<20%).An RFI of 100% was attributed to those tiles from which the native forest had disappeared.

    Table 1 Description of the fragmentation metric parameters analysed according to McGarigal and Marks (1995) and Leautaud Valenzuela (2014)

    Fragmentation patterns

    The fragmentation patterns were described by employing the Getis-Ord Gi* analysis (Ord and Getis 1995) for the years 1990 and 2018 and by considering the RFI values.The resulting Z-scores and p-values indicate a spatial cluster of high or low RFI values.At 5% significance (p ≤0.05),a Z-score greater than 1.96 indicates a hot spot,while a Z-score smaller than ?1.96 indicates a cold spot and the remaining values are classified as not significant (?1.96 0.05),thus suggesting a random spatial process (Feng et al.2018).A transition matrix was created using the categories of the Getis-Ord Gi* analysis for the years 1990 and 2018 to identify the probability of a hexagon disappearing or of its state changing,based on its initial state (1990) upon its categorisation.

    Fragmentation in protected and unprotected areas

    In order to test the trend of the RFI in protected and unprotected areas (Fig.1c),an RFI trend index was calculated as follows (Eq.2):

    We then assigned one of the following three protection categories to each tile:unprotected,protected by the Heritage of Natural Areas of the Ecuadorian State(PANE in Spanish),and Protected forests.The RFI trend index of these three categories was then compared in order to verify whether the degree of protection prevents fragmentation more effectively.

    Statistical analysis

    Wilcoxon paired tests were then used to compare the value of the fragmentation indicators (RFI,NumP,MPS,MedPS,TE,and DE),which were considered as dependent variables,between 1990 and 2018 in each of the three types of forest (deciduous,semi-deciduous and transition).The tiles were employed to pair these tests in order to consider the variations among the same grids between the two periods.A Kruskal-Wallis test was used to compare the RFI values (dependent variable) obtained for the three types of forest in the year 2018.The same type of test was similarly used to compare the RFI trend indices (dependent variable) obtained for the three categories of protection (unprotected,PANE,and Protected forests).In both the Kruskal-Wallis tests,pairwise comparisons (post hoc) were developed in order to verify the differences among the levels of the independent variables.The tiles were the experimental units in all the statistical analyses.InfoStat software was used in all the statistical analyses.

    Results

    Evolution of loss and fragmentation of forest

    Since 1990,2631.91 km2has been lost (Table 2),which signifies a loss of 27.04% of the original surface still remaining in 1990,and an annual deforestation rate of 94 km2(?1.12%).

    This deforestation has changed the degree of fragmentation of the dry forest (Table 3 and Fig.2),since the mean and median RFI values have increased,particularly from 1990 to 2008 (Table 3).In 1990,42.28% of the tiles of seasonal dry forest were classified with a low or very low RFI,while this figure dropped to 29.15% in 2018.Moreover,432 tiles that had some forest patches in 1990 had no patches in 2018.

    Fig.2 Spatial evolution of the RFI of the equatorial dry forest per tile in each of the six time periods

    Spatial evolution of fragmentation

    Deforestation has occurred principally in the north of the study area,which formally contained small forest fragments that have disappeared or been considerably reduced since 1990 (Figs.2 and 3).Other affected areas were located in the Guayas areas and in the centralsouth,where many of the segments had disappeared,leading to a significant increase in RFI.The edge areas of the large forest fragments have been deforested,as has also occurred with the small fragments,which has resulted in the disappearance of those forests throughout the territory analysed (Fig.3).

    Fig.3 Comparison of the central north and central south of the Ecuadorian coast,showing details of forest fragments(green areas)in the years 1990 and 2018

    Fragmentation indexes for the three types of forests

    According with the Wilcoxon paired tests,the RFI value was significantly higher in 2018 than in 1990 for the three types of forests (Table 4),and the NumP,MPS,MedPS and TE values were significantly higher in 1990 than in 2018.The ED was significantly higher in 2018 than in 1990 for the semi-deciduous and transition forests,whereas it was lower in 2018 for deciduous forests(Table 4).The semi-deciduous forest attained the highest increase in RFI from 1990 to 2018 (highest mean difference),with the highest levels of fragmentation occurring in 2018 (Table 4).With regard to the numberof patches (NumP) within the tile,this has not undergone a great variation as regards either the total or forest types.

    Table 2 Surface of equatorial dry native forest and other land uses in each of the periods included in this study

    Table 3 Variation of the mean and median RFI values (%) per tile in each of the 6 periods.SD=standard deviation

    The Kruskal-Wallis test showed differences among the RFI values obtained for the three types of seasonal dry forests in 2018 (H=295.65;p<0.0001),with the highest value being attained for semi-deciduous forest and the lowest for the transition forest(Table 5).

    Fragmentation patterns

    The Getis-Ord Gi* analysis shows the hot and cold fragmentation spots (Fig.4).These results indicate that the hot areas,which had a worse structural connectivity,were more vulnerable to disappearance and had a worse state of conservation.The comparison of 1990 with 2018 highlights this evolution (Fig.4 and Table 6).In 2018,there were 981 (29.70%) tiles catalogued as cold spots,1044 (31.60%) with no significant differences,and 843(25.52%) as hotspots;in 1990,meanwhile,there were 1063 (32.18%),1286 (38.93%) and 863 (26.12%) respectively.The transition matrix (Table 6) shows that of the 863 tiles classified as hotspots in 1990,213 (24.68%) disappeared during the studied period,and 572 (66.28%)remained in the hotspot category.

    Table 4 Variation between 1990 and 2018 as regards the values obtained for the different fragmentation indicators by forest type,showing the Z and p-value from the Wilcoxon paired test.Mean-dif=value 2018–value 1990;SD=standard deviation

    Table 5 RFI values in 2018 for the three types of forests.N=number of tiles;SE=standard error of means.Lower case letters indicate significant differences according to the post hoc test

    Table 6 Transition matrix according to the categorisation provided by the Getis-Ord Gi*analysis,showing the number of tiles whose state changed from 1990 (columns) to 2018 (rows).NS=not significant changes

    Fragmentation in protected and unprotected areas

    Of the 2707 tiles into which the seasonal dry forest was divided,only 7.24% was covered by PANE and 8.32% by protected forests,while 84.45% were unprotected.The RFI trend was significantly different for the three protection categories (H=19.60;p<0.001),with the lowest value for PANE (0.3±0.07),the highest values for Protected forests (0.97±0.17),and the intermediate values for unprotected areas (0.42±0.03).

    Discussion

    Ecuador is undergoing a high rate of deforestation,and the seasonal dry forest is no exception (Sierra 2013).Our results show that 2631.91 km2of seasonal dry forest have been converted to other land uses in the last three decades (87% of the forest that was deforested between 1990 and 2018 had been transformed into agricultural land by 2018,while 7% had been transformed into scrubland),with the extinction of many patches,thus causing a constant increase in fragmentation.This fragmentation has occurred throughout the study area,although we have identified areas with higher fragmentation values and that are spatially aggregated(hotspots),thus suggesting that an important area of the remaining forests runs a high risk of disappearing in the next few years.All these data suggest the urgency of implementing effective conservation measures to preserve the remaining Ecuadorian seasonal dry forest patches and promoting connectivity,with the eventual goal of preventing the disappearance of new areas and ensuring the functional ecology of the remaining forests.

    Deforestation of Ecuadorian seasonal dry forests

    We observed a dramatic level of deforestation of native forests,and consequently assume that this is a threat to the flora and fauna that inhabit these forests.According to our results,the Ecuadorian seasonal dry forests underwent a net loss of 27% from 1990 to 2018,signifying an annual deforestation rate of ?1.12%.This annual deforestation rate was higher than the rates found in other Latin American countries (Brazil:?0.56%;Colombia:?0.31%;Peru:?0.18%),but lower than that of Paraguay (?1.53%) (FAO 2020b).When compared to other dry forests in the region,the deforestation rate was in the same range as that of Paraguay and Chile (between ?1% and ?2%),with lower rates than those found for Argentina and Mexico (>?2%),but greater than those found for Brazil,Costa Rica and Venezuela (

    Fragmentation

    According to our calculations,all the landscape metrics attained worse fragmentation values in 2018 than in 1990 for all three types of forest.The number of patches decreased and the forests had a smaller mean patch size,which led to an overall increase in the fragmentation index (RFI).But if this information is analysed together with the other fragmentation metrics,it will be noted that this is associated with the disappearance of the smaller patches and the fragmentation of large patches,which has kept the number of patches constant,but has increased the fragmentation.The edge density (ED) in the semi-deciduous and transitional forest has probably increased as a consequence of the forest fragments getting smaller (lower mean MPS) and the increase in the number of small fragments (median MPS).However,the ED in the deciduous forest is probably decreasing because the small fragments are disappearing and the largest ones are becoming smaller (thus keeping the MPS constant),and since smaller fragments had higher edge density values,their disappearance may have led to a decrease in the ED value (Hargis et al.1998).This process makes this measure less sensitive because,although these small patches disappear,the landscape fragmentation increases (Whelan and Maina 2005;Tulloch et al.2016).Small patches have been shown to be fundamental to ecosystems,particularly in those that are highly fragmented,and their disappearance may,therefore,have negative consequences for them (Tulloch et al.2016).

    Of the three types of forest considered,the semideciduous forest was the most fragmented.For instance,the average patch size (MPS) of the semi-deciduous forest attained very worrying values (0.58 km2),considering that the tile area is 10 km2.This may be owing to the fact that the more humid forests are more fertile and are,therefore,more prone to the establishment of crops and pastures (Ministerio del ambiente del Ecuador 2012;Lessmann et al.2014).Moreover,many areas of seasonal dry forest have degenerated into savannah,scrub or grasslands owing to the high pressure of livestock and overgrazing,which could be the cause of the disappearance of the small fragments of deciduous forest,thus limiting forest growth and extension (Trejo and Dirzo 2000;Sales et al.2020).Conversely,drier areas are often perceived as areas that are poorer in resources (Siyum 2020),which could explain the lower conversion of the deciduous forest when compared to the semi-deciduous forest.

    Connectivity

    Upon comparing the images from 1990 and 2018(Figs.2,3 and 4),it will be noted that a quarter of the forest fragments classified as hotspots in 1990 had completely disappeared by 2018,which indicates that these areas are more prone to disappearance.Many factors may lead to differences in deforestation among areas,such as the growth rates of the localhuman population,the presence of particular hardwood species,the development of specific types of agriculture,the distance to roads and trails,the distance to rivers or the suitability of the land (e.g.soil features or being steep) for agricultural purposes in general (Andam et al.2008;Barber et al.2014;van Der Hoek 2017).Future studies should,therefore,be carried out to evaluate which factors explain a greater or lesser fragmentation of the landscape in order to identify those forests that are still wellconserved and run the greatest risk of becoming fragmented.

    Deforestation for agricultural and livestock purposes has been identified as one of the main reasons for the loss of seasonal dry forests in Ecuador (Tapia-Armijos et al.2015;Prieto-Torres et al.2018),and this also occurs in other countries,such as in the Brazilian Cerrado(Trigueiro et al.2020) or in the Mexican Yucatan (Smith et al.2019),and in other dry forests in Latin America(Armenteras et al.2017).The spatial analysis of fragmentation indicates two large areas of high concentrations of fragmentation (hotspots) in Ecuadorian seasonal dry forests (Fig.4).The first area is located in the province of Manabí,which is the province with the highest agricultural production,and in which 777,088 ha correspond to cultivated and natural pastures,contributing more than 20% of the country’s agricultural area (INEC 2019).The second fragmented area corresponds to the urban areas of Guayaquil and Machala,the first and third largest cities as regards human population,respectively (Instituto nacional de estadisticas y censos 2010).The population of Ecuador has increased dramatically in the last few decades,since it has grown by 452% in 60 years (Villacís and Carrillo 2012),reaching 17 million inhabitants in 2019 (Instituto Nacional de Estadística 2019).This has led to an increase in demands for food and an increase in the areas devoted to agricultural and livestock production,which are the greatest threats to tropical forests in South America and Africa (Laurance et al.2014).The dry forests are used by the local population,since they have environmental,social and economic importance for various segments of the rural communities (Brice?o et al.2016).The quality of wood products has historically led to interventions in these forests in order to extract wood and food products as a means of subsistence.The seasonal dry forest provides essential ecosystem goods and services,livelihoods and is vital to the well-being of its residents,since it provides supplies(water,wood,food,biofuels) (Nelson et al.2020;Siyum 2020).Population growth consequently also increases the pressure on remnants of forest.In summary,the great importance of agriculture for the economy of the region,together with the growth of the human population,have increased the conversion of natural forests into agricultural land,thus leading to a rise in fragmentation.

    Fig.4 Comparison of the hot and cold spots by means of Getis-Ord Gi*analysis in the years 1990 and 2018

    The analyses of hotspots showed the areas with a worst conservation status,low connectivity and high fragmentation,and these may be priority areas for forest restoration and an increase in connectivity.Furthermore,in areas identified as cold spots,the actions should be focused principally on preventing deforestation.Fragmentation can have negative consequences for populations of wild species that inhabit the dry forest(Solórzano et al.2021),since many remaining patches are becoming isolated and exposed to disappearance(Margules and Pressey 2000).The synergistic effects of fragmentation lead to changes in climate,which can,in turn,change the structure of the vegetation,soil cover and nutrient status,thus affecting the species that inhabit these forest fragments (Margules and Pressey 2000).Changes take place in these isolated fragments,which can lead to the collapse of populations (Laurance et al.2012).

    Conservation implications

    Tropical seasonal dry forests are the ecosystems with the least remaining surface in Ecuador (Ferrer-Paris et al.2018) and are possibly the most threatened tropical ecosystems in the world (Escribano-Avila et al.2017).They are considered an endangered ecosystem owing to the high degree of endemism and species richness;however,less than 10% of their area is protected (Prieto-Torres et al.2018).Protected areas are important for conservation (Barber et al.2014;van Der Hoek 2017),and should be expanded in the case of the Ecuadorian dry forest(see below),which is less protected than other ecosystems (Rivas et al.2020).

    Previous works have shown that the dry ecosystems in the Coastal Region of Ecuador are underrepresented in the PANE (Sierra et al.2002;Lessmann et al.2014;Escribano-Avila et al.2017),thus suggesting that it is necessary to create new protected areas in order to preserve these ecosystems (Lessmann et al.2014;Cuesta et al.2017).As our results show,there has been less fragmentation in the PANE,while it has increased more in unprotected areas and has been particularly dramatic in protected forests.The protected areas included in the PANE have,therefore,been partly effective as regards preventing deforestation,with a smaller increase in the RFI value than in unprotected areas from 1990 to 2018.These results coincide with those of two previous works,which demonstrated that the deforestation rates were lower inside protected areas,although deforestation still took place in those areas (van Der Hoek 2017;Ford et al.2020).Protected forests are not,however,an effective conservation tool for the conservation of seasonal dry forests since,according our results,the RFI increased even more than in unprotected areas.Indeed,more than half of the areas in protected forests were classified as non-forest land use (Rivas et al.2020).Although intensive agriculture and deforestation is prohibited in those forests (Sandoval et al.2017),our results showed that the RFI dramatically increased inside the protected forests,signifying that the current management system of these forests needs to be reviewed with the aim of ensuring their intended conservation goals.Protection measures should,therefore,be implemented,and they should be established in areas of high priority,which would reduce fragmentation and increase structural connectivity.Several scientists have recently highlighted the importance of small forest patches in fragmented landscapes as regards biodiversity conservation (Tulloch et al.2016;Fahrig et al.2019;Volenec and Dobson 2020),thus suggesting that it is necessary to maximize the total amount of habitat conserved,irrespective of its rate of fragmentation (Fahrig et al.2019;Ríos et al.2021).Specific conservation measures,such as the creation of small reserves (including private protected areas:Guerrero-Casado et al.2021),should,therefore,be implemented in order to protect the few remnants of seasonal dry forest.

    Conclusion

    According to our results,the Equatorial seasonal dry forest has undergone a continuous process of deforestation that has led to the loss of more than 2600 km2of native dry forest in the last three decades,which is causing an increase in fragmentation,with semi-deciduous forest being the most affected.Fragmentation has increased since 1990,and the number of patches has decreased as a result of the reduction in the forest area,thus increasing the border and patching the forest in isolated fragments,and consequently making it more exposed to anthropic and natural changes.Fragmentation occurs throughout the entire distribution area of seasonal dry forest,which degrades the ecosystem,increases its vulnerability,reduces the area and decreases its connectivity,thus leading to high values of biodiversity loss.Our results show that many areas of seasonal dry forests run a great risk of disappearing if effective protection is not provided or conservation measures are not taken,and it is,therefore,urgent to establish conservation measures that will avoid the continued fragmentation of these forests.

    Supplementary Information

    The online version contains supplementary material available at https://doi.org/10.1186/s40663-021-00329-5.

    Additional file 1:Figure S1.Percentage of RFI per tile between the years 1990 and 2018 for the following values:disappeared,very high,high,medium,low and very high.

    Acknowledgements

    José Guerrero-Casado is currently supported by the European Regional Development Fund (ERDF) and the Consejería de Transformación Económica,Industria,Conocimiento y Universidades (project reference:1264483-R).Rafael M Navarro Cerrillo is particularly grateful for the support of the ISOPINE(UCO-1265298) and ESPECTRAMED (CGL2017-86161-R) projects.We acknowledge the institutional support of the University of Cordoba-Campus de Excelencia CEIA3.

    Authors’ contributions

    Carlos A.Rivas:data collection,experimental design and writing the original draft.José Guerrero-Casado:statistical analysis,supervision,review and editing the original draft.Rafael M.Navarro-Cerillo:experimental design,conceptualisation,supervision,review and editing the original draft.The author(s)read and approved the final manuscript.

    Funding

    José Guerrero-Casado is currently supported by the European Regional Development Fund (ERDF) and the Consejería de Transformación Económica,Industria,Conocimiento y Universidades(project reference:1264483-R).Rafael M Navarro Cerrillo is particularly grateful for the support of the ISOPINE(UCO-1265298) and ESPECTRAMED (CGL2017–86161-R) projects.

    Availability of data and materials

    The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.

    Declarations

    Ethics approval and consent to participate

    Not applicable.

    Consent for publication

    Not applicable.

    Competing interests

    The authors declare that they have no competing interests.

    Author details

    1Instituto de Ciencias Básicas,Universidad Técnica de Manabí,Portoviejo,Manabí,Ecuador.2Department of Forest Engineering,Laboratory of Dendrochronology,Silviculture and Global Change–DendrodatLab–ERSAF,University of Cordoba,Campus de Rabanales,Crta.IV,km.396,E-14071 Cordoba,Spain.3Departamento de Medicina Veterinaria,Facultad de Ciencias Veterinarias,Universidad Técnica de Manabí,Portoviejo,Manabí,Ecuador.4Department of Zoology,University of Cordoba,Campus de Rabanales,14071 Cordoba,Spain.

    Received:31 January 2021Accepted:14 June 2021

    插阴视频在线观看视频| 亚洲精品乱码久久久久久按摩| 国产精品久久电影中文字幕| 麻豆精品久久久久久蜜桃| 在现免费观看毛片| 亚洲四区av| 啦啦啦啦在线视频资源| 级片在线观看| 免费一级毛片在线播放高清视频| 国产黄片美女视频| 91av网一区二区| 亚洲欧洲日产国产| 精品人妻熟女av久视频| 国产在线男女| 久久精品人妻少妇| 丰满乱子伦码专区| 亚洲国产色片| 国内精品久久久久精免费| 网址你懂的国产日韩在线| 久久99精品国语久久久| 黄色视频,在线免费观看| 亚洲欧美日韩高清专用| 啦啦啦观看免费观看视频高清| 人人妻人人澡人人爽人人夜夜 | 少妇人妻精品综合一区二区 | 国产精品无大码| 99在线人妻在线中文字幕| 一区二区三区免费毛片| 久久久久久久久久久丰满| 特级一级黄色大片| 国产真实伦视频高清在线观看| 长腿黑丝高跟| 国产精品久久久久久精品电影小说 | 99久久无色码亚洲精品果冻| 九九爱精品视频在线观看| 日韩成人伦理影院| 黄片无遮挡物在线观看| 久久久精品欧美日韩精品| 久久久久网色| 啦啦啦韩国在线观看视频| 99视频精品全部免费 在线| 欧美日韩精品成人综合77777| 听说在线观看完整版免费高清| 国产精品久久久久久精品电影小说 | 亚洲天堂国产精品一区在线| 波野结衣二区三区在线| 成人国产麻豆网| 综合色丁香网| 亚洲欧美精品综合久久99| kizo精华| 国产亚洲5aaaaa淫片| 亚洲国产欧美人成| 久久久久性生活片| 亚洲丝袜综合中文字幕| 日韩av不卡免费在线播放| 国产av不卡久久| 高清毛片免费观看视频网站| 国产乱人偷精品视频| 岛国毛片在线播放| 国产一区二区在线观看日韩| 亚洲婷婷狠狠爱综合网| 在线播放国产精品三级| 成人特级黄色片久久久久久久| 热99在线观看视频| 91在线精品国自产拍蜜月| a级毛色黄片| 看片在线看免费视频| 亚洲中文字幕一区二区三区有码在线看| 婷婷六月久久综合丁香| 国产精品久久久久久av不卡| 欧美高清性xxxxhd video| 内射极品少妇av片p| 亚洲激情五月婷婷啪啪| 久久久精品欧美日韩精品| 国产美女午夜福利| 在线播放无遮挡| 91午夜精品亚洲一区二区三区| 国产成人精品婷婷| 久久午夜亚洲精品久久| 欧美bdsm另类| 中国国产av一级| 午夜激情欧美在线| 女的被弄到高潮叫床怎么办| 天堂网av新在线| 亚洲国产欧洲综合997久久,| АⅤ资源中文在线天堂| 在线免费观看的www视频| 黑人高潮一二区| 精品久久久久久久久av| 观看美女的网站| 卡戴珊不雅视频在线播放| 亚洲婷婷狠狠爱综合网| 成人性生交大片免费视频hd| 欧美丝袜亚洲另类| 欧美3d第一页| 国产综合懂色| 日韩欧美在线乱码| av.在线天堂| 2022亚洲国产成人精品| 国产精品av视频在线免费观看| 我的女老师完整版在线观看| 亚洲国产色片| 国产精品无大码| 日韩三级伦理在线观看| 成人av在线播放网站| 99精品在免费线老司机午夜| 国产三级在线视频| 免费搜索国产男女视频| 亚洲在线观看片| 身体一侧抽搐| 国产一区二区激情短视频| 三级毛片av免费| 在线播放无遮挡| 91av网一区二区| 欧美zozozo另类| 精品久久久久久久久av| 午夜福利视频1000在线观看| 国产精品不卡视频一区二区| 国产一区二区三区av在线 | 人人妻人人看人人澡| 只有这里有精品99| 亚洲七黄色美女视频| 欧美成人免费av一区二区三区| 国产毛片a区久久久久| 少妇丰满av| 精品久久久久久久末码| 丝袜喷水一区| 亚洲最大成人中文| 欧洲精品卡2卡3卡4卡5卡区| 免费人成视频x8x8入口观看| 99九九线精品视频在线观看视频| 国产老妇伦熟女老妇高清| av黄色大香蕉| 亚洲成av人片在线播放无| 一区二区三区高清视频在线| 免费不卡的大黄色大毛片视频在线观看 | 一个人免费在线观看电影| 久久久久网色| 99热这里只有是精品50| 天堂中文最新版在线下载 | 国产成年人精品一区二区| 观看美女的网站| 夜夜夜夜夜久久久久| 国产成人一区二区在线| 国产 一区精品| 91久久精品国产一区二区三区| 中文字幕精品亚洲无线码一区| 啦啦啦观看免费观看视频高清| 国产熟女欧美一区二区| 亚洲av.av天堂| 春色校园在线视频观看| 99热这里只有是精品在线观看| 一级毛片aaaaaa免费看小| 两个人的视频大全免费| 亚洲av男天堂| 亚洲高清免费不卡视频| 国产v大片淫在线免费观看| 中文资源天堂在线| 久久精品人妻少妇| 日韩欧美三级三区| 精品久久久久久久末码| 午夜精品国产一区二区电影 | 日本一本二区三区精品| 熟妇人妻久久中文字幕3abv| 自拍偷自拍亚洲精品老妇| 99热网站在线观看| 国产精品爽爽va在线观看网站| 最后的刺客免费高清国语| 精品久久久久久久久久久久久| 亚洲经典国产精华液单| 男人舔女人下体高潮全视频| 久久韩国三级中文字幕| 国产中年淑女户外野战色| 精品国内亚洲2022精品成人| 久久久久国产网址| 村上凉子中文字幕在线| 一边亲一边摸免费视频| 久久九九热精品免费| 精品久久久久久久久久久久久| av卡一久久| 在线免费十八禁| 国产av一区在线观看免费| 欧美不卡视频在线免费观看| 日产精品乱码卡一卡2卡三| 丰满人妻一区二区三区视频av| 日本爱情动作片www.在线观看| 久久精品夜色国产| 国产欧美日韩精品一区二区| 全区人妻精品视频| 99久久中文字幕三级久久日本| 日本一二三区视频观看| 天堂网av新在线| 久久久久久九九精品二区国产| 亚洲一级一片aⅴ在线观看| 精品久久久久久久末码| 欧美色视频一区免费| 婷婷色av中文字幕| 午夜精品一区二区三区免费看| 国产精品.久久久| av在线观看视频网站免费| 亚洲内射少妇av| 亚洲精品国产av成人精品| 成人亚洲精品av一区二区| 国产成人一区二区在线| 日韩一本色道免费dvd| 99热只有精品国产| av视频在线观看入口| 久久精品国产亚洲网站| 十八禁国产超污无遮挡网站| 伊人久久精品亚洲午夜| 九九热线精品视视频播放| 亚洲最大成人av| 欧美成人a在线观看| 麻豆乱淫一区二区| 精品久久久噜噜| 精品不卡国产一区二区三区| 中文欧美无线码| 免费看美女性在线毛片视频| 1024手机看黄色片| 久久久精品大字幕| 在线天堂最新版资源| 最新中文字幕久久久久| 日韩成人av中文字幕在线观看| 三级毛片av免费| 成人欧美大片| 色哟哟哟哟哟哟| 国产私拍福利视频在线观看| 岛国在线免费视频观看| 久久99热6这里只有精品| 啦啦啦观看免费观看视频高清| 亚洲国产精品sss在线观看| 日本一二三区视频观看| 国语自产精品视频在线第100页| 国国产精品蜜臀av免费| 丝袜喷水一区| 草草在线视频免费看| 国产淫片久久久久久久久| 免费看光身美女| 国产精品一区二区三区四区久久| 亚洲精品456在线播放app| 1000部很黄的大片| 蜜桃久久精品国产亚洲av| 两个人的视频大全免费| 性色avwww在线观看| 国产成人一区二区在线| АⅤ资源中文在线天堂| 高清在线视频一区二区三区 | 99久久九九国产精品国产免费| 国产成人精品婷婷| 免费看av在线观看网站| 欧美一区二区精品小视频在线| 日日摸夜夜添夜夜爱| 国产精品野战在线观看| 日韩欧美三级三区| 99久久人妻综合| 国产精品久久久久久亚洲av鲁大| 午夜福利视频1000在线观看| 日韩在线高清观看一区二区三区| 成人三级黄色视频| 性色avwww在线观看| 国产综合懂色| 精品久久久久久久久久免费视频| 插阴视频在线观看视频| 内地一区二区视频在线| 不卡一级毛片| 一级毛片电影观看 | 在线观看美女被高潮喷水网站| 老司机福利观看| 欧美最新免费一区二区三区| 精品午夜福利在线看| 麻豆av噜噜一区二区三区| 九九在线视频观看精品| 一边亲一边摸免费视频| 干丝袜人妻中文字幕| 99视频精品全部免费 在线| 变态另类丝袜制服| 国产v大片淫在线免费观看| 欧美日韩一区二区视频在线观看视频在线 | 免费看光身美女| 人人妻人人澡欧美一区二区| 亚洲最大成人av| 日本撒尿小便嘘嘘汇集6| 亚洲,欧美,日韩| av免费观看日本| 亚洲七黄色美女视频| 久久精品国产99精品国产亚洲性色| 亚洲最大成人中文| 最新中文字幕久久久久| 18禁在线无遮挡免费观看视频| 一本精品99久久精品77| 精品一区二区三区人妻视频| 99在线人妻在线中文字幕| 简卡轻食公司| 人人妻人人看人人澡| 黄片无遮挡物在线观看| 午夜精品在线福利| 久久精品91蜜桃| 毛片一级片免费看久久久久| 两性午夜刺激爽爽歪歪视频在线观看| 深爱激情五月婷婷| 国产精品不卡视频一区二区| 日韩强制内射视频| 中文字幕av在线有码专区| 美女国产视频在线观看| 91精品国产九色| 18+在线观看网站| 国产成人91sexporn| 人妻久久中文字幕网| 男人舔奶头视频| 嫩草影院入口| 天天一区二区日本电影三级| 搡老妇女老女人老熟妇| 国产成人a区在线观看| 精品国产三级普通话版| 国产精品久久电影中文字幕| 久久精品久久久久久噜噜老黄 | 久久久欧美国产精品| 亚洲国产欧美在线一区| 国产成人a区在线观看| 国产成人91sexporn| 亚洲av电影不卡..在线观看| 精华霜和精华液先用哪个| 国产日本99.免费观看| 两个人的视频大全免费| 日韩欧美精品免费久久| 国产精品一区二区在线观看99 | 97人妻精品一区二区三区麻豆| 亚洲人与动物交配视频| 欧美一区二区精品小视频在线| 免费观看的影片在线观看| 三级经典国产精品| 日韩av在线大香蕉| 亚洲色图av天堂| 嫩草影院新地址| 尾随美女入室| 在线观看美女被高潮喷水网站| 国产成人午夜福利电影在线观看| 国产精品爽爽va在线观看网站| 免费人成在线观看视频色| 国产成年人精品一区二区| 天堂影院成人在线观看| 99热网站在线观看| 国产亚洲av片在线观看秒播厂 | h日本视频在线播放| 日韩成人av中文字幕在线观看| 国产精品不卡视频一区二区| 欧美性感艳星| 国产免费男女视频| 91在线精品国自产拍蜜月| 麻豆成人av视频| 91久久精品国产一区二区三区| 国产精品电影一区二区三区| 国产一区二区在线av高清观看| 亚洲三级黄色毛片| 亚洲国产欧美人成| 欧美高清成人免费视频www| 一级毛片电影观看 | 久久精品国产亚洲av涩爱 | 黄色日韩在线| 国模一区二区三区四区视频| 成人特级黄色片久久久久久久| 一边摸一边抽搐一进一小说| 九九在线视频观看精品| 给我免费播放毛片高清在线观看| 三级国产精品欧美在线观看| 麻豆成人午夜福利视频| 又粗又爽又猛毛片免费看| 亚洲国产精品成人久久小说 | 亚洲成a人片在线一区二区| 久久久午夜欧美精品| 午夜激情福利司机影院| 日本与韩国留学比较| 欧美xxxx性猛交bbbb| 国产成人精品久久久久久| 能在线免费看毛片的网站| 成人亚洲欧美一区二区av| 亚洲av熟女| 亚洲欧美日韩东京热| 18禁在线播放成人免费| 亚洲最大成人中文| 毛片一级片免费看久久久久| 伊人久久精品亚洲午夜| 夜夜爽天天搞| 熟女电影av网| 久久九九热精品免费| 国产成人精品久久久久久| 午夜精品在线福利| 国产乱人视频| 非洲黑人性xxxx精品又粗又长| 在线观看av片永久免费下载| 一级黄片播放器| 波多野结衣高清无吗| 久久人人爽人人片av| 午夜福利在线在线| 久久人人爽人人片av| 此物有八面人人有两片| 国产亚洲欧美98| 国产精品福利在线免费观看| 亚洲精品456在线播放app| www日本黄色视频网| 中文字幕免费在线视频6| 少妇被粗大猛烈的视频| 又爽又黄a免费视频| 国产成人91sexporn| 观看免费一级毛片| 床上黄色一级片| 日韩一区二区视频免费看| 国产成年人精品一区二区| 国产精品久久电影中文字幕| av天堂在线播放| 性欧美人与动物交配| 欧美日韩精品成人综合77777| 校园春色视频在线观看| 高清日韩中文字幕在线| 久久久精品大字幕| 久久久久久大精品| 欧美日本视频| 黄色视频,在线免费观看| 欧美+日韩+精品| 观看美女的网站| 黄色一级大片看看| 一本久久精品| 一区二区三区高清视频在线| a级毛色黄片| 亚洲精品日韩av片在线观看| 久久精品国产鲁丝片午夜精品| 国产成人91sexporn| 国产在线精品亚洲第一网站| 久久精品综合一区二区三区| 国产成人一区二区在线| 男女视频在线观看网站免费| 精品熟女少妇av免费看| 乱系列少妇在线播放| 欧美+亚洲+日韩+国产| 99热全是精品| 国产不卡一卡二| 亚洲不卡免费看| 久久精品国产清高在天天线| 黄色欧美视频在线观看| 舔av片在线| 熟妇人妻久久中文字幕3abv| 国产精品久久视频播放| 悠悠久久av| 日韩 亚洲 欧美在线| 波野结衣二区三区在线| 1024手机看黄色片| 亚洲真实伦在线观看| 99热全是精品| 26uuu在线亚洲综合色| 晚上一个人看的免费电影| 国产高清激情床上av| 国产高清三级在线| 日韩强制内射视频| 婷婷亚洲欧美| 在现免费观看毛片| 一夜夜www| 久久精品国产亚洲av涩爱 | 99视频精品全部免费 在线| 99久国产av精品国产电影| 中文字幕av成人在线电影| 国产精品一二三区在线看| 26uuu在线亚洲综合色| 老师上课跳d突然被开到最大视频| 亚洲av免费在线观看| 国产色婷婷99| 国产亚洲91精品色在线| 丝袜喷水一区| 亚洲精品色激情综合| 国产高清三级在线| 日本免费一区二区三区高清不卡| 久久精品人妻少妇| 国产成人a∨麻豆精品| 看黄色毛片网站| 国产精品久久久久久精品电影| 校园春色视频在线观看| 亚洲乱码一区二区免费版| 美女内射精品一级片tv| 亚洲无线观看免费| 日韩三级伦理在线观看| 中文在线观看免费www的网站| 亚洲成人精品中文字幕电影| 国产真实乱freesex| 国产免费一级a男人的天堂| 99riav亚洲国产免费| 国产精品一区二区在线观看99 | 99热这里只有精品一区| 欧美日韩精品成人综合77777| 国产精品一二三区在线看| 日本-黄色视频高清免费观看| 女人十人毛片免费观看3o分钟| 麻豆成人av视频| 激情 狠狠 欧美| 中文字幕av在线有码专区| 国产精品三级大全| 亚洲欧美日韩高清在线视频| 国产精品免费一区二区三区在线| 蜜桃久久精品国产亚洲av| 黄片wwwwww| 嘟嘟电影网在线观看| 欧美日本亚洲视频在线播放| 99热全是精品| 国产精品久久久久久亚洲av鲁大| 午夜福利视频1000在线观看| 九九在线视频观看精品| 一级毛片我不卡| av福利片在线观看| 精品久久久久久久末码| 人人妻人人澡人人爽人人夜夜 | 欧美最新免费一区二区三区| 精品人妻熟女av久视频| 26uuu在线亚洲综合色| 国产在线精品亚洲第一网站| 欧美一区二区亚洲| 小说图片视频综合网站| 欧美色视频一区免费| 欧美日韩在线观看h| 亚洲美女搞黄在线观看| 国产一区二区激情短视频| 亚洲国产精品合色在线| 99在线视频只有这里精品首页| 神马国产精品三级电影在线观看| 99久久人妻综合| 亚洲熟妇中文字幕五十中出| 国产精品野战在线观看| 一个人观看的视频www高清免费观看| 在线天堂最新版资源| 亚洲国产精品成人久久小说 | 国产乱人偷精品视频| 久久久精品94久久精品| 伦理电影大哥的女人| 国产精品久久久久久精品电影小说 | 成人三级黄色视频| 天堂网av新在线| 日本av手机在线免费观看| 日本黄色视频三级网站网址| 国产69精品久久久久777片| 我的女老师完整版在线观看| 久久久久免费精品人妻一区二区| 久久久久九九精品影院| 人人妻人人澡欧美一区二区| 久久人妻av系列| 乱人视频在线观看| 亚洲成人av在线免费| 久久久久性生活片| 丰满人妻一区二区三区视频av| 尾随美女入室| 国产精华一区二区三区| 美女大奶头视频| 精品久久久久久久久久免费视频| 天美传媒精品一区二区| 色哟哟哟哟哟哟| 国产精品蜜桃在线观看 | 国产久久久一区二区三区| 久久久精品94久久精品| 国产高清不卡午夜福利| 久久99热6这里只有精品| 免费电影在线观看免费观看| 91狼人影院| 午夜福利视频1000在线观看| 精品人妻视频免费看| 中文字幕久久专区| 色播亚洲综合网| 男的添女的下面高潮视频| 中文字幕制服av| 亚洲内射少妇av| 色综合色国产| 国产麻豆成人av免费视频| 男女下面进入的视频免费午夜| 国产激情偷乱视频一区二区| 精品久久久久久久末码| 最好的美女福利视频网| 国产成人精品一,二区 | 夜夜看夜夜爽夜夜摸| 国产单亲对白刺激| 长腿黑丝高跟| 网址你懂的国产日韩在线| 99热6这里只有精品| www.色视频.com| 国产午夜精品论理片| 乱码一卡2卡4卡精品| 国产黄a三级三级三级人| 国产高清三级在线| 精品久久久久久久人妻蜜臀av| 婷婷色综合大香蕉| 18禁在线无遮挡免费观看视频| 日本黄大片高清| 亚洲成a人片在线一区二区| 久久99精品国语久久久| 女人被狂操c到高潮| 国产精品av视频在线免费观看| 午夜精品一区二区三区免费看| 美女xxoo啪啪120秒动态图| 一级黄色大片毛片| 精品人妻熟女av久视频| 国产精品无大码| 99久久人妻综合| 少妇的逼好多水| 两性午夜刺激爽爽歪歪视频在线观看| 69av精品久久久久久| 色5月婷婷丁香| 99久久无色码亚洲精品果冻| 高清午夜精品一区二区三区 | 欧美高清成人免费视频www| 国产黄片视频在线免费观看| 亚洲经典国产精华液单| 在现免费观看毛片| 国产av在哪里看| 亚洲国产精品久久男人天堂| 爱豆传媒免费全集在线观看| 少妇的逼水好多| 狠狠狠狠99中文字幕| 欧美日韩一区二区视频在线观看视频在线 | 久久精品国产清高在天天线| a级一级毛片免费在线观看| 精品免费久久久久久久清纯|