Non-stationary frequency analysis of extreme precipitation in South Korea using peaks-over-threshold and annual maxima

The conventional approach to the frequency analysis of extreme precipitation is complicated by non-stationarity resulting from climate variability and change. This study utilized a non-stationary frequency analysis to better understand the time-varying behavior of short-duration (1-, 6-, 12-, and 24-h) precipitation extremes at 65 weather stations scattered across South Korea. Trends in precipitation extremes were diagnosed with respect to both annual maximum precipitation (AMP) and peaks-over-threshold (POT) extremes. Non-stationary generalized extreme value (GEV) and generalized Pareto distribution (GPD) models with model parameters made a linear function of time were applied to AMP and POT respectively. Trends detected using the Mann–Kendall test revealed that the stations showing an increasing trend in AMP extremes were concentrated in the mountainous areas (the northeast and southwest regions) of South Korea. Trend tests on POT extremes provided fairly different results, with a significantly reduced number of stations showing an increasing trend and with some stations showing a decreasing trend. For most of stations showing a statistically significant trend, non-stationary GEV and GPD models significantly outperformed their stationary counterparts, particularly for precipitation extremes with shorter durations. Due to a significant-increasing trend in the POT frequency found at a considerable number of stations (about 10 stations for each rainfall duration), the performance of modeling POT extremes was further improved with a non-homogeneous Poisson model. The large differences in design storm estimates between stationary and non-stationary models (design storm estimates from stationary models were significantly lower than the estimates of non-stationary models) demonstrated the challenges in relying on the stationary assumption when planning the design and management of water facilities. This study also highlighted the need of caution when quantifying design storms from POT and AMP extremes by showing a large discrepancy between the estimates from those two approaches.     

The role of antecedent groundwater heads in controlling transient aquifer storage and flood peak attenuation in karst watersheds

Transient storage of floodwaters in aquifers is known to attenuate peak flows in rivers and drive subsurface dissolution. Transient aquifer storage could be enhanced in watersheds overlying karst aquifers where caves facilitate surface and groundwater exchange. Few studies, however, have examined controls on, or magnitudes of, transient aquifer storage or flood peak attenuation in karstic watersheds. Here we evaluate flood peak attenuation with multiple linear regression analyses of 10 years of river and groundwater data from the Suwannee River, which flows over the karstic upper Floridan aquifer in north-central Florida and experiences frequent flooding. Regressions show antecedent river stage exerts the dominant control on magnitudes of transient aquifer storage, with recharge and time to peak having secondary controls. Specifically, low antecedent stages result in larger magnitudes of transient aquifer storage and thus greater flood attenuation than conditions of elevated antecedent stage. These findings suggest subsurface weathering, including cave formation and enlargement, caused by transient aquifer storage could occur on a more frequent basis in aquifers where groundwater table elevation is lowered due to anthropogenic or climatic influences. Our work also shows that measures of groundwater table elevation prior to an event could be used to improve predictive flood models. © 2018 John Wiley & Sons, Ltd.     

Tide-induced variability in beach surface moisture: Observations and modelling

The moisture content w_s of a beach surface strongly controls the availability of sand for aeolian transport. Our predictive capability of the spatiotemporal variability in w_s, which depends to a large extent on water table depth, is, however, limited. Here we show that water table fluctuations and surface moisture content observed during a 10-day period on a medium-grained (365μm) planar (1:30) beach can be predicted well with the nonlinear Boussinesq equation extended to include run-up infiltration and a soil–water retention curve under the assumption of hydrostatic equilibrium. On the intertidal part of the beach the water table is observed and predicted to continuously fall from the moment the beach surface emerges from the falling tide to just before it is submerged by the incoming tide. We find that on the lower 30% of the intertidal beach the water table remains within 0.1–0.2 m from the surface and that the sand is always saturated (w_s≈20%, by mass). Higher up on the intertidal beach, the surface can dry to about 5% when the water table has fallen to 0.4–0.5 m beneath the surface. Above the high-tide level the water table is always too deep (>0.5 m) to affect surface moisture and, without precipitation, the sand is dry (w_s < 5 − 8%). Because the water table depth on the emerged part of the intertidal beach increases with time irrespective of whether the (ocean) tide falls or rises, we find no need to include hysteresis (wetting and drying) effects in the surface-moisture modelling. Model simulations suggest that at the present planar beach only the part well above mean sea level can dry sufficiently (w_s < 10%) for sand to become available for aeolian transport. ©2018 The Authors. Earth Surface Processes and Landforms published by John Wiley & Sons Ltd.     

CO2 saturation estimates at Sleipner (North Sea) from seismic tomography and rock physics inversion

CO₂ saturations are estimated at Sleipner using a two-step imaging workflow. The workflow combines seismic tomography (full-waveform inversion) and rock physics inversion and is applied to a two-dimensional seismic line located near the injection point at Sleipner. We use baseline data (1994 vintage, before CO₂ injection) and monitor data that was acquired after 12 years of CO₂ injection (2008 vintage). P-wave velocity models are generated using the Full waveform inversion technology and then, we invert selected rock physics parameters using an rock physics inversion methodology. Full waveform inversion provides high-resolution P-wave velocity models both for baseline and monitor data. The physical relations between rock physics properties and acoustic wave velocities in the Utsira unconsolidated sandstone (reservoir formation) are defined using a dynamic rock physics model based on well-known Biot–Gassmann theories. For data prior to injection, rock frame properties (porosity, bulk and shear dry moduli) are estimated using rock physics inversion that allows deriving physically consistent properties with related uncertainty. We show that the uncertainty related to limited input data (only P-wave velocity) is not an issue because the mean values of parameters are correct. These rock frame properties are then used as a priori constraint in the monitor case. For monitor data, the Full waveform inversion results show nicely resolved thin layers of CO₂–brine saturated sandstones under intra-reservoir shale layers. The CO₂ saturation estimation is carried out by plugging an effective fluid phase in the rock physics model. Calculating the effective fluid bulk modulus of the brine–CO₂ mixture (using Brie equation in our study) is shown to be the key factor to link P-wave velocity to CO₂ saturation. The inversion tests are done with several values of Brie/patchiness exponent and show that the CO₂ saturation estimates are varying between 0.30 and 0.90 depending on the rock physics model and the location in the reservoir. The uncertainty in CO₂ saturation estimation is usually lower than 0.20. When the patchiness exponent is considered as unknown, the inversion is less constrained and we end up with values of exponent varying between 5 and 20 and up to 33 in specific reservoir areas. These estimations tend to show that the CO₂–brine mixing is between uniform and patchy mixing and variable throughout the reservoir.     

Leaky Aquifer Models to Simulate the Well Flow in Fractured Aquifers with Linear Interporosity Flow

Following Hemker and Maas (1987) the models of two or three leaky aquifers are applied to simulate the flow to vertical wells operating in the fractured or dual porosity aquifers. The software WellTest (WT) (Székely 2015) is used for calculating the drawdown and discharge rate variation. The comparative analysis with the independent analytical solutions by Boulton and Streltsova-Adams (1978), Warren and Root (1963), Kazemi et al. (1969) concluded with acceptable agreement between the WT simulation and the alternate calculation methods. The selected field tests have been conducted in fractured limestone aquifers. The pumping test west of Copenhagen shows an example of fractured aquifer with considerable negative skin effect at the well face. The flowing well Wafra W1 in Kuwait operates in the two-zone aquifer exhibiting sufficient vertical recharge via leakage beyond a circular domain of estimated radius of 2460 m.     

Evaluation of acid drainage in ore deposits

Water–rock–gas system is simulated using a technique based on the fundamental principles of modern hydrogeochemistry regarding the formation parameters of natural-water chemistry to show that groundwater pollution hazard at a gold-bearing placer is acute for any site of water exchange changing from passive to active. Under passive water exchange, the products of sulfide oxidation are accumulated; whereas under active water exchange, the oxidation products are leached out to form acid drainage flows. The treatment of acid drainage water and the neutralization of acid solutions are more effective under reduction or weakly alkaline conditions.     

Long-term Geomechanical Observation in the Jeroným Mine

Geomechanical measurements have been carried out in the medieval Jeroným Mine since 2001. At first quarterly observation and later also continual monitoring by distributed measurement network were used. The network includes a number of different sensors that have been installed in selected places of underground spaces. There have not been any manifestations documenting dangerous changes of stability during the period of monitoring (2001–2015). The results of accomplished measurements from all convergence profiles confirm the total stability of all underground spaces. However, a number of places have been documented as potentially unstable. Measured values of selected parameters are presented in the paper, namely, movements along fractures, changes of stress tensor and fluctuations in the level of mine water table.     

Channel response to an extreme flood and sediment pulse in a mixed bedrock and gravel‐bed river

We exploit a natural experiment caused by an extreme flood (~500 year recurrence interval) and sediment pulse derived from more than 2500 concurrent landslides to explore the influence of valley‐scale geomorphic controls on sediment slug evolution and the impact of sediment pulse passage and slug deposition and dispersion on channel stability and channel form. Sediment slug movement is a crucial process that shapes gravel‐bed rivers and alluvial valleys and is an important mechanism of downstream bed material transport. Further, increased bed material transport rates during slug deposition can trigger channel responses including increases in lateral mobility, channel width, and alluvial bar dominance. Pre‐ and post‐flood LiDAR and aerial photographs bracketing the 2007 flood on the Chehalis River in south‐western Washington State, USA, document the channel response with high spatial and temporal definition. The sediment slug behaved as a Gilbert Wave, with both channel aggradation and sequestration of large volumes of material in floodplains of headwaters' reaches and reaches where confined valleys enter into broad alluvial valleys. Differences between the valley form of two separate sub‐basins impacted by the pulse highlight the important role channel and channel‐floodplain connectivity play in governing downstream movement of sediment slug material. Finally, channel response to the extreme flood and sediment pulse illustrate the connection between bed material transport and channel form. Specifically, the channel widened, lateral channel mobility increased, and the proportion of the active channel covered by bars increased in all reaches in the study area. The response scaled tightly with the relative amount of bed material sediment transport through individual reaches, indicating that the amount of morphological change caused by the flood was conditioned by the simultaneous introduction of a sediment pulse to the channel network. Copyright © 2015 John Wiley & Sons, Ltd.