• <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高清观看| 国产精品综合久久久久久久免费| 欧美最黄视频在线播放免费| 1000部很黄的大片| 国产精品免费一区二区三区在线| 一本久久中文字幕| 国产精品久久电影中文字幕| 波多野结衣巨乳人妻| 欧美区成人在线视频| 亚洲国产精品成人综合色| 精品无人区乱码1区二区| 国产成人a区在线观看| 99久久99久久久精品蜜桃| 欧美另类亚洲清纯唯美| 特级一级黄色大片| 麻豆久久精品国产亚洲av| 亚洲 国产 在线| 亚洲精品日韩av片在线观看| 亚洲av美国av| 99热这里只有是精品50| 亚洲中文日韩欧美视频| 国产午夜精品论理片| 欧美日韩瑟瑟在线播放| ponron亚洲| 日韩欧美国产一区二区入口| 久久婷婷人人爽人人干人人爱| 一二三四社区在线视频社区8| 亚洲狠狠婷婷综合久久图片| 亚洲五月婷婷丁香| xxxwww97欧美| 真实男女啪啪啪动态图| 赤兔流量卡办理| 两性午夜刺激爽爽歪歪视频在线观看| 成人毛片a级毛片在线播放| 此物有八面人人有两片| 一个人免费在线观看电影| 国产国拍精品亚洲av在线观看| 少妇的逼好多水| 日本 av在线| 国产精品,欧美在线| 日韩欧美国产一区二区入口| 人妻丰满熟妇av一区二区三区| 午夜免费激情av| 一个人免费在线观看电影| eeuss影院久久| 大型黄色视频在线免费观看| 国产视频一区二区在线看| 国产一区二区三区视频了| 校园春色视频在线观看| 午夜日韩欧美国产| 亚洲人成电影免费在线| 国产国拍精品亚洲av在线观看| 久久精品夜夜夜夜夜久久蜜豆| 97超视频在线观看视频| 亚洲自拍偷在线| 在线a可以看的网站| 成人国产一区最新在线观看| 一进一出好大好爽视频| 亚洲av成人av| h日本视频在线播放| 久久精品国产亚洲av涩爱 | 国内毛片毛片毛片毛片毛片| 成年人黄色毛片网站| 国产精品爽爽va在线观看网站| 99久久无色码亚洲精品果冻| 国产欧美日韩一区二区精品| 性欧美人与动物交配| 成年女人毛片免费观看观看9| 成年免费大片在线观看| 久久久国产成人免费| 欧美黄色淫秽网站| aaaaa片日本免费| av女优亚洲男人天堂| 国产亚洲av嫩草精品影院| 51午夜福利影视在线观看| 我要看日韩黄色一级片| 欧美成人一区二区免费高清观看| 老鸭窝网址在线观看| 成人无遮挡网站| 久久天躁狠狠躁夜夜2o2o| 日韩高清综合在线| 亚洲黑人精品在线| 变态另类丝袜制服| 岛国在线免费视频观看| 男女视频在线观看网站免费| 成人午夜高清在线视频| 观看美女的网站| 人人妻人人澡欧美一区二区| 国产色爽女视频免费观看| 日韩 亚洲 欧美在线| 精品一区二区三区视频在线| 无人区码免费观看不卡| 久9热在线精品视频| 最近中文字幕高清免费大全6 | 久久久久久九九精品二区国产| 久久久成人免费电影| 色播亚洲综合网| 99热这里只有是精品在线观看 | 91av网一区二区| 精品乱码久久久久久99久播| 热99re8久久精品国产| 在线国产一区二区在线| 久久久久久九九精品二区国产| a级毛片免费高清观看在线播放| 欧美色欧美亚洲另类二区| 人人妻人人澡欧美一区二区| 欧美区成人在线视频| 国产激情偷乱视频一区二区| 中出人妻视频一区二区| 中文字幕精品亚洲无线码一区| 国产高清视频在线观看网站| 日韩成人在线观看一区二区三区| 久久国产乱子伦精品免费另类| 91久久精品电影网| 男人舔女人下体高潮全视频| 亚洲欧美日韩高清专用| 亚洲人成网站高清观看| 精品一区二区三区视频在线观看免费| 我要看日韩黄色一级片| 51国产日韩欧美| 亚洲天堂国产精品一区在线| 全区人妻精品视频| 亚洲av免费在线观看| 亚洲精品粉嫩美女一区| 久久久久性生活片| 日韩欧美在线乱码| 搡女人真爽免费视频火全软件 | 国产又黄又爽又无遮挡在线| 成人国产综合亚洲| 国产熟女xx| av在线蜜桃| 亚洲真实伦在线观看| 国产精品久久久久久久电影| 国产视频一区二区在线看| 激情在线观看视频在线高清| 亚洲午夜理论影院| 最近最新中文字幕大全电影3| 中文字幕av在线有码专区| 免费av不卡在线播放| 免费看美女性在线毛片视频| 久久久久久国产a免费观看| 久久久久久久午夜电影| 99国产综合亚洲精品| 国产亚洲av嫩草精品影院| 最新在线观看一区二区三区| 人妻丰满熟妇av一区二区三区| 欧美最新免费一区二区三区 | 国产欧美日韩精品亚洲av| 国产中年淑女户外野战色| 日本熟妇午夜| 美女高潮喷水抽搐中文字幕| 久久久久性生活片| 久久精品国产亚洲av涩爱 | 脱女人内裤的视频| 日韩中字成人| 禁无遮挡网站| 国产午夜精品论理片| 亚洲片人在线观看| 午夜精品久久久久久毛片777| 观看免费一级毛片| 两性午夜刺激爽爽歪歪视频在线观看| 成人高潮视频无遮挡免费网站| 在线a可以看的网站| 日本五十路高清| 国产野战对白在线观看| 亚洲自偷自拍三级| 99久久精品热视频| 最近中文字幕高清免费大全6 | 国产精品女同一区二区软件 | 熟女人妻精品中文字幕| 99久久精品一区二区三区| 三级国产精品欧美在线观看| 亚洲自拍偷在线| 日韩欧美精品免费久久 | 亚洲激情在线av| 亚洲国产精品成人综合色| 三级毛片av免费| 床上黄色一级片| 国产精品永久免费网站| 亚洲人成网站在线播| 乱码一卡2卡4卡精品| 特大巨黑吊av在线直播| 青草久久国产| 噜噜噜噜噜久久久久久91| 亚洲国产欧洲综合997久久,| 热99在线观看视频| 一级毛片久久久久久久久女| 精品午夜福利视频在线观看一区| 中亚洲国语对白在线视频| 国产单亲对白刺激| 欧美激情久久久久久爽电影| www.色视频.com| 最近中文字幕高清免费大全6 | 亚洲av成人av| 搡老岳熟女国产| 国产白丝娇喘喷水9色精品| 亚洲av五月六月丁香网| 在线a可以看的网站| 又爽又黄无遮挡网站| 午夜免费男女啪啪视频观看 | 国语自产精品视频在线第100页| 亚洲成人久久爱视频| 深夜精品福利| 男女做爰动态图高潮gif福利片| 国产一区二区激情短视频| 亚洲av免费在线观看| 国产亚洲精品久久久久久毛片| 国产精品乱码一区二三区的特点| 搡老熟女国产l中国老女人| 精品熟女少妇八av免费久了| 国模一区二区三区四区视频| 观看美女的网站| 久久性视频一级片| 欧美日韩福利视频一区二区| 日本a在线网址| 他把我摸到了高潮在线观看| 亚洲av美国av| 国产一区二区激情短视频| 脱女人内裤的视频| 久久久久久久久久成人| 欧美不卡视频在线免费观看| 精品福利观看| 99热精品在线国产| 91在线观看av| 白带黄色成豆腐渣| 一进一出抽搐gif免费好疼| 日本a在线网址| 色综合亚洲欧美另类图片| 精品一区二区三区人妻视频| 精品一区二区三区人妻视频| 日韩人妻高清精品专区| 日韩人妻高清精品专区| 国产三级在线视频| 亚洲精品亚洲一区二区| 亚洲国产精品999在线| 全区人妻精品视频| 免费av不卡在线播放| 亚洲精品一卡2卡三卡4卡5卡| 久久国产乱子伦精品免费另类| 国产免费一级a男人的天堂| ponron亚洲| 色在线成人网| 村上凉子中文字幕在线| 神马国产精品三级电影在线观看| 国产亚洲精品综合一区在线观看| 国产av麻豆久久久久久久| 在线观看66精品国产| av在线观看视频网站免费| 亚洲精品日韩av片在线观看| 国产高清激情床上av| www.999成人在线观看| 久久久久久大精品| 久久久久久久午夜电影| 日本免费a在线| 中亚洲国语对白在线视频| 97超级碰碰碰精品色视频在线观看| 亚洲在线观看片| 日本熟妇午夜| 在线观看66精品国产| 91久久精品电影网| 露出奶头的视频| 成人特级av手机在线观看| a级毛片免费高清观看在线播放| 国产欧美日韩一区二区三| 色吧在线观看| 青草久久国产| 看十八女毛片水多多多| 九色成人免费人妻av| 欧美极品一区二区三区四区| 欧美+亚洲+日韩+国产| 草草在线视频免费看| 最近在线观看免费完整版| 看片在线看免费视频| 成人av一区二区三区在线看| 欧美乱妇无乱码| 99久久99久久久精品蜜桃| 一二三四社区在线视频社区8| 午夜激情福利司机影院| 国产高清视频在线播放一区| 欧美性感艳星| 精品99又大又爽又粗少妇毛片 | 变态另类丝袜制服| 日韩欧美国产在线观看| 99久久九九国产精品国产免费| 日韩中文字幕欧美一区二区| 中出人妻视频一区二区| 简卡轻食公司| 一进一出好大好爽视频| 国产一区二区三区在线臀色熟女| 两人在一起打扑克的视频| 国产精品99久久久久久久久| 国产伦精品一区二区三区视频9| 国产精品亚洲av一区麻豆| 国产精品亚洲美女久久久| 国产主播在线观看一区二区| 成人av在线播放网站| 天堂影院成人在线观看| 熟女人妻精品中文字幕| 欧美性感艳星| 在线观看午夜福利视频| 啪啪无遮挡十八禁网站| 亚洲人成电影免费在线| 成人精品一区二区免费| 国产精品av视频在线免费观看| 久久精品夜夜夜夜夜久久蜜豆| 日本一二三区视频观看| 成年人黄色毛片网站| 波多野结衣巨乳人妻| 亚洲成人精品中文字幕电影| 麻豆av噜噜一区二区三区| 国产一区二区三区视频了| 久久亚洲真实| 一个人看的www免费观看视频| 五月伊人婷婷丁香| 老司机福利观看| 午夜免费成人在线视频| bbb黄色大片| 国产午夜福利久久久久久| 国模一区二区三区四区视频| 久久国产乱子伦精品免费另类| 18禁裸乳无遮挡免费网站照片| 成年版毛片免费区| 日日干狠狠操夜夜爽| 亚洲午夜理论影院| 综合色av麻豆| 波多野结衣巨乳人妻| 亚洲片人在线观看| 日本黄色片子视频| 亚洲无线观看免费| 老熟妇乱子伦视频在线观看| 亚洲电影在线观看av| 五月伊人婷婷丁香| 他把我摸到了高潮在线观看| 丰满乱子伦码专区| 日韩av在线大香蕉| 国产麻豆成人av免费视频| 亚洲 欧美 日韩 在线 免费| 丝袜美腿在线中文| 免费看a级黄色片| av国产免费在线观看| 非洲黑人性xxxx精品又粗又长| 色哟哟哟哟哟哟| a在线观看视频网站| 亚洲色图av天堂| 亚洲七黄色美女视频| 亚洲成a人片在线一区二区| 每晚都被弄得嗷嗷叫到高潮| 久久亚洲真实| 国内精品一区二区在线观看| 国产一区二区激情短视频| 国产精品亚洲av一区麻豆| 人人妻,人人澡人人爽秒播| 成人高潮视频无遮挡免费网站| 男人舔女人下体高潮全视频| 欧美性感艳星| 中文在线观看免费www的网站| 国产一级毛片七仙女欲春2| 国产黄a三级三级三级人| 色视频www国产| 成人国产一区最新在线观看| 99久久精品国产亚洲精品| 成人av一区二区三区在线看| 真人一进一出gif抽搐免费| 成人国产综合亚洲| 亚洲最大成人av| 特大巨黑吊av在线直播| 五月伊人婷婷丁香| 久久久国产成人免费| 欧美绝顶高潮抽搐喷水| 床上黄色一级片| 国产成人a区在线观看| 午夜福利在线在线| 人人妻人人看人人澡| 亚洲av免费在线观看| 色视频www国产| 亚洲真实伦在线观看| 熟女人妻精品中文字幕| 精品国产亚洲在线| 亚洲欧美日韩高清专用| 国产老妇女一区| 有码 亚洲区| 怎么达到女性高潮| 精品一区二区免费观看| 国产av一区在线观看免费| 大型黄色视频在线免费观看| 亚洲欧美清纯卡通| 国产精品爽爽va在线观看网站| 又黄又爽又刺激的免费视频.| 久久国产乱子免费精品| 亚洲成人久久爱视频| 日韩免费av在线播放| 亚洲黑人精品在线| 国产在线男女| 偷拍熟女少妇极品色| 亚洲av成人精品一区久久| 2021天堂中文幕一二区在线观| 久久性视频一级片| 人妻制服诱惑在线中文字幕| 精品免费久久久久久久清纯| 校园春色视频在线观看| 又紧又爽又黄一区二区| 听说在线观看完整版免费高清| 九色国产91popny在线| 亚洲真实伦在线观看| 午夜福利免费观看在线| 国产v大片淫在线免费观看| av在线蜜桃| 国产单亲对白刺激| 麻豆成人av在线观看| 国产精品自产拍在线观看55亚洲| av黄色大香蕉| 亚洲成人中文字幕在线播放| x7x7x7水蜜桃| 黄色日韩在线| 香蕉av资源在线| 高清在线国产一区| 18禁黄网站禁片免费观看直播| 在线观看美女被高潮喷水网站 | 国产午夜精品论理片| or卡值多少钱| 国产亚洲精品综合一区在线观看| 热99在线观看视频| 午夜久久久久精精品| 国产精品影院久久| 亚洲狠狠婷婷综合久久图片| 久久精品国产亚洲av香蕉五月| 精品国产亚洲在线| 精品久久久久久,| 午夜亚洲福利在线播放| 一边摸一边抽搐一进一小说| 日本黄色片子视频| 成人一区二区视频在线观看| 欧美三级亚洲精品| 美女高潮喷水抽搐中文字幕| 日本在线视频免费播放| 欧美黄色片欧美黄色片| 久久精品国产99精品国产亚洲性色| 欧美最黄视频在线播放免费| 久久久久久久午夜电影| 国产色婷婷99| 免费黄网站久久成人精品 | 一进一出抽搐gif免费好疼| 成人性生交大片免费视频hd| 午夜影院日韩av| 欧美日韩亚洲国产一区二区在线观看| 91麻豆精品激情在线观看国产| 国产aⅴ精品一区二区三区波| 99在线人妻在线中文字幕| 日韩欧美在线乱码| av视频在线观看入口| 亚洲经典国产精华液单 | 久久九九热精品免费| 亚洲精品成人久久久久久| 少妇高潮的动态图| 国产伦精品一区二区三区视频9| 黄色视频,在线免费观看| 最后的刺客免费高清国语| 亚洲最大成人中文| 亚洲最大成人av| 日韩大尺度精品在线看网址| 午夜a级毛片| 久久99热6这里只有精品| 99久久久亚洲精品蜜臀av| 久久精品综合一区二区三区| 99国产精品一区二区蜜桃av| 亚洲第一欧美日韩一区二区三区| 欧美一级a爱片免费观看看| 久久精品人妻少妇| 俺也久久电影网| 日本撒尿小便嘘嘘汇集6| 亚洲成人免费电影在线观看| 欧美日韩中文字幕国产精品一区二区三区| 别揉我奶头 嗯啊视频| 亚洲自偷自拍三级| 久久久久精品国产欧美久久久| 久久人人爽人人爽人人片va | 禁无遮挡网站| 欧美绝顶高潮抽搐喷水| 国内毛片毛片毛片毛片毛片| 在线国产一区二区在线| 亚洲aⅴ乱码一区二区在线播放| 在线观看舔阴道视频| 亚洲av成人av| 日韩欧美在线乱码| 床上黄色一级片| 久久久久久久久久成人| 国产精品爽爽va在线观看网站| 欧美一级a爱片免费观看看| a级毛片免费高清观看在线播放| 亚洲avbb在线观看| 亚洲一区二区三区色噜噜| 黄色视频,在线免费观看| 亚洲内射少妇av| 久久伊人香网站| 午夜福利高清视频| 国产三级黄色录像| 最新在线观看一区二区三区| 18禁黄网站禁片免费观看直播| 深爱激情五月婷婷| 中文在线观看免费www的网站| 国产精品av视频在线免费观看| 特大巨黑吊av在线直播| 久久精品久久久久久噜噜老黄 | 亚洲第一欧美日韩一区二区三区| 日韩成人在线观看一区二区三区| 1024手机看黄色片| 国产白丝娇喘喷水9色精品| 午夜福利免费观看在线| 久久人妻av系列| 国产免费av片在线观看野外av| 中文字幕人成人乱码亚洲影| 久久久久久久精品吃奶| 亚洲电影在线观看av| 午夜免费激情av| 婷婷精品国产亚洲av在线| 国产精品久久久久久久电影| 亚洲无线在线观看| 欧美激情久久久久久爽电影| 久久久色成人| 悠悠久久av| 精品午夜福利视频在线观看一区| 18禁黄网站禁片免费观看直播| 国产av一区在线观看免费| 久久久成人免费电影| 国产在线男女| 午夜福利免费观看在线| 国语自产精品视频在线第100页| 两个人视频免费观看高清| 精品一区二区三区视频在线观看免费| 成年人黄色毛片网站| 神马国产精品三级电影在线观看| 麻豆av噜噜一区二区三区| 一进一出抽搐gif免费好疼| 亚洲av电影不卡..在线观看| 国内揄拍国产精品人妻在线| 真人做人爱边吃奶动态| or卡值多少钱| 日日摸夜夜添夜夜添小说| 亚洲av熟女| 亚洲va日本ⅴa欧美va伊人久久| 男女之事视频高清在线观看| 成人欧美大片| 欧美黑人巨大hd| 成人国产综合亚洲| 99久久无色码亚洲精品果冻| 色噜噜av男人的天堂激情| 一进一出抽搐动态| 午夜福利在线在线| 亚洲激情在线av| 国产精品不卡视频一区二区 | 久久天躁狠狠躁夜夜2o2o| 在线观看av片永久免费下载| 欧美成人性av电影在线观看| 亚洲不卡免费看| 亚洲第一电影网av| 天天一区二区日本电影三级| 搞女人的毛片| 日本精品一区二区三区蜜桃| 在线观看午夜福利视频| 又黄又爽又刺激的免费视频.| 国产三级在线视频| 在线看三级毛片| 91麻豆精品激情在线观看国产| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 日韩欧美 国产精品| 简卡轻食公司| 好男人电影高清在线观看| 在线看三级毛片| 日韩国内少妇激情av| 男插女下体视频免费在线播放| 国产午夜福利久久久久久| 在现免费观看毛片| 全区人妻精品视频| 麻豆国产av国片精品| 在线观看免费视频日本深夜| 午夜福利在线观看吧| av天堂在线播放| 久久精品国产99精品国产亚洲性色| 国语自产精品视频在线第100页| 性插视频无遮挡在线免费观看| 午夜视频国产福利| 国产成人影院久久av| 亚洲成人精品中文字幕电影| 啦啦啦韩国在线观看视频| 国产精品一及| 日韩av在线大香蕉| 91在线观看av| 国产亚洲精品久久久com| 搡老岳熟女国产| 国产久久久一区二区三区| 俄罗斯特黄特色一大片| 欧美精品啪啪一区二区三区| 国产人妻一区二区三区在| 午夜福利欧美成人| 欧洲精品卡2卡3卡4卡5卡区| 90打野战视频偷拍视频| 内射极品少妇av片p| 波多野结衣高清作品| 婷婷精品国产亚洲av| 国产精品精品国产色婷婷| 亚洲欧美日韩东京热| 欧美成人a在线观看| 最近最新中文字幕大全电影3| 深夜精品福利| 看片在线看免费视频| 午夜老司机福利剧场| 91在线观看av| 国产精品嫩草影院av在线观看 | 亚洲第一区二区三区不卡| 久久九九热精品免费| 热99在线观看视频| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 成年人黄色毛片网站| 在线观看午夜福利视频| av天堂中文字幕网|