Spatiotemporal analysis of residential flood exposure in the Atlanta,Georgia metropolitan area

Natural Hazards - This research examines changes in residential built-environment flood exposure within the current boundaries of the Atlanta, Georgia metropolitan statistical area, by estimating...     

Diachronous isotopic and sedimentary responses to topographic change as indicators of mid‐Eocene hydrologic reorganization in the western United States

Early Cenozoic terrestrial deposits in the western United States represent well‐preserved archives of climatic and tectonic processes that together shaped the Earth's surface during the demise of a large continental plateau. This study examines a Cenozoic terrestrial sedimentary sequence in the central part of the Cordilleran orogen (Montana) using sedimentologic and geochemical techniques. At ～49 Ma, we observe rapid major shifts in oxygen, carbon and strontium isotope records that are too large to directly reflect changes in meteoric water composition due to simple orographic rainout. The transition to low‐δ¹⁸O values in pedogenic carbonate in concert with changes in the composition of clastic material at ～49 Ma points to the input of evolved meteoric water to the hydrological cycle due to a change in the source of waters reaching Cordilleran intermontane regions in southwestern Montana. This drainage reorganization coincides with the initiation of magmatism and extension to the west in what is now Montana and Idaho. The sedimentological record shows evidence that depositional gradients increased in the study area ～46 Ma, ～3 Myr after the drainage reorganization occurred. This interval is most likely the time it took for extensional deformation to propagate to the study area itself. Evidence of freshening events in Laramide Basins to the southeast suggests that this drainage reorganization diverted waters to progressively fill these basins and highlights the impact of post‐plateau extension‐related landscape reorganization on river networks and lake dynamics. This study also emphasizes the importance of using multiple tools in deciphering topographic history through the study of terrestrial basin deposits, in that interpretation based on any single method employed would have compromised our ability to successfully identify the regional evolution of topography and drainage networks.     

A numerical study of detached shear layers in a rotating sliced cylinder

Abstract

The low Rossby number flow in a rotating cylinder with an inclined bottom, of small slope, is examined when part of the lid of the container is rotating at a slightly different rate. The resulting flow is calculated numerically by solving the governing equations for the two-dimensional geostrophic motion which approximates the flow in most of the fluid including the inertially-modified E ^¼ -layers. The presence of ageostrophic regions, on the container walls and beneath the velocity discontinuity on the lid, is accounted for in the governing equations and their boundary conditions. This study supplements previous work on this configuration, in which the zero Rossby number flow was calculated and experimental results were presented, by enabling a direct comparison to be made between the results of the low Rossby number theory and the experiments. The numerical results for a range of Rossby and Ekman numbers compare well with those from the experiments despite a severe limitation on the size of the Rossby number arising from the analysis in the ageostrophic part of the detached shear layer.     

Kinematic dynamo action in a sphere: Effects of periodic time-dependent flows on solutions with axial dipole symmetry

Choosing a simple class of flows, with characteristics that may be present in the Earth's core, we study the ability to generate a magnetic field when the flow is permitted to oscillate periodically in time. The flow characteristics are parameterised by D, representing a differential rotation, M, a meridional circulation, and C, a component characterising convective rolls. The dynamo action of all solutions with fixed parameters (steady flows) is known from earlier studies. Dynamo action is sensitive to these flow parameters and fails spectacularly for much of the parameter space where magnetic flux is concentrated into small regions, leading to high diffusion. In addition, steady flows generate only steady or regularly reversing oscillatory fields and cannot therefore reproduce irregular geomagnetic-type reversal behaviour. Oscillations of the flow are introduced by varying the flow parameters in time, defining a closed orbit in the space ( D,?M ). When the frequency of the oscillation is small, the net growth rate of the magnetic field over one period approaches the average of the growth rates for steady flows along the orbit. At increased frequency time-dependence appears to smooth out flux concentrations, often enhancing dynamo action. Dynamo action can be impaired, however, when flux concentrations of opposite signs occur close together as smoothing destroys the flux by cancellation. It is possible to produce geomagnetic-type reversals by making the orbit stray into a region where the steady flows generate oscillatory fields. In this case, however, dynamo action was not found to be enhanced by the time-dependence. A novel approach is being taken to solve the time-dependent eigenvalue problem where, by combining Floquet theory with a matrix-free Krylov-subspace method, we can avoid large memory requirements for storing the matrix required by the standard approach.     

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.     

Marine Macrophyte Wrack Inputs and Dissolved Nutrients in Beach Sands

We investigated the role of sandy beaches in nearshore nutrient cycling by quantifying macrophyte wrack inputs and examining relationships between wrack accumulation and pore water nutrients during the summer dry season. Macrophyte inputs, primarily giant kelp Macrocystis pyrifera, exceeded 2.3 kg m⁻¹ day⁻¹. Mean wrack biomass varied 100-fold among beaches (range = 0.41 to 46.43 kg m⁻¹). Mean concentrations of dissolved inorganic nitrogen (DIN), primarily NO_x⁻-N, and dissolved organic nitrogen (DON) in intertidal pore water varied significantly among beaches (ranges = 1 to 6,553 μM and 7 to 2,006 μM, respectively). Intertidal DIN and DON concentrations were significantly correlated with wrack biomass. Surf zone concentrations of DIN were also strongly correlated with wrack biomass and with intertidal DIN, suggesting export of nutrients from re-mineralized wrack. Our results suggest beach ecosystems can process and re-mineralize substantial organic inputs and accumulate dissolved nutrients, which are subsequently available to nearshore waters and primary producers.     

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.