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

    Using Multidisciplinary Analysis to Develop Adaptation Options against Extreme Coastal Floods

    2022-12-09 03:21:14XinmengShanJunWangJiahongWenHengzhiHuLeiWangJieYinMengyaLi

    Xinmeng Shan · Jun Wang · Jiahong Wen · Hengzhi Hu · Lei Wang · Jie Yin · Mengya Li

    Abstract Long-term flood risk adaptation and decision making are complex because the future is full of deep uncertainties.Flexibility and robustness can be used to deal with future uncertainty.This study developed an integrated modeling framework that extends previous studies to the spatial domain to assess the future flood risks and the cost and benefit of three adaptation measures for four types of buildings in Shanghai.Real options analysis (ROA) and dynamic adaptive policy pathways (DAPP) were integrated to develop a dynamic adaptation pathway and identify robust adaptation options.The results show that: (1) Sea level rise and land subsidence will significantly exacerbate the flood risks in Shanghai; (2) Among the three flood control measures, wet-floodproofing has the best economic performance in terms of both the net present value and the benefit/cost ratio, followed by dry-floodproofing, and elevation; (3) Dryfloodproofing can be used at the beginning of the future period (2030—2100), and it can be replaced by wet-floodproofing in 2035—2042; the elevation measure also shows good performance at the beginning of implementation, but its performance will decline after 2041—2045; (4) The combined strategy of dry-and wet-floodproofing in 2044—2046 and a hybrid strategy combining the three measures should be the optimal solution for reducing the flood risks in 2047—2051.The methodology developed in this study can provide insights for coastal cities to formulate cost-effective and feasible adaptation strategies in a deeply uncertain future.

    Keywords Adaptation options · Building-level measures · Climate change · Coastal floods · Cost-benefit analysis · Shanghai

    1 Introduction

    The Intergovernmental Panel on Climate Change’s Sixth Assessment Report shows that global sea level will continue to rise throughout the twenty-f irst century, leading to more frequent and severe f loods in coastal areas (IPCC 2021).The combined risk of extreme f loods due to the sea level rise, land subsidence, and climate change is threatening coastal cities (Abadie et al.2020; Catalao et al.2020; Strauss et al.2021).Under Representative Concentration Pathway (RCP) 2.6 and RCP8.5 carbon emission scenarios, the global comprehensive risk will increase two-and fourfold, respectively (Magnan et al.2021).Risk analysis of extreme storm f loods is crucial for coastal cities to develop strategic plans and adaptation measures (Scussolini et al.2017; Du et al.2020; Wang et al.2020).In particular, extreme f loods with low probability and high magnitude have become the focus of recent studies because of their catastrophic consequences and losses (Shan et al.2021; Tang et al.2021).Many studies have assessed f lood risks in coastal cities in the context of sea level rise, land subsidence, and climate change (Lin et al.2016; Hinkel et al.2018; Shan et al.2019; Ngo et al.2020; Nguyen et al.2021).Due to sea level rise the 1/500-year flood in New York City before the Industrial Revolution, for example, has become the current 1/25-year flood (Reed et al.2015; Garner et al.2017).Hence, how coastal cities adapt to the increasing risk of extreme storm flooding needs extensive and insightful discussions (Barnett et al.2014; Haasnoot et al.2021).

    Flood adaptation refers to the implementation of a series of measures to reduce the occurrence and potential impact of flood events (Tanoue et al.2021).Cost-benefit analysis (CBA) is commonly used to evaluate the effectiveness of flood adaptation measures, and is an effective tool to support making adaptation decisions (Ward et al.2017; de Ruig et al.2019).Many studies have examined the cost and benefit of flood adaptation measures in coastal cities.Studies for Ho Chi Minh City (Scussolini et al.2017), Jakarta (Ward et al.2011), New York (Aerts et al.2014), Los Angeles (de Ruig et al.2020), London (Hall et al.2019), and Osaka Bay (Ha et al.2021) have focused on the impact of coastal flooding on people and assets, and the cost effectiveness of adaptation measures under different sea level rise scenarios.In addition, some studies have focused on the cost-benefit analysis of flood adaptation measures at the household level (Wang et al.2015; Qin and Stewart 2020; Shan et al.2021).For example, Shan et al.( 2021) evaluated the cost-effectiveness of two flood adaptation measures—dry flood-proofing and wet flood-proofing—for flood adaptation of residential buildings in Shanghai.Qin and Stewart ( 2020) assessed the cost-effectiveness of four adaptation measures for reducing the risk of residential buildings.However, in the research of extreme flood risk adaptation, uncertainty exists not only in disaster risk assessment, but also in cost-benefit analysis of adaptation measures, as well as in the decision-making process.Therefore, extreme flood risk adaptation carries deep uncertainty.

    Due to the deep uncertainty of the changes, decisions need to be formulated as flexible and adaptable.Robust decision-making methods such as real options analysis (ROA) and dynamic adaptive policy pathways (DAPP) provide solutions to decision making that involve deep uncertainty and can be applied to research on decision making and planning to deal with climate change.

    Real options analysis is a method to deal with the uncertainty of investment timing by considering management flexibility and volatility (Erfani et al.2018; Kim et al.2019).Decision makers can decide when to invest and maintain the possibility of greater investment in the future using ROA.In recent years, the concept of ROA has been used in flood risk management to help decision makers formulate robust decisions to enhance flood adaptation in the context of climate change.Adaptation can be achieved because ROA provides an optimal sequence of future investment decisions to cope with uncertainty in flood hazards over time (Kind et al.2018).Flexible options can give more time to manage areas where investments may be overwhelmed or future risks are more severe than expected (Wreford et al.2020; Ginbo et al.2021).Kim et al.( 2019) evaluated the adaptive investment timing of Lymington, a flood-prone coastal town in the UK, under the scenario of sea level rise based on the ROA method.Liu et al.( 2018) used ROA to assess adaptation options for flood risks under climate change scenarios in London, and pointed out that this method was more flexible than fixed adaptation (that is, maintaining the current adaptation measures).Hino and Hall ( 2017) showed that a flexible adaptation plan can better deal with deep uncertainty than implementing a single plan at present.Similarly, postponing the timing of decision making under uncertain conditions allows for learning more from the change and making better decisions in the future (Woodward et al.2014).

    The dynamic adaptive policy pathways (DAPP) approach is a method that explores alternative sequences of investment decisions to achieve goals over time in the context of uncertain future development and environmental changes (Haasnoot et al.2013a; Haasnoot et al.2013b).Haasnoot et al.( 2020) proposed an economic evaluation framework to explore adaptation pathways or a series of strategic investment options and took flood risk management in the Netherlands as a case study.Monitoring and reevaluation of the adaptation strategies should be an ongoing process, and investment tipping points should be identified in time to transit or make changes to the adaptation strategies accordingly.

    The DAPP and ROA approaches have commonalities in adjusting options or strategies depending on the realization in the future.At present, these approaches are increasingly used to design or plan long-term adaptation options over multiple decades and longer terms (Woodward et al.2014; Kim et al.2019).Gersonius et al.( 2015) compared ROA (Gersonius et al.2013) with the tipping point analysis (the core idea of DAPP) (Kwadijk et al.2010) and showed that the two methods could be combined.Buurman and Babovic ( 2016) developed a conceptual framework to deal with deep uncertainty in climate change adaptation policy making based on ROA and DAPP, and demonstrated the framework to policy making for pluvial flood protection in Singapore.The combination of ROA and DAPP can be used to investigate the optimal timing of flood adaptation measures under deep uncertainty, and develop effective adaptation pathways for both the short term and the long term.

    A wide range of economic and robust decision-making methods is available to support decision making on coastal flood protection options under changing flood risks for effective adaptation, such as CBA (Aerts et al.2014), ROA (Kim et al.2019), and DAPP (Haasnoot et al.2020).However, few studies have demonstrated how to combine several types of such methods to generate complementary cost-effectiveness information to support coastal flood control decisions.Concerning the research gaps, based on a case study in Shanghai, this study aimed to develop a multidisciplinary analysis framework that encompasses risk analysis, cost-benefit analysis of adaptation measures, and adaptation pathways study using ROA and DAPP.The specific scientific questions to be examined with the data analyses are:

    (1) Will there be a significant increase in the future flood risk of buildings in the case study area?

    (2) What adaptation measures should be taken to reduce storm flood risks and how can the costs and benefits of measures be quantified?

    (3) How can dynamic optimal pathways for the short, medium, and long terms under deep uncertainty be formulated to help decision makers make robust decisions? The methodology to be developed should also provide insights into flood risk assessment and management in other coastal cities.

    2 Data and Methods

    This study used a set of multidisciplinary and comprehensive methods—flood risk analysis, cost-benefit analysis of adaptation measures, and adaptation pathways—to assess the flood risk dynamics, the effect of adaptation measures on reducing risks, and the formulation of dynamic adaptive policy pathways for buildings in 2010, 2030, 2050, and 2100 (Fig.1).

    First, hydrodynamic modeling was conducted to produce flood hazard maps associated with the storm surges in the baseline year of 2010 and future years of 2030, 2050, and 2100.The asset value of buildings was estimated and mapped according to the area, number of floors, and cost of each building type.The baseline year and future extreme storm flood simulation results were overlaid respectively with the asset values of buildings in 2010, 2030, 2050, and 2100 to evaluate the exposed assets.Building losses were obtained by combining exposed assets and vulnerability curves.The risk was quantified as estimated expected annual damage (EAD) (see more details in Sect.2.1.1).

    Second, cost-benefit analysis of adaptation measures was performed through net present value and benefit/cost ratio.This study evaluated the cost and benefit of the buildingcode strategies for individual buildings, including dry-floodproofing, wet-floodproofing, and elevation in Shanghai under the current and two representative concentration pathways (RCP2.6 and RCP8.5) scenarios.Dry-floodproofing measures mean that the flood-prone parts of a building are waterproofed, so that floodwater is prevented from entering inside.Wetfloodproofing measures allow floodwater to enter a building and its contents through adjusting the interior and material of the building.Elevation measures refer to the use of materials such as concrete to raise the height of the steps of the building.

    Third, on the basis of evaluating the performance of flood adaptation measures, by integrating the ROA of the continuation value and termination value of flood adaptation measures and the tipping point analysis of DAPP, the flood dynamic adaptation pathways of short, medium, and long terms were developed.

    2.1 Methods

    This section outlines the methodology employed in the study.The methods specifically include flood risk analysis, costbenefit analysis, real options analysis (ROA), and dynamic adaptive policy pathways (DAPP).

    2.1.1 Flood Risk Analysis

    The risk analysis mainly includes four parts: (1) The total value of the elements at risk can be calculated from the total asset value; (2) Exposure is measured by overlaying flood simulation results with the total assets; (3) The depth-damage functions for Shanghai were developed by Ke ( 2014)—they represent the relationship between flood characteristics and damage rate, and reflect the relationship between inundation depth below 0.5 m to above 3.0 m at an interval of 50 cm and the damage rates of residential, commercial, public, and industrial buildings (Fig.2).The losses can be evaluated by Eq.( 1); (4) Risk can be represented by a curve, which combines the losses under all return periods to obtain the expected annual loss (EAD), and EAD can be calculated using Eq.( 2).

    whereLis the losses of the elements at risk;Eis the exposed asset (Eq.1); andVis the loss rate of the elements at risk under different inundation depths.

    wherexis the losses;f(x)is the corresponding probabilities; and EAD can be used as the basis for flood adaptation costbenefit analysis.

    2.1.2 Cost-Benefit Analysis

    The adaptation measures are economically evaluated through a cost-benefit analysis (Aerts et al.2014; Scussolini et al.2017).The building-code strategies that were used in this study include wet-floodproofing, dry-floodproofing, and building elevation.The costs include both initial investments and maintenance costs (Du et al.2020; Shan et al.2021).The benefits are expressed as the reduction in EAD from the current level.For each strategy, the net present value, the payoffperiod, and the benefit/cost ratio are estimated.Net Present Values: Net present value provides the amount of net economic benefits that an adaptation strategy generates.The discount rate is a key indicator of cost-benefit evaluation, with great uncertainty and time dynamics (van den Bergh and Botzen 2014).This study used three discount rates: the Chinese official discount rate of 8% as a higher bound, a lowerbound discount rate of 4% to reflect the long-term benefits of investments, and an in-between rate of 6%, all of which have been applied in relevant studies (Du et al.2020; Shan et al.2021).It can be calculated using Eq.( 3) (Aerts et al.2014):

    where NPV stands for net present values;B tis the benefit (that is, avoided EAD) of a flood risk adaptation measure in yeart;C tis the cost in yeart;ris the discount rate; and the investment horizon isTyears.

    Benefit/Cost Ratio: The benefit/cost ratio is an evaluation index equivalent to the net present value.Comparing the benefit/cost ratios of different adaptation measures can determine the most efficient flood adaptation measures (Aerts et al.2014).The benefit/cost ratio (BCR) can be calculated using Eq.( 4):

    2.1.3 Real Options Analysis

    Real options analysis (ROA) is a probabilistic decisionmaking process.When irreversibility and uncertainty are the key characteristics of decision-making problems, ROA can assess the flexibility and adaptability of future decisions.Real option analysis evaluates the value of two options: (1) the continuation value, which is an option value when the option is deferred; and (2) the termination value, which is an option value when the option is implemented (Kim et al.2019).The real option value in any given yeartcan be defined by Eq.( 5):

    whereR tis an option value in the yeart, which is the higher of the two values;Rcon,tis a continuation value in yeart; andRter,tis a termination value in yeart.

    The termination value (Rter,t) at any yeartis the net present value of the investment made in yeart.It can be calculated using Eq.( 6):

    where EABiis the expected annual benefit of the project at yeari;ris a discount rate;Iis investment cost; andLis the service life of adaptation measures.

    A continuation value (Rcon,t+1) and a termination value (Rter,t+1) are values that the option holder can expect in yeart+1; the continuation value at yeartis the higher one of these two values discounted from yeart+1to yeart.It can be calculated using Eq.( 7):

    The calculation of the continuation value and the termination value starts from the end of the waiting period (T), using a backward induction method (Kim et al.2019).The backward induction method requires a boundary value at the end of the time horizon for option evaluation.This study set the end year at 2100.

    2.1.4 Dynamic Adaptive Policy Pathways

    Tipping point analysis (TPA) is the core of the dynamic pathways of adaptation measures, and its purpose is to measure when the adaptation plan can no longer meet the specific objectives and requirements.The dynamic adaptive policy pathways (DAPP) method provides a framework for a series of adaptation schemes that can be implemented.After the tipping point analysis is completed, it is displayed in the form of a decision tree or subway route map.Any connectable pathways in the framework is a DAPP (Fig.3) (Haasnoot et al.2020).The adaptation pathways scoreboard includes the implementation difficulty and cost ranking of each adaptation pathway.Generally speaking, this method uses scenario analysis to evaluate the distribution of the effective period among the various adaptation measures.

    2.2 Study Area

    Shanghai had a permanent population of 24.88 million in 2020, making it one of the most populous estuary megacities in the world.Shanghai is located on the coast of the Northwest Pacific, with the East China Sea to the east and Hangzhou Bay to the south.The terrain is low and flat, and the city is extremely sensitive to sea level rise and typhoon storm surges.One of the deadliest storm surge events that killed more than 29,000 people in Shanghai occurred in 1905; Typhoon Winnie killed seven people and flooded more than 5000 households in 1997.In August 2005, Typhoon Maisha hit the urban area of Shanghai and thousands of houses were flooded, followed by Typhoon Fitow in 2013, Typhoon Lekima in 2019, and Typhoon In-fa in 2021.These events exhibited large uncertainty in timing, magnitudes, and damages to Shanghai.The rainstorms brought by the typhoons may cause a rapid increase in river network flow, and land subsidence and insufficient drainage capacity in Shanghai have also exacerbated the rainstorm and waterlogging disasters in low-lying areas.In the future, climate change, sea level rise, land subsidence, and rapid urbanization will increase the risks of extreme storm floods in Shanghai (Wang et al.2018).The study area consists of 16 districts, among which Putuo (PT), Hongkou (HK), Yangpu (YP), Changning (CN), Xuhui (XH), Jingan (JA), and Huangpu (HP) Districts are located inside the Outer Ring, constituting the city center of Shanghai (Fig.4).

    2.3 Data

    The data used in this study mainly include extreme storm floods in Shanghai, land use data, Shanghai’s building data, and building cost data.

    2.3.1 Extreme Storm Floods

    A simplified 2D flood inundation model (FloodMap Inertial) was used to simulate dike failure-induced flooding and to derive inundation maps.For the risk decision of sea level rise, the upper limit scenario is usually considered (Hinkel et al.2015).The representative concentration pathways RCP8.5 and RCP2.6 (the lower limit scenario of sea level rise) were used to estimate the multiple climate change scenarios and probability of future sea level rise.All relative sea level projections by 2100 are given with the baseline year of 2010 (Yin et al.2020).Then, according to the height and structure of the river embankments along the Shanghai coast, the overwhelming and breaching discharges of storm surges in 1/10, 100, and 1000-years are calculated.A 2-day (48 h) tidal cycle including four rising phases and four falling limbs is applied in the flood simulations.The two-dimensional hydrodynamic model (FloodMap-Inertial) was used to simulate the baseline year (2010) and future (2030, 2050, and 2100) flood inundation depths of extreme storm floods under the RCP8.5 and RCP2.6 scenarios, with a spatial resolution of 50 m by 50 m.Figure 5 shows the most extreme storm flood inundation map in Shanghai under the RCP8.5 scenario.The simulation results have been verified and used in previous research (Yin et al.2020; Shan et al.2022).

    2.3.2 Land Use Type

    The land use data of Shanghai in 2010 were provided by the Shanghai Institute of Surveying and Mapping.Shanghai’s land use data in 2030, 2050, and 2100 are based on the 2010 land use data and were simulated using the FLUS model,with a spatial resolution of 50 m by 50 m.The Kappa coefficient and figure of merit index were used to quantitatively assess the simulation results, which show that the simulated land use has high accuracy and meets the research needs (Shan et al.2022).

    2.3.3 Building Data

    The building data of Shanghai in 2017 included the buildings’ footprint and height information, which were obtained from Baidu Map1https://map.baidu.com.using a python-based web crawler and then processed by a GIS software.The resulting total buildings’ area is 1238.2 km2http://tjj.sh.gov.cn/tjnj/index.html..Shan et al.( 2019) selected one of the Shanghai building types—residential buildings—and used Google Earth high-precision remote sensing images to evaluate the accuracy of residential building data, which proved that the accuracy is high.Based on the spatial distribution of buildings in 2017 and according to the residential, commercial, public, and industrial land areas in 2030, 2050, and 2100, the area of residential, commercial, public, and industrial buildings in Shanghai in 2030, 2050, and 2100 can be obtained according to a certain proportion coefficient.The results show that the total buildings’ area of Shanghai in 2030, 2050, and 2100 is 1518.1 km 2 , 1739.4 km 2 , and 1879.4 km 2 , respectively.

    2.3.4 Building Cost Data

    The completed construction floor area and cost of construction by building use type were obtained from the Shanghai Statistical Yearbooks (1995—2019).2 The construction cost per square meter was calculated by dividing the cost of construction of buildings completed by the building floor area completed, which was applied in the study by Wu et al.( 2019).For the prediction of the construction cost of residential buildings, commercial buildings, and public buildings, a polynomial curve was fitted on the future construction cost trend, by assuming that this trend will continue until 2050.Industrial building construction cost was based on the construction cost report of China’s mainland and Hong Kong released by the construction engineering consulting company ARCADIS,33 https://www.arcadis.com/en.which published the construction cost of different building types in Shanghai.According to a certain average annual growth rate, the cost of industrial construction in 2050 is predicted.Due to the uncertainty of future building costs on a long-term scale, this study assumed that the building cost in 2100 will remain the same as in 2050.

    3 Results

    This section presents the results of this study, which include dynamic changes of flood risk; performance of adaptation measures; and robust decision making of adaptation strategies under deep uncertainties.

    3.1 Dynamic Change Analysis of Flood Risks

    Under the RCP2.6 and RCP8.5 scenarios, the exposed area of the four types of buildings in Shanghai in 2100 will be 905.9 km 2 and 1010.9 km 2 respectively—that is 4.6 times and 6.4 times the 1/10-year flood, respectively—and the exposed assets will be 5—9 times larger in value.Of the 16 districts, the exposed building assets are mainly concentrated in the waterfront districts of Pudong and Baoshan and the city center under the two emission scenarios.Pudong has the largest proportion of exposed assets, accounting for about 30—33% of the total exposed assets.For the 1/1,000-year floods, the exposed assets of Shanghai city center amount to 33.3% in 2100.

    The spatial pattern of building losses changed from a scattered distribution at the 10-year return period to a contiguous distribution at the 1000-year return period, with a total of USD 15.93 billion for the 1000-year floods under the RCP8.5 scenario, which is 10.6—17.9 times the 10-year flood losses (Fig.6).Of the 16 districts, the loss of buildings is mainly concentrated in Pudong and Baoshan Districts and the city center.For different building types, residential building losses in 2100 are the largest at USD 6.52 billion, followed by commercial buildings with USD 3.81 billion in losses and public buildings with USD 3.66 billion in losses; industrial buildings have the smallest losses at USD 1.94 billion.

    The annual exceedance probability (AEP)-loss curves of extreme flooding for Shanghai’s buildings under the RCP8.5 and RCP2.6 scenarios are based on the losses under three flood return periods.The expected annual damage (EAD) was calculated by Eq.( 2).In the current scenario, the EAD of Shanghai’s buildings in 2010 was USD 22.35 million/year.Under the RCP2.6 scenario, the EAD of the buildings is 2—10 times that of the current scenario.Under the RCP8.5 scenario, the EAD of the buildings in 2030, 2050, and 2100is USD 40.01 million/year, USD 79.71 million/year, and USD 259.30 million/year, respectively—that is 2—15 times that of the current scenario.Pudong, Baoshan, and Xuhui Districts have the highest EAD.The EAD of the city center in 2010, 2030, 2050, and 2100 accounts for 27.7%, 29.6%, 32.7%, and 40% of the total, respectively.

    3.2 Performance of Adaptation Measures

    Using Eqs.( 3) and ( 4), the net present value and benefit/cost ratio of dry-floodproofing, wet-floodproofing, and elevation are obtained under different discount rates (4%, 6%, 8%) (Fig.7).Under the RCP2.6 scenario, the implementation of the three flood adaptation measures will effectively reduce the risks of Shanghai’s buildings in 2030, 2050, and 2100.If we implement wet-floodproofing in Shanghai in 2030, the net present value is USD 0.59—0.88 billion and the benefit/cost ratio is 6.87—7.54.The net present value and benefit/cost ratio in 2100 are 10.8—14.0 times and 0.8—1.1 times of 2030, respectively.Combining the dry and wet-floodproofing measures to form a “combined” strategy leads to a net present value of USD 0.41—0.64 billion and a benefit/cost ratio of 2.23—2.46 in 2030.It can be concluded that the “combined” strategy cannot generate good benefits in the near future (2030); however, the net present value and benefit/cost ratio in 2100 are 18.6—29.5 times and 3.3—4.3 times the 2030 values, respectively.We also formulated a “hybrid” strategy by combining dry-floodproofing, wet-floodproofing, and elevation, which has a better performance in terms of low residual risk and high benefit/cost ratio.The net present value and benefit/cost ratio in 2100 are USD 6.95—17.81 billion and 4.37—6.36, respectively.Under the RCP8.5 scenario, the net present value and benefit/cost ratio of wetfloodproofing implemented in 2100 will be 9.0—15.0 times and 1.3—1.4 times the 2030 values, respectively.Comparing the net present value and benefit/cost ratio of the three flood adaptation measures, wet-floodproofing measures have the best flood control outcome, followed by dry-floodproofing, and elevation.The net present values of the combined strategy (dry-and wet-floodproofings) and the hybrid strategy (dry-and wet-floodproofings and elevation) in 2100 are USD 7.21—14.20 billion and USD 6.44—13.29 billion, respectively.Correspondingly, the benefit/cost ratios are 4.79—6.19 and 3.08—3.99, respectively.The net present value and benefit/cost ratio of the hybrid strategy are significantly higher than those of the individual measures as well as the combined strategy.However, its cost is higher in the shorter term.Therefore, the hybrid strategy is suitable for long-term planning.

    Under the RCP2.6 and RCP8.5 scenarios, through the implementation of flood adaptation measures, the changes in Shanghai’s building EAD before and after the implementation of the measures were evaluated and compared (Fig.8).Under the RCP8.5 scenario, after the implementation of dry-floodproofing, the annual benefits (that is, the reduction in EAD) of the four types of buildings can reach USD 37.04 million to USD 169.74 million.In addition, after the implementation of wet-floodproofing, the annual benefit is USD 38.61 million to USD 202.88 million; after the implementation of elevation, the annual benefit is USD 26.73—136.45 million.The annual benefits of implementing dry-floodproofing, wet-floodproofing, and elevation measures in 2030 is USD 37.04 million, USD 38.61 million, and USD 26.73 million, respectively.After the implementation of all three measures, the annual benefits of the four types of buildings in 2100 is USD 169.74 million, USD 202.88 million, and USD 136.45 million, respectively.These results show that the adaptations tend to generate better benefits in the long term, indicating that other coastal cities should begin their adaptation plans as early as possible.

    3.3 Robust Decision Making of Adaptation Strategies Under Deep Uncertainties

    Table 1 Results of option values and optimal investment times for different sea level rise (SLR) scenarios in Shanghai

    We used Eqs.( 5)—( 7), the continuation value, and the termination value of the three flood adaptation measures to determine the best time for investment.The optimal investment time was determined by checking whether expected annual benefit (EAB) reaches the critical threshold value (Table 1).We found that the critical threshold of the optimal investment time highly depends on the investment cost and discount rate.The more extreme the sea level rise scenario is (that is, RCP8.5), the earlier the optimal investment time appears.Therefore, a rapid sea level rise requires an immediate investment to obtain the maximum benefit value.

    Under the RCP2.6 and RCP8.5 scenarios, the dry-floodproofing measure can be used in the initial period of 2033 and 2036.With time lapses, the performance of the dryfloodproofing may drop, and wet-floodproofing implemented in 2035 and 2042 will obtain better economic benefits.The elevation measures also show good performance at the beginning of implementation, but the performance will decline after 2041—2045.Dry-floodproofing and wetfloodproofing measures can be transformed into the combined strategy in 2044—2046.The combined strategy can be upgraded to the hybrid strategy in 2047—2051.

    Figure 9 shows the optimal investment time for the three adaptation measures.Under the RCP2.6 and RCP8.5 scenarios, the dry-floodproofing measure is implemented first.Under the RCP2.6 scenario, after the wet-floodproofing measure is implemented for 6 years the elevation measure is implemented, while under the RCP8.5 scenario, after the implementation of the wet-floodproofing measure for 3 years the elevation measure is implemented immediately.

    Meanwhile, in order to enable decision makers to make robust decisions, three considerations were integrated: (1) economic benefits of measures; (2) the continuation value and the termination value; and (3) tipping points.The dynamic pathways for Shanghai’s future buildings to adapt to flood risks were formulated for both the short and long terms and are presented in Fig.10.The figure shows a variety of adaptation pathways, and the tipping point of adaptation can be evaluated from the effectiveness of adaptation measures.A tipping point is the time period during which a given measure no longer meets its goal.These tipping points can be expressed in terms of the progression of the analyzed climate variable (that is, sea level rise) on a temporal scale.Through this map, the decision maker will know when to implement the investment in order to generate the largest net benefit (that is, the overall benefits minus the overall costs).It is necessary to point out that the timely switch between adaptation measures critically depends on periodic monitoring of the sea level rise and revising the estimations of its impacts and the effectiveness of existing measures.

    4 Discussion

    Future uncertainties in climate change and sea level rise pose a major challenge to flood risk management and adaptation in coastal megacities.This study applied an integrated modeling framework to assess the flood risk, provide a cost-benefit analysis of adaptation measures used to reduce risks, and coupled a ROA analysis with the concept of DAPP.Taking Shanghai as a case study, we developed a multidisciplinary approach to understand impacts and adaptation of coastal cities to sea level rise.We have found that sea level rise and land subsidence will significantly exacerbate the flood risks in Shanghai.Under RCP2.6 and RCP8.5, we identified that buildings with high risks are located in the city center and along the banks of the Huangpu River.We also found that with the rise of the sea level caused by climate warming, under the RCP8.5 scenario, the loss of Shanghai’s buildings in the future (2030, 2050, and 2100) accounts for 4.6%, 5.2% and 8.8% of the exposed assets, respectively.Climate change and sea level rise further increase the risk of flood, which are also a signal for other coastal megacities to plan for early adaptation in highly exposed and economically important areas.

    The formulation of climate change adaptation measures needs to fully consider the uncertainty of future climate change in order to find a reliable solution for planning and design of flood control infrastructures.According to Yin et al.( 2020), if no adaptation measures are taken, the risk faced by Shanghai is expected to increase by 3—160 times by 2100, which may repeat the tragedy of New Orleans.The results of this study show that the EAD of Shanghai’s buildings will increase by approximately 15 times in 2100 under the RCP8.5 scenario.Under the current, RCP2.6, and RCP8.5 scenarios, the three flood adaptation measures advocated in this study—through the use of building codes to reduce the exposure and vulnerability of the elements at risk—can effectively reduce the risks faced by Shanghai’s buildings in the future.Compared to the adaptation effects of dry-or wet-floodproofing, the effect of the elevation measure is smaller.Wang et al.( 2015) showed that raising the elevation of old buildings may not bring substantive benefits at all, because the moderate long-term benefits may be offset by a large amount of social costs and high building renovation costs.Lifting new buildings separately is more costeffective, easier to accept by society, and an easier policy to implement.Therefore, in the future risk adaptation in Shanghai, if the elevation measure is adopted, it should consider new buildings, where the maximum benefit will be obtained.In addition to adopting the three measures proposed in this study, for newly built low-lying communities, floor elevation of the lowest floor must exceed the local flood warning line.

    From the perspective of decision support, the future challenge is not only the improvement of technical means, but also the application of risk assessment models to practical decision making, which is reflected in how to deal with uncertainty in the decision making process and quantify adaptation measures, and how to promote communication with decision makers and stakeholders.The combination of ROA and DAPP methods can help decision makers distinguish the advantages and disadvantages of the implementation results of their plans, and formulate the optimal combination of adaptation strategies and flexible adaptation pathways.The ROA and DAPP methods have been increasingly used to design or plan urban development for decades (Woodward et al.2014; Buurman and Babovic 2016).Although Shanghai has established a system of policies and measures to adapt to climate change, there is a lack of a dynamic adaptation pathways map of measures for decision makers and stakeholders.The comprehensive analysis of ROA and DAPP can guide decision makers to make robust decisions, and provide a reference for flood risk adaptation in other coastal cities and delta areas.

    There are some limitations to our study.First, the vulnerability curve established by Ke ( 2014) was used to evaluate the loss of buildings in Shanghai in 2030, 2050, and 2100.With the acceleration of urbanization and the renewal of the city, the flood vulnerability of buildings may change dynamically and need an updated vulnerability curve.Second, due to the current rapid urbanization process and many uncertain factors, there are inevitably limitations in the prediction of the spatial distribution of buildings and costs in the long term, especially for 2050 and 2100.In future research, Internet big data and machine learning methods can be adopted to map the building footprint, height, cost, and economic values, providing higher resolution spatial and economic information of elements at risk.

    5 Conclusion

    This study integrated risk analysis, cost-benefit analysis, real options analysis (ROA), and dynamic adaptive policy pathways (DAPP) methods to analyze the flood risk, adaptation measures, and adaptation decision making of coastal cities.With the case study in Shanghai, we quantified the future risks of buildings under extreme storm flooding scenarios for 2030, 2050, and 2100, analyzed the cost and benefit of three flood adaptation measures to reduce flood risk, and formulated the adaptation dynamic pathways for decision making.The main conclusions are:

    (1) Sea level rise, land subsidence, and socioeconomic development will cause a sharp increase in the potential future flood risk.Under the RCP8.5 scenario, the EAD of the buildings in 2030, 2050, and 2100 will be 2—15 times that of the current scenario.

    (2) Three adaptation measures based on building codes can effectively reduce the risks of Shanghai’s buildings in 2030, 2050, and 2100.Wet-floodproofing has the best flood control effect, followed by dry-floodproofing and building elevation.

    (3) In terms of the adaptation pathways, we found that the flexible transformation of adaptation measures can achieve the highest economic benefits.Dry-floodproofing can be used in the initial period of 2033 and 2036, and should be replaced by wet-floodproofing in 2035 and 2042 under the RCP2.6 and RCP8.5 scenarios.The elevation measures also showed good performance at the beginning of implementation, but the gain in benefit will decline after 2041—2045.Dry-floodproofing and wet-floodproofing measures can be used together to form a combined strategy in 2044—2046.A hybrid strategy of dry-and wet-floodproofing and elevation measures can be adopted in 2047—2051.

    However, future changes in storminess were not considered in the extreme storm flood simulation, and the cost effectiveness of adaptation strategies should be higher if their benefit of reducing indirect flood risk were considered.Moreover, the performance of adaptation strategies should be explained in terms of residual flood risks.Future studies can be improved through incorporating these aspects.

    AcknowledgementsThis research was funded by the National Key Research and Development Program of China (Grant No.2018YFC1508803), the National Social Science Foundation of China (Grant No.18ZDA105), the National Natural Science Foundation of China (Grant No.41971199, 42171080, 42001182), and the Shanghai Science and Technology Support Program (Grant No.19DZ1201505).

    Open AccessThis article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made.The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material.If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder.To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

    超碰97精品在线观看| 在线观看免费视频日本深夜| 国产高清国产精品国产三级| 免费在线观看黄色视频的| 一区二区三区精品91| 电影成人av| 男女免费视频国产| 久久中文字幕一级| 国产一区二区激情短视频| 国产又爽黄色视频| 在线观看免费视频日本深夜| 怎么达到女性高潮| 久久国产精品人妻蜜桃| 最新在线观看一区二区三区| 精品熟女少妇八av免费久了| 黄色片一级片一级黄色片| 三级毛片av免费| 欧美日韩成人在线一区二区| 国内毛片毛片毛片毛片毛片| 日本黄色日本黄色录像| 欧美激情 高清一区二区三区| 久久天堂一区二区三区四区| 村上凉子中文字幕在线| 天天躁夜夜躁狠狠躁躁| 成人影院久久| 亚洲精华国产精华精| 亚洲成人国产一区在线观看| 日本精品一区二区三区蜜桃| 99国产精品免费福利视频| av视频免费观看在线观看| 日韩精品免费视频一区二区三区| 国产精华一区二区三区| 交换朋友夫妻互换小说| 亚洲黑人精品在线| 男女之事视频高清在线观看| 亚洲aⅴ乱码一区二区在线播放 | 真人做人爱边吃奶动态| 怎么达到女性高潮| 真人做人爱边吃奶动态| 免费av中文字幕在线| 日韩有码中文字幕| 免费在线观看亚洲国产| 18禁观看日本| 少妇猛男粗大的猛烈进出视频| 91老司机精品| 亚洲精品国产一区二区精华液| 一区福利在线观看| 国产亚洲精品第一综合不卡| 精品亚洲成a人片在线观看| videos熟女内射| 12—13女人毛片做爰片一| 人妻丰满熟妇av一区二区三区 | 欧美色视频一区免费| 十八禁人妻一区二区| 亚洲精品中文字幕一二三四区| 大型av网站在线播放| 丝袜美腿诱惑在线| 亚洲成人国产一区在线观看| 国产精品久久久久久精品古装| 亚洲午夜精品一区,二区,三区| 国产精品1区2区在线观看. | 美女扒开内裤让男人捅视频| 99香蕉大伊视频| 国产精品久久久久久人妻精品电影| 视频区欧美日本亚洲| 欧美老熟妇乱子伦牲交| 午夜福利乱码中文字幕| 亚洲 国产 在线| av电影中文网址| 夜夜躁狠狠躁天天躁| 50天的宝宝边吃奶边哭怎么回事| 欧美人与性动交α欧美精品济南到| 亚洲国产精品一区二区三区在线| 日韩三级视频一区二区三区| 热99re8久久精品国产| 午夜福利免费观看在线| 国产又色又爽无遮挡免费看| 久久婷婷成人综合色麻豆| av不卡在线播放| 一级毛片精品| 久久午夜综合久久蜜桃| 狂野欧美激情性xxxx| 成人亚洲精品一区在线观看| 精品国产美女av久久久久小说| 国产一卡二卡三卡精品| 国产一区二区三区在线臀色熟女 | 大香蕉久久网| 国产亚洲精品第一综合不卡| 亚洲一区中文字幕在线| 久久国产亚洲av麻豆专区| 女人爽到高潮嗷嗷叫在线视频| 在线十欧美十亚洲十日本专区| av不卡在线播放| 啦啦啦视频在线资源免费观看| 国产在视频线精品| 亚洲精品一卡2卡三卡4卡5卡| 亚洲熟女毛片儿| 亚洲 国产 在线| 久久久久久免费高清国产稀缺| 看黄色毛片网站| 一边摸一边抽搐一进一小说 | 18禁黄网站禁片午夜丰满| 国产成人av教育| www.精华液| 亚洲色图综合在线观看| 一区二区三区国产精品乱码| 欧美中文综合在线视频| 成熟少妇高潮喷水视频| 欧美一级毛片孕妇| 久久久国产一区二区| 欧美久久黑人一区二区| 9191精品国产免费久久| 国产精品免费视频内射| 91字幕亚洲| 国产无遮挡羞羞视频在线观看| 亚洲全国av大片| 国产淫语在线视频| 日本精品一区二区三区蜜桃| 国产主播在线观看一区二区| 91精品三级在线观看| 男女之事视频高清在线观看| 国产精品成人在线| 亚洲成a人片在线一区二区| 久久影院123| 久久精品国产综合久久久| 久久午夜综合久久蜜桃| 99国产精品一区二区蜜桃av | 亚洲色图综合在线观看| 变态另类成人亚洲欧美熟女 | 日韩三级视频一区二区三区| 狠狠狠狠99中文字幕| 中文字幕av电影在线播放| 丁香六月欧美| 国产精品永久免费网站| 制服人妻中文乱码| 欧美乱色亚洲激情| 每晚都被弄得嗷嗷叫到高潮| av免费在线观看网站| 欧美精品啪啪一区二区三区| 亚洲专区字幕在线| 亚洲国产欧美网| 男男h啪啪无遮挡| 亚洲色图综合在线观看| 搡老乐熟女国产| 91麻豆av在线| 亚洲中文日韩欧美视频| 亚洲色图av天堂| 久久久久久人人人人人| 国产精品98久久久久久宅男小说| 一级片'在线观看视频| 在线观看免费午夜福利视频| 国产亚洲欧美精品永久| netflix在线观看网站| 老司机在亚洲福利影院| 男人舔女人的私密视频| 宅男免费午夜| 99国产综合亚洲精品| 国产一区二区三区综合在线观看| 亚洲av日韩在线播放| 99香蕉大伊视频| 精品久久久久久久久久免费视频 | 国产97色在线日韩免费| av一本久久久久| 极品人妻少妇av视频| 母亲3免费完整高清在线观看| 久久国产精品大桥未久av| 老鸭窝网址在线观看| 99精品久久久久人妻精品| 母亲3免费完整高清在线观看| 亚洲,欧美精品.| 日韩欧美三级三区| 久久久水蜜桃国产精品网| 精品国产乱码久久久久久男人| 久久午夜亚洲精品久久| 国产精品 欧美亚洲| 无限看片的www在线观看| 欧美 日韩 精品 国产| 久久精品国产清高在天天线| 大型av网站在线播放| 国产成人av教育| 日本撒尿小便嘘嘘汇集6| 国产在线精品亚洲第一网站| 一边摸一边抽搐一进一出视频| 国产无遮挡羞羞视频在线观看| 伊人久久大香线蕉亚洲五| 91麻豆精品激情在线观看国产 | 岛国毛片在线播放| 交换朋友夫妻互换小说| 女人精品久久久久毛片| 丰满迷人的少妇在线观看| 精品欧美一区二区三区在线| 操出白浆在线播放| 精品欧美一区二区三区在线| 窝窝影院91人妻| 不卡av一区二区三区| 又黄又爽又免费观看的视频| 丝袜在线中文字幕| 中文字幕制服av| 老汉色∧v一级毛片| 精品电影一区二区在线| 国产一区有黄有色的免费视频| 日韩欧美一区二区三区在线观看 | 最近最新免费中文字幕在线| 久9热在线精品视频| 18禁裸乳无遮挡动漫免费视频| 午夜久久久在线观看| av国产精品久久久久影院| 热re99久久精品国产66热6| 1024香蕉在线观看| 欧美日韩中文字幕国产精品一区二区三区 | 天天躁夜夜躁狠狠躁躁| 午夜福利视频在线观看免费| 国产不卡一卡二| 91精品三级在线观看| 日韩三级视频一区二区三区| 国产成人精品久久二区二区91| 日韩一卡2卡3卡4卡2021年| 久久 成人 亚洲| 不卡一级毛片| 少妇被粗大的猛进出69影院| 在线永久观看黄色视频| 一本一本久久a久久精品综合妖精| 捣出白浆h1v1| 久久精品亚洲av国产电影网| 成人18禁高潮啪啪吃奶动态图| av免费在线观看网站| av欧美777| 亚洲av第一区精品v没综合| 大片电影免费在线观看免费| 一本综合久久免费| 老熟妇乱子伦视频在线观看| 亚洲成人免费电影在线观看| 一区二区日韩欧美中文字幕| 精品免费久久久久久久清纯 | 妹子高潮喷水视频| 精品国产美女av久久久久小说| 色94色欧美一区二区| 国产免费现黄频在线看| 亚洲av成人av| 色94色欧美一区二区| 中文欧美无线码| av天堂久久9| 侵犯人妻中文字幕一二三四区| 久久国产乱子伦精品免费另类| 国产精品久久久人人做人人爽| 视频区欧美日本亚洲| 精品国产乱码久久久久久男人| 色精品久久人妻99蜜桃| 女人被狂操c到高潮| 亚洲av成人一区二区三| 一进一出抽搐gif免费好疼 | 午夜福利,免费看| 久久国产精品大桥未久av| 黄片大片在线免费观看| 成熟少妇高潮喷水视频| av电影中文网址| 女人被躁到高潮嗷嗷叫费观| av视频免费观看在线观看| 国产亚洲精品久久久久5区| 久久草成人影院| netflix在线观看网站| 亚洲少妇的诱惑av| 女性生殖器流出的白浆| 天天操日日干夜夜撸| 中亚洲国语对白在线视频| 色婷婷av一区二区三区视频| 韩国av一区二区三区四区| tube8黄色片| 久久精品国产亚洲av高清一级| 免费av中文字幕在线| 久久香蕉精品热| 国产精品一区二区在线观看99| 高清毛片免费观看视频网站 | 一a级毛片在线观看| 亚洲中文av在线| 一区福利在线观看| 美女高潮到喷水免费观看| 高清黄色对白视频在线免费看| 十分钟在线观看高清视频www| 欧美乱码精品一区二区三区| 一区二区三区激情视频| 波多野结衣av一区二区av| 国产欧美日韩精品亚洲av| 亚洲av熟女| 精品无人区乱码1区二区| 久久国产精品人妻蜜桃| 欧美国产精品va在线观看不卡| 成人av一区二区三区在线看| 国产一区在线观看成人免费| 欧美日韩中文字幕国产精品一区二区三区 | 大香蕉久久网| 亚洲精品自拍成人| 99热网站在线观看| 怎么达到女性高潮| 亚洲专区国产一区二区| 午夜福利视频在线观看免费| 十八禁网站免费在线| av超薄肉色丝袜交足视频| 国产精品亚洲av一区麻豆| 成人精品一区二区免费| 国产精品国产高清国产av | 性少妇av在线| 国产精品久久久久成人av| 纯流量卡能插随身wifi吗| 19禁男女啪啪无遮挡网站| 无人区码免费观看不卡| 男男h啪啪无遮挡| 亚洲少妇的诱惑av| 中文字幕人妻丝袜一区二区| 母亲3免费完整高清在线观看| 狂野欧美激情性xxxx| 美女福利国产在线| 新久久久久国产一级毛片| 男女之事视频高清在线观看| 免费久久久久久久精品成人欧美视频| 侵犯人妻中文字幕一二三四区| 精品国产美女av久久久久小说| 夜夜爽天天搞| 制服人妻中文乱码| 美女视频免费永久观看网站| 极品教师在线免费播放| av网站在线播放免费| 欧美成人午夜精品| 岛国毛片在线播放| 99久久精品国产亚洲精品| 黄色毛片三级朝国网站| 亚洲全国av大片| 国产不卡av网站在线观看| 国产成人精品无人区| 无人区码免费观看不卡| 在线播放国产精品三级| 婷婷成人精品国产| 黄色视频,在线免费观看| 校园春色视频在线观看| 丰满饥渴人妻一区二区三| 一个人免费在线观看的高清视频| 777米奇影视久久| 亚洲色图av天堂| 国产野战对白在线观看| 俄罗斯特黄特色一大片| 久久久久久亚洲精品国产蜜桃av| 50天的宝宝边吃奶边哭怎么回事| 丝袜在线中文字幕| 人人妻人人澡人人爽人人夜夜| 久久久久久亚洲精品国产蜜桃av| a级毛片黄视频| 宅男免费午夜| 黄片大片在线免费观看| 国产激情久久老熟女| 一级a爱视频在线免费观看| 91精品三级在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| 国产成人av激情在线播放| 热re99久久精品国产66热6| 19禁男女啪啪无遮挡网站| 成熟少妇高潮喷水视频| 久久国产精品影院| 一二三四社区在线视频社区8| 久久精品国产99精品国产亚洲性色 | 黄色丝袜av网址大全| 欧美乱色亚洲激情| 精品国产一区二区三区四区第35| 男女下面插进去视频免费观看| 精品久久久久久久久久免费视频 | 啦啦啦在线免费观看视频4| 国产精品免费视频内射| 免费看a级黄色片| 免费观看人在逋| 老司机福利观看| 九色亚洲精品在线播放| 午夜福利一区二区在线看| 99riav亚洲国产免费| 国产aⅴ精品一区二区三区波| 久久久精品区二区三区| 日韩欧美一区二区三区在线观看 | 黑人猛操日本美女一级片| 亚洲精品国产区一区二| 他把我摸到了高潮在线观看| a级毛片在线看网站| 男女高潮啪啪啪动态图| 亚洲精品成人av观看孕妇| 国产欧美日韩一区二区三| 亚洲精华国产精华精| 亚洲欧洲精品一区二区精品久久久| 侵犯人妻中文字幕一二三四区| 久久天堂一区二区三区四区| 国产精品久久视频播放| 搡老乐熟女国产| 999久久久精品免费观看国产| 十八禁网站免费在线| 午夜精品久久久久久毛片777| 日日摸夜夜添夜夜添小说| 成人国产一区最新在线观看| 91成年电影在线观看| 免费在线观看亚洲国产| 国产成+人综合+亚洲专区| 啦啦啦免费观看视频1| 精品福利永久在线观看| 欧美日韩亚洲综合一区二区三区_| 久99久视频精品免费| 99久久99久久久精品蜜桃| 久久久国产欧美日韩av| 亚洲七黄色美女视频| 久久久国产一区二区| 变态另类成人亚洲欧美熟女 | 成熟少妇高潮喷水视频| 老司机深夜福利视频在线观看| 成年动漫av网址| 国产av精品麻豆| 99久久综合精品五月天人人| 老熟妇仑乱视频hdxx| 精品福利永久在线观看| 少妇猛男粗大的猛烈进出视频| 少妇被粗大的猛进出69影院| 女警被强在线播放| 又紧又爽又黄一区二区| 99精国产麻豆久久婷婷| 久久久国产一区二区| 女同久久另类99精品国产91| 免费在线观看完整版高清| 操出白浆在线播放| 欧美乱色亚洲激情| 天堂动漫精品| 美女扒开内裤让男人捅视频| 午夜视频精品福利| 天堂中文最新版在线下载| 国产激情欧美一区二区| 一级a爱片免费观看的视频| 久久精品成人免费网站| 我的亚洲天堂| 一级a爱片免费观看的视频| 麻豆国产av国片精品| 不卡av一区二区三区| 日韩三级视频一区二区三区| 成人三级做爰电影| av电影中文网址| 国产精品国产av在线观看| 国产97色在线日韩免费| 老司机亚洲免费影院| 久久人妻熟女aⅴ| 人人妻人人澡人人看| 欧美国产精品va在线观看不卡| 黄片小视频在线播放| 热99re8久久精品国产| 黑人巨大精品欧美一区二区mp4| 国产精品自产拍在线观看55亚洲 | 精品亚洲成国产av| 国产精品成人在线| 欧美精品亚洲一区二区| 亚洲欧美一区二区三区黑人| 一级毛片女人18水好多| 国产真人三级小视频在线观看| 亚洲熟女精品中文字幕| 亚洲国产精品合色在线| 久久精品人人爽人人爽视色| 国产成人啪精品午夜网站| 国产高清视频在线播放一区| 满18在线观看网站| 中文字幕高清在线视频| 精品乱码久久久久久99久播| 婷婷丁香在线五月| 日本五十路高清| 欧美色视频一区免费| 久热爱精品视频在线9| 黄色视频,在线免费观看| 精品第一国产精品| 国产男靠女视频免费网站| 最新的欧美精品一区二区| 咕卡用的链子| 成年人免费黄色播放视频| 成人影院久久| 国产主播在线观看一区二区| 少妇的丰满在线观看| 三上悠亚av全集在线观看| 最新的欧美精品一区二区| 亚洲精品一卡2卡三卡4卡5卡| 欧美日韩乱码在线| 国产精品电影一区二区三区 | 精品国产一区二区久久| 国产视频一区二区在线看| 国产精品av久久久久免费| 久久久久久久久免费视频了| 巨乳人妻的诱惑在线观看| 波多野结衣一区麻豆| 下体分泌物呈黄色| 国产亚洲欧美在线一区二区| 色尼玛亚洲综合影院| 国产男靠女视频免费网站| 国产精华一区二区三区| 亚洲va日本ⅴa欧美va伊人久久| 午夜福利免费观看在线| 国产精品电影一区二区三区 | 免费观看精品视频网站| 午夜免费观看网址| 伦理电影免费视频| 在线永久观看黄色视频| 欧美国产精品一级二级三级| 老熟妇仑乱视频hdxx| 最新的欧美精品一区二区| 国产亚洲精品一区二区www | 久久久国产一区二区| 亚洲全国av大片| 99riav亚洲国产免费| 国产免费男女视频| 黄色视频不卡| 欧美精品av麻豆av| 国产一区有黄有色的免费视频| 久9热在线精品视频| 欧美日韩精品网址| 亚洲性夜色夜夜综合| 50天的宝宝边吃奶边哭怎么回事| 日本欧美视频一区| 视频区图区小说| 亚洲三区欧美一区| 国产成人精品久久二区二区91| 少妇猛男粗大的猛烈进出视频| 宅男免费午夜| 老司机福利观看| 久久精品国产a三级三级三级| 免费在线观看黄色视频的| 欧美 日韩 精品 国产| 波多野结衣一区麻豆| 丰满饥渴人妻一区二区三| 亚洲一区二区三区欧美精品| 国产精品亚洲一级av第二区| 母亲3免费完整高清在线观看| av不卡在线播放| 国产亚洲欧美在线一区二区| 国产高清国产精品国产三级| 国产在线观看jvid| 午夜老司机福利片| 国产高清视频在线播放一区| 国产av又大| 女同久久另类99精品国产91| 国产精品 国内视频| 激情在线观看视频在线高清 | 国产成人精品久久二区二区免费| 91精品三级在线观看| av天堂久久9| 日韩欧美一区二区三区在线观看 | 欧美日韩av久久| 啪啪无遮挡十八禁网站| 曰老女人黄片| 亚洲一码二码三码区别大吗| 热99国产精品久久久久久7| 天堂动漫精品| 伦理电影免费视频| 国产高清激情床上av| 成年女人毛片免费观看观看9 | 黄片小视频在线播放| 法律面前人人平等表现在哪些方面| 久久久国产欧美日韩av| 国产区一区二久久| 国产男靠女视频免费网站| 在线av久久热| 精品亚洲成国产av| 国产单亲对白刺激| 国产精品香港三级国产av潘金莲| 欧洲精品卡2卡3卡4卡5卡区| 久久久精品免费免费高清| 搡老乐熟女国产| 久久中文看片网| 伦理电影免费视频| 日韩欧美一区二区三区在线观看 | 国产亚洲精品一区二区www | 国产99久久九九免费精品| 亚洲av美国av| xxx96com| 久久久久精品人妻al黑| 日日夜夜操网爽| 亚洲av电影在线进入| 欧美日韩国产mv在线观看视频| 精品国产一区二区三区久久久樱花| 天天躁日日躁夜夜躁夜夜| 狠狠婷婷综合久久久久久88av| 欧美乱色亚洲激情| 黄频高清免费视频| 超碰成人久久| 亚洲欧美一区二区三区黑人| 老司机亚洲免费影院| 啪啪无遮挡十八禁网站| 又紧又爽又黄一区二区| 久久精品91无色码中文字幕| 精品福利观看| 亚洲av欧美aⅴ国产| 国产精品1区2区在线观看. | 国产aⅴ精品一区二区三区波| 人人妻人人爽人人添夜夜欢视频| 看免费av毛片| 黑丝袜美女国产一区| 热99国产精品久久久久久7| 久久天躁狠狠躁夜夜2o2o| 精品少妇一区二区三区视频日本电影| 91国产中文字幕| 女人被躁到高潮嗷嗷叫费观| 三级毛片av免费| 免费黄频网站在线观看国产| 可以免费在线观看a视频的电影网站| 日韩成人在线观看一区二区三区| 欧美一级毛片孕妇| 在线观看免费视频网站a站| 欧美日韩亚洲国产一区二区在线观看 | www.自偷自拍.com| 久久草成人影院| 欧美午夜高清在线| 人人妻人人添人人爽欧美一区卜| 精品人妻熟女毛片av久久网站| 精品高清国产在线一区| 最近最新中文字幕大全免费视频| 下体分泌物呈黄色| 超色免费av| 美女视频免费永久观看网站| 别揉我奶头~嗯~啊~动态视频| 国产欧美日韩一区二区三区在线| 午夜老司机福利片|