SEARCH · Search USGS

Results for “Mountain Views”

Search indexed USGS publications on groundwater, aquifers, geologic maps, mineral resources and earthquakes. Explore source records by subject and place.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

170 records · Page 10Linked to original sources

USGS2017

Mineral-deposit model for lithium-cesium-tantalum pegmatites

Lithium-cesium-tantalum (LCT) pegmatites comprise a compositionally defined subset of granitic pegmatites. The major minerals are quartz, potassium feldspar, albite, and muscovite; typical accessory minerals include biotite, garnet, tourmaline, and apatite. The principal lithium ore minerals are spodumene, petalite, and lepidolite; cesium mostly comes from pollucite; and tantalum mostly comes from columbite-tantalite. Tin ore as cassiterite and beryllium ore as beryl also occur in LCT pegmatites, as do a number of gemstones and high-value museum specimens of rare minerals. Individual crystals in LCT pegmatites can be enormous: the largest spodumene was 14 meters long, the largest beryl was 18 meters long, and the largest potassium feldspar was 49 meters long. Lithium-cesium-tantalum pegmatites account for about one-fourth of the world’s lithium production, most of the tantalum production, and all of the cesium production. Giant deposits include Tanco in Canada, Greenbushes in Australia, and Bikita in Zimbabwe. The largest lithium pegmatite in the United States, at King’s Mountain, North Carolina, is no longer being mined although large reserves of lithium remain. Depending on size and attitude of the pegmatite, a variety of mining techniques are used, including artisanal surface mining, open-pit surface mining, small underground workings, and large underground operations using room-and-pillar design. In favorable circumstances, what would otherwise be gangue minerals (quartz, potassium feldspar, albite, and muscovite) can be mined along with lithium and (or) tantalum as coproducts. Most LCT pegmatites are hosted in metamorphosed supracrustal rocks in the upper greenschist to lower amphibolite facies. Lithium-cesium-tantalum pegmatite intrusions generally are emplaced late during orogeny, with emplacement being controlled by pre-existing structures. Typically, they crop out near evolved, peraluminous granites and leucogranites from which they are inferred to be derived by fractional crystallization. In cases where a parental granite pluton is not exposed, one is inferred to lie at depth. Lithium-cesium-tantalum LCT pegmatite melts are enriched in fluxing components including H2O, F, P, and B, which depress the solidus temperature, lower the density, and increase rates of ionic diffusion. This, in turn, enables pegmatites to form thin dikes and massive crystals despite having a felsic composition and temperatures that are significantly lower than ordinary granitic melts. Lithium-cesium-tantalum pegmatites crystallized at remarkably low temperatures (about 350–550 °C) in a remarkably short time (days to years). Lithium-cesium-tantalum pegmatites form in orogenic hinterlands as products of plate convergence. Most formed during collisional orogeny (for example, Kings Mountain district, North Carolina). Specific causes of LCT pegmatite-related magmatism could include: ordinary arc processes; over thickening of continental crust during collision or subduction; slab breakoff during or after collision; slab delamination before, during, or after collision; and late collisional extensional collapse and consequent decompression melting. Lithium-cesium-tantalum pegmatite deposits are present in all continents including Antarctica and in rocks spanning 3 billion years of Earth history. The global age distribution of LCT pegmatites is similar to those of common pegmatites, orogenic granites, and detrital zircons. Peak times of LCT pegmatite genesis at about 2640, 1800, 960, 485, and 310 Ma (million years before present) correspond to times of collisional orogeny and supercontinent assembly. Between these pulses were long intervals when few or no LCT pegmatites formed. These minima overlap with supercontinent tenures at ca. 2450–2225, 1625–1000, 875–725, and 250–200 Ma. Exploration and assessment for LCT pegmatites are guided by a number of observations. In frontier areas where exploration has been minimal at best, the key first-order criteria are an orogenic hinterland setting, appropriate regional metamorphic grades, and the presence of evolved granites and common granitic pegmatites. New LCT pegmatites are most likely to be found near known deposits. Pegmatites tend to show a regional mineralogical and geochemical zoning pattern with respect to the inferred parental granite, with the greatest enrichment in the more distal pegmatites. Mineral-chemical trends in common pegmatites that can point toward an evolved LCT pegmatite include: increasing rubidium in potassium feldspar, increasing lithium in white mica, increasing manganese in garnet, and increasing tantalum and manganese in columbite-tantalite. Most LCT pegmatite bodies show a distinctive internal zonation featuring four zones: border, wall, intermediate (where lithium, cesium, and tantalum are generally concentrated), and core. This zonation is expressed both in cross section and map view; thus, what may appear to be a common pegmatite may instead be the edge of a mineralized body. Neither lithium-cesium-tantalum pegmatites nor their parental granites are likely to cause serious environmental concerns. Soils and country rock surrounding a LCT pegmatite, as well as waste from mining operations, may be enriched in characteristic elements relative to global average soil and bedrock values. These elements may include lithium, cesium, tantalum, beryllium, boron, fluorine, phosphorus, manganese, gallium, rubidium, niobium, tin, and hafnium. Among this suite of elements, however, the only ones that might present a concern for environmental health are beryllium and fluorine, which are included in the U.S. Environmental Protection Agency drinking-water regulations with maximum contaminant levels of 4 micrograms per liter and 4 milligrams per liter, respectively.

Scientific Investigations Report↗

USGS2022

Bedrock geologic map of the Crown Point quadrangle, Essex County, New York, and Addison County, Vermont

The bedrock geology of the 7.5-minute Crown Point quadrangle consists of deformed and metamorphosed Mesoproterozoic gneisses of the Adirondack Highlands unconformably overlain by weakly deformed lower Paleozoic sedimentary rocks of the Champlain Valley. The Mesoproterozoic rocks occur on the eastern edge of the Adirondack Highlands and represent an extension of the Grenville Province of Laurentia. Granulite facies Mesoproterozoic paragneiss, marble, and amphibolite hosted the emplacement of granitic orthogneiss at approximately 1.18–1.15 giga-annum (Ga, billion years before present). The earliest of four phases of deformation (D1) is characterized by gneissosity, rarely preserved F1 isoclinal folds, and migmatite in the host rocks. Subsequent D2 deformation produced a composite penetrative gneissosity, migmatite, and isoclinal F2 folds. Towards the end of D2, felsic magmatism (including the regionally extensive Lyon Mountain Granite Gneiss, abbreviated “LMG”) spread by penetrative migration as semiconcordant alkali feldspar granite sheets subparallel to S2 into previously deformed lithologies. The LMG crystallized at approximately 1.15 Ga and displays synkinematic F2 folds thus constraining the time of D2 deformation. Exhumation during D3 produced F3 folds exhibited in regional domes and basins, such as the Keeney Mountain synform, local reactivation of the S2 foliation, partial melting, metamorphism, metasomatism, iron ore remobilization, and intrusion of magnetite-bearing pegmatite both as layer-parallel sills and crosscutting dikes. D4 created NE- and NW-trending boudinage, local high-grade ductile shear zones, and crosscutting granitic pegmatite dikes. Kilometer (km)-scale lineaments readily observed in lidar data are Ediacaran mafic dikes and Phanerozoic brittle faults. The Paleozoic rocks are part of the Early Cambrian to Late Ordovician great American carbonate bank on the ancient margin of Laurentia. Cambrian-Ordovician stratigraphy records an approximately 1-km-thick section and a transition from synrift clastics to passive margin peritidal carbonate buildups to gradually deeper water subtidal to shelf carbonates during foreland basin development associated with the Taconic orogeny. The Paleozoic rocks are weakly folded and block faulted. Large areas of the Champlain Valley are covered by undifferentiated glacial deposits, some of which contain mapped landslides. The map also shows waste rock piles and tailings from historical mining operations and large areas of artificial fill. This study was undertaken to improve our understanding of the bedrock geology in the Adirondack Highlands, establish a modern framework for 1:24,000-scale bedrock geologic mapping in the Adirondacks, provide a context for historical iron mines in the eastern Adirondacks, and update the stratigraphy of the Champlain Valley in New York and Vermont. This Scientific Investigations Map of the Crown Point 7.5-minute quadrangle consists of a map sheet, an explanatory pamphlet, and a geographic information system database that includes bedrock geologic units, faults, outcrops, and structural geologic information. The map sheet includes a bedrock geologic map, a correlation of map units, a description of map units, an explanation of map symbols, three cross sections, and a simplified surficial geologic map that includes lidar percent slope. The explanatory pamphlet includes a discussion of the geology. The bedrock geologic map on the map sheet is multi-layered and has been designed to enable the user to turn off the surficial map layer to view the concealed bedrock map units.

New York, Vermont↗

USGS2008

Eruptive history and tectonic setting of Medicine Lake Volcano, a large rear-arc volcano in the southern Cascades

Medicine Lake Volcano (MLV), located in the southern Cascades ∼ 55 km east-northeast of contemporaneous Mount Shasta, has been found by exploratory geothermal drilling to have a surprisingly silicic core mantled by mafic lavas. This unexpected result is very different from the long-held view derived from previous mapping of exposed geology that MLV is a dominantly basaltic shield volcano. Detailed mapping shows that < 6% of the ∼ 2000 km 2 of mapped MLV lavas on this southern Cascade Range shield-shaped edifice are rhyolitic and dacitic, but drill holes on the edifice penetrated more than 30% silicic lava. Argon dating yields ages in the range ∼ 475 to 300 ka for early rhyolites. Dates on the stratigraphically lowest mafic lavas at MLV fall into this time frame as well, indicating that volcanism at MLV began about half a million years ago. Mafic compositions apparently did not dominate until ∼ 300 ka. Rhyolite eruptions were scarce post-300 ka until late Holocene time. However, a dacite episode at ∼ 200 to ∼ 180 ka included the volcano's only ash-flow tuff, which was erupted from within the summit caldera. At ∼ 100 ka, compositionally distinctive high-Na andesite and minor dacite built most of the present caldera rim. Eruption of these lavas was followed soon after by several large basalt flows, such that the combined area covered by eruptions between 100 ka and postglacial time amounts to nearly two-thirds of the volcano's area. Postglacial eruptive activity was strongly episodic and also covered a disproportionate amount of area. The volcano has erupted 9 times in the past 5200 years, one of the highest rates of late Holocene eruptive activity in the Cascades. Estimated volume of MLV is ∼ 600 km 3 , giving an overall effusion rate of ∼ 1.2 km 3 per thousand years, although the rate for the past 100 kyr may be only half that. During much of the volcano's history, both dry HAOT (high-alumina olivine tholeiite) and hydrous calcalkaline basalts erupted together in close temporal and spatial proximity. Petrologic studies indicate that the HAOT magmas were derived by dry melting of spinel peridotite mantle near the crust mantle boundary. Subduction-derived H 2 O-rich fluids played an important role in the generation of calcalkaline magmas. Petrology, geochemistry and proximity indicate that MLV is part of the Cascades magmatic arc and not a Basin and Range volcano, although Basin and Range extension impinges on the volcano and strongly influences its eruptive style. MLV may be analogous to Mount Adams in southern Washington, but not, as sometimes proposed, to the older distributed back-arc Simcoe Mountains volcanic field.

Journal of Volcanology and Geothermal Research↗

USGS1986

Effect of height and orientation (microclimate) on geomorphic degradation rates and processes, late-glacial terrace scarps in central Idaho

Skip Nav Destination RESEARCH ARTICLE | JULY 01, 1986 Effect of height and orientation (microclimate) on geomorphic degradation rates and processes, late-glacial terrace scarps in central Idaho KENNETH L. PIERCE ; STEVEN M. COLMAN Author and Article Information GSA Bulletin (1986) 97 (7): 869–885. https://doi.org/10.1130/0016-7606(1986)97<869:EOHAOM>2.0.CO;2 Article history Standard View 2.0.CO;2" data-doctype="contentPdf" data-article-id="188507" data-mce-href="https://pubs.geoscienceworld.org/gsa/gsabulletin/article-pdf/97/7/869/3434736/i0016-7606-97-7-869.pdf"> Open the PDF for in another window Cite Share Icon Share Permissions Abstract Terrace scarps can serve as a nearly ideal natural laboratory for the study of the evolution of slopes. This paper examines the effects of scarp size (height) and orientation (microclimate) by keeping constant variables such as age, lithology, and regional climate. If a scarp degrades as a closed system, and downslope movement is directly proportional to surface gradient , the evolution of the scarp is modeled by the diffusion equation. For a group of scarps of same age and known starting angle, the diffusion-equation model predicts the relation between maximum scarp angle (𝛉) and scarp height ( h ). Late Pleistocene terrace scarps now as steep as 33.25°, as well as measured angles of repose for sand and gravel, require a starting angle as steep as 33.5°. For latest Pleistocene Idaho and Utah scarps, as h increases, 𝛉 is gentler (more degraded) than modeled by the diffusion equation with a constant rate coefficient. The degradation-rate coefficient ( c ) increases tenfold with scarp height; it should not change with scarp height if downslope movement is solely determined by surface gradient (to the first power). Soil wash appears to be responsible for this departure from the diffusion-equation model, for transport rate by soil wash is a function of scarp size (height). South-facing scarps are less vegetated and more degraded than north-facing scarps. For scarps 2 m high, the degradation rate ( c *) on S-facing scarps is 2 times that on N-facing scarps; for 10-m scarps, it is 5 times. The observed dependence of the rate coefficient c * on scarp height can be removed by normalizing c * to values for west-facing scarps of the same height. The residual c * values calculated by this method correlate well with differences in incident solar radiation resulting from the different scarp orientations and maximum gradients. This correlation demonstrates the importance of orientation on slope processes and their rates through the differences in freeze-thaw cycles, soil moisture, and vegetative cover. Scarp morphology may be used to estimate age, if one accounts for the effects of climate and for scarp height, orientation, and lithology. For example, using the dated Bonneville shoreline scarps for calibration and comparing only scarps of equal height, we estimate the Drum Mountains fault scarps to be 9,000 yr old. This age is about twice that produced by previous diffusion-equation calculations that have not accounted for the height as we have here, but it is the same as independent geologic estimates of their age.

Idaho↗

USGS2013

A natural resource condition assessment for Sequoia and Kings Canyon National Parks: Appendix 14: plants of conservation concern

Sequoia and Kings Canyon National Parks are located in the California Floristic Province, which has been named one of world‘s hotspots of endemic biodiversity (Myers et al. 2000). The California Floristic Province is the largest and most important geographic floristic unit in California and extends from the Klamath Mountains of southwestern Oregon to the northwestern portion of Baja California (Hickman 1993). The Sierra Nevada, one of six regions that make up the California Floristic Province, covers nearly 20% of the land in California yet contains over 50% of its flora. Within the Sierra Nevada, the southern Sierra supports more Sierran endemic and rare plant taxa than the central and northern portions of the region (Shevock 1996). Sequoia and Kings Canyon National Parks (SEKI) encompass roughly 20% of the southern Sierra Nevada region. The parks overlap three floristic subregions (central Sierra Nevada High, southern Sierra Nevada High, and southern Sierra Nevada Foothills), and border the Great Basin Floristic Province. The parks support a rich and diverse vascular flora composed of over 1,560 taxa. Of these, 150 taxa are identified as having special status. The term special status is applied here to include taxa that are state or federally listed, rare in California, or at risk because they have a limited distribution. Only one species from these parks is listed under the state or federal Endangered Species Acts ( Carex tompkinsii , Tompkins‘ sedge, is listed as a rare species under the California Endangered Species Act), and one species is under review for federal endangered listing ( Pinus albicaulis , whitebark pine). However, an absence of threatened and endangered species recognized by Endangered Species Acts is not equivalent to an absence of species at risk. There are 83 plant taxa documented as occurring in SEKI that are considered imperiled or vulnerable in the state by the California Department of Fish and Game‘s California Natural Diversity Database (CNDDB 2010a). There are an additional 66 taxa not formally listed by CNDDB that are recognized as having special status because their distribution is restricted to the Sierra Nevada. Special status plants are distributed throughout the two parks and inhabit a wide range of environments along the length of the elevation gradient that characterizes these parks. Ideally, we would assess the condition (status and trends) of each of the taxa on the SEKI special status plant list, documenting current population sizes, demographic rates and demographic trends. We would also hope to quantify the effects of individual stressors on each species based on existing monitoring and research. However, no data are available for most of the species on the special status plant list. For those few species (12 herbaceous species and two tree species) for which we possess some change over time information, the data are not adequate to make a competent assessment. Note that we have not explored the tree demographic information in any detail, as is covered in the NRCA Intact Forest/Five Needle Pines and Sequoia chapters. In general, we are unable to present an ‗integrity‘ metric for special status species in the parks, since the data to quantify the condition of each species in such a manner is not available. In contrast, the park does possess substantial data describing biodiversity in the parks. Therefore, our analysis focuses on describing the distribution and rarity of special status plants within the parks, with a particular focus on assessing the spatial distribution of species richness. We hope that such information will prove useful to park managers in determining which areas in the parks merit the most attention (for example in developing monitoring protocols). We also assess potential vulnerability of special status species to the stressors chosen by the NRCA working group, using both park data and available literature. As a first step, we spent considerable effort updating and refining the criteria for the special status plant list, as this list defines which taxa are considered in our assessment. Observation data of these species was then compiled from all known sources in order to provide a comprehensive view of where special status plants have been documented and, ultimately, to enable the most informed determinations of areas in the parks that potentially support the highest number of rare and endemic taxa. These ‗hot spot‘ analyses are presented by geographic region, vegetation type and elevation. For these and other analyses presented in this report, we place more focus on summarizing findings for the herbaceous and shrub special status taxa than on special status trees. The trees which qualify as special status are the focus of other NRCA chapters, including Giant Sequoia and Intact Forests/Five-needle Pines. We do, however, present their mapped distributions and provide overviews of research related to the special status tree taxa in the Stressors section of this report.

California↗

USGS1971

Structure and petrology of the alpine-type peridotite at Burro Mountain, California, U.S.A.

The alpine-type peridotite at Burro Mountain is a partially serpentinized harzburgite-dunite body approximately 2 km in diameter. It lies in a chaotic mélange derived from the Franciscan Formation (Upper Jurassic to Upper Cretaceous) of the southern Coast Ranges of California. The peridotite is bounded on the east by a vertical fault in the Nacimiento fault zone that brings sedimentary rocks of Taliaferro's (1943 b ) Asuncion Group (Upper Cretaceous) into contact with the peridotite. The peridotite appears to be one of a number of tectonic lenses, having a wide range in size, that make up the mélange . These lenses include metagraywacke, metachert, greenstone, amphibolite, and blueschist, as well as ultramafic rocks, and represent a wide range of pressure-temperature environments. The outer shell of the peridotite is a sheared serpentinite zone 10–15 m thick. The peridotite was tectonically emplaced at its present level as a cold solid mass and had little effect on the mineral assemblages of the Franciscan Formation. Local development of lawsonite and aragonite in shear zones may be related to the peridotite emplacement. Foliated harzburgite forms approximately 60 per cent of the peridotite. It is a lithologically uniform rock that has an olivine: orthopyroxene ratio of approximately 75:25. Accessory clinopyroxene and chromian spinel generally make up less than 5 per cent of the harzburgite. Dunite, composed of olivine, accessory chromian spinel (< 5 per cent), and trace amounts of pyroxene, makes up approximately 40 per cent of the peridotite and occurs as dikes, sills, and irregular bodies in the harzburgite. Olivine and pyroxene show small but significant compositional variations and chromian spinel shows a large range in the cation ratio Cr/(Cr+Al+ Fe 3+ ). The compositional variations in these minerals are related to original differences in bulk chemical composition. The following compositional ranges were determined for minerals in the harzburgite: olivine, Fo 91.1 −Fo 91.4 ; orthopyroxene, En 89.8 −En 91.1 ; clinopyroxene, Ca 47.0 Mg 50.0 Fe 3.0 −Ca 48.7 Mg 48.2 Fe 3.1 ; chromian spinel, Cr/(Cr+Al+Fe 3+ ) 0.37−0.55. The pyroxenes have a range in A1 2 O 3 content of 1.3−3.0 wt per cent. Olivine from dunite ranges from Fo 91 to Fo 92 7 and the chromian spinel has a range in the Cr/(Cr+Al+Fe 3+ ) ratio of 0.30−0.75. Although all the dunites are lithologically similar, three distinct types are recognized on the basis of composition of coexisting olivine and chromian spinel. Structural relations between the three types of dunite suggest three periods of emplacement (possibly overlapping) of dunite into harzburgite. The evidence indicates that the dunite, and probably also the harzburgite crystallized from an ultramafic magma, probably in the upper mantle. After the magmatic episode and crystallization, the peridotite was subjected to a deep-seated plastic deformation and recrystallization. The first phase of the deformation produced a pervasive, planar structural element (S 1 ) that crosscuts many harzburgite-dunite contacts. It is probable that some of the dunite sills were emplaced during this deformation. The foliation, S 1 , is defined by layers of different orthopyroxene content in harzburgite, and by discontinuous layers of chromian spinel in dunite. Flow or slip along S 1 produced slip folds in harzburgite—dunite contacts with axial planes parallel to S 1 . At a later stage, isoclinal folds developed in S 1 , and the present olivine microfabric was probably formed by recrystallization in the stress field that produced the isoclinal folding. In the olivine microfabric, X tends to be perpendicular to the axial planes (S 2 ) of the isoclinal folds and Y and Z tend to form double maxima in S 2 approximately 90° apart. Mg−Fe 2+ distribution between coexisting mineral pairs yields a calculated temperature of formation of approximately 1200 °C. Although this temperature is only a nominal value, it indicates that the mineral pairs equilibrated at a significantly high temperature. In view of the deformation and recrystallization, the calculated temperature possibly represents subsolidus re-equilibration of the minerals during this event. The deformation and recrystallization probably occurred shortly after crystallization while the peridotite was still at a high temperature. A later deep-seated deformation produced small scattered kink folds in S 1 that tend to disrupt the major olivine microfabric. The kink folding was accompanied or followed by the development of kink bands in olivine that reflect intragranular gliding on the system T = [ Okl ], t = [100]. The kink bands probably formed at a minimum temperature of 1000 °C. Following the deep-seated deformation, which probably took place in the mantle, the peridotite mass was tectonically detached and moved upward to its present level in the crust. Cleavages, joints, and faults provided channels for water to pervade the peridotite and allow alteration of the primary minerals.

California↗

USGS1976

Summary appraisals of the nation's ground-water resources – California region

Most people in the California Region live in a semiarid or arid climate, with precipitation less than the potential evapotranspiration- environments of perennial water deficiency. The deficiency becomes most onerous during the characteristically rainless summers and during recurrent droughts that may continue for 10--20 years. However, water from winter rain and snow can be stored for use during the dry summer months, and water stored during a wet climatic period can be used in a succeeding dry period; moreover, perennial deficiency can be overcome by bringing water from areas of perennial surplus. Ground-water reservoirs have especial significance in arid and semiarid regions as repositories where water is stored or can be stored with minimum loss by evaporation. Nearly all the ground-water reservoirs of the California Region are in alluvial sediments of valleys and plains that flank the mountain ranges. The largest, underlying the vast Central Valley, occupies 10 percent of the area of the region, has an estimated usable capacity exceeding 100 million acre-feet (125 cubic kilometres), and has an annual pumpage from wells of about 13 million acre-feet (16 cubic kilometres). Another 10 percent of the region is occupied by 55 developed ground-water reservoirs that are widely distributed; aggregate annual pumpage from them is about 3 1/2 million acre-feet (4 cubic kilometres). In the southeastern desert about 60 ground-water reservoirs occupy still another 10 percent of the region; these have been explored only enough to show that most have some usable water, but current use is negligible. In northeastern California and adjacent Oregon and Nevada, ground-water reservoirs are identified only in valleys and lowlands where wells are feasible, but basaltic rocks of the Cascade Range and Modoc Plateau are excellent aquifers distributed over an area constituting about 15 percent of the region. In sum, slightly less than half the California Region is underlain by ground-water reservoirs, either in valley fill or in volcanic rocks, which can yield significant quantities of water to wells. The rest of the California Region includes the mountains, canyons, slopes, and foothills of the Sierra Nevada, Coast Ranges, and Basin Ranges, whose consolidated rocks and products of their weathering may be permeable locally but are not generally so. Here, the prevailing method of ground-water development is still mostly trial and error, and while in many places a well can yield enough water for a family, some families might have to do without amenities such as flush toilets and automatic washers. For more than half a century the California Region has led all others in North America in pumping of ground water as well as in the area, variety, yield, and export of crops irrigated by water from wells. It has led in the development and use of deep-well turbine pumps for large yield and in the drilling of water wells to great depths. At the same time, such developments have resulted in the elimination of artesian pressures that produced thousands of flowing wells in the 19th century and led to the wide distribution of "falling water tables." Also, California was first to induce encroachment of seawater into wells (in 1906); first to recognize subsidence of land caused by pumping from wells (in 1933), generating news about land sinking in San Jose, Long Beach, and along the Delta-Mendota and Friant-Kern Canals; and first to experience pollution of ground-water reservoirs by brines, chemicals, industrial wastes, and petroleum byproducts including gasoline. The region has led in research in several fields leading to solution of many of these problems. Ground-water problems developed rapidly after World War II with booming population, agriculture, industry, and water demand during several years of regionwide drought. Water levels in wells trended downward almost everywhere as a natural effect of the drought and at accelerated rates in areas of pumping for new enterprises or to supplement subnormal surface-water supplies. The declines in many pumping areas exceeded 100 feet (30 metres), and in some confined aquifers the potentiometric surface was drawn down more than 330 feet (100 metres). The depletion of ground-water storage has had "permanent" side effects, including subsidence of the land exceeding 10 feet (3 metres) in extensive areas, and seawater intrusion that ended the useful lives of many wells along the coast and as much as 6 miles (10 kilometres) inland. Some problems have been solved, but these solutions have at times created other problems. Many ground-water reservoirs have gone through one or more stages - exploration for productive aquifers, exploitation and development for use of the water, restriction to the perennial supply or "safe" yield, importation of surface water, artificial recharge of ground water, conjunctive use of surface and ground water, protection of water quality, and integrated management of use and disposal of water. This evolutionary sequence is unique for each reservoir, and so generalizations become difficult in a regional appraisal; also, the changes with time are significant and varied, and knowledge of prior events is a prerequisite in an appraisal of the resource in a specific year. As of 1970, water levels in many wells had risen significantly from the minimum levels of record reached during the 1960's or earlier; only in areas of new development and in desert areas of long-continued "mining" of nonreplenished water was ground-water storage still being depleted. Land subsidence has continued at diminishing rates and practically has come to a halt in some areas; invading seawater has been stopped or nudged back in most places where problems were significant. The current, favorable situation has been helped by climatic variations, from drought in 1945-52 and exceedingly dry years in 1959 and 1961 to above-normal precipitation in 1969 and 1970; but most of the serious problems have been solved by human efforts, including especially the implementation of the California Water Plan, transporting water from areas of perennial surplus to areas where it is used in lieu of ground water or for ground-water replenishment. All major urban areas now import water to supplement or replace the water pumped from wells. Extensive agricultural areas that formerly were irrigated solely by ground water now obtain some of their water from surface reservoirs and canals, especially in the Central Valley. With surface water available as an alternative supply, well owners can view their ground water with complacency. But complacency can lead to neglect and carelessness and consequent deterioration of the ground-water resource by pollution. Claiming heritage from the English Common Law, the existing California law grants to the landowner (riparian) and private enterprise (appropriator) rights to the water stored in ground-water reservoirs or discharged from them, including the base flow of streams. Ground-water development has been by private and local enterprise, and the California legislature has protected and encouraged local responsibility, control, and management of ground water. As to surface water, a constitutional amendment in 1928 limited riparian rights to the quantities of water that were "reasonably required for the beneficial use to be served." The surpluses have become public waters which are collected, stored, transported, and delivered under various contracts by Federal, State, and other agencies. The agencies have not stored water underground because of uncertainty as to their rights, but some local agencies have been encouraged with favorable pricing schedules to undertake the artificial recharge and management of ground-water reservoirs. Thus, conjunctive use of surface and ground water has become a matter of interagency negotiation. Of all the constraints on effective use of ground-water reservoirs, the most formidable may be the attitudes of people. Assurance of water supply is vital in areas of water deficiency, and Government has assumed increasing responsibility for the welfare of people in these areas. Unfortunately, when Government provides this assurance, most beneficiaries demand continued subsidy to the exclusion of perhaps cheaper private development. Indeed, as the water resources are presently segregated-with private rights predominant in ground water and public interest dominant in surface water-ground-water development has suffered for lack of public concern. The region has the scientific and technologic capability for effective use of groundwater reservoirs, as shown by the achievements and programs of several districts, but many districts are not organized or staffed for such comprehensive management and will need assistance and scientific expertise available from State and Federal agencies. Those agencies, in turn, may not have the scientific data that are essential to prevent haphazard activities and to enable programs to be organized for the most effective and attractive utilization of the water resources. In these days of increasing concern over pollution, existing data are generally inadequate to assess the natural deterioration of ground waters as a basis for defining pollution.

California↗

USGS2006

Evaluation of some software measuring displacements using GPS in real-time

For the past decade, the USGS has been monitoring deformation at various locations in the western United States using continuous GPS. The main focus of these measurements are estimates of displacement averaged over one day. Essentially, these consist of recording at 30 seconds intervals the carrier-frequency phase-data (equivalent to travel-time) between a GPS receiver and the GPS satellite network. In turn, these observations, which are converted to pseudo—ranges, are processed using one of the “research grade” programs (GIPSY, Zumberge et al., or GAMIT, wwwgpsg.mit.edu/~simon/gtgk) to estimate the position of the GPS receiver averaged over 24 hours. However, it is possible and desirable to estimate the position of the receiver (actually the antenna) more frequently and to do this within a few seconds of the time actual measurement (known as real-time). A recent example, the 2004 Magnitude 6, Parkfield, California earthquake, demonstrated that having GPS estimates of position more frequently than simply a daily average is required if one requires discrimination between co-seismic and post-seismic deformation (Langbein et al., 2006). The high-rate estimates of position obtained at Parkfield show that post-seismic deformation started less than one-hour after the mainshock and that this deformation was roughly the same magnitude as the co-seismic deformation. The high-rate solutions for Parkfield were done by others including Yehuda Bock at UCSD and Kristine Larson at U. of Colorado, but not the USGS. The Parkfield experience points out the need for an in-house capability by the USGS to be able to accurately measure co-seismic displacements and other rapid, deformation signals using GPS. This applies to both the Earthquake and Volcano Hazard programs. Although at many locations where we monitor deformation, we have strainmeters and tiltmeters in addition to GPS which, in principle, are far more sensitive to rapid deformation over periods of less than a day (Langbein and Bock, 2004). But, not all locales include strain and tiltmeters. Thus, having the capability to extract signals with periods of less than a day is desirable since the distribution of GPS is more extensive than strain and tilt. At both Parkfield and Long Valley, the USGS has been using other software packages to process the GPS data at sub-daily intervals and in real-time. The underlying goal of these types of measurements is to detect any deformation event as it evolves; the 24 hour processing might not provide timely results if such a deformation event is precursory to a geologic hazard (an earthquake for Parkfield and either a volcanic event or an earthquake for Long Valley). In Long Valley, We use the software package called 3DTracker (http://www.3dtracker.com, http://www.condorearth.com) to estimate the changes of in position of a remote site relative to a “fixed” site. The 3DTracker software uses double difference GPS code measurements and receiversatellite-time triple differences from one epoch to the next of the GPS phase data (a proxy for travel-time measurements) and employs a Kalman filter to obtain stability in the estimate of position. That is, the estimate of the current position depends upon the estimate of the prior position. Hence, a time series of position looks fairly smooth depending upon the coefficient selected for the Kalman filter. With triple differences, the sometimes troublesome initial integer cycle ambiguity terms cancel (number of wavelengths between the receiver and each satellite), but only the incremental change in position is calculated. This triple difference Kalman filter solution is slow to converge and less accurate than a double difference (e.g., RTD, Track) solution, but it is robust and computationally efficient (Remondi and Brown, 2000). 3D-Tracker allows use of various single-frequency and dual-frequency GPS phase and code observables including the ionospheric-free combinations (known as LC or L3 and P(L3)) formed from an linear combination of the L1 and L2 carrier phase and code data. The lowest noise observable is the L1 carrier, but it is biased by ionospheric refraction that has amplitudes of about 1 to 10 ppm. This results in a systematic scale error in the relative positions. The L3 phase noise is about 3 times greater than the L1 phase noise, but it is generally used to solve for all but the shortest baselines (< 5 km). In addition, the software does output the position changes is a standard format that can be used for other analysis. At Parkfield, we use the software package called RTD (http://www.geodetics.com). The RTD software has been described in the literature (Bock et al., 2000) but basically, it estimates the position without the constraint of a Kalman filter. It uses double differences (in our studies the LC or ionospheric free observable is used) and the integer ambiguities are resolved independently for each 1-second measurement; Most GPS software that use double-differences require several epochs of measurements to resolve the integer ambiguities. The data files use a proprietary format and can not be read by me or others; rather, Yehuda Bock at UCSD (and author of RTD) translates these files into a standard format that can be read by me. Recently, Tom Herring of MIT has modified the GAMIT software to process kinematically GPS data (www-gpsg.mit.edu/~simon/gtgk/tutorial/Lecture_13.pdf). At this time, the software, known as TRACK, does not process the observations in real-time. Consequently, the latency between the time of the observation and the time when a position estimate is available depends upon the frequency that the data are downloaded and the speed of actually processing the observations; there could be a delay of an hour or two before the a position estimates are available. Unlike RTD and 3DTracker, TRACK comes with GAMIT (which is distributed freely) and is currently operating in a test mode at the USGS office in Pasadena. The LC or ionosphere free observable is used in our TRACK solutions. JPL has a version of their GIPSY software called “Real-time GIPSY (RTG)” (gipsy.jpl.nasa.gov/orms/rtg), which, like TRACK, can process the pseudo-range data “off—line”. However, this software is not freely distributed. Instead, at least one company, NAVCOM, has teamed with JPL to integrate RTG with GPS receivers and telemetry that yields positions in realtime. Kristine Larson of University of Colorado has modified the original GIPSY to estimate positions kinematically. Again, like TRACK, the positions are estimated off—line. Much of her research is described in Larson et al. (2003), and Choi et al. (2004). For Long Valley, out of the 17 GPS sites, we monitor 5 baselines within the caldera at 5 second intervals relative to the Bald Mountain site at the edge of the caldera using 3DTracker. The baseline measurement using 3DTracker consists of determination of the 3 dimensional positions of the 5 remote points (GPS receivers) relative to a GPS site at Bald. A second, independent system collects and downloads once a day the 30-second data used for the 24-hour solutions for the 12 sites not monitored with 3DTracker. For the sites monitored with 3DTracker, the pseudo—range data are decimated to 30 seconds and converted to a form used for the 24-hour solutions. Both sets of telemetry employ 900 MHz spread spectrum radios which require line of site between all of the links. The telemetry for the 3DTracker sites require a dedicated radios at each end and intermediate repeaters as needed, while the telemetry required for the other sites use a single master radio, repeaters as needed, and a radio at each remote site. (The 5 sites being monitored with 3DTracker require 13 radios.) At Parkfield, RTD is used to measure the position changes all 12 baselines at 1 second intervals relative to a site, Pomm, adjacent to the San Andreas Fault. The complete RTD package (hardware and software) collects all of the data and determines the position of each site relative to Pomm. In addition, the system stores both the 1-second and 30-second pseudo-range data for later downloading which are ultimately used in the 24-hour solutions. To do this, each site has a 2.4 GHz radio and a telemetry buffer. The telemetry buffer holds 24-hours of data (in the event that the telemetry link is broken) and converts the RS232 data stream from the GPS receiver into a form compatible with an IP (Internet protocol) network connection. In contrast with the Long Valley system, the telemetry link for GPS at Parkfield consists of a single radio at each remote sites and a single radio at the central site. Although position estimates are produced within 1-second of the observations, these results are not immediately available because there is no high speed Internet connection to Parkfield. Instead, the data are stored on a removable disk and sent to UCSD once per month. Below, I describe the results of a simple experiment to examine the response of some of these systems to simulated deformation that could be an analogue of a tectonic or volcanic event. In many engineering applications, the system response is tested by inputting a step to the system and measuring the output of the system. Essentially, this is what I've done. The experiment described below moves the GPS antenna from its original position to a new position within 1 second; the software tracks the translation. These measurements were conducted in August 2004 with the RTD software at Parkfield, and twice in Long Valley. The first Long Valley test was conducted in September 2004 using 3DTracker on a single baseline. The test was repeated in September 2005 using 3DTracker on two baselines and, importantly, saving the RINEX files of the data so that the data could be replayed through 3DTracker using other options in the program and, using other software packages including TRACK. In addition, we observed a short-term event at the Three Sisters volcano in Oregon. This event was snow melt at a remote GPS site which gave an apparent 15 cm displacement in vertical in less than one-day. 3DTracker is used to monitor this site, and the event was captured with this software. In addition, with the assistance of others, I got additional estimates of position using other software packages; those results are presented. Finally, the precision of both 3DTracker and RTD are compared using a power spectrum. Those results would suggest that 3DTracker using appropriate Kalman filter coefficients would have better precision than RTD; instead, the lower noise level from 3DTracker is a result of smoothing from the Kalman filter. Given the results described in this report, high-rate GPS is certainly capable of accurately measuring displacements of 1 centimeter with a high degree of statistical confidence. Plotting these results show that the time of the displacement can be visually determined to that of the sampling interval of the data. However, especially with small amplitude signals, any of the software packages can yield erroneous deformation “signals” that are either due excess travel-time of the GPS carrier frequency from multipath or a limitation in the software. Thus, the time series of displacements must be viewed with caution and knowledge of external circumstances that might cause a change in position. The casual reader should continue with the next section describing the methods then jump to the last two sections for the discussion and conclusions. I have made some recommendations there.

Open-File Report↗