Geophysical Finite-Element Simulation Tool (GeoFEST): Algorithms and Validation for Quasistatic Regional Faulted Crust Problems

GeoFEST (Geophysical Finite Element Simulation Tool) is a two- and three-dimensional finite element software package for the modeling of solid stress and strain in geophysical and other continuum domain applications. It is one of the featured high-performance applications of the NASA QuakeSim project. The program is targeted to be compiled and run on UNIX systems, and is running on diverse systems including sequential and message-passing parallel systems. Solution to the elliptical partial differential equations is obtained by finite element basis sampling, resulting in a sparse linear system primarily solved by conjugate gradient iteration to a tolerance level; on sequential systems a Crout factorization for the direct inversion of the linear system is also supported. The physics models supported include isotropic linear elasticity and both Newtonian and power-law viscoelasticity, via implicit quasi-static time stepping. In addition to triangular, quadrilateral, tetrahedral and hexahedral continuum elements, GeoFEST supports split-node faulting, body forces, and surface tractions. This software and related mesh refinement strategies have been validated on a variety of test cases with rigorous comparison to analytical solutions. These include a box-shaped domain with imposed motion on one surface, a pair of strike slip faults in stepover arrangement, and two community-agreed benchmark cases: a strike slip fault in an enclosing box, and a quarter-domain circular fault problem. Scientific applications of the code include the modeling of static and transient co- and post-seismic earth deformation, Earth response to glacial, atmospheric and hydrological loading, and other scenarios involving the bulk deformation of geologic media.     

Dynamic characteristics and horizontal transports of internal solitons generated at the Columbia River plume front

River plume front-generated internal solitons play an important role in the interaction between the plume and coastal waters. The internal solitons drive a non-harmonic velocity field, resulting in a horizontal transport that carries plume water seaward and redistributes nutrients and sediments. In this study, we present observations of internal solitons generated at the Columbia River plume front that separates the new, tidal plume, older plume and coastal waters. Scale analyses suggest that the plume front-generated internal solitons are highly non-linear waves, and their dynamic properties do not conform to any weakly non-linear theory. Thus, a high-order Korteweg–de Vries (KdV) theory is used to analyze the internal solitons. The comparison between theoretical values and cruise data shows that the high-order KdV model is much better than the weakly non-linear theories for prediction of the soliton dynamic parameters. Based on the model, we develop theoretical and numerical solutions of the soliton-induced upper layer horizontal transport and Lagrangian water parcel transport distance, which shows that the water particle drift, during the internal soliton passage, is as far as 1 km, and demonstrates the role of the internal solitons on the exchange between the plume and ambient coastal water. Energy fluxes caused by the internal solitons are estimated using the high-order KdV theory. The leading soliton fluxes 2.0×10³ W m⁻¹ per unit crest length, and carries energy of 4.2×10⁵ J m⁻¹. The total energy carried by the eight internal solitons is 1.6×10⁶ J m⁻¹, about 70% of the total frontal energy.     

Geochemistry,chronology and stratigraphy of Neogene tuffs of the Central Andean region

Tuff layers are vital stratigraphic tools that allow correlations to be made between widely dispersed exposures. Despite their widespread occurrence in the central Andes, tuffs from both natural exposures and sedimentary cores extracted from the region's extensive salars (salt pans) are relatively unstudied. Here we lay the foundation for a tephrostratigraphic framework in the central Andes (14–28°S) by chemically and morphologically characterizing ash shards, and in some cases dating 36 Neogene distal tuffs. These tuffs occur in lacustrine and alluvial deposits from the southern Bolivian Altiplano and adjacent Atacama Desert. All tuffs are calc-akaline rhyolites, consistent with their setting in the Central Andean Volcanic Zone. Five of the older tuffs were ⁴⁰Ar/³⁹Ar dated and yield an age range of 6.63–0.75 Ma. Organic material associated with tuffs deposited into paleolake sediments, paleowetland deposits, or urine-encrusted rodent middens provide constraints on the age of several Late Pleistocene and Holocene tuffs.These tuffs provide key stratigraphic markers and ages for lake cycles and archeological sites on the Bolivian Altiplano and for assessing rates of surficial processes and archeology in both the Atacama and Altiplano. While modern climate, and consequently questions about geomorphic processes and climate change, differs in the hyperarid Atacama and the semi-arid Altiplano, the most extensive air-fall tuffs covered both regions, placing the Atacama and the Bolivian Altiplano in the same tephrostratigraphic province. For example, the Escara B tuff (～1.85 Ma), can be securely identified in both the Altiplano and Atacama. On the Altiplano, dates from the Escara B and E tuffs securely establish the age of the Escara Formation—representing the oldest expansive lake documented on the Bolivian Altiplano. By contrast, the presence of the Escara B tuff below ～6 m of alluvial sediment at the Blanco Encalado site in the Atacama desert yields information about sedimentation rates in this hyperarid region. Indeed, most tuffs from the Atacama Desert are older than 600,000 years, even though they occur within fluvial terraces immediately adjacent to the alluvial fans that are still active. Most of these geomorphic surfaces in the Atacama also possess well-developed saline soils that, when combined with the radiometric ages of the distal tuffs, suggest slow rates of geomorphic change and exceptional landscape stability for this area during the Quaternary.In contrast, younger tuffs are more abundant in the more recent lake records of the Altiplano. The Chita tuff was deposited at ～15,650 cal yr B.P., during the regressive phase of the region's deepest late Quaternary lake cycle—the “Tauca lake cycle”—which spanned 18.1–14.1 cal yr B.P. Two Holocene tuffs, the Sajsi tuff and the Cruzani Cocha tuff, are widespread. The Sajsi tuff was deposited just before 1700 cal yr B.P., whereas the Cruzani Cocha tuff appears to be mid-Holocene in age and shows some chemical affinities to a Holocene tuff (202B) deposited between 4420 and 5460 cal yr B.P. in a urine-encased rodent midden in the Atacama Desert.     

Effects of soil layering on the characteristics of basin-edge induced surface waves

This paper presents the effects of soil layering on the characteristics of basin-edge induced surface waves and associated strain and aggravation factor. The simulated results revealed surface wave generation near the basin-edge. The first mode of induced Love wave was obtained in models having increasing velocity with depth and a large impedance contrast between the soil layers. Amplitude amplification or de-amplification of body waves was proportional to the impedance contrast between the soil layers. The average aggravation factor was inversely proportional to the impedance contrast between the soil layers in case of increasingvelocity models and proportional in case of decreasing-velocity basinedge models. On the other hand, the maximum strain was inversely proportional to the impedance contrast between the soil layers in both cases. On the average, strain was greater in case of increasing-velocity models but the average aggravation factor was greater in case of decreasingvelocity models.     

AMS 14C age constraints on geoglyphs in the lower colorado river region,arizona and california

Giant ground figures are widespread in the lower Colorado River area of southwestern North America, yet their chronology has remained unconstrained by numerical ages. Thirteen AMS ¹⁴C measurements reported here indicate that geoglyphs were made from before ˜A. D. 1200 to before ˜900 B. C. We account for potential contamination from prior organics in weathering rinds. All other potential errors point to ¹⁴C dates being minimum-limiting ages for the manufacturing of geoglyphs. Although these ages indicate considerable chronological complexity among geoglyphs, our data are consistent with the linguistic hypothesis that the Yuman people in the desert of southeastern California migrated from Baja California—rather than from the north. These results must, however, be placed under the cloud of uncertainty that hangs over the entire field of AMS dating of rock art: the untested assumption surrounding contemporeneity of organics in a surface context. © 1995 John Wiley & Sons, Inc.     

Simulation of urban expansion via integrating artificial neural network with Markov chain – cellular automata

Accurate simulations and predictions of urban expansion are critical to manage urbanization and explicitly address the spatiotemporal trends and distributions of urban expansion. Cellular Automata integrated Markov Chain (CA-MC) is one of the most frequently used models for this purpose. However, the urban suitability index (USI) map produced from the conventional CA-MC is either affected by human bias or cannot accurately reflect the possible nonlinear relations between driving factors and urban expansion. To overcome these limitations, a machine learning model (Artificial Neural Network, ANN) was integrated with CA-MC instead of the commonly used Analytical Hierarchy Process (AHP) and Logistic Regression (LR) CA-MC models. The ANN was optimized to create the USI map and then integrated with CA-MC to spatially allocate urban expansion cells. The validated results of kappa and fuzzy kappa simulation indicate that ANN-CA-MC outperformed other variously coupled CA-MC modelling approaches. Based on the ANN-CA-MC model, the urban area in South Auckland is predicted to expand to 1340.55 ha in 2026 at the expense of non-urban areas, mostly grassland and open-bare land. Most of the future expansion will take place within the planned new urban growth zone.     

Using gravity and topography-implied anomalies to assess data requirements for precise geoid computation

Many regions around the world require improved gravimetric data bases to support very accurate geoid modeling for the modernization of height systems using GPS. We present a simple yet effective method to assess gravity data requirements, particularly the necessary resolution, for a desired precision in geoid computation. The approach is based on simulating high-resolution gravimetry using a topography-correlated model that is adjusted to be consistent with an existing network of gravity data. Analysis of these adjusted, simulated data through Stokes’s integral indicates where existing gravity data must be supplemented by new surveys in order to achieve an acceptable level of omission error in the geoid undulation. The simulated model can equally be used to analyze commission error, as well as model error and data inconsistencies to a limited extent. The proposed method is applied to South Korea and shows clearly where existing gravity data are too scarce for precise geoid computation.     

Report of centimeter-sized tektites in Pima County,Arizona, cannot be verified

Abstract— Glassy objects reportedly found in Pima County, Arizona, have been identified as tektites. A field survey of the area where they were said to occur, however, did not produce any other tektites, nor did it reveal any other geologic features that might indicate a nearby impact crater. The major-, minor-, and trace-element composition of one specimen is similar to those measured in indochinites, which suggests the objects reportedly from Pima County were instead transported to southern Arizona from Indochina by people.     

Coronal images from the 1984 solar eclipse

We present digitized photographs of the white-light solar corona taken during the total solar eclipse of 22–23 November, 1984, on both calibrated black-and-white film and on color film. Conditions on site in Hula, Papua New Guinea, were exceptionally clear. The color image was used to produce an isophotal map of the inner corona, from which a flattening coefficient of 0.23 was measured. The black-and-white image was enhanced through a digital radial filter. Our images are the best processed images available from the 1984 eclipse and so provide important data for synoptic observations.