Interactive Visualization to Advance Earthquake Simulation

The geological sciences are challenged to manage and interpret increasing volumes of data as observations and simulations increase in size and complexity. For example, simulations of earthquake-related processes typically generate complex, time-varying data sets in two or more dimensions. To facilitate interpretation and analysis of these data sets, evaluate the underlying models, and to drive future calculations, we have developed methods of interactive visualization with a special focus on using immersive virtual reality (VR) environments to interact with models of Earth’s surface and interior. Virtual mapping tools allow virtual “field studies” in inaccessible regions. Interactive tools allow us to manipulate shapes in order to construct models of geological features for geodynamic models, while feature extraction tools support quantitative measurement of structures that emerge from numerical simulation or field observations, thereby enabling us to improve our interpretation of the dynamical processes that drive earthquakes. VR has traditionally been used primarily as a presentation tool, albeit with active navigation through data. Reaping the full intellectual benefits of immersive VR as a tool for scientific analysis requires building on the method’s strengths, that is, using both 3D perception and interaction with observed or simulated data. This approach also takes advantage of the specialized skills of geological scientists who are trained to interpret, the often limited, geological and geophysical data available from field observations.     

Surficial fractures in the Navajo sandstone,south-western USA: the roles of thermal cycles,rainstorms, granular disintegration,and iterative cracking

Deep (> 5 m) sheeting fractures in the Navajo sandstone are evident at numerous sites in southern Utah and derive from tectonic stresses. Strong diurnal thermal cycles are, however, the likely triggers for shallow (< 0.3 m) sheeting fractures. Data from subsurface thermal sensors reveal that large temperature differences between sensors at 2 and 15 cm depth on clear summer afternoons are as great as those that trigger sheeting fractures in exposed California granite. Extensive polygonal patterns in the Navajo sandstone are composed of surface-perpendicular fractures and were produced by contractile stresses. Numerous studies have shown that porewater diminishes the tensile strength of sandstone. Based on our thermal records, we propose that cooling during monsoonal rainstorms triggers polygonal fracturing of temporarily weakened rock. On steep outcrops, polygonal patterns are rectilinear and orthogonal, with T-vertices. Lower-angle slopes host hexagonal patterns (defined by the dominance of Y-vertices). Intermediate patterns with rectangles and hexagons of similar scale are common. We posit that outcropping fractures are advancing downward by iterative steps, and that hexagons on sandstone surfaces (like prismatic columns of basalt) have evolved from ancestral orthogonal polygons of similar scale. In lava flows, fractures elongate intermittently as they follow a steep thermal gradient (the source of stress) as it rapidly moves through the rock mass. In our model, a steep, surficial thermal gradient descends through unfractured sandstone, but at the slow pace of granular disintegration. Through time, as the friable rock on stable slopes erodes, iterative cracking advances into new space. Hexagonal patterns form as new fractures, imperfectly guided by the older ones, propagate in new directions, and vertices drift into a configuration that minimizes the ratio of fracture length to polygon area. © 2020 John Wiley & Sons, Ltd.     

The use of monitoring data for identifying factors influencing phytoplankton bloom dynamics in the eutrophic Taw Estuary, SW England

Using the Taw Estuary as an example, data routinely collected by the Environment Agency for England and Wales over the period 1990-2004 were interrogated to identify the drivers of excessive algal growth. The estuary was highly productive with chlorophyll concentrations regularly exceeding 100 μg L⁻¹, mostly during periods of low freshwater input from the River Taw when estuarine water residence times were longest. However, algal growth in mid estuary was often inhibited by ammonia inputs from the adjacent sewage treatment works. The reported approach demonstrates the value of applying conventional statistical analyses in a structured way to existing monitoring data and is recommended as a useful tool for the rapid assessment of eutrophication. However, future estuarine monitoring should include the collection of dissolved organic nutrient data and targeted high temporal resolution data because the drivers of eutrophication are complex and often very specific to a particular estuary.     

Spatial Patterns of Sea Level Variability Associated with Natural Internal Climate Modes

Sea level rise (SLR) can exert significant stress on highly populated coastal societies and low-lying island countries around the world. Because of this, there is huge societal demand for improved decadal predictions and future projections of SLR, particularly on a local scale along coastlines. Regionally, sea level variations can deviate considerably from the global mean due to various geophysical processes. These include changes of ocean circulations, which partially can be attributed to natural, internal modes of variability in the complex Earth’s climate system. Anthropogenic influence may also contribute to regional sea level variations. Separating the effects of natural climate modes and anthropogenic forcing, however, remains a challenge and requires identification of the imprint of specific climate modes in observed sea level change patterns. In this paper, we review our current state of knowledge about spatial patterns of sea level variability associated with natural climate modes on interannual-to-multidecadal timescales, with particular focus on decadal-to-multidecadal variability. Relevant climate modes and our current state of understanding their associated sea level patterns and driving mechanisms are elaborated separately for the Pacific, the Indian, the Atlantic, and the Arctic and Southern Oceans. We also discuss the issues, challenges and future outlooks for understanding the regional sea level patterns associated with climate modes. Effects of these internal modes have to be taken into account in order to achieve more reliable near-term predictions and future projections of regional SLR.     

Triaxial induction tool with electrode sleeve: FD modeling in 3D geometries

We have developed a finite-difference (FD) method to model the response of a multiple-spacing triaxial array induction tool with a multiple-electrode type sleeve. The FD software was verified versus two other independent modeling methods. The model response was checked for the FD grid refinement, while attempting to minimize run time. Many details of the electrode-sleeve geometry were studied using this method that led to the final electrode-sleeve design. The electrode sleeve reduces the response of the transverse couplings to eccentering in a conductive borehole. Quantitative agreement was found between the measurements and calculations. Using the FD software we study the tool response to various effects of 3D geometry.     

Sulfur chemistry in the Venus mesosphere from SO2 and SO microwave spectra

First measurements of SO₂ and SO in the Venus mesosphere (70-100 km) are reported. This altitude range is distinctly above the ∼60-70 km range to which nadir-sounding IR and UV investigations are sensitive. Since July 2004, use of ground-based sub-mm spectroscopy has yielded multiple discoveries. Abundance of each molecule varies strongly on many timescales over the entire sub-Earth Venus hemisphere. Diurnal behavior is evident, with more SO₂, and less SO, at night than during the day. Non-diurnal variability is also present, with measured SO₂ and SO abundances each changing by up to 2× or more between observations conducted on different dates, but at fixed phase, hence identical sub-Earth Venus local times. Change as large and rapid as a 5σ doubling of SO on a one-week timescale is seen. The sum of SO₂ and SO abundances varies by an order of magnitude or more, indicating at least one additional sulfur reservoir must be present, and that it must function as both a sink and source for these molecules. The ratio SO₂/SO varies by nearly two orders of magnitude, with both diurnal and non-diurnal components. In contrast to the strong time dependence of molecular abundances, their altitude distributions are temporally invariant, with far more SO₂ and SO at 85-100 km than at 70-85 km. The observed increase of SO₂ mixing ratio with altitude requires that the primary SO₂ source be upper mesospheric photochemistry, contrary to atmospheric models which assert upward transport as the only source of above-cloud SO₂. Abundance of upper mesospheric aerosol, with assumption that it is composed primarily of sulfuric acid, is at least sufficient to provide the maximum gas phase (SO + SO₂) sulfur reported in this study. Sulfate aerosol is thus a plausible source of upper mesospheric SO₂.     

EQUILIBRIUM AND NONEQUILIBRIUM CONDITIONS IN DRYLAND RIVERS

Rivers in drylands typically are characterized by extreme flow variability, with long periods of little or no flow interspersed with occasional large, sometimes extreme, floods. Complete adjustment of river form and process is sometimes inhibited, resulting in a common assumption that equilibrium conditions may rarely, if ever, exist in dryland rivers, and that transient and unstable (nonequilibrium) behavior is the norm. Examples from the Channel Country and the Northern Plains in central Australia challenge that notion. Along the middle reaches of these intermediate and large, low-gradient rivers, where long duration floods generate moderate to low unit stream powers and boundary resistance is high as a result of indurated alluvial terraces, cohesive muds or riparian vegetation, there is evidence that: (1) channels have remained essentially stable despite large floods; (2) sediment transport discontinuities, while present at a catchment scale, are largely insignificant for channel form and process in individual reaches; (3) there are strong correlations between many channel form and process variables; and (4) many rivers appear to be adjusted to maximum sediment transport efficiency under conditions of low gradient, abundant within-channel vegetation and declining downstream discharge. In these middle reaches, rivers are characterized by equilibrium conditions. However, in the aggradational lower reaches of rivers on the Northern Plains, where upstream terraces are buried by younger sediments and channels are less confined, nonequilibrium conditions prevail. Here, channels sometimes undergo sudden and substantial changes in form during large floods, sediment transport discontinuities are readily apparent, and landforms such as splays remain out-of-balance with normal flows. Hence, dryland rivers can exhibit both equilibrium and nonequilibrium conditions, depending on factors such as catchment size, channel gradient, flood duration, unit stream power, channel confinement, sediment cohesion, and bank strength. [Key words: dryland rivers, floods, equilibrium, nonequilibrium, central Australia.]     

Inbank and overbank velocity conditions in an arid zone anastomosing river

Relative to those at sub‐bankfull flow, hydraulic conditions at overbank flow, whether in the channel or on the floodplain, are poorly understood. Here, velocity conditions are analysed over an unusually wide range of flows in the arid zone river of Cooper Creek with its complex system of anastomosing channels and large fluctuations in floodplain width. At‐a‐station hydraulic geometry relationships reveal sharp discontinuities in velocity at the inbank–overbank transition, the nature of the discontinuity varying with the degree of flow confinement and the level of channel–floodplain interaction. However, despite inter‐sectional differences, velocities remain modest throughout the flow range in this low‐gradient river, and the large increases in at‐a‐station discharge are principally accommodated by changes in cross‐sectional area. Velocity distribution plots suggest that within‐channel conditions during overbank flow are characterized by a central band of high velocity which penetrates far toward the bed, helping to maintain already deep cross‐sections. Floodplain resistance along Cooper Creek is concentrated at channel bank tops where vegetation density is highest, and the subsequent flow retardation is transmitted across the surface of the channels over distances as large as 50–70 m. The rough floodplain surface affects flood wave transmission, producing significant decreases in wave speeds downstream. The character of the wave‐speed–discharge relationship also changes longitudinally, from log–linear in the upper reaches to nonlinear where the floodplain broadens appreciably. The nonlinear form is similar in several respects to relationships proposed for more humid rivers, with flood wave speed reaching an intermediate maximum at about four‐fifths bankfull discharge before decreasing to a minimum at approximately Q_2·33. It does not regain the value at the intermediate maximum until the 10 year flood, by which time floodplain depths have become relatively large and broad floodways more active. Copyright © 2002 John Wiley & Sons, Ltd.