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

    Global Rainstorm Disaster Risk Monitoring Based on Satellite Remote Sensing

    2022-03-12 07:52:18SulingRENXiuzhenHANJunYANGShihaoTANGYangZHANGTianchanSHANandChengLIU
    Journal of Meteorological Research 2022年1期

    Suling REN, Xiuzhen HAN, Jun YANG, Shihao TANG, Yang ZHANG, Tianchan SHAN, and Cheng LIU

    1 National Satellite Meteorological Center, China Meteorological Administration, Beijing 100081

    2 Heavy Rain and Drought–Flood Disasters in Plateau and Basin Key Laboratory of Sichuan Province, Chengdu 610072

    ABSTRACT

    Key words: satellite-retrieved precipitation, flood, rainstorm disaster risk monitoring, rainstorm disaster risk index

    1. Introduction

    Rainstorms and brief torrential rain can cause flooding, landslides, debris flow, urban waterlogging, and other natural disasters, causing serious casualties and economic losses. Flash floods caused by excessive rainfall are one of the most severe natural hazards in the world, presenting the most significant potential threat to the natural environment, human life, public health, and property (Kundzewicz et al., 2013). Since 1950, all very large and extreme floods in the Yangtze River basin (including the Huaihe River basin), such as in 1954, 1969, 1975, 1991,1998, 2003, 2007, 2010, and 2020, were caused by rainstorms, which caused huge losses to the country (Wang et al., 2003). Studies on rainstorm disasters include many aspects such as rainstorm flood disaster prediction, risk assessment of rainstorm disaster, zoning of flood disaster risk, and climatic characteristics of rainstorm disaster distribution.

    Because the type of disaster caused by rainstorms is generally regional flooding, flood forecasting technology is an interdisciplinary field. Generally, hydrological models are used in flood forecasting, including empirical models, lumped conceptual hydrologic models, and distributed hydrological models (Freeze and Harlan, 1969).A distributed hydrological model, considering the spatial heterogeneity of precipitation and land surface with digital elevation model (DEM) elevation data, can simulate the impact of human activities and underlying surface factors on the hydrological cycle. With the development of meteorological science and technology and the advancement of observation techniques, radar precipitation and model forecast precipitation have been gradually applied to hydrological models to improve the accuracy (Ren et al.,2003; Bao et al., 2011).

    Operational services for rainstorm disaster monitoring and prediction are in place in many countries. In rainstorm-caused flood monitoring, a hydrologic model is generally used. For example, the Australia Bureau of Meteorology flood forecasting and warning service uses rainfall and streamflow observations, numerical weather predictions, and hydrologic models to forecast possible flood events across Australia (http://www.bom.gov.au/australia/flood/). In India, during the flood season, the Flood Meteorological Offices (FMOs) provide valuable meteorological support to the Central Water Commission for issuing flood warnings in respect to many rivers.The real-time monitoring and statistical analysis of district-wide daily rainfall is one of the important functions of the Hydromet Division of India Meteorological department (IMD; https://mausam.imd.gov.in/). At the University of Maryland, USA, a Global Flood Monitoring System (GFMS) was developed using real-time satellite precipitation, hydrological runoff, and a routing model to create flood detection/intensity estimates, the latest maps of instantaneous precipitation and totals from the past day, three days, and seven days are also displayed (Wu et al., 2012, 2014). According to the above national operational flood monitoring systems, precipitation is a key factor in this kind of disaster.

    Due to the need for disaster prevention and mitigation,risk classification and evaluation of rainstorm and flood disaster losses have been a focus of academic study in recent years (Chou et al., 2013). Since the early 1990s,Chinese scholars have begun to research the risk evaluation of natural disaster losses and proposed evaluation standard, indicators, and methods for quantifying natural disaster losses (Liu, 1994; Wang et. al., 2007; Fu et al.,2016). According to the information on 263 rainstorms and flood disasters from 1989 to 1993, Li and Xu (1995)analyzed the economic losses in different disaster areas and the occurrence times of various grades of rainstorm and flood disaster. Hao et al. (2012) discovered a risk analysis method based on information diffusion theory.Jin et al. (2012) discussed a connection number of risk evaluation of China’s natural disaster systems based on triangular fuzzy numbers and stochastic simulation. Yao et al. (2012) analyzed the temporal and spatial distribution characteristics of rainstorms and flood disasters in Heilongjiang Province since 1990. Zhao et al. (2014) introduced the gray correlation method of risk evaluation in meteorological disaster losses based on historical disaster data from the Chinese mainland, and applied the improved gray relational analysis model (the triangular gray relational model) to the risk evaluation of rainstorm and flood disaster losses. At present, geographic information system (GIS) technology is widely used in the study of rainstorm disaster risk zoning based on rainstorm disaster assessment. The assessment of disasters includes not only precipitation, but also topography and other land surface conditions, as well as the population and economics. According to the operational requirements, a number of national and local standards for rainstorm disasters have been established.

    With global warming, research has shown that the global mean surface temperature has increased by 0.66°C from 1961–1990 to 2000–2019 (Valipour et al., 2021).The research on changes in temperature and precipitation over South Asia during the 21st century shows that there is still a continuous increase in the annual mean temperature and a robust increase in summer precipitation over most of the South Asian region (Almazroui et al., 2020). The continuous increase in temperature and precipitation will increase the frequency and intensity of rainstorm disasters in the future.

    Most of the rainstorm disaster assessment models and rain-caused flooding models are regionally applicable,which means that they need to be optimized separately for different regions and basins, and it is difficult to obtain the monitoring data needed in some models, such as hydrological data and ground observation rainfall. With the increasing demand for global rainstorm disaster monitoring services, it is necessary to build rainstorm disaster risk monitoring (RDRM) model based on regionwide coverage, high temporal resolution, and easily obtained datasets. With the development of satellite remote sensing technology, global real-time rainstorm observation can be realized by satellite retrieval rainfall, which provides the necessary conditions for real-time monitoring of global rainstorm disasters. Using rainfall data retrieved by multiple meteorological satellites, this paper focuses on the near-real-time monitoring of global severe rainstorm disaster risk with reference to rainstorm and flood disaster prediction and disaster assessment methods. In this paper, a worldwide dynamic rainstorm threshold map will be established according to historical data on severe rainstorm-caused flood disasters from the past 20 years. The paper will also establish the calculation method of rainstorm intensity index with multiple precipitation parameters. The factors of hazard-formative environments and hazard-affected bodies will also be considered. Finally, an RDRM model is established, and a rainstorm risk index (RRI) can be obtained in near real time. At the same time, the applicability of the RDRM in rainstorm disaster monitoring in different areas in 2020 is analyzed.

    2. Data and methods

    2.1 GSMaP_Gauge rainfall

    The Japanese Global Satellite Mapping of Precipitation (GSMaP) datasets are provided by the Japan Aerospace Exploration Agency (JAXA). They include a multisatellite global precipitation map under the Global Precipitation Measurement (GPM) Mission, using Dualfrequency Precipitation Radar (DPR) onboard GPM core satellites, other GPM constellation satellites, and geostationary satellites (Hou et al., 2014; Skofronick-Jackson et al., 2017; Deng et al., 2019; Gao et al., 2019). Three different remote sensing precipitation products are produced, GSMaP_NRT, GSMaP_MVK, and GSMaP_Gauge. In this study, the 24-h accumulated and hourly precipitation GSMaP_Gauge data, which include the correction of global surface rainfall station data, are used from 2000 to 2020. The variable is rainfall rate (mm h?1),the domain is global (60°N–60°S, 0–360°), and the grid resolution is 0.1° latitude/longitude.

    2.2 DFO flood disaster

    The global flood disaster dataset is from the Dartmouth Flood Observation (DFO) Center, Boulder, CO,USA (Kundzewicz et al., 2013). A sample of the data is given in Table 1, including the location, area, time, influence, main cause, and severity of the flood. There are 10 main causes of flooding disasters, 5 of them related to precipitation: heavy rain, tropical cyclones (TC), monsoonal rain, brief torrential rain, and extratropical cyclones. In this study, we mainly focus on the rain-induced kinds of disasters. The time range of the data used in this study is 2000 to 2020.

    The flood severity is defined in three classes, 1.0, 1.5,and 2.0, i.e., Class 1.0: large flood events, significant damage to structures or agriculture, fatalities, and/or a 1–2-decade interval since the last similar event; Class 1.5: very large events, with a greater than 20-yr but less than 100-yr recurrence interval, and/or a local recurrence interval of 10–20 yr; and Class 2.0: extreme events, with an estimated recurrence interval greater than 100 yr.

    2.3 ETOPO1 terrain height and river network

    ETOPO1 is a 1 arc-minute global relief model of Earth’s surface that integrates land topography and ocean bathymetry from the NOAA, Silver Spring, MD, USA(Hastings and Dunbar, 1999). Global river network data used in this study were provided by the National Geomatics Center of China, including the river information from level 1 to level 5.

    2.4 Gross domestic product and population

    The estimate of gross domestic product (GDP) is derived from nighttime lights satellite imagery and the LandScan population grid in 2010, which provide an alternative means for measuring economic activity (Ghosh et al., 2010), spatially distributed to generate a spatially disaggregated 1-km2map of total economic activity.

    The global population dataset in 2020 gives the estimated total number of people per grid-cell at a resolution of 1 km from WordPop (www.worldpop.org, https://dx.doi.org/10.5258/SOTON/WP00647; Stevens et al., 2015).

    2.5 Data normalization method

    The input datasets of the rainstorm disaster risk monitoring model built in this study include not only rainfall,but also terrain height, river network, economy (GDP),and population. They have different dimensions and units. In order to eliminate the influence of dimensions among influence factors in the model, data standardization is needed to solve the comparability between datasets. Two common normalization methods are Min–Max and Z-score (Carrino, 2017), Min–Max normalization is also known as deviation standardization, it is a linear transformation of the original data, so the result is mapped to 0–1, which in essence changes a number into a decimal between 0 and 1. The transfer function is (X?Min)/(Max ? Min), whereXis the data value to be normalized and Max and Min are the maximum and minimumvalues, respectively. After the normalization of the original or intermediate data, each influence factor is on the same order of magnitude and can be comprehensively analyzed in the model.

    Table 1. A sample of the flood disaster dataset

    2.6 Establishment of the rainstorm disaster risk monitoring model

    The RDRM model built in this study is an empirical model according to the current operational rainstorm disaster-related standards, including national standards(Chen et al., 2017), meteorological industry standards(Bao et al., 2018), and some provincial standards (Chen et al., 2017) in China. It also references many research results of rainstorm-caused flood disaster assessment and operational flood monitoring methods from many other countries (Wu et al., 2012, 2014).

    The prototype of the model is from the Zhejiang Provincial standard in China (Chen et al., 2017). The rainstorm process risk index is defined as follows:

    In this study, based on the rainstorm process risk assessment model, a rainstorm disaster risk real-time monitoring model will be built with multiple parameters referring to the original model and will be a great improvement in many aspects. Improvements include: (1) a dynamic global rainstorm disaster threshold map is established by using global climate data on rainstorm-induced floods and the precipitation retrieved from satellites in the last 20 years; (2) 7-day accumulated rainfall and the number of days when the daily rainfall exceeds the rainstorm disaster threshold in the first seven days are introduced to the new model; (3) gross domestic product(GDP) and population are considered as the rainstorm hazard-affected bodies influence coefficient; (4) using Min–Max data normalization to get more reasonable parameters.

    3. Results

    3.1 Distribution of rain-induced flood disasters

    According to the global flood disaster dataset of DFO from 2000 to 2020, there were about 3427 floods in the world during this period. Heavy rain caused 2111 floods,accounting for 61.6% of the total, which is the largest proportion. The numbers of brief torrential rain-induced,monsoonal rain-induced, and tropical cyclone/storm raininduced floods were 571, 323, and 312, respectively. The above four types of flooding disasters are most common,and are all related to precipitation, including extratropical cyclones, so the total number of flooding disasters related to precipitation is about 3319 from 2000 to 2020.The main causes of flooding disasters are shown in Fig. 1.

    According to the central location and areas of flooding disasters caused by precipitation (that is, 3319 cases),the distribution of flooding disasters is shown in Fig. 2.It can be seen that the regions of a high frequency of rain-induced flooding disasters in the world include Japan, North Korea, Central and Eastern China, Philippines, South and Southeast Asia, Central and East Africa,East Australia, Central and South Europe, southeastern North America, northwestern and southeastern South America, etc. Most of them belong to monsoon regions or tropical cyclone landing areas, or are areas affected by convective cloud systems in the intertropical convergence zone (ITCZ); and some of them are related to the major river basins in the world. For example, the high frequency of rain-induced flood disasters in South Asia is correlated with the summer monsoonal rainfall, tropical cyclones landing in the Arabian Sea and Bay of Bengal,and the Indus Basin and Ganges Basin.

    Fig. 1. Main causes of flooding disasters from 2000 to 2020.

    Fig. 2. Distribution of the frequency of global rain-induced flooding disasters from 2000 to 2020.

    3.2 Threshold for rainstorm disasters

    In the prototype of the model (Chen et al., 2017), the rainstorm threshold is a constant, that is, 50, which is the rainstorm standard for 24-h rainfall in eastern China.However, rainstorm disasters are affected by many factors, such as geographical location, topography,rivers, and climate events, resulting in obvious regional differences in the rainfall threshold of rainstorm disasters. Zhai et al. (2018) and Yuan et al. (2019) performed a comprehensive assessment and rechecking of the rainfall threshold for floods based on the disaster information, offering a valuable reference for the analysis,and calculation of early warning indicators of a flash flood. In this section, the distribution characteristics of rainfall threshold in different areas according to the flood disasters in 2000–2020 will be analyzed. This will provide a dynamic parameter of precipitation for the establishment of a rainstorm disaster monitoring model.

    The corrected distribution of global rainstorm threshold is shown in Fig. 3. It can be seen that the areas with a high rainfall threshold appear in the Korean Peninsula,Japan, Philippines, Malaysia, India, Madagascar, western and eastern Australia, the west coast of the United States, Northwest Mexico, and so on. When the flood disasters caused by TC rainfall are in the majority in a region, it is more common for there to be a large rainstorm threshold because the TC landing is characterized by strong rainfall. In the southwest of India, western Australia, western USA, and North Madagascar, the large rainstorm threshold is related to tropical storms. If continuous frontal precipitation-caused floods make up the majority, there is more likely to be a relatively small rainstorm threshold such as in Pakistan, Afghanistan,central Africa, and eastern Australia.

    3.3 Satellite data-based RDRM model

    Fig. 3. Distribution of the rainfall threshold for rainstorm disasters based on the flooding disasters during 2000–2020.

    Similar weather events may have very different impacts in different places. Thus, in the establishment of a rainstorm disaster monitoring model, not only the hazard-formative factors (that is, strong precipitation) but also hazard-formative environments and hazard-affected bodies are analyzed. Based on relevant national and industry standards of rainstorm disaster in China and monitoring methods of rainstorm and flood disasters (Chen H. Y. et al., 2017; Chen X. Y. et al., 2017; Bao et al.,2018), in this study, the rainstorm disaster risk index, including these three factors, is defined as follows:threshold is shown in Fig. 3. The rainstorm intensity index obtained by Eq. (3) is shown in Fig. 4d.

    Fig. 4. Distributions of precipitation parameters during the rainstorm disaster in Vietnam in October 2020: (a) precipitation ( R) on 6 October,(b) accumulated precipitation ( Rall) from 29 September to 5 October, (c) the number of days when daily precipitation ( Rd) exceeded the threshold from 29 September to 6 October, and (d) rainstorm intensity index ( If) on 6 October.

    The daily precipitation causing the rainstorm and flooding disaster occurred in central Vietnam, and was more than 50 mm in most areas. The accumulated precipitation in the first seven days in central and western Vietnam was more than 100 mm. The number of days when the precipitation exceeded the threshold was 1–2. The rainstorm intensity index in central Vietnam was greater than 1.2 when the rainstorm disaster occurred. In addition, the rainstorm intensity index on 6 October also showed that there was a high index greater than 1.2 in south central Cambodia, southern Thailand, and southern Myanmar. Disaster information was also reported on the internet (http://floodlist.com). It is said that the first of the recent spate of flooding began in early October 2020, when days of heavy rainfall in South East Asia caused flooding in Vietnam, Cambodia, Laos, and Thailand.

    Table 2. Values of the terrain influence factor ph drived based on the elevation h and its standard deviationSh

    Fig. 5. Distribution of terrain influence factor ph.

    whereSrmaxandSrminare the maximum and minimum ofSr, respectively. Thus, the value ofpris 0–1. The distribution of the global river network influence factor is shown in Fig. 6. It can be seen that the river influence coefficients in eastern and northeastern China, South Asia, Southeast Asia, and northern North America are high. The corresponding main rivers are the Yangtze River, Yellow River, Songhua River, Indus River,Ganges River, Mekong River, and Mississippi River.Compared with Fig. 2, the flood disaster frequency is also high in areas with a high value of the river influence factor.

    3.3.3Hazard-affected bodies influence coefficientIb′

    Severe rainstorm disaster will cause casualties or population transfer and economic losses. Therefore, two factors, population and GDP, are considered in the monitoring model of rainstorm disasters. If the impact of disaster prevention and reduction on disaster losses is not considered, it is generally believed that the more densely populated and economically developed areas are, the greater the losses caused by the same rainstorm. Based on this, the influence coefficient of rainstorm hazard-affected bodies is calculated as follows:

    whereppandwpare the population influence factor and its coefficient;pgandwgare the GDP influence factor and its coefficient. According to Chen et al. (2017), the sum of the coefficients of these two factors should be set to 1. If we consider which factor accounts for a larger proportion, then the coefficient will have a higher response. The distribution of the GDP influence factor and population influence factor are consistent with each other in most areas (see Sections 3.4.1 and 3.4.2). In this paper,wp= 0.5 andwg= 0.5 are selected.

    a. Population influence factorpp

    According to the global population density data, the global population distribution is extremely inhomogeneous, the maximum population density is 122,208 km–2in very few areas, while in low-population-density areas the number of people is less than 1 km–2. Therefore,when calculating the population influence factor through the Min–Max normalization method to deal with values between 0 and 1, the maximum value needs to be given a reasonable number in order to inform a more reasonable population influence coefficient. By analyzing the population distribution map and many experiments, finally,the maximum value (Max) in (X– Min)/(Max – Min) is selected as 800, that is, Max = 800 and Min = 0. The calculation method of population influence factor is as follows:

    wherepdis the population density. According to this relationship, the distribution of population influence coefficient can be calculated and is shown in Fig. 7. It can be seen that a population influence coefficient greater than 0.6 is mainly seen in east, southeast, and south Asia,western Europe, eastern and southern North America,central and eastern Africa, etc. In particular, northern India, Bangladesh, and eastern Pakistan in South Asia are the most densely populated areas with a population influence coefficient greater than 0.8. Rainstorm disasters occur frequently in these areas, which can be seen from the climate distribution of rainstorms and flooding disasters(Fig. 2).

    b. GDP influence factorpg

    There is a close relationship between economic development and population density. Studies on the relationship between economy and the population show that population density is one of the necessary conditions for economic activities and economic development. The demand pressure generated by population growth contributes to the generation of innovative behavior, thus promoting economic development and technological pro-

    Fig. 6. Distribution of the river network influence factor pr.

    Fig. 7. Distribution of population influence factor .

    pd gress. Generally speaking, the economy of densely populated areas will be more developed. However, the economic development is also affected by other factors. In the analysis of the GDP influence coefficient, the method is similar to that for the population impact coefficient.The maximum value of GDP in (X– Min)/(Max – Min)is selected as 50 (the maximum GDP is 723 in the dataset), that is, Max = 50 and Min = 0. The calculation method of the GDP influence factor is as follows:

    3.4 Evaluation and application in 2020

    The most important input data to the rainstorm disaster risk monitoring model concern rainfall; the main output of the model is the RRI, with which we hope to provide monitoring and early warning information for severe rainstorm disasters worldwide based on real-time satellite-retrieved rainfall. In order to analyze its reliability, the floods caused by rainstorms in 2020 were examined. We used flood information to study the RRI in corresponding areas and the active periods.

    3.4.1Evaluation

    In 2020, there were 172 floods worldwide in the DFO dataset, 146 of which were caused by rainfall (Fig. 9;Kundzewicz et al., 2013). The flood severity was classified into three categories in the DFO flood list (Table 1),that is, 1.0, 1.5, and 2.0, which are represented by green,blue, and red, respectively (Fig. 9). There were 39, 85,and 22 floods with a severity of 1.0, 1.5, and 2.0, respectively. The dots in Fig. 9 are the central locations of the 146 floods in 2020.

    Using the RDRM model established in this study, a daily RRI map was created for 2020. By examining the 146 floods and the corresponding RRI near the central location of each flood one by one in different flood severity classes, we obtained the results shown in Fig.10. The average maximum RRI for the 39 floods with severity 1.0 is about 1.68, which was the smallest among the three classes, and in most cases the RRI was between 1 and 2. There were five cases in which the maximum RRI were greater than 3, and one flood was caused by TC rainfall. The average maximum RRI were 2.16 and 3.63 for floods with severity 1.5 and 2.0, respectively.

    Fig. 8. Distribution of GDP influence factor pg.

    Fig. 9. Distribution of flood disaster central locations caused by rainstorms in 2020.

    It can be seen that, on average, the rainstorm risk index is consistent with the flood disaster severity for 2020(Fig. 10). However, there were some cases that do not match exactly. For example, in case No. 8 the flood severity was 2.0 but the maximum RRI was about 2.0,which is much smaller than the average of 3.63 (Fig.10c). The flood occurred in the United States in May; the main causes were a dam break and heavy rain. It may be that the dam break aggravated the influence of precipitation, leading a relatively small RRI to be associated with a stronger flood disaster. Another similar example is case No. 14 for flood severity 2.0 (Fig. 10c); the flood occurred in Chad, Africa, from August 11 to 28. Although the daily precipitation intensity was relatively small (the maximum RRI is about 2.0), long-lasting rainfall can still bring about an extreme flood event.

    The formation of a flood disaster is the result of a complex multifactor synthesis, although in the RDRM model, the effect of precipitation intensity on the disaster under the same daily precipitation is considered to some extent. There is still a lot of work to do in terms of detailed classification according to the disaster information to build a better RDRM model. However, this model reflects the disaster level of rainstorms to a certain extent.

    3.4.2Application in 2020

    It can be seen from the above analysis that RDRM model includes precipitation, topography and rivers, population, and economy, which can quantitatively represent the potential hazards caused by rainstorms. In order to analyze the reasonableness of the RRI in operational rainstorm service, two rainstorm-caused floods in 2020 were selected as examples.

    The water level of the Yangtze River and Poyang lake reservoirs increased, and the water area increased due to the heavy rainfall. According to the daily accumulated precipitation, since 29 June, the northern part of Jiangxi Province has been hit by continuous rainstorms, resulting in serious floods. During this period, the rainstorm from 7 to 8 July in Jiangxi Province, which is the most severe since complete meteorological records were made available in 1961, has had a wide range of influence.Based on satellite remote sensing and meteorological observation data, monitoring and evaluation of the changes in the main body and nearby water area of Poyang Lake in Jiangxi Province were carried out (Fig. 12). The result shows that the area of Poyang Lake is the largest it has been in the last 10 years (figure omitted).

    Fig. 10. The maximum RRI during each of the 146 floods caused by rainfall near the central location of the floods in 2020 for (a) the 39 floods with severity of 1.0, (b) the 85 floods with severity of 2.0, and (c) the 22 floods with severity of 2.0 based on the DFO flood disaster dataset.

    Using the RDRM model in this study, the RRI distribution was obtained as shown in Fig. 13. It can be seen that, affected by the continuous heavy rainfall, the regions with a RRI greater than 1.0 on 9 July were distributed in eastern Hubei, eastern Hunan, southern Anhui,and north–central Jiangxi provinces. The maximum center of the RRI, which was greater than 4.5, was located in the north–central part of Jiangxi Province. Compared with the precipitation (Fig. 11), the high risk index in southern Anhui and northern Jiangxi provinces had a close relationship with the strong 7-day accumulated precipitation, while the high RRI in central Jiangxi Province was related to the 24-h heavy rainfall. Figure 13 also shows that the value of the RRI in western Japan is 1–3 during to the eastward extension of the Meiyu front cloud belt. At the same time, affected by heavy rainfall, the rainstorm disaster risk index in Northeast India is also higher.

    The rainstorm disaster monitoring model is based on the historical data of rainstorm disasters, and a global dynamic threshold of rainstorm disaster is established;therefore, the rainstorm disaster risk index can be applied worldwide. The applicability of the index to other regions is illustrated by an example of monsoonal rainfall disaster in South Asia from late August 2020.

    In late August 2020, Pakistan and India suffered severe rainstorms due to the strong summer monsoonal rainfall. The rainstorm and flood disaster dataset (Table 1)shows that the rainstorm disaster in Pakistan occurred from 24 to 27 August, with the center located near 26.9573°N, 68.2157°E—that is, Southeast Pakistan, with an affected area of about 97,000 km2. The rainstorm disaster process in India occurred from 21 to 27 August, and the center is located near 24.7095°N, 77.4698°E, which is in the north central part of India, with an area of about 61,000 km2.

    The rainstorm disaster risk index (Fig. 14) shows that there were two high-risk centers, one in Sindh Province in the south of Pakistan and the other in Orissa state in the northeast of India, with rainstorm risk greater than 2.0. The maximum value in Southeast Pakistan was more than 4.5, while the maximum value of the RRI in Northeast India was 3–4.

    Form the comparison of water body monitoring before (5 August) and after (6 September) the floods in southern Pakistan byFY-3Dsatellite, it can be seen in Fig. 15 that, on 6 September, affected by the rainstorm,the Indus River, which flows through Punjab and Sindh,widened significantly, the water body area of Manchar Lake increased, and some areas such as Hyderabad, Mirpur Khas, and Umerkot in southern Sindh were affected by the rainstorm and floods.

    The rainstorm disaster risk index on 25 August was consistent with that of the rainstorm flood disaster dataset, and the surface water body monitoring by meteorological satellite also provided further evidence of the flooding disaster caused by the rainstorm. Therefore, the rainstorm disaster monitoring model established in this paper effectively monitored this rainstorm-induced disaster in South Asia.

    4. Conclusions and discussion

    Fig. 11. Distributions of precipitation during the rainstorm disasters in China in July 2020: (a) 7-day (2–8 July) accumulated precipitation,(b) 24-h precipitation on 9 July, and (c) R = max(×R3,×R12, R24)precipitation on 9 July 2020.

    In this study, an RDRM empirical model is built based on the remote sensing rainfall dataset, which can be used worldwide. Based on global rainstorm and flood information from the past 20 years, the dynamic threshold of rainstorm is established. Therefore, the difference in the precipitation intensity that can cause flooding disasters in different regions of the world is considered. Moreover,rainstorm disaster-causing factors include many parameters, and the hazard-formative environment factor and hazard-affected body factor of rainstorm disaster are built on terrain and river datasets, and population and economic datasets, respectively. The combination of the above factors makes the rainstorm disaster monitoring results have regional differences under similar precipitation conditions. The application of the rainstorm disaster monitoring model in 2020 in China, India, and Pakistan also shows the applicability of the model in different regions.The main contributions of this study are summarized as follows:

    Fig. 12. Map of the water body change (difference between 8 and 2 July) in and around Poyang Lake monitored by Fengyun satellites.

    Fig. 13. Distribution of the rainstorm disaster risk index over China on 9 July 2020.

    According to the global flooding disaster data from the past 20 years, the top four flood disaster causes (accounting for 96.8% in total) related to precipitation are heavy rain (accounting for 61.6%), brief torrential rain (16.7%),monsoonal rain (9.4%), and tropical cyclone/storm rain(9.1%). The distribution of flooding disasters is inhomogeneous; high-frequency areas are distributed in Japan,North Korea, central and eastern China, Philippines,South and Southeast Asia, Central and East Africa, East Australia, Central and South Europe, southeastern North America, northwestern and southeastern South America,etc.

    Based on data on the historical flood disasters caused by rainstorms and daily precipitation in the past 20 years,the worldwide dynamic threshold of precipitation during the floods is established, and the distribution of the precipitation threshold also shows regional non-uniformity.In the calculation of the rainstorm intensity index, a dynamic threshold is introduced, and the cumulative precipitation of the previous seven days, the daily precipitation,and the number of days when the precipitation exceeds the rainstorm disaster threshold are used as the hazardformative factors.

    Fig. 14. Distributions of the rainstorm disaster risk index on 25 August 2020: (a) global view and (b) regional view.

    Fig. 15. Water body change of the Indus River in Sindh, Pakistan monitored by FY-3D: (a) normal water level before the rainstorm on 5 August 2020 and (b) river widening due to the rainstorms on 6 September 2020.

    Two variables are selected for the hazard-formative environments (i.e., terrain and river networks) and two for the hazard-affected bodies (i.e., population distribution and GDP). By adding the influence factor into the rainstorm disaster monitoring model, the rainstorm intensity index can be modified to obtain a more reasonable RRI.

    By using the rainstorm disaster monitoring model established in this study, the long-lasting regional flood process in Southeast China and the monsoonal rain flood process in South Asia in 2020 are monitored. The rainstorm risk index monitors these two disaster processes well.

    The goal of this study was to establish a real-time rainstorm-caused flood monitoring model that can be used in daily operational services. Thus, the model was established according to the current operational rainstorm disaster-related standards, including national standards, meteorological industry standards, and some provincial standards of China. The design of the model also referenced many research results of rainstorm-caused flood disaster assessment and operational flood monitoring methods in other countries such as Australia, India,Pakistan, the USA, and so on. In the operational service of rainstorm disaster monitoring and forecasting, many parameters in the standard are set based on operational experience. Therefore, the parameters of the global rainstorm and flood disaster empirical model built in this study can be adjusted according to the specific needs.Because the rainstorm threshold in the model established in this study is based on information from the past 20 years on rainstorm-caused flooding disasters around the world, the model has a certain applicability worldwide.

    The GSMaP_Gauge daily rainfall data used in this study from the website have a three-day delay at present.In order to realize the real-time monitoring of global rainstorm disasters based on the meteorological satellite,there are higher requirements for the effectiveness of satellite precipitation retrieval products. Therefore, a precision assessment of satellite precipitation estimation products such as FY-2/3/4 satellites will be carried out later, determining the relationship between the precipitation estimation products of FY-2/3/4 and the GSMaP_Gauge, and the FY satellite precipitation data will be corrected in its scanning areas (Singh et al., 2018). According to the rainstorm disaster monitoring model established in this study, the RRI monitoring products of global key regions (such as in the Belt and Road regions) will be established based on the FY satellites, which can be applied to the monitoring and early warning services of heavy rain disasters in the key regions of the world.

    亚洲aⅴ乱码一区二区在线播放| 身体一侧抽搐| 黄片无遮挡物在线观看| 成人特级av手机在线观看| 性色avwww在线观看| 欧美日韩亚洲高清精品| 日韩国内少妇激情av| 亚洲,欧美,日韩| 夫妻性生交免费视频一级片| 看非洲黑人一级黄片| 国产一区二区亚洲精品在线观看| 麻豆国产97在线/欧美| 校园人妻丝袜中文字幕| 97热精品久久久久久| 日本黄色片子视频| 91精品一卡2卡3卡4卡| 纵有疾风起免费观看全集完整版| 又粗又硬又长又爽又黄的视频| 日日摸夜夜添夜夜爱| 韩国高清视频一区二区三区| 国产伦在线观看视频一区| 欧美精品人与动牲交sv欧美| 有码 亚洲区| 久久国内精品自在自线图片| 国产免费又黄又爽又色| 边亲边吃奶的免费视频| 久久久久久国产a免费观看| 国产片特级美女逼逼视频| 涩涩av久久男人的天堂| 国产精品久久久久久久电影| 菩萨蛮人人尽说江南好唐韦庄| 日本av手机在线免费观看| 久久久久精品性色| 99热这里只有是精品在线观看| 日韩制服骚丝袜av| 日韩一本色道免费dvd| 亚洲国产精品专区欧美| 一边亲一边摸免费视频| 一区二区三区免费毛片| 国产人妻一区二区三区在| 国产精品国产三级国产专区5o| 又爽又黄无遮挡网站| 女的被弄到高潮叫床怎么办| tube8黄色片| 免费av不卡在线播放| 国产一区二区在线观看日韩| 一级毛片我不卡| 99热这里只有精品一区| 熟女电影av网| 高清av免费在线| 在现免费观看毛片| 成人午夜精彩视频在线观看| 蜜臀久久99精品久久宅男| 97在线视频观看| 亚洲国产色片| 午夜免费男女啪啪视频观看| 国产 一区精品| 一级毛片久久久久久久久女| 成年av动漫网址| 国产精品一区二区三区四区免费观看| 日韩精品有码人妻一区| 国产精品久久久久久av不卡| 观看免费一级毛片| 中文天堂在线官网| 亚洲最大成人av| 好男人在线观看高清免费视频| 99热这里只有是精品50| 亚洲精品亚洲一区二区| 少妇裸体淫交视频免费看高清| 三级国产精品片| 中文资源天堂在线| 26uuu在线亚洲综合色| 91在线精品国自产拍蜜月| 色视频www国产| 国产伦理片在线播放av一区| 日韩伦理黄色片| 久久99热这里只频精品6学生| 国产精品蜜桃在线观看| 欧美区成人在线视频| 久久久成人免费电影| 六月丁香七月| 乱码一卡2卡4卡精品| 99视频精品全部免费 在线| 国产白丝娇喘喷水9色精品| 一级毛片黄色毛片免费观看视频| 国产国拍精品亚洲av在线观看| 亚洲精品日本国产第一区| 日韩视频在线欧美| 日本与韩国留学比较| 99九九线精品视频在线观看视频| 熟妇人妻不卡中文字幕| 亚洲电影在线观看av| 免费电影在线观看免费观看| 欧美高清性xxxxhd video| 日本一本二区三区精品| 久久精品人妻少妇| 亚洲av电影在线观看一区二区三区 | 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 联通29元200g的流量卡| 国产永久视频网站| 成人漫画全彩无遮挡| 麻豆精品久久久久久蜜桃| 麻豆成人av视频| 一区二区三区精品91| 国产亚洲最大av| 国产精品99久久99久久久不卡 | 久久人人爽av亚洲精品天堂 | 精品酒店卫生间| 亚洲怡红院男人天堂| 亚洲高清免费不卡视频| 中文在线观看免费www的网站| 国产伦理片在线播放av一区| 自拍偷自拍亚洲精品老妇| 综合色丁香网| 又爽又黄a免费视频| 国产午夜福利久久久久久| 国产国拍精品亚洲av在线观看| 午夜激情福利司机影院| 大香蕉久久网| 国产成人免费观看mmmm| 国产黄色视频一区二区在线观看| 麻豆成人av视频| 免费播放大片免费观看视频在线观看| 日韩亚洲欧美综合| 国产精品99久久久久久久久| 日日撸夜夜添| 国产黄色视频一区二区在线观看| 岛国毛片在线播放| 国产精品人妻久久久影院| 亚洲国产色片| www.av在线官网国产| 七月丁香在线播放| 三级国产精品片| 色综合色国产| 一级片'在线观看视频| 少妇 在线观看| 色5月婷婷丁香| 婷婷色综合大香蕉| 麻豆成人午夜福利视频| 不卡视频在线观看欧美| 国产精品一及| 久久精品久久久久久久性| 看免费成人av毛片| 国产黄a三级三级三级人| 国产精品一二三区在线看| 国产成人freesex在线| 亚洲伊人久久精品综合| 亚洲精品日本国产第一区| 国产伦理片在线播放av一区| 18禁裸乳无遮挡动漫免费视频 | 26uuu在线亚洲综合色| 亚洲精品国产av成人精品| 国产在线一区二区三区精| 在线播放无遮挡| 大香蕉久久网| 精品人妻偷拍中文字幕| 看非洲黑人一级黄片| 国产男女内射视频| 欧美性猛交╳xxx乱大交人| 中文乱码字字幕精品一区二区三区| 亚洲国产精品国产精品| 九草在线视频观看| eeuss影院久久| 超碰av人人做人人爽久久| 狂野欧美激情性bbbbbb| 中文在线观看免费www的网站| 久久久久国产网址| 自拍偷自拍亚洲精品老妇| av一本久久久久| 午夜福利视频1000在线观看| 国产乱人偷精品视频| 久久久久久久久久久免费av| 午夜福利在线在线| av又黄又爽大尺度在线免费看| 啦啦啦中文免费视频观看日本| 国产乱来视频区| 成人综合一区亚洲| 国产白丝娇喘喷水9色精品| 97人妻精品一区二区三区麻豆| av一本久久久久| 男男h啪啪无遮挡| 亚洲国产成人一精品久久久| 国产大屁股一区二区在线视频| 欧美老熟妇乱子伦牲交| 国产熟女欧美一区二区| 国产精品久久久久久久电影| 只有这里有精品99| 亚州av有码| 国产精品精品国产色婷婷| 一本一本综合久久| 国产亚洲一区二区精品| 国产精品一二三区在线看| 精品视频人人做人人爽| 夫妻性生交免费视频一级片| 神马国产精品三级电影在线观看| 国产午夜精品久久久久久一区二区三区| 久久人人爽人人爽人人片va| 能在线免费看毛片的网站| 亚洲人成网站在线观看播放| 亚州av有码| 又爽又黄无遮挡网站| 你懂的网址亚洲精品在线观看| 亚洲精品国产色婷婷电影| 成人亚洲精品av一区二区| 国产熟女欧美一区二区| 国产亚洲最大av| 国产在线一区二区三区精| 久久精品人妻少妇| 美女cb高潮喷水在线观看| 国产精品不卡视频一区二区| 精品午夜福利在线看| 美女高潮的动态| 国产老妇女一区| 青春草亚洲视频在线观看| 久久精品国产鲁丝片午夜精品| 日韩av在线免费看完整版不卡| 国产免费视频播放在线视频| 在线播放无遮挡| 丰满人妻一区二区三区视频av| 日韩亚洲欧美综合| 99热这里只有是精品50| 国产精品精品国产色婷婷| 国产亚洲最大av| 99热全是精品| kizo精华| 少妇人妻一区二区三区视频| 国产高清国产精品国产三级 | 色播亚洲综合网| 国产美女午夜福利| 联通29元200g的流量卡| 大话2 男鬼变身卡| 丰满人妻一区二区三区视频av| 亚洲国产精品成人综合色| 蜜桃亚洲精品一区二区三区| 99热这里只有是精品50| 精品熟女少妇av免费看| 免费av毛片视频| 午夜福利视频精品| 精品久久久久久久人妻蜜臀av| 精品久久久久久久久亚洲| 超碰av人人做人人爽久久| av线在线观看网站| 欧美精品人与动牲交sv欧美| 超碰av人人做人人爽久久| 国产乱人偷精品视频| 老师上课跳d突然被开到最大视频| 看十八女毛片水多多多| 在线观看国产h片| 高清在线视频一区二区三区| 亚洲精品日韩av片在线观看| 国产精品一区www在线观看| 亚洲av二区三区四区| 国产精品爽爽va在线观看网站| av黄色大香蕉| 亚洲欧美精品专区久久| 国产成人aa在线观看| 黄色日韩在线| 嫩草影院新地址| 欧美潮喷喷水| 国产亚洲91精品色在线| 国产乱人偷精品视频| 黄片无遮挡物在线观看| 精品一区二区三区视频在线| 精品国产露脸久久av麻豆| 青青草视频在线视频观看| 青春草视频在线免费观看| 最新中文字幕久久久久| 狂野欧美激情性bbbbbb| 麻豆久久精品国产亚洲av| 亚洲色图综合在线观看| 国产毛片在线视频| 亚洲最大成人av| 国产一区有黄有色的免费视频| 99热国产这里只有精品6| 日韩 亚洲 欧美在线| 亚洲成人av在线免费| 国产亚洲最大av| 精品久久久噜噜| 精品久久久久久久久亚洲| 免费在线观看成人毛片| 天天躁夜夜躁狠狠久久av| 亚洲国产最新在线播放| 欧美激情国产日韩精品一区| 欧美少妇被猛烈插入视频| videos熟女内射| 国产 一区精品| 在线免费观看不下载黄p国产| 国产熟女欧美一区二区| 精品久久久精品久久久| 国产白丝娇喘喷水9色精品| 青青草视频在线视频观看| 成人午夜精彩视频在线观看| 亚洲欧洲日产国产| videos熟女内射| 久久精品久久久久久噜噜老黄| 免费在线观看成人毛片| 成人毛片a级毛片在线播放| 一区二区三区免费毛片| av福利片在线观看| 少妇猛男粗大的猛烈进出视频 | 国模一区二区三区四区视频| 波野结衣二区三区在线| 丝袜美腿在线中文| 国产一区二区三区综合在线观看 | 国产精品一二三区在线看| 亚洲精品aⅴ在线观看| 欧美成人一区二区免费高清观看| 伦理电影大哥的女人| 听说在线观看完整版免费高清| 精品一区二区三卡| av又黄又爽大尺度在线免费看| 尤物成人国产欧美一区二区三区| 成人黄色视频免费在线看| 麻豆乱淫一区二区| 久久99精品国语久久久| 99热全是精品| 国产精品国产三级专区第一集| 在线 av 中文字幕| 在线看a的网站| 久久久久九九精品影院| 精品国产一区二区三区久久久樱花 | 3wmmmm亚洲av在线观看| 中文字幕亚洲精品专区| 五月伊人婷婷丁香| 久热久热在线精品观看| 日韩欧美精品v在线| 久久久久国产精品人妻一区二区| 午夜日本视频在线| 最近的中文字幕免费完整| 大话2 男鬼变身卡| 在线免费十八禁| 肉色欧美久久久久久久蜜桃 | 国内少妇人妻偷人精品xxx网站| 美女cb高潮喷水在线观看| 大陆偷拍与自拍| 国产精品久久久久久av不卡| 国产精品一区二区在线观看99| 亚洲aⅴ乱码一区二区在线播放| 人人妻人人爽人人添夜夜欢视频 | 欧美日韩在线观看h| a级毛色黄片| 亚洲精品日韩在线中文字幕| 亚洲精品第二区| 男的添女的下面高潮视频| 中文乱码字字幕精品一区二区三区| 波野结衣二区三区在线| 午夜老司机福利剧场| 国产白丝娇喘喷水9色精品| 欧美zozozo另类| 精品少妇黑人巨大在线播放| 久久99热这里只频精品6学生| 久久亚洲国产成人精品v| 久久ye,这里只有精品| 美女被艹到高潮喷水动态| 亚洲自偷自拍三级| 一级毛片aaaaaa免费看小| 九色成人免费人妻av| av在线app专区| 亚洲欧美成人综合另类久久久| 极品教师在线视频| 黄片wwwwww| 美女视频免费永久观看网站| 最新中文字幕久久久久| 又爽又黄a免费视频| 美女国产视频在线观看| 久久久久久久亚洲中文字幕| 中文精品一卡2卡3卡4更新| 又黄又爽又刺激的免费视频.| 高清日韩中文字幕在线| 国产成人a区在线观看| 亚洲一级一片aⅴ在线观看| 高清在线视频一区二区三区| 亚洲av免费在线观看| 国产亚洲av片在线观看秒播厂| 26uuu在线亚洲综合色| 下体分泌物呈黄色| 一级毛片aaaaaa免费看小| av国产久精品久网站免费入址| 寂寞人妻少妇视频99o| 男女边吃奶边做爰视频| 五月天丁香电影| 插逼视频在线观看| 综合色丁香网| 亚洲三级黄色毛片| 亚洲av成人精品一二三区| 国产精品嫩草影院av在线观看| 国产又色又爽无遮挡免| 我的女老师完整版在线观看| 日韩一区二区三区影片| 欧美精品人与动牲交sv欧美| 日韩免费高清中文字幕av| 十八禁网站网址无遮挡 | 久久精品综合一区二区三区| 欧美zozozo另类| 国产欧美亚洲国产| 黄色怎么调成土黄色| 秋霞伦理黄片| 亚洲av男天堂| 久久人人爽av亚洲精品天堂 | 99久国产av精品国产电影| 国产 一区 欧美 日韩| freevideosex欧美| 欧美日韩国产mv在线观看视频 | 极品少妇高潮喷水抽搐| 亚洲精品亚洲一区二区| 嫩草影院精品99| av黄色大香蕉| 久久99热这里只有精品18| 日韩一区二区三区影片| 久久久精品免费免费高清| 我的女老师完整版在线观看| 嘟嘟电影网在线观看| 中文字幕av成人在线电影| 亚洲欧洲日产国产| 国产伦精品一区二区三区视频9| 亚洲欧美清纯卡通| 美女xxoo啪啪120秒动态图| 国内揄拍国产精品人妻在线| 狂野欧美激情性bbbbbb| 亚洲成色77777| 日本av手机在线免费观看| 亚洲在线观看片| 成人二区视频| 日本熟妇午夜| 久久久久久伊人网av| 亚洲av不卡在线观看| 亚洲精品,欧美精品| 久久国内精品自在自线图片| 国产伦在线观看视频一区| 日韩三级伦理在线观看| 99久久人妻综合| 国产成人精品福利久久| av天堂中文字幕网| 精品人妻视频免费看| 国产高清国产精品国产三级 | 久久97久久精品| 一级毛片我不卡| 大香蕉久久网| 久久久久久国产a免费观看| 丝袜喷水一区| 国产伦精品一区二区三区视频9| 晚上一个人看的免费电影| 在线精品无人区一区二区三 | 3wmmmm亚洲av在线观看| 亚洲精品国产av成人精品| 国产黄片美女视频| 亚洲精品久久午夜乱码| 2018国产大陆天天弄谢| 精品亚洲乱码少妇综合久久| 神马国产精品三级电影在线观看| 国内精品宾馆在线| 日韩成人伦理影院| 国产成人精品婷婷| 亚洲精品一区蜜桃| 久久国产乱子免费精品| 51国产日韩欧美| 欧美+日韩+精品| 色哟哟·www| 97在线人人人人妻| 国产永久视频网站| 看非洲黑人一级黄片| 最近中文字幕高清免费大全6| 99久久人妻综合| 久久久午夜欧美精品| 一区二区三区乱码不卡18| kizo精华| av又黄又爽大尺度在线免费看| 国产乱人偷精品视频| 22中文网久久字幕| 日本一二三区视频观看| 三级经典国产精品| 青青草视频在线视频观看| 午夜视频国产福利| 亚洲av.av天堂| 国产v大片淫在线免费观看| 日韩一本色道免费dvd| 成人无遮挡网站| 欧美+日韩+精品| 国产精品一及| 国产精品av视频在线免费观看| 亚洲av不卡在线观看| 男女下面进入的视频免费午夜| 青春草国产在线视频| 天天躁夜夜躁狠狠久久av| 亚洲人与动物交配视频| 欧美亚洲 丝袜 人妻 在线| 久久久久性生活片| av在线观看视频网站免费| 热re99久久精品国产66热6| 欧美日韩视频精品一区| 中文字幕制服av| 亚洲国产av新网站| 又爽又黄a免费视频| av又黄又爽大尺度在线免费看| 精品视频人人做人人爽| 日韩三级伦理在线观看| 国产免费又黄又爽又色| www.色视频.com| 老师上课跳d突然被开到最大视频| 一级毛片 在线播放| 免费黄频网站在线观看国产| 一级片'在线观看视频| 亚洲av成人精品一区久久| 麻豆国产97在线/欧美| 欧美日本视频| 国产综合精华液| 免费看不卡的av| 日韩av免费高清视频| 久久久久久久亚洲中文字幕| 久久久成人免费电影| 中文在线观看免费www的网站| 3wmmmm亚洲av在线观看| 国产精品久久久久久久久免| 免费观看av网站的网址| 国产免费视频播放在线视频| 国产探花极品一区二区| 91久久精品国产一区二区三区| 精品一区二区三卡| 亚洲精品乱久久久久久| 国产乱人偷精品视频| 精品国产三级普通话版| 成人鲁丝片一二三区免费| 91午夜精品亚洲一区二区三区| 午夜福利视频1000在线观看| 一个人看的www免费观看视频| 欧美性猛交╳xxx乱大交人| 又黄又爽又刺激的免费视频.| 国产成人91sexporn| 中文字幕免费在线视频6| 亚洲自拍偷在线| 建设人人有责人人尽责人人享有的 | 国产真实伦视频高清在线观看| 少妇 在线观看| 麻豆久久精品国产亚洲av| 国产亚洲av嫩草精品影院| 老司机影院毛片| 美女内射精品一级片tv| 成人免费观看视频高清| 国产高清国产精品国产三级 | 日韩一区二区视频免费看| 在线观看三级黄色| 亚洲av成人精品一区久久| xxx大片免费视频| 成人美女网站在线观看视频| 亚洲av免费在线观看| 精品久久久噜噜| 亚洲内射少妇av| 午夜视频国产福利| 国产亚洲最大av| 高清视频免费观看一区二区| 欧美性感艳星| 永久网站在线| 干丝袜人妻中文字幕| 亚洲av福利一区| 国产精品国产av在线观看| 亚洲aⅴ乱码一区二区在线播放| 国产女主播在线喷水免费视频网站| 国产精品伦人一区二区| 中文精品一卡2卡3卡4更新| 国产毛片a区久久久久| 最近2019中文字幕mv第一页| 久久久久久久久大av| 国产又色又爽无遮挡免| 99九九线精品视频在线观看视频| 人妻制服诱惑在线中文字幕| 99久久精品国产国产毛片| 2022亚洲国产成人精品| 久久久久久久大尺度免费视频| 日韩av免费高清视频| 人人妻人人澡人人爽人人夜夜| 久久久久性生活片| 国产成年人精品一区二区| 久久鲁丝午夜福利片| 国产美女午夜福利| 国产成人免费观看mmmm| 人妻一区二区av| 亚洲va在线va天堂va国产| 久久国产乱子免费精品| 黄色欧美视频在线观看| 天美传媒精品一区二区| 国产永久视频网站| av黄色大香蕉| 亚洲精品久久午夜乱码| 久久久亚洲精品成人影院| 尤物成人国产欧美一区二区三区| 国产综合精华液| 亚洲精品一区蜜桃| 免费大片18禁| 亚洲精品aⅴ在线观看| 国产精品国产三级国产专区5o| 久久国产乱子免费精品| 久久精品人妻少妇| 水蜜桃什么品种好| 亚洲真实伦在线观看| 免费看日本二区| 少妇被粗大猛烈的视频| 日本三级黄在线观看| 97超视频在线观看视频| 国产精品麻豆人妻色哟哟久久| 啦啦啦在线观看免费高清www| 成人美女网站在线观看视频| 国产亚洲午夜精品一区二区久久 | 国产女主播在线喷水免费视频网站| 狂野欧美激情性bbbbbb| 亚洲av国产av综合av卡| 干丝袜人妻中文字幕| 亚洲av成人精品一二三区| 国产探花在线观看一区二区| 搡老乐熟女国产| 国产成人freesex在线| 久久人人爽av亚洲精品天堂 | 色网站视频免费| 高清欧美精品videossex| 国内少妇人妻偷人精品xxx网站| 国产乱来视频区|