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

    Impacts of climate change on water quantity,water salinity,food security,and socioeconomy in Egypt

    2021-05-02 13:45:12MohieElDinMohameOmarAhmeMoustafaAhmeMoussaReinharHinkelmann
    Water Science and Engineering 2021年1期

    Mohie El Din Mohame Omar ,Ahme Moustafa Ahme Moussa ,Reinhar Hinkelmann

    a National Water Research Center,Fum Ismailiya Canal,Shoubra El-Kheima 13411,Egypt

    b International Center for Agricultural Research in the Dry Areas,Cairo 2416,Egypt

    c Nile Research Institute(NRI),National Water Research Center,El-Qanater 13621,Egypt

    d Department of Water Engineering,Technische Universit¨at Berlin(TU Berlin)Campus El Gouna,El Gouna 84513,Egypt

    Received 16 April 2020;accepted 30 August 2020 Available online 23 March 2021

    Abstract Climate change might have direct impacts on water quantity in Egypt and lead to indirect effects on Mediterranean saltwater intrusion to groundwater,which exposes agriculture to vulnerability.This study investigated impacts of climate change on agriculture,with particular regard to food security and socioeconomy,and quantified the effectiveness of cropping pattern adaptation measures by integrating three mathematical models.The BlueM model was used for hydrological simulations of Nasser Lake under flooding scenarios to predict the water supply from the High Aswan Dam.The water and salinity balance(WB-SAL)model was adopted to estimate the water salinity in the Nile Delta.The simulated results from the BlueM and WB-SAL models were integrated with the agricultural simulation model for Egypt(ASME)to project cropping patterns,food security,and socioeconomy throughout the country.The results showed that future climate change will directly affect the total crop area;crop areas for 13 crop types;the self-sufficiency of wheat,rice,cereal,and maize supplies;and socioeconomic indicators.The proposed cropping pattern adaptation measures focus on fixing the crop areas of rice and orchards and providing half of the population with lentils,maize,onion,vegetables,milk,and meat.The adaptation measures have the potential to promote food security without causing deterioration of the socioeconomic situation.However,water availability has much more significant effects on food security and socioeconomy than cropping pattern adaptation measures do.Accordingly,the country should rationalize water use efficiency and increase water supply.? 2021 Hohai University.Production and hosting by Elsevier B.V.This is an open access article under the CC BY-NC-ND license(http://creativecommons.org/licenses/by-nc-nd/4.0/).

    Keywords:Climate change;Water quantity;Salinity;Cropping pattern adaptation measures;Food security;Socioeconomy

    1.Introduction

    Egypt is almost entirely dependent on the water of the Nile River,which supplies an annual water volume of approximately 55.5 km3,accounting for 93% of the country"s conventional water resources.The annual total water demand is 81.3 km3,with approximately 86% for agriculture,2.5% for industry,and 11.5% for the domestic sector(Omar and Moussa,2016).There is a gap between water demand and water availability,which is compensated for by the reuse of drainage water,wastewater,and shallow groundwater.The uncertainty of climate change impacts is another challenge to the water resources system in Egypt.

    Climate change affects air temperature and precipitation around the world.It tends to result in sea-level rise,thereby influencing groundwater hydraulics and causing greater seawater intrusion in many coastal aquifers(Ketabchi et al.,2016),including the coastal aquifers in Egypt.The global climate has become 0.5°C warmer over the last 100 years due to greenhouse gas emissions caused by human activities(WMO,2017).Air temperature rise and precipitation variation have the potential to change water flows and might lead to intensified extreme hydrological events.However,previous studies focusing on the uncertainty of water flow variation have produced contrasting results.Elshamy and Wheater(2009)used the bias-corrected statistical downscaling approach to downscale 17 general circulation models(GCMs)and projected that by the end of this century the Blue Nile flow would have a change rate ranging from-60% to 45%.Strzepek and McCluskey(2007)used five GCMs based on two emission scenarios and estimated the variations of the Nile flows entering Nasser Lake in 2050 and 2100 under 20 scenarios.Twelve reduced flow scenarios and eight increased flow scenarios were projected.Kotb(2015)presented six scenarios using the Quantifying Uncertainty in Model Projections(QUMP)model and the regional climate model(RCM),which showed a precipitation increase of 5%-11% and a flow increase of 7.15%-23.21% at the Dongola Station at the Nasser Lake entrance.Regarding the sea-level rise,Eissa et al.(2017)used an analytical model combining hydrogeological and geochemical characterization of aquifer geometry to provide saltwater and freshwater balances at the Ras El Hekma site on the Egyptian Mediterranean Sea coast.Sefelnasr and Sherif(2014)used the finite element subsurface flow system model and found that large areas in the coastal zone of the Nile Delta would be submerged in seawater.UNESCO(2013)quantified the regions vulnerable to sea-level rise in Egypt.The increase in shallow groundwater salinity and drainage water reuse would lead to a gradual increase in the salt content of irrigation water in the Nile Delta.In Egypt,possible risks induced by climate change include the following:(1)water scarcity due to the reduction of water supply from the Nile River and high water requirement for crops,and(2)salinity increase in irrigation water in the Nile Delta resulting from saltwater intrusion to shallow groundwater and water reuse.

    Agriculture consumes approximately 86% of total water resources in Egypt.Irrigated crops contribute mostly to local food production and exports and therefore to the prevention of saltwater intrusion in the Northern Nile Delta by rice cultivation.Climate change in Egypt would decrease crop yields for most crops,with wheat yields expected to be reduced by up to 9% in 2030 and by close to 20% in 2060(Smith et al.,2013).Hence,it is necessary to comprehensively investigate future climate change impacts on the Nile River flow,evapotranspiration,and shallow groundwater salinity,and how those impacts will affect the agricultural sector,with regard to the vulnerability of food security and the socioeconomy.It is also important to quantify the effectiveness of adaptation measures.

    The most direct and simplest adaptation measure is to encourage farmers to change cropping patterns.However,the consequences of changing cropping patterns might have adverse impacts on food security and the socioeconomy.Adhikari(2018)reported that farmers in Nepal abandoned rice,cereal crops,tubers,and sweet potatoes,and even abandoned agriculture on their lands throughout 15 drought events owing to climatic conditions in the period from 1972 to 2015.This negatively affected agricultural production and the socioeconomic situation throughout the country.Alabdulkader et al.(2016)reported adverse effects of climate change on the agricultural sector in Saudi Arabia,and presented changes in the date palm cropping pattern as an adaptation measure,which could increase net annual return and decrease water demand.Chebil et al.(2019)investigated the impacts of climate change on 21 strategic crops in Tunisia and found that the total cropped area and agricultural employment in Tunisia would be negatively affected.However,the study did not present the effectiveness of changing cropping patterns as an adaptation measure.

    In Egypt,few previous studies have quantified the effectiveness of changing cropping patterns as climatic adaptation measures.To adapt to the possible impacts of climate change on agriculture in Egypt,including reduction in Nile River flows,air temperature rise,Mediterranean sealevel rise,and salinity increase,changing the cropping structure is required.Therefore,the main objective of this study was to investigate impacts of climate change on water quantity and salinity in Egypt and their consequent impacts on food security and the socioeconomy.In addition,the effectiveness of cropping pattern adaptation measures was quantified.To achieve these goals,climate change impact projections were conducted,with regard to the Nile River flow,air temperature,open water evaporation,evapotranspiration,precipitation,and aquifer salinity in the Nile Delta.Afterward,the impacts of climate projections and adaptation measures on food security indicators,self-sufficiency,and cropping intensity were investigated.Finally,the impacts of climate projections and adaptation measures on socioeconomic indicators,agricultural productivity,consumerproducer surplus,number of jobs,and water marginal value were quantified.

    2.Materials and methods

    In this study,data of climate change projections and demographic changes were collected.Afterward,the BlueM river basin management model was used to predict the water supply from the High Aswan Dam(HAD)based on future flooding scenarios and constraints of reservoir storage volume.Meanwhile,the water and salinity balance(WB-SAL)model was adopted to predict salinity changes in the Nile Delta region according to future HAD release from the BlueM model and the expected demographic and climatic conditions.Finally,the agricultural simulation model for Egypt(ASME)was used to predict the optimal cropping pattern,food security,and socioeconomic indicators in the future based on the results of the BlueM and WB-SAL models.Fig.1 shows the study area.

    2.1.Data acquisition

    Fig.1.Study area.

    Three data sets were obtained.The first was the input data for the BlueM model,which included climate change projections of the Nile River flow at the entrance of Nasser Lake,meteorological conditions in the lake region,seepage losses,and constraints controlling the lake operation processes.This BlueM-based simulation estimated the water release from the HAD.The output data in addition to the second data set were used as input data for the WB-SAL model.The second data set mainly focused on the Nile Delta region and consisted of meteorological conditions,evapotranspiration,the Mediterranean seawater intrusion rate to shallow groundwater,volumes and salinities of water supplies,and salt emissions from agricultural crops and industrial sectors.The outputs of the WB-SAL model were the salinities of irrigation water and drainage water.The third data set included the data affecting agriculture throughout the country,including water resources,demographic changes,agricultural areas,and socioeconomic conditions.These data,in addition to the results of the BlueM and WB-SAL models,were used as the input data for the ASME model.

    2.1.1.Climate change projections

    Many studies have combined GCMs and RCMs with hydrological models to investigate the link between climate and water resources.GCMs estimate the effects of emissions on global climate and describe the physical processes in the atmosphere,in the oceans,and on the land surface.By contrast,RCMs add detailed information about future climate change on fine scales to the GCMs’large-scale projections.At the Dongola Station at the Nasser Lake entrance,Strzepek and McCluskey(2007)used GCMs to predict water flow in 2050.In terms of RCMs,Kotb(2015)used six perturbed physical ensembles developed by the QUMP model to downscale climate projections over the Nile River Basin.Because RCMs account for the regions with mountains,coastlines,and islands on the scale of 100 km or less,and RCMs are more realistic than GCMs,this study utilized the results of Kotb(2015)that all future flows would increase.This study also adopted the findings of Strzepek and McCluskey(2007)that in six out of ten cases,river flow would decrease.Additionally,flow reductions projected from GCMs were considered,because the risks of flow reduction in arid countries such as Egypt cannot be ignored.Therefore,this study constructed one scenario that represented an average river flow increase of 14.3% according to RCMs and another scenario with an average flow decrease of 11.8% according to GCMs.Table 1 presents various climate change projections for the Nile River flow,air temperature,precipitation,open water evaporation,evapotranspiration,and Mediterranean saltwater intrusion rate.

    2.1.2.Future planning alternatives

    The present study collected other data affecting water demand.The population increase,expansion of agricultural areas,and establishment of new communities were combined with different climate scenarios.In 2050,population would be expected to reach 159 956 808 in Egypt(UN,2019).Due to the agriculture expansion plan,Egypt is expected to increase agricultural area by 315 000 hm2(Omar and Moussa,2016).

    2.2.Mathematical modeling

    2.2.1.BlueM model

    The BlueM model was used in this study to simulate the operational process of Nasser Lake considering both the decrease and increase of river flow entering the lake.The BlueM model was developed by the Darmstadt University of Technology,in Germany,for river basin management(Bach et al.,2009).The operation of Nasser Lake was described by the water balance equation under various constraints concerning storage volume,outflow from the lake,and water losses,and water balance calculation was performed on a monthly basis,as follows:

    Table 1Changes in different parameters in Egypt in 2050 due to climate change.

    whereItis the mean inflow in montht(m3);Qtis the outflow from the dam(m3);Mtis the flow released from the emergency spillway(m3);Dtis the water demand of the Toshka Project for reclamation of 226 800 hm2from Nasser Lake(m3);Ttis the volume of water released from the Toshka Spillway to empty the reservoir down to the level of 175 m before floods(m3);Stis the seepage losses from the lake(m3),which was computed using the approach of Ebaid and Ismail(2010);andEtis the mean open water evaporation from the lake(m3),which is expressed as follows:

    whereAtandAt+1are the lake areas at the beginning and end of montht,respectively;andCtis the evaporation coefficient in montht,which was calculated using the method of Ebaid and Ismail(2010).Model calibration was conducted by comparing the simulated daily water level upstream of the HAD with the observations,with the absolute error as the error index.

    2.2.2.WB-SAL model

    Owing to the intensive drainage water reuse,salt fluxes from the Nile Delta to the Mediterranean Sea and coastal lakes are higher than the fluxes entering the delta.As climate change might result in severe water shortage and saltwater intrusion,the salinity of drainage,surface irrigation,and shallow groundwater in the delta region may increase in the future.This would threaten agricultural production and its related socioeconomic aspects in the delta region.Therefore,it is necessary to project the future salinity of water bodies in the delta region.In general,traditional water quality models simulate the concentration of materials in any water system as a function of space and time based on different water volumes and loads entering and exiting the system.This means that the salt loads or concentrations of agricultural,municipal,and industrial disposals should be collected,measured,or assumed as the model input.However,it is a difficult,time-consuming,and high-cost process to collect these data in large basins such as the delta region.Hence,the WB-SAL model was used to project future salt emissions of agricultural,municipal,and industrial disposals to water bodies.The loads of agricultural drainage were estimated according to cropping patterns.The loads of industrial wastewater were based on industry categories and their share percentages.Afterward,the model was used to estimate the salinity of mixed drainage and irrigation water,with the total dissolved solids(TDS)as the salinity indicator.In the WB-SAL model,the salt balance is calculated based on the fact that dissolved salts are conservative and not subject to transformation,sedimentation,or biomass removal.When water is evaporated,the dissolved salts are left.The salt balance equation in the WB-SAL model is expressed as follows:

    whereVoandCoare the volume and salinity of drainage water that discharges to the Mediterranean Sea,respectively;ViandCiare the volume and salinity of the Nile River flow entering the delta region,respectively;VpandCpare the volume and salinity of precipitation,respectively;VRandCRare the volume and salinity of reused drainage,respectively;andVDandCDare the volume and salinity of water diverted for different demands,respectively,and are expressed as follows:

    whereVAandCAare the volume and salinity of diverted water for agricultural demand,respectively;VMandCMare the volume and salinity of diverted water for municipal demand,respectively;VIandCIare the volume and salinity of diverted water for industrial demand,respectively;andVAqandCAqare the volume and salinity of diverted water for aquacultural demand,respectively.The salinity of irrigation water was estimated using the following equation:

    whereVSandCSare the volume and salinity of surface water,respectively;andVSGandCSGare the volume and salinity of shallow groundwater,respectively.

    In this study,drainage salinity was calculated as a mixture of agricultural drainage,industrial wastewater,treated wastewater,untreated wastewater,and aquaculture drainage.Some drainage water is reused,and the remaining drainage water is discharged to the Mediterranean Sea and northern lakes via drainage systems.Hence,the salinity of reused water was defined as drainage salinity.Salt loads from agricultural drainage and industrial wastewater were estimated from emission factors.For agricultural drainage,emission factors for different crops were adopted according to those from Deltares(2012).Based on current cropping patterns,salt loads were separately estimated using the WBSAL model.Similarly,loads from industrial wastewater were derived from the emission factors of various sectors and their percentage of the total industry.Industrial emission factors were adopted according to those developed by Wahab and Badawy(2004).The share percentages of various industrial sectors in each governorate were taken from Deltares(2012).The share percentages of various industrial sectors for the entire delta region were calculated.The emission loads were added to the original loads in the supplied water.The salinity from agriculture and industry was calculated by dividing the salt loads by the volumes of agricultural drainage and industrial wastewater,respectively.Salt emission loads were calculated separately and were linked to water and salt pathways as shown in Fig.2.

    The salinities of treated and untreated wastewater were estimated according to the findings of El Gammal(2011),who collected water quality data of raw wastewater at seven locations in Egypt.This study defined the salinity required by the WB-SAL model as the average salinity at the seven locations in El Gammal(2011).The salinity of aquacultural drainage water was adopted from the findings of Emara et al.(2016),who analyzed the effluent of an aquacultural farm discharging to Burullus Lake in the northern Nile Delta on a seasonal basis,and found an average TDS of 3 000 mg/L.Based on the predictions from the BlueM model,factors affecting Nile River flows were altered.Due to both cultivation expansion and urbanization,the agricultural area will change.Total irrigation demand was estimated from the agricultural area as well as from water requirements for crops,which will increase due to climate change.Municipal water demand was dependent on population and the water consumption rate per capita.Treated and untreated wastewater volumes were estimated from wastewater discharge per capita and the ratio of the volume of treated wastewater to that of total wastewater.The calibration process was performed by comparing the simulated salinity of drainage water to the coastal lakes and the Mediterranean Sea with the observations.

    Fig.2 shows the water and salt balances calculated by the WB-SAL model.The current annual conventional water resources in the delta region are 31.7 km3from the Nile River and 1.0 km3from effective rainfall.By contrast,the current annual water demand is 38.39 km3for agriculture,4.63 km3for municipal use,1.16 km3for industry,and 0.80 km3for aquaculture.The gap between water demand and supply is approximately 12.67 km3,which is compensated for by the recycling of drainage water and wastewater as well as utilization of shallow groundwater.

    Fig.2.Water and salt balances calculated by WB-SAL model.

    2.2.3.TDS measurements

    Water sampling was conducted for the Nile River flow entering the delta.The measured TDS was used as the input data for the WB-SAL model.Drainage water sampling was also conducted at the drainage sites of Bahr Hadus,El Serw,Bahr El Baqar,El-Omoum,Nasser,East Burullus,and West Burullus before drainage water flowed into the coastal lakes and the Mediterranean Sea.The measured TDS concentrations at the drainage sites were used for comparison with the simulated TDS from the WB-SAL model.Water sampling was performed in both summer and winter.

    2.2.4.ASME model

    ASME is an optimization model that determines the cropping patterns needed to produce the highest agriculture-related welfare with a certain availability of water.The optimization mode of ASME has no restrictions on crop area but merely has rotational restrictions to ensure feasible cropping patterns.In this study,several restrictions were set up,including fixed crop areas for specific crops,controlled import rates,and securing part of the population with specific crops.These restrictions enabled ASME to predict the optimal cropping pattern to achieve optimal food security and socioeconomy under different climatic scenarios.In addition,they enabled ASME to quantify effective adaptation measures.

    ASME is programmed with the General Algebraic Modeling System.The mathematical model was developed and applied at a governorate level using linear programming.The model is static when irrigation water is allocated for a year among different users.The objective function of ASME is to maximize the welfare of producers by aggregating gross margins from crop production with limited land and water resources.The objective function(Z)maximizes the total gross margins from all crops by selecting the optimal crop grouping subject to a set of constraints.In addition,the maximization of net return in the short run is equivalent to the maximum gross margin.The objective function at each level of mathematical analysis is as follows:

    whereZrepresents the value of the objective function(USD),Pijis the unit price of cropjin governorateipaid to the producers(USD/t),Yijis the yield per unit area for cropjin governoratei(t/hm2),Cijisthe investment costper unitareafor cropjin governoratei(USD/hm2),andXijis the crop area for cropjin governoratei(hm2).Theoptimalareaofeachcropdependsonthe total amount of water and crop water demand.The maximization of gross margin per unit area is equivalent to the maximum gross margin per unit water.This helps to estimate the total quantity of water that should be used for a given crop:

    whereWtis the total amount of water available for irrigation(m3),andWijis the crop water demand for cropjin governorateiin monthm(m3).

    The ASME model decomposes the commodities into directuse crop commodities,processing-input crop commodities,processing-output crop commodities,processing-by-product crop commodities,and livestock commodities.The water system in ASME consists of the Nile River water from the HAD and irrigation canals.Variation of future water supply from the HAD due to climate change estimated from the BlueM model was considered in future scenarios.Evaporation from open water and fallow land was considered.Agricultural water use refers to evapotranspiration,part of which is compensated for by effective rainfall.Net of rain and losses abstracted from the Nile system werecalculatedatagovernoratelevel.Evaporationwasestimated at the national level and was proportionally distributed over the study area.Changes in rainfall,evaporation,and evapotranspiration arising from climate change were obtained and considered in future scenarios.Future salinity alteration in the delta region calculated by the WB-SAL model was considered in the delta governorates.Municipal,industrial,and agricultural drainage water returns to the Nile River,canals,sea,or lakes via gravity or pumping stations.Additionally,ASME calculated the deep groundwater balance in New Valley,Matruh,Red Sea,and South Sinai governorates.To calibrate the ASME model,the estimated volume of reused drainage was compared with the observations.

    2.2.5.Measures in adaptation to climate change

    Mitigation measures for emission reduction are not enough,because certain greenhouse gases will still exist in the atmosphere.Adaptation is essential in reducing these impacts.We enacted a set of cropping pattern adaptation measures as follows:

    (1)Only 294 000 hm2of rice was cultivated in coastal governorates,such as Beheira,Kafr El Sheik,Dakhalia,Damietta,Sharkia,and Port Said.

    (2)Orchard area was fixed at the 2017 level based on a questionnaire in this study in the five governorates of Fayoum,Qena,Damietta,Behira,and Sharkia.It was found that 96%of the farmers prioritized orchards in the case of water shortage.

    (3)Half of the population demand for lentils,maize,onions,vegetables,meat,and milk was produced in each governorate to guarantee the food self-sufficiency of the rural population.

    (4)According to FAO(2011),Egypt lost 13%-15% of available cereal between harvesting and final consumption.Hence,a loss of 14% was considered in this study.

    2.2.6.Modeling scenarios

    The current scenario and three future scenarios were established(Table 2).Due to the uncertainty of climate change projections of Nile River flows and water salinity in the delta region,the three future scenarios were set to investigate the impacts of all possible projections and the suggested cropping pattern adaptation measures on cropping patterns,food security,and the socioeconomy in Egypt.

    Scenario 1 investigated the effects of the reduction in water supply from the HAD and the consequent increase of water salinity in the delta region.The reduced water supply from theHAD was estimated from the BlueM model based on the GCM-predicted average reduction percentage of the Nile River flow at the Nasser Lake entrance(11.8%)(Table 1).The change of water salinity in the delta region was estimated by the WB-SAL model,based on the average reduction rate of the Nile River flow at the Nasser Lake entrance and such parameters as air temperature,precipitation,open water evaporation,evapotranspiration,and shallow groundwater salinity in the delta region(Table 1).This scenario did not assess adaptation measures,but only assumed that farmers would adapt to the reduction of cultivated area and livestock as well as the same 30 crop types currently produced.A linear relationship was assumed between water availability and cultivation area.Additionally,this scenario assumed that water shortage would be proportionally distributed across governorates.The change rate of salinity in irrigation water was considered in coastal governorates such as Beheira,Kafr El Sheik,Dakhalia,Damietta,Sharkia,and Port Said.

    Table 2Basic information for current status and three future scenarios.

    Scenario 2 evaluated the impacts of water supply increase from the HAD and the consequent water salinity reduction in the delta region.The increase in water supply from the HAD was estimated by the BlueM model based on the GCMprojected average increase percentage of the Nile River flow at the Nasser Lake entrance(14.3%)(Table 1).The results of Scenarios 1 and 2 were compared to estimate the impacts of the Nile River flow and water salinity in the delta region on cropping pattern,food security,and the socioeconomy in Egypt.

    Scenario 3 evaluated the effectiveness of the cropping pattern adaptation measures as described in Section 2.2.5.The adaptation measures were utilized in the case of the reduction in Nile River flows.This scenario used the same reduction rate of the Nile River flow and growth rate of salinity as used in Scenario 1.The results of Scenarios 1 and 3 were compared to indicate the possible impacts of adaptation measures on cropping patterns,food security,and the socioeconomy in Egypt.

    3.Results

    3.1.Results of model calibration

    With regard to the BlueM model,the average monthly water level data upstream of the HAD in the period from 2013 to 2017 were used for model calibration.Fig.3 shows that the BlueM model provided a satisfactory water level simulation with an absolute error index ranging from-0.23 m to 0.17 m.As for the WB-SAL model,the relative error for the salinity of water drainage to the Mediterranean Sea and coastal lakes was 1.45%,indicating a strong simulation performance(Fig.4).The ASME model obtained a relative error of 3.38% for the volume of reused drainage water(Fig.4),and the modeling performance was satisfactory.

    3.2.Vulnerability of HAD water supply and salinity of water in delta region

    Climate change in the Nile River Basin and in Egypt will expose the water balance,water salinity,food security,and the socioeconomy in Egypt to vulnerability in 2050.Climate change will alter the flow entering Nasser Lake.Based on the BlueM model,the mean annual water supply from the HAD under drying climate scenarios would be 64.4 km3and 48.8 km3with an average inflow change rate of 14.3% and-11.8%,respectively.Based on the WB-SAL model,the salinity of drainage water to the sea would increase from its current level of 1 859 mg/L to 2 073 mg/L under Scenario 1 with a growth rate of 11.5%,and the salinity of irrigation water would rise from 361.7 mg/L to 490.2 mg/L with a growth rate of 35.5%.By contrast,the drainage salinity and irrigation salinity would decrease to 1 561 mg/L and 305.4 mg/L,respectively,in Scenario 2.Notably,the salinity changes in the delta region under future scenarios could be attributed to the alteration of the five factors:water supply from the Nile River,rate of seawater intrusion to shallow groundwater,population,evaporation,and agricultural area.Given that all these factors,except for the Nile River flow,were assumed to be the same under future scenarios,the water salinity in the delta region would be inversely proportional to the Nile River flow.

    Fig.3.Simulated mean monthly water level upstream of HAD in period from 2013 to 2017 using BlueM model.

    Fig.4.TDS of drainage water flowing to sea simulated by WB-SAL model and volume of reused drainage water simulated by ASME model.

    3.2.1.Vulnerability of cropping patterns

    The simulated results from the BlueM and WB-SAL models were integrated with the ASME model to project future changes in cropping patterns.For the entire country,crop area for each crop type would significantly change in the three future scenarios(Fig.5).In comparison with the current situation,the crop areas for 13 crop types would decrease in Scenario 1 and increase in Scenario 2.These crop types are faba bean,groundnut,long berseem,lentils,Nili maize,summer onion,soybeans and sunflower,short berseem,sugar cane,cotton,Nili sorghum,and Nili tomato.Only the crop area for sugar beet would rise in Scenario 1 and drop in Scenario 2.In both scenarios,the crop areas for barley,summer maize,and rice would decrease.For the eight crop types of orchard,other legumes,summer tomato,winter tomato,Nili vegetables,summer vegetables,winter vegetables,and wheat,the crop areas were projected to increase in both scenarios.The total crop area would decrease from the current amount of 5 516 700 hm2to 4 146 660 hm2in Scenario 1 and 5 907 300 hm2in Scenario 2.

    3.2.2.Vulnerability of food security

    The projected cropping patterns in future scenarios will directly impact future food security(Table 3).Cropping intensity,which refers to the fraction of cultivated area that is harvested,is an important indicator of food security.A cropping intensity value exceeding 100%indicates that more than one cropping cycle is permitted in a year in the same region.As shown in Table 3,the cropping intensity would decrease from its current level of 161.07%-124.85% in Scenario 1.By contrast,it would increase to 169.05% in Scenario 2.In Scenario 1,the self-sufficiencies for various crops were projected to decrease,and the decrease magnitude for wheat,rice,cereal,and maize were 7.56%,83.8%,19.55%,and 16.80%,respectively.The self-sufficiencies in Scenario 2 were higher than those in Scenario 1,with those for wheat,rice,cereal,and maize in Scenario 2 being 19.71%,19.31%,26.31%,and 27.7% higher,respectively,than those in Scenario 1.

    Fig.5.Cropping patterns in Egypt in current situation and three future scenarios.

    Table 3Food security at current level and in three future scenarios.

    3.2.3.Socioeconomic vulnerability

    Evaluation of the socioeconomic situation in the agricultural sector depends on many criteria,including agricultural productivity at the farm-gate price,consumer-producer surplus,number of jobs,and water marginal value(Table 4).The farm-gate value of an agricultural product is the net price of this product when it leaves the farm after marketing costs are subtracted.Compared with the current status,the net agricultural productivity at the farm-gate price was projected to decrease in Scenario 1.It would be 2.49×109USD higher in Scenario 2 than in Scenario 1.Consumer-producer surplus is an important indicator for assessment of the economic progress in the agricultural sector,and it was calculated as the difference between the revenue of total crops and the combined costs for crop harvesting and irrigation.The consumerproducer surplus in Scenario 2 was projected to be 1.086×1010USD higher than that in Scenario 1.Water marginal value measures the willingness of water users to pay for additional water when water demand is not sustained.Water marginal value is related to water scarcity,and it was reformulated into a scarcity cost indicator.The reduction in Nile River flows in Scenario 1 led to an increase in water marginal value,indicating that water would become scarcer.The marginal value in Scenario 2 was projected be 52.17%lower than that in Scenario 1.Meanwhile,the number of jobs for crop production,animal husbandry,and agro-processing were projected to be higher in Scenario 2 than in Scenario 1.

    3.3.Effectiveness of adaptation measures

    The cropping pattern adaptation measures in Scenario 3 included fixing rice cultivation areas in coastal governorates and fixing the cultivation areas of orchards,lentils,maize,onions,and vegetables in all governorates.The adaptation measures have been suggested in the case of the reduction in Nile River flows and water salinity increase.Hence,the values of water supply from the HAD and water salinity in the delta region in Scenario 3 were defined as same as those in Scenario 1.The comparison of Scenarios 1 and 3 helps to clarify the effectiveness of adaptation measures,because all other factors were the same in these two scenarios.In Scenario 3,the cropping areas of faba beans,long berseem,other legumes,soybeans,sunflowers,sugar beets,short berseem,cotton,and Nili sorghum were predicted to increase.However,those for groundnuts,sugar cane,and sesame were projected to decrease.The total crop area would increase from 4 146 660 hm2in Scenario 1 to 4 373 040 hm2in Scenario 3(Fig.5).

    With regard to food security,the adaptation measures in Scenario 3 increased the cropping intensity by 5.02% in comparison with those in Scenario 1(Table 3).These measures also augmented the self-sufficiencies of wheat,rice,cereal,and maize by 9.57%,22.23%,0.54%,and 1.20%,respectively(Table 3).Clearly,the cropping pattern adaptation measures had a significant impact on self-sufficiencies of rice and wheat but played a weak role in self-sufficiencies of cereal and maize.The self-sufficiencies of sugar and meat were projected to decrease in both Scenarios 1 and 3,and milk production would fully sustain the demand in these two scenarios.

    Although the cropping pattern adaptation measures increased the total cropped area and improved food security in comparison with Scenario 1,they did not influence the socioeconomic situation.As shown in Table 4,the current net agricultural productivity at the farm-gate price,consumerproducer surplus,water marginal value,and number of jobs for crop production,animal husbandry,and agro-processing in both Scenarios 1 and 3 had similar values.

    4.Discussion

    This study evaluated the vulnerability of the water supply from the HAD,water salinity in the Nile Delta region,foodsecurity,and the socioeconomy in Egypt in response to climate change by the year 2050.The mean annual water supply from HAD release showed a large uncertainty,either with an average increase of 14.30% or a decrease of-11.80%.The salinity of irrigation water in the Nile Delta would increase by 35.5% in the case of the reduction in Nile River flows.

    Table 4Socioeconomic situation currently and in three future scenarios.

    It was concluded that crop area for each crop type throughout country would be vulnerable to climate change in 2050.A clear and direct relationship existed between Nile River flows and crop areas for crop types such as faba beans,groundnuts,long berseem,lentils,Nili maize,summer onions,soyabeans,sunflowers,short berseem,sugar cane,cotton,Nili sorghum,and Nili tomatoes.Meanwhile,only the sugar beet area was inversely proportional to the Nile River flow.The Nile River flow had an insignificant relationship with the crop areas for barley,summer maize,rice,orchards,other legumes,summer tomatoes,winter tomatoes,Nili vegetables,summer vegetables,winter vegetables,and wheat.In general,the projected total crop area was proportional to the Nile River flow.

    This study found that the reduction in Nile River flows and the increased water salinity in the delta region tended to reduce the self-sufficiencies of wheat,rice,cereal,and maize,and vice versa.Therefore,food security in Egypt would be vulnerable to climate change.In addition,this study quantified the net agricultural productivity at the farm-gate price,consumer-producer surplus,water marginal value,and employment for crop production,animal husbandry,agroprocessing,and the agricultural sector.It was found that the reduced Nile River flow and the increased salinity in the delta region would deteriorate the socioeconomy of the agricultural sector.

    This study proposed a set of cropping pattern adaptation measures,including fixing rice crop areas in coastal governorates and fixing the cropping areas of orchards,lentils,maize,onions,and vegetables in all governorates.These adaptation measures would be enacted only in the case of the reduction in the Nile River flow.The measures would result in an increase in the crop areas of faba beans,long berseem,other legumes,soybeans,sunflower,sugar beets,short berseem,cotton,and Nili sorghum but would lead to a decrease in the crop areas of groundnuts,sugar cane,and sesame.As a result,the total crop area would increase.The adaptation measures would be effective in improving the food security status in the country owing to the increase of cropping intensity and the self-sufficiencies of wheat,rice,cereal,and maize.This study found that when these measures were adopted,the net agricultural productivity at the farm gateprice,consumer-producer surplus,water marginal value,and employment for crop production,animal husbandry,and agroprocessing would not significantly change in comparison with the situation in which no measures were considered.Therefore,the improved food security status would come at the expense of socioeconomic status.

    These results agree with the findings of Chebil et al.(2019),who adopted the agricultural supply model for Tunisia(ASMOT)using the data for 21 strategic crops.They found that the total crop area and agricultural employment in Tunisia were negatively affected due to the decreased irrigation water availability arising from climate change.The findings of this study are also in line with those of Adhikari(2018),who reported that 15 drought events took place in Nepal from 1972 to 2015 in climate change conditions.The adaptation measures included abandoning the crop types of rice,cereal crops,tubers,and sweet potatoes,which negatively affected agricultural production and the socioeconomic situation.

    Comparison of the results in the three future scenarios clearly demonstrated that the Nile River flow would be largely impacted by climate change,which would subsequently affect the water salinity in the delta region,food security,and the socioeconomy.The Nile River flow reduction of 11.80% at the Dongola Station at the Nasser Lake entrance would reduce the average self-sufficiency of strategic crops by 31.08% and net agricultural productivity by 2.316×109USD,and lead to a loss in daily part-time income for 206 000 laborers.The reduction in the Nile River flow combined with cropping pattern adaptation measures would lower the average self-sufficiency of strategic crops by 21.29% and net agricultural productivity by 2.357×109USD,and result in a loss in daily part-time income for 203 000 laborers.However,the Nile River flow increase of 14.30% would augment the average selfsufficiency of strategic crops by 6.06% and net agricultural productivity by 1.73×108USD.It would lead to a gain in daily part-time income for 23 000 laborers.Therefore,it can be concluded that an increase in the Nile River flow would have more positive effects on food security and the socioeconomy than the cropping pattern adaptation measures.

    5.Conclusions

    In this study,the impacts of climate change on food security and the socioeconomy in Egypt were investigated,and the effectiveness of cropping pattern adaptation measures was evaluated by integrating three mathematical models.The main conclusions of this study can be summarized as follows:

    (1)The vulnerability analysis in this study demonstrated that climate change would significantly affect the Nile River flow in Egypt.The Nile River flow would have an inverse relationship with the water salinity in the delta region.A reduction in the Nile River flow would reduce the average selfsufficiency of strategic crops,net agricultural productivity,and consumer-producer surplus.It would increase the number of laborers who lose their daily part-time income.In contrast,an increase in the Nile River flow would augment the average self-sufficiency of strategic crops,net agricultural productivity,and consumer-producer surplus and reduce the number of laborers with decreased daily part-time income.

    (2)This study proposed to maintain a fixed crop area for rice in nine coastal governorates and to maintain fixed crop areas for orchards,lentils,maize,onion,tomatoes,and vegetables.These cropping pattern adaptation measures would increase the total crop area and improve food security without causing a deterioration of socioeconomic indicators.(3)Water availability would have more significant impacts on food security and socioeconomic status than cropping pattern adaptation measures.Accordingly,the country should increase water use efficiency and water supply to preserve agricultural productivity,food security,and socioeconomic health under different climate change scenarios.

    Declaration of competing interest

    The authors declare no conflicts of interest.

    黄片wwwwww| 最新中文字幕久久久久| 成人午夜精彩视频在线观看| 国产精品一区www在线观看| 夫妻午夜视频| 国产又色又爽无遮挡免| 色婷婷久久久亚洲欧美| 伦理电影免费视频| 在线看a的网站| 亚洲av二区三区四区| 在线观看一区二区三区激情| 春色校园在线视频观看| 欧美一区二区亚洲| 精品国产三级普通话版| 女性生殖器流出的白浆| 免费看日本二区| 国产精品熟女久久久久浪| 亚洲精品日韩在线中文字幕| 小蜜桃在线观看免费完整版高清| 秋霞伦理黄片| 国产精品国产三级国产专区5o| 久久99热这里只频精品6学生| 亚洲精品视频女| 国产午夜精品一二区理论片| 亚洲精品中文字幕在线视频 | 七月丁香在线播放| 久久热精品热| 亚洲第一区二区三区不卡| 熟女电影av网| 国产高潮美女av| 在线看a的网站| 国产av精品麻豆| 色综合色国产| 精品国产露脸久久av麻豆| 国产精品麻豆人妻色哟哟久久| 精品久久久久久久末码| 精品久久久久久电影网| 亚洲av日韩在线播放| 国产午夜精品一二区理论片| 如何舔出高潮| 国语对白做爰xxxⅹ性视频网站| 午夜激情久久久久久久| 一区二区三区精品91| 国产成人午夜福利电影在线观看| 欧美丝袜亚洲另类| 欧美一区二区亚洲| 久久久久久伊人网av| 久久久久久九九精品二区国产| 国产成人一区二区在线| 美女中出高潮动态图| 亚洲精品成人av观看孕妇| 国产精品久久久久久av不卡| 少妇的逼好多水| 亚洲av成人精品一二三区| av在线蜜桃| 女性生殖器流出的白浆| 国产淫语在线视频| 免费看日本二区| 最近的中文字幕免费完整| 欧美xxxx黑人xx丫x性爽| 丰满乱子伦码专区| av免费在线看不卡| 国产高潮美女av| 99精国产麻豆久久婷婷| 日韩一区二区视频免费看| 国产亚洲av片在线观看秒播厂| 久久久久性生活片| 高清欧美精品videossex| 国产高清国产精品国产三级 | 人妻系列 视频| 2021少妇久久久久久久久久久| 美女xxoo啪啪120秒动态图| 五月玫瑰六月丁香| av免费在线看不卡| 久久精品国产自在天天线| 亚洲一区二区三区欧美精品| 成人亚洲欧美一区二区av| 精品一区二区三区视频在线| 校园人妻丝袜中文字幕| 国产精品精品国产色婷婷| 久久人人爽人人爽人人片va| 18禁裸乳无遮挡动漫免费视频| 国产永久视频网站| 久久精品国产亚洲av涩爱| 少妇 在线观看| 久久 成人 亚洲| 亚洲美女视频黄频| 国产爱豆传媒在线观看| 黄色一级大片看看| tube8黄色片| 欧美日韩国产mv在线观看视频 | 亚洲精品国产av蜜桃| 丰满少妇做爰视频| 亚洲国产欧美在线一区| 欧美日韩视频精品一区| 1000部很黄的大片| 80岁老熟妇乱子伦牲交| 中文在线观看免费www的网站| 午夜激情福利司机影院| 久久韩国三级中文字幕| 大片电影免费在线观看免费| 国产精品偷伦视频观看了| 18+在线观看网站| 最近最新中文字幕大全电影3| 少妇熟女欧美另类| 欧美xxxx性猛交bbbb| 热99国产精品久久久久久7| 噜噜噜噜噜久久久久久91| 亚洲精品乱久久久久久| 亚洲一区二区三区欧美精品| 日本免费在线观看一区| 国产高清不卡午夜福利| 国产精品一区二区在线观看99| 亚洲色图av天堂| 最近中文字幕2019免费版| 日韩,欧美,国产一区二区三区| 免费看不卡的av| 国产伦理片在线播放av一区| 午夜视频国产福利| 熟女av电影| 亚洲三级黄色毛片| 97在线人人人人妻| 少妇熟女欧美另类| 少妇 在线观看| 日本黄色片子视频| 一级黄片播放器| 久久久色成人| 青青草视频在线视频观看| 亚洲内射少妇av| 99热网站在线观看| 狂野欧美激情性xxxx在线观看| 国模一区二区三区四区视频| 一级黄片播放器| av在线app专区| 在线观看一区二区三区激情| 国产又色又爽无遮挡免| 人妻制服诱惑在线中文字幕| 日本-黄色视频高清免费观看| 免费人妻精品一区二区三区视频| 哪个播放器可以免费观看大片| 最近最新中文字幕免费大全7| 爱豆传媒免费全集在线观看| 精品一区二区免费观看| 少妇人妻 视频| 国产精品成人在线| 亚洲av中文av极速乱| 亚洲精品自拍成人| 亚洲美女黄色视频免费看| 亚洲欧美一区二区三区黑人 | 在线观看三级黄色| 日韩不卡一区二区三区视频在线| 成人手机av| 嫁个100分男人电影在线观看 | 亚洲精品av麻豆狂野| 中文字幕最新亚洲高清| 精品一区在线观看国产| 黑人猛操日本美女一级片| 男女边吃奶边做爰视频| 人人妻人人添人人爽欧美一区卜| 色视频在线一区二区三区| 丝袜美足系列| 国产1区2区3区精品| 天堂8中文在线网| 久久人人爽av亚洲精品天堂| 青青草视频在线视频观看| 黄色毛片三级朝国网站| 黄色怎么调成土黄色| 丝袜美足系列| 午夜激情久久久久久久| 国产成人av教育| 国产欧美日韩综合在线一区二区| 91老司机精品| 久久国产精品人妻蜜桃| 亚洲国产精品一区三区| 91精品伊人久久大香线蕉| 欧美成人午夜精品| 在线亚洲精品国产二区图片欧美| 国产亚洲一区二区精品| 国产精品一二三区在线看| 晚上一个人看的免费电影| 亚洲 国产 在线| 国产日韩欧美视频二区| 天堂俺去俺来也www色官网| 亚洲五月色婷婷综合| 亚洲精品国产区一区二| 午夜两性在线视频| 悠悠久久av| 国产成人系列免费观看| 亚洲视频免费观看视频| 欧美日本中文国产一区发布| 国产精品国产三级专区第一集| 国产激情久久老熟女| 亚洲情色 制服丝袜| 欧美 亚洲 国产 日韩一| 欧美亚洲日本最大视频资源| 国产在线一区二区三区精| 欧美在线黄色| 人妻 亚洲 视频| 99精国产麻豆久久婷婷| 亚洲国产欧美一区二区综合| 日韩熟女老妇一区二区性免费视频| 性高湖久久久久久久久免费观看| 婷婷色综合大香蕉| 国产激情久久老熟女| www.精华液| xxxhd国产人妻xxx| 天堂8中文在线网| 欧美乱码精品一区二区三区| www日本在线高清视频| 欧美成人午夜精品| 无遮挡黄片免费观看| 一二三四社区在线视频社区8| 精品免费久久久久久久清纯 | 国产精品偷伦视频观看了| 女人精品久久久久毛片| 看免费av毛片| 国产精品一区二区精品视频观看| 黑人巨大精品欧美一区二区蜜桃| kizo精华| 大型av网站在线播放| 两性夫妻黄色片| 亚洲欧洲国产日韩| 制服人妻中文乱码| 欧美日韩国产mv在线观看视频| 曰老女人黄片| 亚洲国产欧美日韩在线播放| av在线老鸭窝| 国产主播在线观看一区二区 | xxx大片免费视频| 国产精品一区二区在线不卡| 在线观看免费高清a一片| 欧美国产精品va在线观看不卡| 女性被躁到高潮视频| 晚上一个人看的免费电影| av在线老鸭窝| 尾随美女入室| 美女脱内裤让男人舔精品视频| 精品福利永久在线观看| 婷婷色麻豆天堂久久| 免费在线观看黄色视频的| 一区二区av电影网| 国产亚洲欧美精品永久| 久久热在线av| 久久 成人 亚洲| 亚洲国产精品成人久久小说| 婷婷色av中文字幕| 涩涩av久久男人的天堂| 久久久欧美国产精品| www.av在线官网国产| 久久精品国产亚洲av高清一级| 日韩大码丰满熟妇| 少妇人妻久久综合中文| 久久人人爽av亚洲精品天堂| 在线 av 中文字幕| 乱人伦中国视频| 黄色一级大片看看| 18在线观看网站| 狂野欧美激情性xxxx| 免费在线观看日本一区| 曰老女人黄片| 国产一区有黄有色的免费视频| 国产一区二区激情短视频 | 国产成人精品久久二区二区91| 亚洲美女黄色视频免费看| 色综合欧美亚洲国产小说| 亚洲情色 制服丝袜| 日韩 亚洲 欧美在线| 丝袜喷水一区| 国产一卡二卡三卡精品| 高清视频免费观看一区二区| 欧美日韩综合久久久久久| 中文欧美无线码| 亚洲国产精品一区三区| 欧美精品啪啪一区二区三区 | 日韩制服丝袜自拍偷拍| 国产欧美日韩一区二区三区在线| 亚洲第一av免费看| 久久精品人人爽人人爽视色| 国产精品国产三级专区第一集| 亚洲综合色网址| 好男人视频免费观看在线| 丰满饥渴人妻一区二区三| 国产成人系列免费观看| 国产亚洲欧美在线一区二区| 视频在线观看一区二区三区| 中文精品一卡2卡3卡4更新| 午夜福利影视在线免费观看| 少妇粗大呻吟视频| 赤兔流量卡办理| 国产精品久久久久久人妻精品电影 | 日本五十路高清| 人妻一区二区av| 亚洲欧美精品自产自拍| 国产成人91sexporn| 校园人妻丝袜中文字幕| 国产日韩欧美视频二区| 午夜免费男女啪啪视频观看| 男女下面插进去视频免费观看| 人人妻人人澡人人看| 校园人妻丝袜中文字幕| 久久精品久久久久久久性| 国产精品国产av在线观看| 性色av乱码一区二区三区2| 中文字幕人妻熟女乱码| 一边亲一边摸免费视频| 丝袜脚勾引网站| 国产又色又爽无遮挡免| 亚洲国产精品一区二区三区在线| 国产成人91sexporn| 在现免费观看毛片| 午夜福利一区二区在线看| 日本a在线网址| 爱豆传媒免费全集在线观看| 精品卡一卡二卡四卡免费| 18禁观看日本| 精品久久久久久久毛片微露脸 | 欧美黑人精品巨大| 夫妻性生交免费视频一级片| 视频区图区小说| 9191精品国产免费久久| 999精品在线视频| 中文字幕人妻丝袜一区二区| 国产精品久久久久久人妻精品电影 | 少妇猛男粗大的猛烈进出视频| 亚洲精品第二区| 国产精品久久久av美女十八| 亚洲久久久国产精品| 日韩av不卡免费在线播放| av网站在线播放免费| 亚洲色图综合在线观看| 老汉色∧v一级毛片| 国产免费又黄又爽又色| 亚洲一码二码三码区别大吗| 亚洲欧美成人综合另类久久久| 黄色毛片三级朝国网站| 成人国产一区最新在线观看 | 制服人妻中文乱码| 亚洲熟女毛片儿| 国产又爽黄色视频| 男女下面插进去视频免费观看| 亚洲国产av新网站| 人妻人人澡人人爽人人| 国产亚洲av片在线观看秒播厂| 男女免费视频国产| 777米奇影视久久| 一级a爱视频在线免费观看| 在线观看免费视频网站a站| cao死你这个sao货| 国产成人欧美在线观看 | 国产老妇伦熟女老妇高清| 精品国产一区二区三区四区第35| 亚洲国产日韩一区二区| 日韩中文字幕视频在线看片| 久久 成人 亚洲| 午夜免费成人在线视频| 国产精品二区激情视频| 一区二区三区精品91| 国产一区亚洲一区在线观看| 国产精品久久久久久精品电影小说| 美女中出高潮动态图| 中国美女看黄片| 国产一区二区在线观看av| 亚洲欧洲精品一区二区精品久久久| 黄片播放在线免费| 国产日韩欧美视频二区| 一区福利在线观看| 欧美激情高清一区二区三区| 精品熟女少妇八av免费久了| 久久青草综合色| 黑人巨大精品欧美一区二区蜜桃| 久久久久精品人妻al黑| 国产色视频综合| 黄色怎么调成土黄色| 久久人人97超碰香蕉20202| 亚洲久久久国产精品| 久久九九热精品免费| 久久热在线av| 亚洲精品成人av观看孕妇| 日韩视频在线欧美| av网站免费在线观看视频| 亚洲国产中文字幕在线视频| 欧美精品亚洲一区二区| 亚洲男人天堂网一区| 精品熟女少妇八av免费久了| av国产精品久久久久影院| 18禁观看日本| 99国产精品一区二区蜜桃av | 国产激情久久老熟女| 天天躁日日躁夜夜躁夜夜| 国产精品久久久久久精品电影小说| 久久天堂一区二区三区四区| 只有这里有精品99| 美女福利国产在线| 欧美97在线视频| 成年人黄色毛片网站| 久久久久久久精品精品| 欧美日韩亚洲高清精品| 夫妻性生交免费视频一级片| 老司机午夜十八禁免费视频| 啦啦啦中文免费视频观看日本| 飞空精品影院首页| 美女高潮到喷水免费观看| 99香蕉大伊视频| 国精品久久久久久国模美| 欧美在线黄色| 亚洲色图 男人天堂 中文字幕| 99国产精品99久久久久| 99久久人妻综合| 老鸭窝网址在线观看| 自拍欧美九色日韩亚洲蝌蚪91| 色婷婷久久久亚洲欧美| 久久精品熟女亚洲av麻豆精品| 十八禁人妻一区二区| 国产视频首页在线观看| 少妇粗大呻吟视频| 亚洲成色77777| 国产精品一区二区在线不卡| 极品少妇高潮喷水抽搐| 国产黄色免费在线视频| 人妻人人澡人人爽人人| 亚洲成av片中文字幕在线观看| 成人亚洲欧美一区二区av| 免费av中文字幕在线| 精品亚洲乱码少妇综合久久| 精品福利永久在线观看| 天堂8中文在线网| 精品一区二区三卡| bbb黄色大片| www日本在线高清视频| h视频一区二区三区| 免费在线观看黄色视频的| 美国免费a级毛片| 久久国产精品男人的天堂亚洲| 免费女性裸体啪啪无遮挡网站| 成人国语在线视频| 亚洲国产精品成人久久小说| 人妻 亚洲 视频| 日本色播在线视频| 亚洲欧美清纯卡通| 午夜福利,免费看| 日韩伦理黄色片| 婷婷色综合www| 国产淫语在线视频| 十八禁人妻一区二区| 人人澡人人妻人| 欧美日韩亚洲高清精品| 免费观看人在逋| 免费不卡黄色视频| 亚洲国产最新在线播放| 国产一区有黄有色的免费视频| 成年av动漫网址| 97在线人人人人妻| 电影成人av| 人人澡人人妻人| 男人添女人高潮全过程视频| 免费观看a级毛片全部| 精品高清国产在线一区| 亚洲图色成人| 亚洲中文av在线| av片东京热男人的天堂| 一级毛片女人18水好多 | 精品人妻熟女毛片av久久网站| 在线av久久热| 亚洲精品一卡2卡三卡4卡5卡 | 超碰97精品在线观看| 天堂中文最新版在线下载| 丰满迷人的少妇在线观看| 国产在线观看jvid| 如日韩欧美国产精品一区二区三区| 大片免费播放器 马上看| 1024香蕉在线观看| 两性夫妻黄色片| 99九九在线精品视频| 91麻豆精品激情在线观看国产 | 少妇精品久久久久久久| 亚洲五月婷婷丁香| 亚洲欧美一区二区三区久久| 日韩精品免费视频一区二区三区| 欧美日韩一级在线毛片| 国产精品熟女久久久久浪| 五月开心婷婷网| av天堂久久9| 国产精品麻豆人妻色哟哟久久| 免费在线观看黄色视频的| 成年人午夜在线观看视频| 成人免费观看视频高清| 欧美大码av| 乱人伦中国视频| 亚洲五月婷婷丁香| 好男人视频免费观看在线| 亚洲人成网站在线观看播放| 国产精品免费大片| 老熟女久久久| 中文欧美无线码| 手机成人av网站| 女性生殖器流出的白浆| 欧美另类一区| 久久久国产欧美日韩av| 日本av免费视频播放| 一级a爱视频在线免费观看| 日韩,欧美,国产一区二区三区| 黑人欧美特级aaaaaa片| 国产精品香港三级国产av潘金莲 | 精品国产超薄肉色丝袜足j| 你懂的网址亚洲精品在线观看| 午夜福利视频精品| 亚洲国产日韩一区二区| 人人妻人人添人人爽欧美一区卜| 亚洲视频免费观看视频| e午夜精品久久久久久久| 婷婷色av中文字幕| 免费高清在线观看视频在线观看| 国产精品久久久久久精品电影小说| 肉色欧美久久久久久久蜜桃| 日本av手机在线免费观看| 午夜精品国产一区二区电影| 婷婷色麻豆天堂久久| 美女国产高潮福利片在线看| 精品福利永久在线观看| 国产精品久久久久成人av| 看十八女毛片水多多多| 亚洲av电影在线观看一区二区三区| 国产又爽黄色视频| 观看av在线不卡| 亚洲三区欧美一区| 欧美激情高清一区二区三区| 人人澡人人妻人| 一区福利在线观看| bbb黄色大片| 亚洲五月色婷婷综合| 51午夜福利影视在线观看| 国产精品一区二区免费欧美 | 欧美日韩av久久| 国产高清不卡午夜福利| 黑人欧美特级aaaaaa片| 久热这里只有精品99| e午夜精品久久久久久久| 侵犯人妻中文字幕一二三四区| 午夜福利在线免费观看网站| 免费久久久久久久精品成人欧美视频| 久久久久久人人人人人| 伊人亚洲综合成人网| 无限看片的www在线观看| 99国产精品免费福利视频| 一级毛片电影观看| 国产成人系列免费观看| 日本av免费视频播放| 久热爱精品视频在线9| 亚洲男人天堂网一区| 一级片'在线观看视频| 亚洲av成人不卡在线观看播放网 | 国产亚洲精品久久久久5区| 久久青草综合色| 国产又色又爽无遮挡免| 熟女少妇亚洲综合色aaa.| 免费在线观看影片大全网站 | 日本猛色少妇xxxxx猛交久久| 一个人免费看片子| 日本vs欧美在线观看视频| 天天躁夜夜躁狠狠躁躁| 只有这里有精品99| 亚洲精品成人av观看孕妇| 婷婷色麻豆天堂久久| 亚洲欧美激情在线| 伊人亚洲综合成人网| 在线 av 中文字幕| 大型av网站在线播放| 91九色精品人成在线观看| 黄色视频不卡| 一级毛片女人18水好多 | www.熟女人妻精品国产| 国产精品久久久av美女十八| 精品国产一区二区久久| 男人舔女人的私密视频| 欧美精品一区二区免费开放| 国产精品一二三区在线看| 亚洲av国产av综合av卡| 一边摸一边做爽爽视频免费| 久久精品亚洲av国产电影网| 久久久久国产精品人妻一区二区| 亚洲av日韩在线播放| a 毛片基地| 亚洲精品自拍成人| 一级黄片播放器| 国产黄频视频在线观看| 女人被躁到高潮嗷嗷叫费观| 欧美日韩一级在线毛片| 18禁观看日本| 国产一区二区在线观看av| 丁香六月天网| 大码成人一级视频| 亚洲精品久久成人aⅴ小说| av又黄又爽大尺度在线免费看| 欧美精品高潮呻吟av久久| 色综合欧美亚洲国产小说| 欧美av亚洲av综合av国产av| 操美女的视频在线观看| 国产成人欧美| 后天国语完整版免费观看| 丰满饥渴人妻一区二区三| 波多野结衣一区麻豆| 国产精品久久久久久精品古装| 99热国产这里只有精品6| 99久久人妻综合| 欧美黄色淫秽网站| 精品亚洲成a人片在线观看| 亚洲成国产人片在线观看| 国产日韩欧美亚洲二区| 国产精品人妻久久久影院| 两个人免费观看高清视频| 久久热在线av| 亚洲国产精品999| 日日夜夜操网爽| 国产在线免费精品| 国产成人a∨麻豆精品|