Engineering Papers⌕ Search

Engineering topics

Donnellan, Andrea

Publications and source records attributed to Donnellan, Andrea.

67 records · Page 4
All sourcesNASA NTRSDOE OSTI
NASA NTRS2004

Using the GeoFEST Faulted Region Simulation System

GeoFEST (the Geophysical Finite Element Simulation Tool) simulates stress evolution, fault slip and plastic/elastic processes in realistic materials, and so is suitable for earthquake cycle studies in regions such as Southern California. Many new capabilities and means of access for GeoFEST are now supported. New abilities include MPI-based cluster parallel computing using automatic PYRAMID/Parmetis-based mesh partitioning, automatic mesh generation for layered media with rectangular faults, and results visualization that is integrated with remote sensing data. The parallel GeoFEST application has been successfully run on over a half-dozen computers, including Intel Xeon clusters, Itanium II and Altix machines, and the Apple G5 cluster. It is not separately optimized for different machines, but relies on good domain partitioning for load-balance and low communication, and careful writing of the parallel diagonally preconditioned conjugate gradient solver to keep communication overhead low. Demonstrated thousand-step solutions for over a million finite elements on 64 processors require under three hours, and scaling tests show high efficiency when using more than (order of) 4000 elements per processor. The source code and documentation for GeoFEST is available at no cost from Open Channel Foundation. In addition GeoFEST may be used through a browser-based portal environment available to approved users. That environment includes semi-automated geometry creation and mesh generation tools, GeoFEST, and RIVA-based visualization tools that include the ability to generate a flyover animation showing deformations and topography. Work is in progress to support simulation of a region with several faults using 16 million elements, using a strain energy metric to adapt the mesh to faithfully represent the solution in a region of widely varying strain.

Geophyical Finite Element Simulation Tool (GeoFEST↗
NASA NTRS2004

QuakeSim and the Solid Earth Research Virtual Observatory

We are developing simulation and analysis tools in order to develop a solid Earth science framework for understanding and studying active tectonic and earthquake processes. The goal of QuakeSim and its extension, the Solid Earth Research Virtual Observatory (SERVO), is to study the physics of earthquakes using state-of-the-art modeling, data manipulation, and pattern recognition technologies. We are developing clearly defined accessible data formats and code protocols as inputs to simulations, which are adapted to high-performance computers. The solid Earth system is extremely complex and nonlinear resulting in computationally intensive problems with millions of unknowns. With these tools it will be possible to construct the more complex models and simulations necessary to develop hazard assessment systems critical for reducing future losses from major earthquakes. We are using Web (Grid) service technology to demonstrate the assimilation of multiple distributed data sources (a typical data grid problem) into a major parallel high-performance computing earthquake forecasting code. Such a linkage of Geoinformatics with Geocomplexity demonstrates the value of the Solid Earth Research Virtual Observatory (SERVO) Grid concept, and advances Grid technology by building the first real-time large-scale data assimilation grid.

virtual observatory↗
NASA NTRS2004

The QuakeSim Project: Numerical Simulations for Active Tectonic Processes

In order to develop a solid earth science framework for understanding and studying of active tectonic and earthquake processes, this task develops simulation and analysis tools to study the physics of earthquakes using state-of-the art modeling, data manipulation, and pattern recognition technologies. We develop clearly defined accessible data formats and code protocols as inputs to the simulations. these are adapted to high-performance computers because the solid earth system is extremely complex and nonlinear resulting in computationally intensive problems with millions of unknowns. With these tools it will be possible to construct the more complex models and simulations necessary to develop hazard assessment systems critical for reducing future losses from major earthquakes.

tectonic processes↗
NASA NTRS2004

Complexity Computational Environment: Data Assimilation SERVOGrid

We are using Web (Grid) service technology to demonstrate the assimilation of multiple distributed data sources (a typical data grid problem) into a major parallel high-performance computing earthquake forecasting code. Such a linkage of Geoinformatics with Geocomplexity demonstrates the value of the Solid Earth Research Virtual Observatory (SERVO) Grid concept, and advance Grid technology by building the first real-time large-scale data assimilation grid Here we develop the next steps for both the SERVO concept and the identified need for a Solid Earth problem-solving environment. We use a challenging motivating problem of importance to NASA namely integrating NASA space geodetic observations with numerical simulations of a changing earth.

data assimiliation↗
NASA NTRS1999

Deployment of Autonomous GPS Stations in Marie Byrd Land, Antarctica

During the 1998-1999 Antarctic field season, we installed three autonomous GPS stations in Marie Byrd Land, West Antarctica to measure glacio-isostatic rebound and rates of spreading across the West Antarctic Rift System. The systems collect data throughout the entire year and therefore must function during the warm, relatively mild summer, and cold, harsh winters. They are powered by gel cell batteries that are charged by wind and solar power. The system includes dual data logging capability. We log data at 5 minute intervals within the receiver and at 30 second intervals to a serial data logger. We do not require 365 days of continuous data for well determined crustal velocities, but rather long periods (greater than 24 hours) of continuous data distributed throughout the year. Therefore, for simplicity, we designed the system to accept occasional data interruptions. The batteries, in addition to supplying power, act as a thermal capacitive heat storage device to help regulate the temperatures within the system. This storage system absorbs the majority of the 10-15 watts of power from the receiver and 5 watts from the data logger which helps to maintain temperature for long periods of time. Power is switched off when the temperature within the system enclosure reaches 50'C and is reconnected at 20 C. If battery voltage drops too low the batteries will freeze. Therefore, we cut the power off when the batteries drop to a low voltage of 12.45V. Power is restored at 13.2V. The temperature and power hysteresis allows for a minimum of several days of data to be collected before system shutdowns. A check of all three stations in late January 1999 indicated that the thermal and power control systems are performing as expected. We plan to implement satellite telemetry to the systems during the 2000-2001 season following a year of development.

Donnellan, Andrea↗
NASA NTRS1997

Rate Change Observed at JPLM after the Northridge Earthquake

Geodetic time series detemined with the Global Positioning System indicate that the geodetic rate of a permanent site in Pasadena, Califonia (JPLM) changed significantly after the 17 January 1994 Northridge California earthquake.

Northridge↗
NASA NTRS1994

Hydro-isostatic deflection and tectonic tilting in the central Andes: Initial results of a GPS survey of Lake Minchin shorelines

Sufficiently large lake loads provide a means of probing rheological stratification of the crust and upper mantle. Lake Minchin was the largest of the late Pleistocene pluvial lakes in the central Andes. Prominent shorelines, which formed during temporary still-stands in the climatically driven lake level history, preserve records of lateral variations in subsequent net vertical motions. At its maximum extent the lake was 140 m deep and spanned 400 km N-S and 200 km E-R. The load of surficial water contained in Lake Minchin was sufficient to depress the crust and underlying mantle by 20-40 m, depending on the subjacent rheology. Any other differential vertical motions will also be recorded as departures from horizontality of the shorelines. We recently conducted a survey of shoreline elevations of Lake Minchin with the express intent of monitoring the hydro-isostatic deflection and tectonic tilting. Using real-time differential Global Positioning System (GPS), we measured topographic profiles across suites of shorelines at 15 widely separated locations throughout the basin. Horizontal and vertical accuracies attained are roughly 30 and 70 cm, respectively. Geomorphic evidence suggests that the highest shoreline was occupied only briefly (probably less than 200 years) and radiocarbon dates on gastropod shells found in association with the shore deposits constrain the age to roughly 17 kyr. The basin-side pattern of elevations of the highest shoreline is composed of two distinct signals: (27 +/- 1) m of hydro-isostatic deflection due to the lake load, and a planar tilt with east and north components of (6.8 +/- 0.4) 10(exp -5) and 9-5.3 +/- 0.3) 10(exp -5). This rate of tilting is too high to be plausibly attributed to steady tectonism, and presumably reflects some unresolved combination of tectonism plus the effects of oceanic and lacustrine loads on a laterally heterogeneous substrate. The history of lake level fluctuations is still inadequately known to allow detailed inferences of crust and mantle rheology. However, it is already clear that the effective elastic plate thickness is closer to 40 km than the 60-70 km crustal thickness in the central Andes and the effective viscosity is less than 5 10(exp 20) Pa s.

Bills, Bruce G.↗
NASA NTRS1993

Space geodetic measurement of crustal deformation in central and southern California, 1984-1992

We estimate the velocity field in central and southern California using Global Positioning System (GPS) observations from 1986 to 1992 and very long baseline interferometry (VLBI) observations from 1984 to 1991. Our core network includes 12 GPS sites spaced approximately 50 km apart, mostly in the western Transverse Ranges and the coastal Borderlands. The precision and accuracy of the relative horizontal velocities estimated for these core stations are adequately described by a 95% confidence ellipse with a semiminor axis of approximately 2 mm/yr oriented roughly north-south, and a semimajor axis of approximately 3 mm/yr oriented east-west.

Feigl, Kurt L.↗
NASA NTRS1993

Geodetic measurement of deformation in the Ventura basin region, southern California

We have measured the deformation in the Ventura basin region, southern California, with Global Positioning System (GPS) measurements carried out over 4.6 years between 1987 and 1992. The deformation within our network is spatically variable on scales of tens of kilometers, with strain rates reaching 0.6 +/- 1 micro-rad/yr in the east-central basin. Blocklike rotations are observed south and northwest of the basin where the maximum shear strain rates are an order of magnitude lower (0.06 +/- 1 micro-rad/yr to the south). We also observed clockwise rotations of 1 deg - 7 deg/m.y. Shear strain rates determined by comparing angle changes from historical triangulation spanning several decades and GPS measurements give consistent, though less precise, results. The geodetic rates of shortening across the basin and Western Transverse Ranges are lower than those estimated from geological observations, but the patterns of deformation from the two methods agree qualitatively.

Donnellan, Andrea↗
NASA NTRS1988

Determination of convergence rates across the Ventura Basin, Southern California, using GPS and historical triangulation

Comparison of angles from historical triangulation observations dating as far back as 1932 with Global Positions System (GPS) measurements taken in 1987 indicates that rapid convergence may be taking place on decade timescales in the central and eastern part of the Ventura basin, an east-west trending trough bounded by thrust faults. Changes in angles over this time were analyzed using Prescott's modified Frank's method and in terms of a model which assumes that the regions to the north and south of the basin are rigid blocks undergoing relative motion. For the two block model, inversion of the observed angle changes over the last 28 years for the relative motion vector leads to north-south convergence across the basin of 30 + or - 5 mm/yr, with a left lateral component of 10 + or - 1 mm/yr in the Fillmore-Santa Paula area in the central part of the basin. The modified Frank's method yields strain rates of approximately 2 microrad/yr in both the east and central parts of the basin for measurements spanning the 1971 San Fernando earthquake. Assuming no east-west strain yeilds north-south compression of approximately 3.5 + or - .2 cm/yr. Comparison of triangulation data prior to the earthquake shows no strain outside the margin of error. The convergence rates determined by geodetic techniques are consistent with geologic observations in the area. Such large geodetic deformation rates, with no apparent near-surface creep on the major thrust, can be understood if these faults become subhorizontal at relatively shallow depths and if the subhorizontal portions of the faults are creeping. An alternative explanation of the large displacement rates might be that the pumping of oil in the vicinity of the benchmarks caused large horizontal motions, although it is unlikely that meter scale horizontal motions are due to oil withdrawal. These and other hypotheses are evaluated to better constrain the tectonics of this active region.

Donnellan, Andrea↗