Cassini imaging of Saturn's rings: II. A wavelet technique for analysis of density waves and other radial structure in the rings

We describe a powerful signal processing method, the continuous wavelet transform, and use it to analyze radial structure in Cassini ISS images of Saturn's rings. Wavelet analysis locally separates signal components in frequency space, causing many structures to become evident that are difficult to observe with the naked eye. Density waves, generated at resonances with saturnian satellites orbiting outside (or within) the rings, are particularly amenable to such analysis. We identify a number of previously unobserved weak waves, and demonstrate the wavelet transform's ability to isolate multiple waves superimposed on top of one another. We also present two wave-like structures that we are unable to conclusively identify. In a multi-step semi-automated process, we recover four parameters from clearly observed weak spiral density waves: the local ring surface density, the local ring viscosity, the precise resonance location (useful for pointing images, and potentially for refining saturnian astrometry), and the wave amplitude (potentially providing new constraints upon the masses of the perturbing moons). Our derived surface densities have less scatter than previous measurements that were derived from stronger non-linear waves, and suggest a gentle linear increase in surface density from the inner to the mid-A Ring. We show that ring viscosity consistently increases from the Cassini Division outward to the Encke Gap. Meaningful upper limits on ring thickness can be placed on the Cassini Division (3.0 m at r∼118,800 km, 4.5 m at r∼120,700 km) and the inner A Ring (10-15 m for r<127,000 km).     

Does Magnetic Flux Submerge at Flux Cancelation Sites?

Simultaneous measurements of the magnetic fields in the photosphere and chromosphere were used to investigate if magnetic flux is submerging at sites between adjacent opposite polarity magnetic network elements in which the flux is observed to decrease or `cancel'. These data were compared with chromospheric and coronal intensity images to establish the timing of the emission structures associated with these magnetic structures as a function of height. We found that most of the cancelation sites show either that the bipole is observed longer in the photosphere than in the chromosphere and corona (44%) or that the timing difference of the disappearance of the bipole between these levels of the atmosphere is unresolved. The magnetic axis lengths of the structures associated with the cancelation sites are on average slightly smaller in the chromosphere than the photosphere. These observations suggest that magnetic flux is retracting below the surface for most, if not all, of the cancelation sites studied.     

Protracted thrusting followed by late rapid cooling of the Greater Himalayan Sequence,Annapurna Himalaya,Central Nepal: Insights from titanite petrochronology

The Greater Himalayan Sequence (GHS) has commonly been treated as a large coherently deforming high‐grade tectonic package, exhumed primarily by simultaneous thrust‐ and normal‐sense shearing on its bounding structures and erosion along its frontal exposure. A new paradigm, developed over the past decade, suggests that the GHS is not a single high‐grade lithotectonic unit, but consists of in‐sequence thrust sheets. In this study, we examine this concept in central Nepal by integrating temperature–time (T–t) paths, based on coupled Zr‐in‐titanite thermometry and U–Pb geochronology for upper GHS calcsilicates, with traditional thermobarometry, textural relationships and field mapping. Peak Zr‐in‐titanite temperatures are 760–850°C at 10–13 kbar, and U–Pb ages of titanite range from c. 30 to c. 15 Ma. Sector zoning of Zr and distribution of U–Pb ages within titanite suggest that diffusion rates of Zr and Pb are slower than experimentally determined rates, and these systems remain unaffected into the lower granulite facies. Two types of T–t paths occur across the Chame Shear Zone (CSZ). Between c. 25 and 17–16 Ma, hangingwall rocks cool at rates of 1–10°C/Ma, while footwall rocks heat at rates of 1–10°C/Ma. Over the same interval, temperatures increase structurally upwards through the hangingwall, but by 17–16 Ma temperatures converge. In contrast, temperatures decrease upwards in footwall rocks at all times. While the footwall is interpreted as an intact, structurally upright section, the thermometric inversion within the hangingwall suggests thrusting of hotter rocks over colder from c. 25 to c. 17–16 Ma. Retrograde hydration that is restricted to the hangingwall, and a lithological repetition of orthogneiss are consistent with thrust‐sense shear on the CSZ. The CSZ is structurally higher than previously identified intra‐GHS thrusts in central Nepal, and thrusting duration was 3–6 Ma longer than proposed for other intra‐GHS thrusts in this region. Cooling rates for both the hangingwall and footwall of the CSZ are comparable to or faster than rates for other intra‐GHS thrust sheets in Nepal. The overlap in high‐T titanite U–Pb ages and previously published muscovite ⁴⁰Ar/³⁹Ar cooling ages imply cooling rates for the hangingwall of ≥200°C/Ma after thrusting. Causes of rapid cooling include passive exhumation driven by a combination of duplexing in the Lesser Himalayan Sequence, and juxtaposition of cooler rocks on top of the GHS by the STDS. Normal‐sense displacement does not appear to affect T–t paths for rocks immediately below the STDS prior to 17–16 Ma.     

Multi-scale simulation of wave propagation and liquefaction in a one-dimensional soil column: hybrid DEM and finite-difference procedure

The paper describes a multi-phase, multi-scale rational method for modeling and predicting free-field wave propagation and the weakening and liquefaction of near-surface soils. The one-dimensional time-domain model of a soil column uses the discrete element method (DEM) to track stress and strain within a series of representative volume elements (RVEs), driven by seismic rock displacements at the column base. The RVE interactions are accomplished with a time-stepping finite-difference algorithm. The method applies Darcy’s principle to resolve the momentum transfer between a soil’s solid matrix and its interstitial pore fluid. Different algorithms are described for the dynamic period of seismic shaking and for the post-shaking consolidation period. The method can analyze numerous conditions and phenomena, including site-specific amplification, down-slope movement of sloping ground, dissolution or cavitation of air in the pore fluid, and drainage that is concurrent with shaking. Several refinements of the DEM are described for realistically simulating soil behavior and for solving a range of propagation and liquefaction factors, including the poromechanic stiffness of the pore fluid and the pressure-dependent drained stiffness of the grain matrix. The model is applied to four sets of well-documented centrifuge studies. The verification results are favorable and highlight the importance of the pore fluid conditions, such as the amount of dissolved air within the pore water.

     

Role of aeolian sediment accretion in the formation of heuweltjie earth mounds,western South Africa

The action of organisms in shaping landforms is increasingly recognized; the field of biogeomorphology and the conceptual framework of ecosystem engineering have arisen in response to the need for integrated studies of the interactions between biotic and abiotic components of landscapes. Pathways by which organisms influence landscape development may be complex. For example, primary change initiated by one biotic element may initiate a cascade of other changes that eventually produce a significant landscape modification. Mound‐like landforms in North America and southern Africa are widely cited examples of biogenic structures, yet there is considerable controversy regarding the processes responsible for their formation. Heuweltjies (Afrikaans for little hills) are circular mounds ranging from 10–30 m diameter and 0.5–2 m height and are widespread in western South Africa. Colonies of the termite (Microhodotermes viator) are typically associated with heuweltjies and some investigators have attributed heuweltjie formation to the direct action of termites in redistributing earthen materials. However, rather than being directly responsible in this way, termites simply create nutrient‐rich islands, which support denser vegetation, thereby inducing the localized accretion of aeolian sediments and upward growth of mounds. Contrasting soil features in heuweltjies in one locale indicate these processes have occurred throughout the late Quaternary. Geographic variation in sizes of mounds is explained in part by the local availability of sediments that can be mobilized and redistributed by the wind. Recognition of the operation of aeolian processes in the formation of heuweltjies has important implications for conservation. Any land use that diminishes the sediment‐trapping effect of vegetation on heuweltjies truncates the very process by which new aeolian materials can accrue and may promote irreversible erosion and landscape degradation. Copyright © 2014 John Wiley & Sons, Ltd.     

Data reporting norms for 40Ar/39Ar geochronology

Data reported in ⁴⁰Ar/³⁹Ar geochronology studies are commonly insufficient to allow computation of ages. This deficiency renders it difficult to compare ages based on different standards or constants, and often hinders critical evaluation of the results. Herein are presented an enumeration of the data that should be reported in all ⁴⁰Ar/³⁹Ar studies, including a discussion in support of these requirements. The minimum required data are identified and distinguished from parameters that are useful but may be derived from them by calculation. Finally, recommendations are made for metadata needed to document age calculations (e.g., from age spectrum or isochron analyses).     

SSBUV and NOAA-11 SBUV/2 Solar Variability Measurements

The Shuttle SBUV (SSBUV) and NOAA-11 SBUV/2 instruments measured solar spectral UV irradiance during the maximum and declining phase of solar cycle 22. The SSBUV data accurately represent the absolute solar UV irradiance between 200–405 nm, and also show the long-term variations during eight flights between October 1989 and January 1996. These data have been used to correct long-term sensitivity changes in the NOAA-11 SBUV/2 data, which provide a near-daily record of solar UV variations over the 170–400 nm region between December 1988 and October 1994. The NOAA-11 data demonstrate the evolution of short-term solar UV activity during solar cycle 22.     

A Comparative Plume Study of DRO,GRO, Benzene,and MTBE: Implications for Risk Management

Groundwater remediation and no-further action decision making at petroleum underground storage tank (UST) sites has largely been based on an understanding of plume length, plume stability, and attenuation rates for key hydrocarbon constituents. Regulatory guidance to support and guide such decisions is based in part on plume studies involving individual hydrocarbon constituents, namely benzene and methyl tert-butyl ether (MTBE). Questions remain regarding whether current guidance is applicable to chemical mixtures such as gasoline range organics (GRO), diesel range organics (DRO), and oxygen containing organic compounds (OCOCs) resulting from hydrocarbon biodegradation. To help address this concern, data from California's GeoTracker database were used to estimate maximum plume lengths, plume stability, and attenuation rates of DRO (which can be used as an analytical surrogate for OCOCs) and GRO relative to benzene and MTBE. The distributions of maximum plume lengths were similar for the four constituents with medians ranging from 27 to 32 m. The fraction of monitoring wells with a decreasing concentration trend ranged from 19% for DRO to 40% for MTBE, while fewer than 7% of the wells had an increasing concentration trend for any of the constituents. Median attenuation rates ranged from 0.10% day⁻¹ for DRO to 0.17% day⁻¹ for MTBE. The results suggest attenuation based risk management is appropriate for DRO and GRO plumes at most petroleum UST sites.     

A spatially and temporally adaptive solution of Richards’ equation

Efficient, robust simulation of groundwater flow in the unsaturated zone remains computationally expensive, especially for problems characterized by sharp fronts in both space and time. Standard approaches that employ uniform spatial and temporal discretizations for the numerical solution of these problems lead to inefficient and expensive simulations. In this work, we solve Richards’ equation using adaptive methods in both space and time. Spatial adaption is based upon a coarse grid solve and a gradient error indicator using a fixed-order approximation. Temporal adaption is accomplished using variable order, variable step size approximations based upon the backward difference formulas up to fifth order. Since the advantages of similar adaptive methods in time are now established, we evaluate our method by comparison with a uniform spatial discretization that is adaptive in time for four different one-dimensional test problems. The numerical results demonstrate that the proposed method provides a robust and efficient alternative to standard approaches for simulating variably saturated flow in one spatial dimension.     

Rockfall susceptibility and runout in the Valley of the Kings

Natural Hazards - The UNESCO world heritage site Valley of the Kings or Wadi el-Moluk (???? ??????) near Luxor, Egypt, hosts unique...