Modelling sediment deposition and phosphorus retention in a river floodplain

Phosphorus (P) is one of the major limiting nutrient in many freshwater ecosystems. During the last decade, attention has been focused on the fluxes of suspended sediment and particulate P through freshwater drainage systems because of severe eutrophication effects in aquatic ecosystems. Hence, the analysis and prediction of phosphorus and sediment dynamics constitute an important element for ecological conservation and restoration of freshwater ecosystems. In that sense, the development of a suitable prediction model is justified, and the present work is devoted to the validation and application of a predictive soluble reactive phosphorus (SRP) uptake and sedimentation models, to a real riparian system of the middle Ebro river floodplain. Both models are coupled to a fully distributed two‐dimensional shallow‐water flow numerical model. The SRP uptake model is validated using data from three field experiments. The model predictions show a good accuracy for SRP concentration, where the linear regressions between measured and calculated values of the three experiments were significant (r² ≥ 0.62; p ≤ 0.05), and a Nash–Sutcliffe coefficient (E) that ranged from 0.54 to 0.62. The sedimentation model is validated using field data collected during two real flooding events within the same river reach. The comparison between calculated and measured sediment depositions showed a significant linear regression (p ≤ 0.05; r² = 0.97) and an E that ranged from 0.63 to 0.78. Subsequently, the complete model that includes flow dynamics, solute transport, SRP uptake and sedimentation is used to simulate and analyse floodplain sediment deposition, river nutrient contribution and SRP uptake. According to this analysis, the main SRP uptake process appears to be the sediment sorption. The analysis also reveals the presence of a lateral gradient of hydrological connectivity that decreases with distance from the river and controls the river matter contribution to the floodplain. Copyright © 2014 John Wiley & Sons, Ltd.     

Tropical precipitation anomalies and d‐excess evolution during El Niño 2014‐16

The last 2014‐16 El Niño event was among the three strongest episodes on record. El Niño considerably changes annual and seasonal precipitation across the tropics. Here, we present a unique stable isotope data set of daily precipitation collected in Costa Rica prior to, during, and after El Niño 2014‐16, in combination with Lagrangian moisture source and precipitation anomaly diagnostics. δ²H composition ranged from ‐129.4 to +18.1 (‰) while δ¹⁸O ranged from ‐17.3 to +1.0 (‰). No significant difference was observed among δ¹⁸O (P=0.186) and δ²H (P=0.664) mean annual compositions. However, mean annual d‐excess showed a significant decreasing trend (from +13.3 to +8.7 ‰) (P<0.001) with values ranging from +26.6 to ‐13.9 ‰ prior to and during the El Niño evolution. The latter decrease in d‐excess can be partly explained by an enhanced moisture flux convergence across the southeastern Caribbean Sea coupled with moisture transport from northern South America by means of an increased Caribbean Low Level Jet regime. During 2014‐15, precipitation deficit across the Pacific domain averaged 46% resulting in a very severe drought; while a 94% precipitation surplus was observed in the Caribbean domain. Understanding these regional moisture transport mechanisms during a strong El Niño event may contribute to a) better understanding of precipitation anomalies in the tropics and b) re‐evaluate past stable isotope interpretations of ENSO events in paleoclimatic archives within the Central America region.     

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.     

In vitro screening for endocrine disruptive activity in selected South African harbours and river mouths

Various waterborne anthropogenic contaminants disrupt the endocrine systems of wildlife and humans, targeting reproductive pathways, among others. Very little is known, however, regarding the occurrence of endocrine disruptive activity in South African freshwater ecosystems, and coastal ecosystems have not been studied in this regard. In a first attempt to investigate endocrine disruptive activity in South African coastal waters, surface water samples collected from harbours, river mouths and estuaries in three metropolitan municipalities, eThekwini (which includes Durban), Nelson Mandela (specifically Port Elizabeth Harbour) and City of Cape Town, were screened for (anti) oestrogenicity and (anti)androgenicity using recombinant yeast bioassays. Moreover, levels of the female hormone 17β-(o)estradiol (E2) were determined in all samples, as well as a selection of hydrocarbons in the eThekwini samples. A high proportion of samples collected from eThekwini were oestrogenic, whereas none from Port Elizabeth Harbour and only a single river mouth sampled in the City of Cape Town were oestrogenic. E2 was detected in all the samples tested, but at higher concentrations at the eThekwini and City of Cape Town localities than Port Elizabeth Harbour. In addition, the recombinant yeast assays revealed that anti-androgenicity was widespread, being detected in the majority of samples screened apart from those representing Port Elizabeth Harbour. Conversely, no anti-oestrogenic or androgenic activity was detected. Anti-androgenicity did not associate with hydrocarbon loads, providing evidence that other anti-androgens were responsible for the observed activity. The present data suggest potential reproductive disruption in marine and estuarine fauna inhabiting the eThekwini and City of Cape Town regions.     

Spatially distributed tracer‐aided modelling to explore water and isotope transport,storage and mixing in a pristine,humid tropical catchment

Rapidly transforming headwater catchments in the humid tropics provide important resources for drinking water, irrigation, hydropower, and ecosystem connectivity. However, such resources for downstream use remain unstudied. To improve understanding of the behaviour and influence of pristine rainforests on water and tracer fluxes, we adapted the relatively parsimonious, spatially distributed tracer‐aided rainfall–runoff (STARR) model using event‐based stable isotope data for the 3.2‐km² San Lorencito catchment in Costa Rica. STARR was used to simulate rainforest interception of water and stable isotopes, which showed a significant isotopic enrichment in throughfall compared with gross rainfall. Acceptable concurrent simulations of discharge (Kling–Gupta efficiency [KGE] ~0.8) and stable isotopes in stream water (KGE ~0.6) at high spatial (10 m) and temporal (hourly) resolution indicated a rapidly responding system. Around 90% of average annual streamflow (2,099 mm) was composed of quick, near‐surface runoff components, whereas only ~10% originated from groundwater in deeper layers. Simulated actual evapotranspiration (ET) from interception and soil storage were low (~420 mm/year) due to high relative humidity (average 96%) and cloud cover limiting radiation inputs. Modelling suggested a highly variable groundwater storage (~10 to 500 mm) in this steep, fractured volcanic catchment that sustains dry season baseflows. This groundwater is concentrated in riparian areas as an alluvial–colluvial aquifer connected to the stream. This was supported by rainfall–runoff isotope simulations, showing a “flashy” stream response to rainfall with only a moderate damping effect and a constant isotope signature from deeper groundwater (~400‐mm additional mixing volume) during baseflow. The work serves as a first attempt to apply a spatially distributed tracer‐aided model to a tropical rainforest environment exploring the hydrological functioning of a steep, fractured‐volcanic catchment. We also highlight limitations and propose a roadmap for future data collection and spatially distributed tracer‐aided model development in tropical headwater catchments.     

The effect of possible climate change on natural groundwater recharge based on a simple model: a study of four karstic aquifers in SE Spain

The study presented in this paper constitutes an initial approach to the problematic task of evaluating the effects of possible climate change on natural water recharge to aquifers. To estimate such effects, a purpose-designed mathematical model termed Estimation of Recharge in Over-exploited Aquifers (ERAS) has been used. It enables to simulate the monthly water recharge to an aquifer, provided that prior knowledge of the exploitation to which it is subjected and the variation caused by these two actions on the piezometric level of the aquifer is available. The basic data required for its application are: precipitation, temperature, groundwater extraction, stored groundwater surface and storage coefficient. The main advantage presented by this model is its independence of the mechanism by which water is displaced through the ground and within the unsaturated zone. The ERAS code was applied to four over-exploited karstic aquifers in Alto Vinalopó (Alicante, Spain) with the goal of generating a synthesized series of values for natural groundwater recharge in each of the aquifers for the 100 years of the twentieth century. Each series thus obtained after being grouped into decades was subjected to statistical processing, which revealed that in every case a logarithmically decreasing trend was present.     

Post-injection trapping of mobile CO<Subscript><Emphasis Type="Bold">2</Emphasis></Subscript> in deep aquifers: Assessing the importance of model and parameter uncertainties

Predicting the fate of the injected CO₂ is crucial for the safety of carbon storage operations in deep saline aquifers: especially the evolution of the position, the spreading and the quantity of the mobile CO₂ plume during and after the injection has to be understood to prevent any loss of containment. Fluid flow modelling is challenging not only given the uncertainties on subsurface formation intrinsic properties (parameter uncertainty) but also on the modelling choices/assumptions for representing and numerically implementing the processes occurring when CO₂ displaces the native brine (model uncertainty). Sensitivity analysis is needed to identify the group of factors which contributes the most to the uncertainties in the predictions. In this paper, we present an approach for assessing the importance of model and parameter uncertainties regarding post-injection trapping of mobile CO₂. This approach includes the representation of input parameters, the choice of relevant simulation outputs, the assessment of the mobile plume evolution with a flow simulator and the importance ranking for input parameters. A variance-based sensitivity analysis is proposed, associated with the ACOSSO-like meta-modelling technique to tackle the issues linked with the computational burden posed by the use of long-running simulations and with the different types of uncertainties to be accounted for (model and parameter). The approach is tested on a potential site for CO₂ storage in the Paris basin (France) representative of a project in preliminary stage of development. The approach provides physically sound outcomes despite the challenging context of the case study. In addition, these outcomes appear very helpful for prioritizing the future characterisation efforts and monitoring requirements, and for simplifying the modelling exercise.     

The dust distribution within the inner coma of comet P/Halley 1982i: encounter by Giotto's impact detectors

Analysis of the data from Giotto's Dust Impact Detection System experiment (DIDSY) is presented. These data represent measurement of the size of dust grains incident on the Giotto dust shield along its trajectory through the coma of comet P/Halley on 1986 March 13/14. First detection occurred at some 287000 km distance from the nucleus on the inbound leg; the majority of the DIDSY subsystems remained operational after closest approach (604 km) yielding the last detection at about 202000 km from the nucleus. In order to improve the data coverage (and especially for the smallest grains, to approximately 10(-19) kg particle mass), data from the PIA instrument has been combined with DIDSY data. Flux profiles are presented for the various mass channels showing, to a first approximation, a 1/R2 flux dependence, where R is the distance of the detection point from the cometary nucleus, although significant differences are noted. Deviations from this dependence are observed, particularly close to the nucleus. From the flux profiles, mass and geometrical area distributions for the dust grains are derived for the trajectory through the coma. Groundbased CCD imaging of the dust continuum in the inner coma at the time of encounter is also used to derive the area of grains intercepted by Giotto. The results are consistent with the area functions derived by Giotto data and the low albedo of the grains deduced from infrared emission. For the close encounter period (-5 min to +5 min), the cumulative mass distribution function has been investigated, initially in 20 second periods; there is strong evidence from the data for a steepening of the index of the mass distribution for masses greater than 10(-13) kg during passage through dust jets which is not within the error limits of statistical uncertainty. The fluences for dust grains along the entire trajectory is calculated; it is found that extrapolation of the spectrum determined at intermediate masses (cumulative mass index alpha = 0.85) is not able to account for the spacecraft deceleration as observed by the Giotto Radio Science Experiment and by ESOC tracking operations. Data at large masses (>10(-8) kg) recently analysed from the DIDSY data set show clear evidence of a decrease in the mass distribution index at these masses within the coma, and it is shown that such a value of the mass index can provide sufficient mass for consistency with the observed deceleration. The total particulate mass output from the nucleus of comet P/Halley at the time of encounter would be dependent on the maximum mass emitted if this change in slope observed in the coma were also applicable to the emission from the nucleus; this matter is discussed in the text. The flux time profiles have been converted through a simple approach to modeling of the particle trajectories to yield an indication of nucleus surface activity. There is indication of an enhancement in flux at t approximately -29 s corresponding to crossing of the dawn terminator, but the flux detected prior to crossing of the dawn terminator is shown to be higher than predicted by simple modelling. Further enhancements corresponding to jet activity are detected around +190 s and +270 s.