Critical rainfall statistics for predicting watershed flood responses: rethinking the design storm concept

Recent advances have been made to modernize estimates of probable precipitation scenarios; however, researchers and engineers often continue to assume that rainfall events can be described by a small set of event statistics, typically average intensity and event duration. Given the easy availability of precipitation data and advances in desk‐top computational tools, we suggest that it is time to rethink the ‘design storm’ concept. Design storms should include more holistic characteristics of flood‐inducing rain events, which, in addition to describing specific hydrologic responses, may also be watershed or regionally specific. We present a sensitivity analysis of nine precipitation event statistics from observed precipitation events within a 60‐year record for Tompkins County, NY, USA. We perform a two‐sample Kolmogorov–Smirnov (KS) test to objectively identify precipitation event statistics of importance for two related hydrologic responses: (1) peak outflow from the Six Mile Creek watershed and (2) peak depth within the reservoir behind the Six Mile Creek Dam. We identify the total precipitation depth, peak hourly intensity, average intensity, event duration, interevent duration, and several statistics defining the temporal distribution of precipitation events to be important rainfall statistics to consider for predicting the watershed flood responses. We found that the two hydrologic responses had different sets of statistically significant parameters. We demonstrate through a stochastic precipitation generation analysis the effects of starting from a constrained parameter set (intensity and duration) when predicting hydrologic responses as opposed to utilizing an expanded suite of rainfall statistics. In particular, we note that the reduced precipitation parameter set may underestimate the probability of high stream flows and therefore underestimate flood hazard. Copyright © 2016 John Wiley & Sons, Ltd.     

Fingerprinting the 8.2 ka event climate response in a coupled climate model

Using results from coupled climate model simulations of the 8.2 ka climate event that produced a cold period over Greenland in agreement with the reconstructed cooling from ice cores, we investigate the typical pattern of climate anomalies (fingerprint) to provide a framework for the interpretation of global proxy data for the 8.2 ka climate event. For this purpose we developed an analysis method that isolates the forced temperature response and provides information on spatial variations in magnitude, timing and duration that characterise the detectable climate event in proxy archives. Our analysis shows that delays in the temperature response to the freshwater forcing are present, mostly in the order of decades (30 a over central Greenland). The North Atlantic Ocean initially cools in response to the freshwater perturbation, followed in certain parts by a warm response. This delay, occurring more than 200 a after the freshwater pulse, hints at an overshoot in the recovery from the freshwater perturbation. The South Atlantic and the Southern Ocean show a warm response reflecting the bipolar seesaw effect. The duration of the simulated event varies for different areas, and the highest probability of recording the event in proxy archives is in the North Atlantic Ocean area north of 40° N. Our results may facilitate the interpretation of proxy archives recording the 8.2 ka event, as they show that timing and duration cannot be assumed to correspond with the timing and duration of the event as recorded in Greenland ice cores. Copyright © 2011 John Wiley & Sons, Ltd.     

Severe soil erosion during a 3‐day exceptional rainfall event: combining modelling and field data for a fallow cereal field

Exceptional rainfall events cause significant losses of soil, although few studies have addressed the validation of model predictions at field scale during severe erosive episodes. In this study, we evaluate the predictive ability of the enhanced Soil Erosion and Redistribution Tool (SERT‐2014) model for mapping and quantifying soil erosion during the exceptional rainfall event (~235 mm) that affected the Central Spanish Pyrenees in October 2012. The capacity of the simulation model is evaluated in a fallow cereal field (1.9 ha) at a high spatial scale (1 × 1 m). Validation was performed with field‐quantified rates of soil loss in the rills and ephemeral gullies and also with a detailed map of soil redistribution. The SERT‐2014 model was run for the six rainfall sub‐events that made up the exceptional event, simulating the different hydrological responses of soils with maximum runoff depths ranging between 40 and 1017 mm. Predicted average and maximum soil erosion was 11 and 117 Mg ha^?1 event^?1, respectively. Total soil loss and sediment yield to the La Reina gully amounted to 16.3 and 9.0 Mg event^?1. These rates are in agreement with field estimations of soil loss of 20.0 Mg event^?1. Most soil loss (86%) occurred during the first sub‐event. Although soil accumulation was overestimated in the first sub‐event because of the large amount of detached soil, the enhanced SERT‐2014 model successfully predicted the different spatial patterns and values of soil redistribution for each sub‐event. Further research should focus on stream transport capacity. Copyright © 2014 John Wiley & Sons, Ltd.     

Late Quaternary sea‐level change and evolution of Belfast Lough,Northern Ireland: new offshore evidence and implications for sea‐level reconstruction

The interplay of eustatic and isostatic factors causes complex relative sea‐level (RSL) histories, particularly in paraglacial settings. In this context the past record of RSL is important in understanding ice‐sheet history, earth rheology and resulting glacio‐isostatic adjustment. Field data to develop sea‐level reconstructions are often limited to shallow depths and uncertainty exists as to the veracity of modelled sea‐level curves. We use seismic stratigraphy, 39 vibrocores and 26 radiocarbon dates to investigate the deglacial history of Belfast Lough, Northern Ireland, and reconstruct past RSL. A typical sequence of till, glacimarine and Holocene sediments is preserved. Two sea‐level lowstands (both max. ?40 m) are recorded at c. 13.5 and 11.5k cal a bp . Each is followed by a rapid transgression and subsequent periods of RSL stability. The first transgression coincides temporally with a late stage of Meltwater Pulse 1a and the RSL stability occurred between c. 13.0 and c. 12.2k cal a bp (Younger Dryas). The second still/slowstand occurred between c. 10.3 and c. 11.5k cal a bp . Our data provide constraints on the direction and timing of RSL change during deglaciation. Application of the Depth of Closure concept adds an error term to sea‐level reconstructions based on seismic stratigraphic reconstructions.     

Stormflow concentration–discharge dynamics of suspended sediment and dissolved phosphorus in an agricultural watershed

Concentration–discharge (C-Q) relationships are an effective tool for identifying watershed biogeochemical source and transport dynamics over short and long timescales. We examined stormflow C-Q, hysteresis, and flushing patterns of total suspended sediment (TSS) and soluble reactive phosphorus (SRP) in two stream reaches of a severely impaired agricultural watershed in northeastern Wisconsin, USA. The upper watershed reach—draining a relatively flat, row crop-dominated contributing area—showed predominantly anti-clockwise TSS hysteresis during storms, suggesting that particulate materials were mobilized more from distal upland sources than near- and in-channel areas. In contrast, the incised lower watershed reach produced strong TSS flushing responses on the rising limb of storm hydrographs and clockwise hysteresis, signalling rapid mobilization of near- and in-channel materials with increasing event flows. C-Q relationships for SRP showed complex patterns in both the upper and lower reaches, demonstrating largely non-linear chemodynamic C-Q behaviour during events. As with TSS, anti-clockwise SRP hysteresis in the upper reach suggested a delay in the hydrologic connectivity between SRP sources and the stream, with highly variable SRP concentrations during some events. A broad range of clockwise, anti-clockwise, and complex SRP hysteresis patterns occurred in the lower watershed, possibly influenced by in-channel legacy P stores and connection to tile drainage networks in the lower watershed area. Total suspended sediment and SRP responses were also strongly related to precipitation event characteristics including antecedent precipitation, recovery period, and precipitation intensity, highlighting the complexity of stormflow sediment and phosphorus responses in this severely impaired agricultural stream.     

The Upstream “Strong Signals” of the Water Vapor Transport over the Tibetan Plateau during a Heavy Rainfall Event in the Yangtze River Basin

A heavy rainfall event that occurred over the middle and lower reaches of the Yangtze River Basin(YRB) during July11–13 2000 is explored in this study. The potential/stream function is used to analyze the upstream "strong signals" of the water vapor transport in the Tibetan Plateau(TP). The studied time period covers from 2000 LST 5 July to 2000 LST 15 July(temporal resolution: 6 hours). By analyzing the three-dimensional structure of the water vapor flux, vorticity and divergence prior to and during the heavy rainfall event, the upstream "strong signals" related to this heavy rainfall event are revealed. A strong correlation exists between the heavy rainfall event in the YRB and the convective clouds over the TP. The "convergence zone" of the water vapor transport is also identified, based on correlation analysis of the water vapor flux two days and one day prior to, and on the day of, the heavy rainfall. And this "convergence zone" coincides with the migration of the maximum rainfall over the YRB. This specific coupled structure actually plays a key role in generating heavy rainfall over the YRB. The eastward movement of the coupled system with a divergence/convergence center of the potential function at the upper/lower level resembles the spatiotemporal evolution of the heavy rainfall event over the YRB. These upstream "strong signals" are clearly traced in this study through analyzing the three-dimensional structure of the potential/stream function of upstream water vapor transport.     

敦密断裂带白垩纪两期重要的变形事件

本文报道了敦密断裂带糜棱岩中黑云母~(40)Ar/~(39)Ar定年结果和大规模走滑-逆冲断裂的几何学、运动学特征及其形成时代,以便揭示断裂带两期变形事件的构造属性。黑龙江省密山市花岗质糜棱岩中黑云母~(40)Ar/~(39)Ar加权平均年龄为132.2±1.2Ma,它是敦密断裂带经历伸展事件的冷却年龄,也是东北亚大陆边缘在早白垩世Hauterivian期-Albian期发生强烈区域伸展作用的产物。密山市至辽宁省清原县系列大型走滑-逆冲断层和断层相关褶皱揭示出在晚白垩世晚期-末期发生右旋走滑-逆冲事件,该事件规模大,影响范围广,导致整个断裂带遭受到强烈改造,形成对冲式断裂系统。将研究区走滑-逆冲断裂与山东省郯庐断裂带中段挤压构造对比,认为郯庐断裂带北段和中段在晚白垩世末期都发生了强烈的走滑-逆冲事件,它们具有相同的构造特征和构造属性。     

Oil dispersal modelling: reanalysis of the Rena oil spill using open-source modelling tools

ABSTRACT

Oil spill forecast modelling is typically used immediately after a spill to predict oil dispersal and promote mobilisation of more effective response operations. The aim of this work was to map oil dispersal after the grounding of the MV Rena on Astrolabe Reef and to verify the results against observations. Model predictions were broadly consistent with observed distribution of oil contamination. However, some hot spots of oil accumulation, likely due to surf-zone and rip current circulation, were not well represented. Additionally, the model was run with 81 differing wind conditions to show that the events occurring during the grounding represented the typical likely behaviour of an oil spill on Astrolabe Reef. Oil dispersal was highly dependent on prevailing wind patterns; more accurate prediction would require better observations of local wind patterns. However, comparison of predictions with observations indicated that the GNOME model was an effective low-cost approach.     

Occurrence and Chemical Compositions of Amphiboles in Altered Dioritic Rocks of Laiwu Skarn‐Type Iron Deposit in West Shandong Area,China

The Laiwu Fe deposit is the largest skarn‐type deposit in West Shandong in China, with an estimated reserve of 290 Mt of iron ore. Here, we investigate the occurrence and chemical composition of amphiboles in altered dioritic rocks of this deposit. Three generations of amphibole group minerals were identified in these rocks. The first type (Type 1) is a hornblende with a generally euhedral crystal shape, 150–500 μm in size, and dominantly occurs in the host monzonite. The second type (Type 2) of hornblende occurs as fine grains (<100 μm) or concentric rims, less than 50 μm in thickness, around the Type 1 hornblende. The third type (Type 3) of amphibole is anhedral actinolite occurring along the boundary or the cracks of the Type 1 or 2 hornblende grains. The Type 1 hornblende has higher FeO (12.5–15.6 wt.%) and lower MgO contents (11.2–14.6 wt.%) than the Type 2 hornblende, which has an FeO content ranging from 8.8 to 10.5 wt.% and MgO content ranging from 14.2 to 17.1 wt.%. The Type 3 actinolite shows the lowest FeO (4.6–6.4 wt.%) and highest MgO (19.8–21.2 wt.%) contents. We infer that the Type 2 hornblende formed from Fe‐rich hydrothermal fluids released during rapid upwelling of the crystallizing magma. The fluids were relatively reduced and enriched in Fe. The Type 3 actinolite formed by coupled dissolution and reprecipitation, and its formation is one of the features denoting Fe‐enrichment events. Our study of these types of amphibole provides insights into the ore formation process.