Thursday, June 25, 2009

Reconnaissance Geology of the State of Baja California

Reconnaissance Geology of the State of Baja California
Geological Society of America Memoir 140
Prepared by
R. Gordon Gastil, Richard P. Phillip, and Edwin C. Allison

This map was compiled in 1971 at a scale of 1:250,000 by the students and staff of Iniversidad Autonoma de Baja California and San Diego State University. It was published in 1975 as part of the Geological Society of America Memoir 140. Later in 2003 these maps were made available digitally in the GSA special paper Tectonic Evolution of Northwestern Mexico and the Southwestern United States. This volume was assembled as a dedication to the career achievements of R. Gordon Gastil. The papers focus primarily on the Mesozoic and Tertiary tectonic evolution of Baja California and nearby regions of mainland México and southern California. The volume begins with an extensive tectonic overview with the goals pf summarizing our present understanding of the geology of Baja California, addressing several controversial issues and providing a framework for more detailed discussions in some of the remaining papers. The volume contains a wealth of modern analytical data and authoritatively explores regional issues. The publications accompanying CD-ROM contained valuable geochronological and geochemical data, color images, and out-of-print Geological Society of America maps. In particular, the CD-ROM contains the highly sought after geologic maps of northern Baja California published by Gordon Gastil and colleagues in 1975 as part of the Geological Society of America Memoir 140, Reconnaissance Geology of the State of Baja California.

Monday, June 8, 2009

Geology Map of the San Diego 30' x 60' Quadrangle

Geology pf the San Diego 30' x 60' Quadrangle

Preliminary Geologic Map at 1:100,000 Scale to provide the public timely access to digital geology, prepared by the California Geological Survey Regional Geologic Mapping Project.

The San Diego 30' X 60' quadrangle was prepared by the Department of Conservation, California Geological Survey pursuant to a U.S. Geological Survey STATEMAP cooperative mapping award (# 1434-94-A-1224). It is a product of the Southern California Areal Mapping Project (SCAMP), a cooperative U.S. Geological Survey-California Geological Survey mapping project, This map is a compilation of published geological mapping (Fig. 1). The published mapping has been modified only to the extent necessary to integrate variables in nomenclature and scale. The onshore part of the map was digitized at a scale of 1:24,000 and the offshore part has been enlarged from 1:250,000. The quadrangle is between 32.5° and 33.0° N. latitude and 117.0° and 118.0° W. longitude. It encompasses the greater San Diego area, the second largest metropolitan area of California.
The area is tectonically active and is dissected by four major northwest-trending, oblique right slip faults that lie within the western part of the Pacific/North American Plate boundary. They include the Rose Canyon-Newport-Inglewood Fault Zone along the coastal margin, the Palos Verdes-Coronado Bank Fault Zone on the inner shelf, the San Diego Trough Fault Zone (origin of the 1986, ML=5.3, Oceanside earthquake) in the central offshore and the San Clemente Fault Zone on the outer offshore margin. Within the greater San Diego metropolitan area, the Rose Canyon Fault Zone as depicted by Kennedy and others (1975), Moore and Kennedy (1975), Kennedy and Welday (1980), Clarke and others (1987), Treiman (1993) and Kennedy and Clarke (2001) includes the Mount Soledad, Old Town, Point Loma, Silver Strand, Coronado and Spanish Bight faults. The Rose Canyon Fault Zone displaces Holocene sediment in Rose Canyon 7 km north of San Diego Bay where a late Pleistocene slip rate of 1-2 mm/yr has been estimated (Lindvall and Rockwell, 1995). A study of the recency and character of faulting in the greater San Diego metropolitan area suggests a long-term Tertiary slip rate for the Rose Canyon Fault Zone of about 1-2 mm/yr (Kennedy and others, 1975). Although there is significant late Quaternary deformation in the San Diego region the seismicity is relatively low (Simons, 1977).
The San Diego quadrangle is underlain by a thick sequence (>5 km) of Mesozoic fore-arc and fore-arc basin andesitic flows and coarse-grained volcaniclastic breccias that have been in large part metamorphosed to low-grade greenshist facies and are pervasively penetratively deformed. However, in the upper part of the section these rocks are not metamorphosed and are only moderately deformed. Marine sedimentary interbeds in Penasquitos Canyon, near Del Mar, contain the fossil Buchia piochii, which indicates a Late Jurassic (Tithonian) age for these strata (Fife and others, 1967; Jones and Miller, 1982). Zircon U/Pb ages from the metavolcanic rocks are reported to range from 137 Ma to 119 Ma (Anderson, 1991) indicating that they are coeval with the surrounding plutonic rocks of the western Peninsular Ranges batholith. The batholithic rocks are mostly granodiorite and tonalite and based on U-Pb isotopic ages range from 140 Ma to 105 Ma (Silver and Chappell, 1988). Much of the basement rock has been deeply weathered and altered. The weathered bedrock and Quaternary alluvial deposits derived from them contain expansible clays, mostly smectite.
The western part of the quadrangle is underlain by a relatively thick (>1,000 m) succession of Upper Cretaceous, Tertiary and Quaternary sedimentary rocks that unconformably overlie basement rocks. They consist of marine, paralic, and continental claystone, siltstone, sandstone and conglomerate. The Upper Cretaceous rocks are composed of marine turbidites and continental fan deposits assigned to the Rosario Group (Kennedy and Moore, 1971). The Lusardi Formation, the basal formation of the Rosario Group is a nonmarine boulder fanglomerate deposited along the western margin of a tectonic highland upon a deeply weathered surface of the older Cretaceous and Jurassic plutonic and metamorphic basement rocks. Clasts within the Lusardi Formation are composed exclusively of these weathered basement rocks. The Lusardi Formation is overlain by the Point Loma Formation, the middle part of the Rosario Group. It is composed mostly of marine sandstone, siltstone and conglomerate sequences that together form massive turbidite deposits. The Point Loma Formation is Campanian and Maestrichtian in age (Sliter, 1968; Bukry and Kennedy, 1969) and underlies most of the Point Loma Peninsula and the hills southeast of La Jolla. It is conformably overlain by the uppermost part of the Rosario Group, marine sandstone and conglomerate of the Maestrichtian (Sliter, 1968; Bukry and Kennedy, 1969) Cabrillo Formation. Following the deposition of the Rosario Group, the San Diego coastal margin underwent uplift and erosion until the middle Eocene when nine partially intertonguing middle and upper Eocene sequences composed of siltstone, sandstone, and conglomerate were deposited during several major transgressive-regressive cycles. The succession is over 700 meters thick and grades from nonmarine fan and dune deposits on the east through lagoonal and nearshore beach and beach-bar deposits to marine continental shelf deposits on the west near the present-day coastline. The age and environmental interpretation of the Eocene sequence is based on the mapped distribution of lithofacies coupled with the presence of a pelagic fossil calcareous nannoplankton flora in the continental shelf facies (e.g., Bukry and Kennedy, 1969), a shallow water molluscan fauna in the nearshore facies (e.g., Givens and Kennedy, 1979), and a fossil terrestrial vertebrate mammal fauna in the paralic facies (e.g., Golz, 1973). Cross bedding, cobble imbrications, paleo-stream gradients and clast petrology indicate a local eastern source for these rocks. The nonmarine facies of the Eocene formations are typically well indurated and cemented whereas the lagoonal facies are soft and friable. The nearshore facies are well indurated, well sorted, and locally concretionary. The marine deposits are typically fine-grained, indurated, and cemented. Following the deposition of Eocene rocks the San Diego margin was again elevated and eroded. During the Oligocene, continental and shallow water lagoonal deposits of the Otay Formation, were deposited. The Otay Formation is light-gray and light- brown, medium- and coarse-grained, arkosic sandstone intertongued with light-brown siltstone and light-gray claystone. Much of the claystone is composed of lightgray bentonite in beds up to 1 m in thickness. Following Oligocene time the San Diego coastal margin underwent uplift and extensive erosion. The next major marine transgression did not occur until Pliocene time when the strata of the San Diego Formation were deposited. The San Diego Formation rests unconformably upon Oligocene, Eocene and Upper Cretaceous beds across its outcrop from Pacific Beach to the International border with Mexico. The San Diego Formation is late Pliocene in age and contains a rich molluscan fauna (e.g., Arnold, 1903; Demere, 1983). It consists mostly of yellowish-brown and gray, fine- to mediumgrained, marine sandstone and reddish-brown, transitional marine and nonmarine pebble and cobble conglomerate. Following the deposition of the San Diego Formation and continuing to the present time, the San Diego coastal margin has undergone relatively steady uplift (Fig. 2). A series of continually evolving marine abrasion platforms have been carved and uplifted during this time and are manifest in the marine terraces and their deposits that are ubiquitous to the San Diego coastal region. The deposits consist of nearshore marine, beach, estuarine, lagoonal and continental dune facies that were deposited across a marine/nonmarine transition zone and along a coastal strandline. Changes in sea level coupled with regional uplift give rise to the preservation and/or obliteration of both the abrasion platforms and their overlying deposits (e.g. Lajoie, and others, 1991; Kern and Rockwell, 1992; Kern, 1996a, 1996b).
The authors appreciate very helpful reviews by Victoria R.Todd and J. Philip Kern

Friday, May 1, 2009

Geology of San Diego

A collection of geologic maps from the San Diego Metropolitan Area. Maps were created using Super Image Overlays.

Saturday, April 18, 2009

Geology Map Symbol KML Generator

Geologic maps use a combination of colors, lines, and symbols to give a a third dimension to the rock layers on a flat map surface.

With the Geology Map Symbol Generator you can add geologic map symbols to google earth using a web-based symbol generator.

This is still a Beta version which simply posts html form data to php script.

Sunday, March 29, 2009

USGS Digital Data Series DDS-9

National Geophysical Data Grids:
Gamma-Ray, Gravity, Magnetic, and Topographic
Data for the Conterminous United States


The data published in the \ASCII directory on this CD-ROM consist of regular grids of ASCII values.Each grid has a different origin and sample spacing as defined in its first record, the header.The grid format is explained in detail below.Each grid represents data that have been projected from latitude and longitude coordinates into map coordinates of kilometers.Thus, the x-origin and y-origin values given in the header represent the distance in kilometers from the central meridian and base latitude of the geographic projection (discussed below) to the lower left corner of the grid.
Binary versions of these grid files, suitable for processing on IBM or compatible personal computers using the potential-field software contained on this CD-ROM, are located in the \DOSBIN directory on this CD-ROM.Use the EXTRACT BINARY GRID FILES submenu to access these grid files.For more information on the binary grid file format and the potential-field software, access the POTENTIAL-FIELD SOFTWARE submenu from the MAIN menu.
Projection information

All grids on this CD-ROM have been projected by using an Albers equal-area conic projection with standard parallels of 29.5 degrees and 45.5 degrees north.A central meridian of 96 degrees west and a base latitude of 0 degrees were used.The projection is referenced to the Clarke 1866 ellipsoid, which has an equatorial radius of 6378.2064 km and a polar radius of 6356.5838 km.
Grid format

Each ASCII grid consists of two header records followed by a series of data records.The first header record contains 56 alphanumeric characters of data identification, 8 alphanumeric characters containing the name of the computer program that created the grid, the integer number of columns in the grid, and the integer number of rows in the grid.The second header record contains a dummy integer value of 1, the x-coordinate in kilometers of the first (leftmost) column, the distance in kilometers between columns of the grid, the y-coordinate in kilometers of the first (bottommost) row of the grid, and the distance in kilometers between rows of the grid.

Each data record contains five values in scientific notation.The bottommost row is presented first, starting from the leftmost column.The first value in each row is a dummy value, usually zero (0.000000000E+00). If a row ends in the middle of a record, the record is padded with zeros. Areas of the grid containing no data are represented by a special value, 0.999999968E+38.

The gamma-ray data grids

Aerial gamma-ray surveys measure the gamma-ray flux produced by the radioactive decay of the naturally occurring elements K-40, U-238, and Th-232 in the top few centimeters of rock or soil (Duval, Cook, and Adams, 1971).If the gamma-ray system is properly calibrated (for example, see Grasty and Darnley, 1971), the data can be expressed in terms of the estimated concentrations of the radioactive elements.The potassium concentration data are usually expressed in units of percent potassium (percent K), uranium as parts per million equivalent uranium (ppm eU), and the thorium as parts per million equivalent thorium(ppm eTh).The term equivalent is used because the technique actually measures the gamma-ray flux from the decay of bismuth (Bi-214), which is a decay product of U-238, and from the decay of thallium (Tl-208), which is a decay product of Th-232.Radioactive disequilibrium in the thorium decay series may cause the measured equivalent uranium and equivalent thorium to differ from the actual uranium andthorium present in the surface rocks and soils. Because Rn-222 is a daughter product of the U-238 decay series, the U-238 concentrations can also be used to estimate the amounts of Rn-222 in the near-surface soil gas.

During the period 1975-83, the U.S. Department of Energy carried out the National Uranium Resource Evaluation (NURE) Program, which included aerial gamma-ray surveys of most of the conterminous United States.Although many of the airborne gamma-ray systems used to make these surveys were calibrated, many of the earlier surveys were done without calibration and conversion to the concentrations of the radioactive elements.Detailed examinations of the digital data available on magnetic tape also showed that many of the "calibrated" surveys do not match the data from other "calibrated" surveys of adjacent areas.For these reasons, the data must be corrected to obtain a consistent data base for the conterminous United States.
Because uranium, thorium, and potassium concentration data are useful in geologic studies and because the NURE data are the only nationwide data base on the natural radiation environment, the U.S. Geological Survey (USGS) reprocessed the aerial gamma-ray data to produce maps showing surface concentrations of potassium, uranium, and thorium for the conterminous United States.These maps have been released as USGS Open- File Reports (Duval and others, 1989, 1990).Some of the reprocessed data have also been released in profile form (Duval, 1995).

The magnetic anomaly data grid

Magnetic anomalies are produced by variations in the distribution of iron minerals, usually magnetite, in the rocks of the Earth's crust.Igneous and metamorphic rocks can be very magnetic.By comparison, sedimentary rocks are usually nonmagnetic.Magnetic anomalies therefore provide a way of mapping exposed and buried crystalline rocks.
The grid of magnetic anomaly data for the conterminous United States and adjacent marine areas (Godson, 1986) was created from digitized contours of the east half of the Composite Magnetic Anomaly Map of the United States, Part A (U.S. Geological Survey, 1982), and the Composite Magnetic Anomaly Map of the Conterminous United States West of 96 Degrees Longitude (Bond and Zietz, 1987),with additional data used in the compilation of the Magnetic Anomaly Map of North America (Geological Society of America, Committee for the Magnetic Anomaly Map of North America, 1987).A regional gradient present in the 1982 map was removed by using a corrected geomagnetic reference field (Godson, 1986).The data, originally gridded on a 2-km interval using the spherical Transverse Mercator projection of the Magnetic Anomaly Map of North America, were reprojected to the Albers projection used on this CD-ROM and regridded on a 2-km interval using a minimum curvature gridding program (Webring, 1981).An interpretation of the 1982 anomaly map was presented by Hinze and Zietz (1985).

The Bouguer gravity anomaly data grid

Gravity anomalies are produced by density variations within the rocks of the Earth's crust and upper mantle.Mapping of these density variations is the primary use of gravity anomalies.
Gravity measurements made on the surface of the Earth must be corrected in various ways before they can be made into an anomaly map.The free-air correction reduces the measurement to sea level by assuming that there is no intervening mass.The simple Bouguer correction accounts for the intervening mass as a uniform slab of constant density, and the complete Bouguer correction includes the effects of constant density topography within 166.7 km of the measurement location.A gravity reference field is subtracted from the corrected measurements to produce the free-air, simple Bouguer, or complete Bouguer anomaly.
The grid of gravity data for the conterminous United States and adjacent marine areas (Godson and Scheibe, 1982; Godson, 1985) was constructed from Defense Mapping Agency gravity data files.The onshore data consisted of nearly one million Bouguer gravity anomaly values computed by using a reduction density of 2.67 grams per cubic centimeter.The offshore data consisted of approximately 800,000 free-air gravity anomaly values. Because the Bouguer anomaly equals the free-air anomaly at sea level, there is no discontinuity in the gridded data at the shoreline.All computations were performed by using the International Gravity Standardization Net of 1971 (International Association of Geodesy, 1974) and the 1967 Geodetic Reference System formula for theoretical gravity (International Association of Geodesy, 1971).
In areas of substantial relief, terrain corrections were computed about each station location at radial distances of 0.895 km to 166.7 km by using a density of 2.67 grams per cubic centimeter.The data were projected and gridded on a 4-km interval using minimum curvature (Webring, 1981).The gridding procedure resulted in the extrapolation of grid values up to 40 km beyond the limits of the data; therefore, values around the edges should be viewed with caution.These gridded data were published in map form as the Gravity Anomaly Map of the United States (Society of Exploration Geophysicists, 1982).This map was further discussed by O'Hara and Lyons (1985) and Kane and Godson (1985, 1989).

The isostatic residual gravity anomaly data grid

Isostatic residual gravity anomaly maps are produced by subtracting long- wavelength anomalies produced by masses deep within the crust or mantle from the Bouguer anomalymap.The long-wavelength anomalies are assumed to result from isostatic compensation of topographic loads.Isostatic residual gravity anomaly maps therefore reveal more clearly than Bouguer anomaly maps the density distributions within the upper crust that are of interest in many geologic and tectonic studies.
The grid of isostatic residual gravity anomaly data (Simpson and others, 1986) was produced from the grid of Bouguer gravity anomaly data (Godson and Scheibe, 1982) by using an Airy-Heiskanen compensation model (Heiskanen and Moritz, 1967) with three parameters.The depth to the compensating root at sea level was chosen to be 30 km.The density contrast across the root was chosen to be 0.35 grams per cubic centimeter, and the density of the topography was chosen to be 2.67 grams per cubic centimeter.Other reasonable choices of these parameters would produce similar-looking residual maps.
The computer program and topographic data sets used to produce the data
grid were described by Simpson and others (1983a,b). The data were
published in map form by Jachens and others (1985).Interpretations of the
isostatic residual gravity anomaly map were presented by Simpson and others
(1986) and by Jachens and others (1989).

The topographic data grid

The topographic data grid for the conterminous United States and adjacent areas was constructed from 30x30 second digital terrain files used by the U.S. Geological Survey for the reduction of gravity data.Elevations are in meters; sea level elevations are listed as 1 meter.

The topographic-bathymetric data grid

The topographic and bathymetric data grid for the conterminous United States and adjacent areas was constructed from 5x5 minute North American topographic data and 5x5 minute Synthetic Bathymetric Profiling System data available from the National Oceanic and Atmospheric Administration.The data were interpolated onto an 8x8 km grid in the Albers projection (Simpson and others, 1983b).Elevation units are in meters relative to sea level.