Characterisation of conduction phenomena in soils at the particle-scale: Finite element analyses in conjunction with synthetic 3D imaging

This paper describes an innovative method to characterise conduction parameters in geomaterials at the particle-scale. The technique is exemplified using 3D synthetic grain packing generated with discrete element approaches. This creates a geo-mechanically viable user-defined 3D granular image through which the particle skeleton and the corresponding pore network are constructed. Images are then imported into the finite element analyses to solve the governing equations of hydraulic and thermal conduction. Navier–Stokes equation is uniquely upscaled to Darcy’s law to assess hydraulic conductivity in soils, while a similar approach implements the Fourier equation to evaluate thermal conduction through grain chains and pore network. High performance computing is used to meet demanding numerical calculations of 3D meshed geometries. Packing density (i.e., porosity) and inter-particle contact areas are explored as variables to highlight the effects of pore volume and inter-particle contact condition in hydraulic and thermal conduction. This emerging technique allows not only characterising the macro-scale behaviour of conduction phenomena in soils but also quantifying and visualising the preferential and local conduction behaviour at the particle-scale. Laboratory measurements of hydraulic and thermal conductivities support numerically obtained results and validate the viability of the new methods used herein. This study introduces an alternative way to determine physical parameters of soils using emerging technology of rigorous numerical simulations in conjunction with 3D images, and to enable fundamental observation of particle-scale mechanisms of macro-scale manifestation.     

Proterozoic low-sulfidation epithermal Au-Ag mineralization in the Mallery Lake area, Nunavut, Canada

The Mallery Lake area contains pristine examples of ancient precious metal-bearing low-sulfidation epithermal deposits. The deposits are hosted by rhyolitic flows of the Early Proterozoic Pitz Formation, but are themselves apparently of Middle Proterozoic age. Gold mineralization occurs in stockwork quartz veins that cut the rhyolites, and highest gold grades (up to 24 g/t over 30 cm) occur in the Chalcedonic Stockwork Zone. Quartz veining occurs in two main types: barren A veins, characterized by fine- to coarse-grained comb quartz, with fluorite, calcite, and/or adularia; and mineralized B veins, characterized by banded chalcedonic silica and fine-grained quartz, locally intergrown with fine-grained gold or electrum. A third type of quartz vein (C), which crosscuts B veins at one locality, is characterized by microcrystalline quartz intergrown with fine-grained hematite and rare electrum. Fluid inclusions in the veins occur in two distinct assemblages. Assemblage 1 inclusions represent a moderate temperature (Th=150 to 220 °C), low salinity (~1 eq. wt% NaCl, with trace CO₂), locally boiling fluid; this fluid type is found in both A and B veins and is thought to have been responsible for Au-Ag transport and deposition. Assemblage 2 inclusions represent a lower temperature (Th=90 to 150 °C), high salinity calcic brine (23 to 31 wt% CaCl₂-NaCl), which occurs as primary inclusions only in the barren A veins. Assemblage 1 and 2 inclusions occur in alternating quartz growth bands in the A-type veins, where they appear to represent alternating fluxes of dilute fluid and local saline groundwater. No workable primary fluid inclusions were observed in the C veins. The A-vein quartz yields '¹⁸O values from 8.3 to 14.5‰ (average=10.9ǃ.7‰ [1C], n=30), whereas '¹⁸O values for B-vein quartz range from 11.2 to 14.0‰ (average=13.0ǂ.9‰, n=12). Calculated '¹⁸O_H2O values for the dilute mineralizing fluid from B veins range from -2.6 to 0.2‰ (average=-0.8ǂ.9‰, n=12) and are consistent with a dominantly meteoric origin. No values could be calculated for the brine, however, because all A-vein quartz samples contain mixed fluid inclusion populations. However, the fact that A-vein quartz samples extend to lower '¹⁸O values than the B veins suggests that the brine had a lighter isotopic signature relative to the dilute fluid. Hydrogen isotopic ratios of fluid inclusion waters extracted from eleven quartz samples of both vein types range from 'D_FI=-56 to -134‰, but show no particular correlation with vein type. In most respects, the mineralogical and fluid characteristics of the Mallery Lake system are comparable to those of Phanerozoic low-sulfidation deposits, and although the presence of high salinity brines is unusual in such deposits, it is not unknown (e.g., Creede, Colorado). In addition, one of the few other examples of well-preserved, Precambrian, low-sulfidation epithermal deposits, from the Central Pilbara tectonic zone, Australia, contains a similarly bimodal fluid assemblage. The significance of these saline brines is not clear, but from this study we infer that they were not directly involved with Au-Ag transport or deposition.     

Intra-annual variability of urban effects on streamflow

While considerable research has established the impacts of urbanization on streamflow, there has been little emphasis on how intra-annual variations in streamflow can deepen the understanding of hydrological processes in urban watersheds. This study fills this critical research gap by examining, at the monthly scale, correlations between land-cover and streamflow, differences in streamflow metrics between urban and rural watersheds, and the potential for the inflow and infiltration (I&I) of extraneous water into sewers to reduce streamflow. We use data from 90 watersheds in the Atlanta, GA region over the 2013–2019 period to accomplish our objectives. Similar to other urban areas in temperate climates, Atlanta has a soil-water surplus in winter and a soil-water deficit in summer. Our results show urban watersheds have less streamflow seasonality than do rural watersheds. Compared to rural watersheds, urban watersheds have a much larger frequency of high-flow days during July–October. This is caused by increased impervious cover decreasing the importance of antecedent soil moisture in producing runoff. Urban watersheds have lower baseflows than rural watersheds during December–April but have baseflows equal to or larger than baseflows in rural watersheds during July–October. Intra-annual variations in effluent data from wastewater treatment plants provide evidence that I&I is a major cause of the relatively low baseflows during December–April. The relatively high baseflows in urban watersheds during July–October are likely caused by reduced evapotranspiration and the inflow of municipal water. The above seasonal aspects of urban effects on streamflow should be applicable to most urban watersheds with temperate climates.     

An examination of the rheology of flocculated clay suspensions

A dense cohesive sediment suspension, sometimes referred to as fluid mud, is a thixotropic fluid with a true yield stress. Current rheological formulations struggle to reconcile the structural dynamics of cohesive sediment suspensions with the equilibrium behaviour of these suspensions across the range of concentrations and shear. This paper is concerned with establishing a rheological framework for the range of sediment concentrations from the yield point to Newtonian flow. The shear stress equation is based on floc fractal theory, put forward by Mills and Snabre (1988). This results in a Casson-like rheology equation. Additional structural dynamics is then added, using a theory on the self-similarity of clay suspensions proposed by Coussot (1995), giving an equation which has the ability to match the equilibrium and time-dependent viscous rheology of a wide range of suspensions of different concentration and mineralogy.     

Strengthened Regulation of the Onset of the South China Sea Summer Monsoon by the Northwest Indian Ocean Warming in the Past Decade

Traditionally,a delayed(early)onset of the South China Sea summer monsoon(SCSSM)has been observed to follow a warm(cold)El Ni?o-Southern Oscillation(ENSO)event in winter,supporting high seasonal predictability of SCSSM onset.However,the empirical seasonal forecasting skill of the SCSSM onset,solely based on ENSO,has deteriorated since 2010.Meanwhile,unexpected delayed onsets of the SCSSM have also occurred in the past decade.We attribute these changes to the Northwest Indian Ocean(NWIO)warming of the sea surface.The NWIO warming has teleconnections related to(1)suppressing the seasonal convection over the South China Sea,which weakens the impacts of ENSO on SCSSM onset and delays the start of SCSSM,and(2)favoring more high-frequency,propagating moist convective activities,which enhances the uncertainty of the seasonal prediction of SCSSM onset date.Our results yield insight into the predictability of the SCSSM onset under the context of uneven ocean warming operating within the larger-scale background state of global climate change.     

Cytochrome P450 1A1 expression in cetacean skin biopsies from the Indian Ocean

The study describes cytochrome P450 1A1 (CYPA1) expression in the skin of different cetacean species (Megaptera novaeangliae, n = 15; Stenella attenuata, n = 7 and Stenella longirostris, n = 24) from the Mozambique Channel island of Mayotte. Immunohistochemical examination was performed with a monoclonal antibody against scup cytochrome CYPA1. The sex was determined using a molecular approach consisting in the genotyping sex-specific genes. CYPA1 was detected at the junction between epidermis and blubber on dolphins only, mostly in the endothelial cells. Similar observation was obtained in the dermis of one M. novaeangliae. Immunohistochemical slides were scored to evaluate the expression of the CYPA1 and a higher expression was observed in S. longirostris, suggesting a higher exposure to pollutants for this species. The difference of expression between sexes was not significant.     

Sensitivity of Simulated Hyporheic Exchange to River Bathymetry: The Steinlach River Test Site

This study determines the aspects of river bathymetry that have the greatest influence on the predictive biases when simulating hyporheic exchange. To investigate this, we build a highly parameterized HydroGeoSphere model of the Steinlach River Test Site in southwest Germany as a reference. This model is then modified with simpler bathymetries, evaluating the changes to hyporheic exchange fluxes and transit time distributions. Results indicate that simulating hyporheic exchange with a high-resolution detailed bathymetry using a three-dimensional fully coupled model leads to nested multiscale hyporheic exchange systems. A poorly resolved bathymetry will underestimate the small-scale hyporheic exchange, biasing the simulated hyporheic exchange towards larger scales, thus leading to overestimates of hyporheic exchange residence times. This can lead to gross biases in the estimation of a catchment's capacity to attenuate pollutants when extrapolated to account for all meanders along an entire river within a watershed. The detailed river slope alone is not enough to accurately simulate the locations and magnitudes of losing and gaining river reaches. Thus, local bedforms in terms of bathymetric highs and lows within the river are required. Bathymetry surveying campaigns can be more effective by prioritizing bathymetry measurements along the thalweg and gegenweg of a meandering channel. We define the gegenweg as the line that connects the shallowest points in successive cross-sections along a river opposite to the thalweg under average flow conditions. Incorporating local bedforms will likely capture the nested nature of hyporheic exchange, leading to more physically meaningful simulations of hyporheic exchange fluxes and transit times.     

A new laboratory technique for determining the compressional wave properties of marine sediments at sonic frequencies and in situ pressures

We describe a new laboratory technique for measuring the compressional wave velocity and attenuation of jacketed samples of unconsolidated marine sediments within the acoustic (sonic) frequency range 1–10 kHz and at elevated differential (confining – pore) pressures up to 2.413 MPa (350 psi). The method is particularly well suited to attenuation studies because the large sample length (up to 0.6 m long, diameter 0.069 m) is equivalent to about one wavelength, thus giving representative bulk values for heterogeneous samples. Placing a sediment sample in a water‐filled, thick‐walled, stainless steel Pulse Tube causes the spectrum of a broadband acoustic pulse to be modified into a decaying series of maxima and minima, from which the Stoneley and compressional wave, velocity and attenuation of the sample can be determined. Experiments show that PVC and copper jackets have a negligible effect on the measured values of sediment velocity and attenuation, which are accurate to better than ± 1.5% for velocity and up to ± 5% for attenuation. Pulse Tube velocity and attenuation values for sand and silty‐clay samples agree well with published data for similar sediments, adjusted for pressure, temperature, salinity and frequency using standard equations. Attenuation in sand decreases with pressure to small values below Q^?1 = 0.01 (Q greater than 100) for differential pressures over 1.5 MPa, equivalent to sub‐seafloor depths of about 150 m. By contrast, attenuation in silty clay shows little pressure dependence and intermediate Q^?1 values between 0.0206–0.0235 (Q = 49–43). The attenuation results fill a notable gap in the grain size range of published data sets. Overall, we show that the Pulse Tube method gives reliable acoustic velocity and attenuation results for typical marine sediments.     

Marine debris collects within the North Pacific Subtropical Convergence Zone

Floating marine debris, particularly derelict fishing gear, is a hazard to fish, marine mammals, turtles, sea birds, coral reefs, and even human activities. To ameliorate the economic and environmental impact of marine debris, we need to efficiently locate and retrieve dangerous debris at sea. Guided by satellite-derived information, we made four flights north of Hawaii in March and April 2005. During these aerial surveys, we observed over 1800 individual pieces of debris, including 122 derelict fishing nets. The largest debris concentrations were found just north of the North Pacific Transition Zone Chlorophyll Front (TZCF) within the North Pacific Subtropical Convergence Zone (STCZ). Debris densities were significantly correlated with sea-surface temperature (SST), chlorophyll-a concentration (Chla), and the gradient of Chla. A Debris Estimated Likelihood Index (DELI) was developed to predict where high concentrations of debris would be most likely in the North Pacific during spring and early summer.     

Circulation of bubbly magma and gas segregation within tunnels of the potential Yucca Mountain repository

Following an intersection of rising magma with drifts of the potential Yucca Mountain nuclear waste repository, a pathway is likely to be established to the surface with magma flowing for days to weeks and affecting the performance of engineered structures located along or near the flow path. In particular, convective circulation could occur within magma-filled drifts due to the exsolution and segregation of magmatic gas. We investigate gas segregation in a magma-filled drift intersected by a vertical dyke by means of analogue experiments, focusing on the conditions of sustained magma flow. Degassing is simulated by electrolysis, producing micrometric bubbles in viscous mixtures of water and golden syrup, or by aerating golden syrup, producing polydisperse bubbly mixtures with 40% of gas by volume. The presence of exsolved bubbles induces a buoyancy-driven exchange flow between the dyke and the drift that leads to gas segregation. Bubbles segregate from the magma by rising and accumulating as a foam at the top of the drift, coupled with the accumulation of denser degassed magma at the base of the drift. Steady-state influx of bubbly magma from the dyke into the drift is balanced by outward flux of lighter foam and denser degassed magma. The length and time scales of this gas segregation are controlled by the rise of bubbles in the horizontal drift. Steady-state gas segregation would be accomplished within hours to hundreds of years depending on the viscosity of the degassed magma and the average size of exsolved gas bubbles, and the resulting foam would only be a few cm thick. The exchange flux of bubbly magma between the dyke and the drift that is induced by gas segregation ranges from 1 m³ s⁻¹, for the less viscous magmas, to 10⁻⁸ m³ s⁻¹, for the most viscous degassed magmas, with associated velocities ranging from 10⁻¹ to 10⁻⁹ m s⁻¹ for the same viscosity range. This model of gas segregation also predicts that the relative proportion of erupted degassed magma, that could potentially carry and entrain nuclear waste material towards the surface, would depend on the value of the dyke magma supply rate relative to the value of the gas segregation flux, with violent eruption of gassy as well as degassed magmas at relatively high magma supply rates, and eruption of mainly degassed magma by milder episodic Strombolian explosions at relatively lower supply rates.