Fault Roughness at Seismogenic Depths from LIDAR and Photogrammetric Analysis

Fault surface roughness is a principal factor influencing earthquake mechanics, and particularly rupture initiation, propagation, and arrest. However, little data currently exist on fault surfaces at seismogenic depths. Here, we investigate the roughness of slip surfaces from the seismogenic strike-slip Gole Larghe Fault Zone, exhumed from ca. 10 km depth. The fault zone exploited pre-existing joints and is hosted in granitoid rocks of the Adamello batholith (Italian Alps). Individual seismogenic slip surfaces generally show a first phase of cataclasite production, and a second phase with beautifully preserved pseudotachylytes of variable thickness. We determined the geometry of fault traces over almost five orders of magnitude using terrestrial laser-scanning (LIDAR, ca. 500 to <1 m scale), and 3D mosaics of high-resolution rectified digital photographs (10 m to ca. 1 mm scale). LIDAR scans and photomosaics were georeferenced in 3D using a Differential Global Positioning System, allowing detailed multiscale reconstruction of fault traces in Gocad^®. The combination of LIDAR and high-resolution photos has the advantage, compared with classical LIDAR-only surveys, that the spatial resolution of rectified photographs can be very high (up to 0.2 mm/pixel in this study), allowing for detailed outcrop characterization. Fourier power spectrum analysis of the fault traces revealed a self-affine behaviour over 3–5 orders of magnitude, with Hurst exponents ranging between 0.6 and 0.8. Parameters from Fourier analysis have been used to reconstruct synthetic 3D fault surfaces with an equivalent roughness by means of 2D Fourier synthesis. Roughness of pre-existing joints is in a typical range for this kind of structure. Roughness of faults at small scale (1 m to 1 mm) shows a clear genetic relationship with the roughness of precursor joints, and some anisotropy in the self-affine Hurst exponent. Roughness of faults at scales larger than net slip (>1–10 m) is not anisotropic and less evolved than at smaller scales. These observations are consistent with an evolution of roughness, due to fault surface processes, that takes place only at scales smaller or comparable to the observed net slip. Differences in roughness evolution between shallow and deeper faults, the latter showing evidences of seismic activity, are interpreted as the result of different weakening versus induration processes, which also result in localization versus delocalization of deformation in the fault zone. From a methodological point of view, the technique used here is advantageous over direct measurements of exposed fault surfaces in that it preserves, in cross-section, all of the structures which contribute to fault roughness, and removes any subjectivity introduced by the need to distinguish roughness of original slip surfaces from roughness induced by secondary weathering processes. Moreover, offsets can be measured by means of suitable markers and fault rocks are preserved, hence their thickness, composition and structural features can be characterised, providing an integrated dataset which sheds new light on mechanisms of roughness evolution with slip and concomitant fault rock production.     

Assessing Nitrogen Dynamics Throughout the Estuarine Landscape

Assessing nitrogen dynamics in the estuarine landscape is challenging given the unique effects of individual habitats on nitrogen dynamics. We measured net N₂ fluxes, sediment oxygen demand, and fluxes of ammonium and nitrate seasonally from five major estuarine habitats: salt marshes, seagrass beds (SAV), oyster reefs, and intertidal and subtidal flats. Net N₂ fluxes ranged from 332?±?116 μmol?N-N₂?m^?2?h^?1 from oyster reef sediments in the summer to ?67?±?4 μmol?N-N₂?m^?2?h^?1 from SAV in the winter. Oyster reef sediments had the highest rate of N₂ production of all habitats. Dissimilatory nitrate reduction to ammonium (DNRA) was measured during the summer and winter. DNRA was low during the winter and ranged from 4.5?±?3.0 in subtidal flats to 104?±?34 μmol?¹⁵NH ₄ ⁺ ?m^?2?h^?1 in oyster reefs during the summer. Annual denitrification, accounting for seasonal differences in inundation and light, ranged from 161.1?±?19.2 mmol?N-N₂?m^?2?year^?1 for marsh sediments to 509.9?±?122.7 mmol?N-N₂?m^?2?year^?1 for SAV sediments. Given the current habitat distribution in our study system, an estimated 28.3?×?10⁶?mol of N are removed per year or 76 % of estimated watershed nitrogen load. These results indicate that changes in the area and distribution of habitats in the estuarine landscape will impact ecosystem function and services.     

The application of lead isotope ratios in the Antarctic macroalga Iridaea cordata as a contaminant monitoring tool

Lead (Pb) isotope ratios were measured in the marine macroalga Iridaea cordata collected from four locations in the Windmill Islands, East Antarctica. Based on the masses of thalli collected, samples analysed in this study were likely to be a mixture of one and two year old thalli. For a sample of thalli of various ages (<12 months to 2 years old) from the same site there was no apparent variation in Pb concentration or Pb isotope ratio with thallus mass/age, indicating that contaminant sources had been constant over the lifetime of the thalli sampled. I.cordata samples close to the Thala Valley waste disposal site (Brown Bay Inner) near the Australian Station, Casey, displayed isotopic signatures (²⁰⁸Pb/²⁰⁴Pb 35.99; ²⁰⁶Pb/²⁰⁷Pb 1.066; n = 3; average values shown) trending towards that possessed by major Australian Pb sources (Broken Hill and Mt Isa, ²⁰⁸Pb/²⁰⁴Pb 35.60; ²⁰⁶Pb/²⁰⁷Pb 1.041) suggesting that these samples had been exposed to anthropogenic Pb originating from the Thala Valley waste disposal site. Material collected hundreds of metres from the tip location at Brown Bay Outer had isotopic values (²⁰⁸Pb/²⁰⁴Pb 36.32; ²⁰⁶Pb/²⁰⁷Pb 1.088; n = 10) intermediate between Brown Bay Inner and sites further from the contaminant source at Sparkes Bay and Wilkes (²⁰⁸Pb/²⁰⁴Pb 36.46; ²⁰⁶Pb/²⁰⁷Pb 1.094; n = 4) showing that contaminant transport was predominantly restricted to Brown Bay Inner. This study demonstrates that the isotope ratios of Pb in marine macroalgae can provide valuable information as to the origin and extent of heavy metal flux in a marine environment.     

Influence of Outcrop Scale Fractures on the Effective Stiffness of Fault Damage Zone Rocks

We combine detailed mapping and microstructural analyses of small fault zones in granodiorite with numerical mechanical models to estimate the effect of mesoscopic (outcrop-scale) damage zone fractures on the effective stiffness of the fault zone rocks. The Bear Creek fault zones were active at depths between 4 and 15 km and localize mesoscopic off-fault damage into tabular zones between two subparallel boundary faults, producing a fracture-induced material contrast across the boundary faults with softer rocks between the boundary faults and intact granodiorite outside the boundary faults. Using digitized fault zone fracture maps as the modeled fault geometries, we conduct nonlinear uniaxial compression tests using a novel finite-element method code as the experimental “laboratory” apparatus. Map measurements show that the fault zones have high nondimensional facture densities (>1), and damage zone fractures anastamose and intersect, making existing analytical effective medium models inadequate for estimation of the effective elastic properties. Numerical experiments show that the damage zone is strongly anisotropic and the bulk response of the fault zone is strain-weakening. Normal strains as small as 2% can induce a reduction of the overall stiffness of up to 75%. Fracture-induced effective stiffness changes are large enough to locally be greater than intact modulus changes across the fault due to juxtaposition of rocks of different lithologies; therefore mesoscopic fracturing is as important as rock type when considering material or bimaterial effects on earthquake mechanics. These results have important implications for earthquake rupture mechanics models, because mesoscopic damage zone fractures can cause a material contrast across the faults as large as any lithology-based material contrast at seismogenic depths, and the effective moduli can be highly variable during a single rupture event.     

Harvest patterns and effort dynamics of indigenous and non-indigenous commercial sectors of the eastern Torres Strait reef line fishery

The reef line fishery (RLF) in eastern Torres Strait (ETS) is unique in that it has both a commercial indigenous sector and a commercial non-indigenous sector. Recently, concerns have been expressed by all stakeholders about the long-term sustainability of the fishery. These concerns have been exacerbated by the lack of detailed catch and effort information from both sectors, which has precluded any formal assessment of the fishery. In this paper, we characterise the harvest patterns and effort dynamics of the indigenous and non-indigenous commercial sectors of the ETS RLF using a range of data sources including commercial logbooks, community freezer records, voluntary logbooks and observer surveys. We demonstrate that bycatch is a significant component of the catch for both sectors and identify substantial differences in harvest patterns and effort dynamics between the sectors. Differences between sectors were observed in species composition and spatial and temporal patterns in catch, effort and catch per unit effort. These results highlight the inherent variation in catch and effort dynamics between the two commercial sectors of the ETS RLF and provide valuable information for the development of future assessments and appropriate management strategies for the fishery. The more reliable estimates of harvest patterns and effort dynamics for both sectors obtained from observer surveys will also assist in resolving issues relating to allocation of reef fish resources in Torres Strait.     

Determination of High Field Strength Elements, Rb, Sr, Mo, Sb, Cs, Tl and Bi at ng g−1 Levels in Geological Reference Materials by Magnetic Sector ICP-MS after HF/HClO4 High Pressure Digestion

Six low abundance rock reference materials (basalt BIR-1, dunite DTS-1, dolerite DNC-1, peridotite PCC-1, serpentine UB-N and basalt TAFAHI) have been analysed for high field strength elements (Zr, Nb, Hf, Ta, Th and U), Rb, Sr, Mo, Sb, Cs, Tl and Bi at ng g⁻¹ levels (in rock) by magnetic sector inductively coupled plasma-mass spectrometry after HF/HClO₄ high pressure decomposition. The adopted method uses only indium as an internal standard. Detection limits were found to be in the range of 0.08 to 16.2 pg ml⁻¹ in solution (equivalent to 0.08 to 16.2 ng g⁻¹ in rock). Our data for high field strength elements, Rb, Sr, Mo, Sb, Cs, Tl and Bi for the six selected low abundance geological reference materials show general agreement with previously published data. Our Ta values in DTS-1 and PCC-1 (1.3 and 0.5 ng g⁻¹) are lower than in previously published studies, providing smooth primitive mantle distribution patterns. Lower values were also found for Tl in BIR-1, DTS-1 and PCC-1 (2, 0.4 and 0.8 ng g⁻¹). Compared with quadrupole ICP-MS studies, the proposed magnetic sector ICP-MS method can generally provide better detection limits, so that the measurement of high field strength elements, Rb, Sr, Mo, Sb, Cs, Tl and Bi at ng g⁻¹ levels can be achieved without pre-concentration, ion exchange separation or other specialised techniques.     

The microscale cooling effects of water sensitive urban design and irrigation in a suburban environment

There is still considerable uncertainty about precipitation at high elevation in mountain terrain due to the relatively few in situ measurements available and to the particular variability of the parameter. In this study, several spatialization techniques were tested, some for climatological time scale and others for daily fields, for precipitation over the western Alps for the period of 1990–2012. The study domain and period were chosen for the quality of available in situ observations and density of the network. First, a weather-type classification was established with a technique based on canonical correlation analysis combining large- and regional-scale data. The spatialization techniques applied for the climatological time scale were adapted from the Aurelhy method which uses elevation and principal components of the topography as predictors. The spatialization techniques applied to daily fields were based on kriging of daily rain gauges and used the climatological fields as predictors. This study aims to validate the advantage of using the climatology of the weather type of the day as predictor for daily fields over a monthly climatology. The climatology of the weather type of the day seems to demonstrate some small improvement.Finally, annual means over the period of 1990–2012 were produced using several methods, including some from accumulation of daily fields and others from the spatialization of in situ station means. Precipitation at high elevations and vertical climatological gradients were particularly scrutinized. Annual means based on sums of daily fields seem to have better performances.This paper only presents results for precipitation but temperature was also analysed.