An overview of toxicant identification in sediments and dredged materials

The identification of toxicants affecting aquatic benthic systems is critical to sound assessment and management of our nation's waterways. Identification of toxicants can be useful in designing effective sediment remediation plans and reasonable options for sediment disposal. Knowledge of which contaminants affect benthic systems allows managers to link pollution to specific dischargers and prevent further release of toxicant(s). In addition, identification of major causes of toxicity in sediments may guide programs such as those developing environmental sediment guidelines and registering pesticides, while knowledge of the causes of toxicity which drive ecological changes such as shifts in benthic community structure would be useful in performing ecological risk assessments. To this end, the US Environmental Protection Agency has developed tools (toxicity identification and evaluation (TIE) methods) that allow investigators to characterize and identify chemicals causing acute toxicity in sediments and dredged materials. To date, most sediment TIEs have been performed on interstitial waters. Preliminary evidence from the use of interstitial water TIEs reveals certain patterns in causes of sediment toxicity. First, among all sediments tested, there is no one predominant cause of toxicity; metals, organics, and ammonia play approximately equal roles in causing toxicity. Second, within a single sediment there are multiple causes of toxicity detected; not just one chemical class is active. Third, the role of ammonia is very prominent in these interstitial waters. Finally, if sediments are divided into marine or freshwater, TIEs perforMed on interstitial waters from freshwater sediments indicate a variety of toxicants in fairly equal proportions, while TIEs performed on interstitial waters from marine sediments have identified only ammonia and organics as toxicants, with metals playing a minor role. Preliminary evidence from whole sediment TIEs indicates that organic compounds play a major role in the toxicity of marine sediments, with almost no evidence for either metal or ammonia toxicity. However, interpretation of these results may be skewed because only a small number of interstitial water (n = 13) and whole sediment (n = 5) TIEs have been completed. These trends may change as more data are collected.     

Mbosi: An anomalous iron with unique silicate inclusions

Abstract— The Mbosi iron meteorite contains millimeter size silicate inclusions. Mbosi is an ungrouped iron meteorite with a Ge/Ga ratio >10, which is an anomalous property shared with the five-member IIF iron group, the Eagle Station pallasites and four other ungrouped irons. Neither the IIF group nor the four other ungrouped irons are known to have silicate inclusions. Chips from three Mbosi inclusions were studied, but most of the work concentrated on a whole 3.1 mm circular inclusion. This inclusion consists of a mantle and a central core of different mineralogies. The mantle is partially devitrified quartz-normative glass, consisting of microscopic crystallites of two pyroxenes and plagioclase, which are crystalline enough to give an x-ray powder diffraction pattern but not coarse enough to permit analyses of individual minerals. The core consists of silica. The bulk composition does not match any known meteorite type, although there is a similarity in mode of occurrence to quartz-normative silicate inclusions in some HE irons. Mbosi silicate appears to be unique. The bulk rare earth element (REE) pattern of the mantle is flat at ? 7×C1; the core is depleted in REE but shows a small positive Eu anomaly. The O-isotope composition of bulk silicate lies on a unit slope mixing line (parallel and close to the C3 mixing line) that includes the Eagle Station pallasites and the iron Bocaiuva (related to the IIF irons); all of these share the property of having Ge/Ga ratios >10. It is concluded that Mbosi silicate represents a silica-bearing source rock that was melted and injected into metal. Melting occurred early in the history of the parent body because the metal now shows a normal Widmanstätten structure with only minor distortion that was caused when the parent body broke up and released meteorites into interplanetary space. The cause of Ge/Ga ratios being >10 in these irons is unknown. The fact that silicates in Mbosi, Bocaiuva (related to IIF irons) and the Eagle Station trio of pallasites, all characterized by a Ge/Ga ratio >10, lie on a unit slope mixing line in the O-isotope diagram suggests that their origins are closely related. The C3 chondrites appear to be likely precursors for silicates in Mbosi, Bocaiuva and the Eagle Station pallasites.     

Trends in Canadian Short‐Duration Extreme Rainfall: Including an Intensity–Duration–Frequency Perspective

Short-duration (5 minutes to 24 hours) rainfall extremes are important for a number of purposes, including engineering infrastructure design, because they represent the different meteorological scales of extreme rainfall events. Both single location and regional analyses of the changes in short-duration extreme rainfall amounts across Canada, as observed by tipping bucket rain gauges from 1965 to 2005, are presented. The single station analysis shows a general lack of a detectable trend signal, at the 5% significance level, because of the large variability and the relatively short period of record of the extreme short-duration rainfall amounts. The single station 30-minute to 24-hour durations show that, on average, 4% of the total number of stations have statistically significant increasing amounts of rainfall, whereas 1.6% of the cases have significantly decreasing amounts. However, regional spatial patterns are apparent in the single station trend results. Thus, for the same durations regional trends are presented by grouping the single station trend statistics across Canada. This regional trend analysis shows that at least two-thirds of the regions across Canada have increasing trends in extreme rainfall amounts, with up to 33% being significant (depending on location and duration). Both the southwest and the east (Newfoundland) coastal regions generally show significant increasing regional trends for 1- and 2-hour extreme rainfall durations. For the shortest durations of 5–15 minutes, the general overall regional trends in the extreme amounts are more variable, with increasing and decreasing trends occurring with similar frequency; however, there is no evidence of statistically significant decreasing regional trends in extreme rainfall amounts. The decreasing regional trends for the 5- to 15-minute duration amounts tend to be located in the St. Lawrence region of southern Quebec and in the Atlantic provinces. Additional analysis using criteria specified for traditional water management practice (e.g., Intensity-Duration-Frequency (IDF)) shows that fewer than 5.6% and 3.4% of the stations have significant increasing and decreasing trends, respectively, in extreme annual maximum single location observation amounts. This indicates that at most locations across Canada the traditional single station IDF assumption that historical extreme rainfall observations are stationary (in terms of the mean) over the period of record for an individual station is not violated. However, the trend information is still useful complementary information that can be considered for water management purposes, especially in terms of regional analysis.     

Including spatial distribution in a data‐driven rainfall‐runoff model to improve reservoir inflow forecasting in Taiwan

Multi‐step ahead inflow forecasting has a critical role to play in reservoir operation and management in Taiwan during typhoons as statutory legislation requires a minimum of 3‐h warning to be issued before any reservoir releases are made. However, the complex spatial and temporal heterogeneity of typhoon rainfall, coupled with a remote and mountainous physiographic context, makes the development of real‐time rainfall‐runoff models that can accurately predict reservoir inflow several hours ahead of time challenging. Consequently, there is an urgent, operational requirement for models that can enhance reservoir inflow prediction at forecast horizons of more than 3 h. In this paper, we develop a novel semi‐distributed, data‐driven, rainfall‐runoff model for the Shihmen catchment, north Taiwan. A suite of Adaptive Network‐based Fuzzy Inference System solutions is created using various combinations of autoregressive, spatially lumped radar and point‐based rain gauge predictors. Different levels of spatially aggregated radar‐derived rainfall data are used to generate 4, 8 and 12 sub‐catchment input drivers. In general, the semi‐distributed radar rainfall models outperform their less complex counterparts in predictions of reservoir inflow at lead times greater than 3 h. Performance is found to be optimal when spatial aggregation is restricted to four sub‐catchments, with up to 30% improvements in the performance over lumped and point‐based models being evident at 5‐h lead times. The potential benefits of applying semi‐distributed, data‐driven models in reservoir inflow modelling specifically, and hydrological modelling more generally, are thus demonstrated. Copyright © 2012 John Wiley & Sons, Ltd.     

Three-dimensional subsatellite motion under air drag and oblateness perturbations

The three-dimensional relative motion of a subsatellite with respect to a reference station in an elliptical orbit is studied. A general theory based on the variation of the relative elements, i.e. the instantaneous differences between the orbital parameters of the subsatellite and those of the station, is formulated in order to incorporate arbitrary perturbing forces acting on both satellites. The loss of precision inherent in the subtraction of almost identical quantities is avoided by the consistent use of difference variables. In the absence of perturbations exact analytical representations can be obtained for the relative state parameters. The influences of air drag and Earth's oblateness on the relative motion trajectories are investigated and illustrated graphically for a number of cases.     

Numerical Simulation of Strongly Stratified Flow Over a Three-Dimensional Hill

A numerical study of stably stratified flow over a three-dimensional hill is presented. Large-eddy simulation is used here to examine in detail the laboratory experimental flows described in the landmark work of Hunt and Snyder about stratified flow over a hill. The flow is linearly stratified and U/Nh is varied from 0.2 to 1.0. Here N and U are the buoyancy frequency and freestream velocity respectively, and h is the height of the hill. The Reynolds number based on the hill height is varied from 365 to 2968. The characteristic flow patterns at various values of U/Nh have been obtained and they are in good agreement with earlier theoretical and experimental results. It is shown that the flow field cannot be predicted by Drazin's theory when recirculation exists at the leeside of the hill even at UNh 1. The wake structure agrees well with a two-dimensional wake assumption when U/Nh 1 but lee waves start to influence the wake structure as U/Nh increases. The dividing-streamline heights obtained in the simulation are in accordance with experimental results and Sheppard's formula. The energy loss along the dividing streamline due to friction/turbulence approximately offsets the energy gained from pressure field. When lee waves are present, linear theory always underestimates the amplitude and overestimates the wavelength of three-dimensional lee waves. The simulated variations of drag coefficients with the parameterK (=ND/ U) are qualitatively consistent with experimental data and linear theory. Here D is the depth of the tank.     

Letters to the editor

A preliminary scuba survey of oligotrophic Lake Rotoma in 1972 revealed a vegetation mainly composed of native hydrophytes in which exotics were at an early stage of colonisation. In 1973 the presence of species was recorded in 5708 quadrats (625 cm²) at 1 m intervals along a total of 50 line transects placed systematically around the lake. Water depth was measured, and quadrat cover and substrate type were subjectively estimated. Species frequency calculations showed that the dominant vegetation pattern was a characean meadow of Chara fibrosa f. acanthopitys (A.Br.) R.D.W., Nitella leptostachys var. leonhardii (R.D.W.) R.D.W., and N. pseudoflabellata var. mucosa (Nordst.) Bailey. The charophytes extended over a depth range of 1–17.5 m on a wide variety of substrates and gradients. Native vascular plants were absent from many transects, and had a depth range only from 0 to 4.5 m, with most occurring above 3.5 m. The Low Mixed Community, found in shallow water less than 1.25 m in depth at the northeast end of the lake, provided this area with a high species diversity. Exotic hydrophytes had established in many areas around the lake. The distribution of Lagarosiphon major (Ridley) Moss and Elodea canadensis Michx. appeared to coincide with boating access and fallen submerged trees over a depth range of 0–6.0 m, although much of the available habitat had not yet been exploited. Emergent species were most abundant within the southwest inlet and also in the lagoons surrounding the lake where sheltered conditions and shallow gradients prevail.     

Flow resistance and hydraulic geometry in bedrock rivers with multiple roughness length scales

Many models of incision by bedrock rivers predict water depth and shear stress from discharge; conversely, palaeoflood discharge is sometimes reconstructed from flow depth markers in rock gorges. In both cases, assumptions are made about flow resistance. The depth–discharge relation in a bedrock river must depend on at least two roughness length scales (exposed rock and sediment cover) and possibly a third (sidewalls). A conceptually attractive way to model the depth–discharge relation in such situations is to partition the total shear stress and friction factor, but it is not obvious how to quantify the friction factor for rough walls in a way that can be used in incision process models. We show that a single flow resistance calculation using a spatially averaged roughness length scale closely approximates the partitioning of stress between sediment and rock, and between bed and walls, in idealized scenarios. Both approaches give closer fits to the measured depth–discharge relations in two small bedrock reaches than can be achieved using a fixed value of Manning's n or the Chézy friction factor. Sidewalls that are substantially rougher or smoother than the bed have a significant effect on the partitioning of shear stress between bed and sidewalls. More research is needed on how best to estimate roughness length scales from observable or measurable channel characteristics. © 2019 John Wiley & Sons, Ltd.