Incorporating elevation in rainfall interpolation in Tunisia using geostatistical methods

Abstract

This paper compares the performance of three geostatistical algorithms, which integrate elevation as an auxiliary variable: kriging with external drift (KED); kriging combined with regression, called regression kriging (RK) or kriging after detrending; and co-kriging (CK). These three methods differ by the way by in which the secondary information is introduced into the prediction procedure. They are applied to improve the prediction of the monthly average rainfall observations measured at 106 climatic stations in Tunisia over an area of 164 150 km² using the elevation as the auxiliary variable. The experimental sample semivariograms, residual semivariograms and cross-variograms are constructed and fitted to estimate the rainfall levels and the estimation variance at the nodes of a square grid of 20 km?×?20 km resolution and to develop corresponding contour maps. Contour diagrams for KED and RK were similar and exhibited a pattern corresponding more closely to local topographic features when (a) the network is sparse and (b) the rainfall–elevation correlation is poor, while CK showed a smooth zonal pattern. Smaller prediction variances are obtained for the RK algorithm. The cross-validation showed that the RMSE obtained for CK gave better results than for KED or RK.

Editor D. Koutsoyiannis; Associate editor C. Onof

Citation Feki, H., Slimani, M., and Cudennec, C., 2012. Incorporating elevation in rainfall interpolation in Tunisia using geostatistical methods. Hydrological Sciences Journal, 57 (7), 1294–1314.     

A study on the spatial variations in stable isotopic composition of precipitation in a semiarid region of Southern India

The stable isotope analysis of all major rain events from Moinabad (MB), Rajendranagar (RN) and Osmanasagar (OS) reservoir, three closely placed locations in Hyderabad, India, were carried out during the 2005 to 2008 period. The OS station recorded the highest amount of rainfall with an average value of 1000 mm, whereas the MB station recorded the lowest average rainfall of 790 mm. The stable isotope (δ¹⁸O) values of the precipitation samples during these period varied from ?11.43‰ to ?0.03‰ for the MB station, ?8.21‰ to 0.54‰ for the RN station and ?11.47‰ to 0.72‰ for the OS station. The d‐excess of precipitation at the three stations also showed considerable variations and revealed that the precipitation in the region undergoes significant modification through secondary evaporation during its fall. The possible causes for these observed spatial and temporal variations in amount and the isotopic composition of precipitation in a small geographical area within the city were studied. The observed variations may be attributed to the regional scale differences in water budget induced by rapid urbanisation activities in the city coupled with the differences in secondary effects undergone by the falling drops. This study elucidating changes in precipitation patterns in the city and its possible causes may largely help in its water balance calculation. Copyright © 2011 John Wiley & Sons, Ltd.     

A temperature‐precipitation‐based model of thirty‐year mean snowpack accumulation and melt in Oregon,USA

High‐resolution, spatially extensive climate grids can be useful in regional hydrologic applications. However, in regions where precipitation is dominated by snow, snowmelt models are often used to account for timing and magnitude of water delivery. We developed an empirical, nonlinear model to estimate 30‐year means of monthly snowpack and snowmelt throughout Oregon. Precipitation and temperature for the period 1971–2000, derived from 400‐m resolution PRISM data, and potential evapotranspiration (estimated from temperature and day length) drive the model. The model was calibrated using mean monthly data from 45 SNOTEL sites and accurately estimated snowpack at 25 validation sites: R² = 0·76, Nash‐Sutcliffe Efficiency (NSE) = 0·80. Calibrating it with data from all 70 SNOTEL sites gave somewhat better results (R² = 0·84, NSE = 0·85). We separately applied the model to SNOTEL stations located < 200 and ≥ 200 km from the Oregon coast, since they have different climatic conditions. The model performed equally well for both areas. We used the model to modify moisture surplus (precipitation minus potential evapotranspiration) to account for snowpack accumulation and snowmelt. The resulting values accurately reflect the shape and magnitude of runoff at a snow‐dominated basin, with low winter values and a June peak. Our findings suggest that the model is robust with respect to different climatic conditions, and that it can be used to estimate potential runoff in snow‐dominated basins. The model may allow high‐resolution, regional hydrologic comparisons to be made across basins that are differentially affected by snowpack, and may prove useful for investigating regional hydrologic response to climate change. Published in 2011 by John Wiley & Sons, Ltd.     

Geostatistical radar–raingauge merging: A novel method for the quantification of rain estimation accuracy

Compared to other estimation techniques, one advantage of geostatistical techniques is that they provide an index of the estimation accuracy of the variable of interest with the kriging estimation standard deviation (ESD). In the context of radar–raingauge quantitative precipitation estimation (QPE), we address in this article the question of how the kriging ESD can be transformed into a local spread of error by using the dependency of radar errors to the rain amount analyzed in previous work. The proposed approach is implemented for the most significant rain events observed in 2008 in the Cévennes-Vivarais region, France, by considering both the kriging with external drift (KED) and the ordinary kriging (OK) methods. A two-step procedure is implemented for estimating the rain estimation accuracy: (i) first kriging normalized ESDs are computed by using normalized variograms (sill equal to 1) to account for the observation system configuration and the spatial structure of the variable of interest (rainfall amount, residuals to the drift); (ii) based on the assumption of a linear relationship between the standard deviation and the mean of the variable of interest, a denormalization of the kriging ESDs is performed globally for a given rain event by using a cross-validation procedure. Despite the fact that the KED normalized ESDs are usually greater than the OK ones (due to an additional constraint in the kriging system and a weaker spatial structure of the residuals to the drift), the KED denormalized ESDs are generally smaller the OK ones, a result consistent with the better performance observed for the KED technique. The evolution of the mean and the standard deviation of the rainfall-scaled ESDs over a range of spatial (5–300 km²) and temporal (1–6 h) scales demonstrates that there is clear added value of the radar with respect to the raingauge network for the shortest scales, which are those of interest for flash-flood prediction in the considered region.     

Spatial variation in soil resistance to flowing water erosion along a regional transect in the Loess Plateau

The factors influencing soil erosion may vary with scale. It remains unclear whether the spatial variation in soil erosion resistance is controlled by regional variables (e.g. precipitation, temperature, and vegetation zone) or by local specific variables (e.g. soil properties, root traits, land use, and farming operations) when the study area enlarges from a hillslope or catchment to the regional scale. This study was performed to quantify the spatial variations in soil erosion resistance to flowing water under three typical land uses along a regional transect on the Loess Plateau and to identify whether regional or local specific variables are responsible for these changes. The results indicated that the measured soil detachment capacities (D_c) of cropland exhibited an irregular trend along the regional transect. The D_c of grassland increased with mean annual precipitation, except for two sites (Yijun and Erdos). The measured D_c of woodland displayed an inverted ‘U’ shape. The changes in rill erodibility (K_r) of three land uses were similar to D_c, whereas no distinguishable trend was found for critical shear stress (τ_c). No significant correlation was detected between D_c, K_r and τ_c, and the regional variables. The spatial variation in soil erosion resistance could be explained reasonably by changes in soil properties, root traits, land use, and farming operations, rather than regional variables. The adjustment coefficient of K_r for grassland and woodland could be well simulated by soil cohesion and root mass density (R² = 0.70, P < 0.01), and the adjustment coefficient of critical shear stress could be estimated with aggregate stability (R² = 0.57, P < 0.01). The results are helpful for quantifying the spatial variation in soil detachment processes by overland flow and to develop process‐based erosion model at a regional scale. Copyright © 2015 John Wiley & Sons, Ltd.     

Water and nutrient predictions in ungauged basins: set-up and evaluation of a model at the national scale

Abstract

A dynamic water quality model, HYPE, was applied to a large, data-sparse region to study whether reliable information on water quantity and water quality could be obtained for both gauged and ungauged waterbodies. The model (called S-HYPE) was set up for all of Sweden (～450 000 km²), divided into sub-basins with an average area of 28 km². Readily available national databases were used for physiographic data, emissions and agricultural practices, fixed values for representative years were used. Daily precipitation and temperature were used as the dynamic forcing of the model. Model evaluation was based on data from several hundred monitoring sites, of which approximately 90% had not been used in calibration on a daily scale. Results were evaluated using the Nash-Sutcliffe efficiency (NSE), correlation and relative errors: 92% of the spatial variation was explained for specific water discharge, and 88% and 59% for total nitrogen and total phosphorus concentrations, respectively. Day-to-day variations were modelled with satisfactory results for water discharge and the seasonal variation of nitrogen concentrations was also generally well captured. In 20 large, unregulated rivers the median NSE for water discharge was 0.84, and the corresponding number for 76 partly-regulated river basins was 0.52. In small basins, the NSE was typically above 0.6. These major achievements relative to previous similar experiments were ascribed to the step-wise calibration process using representative gauged basins and the use of a modelling concept, whereby coefficients are linked to physiographic variables rather than to specific sites.

Editor D. Koutsoyiannis

Citation Strömqvist, J., Arheimer, B., Dahné, J., Donnelly, C. and Lindström, G., 2012. Water and nutrient predictions in ungauged basins: set-up and evaluation of a model at the national scale. Hydrological Sciences Journal, 57 (2), 229–247.     

BIBLIOGRAPHIES

Abstract

On a regional scale (from about 10 to 10⁵km²) estimations of actual evapotranspiration ET for water balance equations are generally missing. The various methods available at present are briefly described and discussed. The use of remote sensing, in particular infra-red thermometry by satellite, in conjunction with soil energy budget data, is the only method to provide satisfactory data on the same scale in space and lime for water balance models. This method is being studied at the INRA Bioclimatology Station at Avignon in the framework of the European ‘Tellus’ project for using HCMM satellite data. Systems for continuously recording the components of the energy budget and ET have been installed in the Crau plain, north of Marseilles, in both a dry area and a large area of irrigated pastures. By combining the recorded measurements with the satellite data and starting from values for elements of 25 ha (500 m × 500 m), the ET for the whole of the Crau plain can be estimated. First results are analysed and the possibilities and limitations of the method are discussed.     

Comparing various approaches for assessing groundwater recharge at a regional scale in the Canadian Shield

ABSTRACT

Estimating groundwater recharge is crucial to ensuring the proper management of aquifers. In this study, net regional recharge and spatial potential recharge are estimated at four watersheds within the Charlevoix–Haute-Côte-Nord (CHCN) regions, Quebec Province, Canada. Four methods are applied based on available data. The first two approaches are regional water budget methods. These two methods differ in their estimation of vertical inflow (VI), which is estimated from two hydrological models: GR4J and HYDROTEL. The third method estimates potential recharge spatially over the study area. Finally, the streamflow data are analysed using the Eckhardt baseflow separation method to obtain an estimation of recharge, assuming that discharge is equal to recharge. According to the results of all investigated methods, the mean annual recharge for the CHCN region is approximately 183 mm, which is 18% of the total annual precipitation (P). The discussion section highlights uncertainties due to the assumptions of each method and the reliability of the data.     

An integrated approach of flash flood analysis in ungauged Mediterranean watersheds using post-flood surveys and unmanned aerial vehicles

This study analyzes the flash flood event of two ungauged ephemeral streams in Olympiada region (Chalkidiki, North Greece), which occurred at the 21–22 of November 2019. Aim of the study is to reconstruct the specific flash flood event, investigate the causes of flood generation mechanisms, evaluate the performance of SCS-CN hydrological and HEC-RAS hydraulic models, investigate the relation between extreme flash floods and human intervention, using the combination of ground and aerial observations obtained from the field survey and unmanned aerial vehicles (UAVs), respectively. The results of the specific discharge ranged between 9 and 11 m³ s⁻¹ km², values that are typical for flash flood events in Mediterranean region. The comparison between the observed and simulated values of flood extent showed sufficiently good performance of the hydraulic model (CSI = 82%). However, the statistical analysis of the observed and simulated flood depths displayed a flood depth overestimation by the applied model, despite that the values of the used statistic indexes are acceptable (RMSE = 0.35 m, SD = 0.53, NSE = 0.56, PBIAS = 11.26%). The model overestimation of flood depth was attributed to the DEM low resolution and quality. Ground and aerial observations depicted the alluvial fan activation, the alternation of flow paths and the huge sediment transport. Human intervention in main streams, urban sprawl, wet AMC and sediment transport were among the main factors that contributed to the flash flood generation. This integrated approach revealed the necessity of the constant evaluation and validation of hydrological and hydraulic models in small ungauged Mediterranean watersheds and ephemeral streams. The use of UAVs in combination with ground observations and hydraulic simulation could significantly contribute to the enhanced understanding of flash flood mechanisms, in the direction of flood risk mitigation, improvement of the planning efficiency of flood prevent measures, flood hazard estimation, evolution of flood warning systems and floodplain geomorphology analysis.     

Links, comparisons and extensions of the geographically weighted regression model when used as a spatial predictor

In this study, we link and compare the geographically weighted regression (GWR) model with the kriging with an external drift (KED) model of geostatistics. This includes empirical work where models are performance tested with respect to prediction and prediction uncertainty accuracy. In basic forms, GWR and KED (specified with local neighbourhoods) both cater for nonstationary correlations (i.e. the process is heteroskedastic with respect to relationships between the variable of interest and its covariates) and as such, can predict more accurately than models that do not. Furthermore, on specification of an additional heteroskedastic term to the same models (now with respect to a process variance), locally-accurate measures of prediction uncertainty can result. These heteroskedastic extensions of GWR and KED can be preferred to basic constructions, whose measures of prediction uncertainty are only ever likely to be globally-accurate. We evaluate both basic and heteroskedastic GWR and KED models using a case study data set, where data relationships are known to vary across space. Here GWR performs well with respect to the more involved KED model and as such, GWR is considered a viable alternative to the more established model in this particular comparison. Our study adds to a growing body of empirical evidence that GWR can be a worthy predictor; complementing its more usual guise as an exploratory technique for investigating relationships in multivariate spatial data sets.     

Multi‐temporal scale changes of streamflow and sediment load in a loess hilly watershed of China

Global climate change and diverse human activities have resulted in distinct temporal–spatial variability of watershed hydrological regimes, especially in water‐limited areas. This study presented a comprehensive investigation of streamflow and sediment load changes on multi‐temporal scales (annual, flood season, monthly and daily scales) during 1952–2011 in the Yanhe watershed, Loess Plateau. The results indicated that the decreasing trend of precipitation and increasing trend of potential evapotranspiration and aridity index were not significant. Significant decreasing trends (p < 0.01) were detected for both the annual and flood season streamflow, sediment load, sediment concentration and sediment coefficient. The runoff coefficient exhibited a significantly negative trend (p < 0.01) on the flood season scale, whereas the decreasing trend on the annual scale was not significant. The streamflow and sediment load during July–August contributed 46.7% and 86.2% to the annual total, respectively. The maximum daily streamflow and sediment load had the median occurrence date of July 31, and they accounted for 9.7% and 29.2% of the annual total, respectively. All of these monthly and daily hydrological characteristics exhibited remarkable decreasing trends (p < 0.01). However, the contribution of the maximum daily streamflow to the annual total progressively decreased (?0.07% year^?1), while that of maximum daily sediment load increased over the last 60 years (0.08% year^?1). The transfer of sloping cropland for afforestation and construction of check‐dams represented the dominant causes of streamflow and sediment load reductions, which also made the sediment grain finer. Copyright © 2015 John Wiley & Sons, Ltd.     

Multi‐site evaluation of hydrology component of SWAT in the coastal plain of southwest Georgia

Simulation of watershed scale hydrologic and water quality processes is important for watershed assessments. Proper characterization of the accuracy of these simulations, particularly in cases with limited observed data, is critical. The Soil & Water Assessment Tool (SWAT) is frequently used for watershed scale simulation. The accuracy of the model was assessed by extrapolating calibration results from a well studied Coastal Plain watershed in Southwest Georgia, USA, to watersheds within the same geographic region without further calibration. SWAT was calibrated and validated on a 16.7‐km² subwatershed within the Little River Experimental Watershed by varying six model parameters. The optimized parameter set was then applied to a watershed of similar land use and soils, a smaller watershed with different land use and soils and three larger watersheds within the same drainage system without further calibration. Simulation results with percent bias (PB) ±15% ≤ PB < ±25% and Nash–Sutcliffe efficiency (NSE) 0.50 < NSE ≤ 0.65 were considered to be satisfactory, whereas those with PB < ±10% and 0.75 < NSE ≤ 1.00 were considered very good. With these criteria, simulation results for the five non‐calibration watersheds were satisfactory to very good. Differences across watersheds were attributed to differences in soils, land use, and surficial aquifer characteristics. These results indicate that SWAT can be a useful tool for predicting streamflow for ungauged watersheds with similar physical characteristics to the calibration watershed studied here and provide an indication of the accuracy of hydrologic simulations for ungauged watersheds. Copyright © 2012 John Wiley & Sons, Ltd.     

Reconstruction of a flash flood with large wood transport and its influence on hazard patterns in an ungauged mountain basin

The reconstruction of past flash floods in ungauged basins leads to a high level of uncertainty, which increases if other processes are involved such as the transport of large wood material. An important flash flood occurred in 1997 in Venero Claro (Central Spain), causing significant economic losses. The wood material clogged bridge sections, raising the water level upstream. The aim of this study was to reconstruct this event, analysing the influence of woody debris transport on the flood hazard pattern. Because the reach in question was affected by backwater effects due to bridge clogging, using only high water mark or palaeostage indicators may overestimate discharges, and so other methods are required to estimate peak flows. Therefore, the peak discharge was estimated (123 ± 18 m³ s^–1) using indirect methods, but one‐dimensional hydraulic simulation was also used to validate these indirect estimates through an iterative process (127 ± 33 m³ s^–1) and reconstruct the bridge obstruction to obtain the blockage ratio during the 1997 event (~48%) and the bridge clogging curves. Rainfall–Runoff modelling with stochastic simulation of different rainfall field configurations also helped to confirm that a peak discharge greater than 150 m³ s^–1 is very unlikely to occur and that the estimated discharge range is consistent with the estimated rainfall amount (233 ± 27 mm). It was observed that the backwater effect due to the obstruction (water level ~7 m) made the 1997 flood (~35‐year return period) equivalent to the 50‐year flood. This allowed the equivalent return period to be defined as the recurrence interval of an event of specified magnitude, which, where large woody debris is present, is equivalent in water depth and extent of flooded area to a more extreme event of greater magnitude. These results highlight the need to include obstruction phenomena in flood hazard analysis. Copyright © 2012 John Wiley & Sons, Ltd.     

Comparison of five snow water equivalent estimation methods across categories

Snow water equivalent (SWE) is an important indicator used in hydrology, water resources, and climate change impact. There are various methods of estimating SWE (falling in 3 categories: indirect sensors, empirical models, and process‐based models), but few studies that provide comparison across these different categories to help users make decisions on monitoring site design or method selection. Five SWE estimation methods were compared against manual snow course data collected over 2 years (2015–2016) from the Dorset Environmental Science Centre, including the gamma‐radiation‐based CS725 sensor, 3 empirical estimation models (Sexstone snow density model, McCreight & Small snow density model, and a meteorology‐based model), and the University of British Columbia Watershed Model snow energy‐balance model. Snow depth, density, and SWE were measured at the Dorset Environmental Science Centre weather station in south‐central Ontario, on a daily basis over 6 winters from 2011 to 2016. The 2 snow density‐based models, requiring daily snow depth as input, gave the best performance (R² of .92 and .92 for McCreight & Small and Sexstone models, respectively). The CS725 sensor that receives radiation coming from soil penetrating the snowpack provided the same performance (R² = .92), proving that the sensor is an applicable method, although it is expensive. The meteorology‐based empirical model, requiring daily climate data including temperature, precipitation and solar radiation, gave the poorest performance (R² = .77). The energy‐balance‐based University of British Columbia Watershed Model snow module, only requiring climate data, worked better than the empirical meteorology‐based model (R² = .9) but performed worse than the density models or CS725 sensor. Given differences in application objectives, site conditions, and budget, this comparison across SWE estimation methods may help users choose a suitable method. For ongoing and new monitoring sites, installation of a CS725 sensor coupled with intermittent manual snow course measurements (e.g., weekly) is recommended for further SWE method estimation testing and development of a snow density model.     

INTERNATIONAL HYDROLOGICAL DECADE / WORKING GROUP ON “THE INFLUENCE OF MAN ON THE HYDROLOGICAL CYCLE”

Abstract

Floods, as extreme hydrological phenomena, can be described by more than one correlated characteristic, such as peak, volume and duration. These characteristics should be jointly considered since they are generally not independent. For an ungauged site, univariate regional flood frequency analysis (FA) provides a limited assessment of flood events. A recent study proposed a procedure for regional FA in a multivariate framework. This procedure represents a multivariate version of the index-flood model and is based on copulas and multivariate quantiles. The performance of the proposed procedure was evaluated by simulation. However, the model was not tested on a real-world case study data. In the present paper, practical aspects are investigated jointly for flood peak (Q) and volume (V) of a dataset from the Côte-Nord region in the province of Quebec, Canada. The application of the proposed procedure requires the identification of the appropriate marginal distribution, the estimation of the index flood and the selection of an appropriate copula. The results of the case study show that the regional bivariate FA procedure performed well. This performance depends strongly on the performance of the two univariate models and, more specifically, the univariate model of Q. The results show also the impact of the homogeneity of the region on the performance of the univariate and bivariate models.

Editor D. Koutsoyiannis     

Prediction of rainfall‐induced shallow landslides in the Loess Plateau,Yan'an,China, using the TRIGRS model

In this work, a transient rainfall infiltration and grid‐based regional slope‐stability model (TRIGRS) was implemented in a case study of Yan'an City, Northwest China. In this area, widespread shallow landslides were triggered by the 12 July 2013 exceptional rainstorm event. A high‐resolution DEM, soil parameters from in‐situ and laboratory measurements, water table depths, the maximum depth of precipitation infiltration and rain‐gauge‐corrected precipitation of the event, were used as inputs in the TRIGRS model. Shallow landslides triggered on the same day were used to evaluate the modeling results. The summarized results are as follows: (i) The characteristics and distribution of thirty‐five shallow landslides triggered by the 12 July 2013 rainfall event were identified in the study area and all were classified as shallow landslides with the maximum depth, area and volume less than 3 m, 200 m² and 1000 m³, respectively, (ii) Four intermediate factor of safety (F_S) maps were generated using the TRIGRS model to represent the scenarios 6, 12, 18 and 24 hours after the storm event. The area with F_S < 1 increased with the rainfall duration. The percentage of the area with F_S < 1 was 0.2%, 3.3%, 3.8% and 5.1% for the four stages, respectively. Twenty‐four hours after the rainstorm, TRIGRS predicted that 1255 grid cells failed, which is consistent with the field data. (iii) TRIGRS generated more satisfactory results at a given precipitation threshold than SINMAP, which is ideal for landslide hazard zoning for land‐use planning at the regional scale. Comparison results showed that TRIGRS is more useful for landslide prediction for a certain precipitation threshold, also in the regional scale. (iv) Analysis of the responses of loess slope prone to slope failure after different precipitation scenarios revealed that loess slopes are particularly sensitive to extended periods of heavy precipitation. Copyright © 2016 John Wiley & Sons, Ltd.     

Flood frequency analysis of Turkey using L‐moments method

The index flood procedure coupled with the L‐moments method is applied to the annual flood peaks data taken at all stream‐gauging stations in Turkey having at least 15‐year‐long records. First, screening of the data is done based on the discordancy measure (D_i) in terms of the L‐moments. Homogeneity of the total geographical area of Turkey is tested using the L‐moments based heterogeneity measure, H, computed on 500 simulations generated using the four parameter Kappa distribution. The L‐moments analysis of the recorded annual flood peaks data at 543 gauged sites indicates that Turkey as a whole is hydrologically heterogeneous, and 45 of 543 gauged sites are discordant which are discarded from further analyses. The catchment areas of these 543 sites vary from 9·9 to 75121 km² and their mean annual peak floods vary from 1·72 to 3739·5 m³ s^?1. The probability distributions used in the analyses, whose parameters are computed by the L‐moments method are the general extreme values (GEV), generalized logistic (GLO), generalized normal (GNO), Pearson type III (PE3), generalized Pareto (GPA), and five‐parameter Wakeby (WAK). Based on the L‐moment ratio diagrams and the |Z_dist|‐statistic criteria, the GEV distribution is identified as the robust distribution for the study area (498 gauged sites). Hence, for estimation of flood magnitudes of various return periods in Turkey, a regional flood frequency relationship is developed using the GEV distribution. Next, the quantiles computed at all of 543 gauged sites by the GEV and the Wakeby distributions are compared with the observed values of the same probability based on two criteria, mean absolute relative error and determination coefficient. Results of these comparisons indicate that both distributions of GEV and Wakeby, whose parameters are computed by the L‐moments method, are adequate in predicting quantile estimates. Copyright © 2011 John Wiley & Sons, Ltd.     

基于随机森林模型的GRACE数据3种空间降尺度对比

陆地水储量是赋存在陆地上各种形式水的综合体现,研究其时空变化对认识区域水循环过程和水资源调控等具有重要意义。然而现有陆地水储量变化数据实际分辨率较低,限制了其在中小流域或地区中的应用。针对这一问题,本文基于GRACE重力卫星和其后续卫星GRACE-FO反演的陆地水储量变化数据,首先采用随机森林模型,分别基于格点、区域(流域)和区域(全国)3种空间降尺度思路将GRACE数据降尺度至0.25°×0.25°,后结合GLDAS模型数据,基于水量平衡原理计算得到地下水储量变化数据,最后基于降尺度模型模拟效果和实测地下水位数据评估3种降尺度思路在全国的适用性。结果表明:随机森林模型能够较好地模拟驱动数据(降水、气温、植被条件指数和土壤水储量)与GRACE数据的统计关系,验证期格点降尺度思路的平均相关系数总体在0.6左右,区域降尺度思路的平均纳什效率系数、相关系数和均方根误差分别>0.5、>0.75和<6.6 cm,3种空间降尺度思路的模拟精度均满足基本要求;2003—2021年间,GRACE数据、格点降尺度、区域降尺度(流域)和区域降尺度(全国)得到的我国陆地水储量亏缺量分别约为...     

Classification of yearly extreme precipitation events and associated flood risk in the Yangtze-Huaihe River Valley

Fifty cases of regional yearly extreme precipitation events(RYEPEs)were identified over the Yangtze-Huaihe River Valley(YHRV)during 1979–2016 applying the statistical percentile method.There were five types of RYEPEs,namely Yangtze Meiyu(YM-RYEPE),Huaihe Meiyu(HM-RYEPE),southwest-northeast-oriented Meiyu(SWNE-RYEPE)and typhoon I and II(TC-RYEPE)types of RYEPEs.Potential vorticity diagnosis showed that propagation trajectories of the RYEPEs along the Western Pacific Subtropical High and its steering flow were concentrated over the southern YHRV.As a result,the strongest and most frequently RYEPEs events,about 16–21 cases with average rainfall above 100 mm,occurred in the southern YHRV,particularly in the Nanjing metropolitan area.There have been 14 cases of flood-inducing RYEPEs since 1979,with the submerged area exceeding 120 km~2as simulated by the Flood Area hydraulic model,comprising six HM-RYEPEs,five YM-RYEPEs,two TC-RYEPEs,and one SWNE-RYEPE.The combination of evolving RYEPEs and rapid expansion of urban agglomeration is most likely to change the flood risk distribution over the Nanjing metropolitan area in the future.In the RCP6.0(RCP8.5)scenario,the built-up area increases at a rate of about 10.41 km~2(10 yr)~(-1)(24.67 km~2(10 yr)~(-1))from 2010 to 2100,and the area of high flood risk correspondingly increases from 3.86 km~2(3.86 km~2)to 9.00 km~2(13.51 km~2).Areas of high flood risk are mainly located at Chishan Lake in Jurong,Lukou International Airport in Nanjing,Dongshan in Jiangning District,Lishui District and other low-lying areas.The accurate simulation of flood scenarios can help reduce losses due to torrential flooding and improve early warnings,evacuation planning and risk analysis.More attention should be paid to the projected high flood risk because of the concentrated population,industrial zones and social wealth throughout the Nanjing metropolitan area.     

Hydrological and meteorological extremes derived from taxation records: the estates of Brtnice,Třebíč and Velké Meziříčí, 1706–1849

Abstract

This paper addresses the hydrological and meteorological extremes that may be deduced from the taxation records of the estates of Brtnice, T?ebí? and Velké Mezi?í?í, all in the Moravian-Bohemian Highlands of the Czech Republic, for the years 1706–1849. At that time, damage to agricultural crops constituted grounds for tax remission for individual farmers and landowners. Where it survives, the relevant administrative documentation generally includes a statement from the applicant, a report by the official commission tasked with checking the contents of it, and any decisions made by taxation authorities at regional and “land” level (the Jihlava regional office and the Moravian Land Administration (“Gubernium”) respectively). Data extracted may include the type of event, dating, places of occurrence and damage done. The chronology of hydrological and meteorological extremes (torrential rain, flash flood, flood, hailstorm, lightning, frost) covers the period 1706–1849, but only four events are evident before 1748 and there is a gap in records between 1757 and 1789. Extremes are analysed from a spatio-temporal point of view. A total of 97 extreme events (171 extremes of particular type) were identified for the region studied. Torrential rain, hailstorm and flash flood were the major devastating phenomena, and occurred mainly from May to August. Torrential rain and hailstorm are clearly attributable to thunderstorms with very intense convection. Five outstanding events and their impacts upon individual farmers are described in detail. The results are discussed with respect to uncertainties in the basic data and in the context of the Czech Lands, because only some of the extremes disclosed are known and confirmed by other documentary data.

Editor Z.W. Kundzewicz

Citation Dolák, L., Brázdil, R., and Valá?ek, H., 2013. Hydrological and meteorological extremes derived from taxation records: the estates of Brtnice, T?ebí? and Velké Mezi?í?í, 1706–1849. Hydrological Sciences Journal, 58 (8), 1620–1634.