Ten best practices to strengthen stewardship and sharing of water science data in Canada

Water science data are a valuable asset that both underpins the original research project and bolsters new research questions, particularly in view of the increasingly complex water issues facing Canada and the world. Whilst there is general support for making data more broadly accessible, and a number of water science journals and funding agencies have adopted policies that require researchers to share data in accordance with the findable, accessible, interoperable, reusable (FAIR) principles, there are still questions about effective management of data to protect their usefulness over time. Incorporating data management practices and standards at the outset of a water science research project will enable researchers to efficiently locate, analyse and use data throughout the project lifecycle, and will ensure the data maintain their value after the project has ended. Here, some common misconceptions about data management are highlighted, along with insights and practical advice to assist established and early career water science researchers as they integrate data management best practices and tools into their research. Freely available tools and training opportunities made available in Canada through Global Water Futures, The Gordon Foundation DataStream, the Digital Research Alliance of Canada Portage Network, Compute Canada, and university libraries, among others are compiled. These include webinars, training videos, and individual support for the water science community that together enable researchers to protect their data assets and meet the expectations of journals and funders. The perspectives shared here have been developed as part of the Global Water Futures programme's efforts to improve data management and promote the use of common data practices and standards in the context of water science in Canada. Ten best practices are proposed that may be broadly applicable to other disciplines in the natural sciences and can be adopted and adapted globally.     

Absolute geomagnetic paleointensity as recorded by ∼1.09 Ga Lake Shore Traps (Keweenaw Peninsula, Michigan)

Absolute geomagnetic paleointensity measurements were made on 255 samples from 38 lava flows of the ～1.09 Ga Lake Shore Traps exposed on the Keweenaw Peninsula (Michigan, USA). Samples from the lava flows yield a well-defined characteristic remanent magnetization (ChRM) component within a ～375°C–590°C unblocking temperature range. Detailed rock magnetic analyses indicate that the ChRM is carried by nearly stoichiometric pseudo-single-domain magnetite and/or low-Ti titanomagnetite. Scanning electron microscopy reveals that the (titano)magnetite is present in the form of fine intergrowths with ilmenite, formed by oxyexsolution during initial cooling. Paleointensity values were determined using the Thellier double-heating method supplemented by low-temperature demagnetization in order to reduce the effect of magnetic remanence carried by large pseudosingle-domain and multidomain grains. One hundred and two samples from twenty independent cooling units meet our paleointensity reliability criteria and yield consistent paleofield values with a mean value of 26.3 ± 4.7μT, which corresponds to a virtual dipole moment of 5.9 ± 1.1×10²² Am². The mean and range of paleofield values are similar to those of the recent Earth’s magnetic field and incompatible with a “Proterozoic dipole low”. These results are consistent with a stable compositionally-driven geodynamo operating by the end of Mesoproterozoic.     

Analytical studies on transient groundwater flow induced by land reclamation using different fill materials

Land reclamation may have a significant influence on groundwater regimes. Analytical solutions have been developed in the past to study the impact of land reclamation on a steady‐state groundwater flow and transient flow in fill materials, assuming that the reclamation site consists of a single zone of uniform hydraulic parameters. In this paper, we derive analytical solutions to describe the transient water table change in response to multi‐stage land reclamation where the fill material is uniform in each stage but the hydraulic conductivity of the fill material varies from stage to stage. By introducing the method of separation of variables, we develop a transient analytical solution to study the impact of land reclamation consisting of fill material with different hydraulic properties on groundwater dynamics. The results show that the water table first increases significantly into the reclaimed zone following the fill material deposition, and then the increase gradually propagates into the original aquifer. The change of water table in the original aquifer mainly depends on the value of hydraulic conductivity of the fill materials. Examples in this paper illustrate how the aquifer system experiences a long time unsteady‐state flow as a result of the reclamation, and it takes at least tens of years for the system to approach a new equilibrium. It is suggested that for a large‐scale reclamation project, the response of the groundwater regime to reclamation should be carefully studied. Copyright © 2013 John Wiley & Sons, Ltd.     

Impact of coastal land reclamation on ground water level and the sea water interface

Land reclamation in coastal areas may have a significant effect on local ground water systems. Steady-state analytic solutions based on Dupuit and Ghyben-Herzberg assumptions are derived to evaluate this effect. Two situations are considered, both with ground water flow resulting from precipitation recharge: the coastal aquifer of an extensive landmass and an island. The results show that after reclamation, the water table rises and the salt water-fresh water interface moves seaward. The degree of these changes depends on the extent of reclamation and the hydraulic conductivity of the fill material. For the island situation, the reclamation displaces the ground water divide and changes the ground water conditions in the entire island. An unintended advantage of the reclamation is an increase of fresh ground water resource because the reclaimed land can be an additional aquifer and rain recharge takes place over a larger area.     

Macrobenthic Community in Tolo Harbour,Hong Kong and its Relations with Heavy Metals

The present study investigated the macrobenthic community in Tolo Harbour, Hong Kong, aiming at linking heavy metal concentrations to differences in macrobenthic community. The stations investigated in Tolo Harbour have widely contrasting features, with some areas located in the Plove Cove displaying both high species richness and abundance while other areas displaying quite impoverished or even void of macrobenthos. High diversity and abundance of macrobenthos in areas with low heavy metal concentrations were recorded. Strong negative correlation between macrobenthic diversity and heavy metal concentrations was found, and this implicated the pollution-induced degradation of macrobenthos in some locations in Tolo Harbour. These results support the Pearson–Rosenberg model for succession along a pollution gradient.     

Estimation of fault-zone conductance by calibration of a regional groundwater flow model: Desert Hot Springs,California

The hydraulic conductance of a large fault zone has been estimated by calibrating a regional groundwater flow model. Drops in groundwater elevations of over 80 m have been observed along a 15-km length of the Mission Creek fault, California, USA. The large drops in elevation are attributed to the reduced hydraulic conductivity of the fault materials. A conceptual and numerical model of the two hydrologic subbasins in Desert Hot Springs, separated by the Mission Creek fault, was developed. The model was used to estimate the hydraulic conductance along the fault. The parameter estimation involved calibrating the model with observed groundwater elevations from over 40 locations over a 60-year period. The fault hydraulic conductances were estimated assuming a linear trend in the fault length, yielding variations in the fault hydraulic conductance of about an order of magnitude along the fault length (2?×?10^?11–4?×?10^?10 1/s). When an average fault thickness of 35 m is assumed, the fault hydraulic conductivity values are estimated to be from three to five orders of magnitude lower than the surrounding materials. A sensitivity analysis indicated that assumptions made in the conceptual model do not significantly affect estimated fault hydraulic conductances.     

The influence of foliation on the fragility of granitic rocks, image analysis and quantitative microscopy

Fabric-dependent anisotropy is important in understanding the mechanical behaviour of foliated rocks because the foliation creates mechanically weak discontinuities. Using optical microscopy, a foliation index (FIX) was calculated for a group of granitic rocks in southern Sweden that displays a range of textures from foliated to nonfoliated. Image analysis from scanning electron microscope with backscattered detector (SEM/BSE) images was used in order to measure the rock texture of these samples. The analyses were compared with the Los Angeles (LA) test, which measures resistance of rock aggregate to fragmentation. The results demonstrate that the measured grain size, shape and spatial arrangement are important parameters controlling the fragility of the rock. Both the texture and foliation must be taken into account to obtain a significant relationship with the mechanical analyses. We suggest that this alternative method could be a helpful tool for predicting a rock's suitability as an aggregate as it gives a better understanding of the material properties compared to mechanical testing.     

Paleomagnetic constraints on eruption patterns at the Pacaya composite volcano,Guatemala

Pacaya volcano is an active composite volcano located in the volcanic highlands of Guatemala about 40 km south of Guatemala City. Volcanism at Pacaya alternates between Strombolian and Vulcanian, and during the past five years there has been a marked increase in the violence of eruptions. The volcano is composed principally of basalt flows interbedded with thin scoria fall units, several pyroclastic surge beds, and at least one welded tuff. Between 400 and 2000 years BP the W-SW sector of the volcano collapsed producing a horseshoeshaped amphitheater (0.65 km³) and providing a window into the cone's infrastructure. Lava flows and tephra exposed in the amphitheater are more then 200 m thick and when combined with flows erupted recently represent between 30 and 40% of the cone's history. Pacaya is ideally suited for a paleomagnetic study into the timing and duration of eruption episodes at a large, composite volcano. We drilled 27 paleomagnetic sites (25 aa flows, 1 dike, and 1 welded tuff) from four lava-flow sequences with between 4 and 14 sites per sequence. The four sequences represent initial through historic activity at Pacaya. We resolved, what appear to be, 22 time-independent paleomagnetic sites by averaging together directions from successive sites where the sitemean directions were indistinguishable at the 95% level of confidence. However, mean-sequence directions of individual lava-flow sequences yielded unusually high Fisher precision parameters (k=44–224) and small circles of 63% confidence (a₆₃=1.6–6.1°) suggesting as few as three or four time-independent sites were collected. This indicates that activity as Pacaya is strongly episodic and that episodes are characterized by voluminous outpouring of lavas. Modelling the data using Holocene PSV rates confirms this and shows that differences in within-sequence directions (6–11.5°) are consistent with emplacement of lava-flow sequences in less than 100 years to as many as 300 years. Relatively larger differences in directions (18–23°) between subjacent lava-flow sequences indicates that repose is at least 300–500 years and could be even longer.     

Optimal weight absorber designs for vibrating structures exposed to random excitations

Optimal mass ratios that minimize the response of a laminated beam with an attached absorber are tabulated for various values of beam damping. The beam is treated as an equivalent one degree of freedom (1DOF) main system vibrating in the fundamental mode. The beam is subjected to Gaussian white noise force and Gaussian white noise base frame acceleration. Optimal absorber frequency ratios and absorber damping ratios have been tabulated by others; the results for the classical 1DOF main system with attached absorber suggest that the optimized non-dimensional response decreases monotonically as the mass ratio increases. However, to generalize this monotonic relation may lead to inappropriate conclusions. If we define a constraint such that an increase in absorber mass leads to a proportional decrease in available beam construction material, i.e. effectively the combined mass of the beam and absorber is minimized, then variations in the mass ratio will affect the beam's parameters such as mass, stiffness and damping. Since some of these parameters are used for non-dimensionalising the response, inspection of non-dimensional responses may in some cases lead to inappropriate conclusions. This paper shows the optimal mass ratios for minimizing the response of a structure exposed to earthquake or fluid flow type random excitations.