An enhanced depth-averaged tidal model for morphological studies in the presence of rotary currents

A simple and efficient method to improve morphological predictions using depth-averaged tidal models is presented. The method includes the contribution of secondary flows in sediment transport using the computed flow field from a depth-averaged model. The method has been validated for a case study using the 3D POLCOMS model and ADCP data. The enhanced depth-averaged tidal model along with the SWAN wave model are applied to morphological prediction around the Lleyn Peninsula and Bardsey Island as a case study in the Irish Sea. Due to the presence of a headland in this area two asymmetrical tidal eddies are developed in which the cyclonic eddy is stronger as a result of Coriolis effects. The results show that the enhanced model can effectively predict formation of sand banks at the centre of cyclonic eddies, while the depth-averaged model, due to its inability to accommodate secondary flow, is inadequate in this respect.     

Groundwater recharge in a confined paleovalley setting,Northeast British Columbia,Canada

Hydrogeology Journal - Ancient river channels or subglacial drainage networks infilled with younger sediments can include significant deposits of highly permeable sands and gravels. Despite being...     

A model for assessing iceberg hazard

Natural Hazards - With the polar regions opening up to more marine activities but iceberg numbers more likely to increase than decline as a result of global warming, the risk from icebergs to...     

Nanogeochemistry of hydrothermal magnetite

Magnetite from hydrothermal ore deposits can contain up to tens of thousands of parts per million (ppm) of elements such as Ti, Si, V, Al, Ca, Mg, Na, which tend to either structurally incorporate into growth and sector zones or form mineral micro- to nano-sized particles. Here, we report micro- to nano-structural and chemical data of hydrothermal magnetite from the Los Colorados iron oxide–apatite deposit in Chile, where magnetite displays both types of trace element incorporation. Three generations of magnetites (X–Z) were identified with concentrations of minor and trace elements that vary significantly: SiO₂, from below detection limit (bdl) to 3.1 wt%; Al₂O₃, 0.3–2.3 wt%; CaO, bdl–0.9 wt%; MgO, 0.02–2.5 wt%; TiO₂, 0.1–0.4 wt%; MnO, 0.04–0.2 wt%; Na₂O, bdl–0.4 wt%; and K₂O, bdl–0.4 wt%. An exception is V₂O₃, which is remarkably constant, ranging from 0.3 to 0.4 wt%. Six types of crystalline nanoparticles (NPs) were identified by means of transmission electron microscopy in the trace element-rich zones, which are each a few micrometres wide: (1) diopside, (2) clinoenstatite; (3) amphibole, (4) mica, (5) ulvöspinel, and (6) Ti-rich magnetite. In addition, Al-rich nanodomains, which contain 2–3 wt% of Al, occur within a single crystal of magnetite. The accumulation of NPs in the trace element-rich zones suggest that they form owing to supersaturation from a hydrothermal fluid, followed by entrapment during continuous growth of the magnetite surface. It is also concluded that mineral NPs promote exsolution of new phases from the mineral host, otherwise preserved as structurally bound trace elements. The presence of abundant mineral NPs in magnetite points to a complex incorporation of trace elements during growth, and provides a cautionary note on the interpretation of micron-scale chemical data of magnetite.     

Cloud fraction at the ARM SGP site

The Atmospheric Radiation Measurement (ARM) Program Southern Great Plains (SGP) site has a rich history of actively sensed cloud observations. Fourteen years (1997–2010) of observations from the Millimeter Cloud Radar (MMCR), Micropulse Lidar (MPL), and Belfort/Vaisala Ceilometers are used to understand how instrument selection and sampling impacts estimates of Cloud Fraction (CF) at this location. Although all instruments should be used in combination for the best estimates of CF, instrument downtime limits available samples and increases observational errors, demanding that users make sacrifices when calculating CF at longer intervals relevant to climate studies. Selection of MMCR or MMCR + MPL cloud masks changes very little in the overall understanding of total CF. Addition of the MPL increases the 14-year average CF by 9 %, mainly through an increase in optically thin high clouds year-round, and mid-level clouds during the summer months. Splitting the period into two equal 7-year periods reveals negligible change in MMCR + MPL CF. For the MMCR, however, CF deceases by 6.1 %. This sudden change in CF occurs around the time the radar was upgraded, suggesting that this decrease is tied to hardware sensitivity or scanning strategy changes. Users must be cognizant of this and other issues when calculating CF from the variety of observations available at the ARM SGP site.     

Constraining the age and origin of the seamount province in the Northeast Indian Ocean using geophysical techniques

The breakup of western margin of Australia from Greater India started around 155 Ma and progressed southwards. After the separation, the interceding intraplate region experienced large volumes of submarine volcanism, extending over 100 Myrs. The Christmas Island Seamount Province (CHRISP, as it has been dubbed) lies south of the Java-Sunda Trench, and contains numerous submerged volcanic seamounts, and two sub-aerially exposed island groups—Cocos (Keeling) Islands, and Christmas Island. While recent geochronological investigations have shed light on the diverse eruption ages of the volcanics of this region, some islands/seamounts have demonstrated protracted volcanic histories, and it is not clear how the volcanic loading, tectonic subsidence, and subsequent emergence history of the islands relates to these discrete volcanic episodes. This study utilises a number of geophysical techniques to determine the crustal structure, loading and subsidence history, and last sub-aerial exposure age for the CHRISP. The study shows that flexural and subsidence modelling are reliable techniques in constraining the age of the seamounts when geochronological techniques are not possible. Utilising regional gravity signatures, we model the crustal structure underneath the Cocos (Keeling) Island, and constrain the thickness of the limestone cover between 900 and 2,100 m. Using age-depth subsidence curves for oceanic lithosphere the time since these seamounts were exposed above sea-level was determined, and a trend in exposure ages that youngs towards the west is observed. Two episodes of volcanism have been recorded at Christmas Island and they are of different origin. The younger phase in the Pliocene is a manifestation of flexure induced cracks produced in the lithosphere as it rides the subduction fore-bulge, whereas a low velocity seismic zone rising from the lower mantle, and tectonic reorganization, may be associated with the older Eocene volcanic phase, as well as much of the rest of the province. Our modelling also supports the existence of an older, undated volcanic core to Christmas Island, based on the loading ages from flexural modelling.     

Li isotopes and trace elements as a petrogenetic tracer in zircon: insights from Archean TTGs and sanukitoids

We report δ⁷Li, Li abundance ([Li]), and other trace elements measured by ion probe in igneous zircons from TTG (tonalite, trondhjemite, and granodiorite) and sanukitoid plutons from the Superior Province (Canada) in order to characterize Li in zircons from typical Archean continental crust. These data are compared with detrital zircons from the Jack Hills (Western Australia) with U–Pb ages greater than 3.9 Ga for which parent rock type is not known. Most of the TTG and sanukitoid zircon domains preserve typical igneous REE patterns and CL zoning. [Li] ranges from 0.5 to 79 ppm, typical of [Li] in continental zircons. Atomic ratios of (Y + REE)/(Li + P) average 1.0 ± 0.7 (2SD) for zircons with magmatic composition preserved, supporting the hypothesis that Li is interstitial and charge compensates substitution of trivalent cations. This substitution results in a relatively slow rate of Li diffusion. The δ⁷Li and trace element data constrain the genesis of TTGs and sanukitoids. [Li] in zircons from granitoids is significantly higher than from zircons in primitive magmas in oceanic crust. TTG zircons have δ⁷Li (3 ± 8‰) and δ¹⁸O in the range of primitive mantle-derived magmas. Sanukitoid zircons have average δ⁷Li (7 ± 8‰) and δ¹⁸O higher than those of TTGs supporting genesis by melting of fluid-metasomatized mantle wedge. The Li systematics in sanukitoid and TTG zircons indicate that high [Li] in pre-3.9-Ga Jack Hills detrital zircons is a primary igneous composition and suggests the growth in proto-continental crust in magmas similar to Archean granitoids.     

The importance of defining chemical potentials,substitution mechanisms and solubility in trace element diffusion studies: the case of Zr and Hf in olivine

The diffusion, substitution mechanism and solubility limits of Zr and Hf in synthetic forsterite (Mg₂SiO₄) and San Carlos olivine (Mg_0.9Fe_0.1)₂SiO₄ have been investigated between 1,200 and 1,500 °C as a function of the chemical potentials of the components in the system MgO(FeO)–SiO₂–ZrO₂(HfO₂). The effect of oxygen fugacity and crystallographic orientation were also investigated. The solubilities of Zr in forsterite are highest and diffusion fastest when the coexisting three-phase source assemblage includes ZrSiO₄ (zircon) or HfSiO₄ (hafnon), and lower and slower, respectively, when the source assemblage includes MgO (periclase). This indicates that Zr and Hf substitute on the octahedral sites in olivine, charge balanced by magnesium vacancies. Diffusion is anisotropic, with rates along the crystal axes increasing in the order a < b < c. The generalized diffusion relationship as a function of chemical activity (as \(a_{{{\text{SiO}}_{2} }}\)), orientation and temperature is: \(logD_{\text{Zr}} = \frac{1}{4}loga_{{{\text{SiO}}_{2} }} + logD_{0} - \left( {\frac{{368 \pm 17\;{\text{kJ}}\;{\text{mol}}^{ - 1} }}{{2.303\;{\text{RT}}}}} \right)\) where the values of log D ₀ are ?3.8(±0.5), ?3.4(±0.5) and ?3.1(±0.5) along the a, b and c axes, respectively. Most experiments were conducted in air (fO₂ = 10^?0.68 bars), but one at fO₂ = 10^?11.2 bars at 1,400 °C shows no resolvable effect of oxygen fugacity on Zr diffusion. Hf is slightly more soluble in olivine than Zr, but diffuses slightly slower. Diffusivities of Zr in experiments in San Carlos olivine at 1,400 °C, fO₂ = 10^?6.6 bars are similar to those in forsterite at the same conditions, showing that the controls on diffusivities are adequately captured by the simple system (nominally iron-free) experiments. Diffusivities are in good agreement with those measured by Spandler and O’Neill (Contrib Miner Petrol 159:791–818, 2010) in San Carlos olivine using silicate melt as the source at 1,300 °C, and fall within the range of most measurements of Fe–Mg inter-diffusion in olivine at this temperature. Forsterite–melt partitioning experiments in the CaO–MgO–Al₂O₃–SiO₂–ZrO₂/HfO₂ show that the interface concentrations from the diffusion experiments represent true equilibrium solubilities. Another test of internal consistency is that the ratios of the interface concentrations between experiments buffered by Mg₂SiO₄ + Mg₂Si₂O₆ + ZrSiO₄ or Mg₂SiO₄ + ZrSiO₄ + ZrO₂ (high silica activity) to those buffered by Mg₂SiO₄ + MgO + ZrO₂ (low silica activity) agree well with the ratios calculated from thermodynamic data. This study highlights the importance of buffering chemical potentials in diffusion experiments to provide constraints on the interface diffusant concentrations and hence validate the assumption of interface equilibrium.     

Structural controls on groundwater flow in a fractured bedrock aquifer underlying an agricultural region of northwestern New Brunswick,Canada

A hydrogeological study was conducted in northwestern New Brunswick, Canada, to improve the predictability of fracture-dominated groundwater flow within folded bedrock composed of fine-grained turbidites. Borehole televiewer logging and outcrop mapping, integrated with hydraulic packer tests revealed enhanced hydraulic conductivity associated with northeasterly striking bedding-plane fractures formed during folding and flexural slip. These fractures impart azimuthal anisotropy to the aquifer because of moderately dipping fold limbs. High-angle fractures form a well-developed non-stratabound network, comprising two open fracture sets striking NNE parallel to the current direction of principal stress, and WNW parallel to the direction of principal stress that dominated during the Acadian orogeny. The subset of fractures showing significant oxidation, deemed most important to the groundwater flow system, is dominated by bedding-plane and high-angle fractures striking near-parallel to the maximum principal stress direction, resulting in extensional opening and enhanced hydraulic conductivities. An equivalent porous media model, incorporating anisotropy and varying hydraulic conductivity with depth, indicates that horizontal flow dominates the aquifer with relatively minor exchange between different model layers. These findings have implications for understanding flow directions in the Black Brook Watershed and elsewhere in the Matapédia Basin where fractures formed under similar stress conditions.     

Water quality and thermal regime of the Motueka River: Influences of land cover,geology and position in the catchment

We examined the effects of land use, geology, and longitudinal position within the river network on water quality and thermal regime at 23 sites within the Motueka River catchment. The concentrations of suspended solids, nitrate nitrogen, total nitrogen, Escherichia coli, and Campylobacter were higher at sites draining pastoral and horticultural land than in similar‐sized native or plantation forest streams. Average daily mean temperature and minimum temperature in summer and maximum winter temperature were higher in unshaded pastoral and horticultural streams than at native forest sites. Differences in water quality and thermal regime were also observed among sites with contrasting geology. Conductivity, pH, and minimum winter temperatures were highest at sites draining marble terrain. In contrast, longitudinal patterns in water quality and temperature regime along the 120‐km length of the river were relatively weak, although longitudinal patterns in amplitude of daily temperature fluctuation matched theoretical predictions. In this study, differences in land use appeared to have the strongest influence on most water quality and thermal variables examined. However, geology was an important factor explaining variation in certain variables (e.g., pH and conductivity). Longitudinal patterns in water quality and temperature were relatively weak and in many instances were linked with longitudinal patterns in land use and geology rather than catchment location alone.