Kinematic analysis of mélange fabrics: examples and applications from the McHugh Complex, Kenai Peninsula, Alaska

Permian to Cretaceous mélange of the McHugh Complex on the Kenai Peninsula, south-central Alaska includes blocks and belts of graywacke, argillite, limestone, chert, basalt, gabbro, and ultramafic rocks, intruded by a variety of igneous rocks. An oceanic plate stratigraphy is repeated hundreds of times across the map area, but most structures at the outcrop scale extend lithological layering. Strong rheological units occur as blocks within a matrix that flowed around the competent blocks during deformation, forming broken formation and mélange. Deformation was noncoaxial, and disruption of primary layering was a consequence of general strain driven by plate convergence in a relatively narrow zone between the overriding accretionary wedge and the downgoing, generally thinly sedimented oceanic plate. Soft-sediment deformation processes do not appear to have played a major role in the formation of the mélange. A model for deformation at the toe of the wedge is proposed in which layers oriented at low angles to σ₁ are contracted in both the brittle and ductile regimes, layers at 30–45° to σ₁ are extended in the brittle regime and contracted in the ductile regime, and layers at angles greater than 45° to σ₁ are extended in both the brittle and ductile regimes. Imbrication in thrust duplexes occurs at deeper levels within the wedge. Many structures within mélange of the McHugh Complex are asymmetric and record kinematic information consistent with the inferred structural setting in an accretionary wedge. A displacement field for the McHugh Complex on the lower Kenai Peninsula includes three belts: an inboard belt of Late Triassic rocks records west-to-east-directed slip of hanging walls, a central belt of predominantly Early Jurassic rocks records north–south directed displacements, and Early Cretaceous rocks in an outboard belt preserve southwest–northeast directed slip vectors. Although precise ages of accretion are unknown, slip directions are compatible with inferred plate motions during the general time frame of accretion of the McHugh Complex. The slip vectors are interpreted to preserve the convergence directions between the overriding and underriding plates, which became more oblique with time. They are not considered indicative of strain partitioning into belts of orogen-parallel and orogen-perpendicular displacements, because the kinematic data are derived from the earliest preserved structures, whereas fabrics related to strain partitioning would be expected to be superimposed on earlier accretion-related fabrics.     

Strong slope-aspect control of regolith thickness by bedrock foliation

The porous near-surface layer of the Earth's crust – the critical zone – constitutes a vital reservoir of water for ecosystems, provides baseflow to streams, guides recharge to deep aquifers, filters contaminants from groundwater, and regulates the long-term evolution of landscapes. Recent work suggests that the controls on regolith thickness include climate, tectonics, lithology, and vegetation. However, the relative paucity of observations of regolith structure and properties at landscape scales means that theoretical models of critical zone structure are incompletely tested. Here we present seismic refraction and electrical resistivity surveys that thoroughly characterize subsurface structure in a small catchment in the Santa Catalina Mountains, Arizona, USA, where slope-aspect effects on regolith structure are expected based on differences in vegetation. Our results show a stark contrast in physical properties and inferred regolith thickness on opposing slopes, but in the opposite sense of that expected from environmental models and observed vegetation patterns. Although vegetation (as expressed by normalized difference vegetation index [NDVI]) is denser on the north-facing slope, regolith on the south-facing slope is four times thicker (as indicated by lower seismic velocities and resistivities). This contrast cannot be explained by variations in topographic stress or conventional hillslope morphology models. Instead, regolith thickness appears to be controlled by metamorphic foliation: regolith is thicker where foliation dips into the topography, and thinner where foliation is nearly parallel to the surface. We hypothesize that, in this catchment, hydraulic conductivity and infiltration capacity control weathering: infiltration is hindered and regolith is thin where foliation is parallel to the surface topography, whereas water infiltrates deeper and regolith is thicker where foliation intersects topography at a substantial angle. These results suggest that bedrock foliation, and perhaps by extension sedimentary layering, can control regolith thickness and must be accounted for in models of critical zone development. © 2020 John Wiley & Sons, Ltd.     

Volcano monitoring in NZ and links to SW Pacific via the Wellington VAAC

Monitoring for natural hazards like earthquakes, volcanic eruptions and landslides in New Zealand is under taken by GNS Science through the GeoNet project, funded by the Earthquake Commission (EQC). The volcano monitoring function, which includes seismograph and GPS networks, visual observations, webcams, deformation, gas and chemical measurements, assesses the status of New Zealand’s 12 active volcanoes and disseminates warnings as necessary. Interaction with the Wellington VAAC addresses aviation concerns. Near-real-time data is obtained from the monitoring networks and processed via two data centres, with two Duty Officers who respond to pager alerts. Through the National Civil Defence Plan, the volcano monitoring function is delegated to GNS Science, which is responsible for setting the volcanic alert level at each volcano. The Alert level is then used by responding agencies, industries, utility providers and the public to set their response. The New Zealand alert level system is based on the current status of the volcano, and is not necessarily predictive. The aviation industry is one of these responding industries. As a consequence of the 1995 and 1996 eruptions of Ruapehu volcano, it became apparent that the aviation industry required more information than the standard VAAC data and an additional service known as VAAS (NZ Volcanic Ash Advisory System) was established by MetService NZ and GNS Science. This arrangement provides enhanced services to the airlines and ensures rapid dissemination and quantification of information like SIGMET and NOTAM. GNS Science has informal relationships with similar geological based organisations in SW Pacific countries, especially Vanuatu and the Solomon Islands, and is able to quantify information about eruptions in those countries. Similarly, MetService NZ has relationships and is able to obtain pilot reports and updates from the meteorological and air traffic control organisations.     

Dissolved oxygen as an indicator of bioavailable dissolved organic carbon in groundwater

Concentrations of dissolved oxygen (DO) plotted vs. dissolved organic carbon (DOC) in groundwater samples taken from a coastal plain aquifer of South Carolina (SC) showed a statistically significant hyperbolic relationship. In contrast, DO-DOC plots of groundwater samples taken from the eastern San Joaquin Valley of California (CA) showed a random scatter. It was hypothesized that differences in the bioavailability of naturally occurring DOC might contribute to these observations. This hypothesis was examined by comparing nine different biochemical indicators of DOC bioavailability in groundwater sampled from these two systems. Concentrations of DOC, total hydrolysable neutral sugars (THNS), total hydrolysable amino acids (THAA), mole% glycine of THAA, initial bacterial cell counts, bacterial growth rates, and carbon dioxide production/consumption were greater in SC samples relative to CA samples. In contrast, the mole% glucose of THNS and the aromaticity (SUVA(254)) of DOC was greater in CA samples. Each of these indicator parameters were observed to change with depth in the SC system in a manner consistent with active biodegradation. These results are uniformly consistent with the hypothesis that the bioavailability of DOC is greater in SC relative to CA groundwater samples. This, in turn, suggests that the presence/absence of a hyperbolic DO-DOC relationship may be a qualitative indicator of relative DOC bioavailability in groundwater systems.     

MTBE, TBA, and TAME Attenuation in Diverse Hyporheic Zones

Groundwater contamination by fuel-related compounds such as the fuel oxygenates methyl tert -butyl ether (MTBE), tert -butyl alcohol (TBA), and tert -amyl methyl ether (TAME) presents a significant issue to managers and consumers of groundwater and surface water that receives groundwater discharge. Four sites were investigated on Long Island, New York, characterized by groundwater contaminated with gasoline and fuel oxygenates that ultimately discharge to fresh, brackish, or saline surface water. For each site, contaminated groundwater discharge zones were delineated using pore water geochemistry data from 15 feet (4.5 m) beneath the bottom of the surface water body in the hyporheic zone and seepage-meter tests were conducted to measure discharge rates. These data when combined indicate that MTBE, TBA, and TAME concentrations in groundwater discharge in a 5-foot (1.5-m) thick section of the hyporheic zone were attenuated between 34% and 95%, in contrast to immeasurable attenuation in the shallow aquifer during contaminant transport between 0.1 and 1.5 miles (0.1 to 2.4 km). The attenuation observed in the hyporheic zone occurred primarily by physical processes such as mixing of groundwater and surface water. Biodegradation also occurred as confirmed in laboratory microcosms by the mineralization of U- ¹⁴C-MTBE and U-¹⁴C-TBA to ¹⁴CO₂ and the novel biodegradation of U- ¹⁴C-TAME to ¹⁴CO₂ under oxic and anoxic conditions. The implication of fuel oxygenate attenuation observed in diverse hyporheic zones suggests an assessment of the hyporheic zone attenuation potential (HZAP) merits inclusion as part of site assessment strategies associated with monitored or engineered attenuation.     

Component correlations in structure‐specific seismic loss estimation

This paper addresses correlations between multiple components in structure‐specific seismic loss estimation. To date, the consideration of such correlations has been limited by methodological tractability, increased computational demand, and a paucity of data for their computation. The effect of component correlations, which arises in various forms, is however a significant factor affecting the results of structure‐specific seismic loss estimation and therefore it is prudent that adequate consideration be given to their effect. This paper provides the details of a tractable and computationally efficient seismic loss estimation methodology in which correlations can be considered. Methods to determine the necessary correlations are discussed, particularly those that can be used in the absence of sufficient empirical data, for which values are suggested based on the judgement. The effects of various assumptions regarding correlations are illustrated via application to a case‐study office structure. It is observed that certain correlation assumptions can lead to errors in excess of 50% in the lognormal standard deviation in the loss given intensity and loss hazard relationships, while full consideration of partial correlations is 50 times more computationally expensive than other assumptions. Copyright © 2009 John Wiley & Sons, Ltd.     

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.     

Multidimensional upstream weighting for multiphase transport on general grids

The governing equations for multiphase flow in porous media have a mixed character, with both nearly elliptic and nearly hyperbolic variables. The flux for each phase can be decomposed into two parts: (1) a geometry- and rock-dependent term that resembles a single-phase flux; and (2) a mobility term representing fluid properties and rock–fluid interactions. The first term is commonly discretized by two- or multipoint flux approximations (TPFA and MPFA, respectively). The mobility is usually treated with single-point upstream weighting (SPU), also known as dimensional or donor cell upstream weighting. It is well known that when simulating processes with adverse mobility ratios, SPU suffers from grid orientation effects. An important example of this, which will be considered in this work, is the displacement of a heavy oil by water. For these adverse mobility ratio flows, the governing equations are unstable at the modeling scale, rendering a challenging numerical problem. These challenges must be addressed in order to avoid systematic biasing of simulation results. In this work, we present a framework for multidimensional upstream weighting for multiphase flow with buoyancy on general two-dimensional grids. The methodology is based on a dual grid, and the resulting transport methods are provably monotone. The multidimensional transport methods are coupled with MPFA methods to solve the pressure equation. Both explicit and fully implicit approaches are considered for time integration of the transport equations. The results show considerable reduction of grid orientation effects compared to SPU, and the explicit multidimensional approach allows larger time steps. For the implicit method, the total number of non-linear iterations is also reduced when multidimensional upstream weighting is used.     

Submarine Groundwater Discharge-Derived Nutrient Loads to San Francisco Bay: Implications to Future Ecosystem Changes

Submarine groundwater discharge (SGD) was quantified at select sites in San Francisco Bay (SFB) from radium (²²³Ra and ²²⁴Ra) and radon (²²²Rn) activities measured in groundwater and surface water using simple mass balance box models. Based on these models, discharge rates in South and Central Bays were 0.3?C7.4?m³?day^?1?m^?1. Although SGD fluxes at the two regions (Central and South Bays) of SFB were of the same order of magnitude, the dissolved inorganic nitrogen (DIN) species associated with SGD were different. In the South Bay, ammonium (NH ₄ ⁺ ) concentrations in groundwater were three-fold higher than in open bay waters, and NH ₄ ⁺ was the primary DIN form discharged by SGD. At the Central Bay site, the primary DIN form in groundwater and associated discharge was nitrate (NO ₃ ^? ). The stable isotope signatures (??¹⁵N_NO3 and ??¹⁸O_NO3) of NO ₃ ^? in the South Bay groundwater and surface waters were both consistent with NO ₃ ^? derived from NH ₄ ⁺ that was isotopically enriched in ¹⁵N by NH ₄ ⁺ volatilization. Based on the calculated SGD fluxes and groundwater nutrient concentrations, nutrient fluxes associated with SGD can account for up to 16?% of DIN and 22?% of DIP in South and Central Bays. The form of DIN contributed to surface waters from SGD may impact the ratio of NO ₃ ^? to NH ₄ ⁺ available to phytoplankton with implications to bay productivity, phytoplankton species distribution, and nutrient uptake rates. This assessment of nutrient delivery via groundwater discharge in SFB may provide vital information for future bay ecological wellbeing and sensitivity to future environmental stressors.     

TAZ-level variation in work trip mode choice between 1990 and 2000 and the presence of rail transit

Rail transit continues to be a popular alternative for cities as a tool for alleviating automobile congestion and for redeveloping areas into transit and pedestrian-friendly environments. Ideally, rail transit will draw trips away from cars, but the quantitative work that tests this notion has often been either case studies of neighborhoods, in which conclusions are tough to generalize, or citywide comparisons where important spatial variation is lost in aggregation. This study seeks to narrow this gap in the research by using multivariate analysis of covariance to isolate the effect of covariates and cities on changes in work trip mode choice at the traffic analysis zone (TAZ) level for nine cities between 1990 and 2000. The results suggest differences by city in the change at the TAZ level of the proportion of people driving alone and taking transit. Increases in transit usage were associated with cities that built rail transit, while increases in automobile commuting and decreases in transit usage were associated with cities that did not.