USGS ScienceSearch

SEARCH · USGS Science

Results for “Journal of Geophysical Research - Atmospheres”

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

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

At least 379 records · Page 21Linked to original sources

Thermal signature, eruption style, and eruption evolution at Pele and Pillan on Io

The Galileo spacecraft has been periodically monitoring volcanic activity on Io since June 1996, making it possible to chart the evolution of individual eruptions. We present results of coanalysis of Near-Infrared Mapping Spectrometer (NIMS) and solid-state imaging (SSI) data of eruptions at Pele and Pillan, especially from a particularly illuminating data set consisting of mutually constraining, near-simultaneous NIMS and SSI observations obtained during orbit C9 in June 1997. The observed thermal signature from each hot spot, and the way in which the thermal signature changes with time, tightly constrains the possible styles of eruption. Pele and Pillan have very different eruption styles. From September 1996 through May 1999, Pele demonstrates an almost constant total thermal output, with thermal emission spectra indicative of a long-lived, active lava lake. The NIMS Pillan data exhibit the thermal signature of a “Pillanian” eruption style, a large, vigorous eruption with associated open channel, or sheet flows, producing an extensive flow field by orbit C10 in September 1997. The high mass eruption rate, high liquidus temperature (at least 1870 K) eruption at Pillan is the best candidate so far for an active ultramafic (magnesium-rich, “komatiitic”) flow on Io, a style of eruption never before witnessed. The thermal output per unit area from Pillan is, however, consistent with the emplacement of large, open-channel flows. Magma temperature at Pele is ≥1600 K. If the magma temperature is 1600 K, it suggests a komatiitic-basalt composition. The power output from Pele is indicative of a magma volumetric eruption rate of ∼250 to 340 m 3 s −1 . Although the Pele lava lake is considerably larger than its terrestrial counterparts, the power and mass fluxes per unit area are similar to active terrestrial lava lakes.

Journal of Geophysical Research E: Planets

Imaging of volcanic activity on Jupiter's moon Io by Galileo during the Galileo Europa Mission and the Galileo Millennium Mission

The Solid-State Imaging (SSI) instrument provided the first high- and medium-resolution views of Io as the Galileo spacecraft closed in on the volcanic body in late 1999 and early 2000. While each volcanic center has many unique features, the majority can be placed into one of two broad categories. The “Promethean” eruptions, typified by the volcanic center Prometheus, are characterized by long-lived steady eruptions producing a compound flow field emplaced in an insulating manner over a period of years to decades. In contrast, “Pillanian” eruptions are characterized by large pyroclastic deposits and short-lived but high effusion rate eruptions from fissures feeding open-channel or open-sheet flows. Both types of eruptions commonly have ∼100-km-tall, bright, SO 2 -rich plumes forming near the flow fronts and smaller deposits of red material that mark the vent for the silicate lavas.

Journal of Geophysical Research E: Planets

Rock types of South Pole-Aitken basin and extent of basaltic volcanism

The enormous pre-Nectarian South Pole-Aitken (SPA) basin represents a geophysically and compositionally unique region on the Moon. We present and analyze the mineralogical diversity across this basin and discuss the implications for basin evolution. Rock types are derived from Clementine multispectral data based on diagnostic characteristics of ferrous absorptions in fresh materials. Individual areas are characterized as noritic (dominated by low-Ca pyroxene), gabbroic/basaltic (dominated by high-Ca pyroxene), feldspathic (<3-6% FeO), and olivine-gabbro (dominated by high-Ca pyroxene and olivine). The anorthositic crust has effectively been removed from the interior of the basin. The style of volcanism within the basin extends over several 100 Myr and includes mare basalt and pyroclastic deposits. Several areas of ancient (pre-Orientale) volcanism, or cryptomaria, have also been identified. The nonmare mafic lithology that occurs across the basin is shown to be noritic in composition and is pervasive laterally and vertically. We interpret this to represent impact melt/breccia deposits derived from the lower crust. A few localized areas are identified within the basin that contain more diverse lithologies (gabbro, olivine-gabbro), some of which may represent material from the deepest part of the lower crust and perhaps uppermost mantle involved in the SPA event.

Journal of Geophysical Research E: Planets

The roughness of natural terrain: A planetary and remote sensing perspective

We examine the various methods and parameters in common use for quantifying and reporting surface topographic "roughness." It is shown that scale-dependent roughness parameters are almost always required, though not widely used. We suggest a method of standardizing the parameters that are computed and reported so that topographic data gathered by different workers using different field techniques can be directly and easily intercompared. We illustrate the proposed method by analyzing topographic data from 60 different surfaces gathered by five different groups and examine the information for common features. We briefly discuss the implications of our analysis for studies of planetary surface roughness, lander safety, and radar remote sensing modeling and analysis.

Journal of Geophysical Research E: Planets

Lunar prospector epithermal neutrons from impact craters and landing sites: Implications for surface maturity and hydrogen distribution

Initial studies of neutron spectrometer data returned by Lunar Prospector concentrated on the discovery of enhanced hydrogen abundances near both lunar poles. However, the nonpolar data exhibit intriguing patterns that appear spatially correlated with surface features such as young impact craters (e.g., Tycho). Such immature crater materials may have low hydrogen contents because of their relative lack of exposure to solar wind-implanted volatiles. We tested this hypothesis by comparing epithermal* neutron counts (i.e., epithermal −0.057 × thermal neutrons) for Copernican-age craters classified as relatively young, intermediate, and old (as determined by previous studies of Clementine optical maturity variations). The epithermal* counts of the crater and continuous ejecta regions suggest that the youngest impact materials are relatively devoid of hydrogen in the upper 1 m of regolith. We also show that the mean hydrogen contents measured in Apollo and Luna landing site samples are only moderately well correlated to the epithermal* neutron counts at the landing sites, likely owing to the effects of rare earth elements. These results suggest that further work is required to define better how hydrogen distribution can be revealed by epithermal neutrons in order to understand more fully the nature and sources (e.g., solar wind, meteorite impacts) of volatiles in the lunar regolith.

Journal of Geophysical Research E: Planets

Io in the near infrared: Near-Infrared Mapping Spectrometer (NIMS) results from the Galileo flybys in 1999 and 2000

Galileo's Near-Infrared Mapping Spectrometer (NIMS) observed Io during the spacecraft's three flybys in October 1999, November 1999, and February 2000. The observations, which are summarized here, were used to map the detailed thermal structure of active volcanic regions and the surface distribution of SO 2 and to investigate the origin of a yet unidentified compound showing an absorption feature at ∼1 μm. We present a summary of the observations and results, focusing on the distribution of thermal emission and of SO 2 deposits. We find high eruption temperatures, consistent with ultramafic volcanism, at Pele. Such temperatures may be present at other hot spots, but the hottest areas may be too small for those temperatures to be detected at the spatial resolution of our observations. Loki is the site of frequent eruptions, and the low thermal emission may represent lavas cooling on the caldera's surface or the cooling crust of a lava lake. High-resolution spectral observations of Emakong caldera show thermal emission and SO 2 within the same pixels, implying that patches of SO 2 frost and patches of cooling lavas or sulfur flows are present within a few kilometers from one another. Thermal maps of Prometheus and Amirani show that these two hot spots are characterized by long lava flows. The thermal profiles of flows at both locations are consistent with insulated flows, with the Amirani flow field having more breakouts of fresh lava along its length. Prometheus and Amirani each show a white ring at visible wavelengths, while SO 2 distribution maps show that the highest concentration of SO 2 in both ring deposits lies outside the white portion. Visible measurements at high phase angles show that the white deposit around Prometheus extends into the SO 2 ring. This suggests that the deposits are thin and that compositional or grain size variations may occur in the radial direction. SO 2 mapping of the Chaac region shows that the interior of a caldera adjacent to Chaac has almost pure SO 2 . The deposit appears to be topographically controlled, suggesting a possible origin by liquid flow.

Journal of Geophysical Research E: Planets

Ancient drainage basin of the Tharsis region, Mars: Potential source for outflow channel systems and putative oceans or paleolakes

Paleotopographic reconstructions based on a synthesis of published geologic information and high-resolution topography, including topographic profiles, reveal the potential existence of an enormous drainage basin/aquifer system in the eastern part of the Tharsis region during the Noachian Period. Large topographic highs formed the margin of the gigantic drainage basin. Subsequently, lavas, sediments, and volatiles partly infilled the basin, resulting in an enormous and productive regional aquifer. The stacked sequences of water-bearing strata were then deformed locally and, in places, exposed by magmatic-driven uplifts, tectonic deformation, and erosion. This basin model provides a potential source of water necessary to carve the large outflow channel systems of the Tharsis and surrounding regions and to contribute to the formation of putative northern-plains ocean(s) and/or paleolakes.

Journal of Geophysical Research E: Planets

Remotely sensed geology from lander-based to orbital perspectives: Results of FIDO rover May 2000 field tests

Blind field tests of the Field Integration Design and Operations (FIDO) prototype Mars rover were carried out 7-16 May 2000. A Core Operations Team (COT), sequestered at the Jet Propulsion Laboratory without knowledge of test site location, prepared command sequences and interpreted data acquired by the rover. Instrument sensors included a stereo panoramic camera, navigational and hazard-avoidance cameras, a color microscopic imager, an infrared point spectrometer, and a rock coring drill. The COT designed command sequences, which were relayed by satellite uplink to the rover, and evaluated instrument data. Using aerial photos and Airborne Visible and Infrared Imaging Spectrometer (AVIRIS) data, and information from the rover sensors, the COT inferred the geology of the landing site during the 18 sol mission, including lithologic diversity, stratigraphic relationships, environments of deposition, and weathering characteristics. Prominent lithologic units were interpreted to be dolomite-bearing rocks, kaolinite-bearing altered felsic volcanic materials, and basalt. The color panoramic camera revealed sedimentary layering and rock textures, and geologic relationships seen in rock exposures. The infrared point spectrometer permitted identification of prominent carbonate and kaolinite spectral features and permitted correlations to outcrops that could not be reached by the rover. The color microscopic imager revealed fine-scale rock textures, soil components, and results of coring experiments. Test results show that close-up interrogation of rocks is essential to investigations of geologic environments and that observations must include scales ranging from individual boulders and outcrops (microscopic, macroscopic) to orbital remote sensing, with sufficient intermediate steps (descent images) to connect in situ and remote observations.

Journal of Geophysical Research E: Planets

Upper crustal structure in Puget Lowland, Washington: Results from the 1998 Seismic Hazards Investigation in Puget Sound

A new three-dimensional (3-D) model shows seismic velocities beneath the Puget Lowland to a depth of 11 km. The model is based on a tomographic inversion of nearly one million first-arrival travel times recorded during the 1998 Seismic Hazards Investigation in Puget Sound (SHIPS), allowing higher-resolution mapping of subsurface structures than previously possible. The model allows us to refine the subsurface geometry of previously proposed faults (e.g., Seattle, Hood Canal, southern Whidbey Island, and Devils Mountain fault zones) as well as to identify structures (Tacoma, Lofall, and Sequim fault zones) that warrant additional study. The largest and most important of these newly identified structures lies along the northern boundary of the Tacoma basin; we informally refer to this structure here as the Tacoma fault zone. Although tomography cannot provide information on the recency of motion on any structure, Holocene earthquake activity on the Tacoma fault zone is suggested by seismicity along it and paleoseismic evidence for abrupt uplift of tidal marsh deposits to its north. The tomography reveals four large, west to northwest trending low-velocity basins (Tacoma, Seattle, Everett, and Port Townsend) separated by regions of higher velocity ridges that are coincident with fault-bounded uplifts of Eocene Crescent Formation basalt and pre-Tertiary basement. The shapes of the basins and uplifts are similar to those observed in gravity data; gravity anomalies calculated from the 3-D tomography model are in close agreement with the observed anomalies. In velocity cross sections the Tacoma and Seattle basins are asymmetric: the basin floor dips gently toward a steep boundary with the adjacent high-velocity uplift, locally with a velocity "overhang" that suggests a basin vergent thrust fault boundary. Crustal fault zones grow from minor folds into much larger structures along strike. Inferred structural relief across the Tacoma fault zone increases by several kilometers westward along the fault zone to Lynch Cove, where we interpret it as a zone of south vergent faulting overthrusting Tacoma basin. In contrast, structural relief along the Seattle fault zone decreases west of Seattle, which we interpret as evidence that the N-S directed compression is being accommodated by slip transfer between the Seattle and Tacoma fault zones. Together, the Tacoma and Seattle fault zones raise the Seattle uplift, one of a series of east-west trending, pop-up structures underlying Puget Lowland from the Black Hills to the San Juan Islands.

Washington

Strain accumulation near Yucca Mountain, Nevada, 1993-1998

A 50‐km aperture geodetic network centered on the proposed high‐level radioactive waste disposal site at Yucca Mountain, Nevada, was surveyed with GPS in 1993 and 1998. The average deformation rate across the area is described by the principal strain rates 22.8±8.8 nstrain yr N77.6°W±13.5° and −8.8±11.9 nstrain yr N12.5°E±13.5° (extension reckoned positive) and a clockwise rotation rate about a vertical axis of 9.6±7.4 nrad yr relative to fixed North America. Quoted uncertainties are standard deviations. Those strain rates are consistent with the geodetic strain rates (2±12 nstrain yr N87°±12°W and −22±12 nstrain yr N03°±12°E) previously reported by . [1999] for the 1983–1998 interval and with the low extension rate (5–20 nstrain yr) [., 1998] inferred from the geologic record. None of those strain rates is consistent with the 50±9 nstrain yr N65°W extension rate for the area reported by Wernicke et al.

Nevada

Response of the San Andreas fault to the 1983 Coalinga-Nuñez earthquakes: an application of interaction-based probabilities for Parkfield

The Parkfield-Cholame section of the San Andreas fault, site of an unfulfilled earthquake forecast in 1985, is the best monitored section of the world's most closely watched fault. In 1983, the M = 6.5 Coalinga and M = 6.0 Nu&ntilde;ez events struck 25 km northeast of Parkfield. Seismicity rates climbed for 18 months along the creeping section of the San Andreas north of Parkfield and dropped for 6 years along the locked section to the south. Right-lateral creep also slowed or reversed from Parkfield south. Here we calculate that the Coalinga sequence increased the shear and Coulomb stress on the creeping section, causing the rate of small shocks to rise until the added stress was shed by additional slip. However, the 1983 events decreased the shear and Coulomb stress on the Parkfield segment, causing surface creep and seismicity rates to drop. We use these observations to cast the likelihood of a Parkfield earthquake into an interaction-based probability, which includes both the renewal of stress following the 1966 Parkfield earthquake and the stress transfer from the 1983 Coalinga events. We calculate that the 1983 shocks dropped the 10-year probability of a M &sim; 6 Parkfield earthquake by 22% (from 54 &plusmn; 22% to 42 &plusmn; 23%) and that the probability did not recover until about 1991, when seismicity and creep resumed. Our analysis may thus explain why the Parkfield earthquake did not strike in the 1980s, but not why it was absent in the 1990s. We calculate a 58 &plusmn; 17% probability of a M &sim; 6 Parkfield earthquake during 2001&ndash;2011.

Journal of Geophysical Research B: Solid Earth

Aftershock triggering by complete Coulomb stress changes

We examine the correlation between seismicity rate change following the 1992, M 7.3, Landers, California, earthquake and characteristics of the complete Coulomb failure stress (CFS) changes (ΔCFS( t )) that this earthquake generated. At close distances the time-varying “dynamic” portion of the stress change depends on how the rupture develops temporally and spatially and arises from radiated seismic waves and from permanent coseismic fault displacement. The permanent “static” portion (ΔCFS) depends only on the final coseismic displacement. ΔCFS diminishes much more rapidly with distance than the transient, dynamic stress changes. A common interpretation of the strong correlation between ΔCFS and aftershocks is that load changes can advance or delay failure. Stress changes may also promote failure by physically altering properties of the fault or its environs. Because it is transient, ΔCFS( t ) can alter the failure rate only by the latter means. We calculate both ΔCFS and the maximum positive value of ΔCFS( t ) (peak ΔCFS( t )) using a reflectivity program. Input parameters are constrained by modeling Landers displacement seismograms. We quantify the correlation between maps of seismicity rate changes and maps of modeled ΔCFS and peak ΔCFS( t ) and find agreement for both models. However, rupture directivity, which does not affect ΔCFS, creates larger peak ΔCFS( t ) values northwest of the main shock. This asymmetry is also observed in seismicity rate changes but not in ΔCFS. This result implies that dynamic stress changes are as effective as static stress changes in triggering aftershocks and may trigger earthquakes long after the waves have passed.

Journal of Geophysical Research B: Solid Earth

Inferring rate and state friction parameters from a rupture model of the 1995 Hyogo-ken Nanbu (Kobe) Japan earthquake

We consider the applicability of laboratory-derived rate-and state-variable friction laws to the dynamic rupture of the 1995 Kobe earthquake. We analyze the shear stress and slip evolution of Ide and Takeo's [1997] dislocation model, fitting the inferred stress change time histories by calculating the dynamic load and the instantaneous friction at a series of points within the rupture area. For points exhibiting a fast-weakening behavior, the Dieterich-Ruina friction law, with values of d c = 0.01–0.05 m for critical slip, fits the stress change time series well. This range of d c is 10–20 times smaller than the slip distance over which the stress is released, D c , which previous studies have equated with the slip-weakening distance. The limited resolution and low-pass character of the strong motion inversion degrades the resolution of the frictional parameters and suggests that the actual d c is less than this value. Stress time series at points characterized by a slow-weakening behavior are well fitted by the Dieterich-Ruina friction law with values of d c ≥0.01–0.05 m. The apparent fracture energy G c can be estimated from waveform inversions more stably than the other friction parameters. We obtain a G c ≈ 1.5×l0 6 J m −2 for the 1995 Kobe earthquake, in agreement with estimates for previous earthquakes. From this estimate and a plausible upper bound for the local rock strength we infer a lower bound for D c of about 0.008 m.

Kobe, Osaka

Crustal structure of central Lake Baikal: Insights into intracontinental rifting

The Cenozoic rift system of Baikal, located in the interior of the largest continental mass on Earth, is thought to represent a potential analog of the early stage of breakup of supercontinents. We present a detailed P wave velocity structure of the crust and sediments beneath the Central Basin, the deepest basin in the Baikal rift system. The structure is characterized by a Moho depth of 39–42.5 km; an 8-km-thick, laterally continuous high-velocity (7.05–7.4 km/s) lower crust, normal upper mantle velocity (8 km/s), a sedimentary section reaching maximum depths of 9 km, and a gradual increase of sediment velocity with depth. We interpret the high-velocity lower crust to be part of the Siberian Platform that was not thinned or altered significantly during rifting. In comparison to published results from the Siberian Platform, Moho under the basin is elevated by <3 km. On the basis of these results we propose that the basin was formed by upper crustal extension, possibly reactivating structures in an ancient fold-and-thrust belt. The extent and location of upper mantle extension are not revealed by our data, and it may be offset from the rift. We believe that the Baikal rift structure is similar in many respects to the Mesozoic Atlantic rift system, the precursor to the formation of the North Atlantic Ocean. We also propose that the Central Baikal rift evolved by episodic fault propagation and basin enlargement, rather than by two-stage rift evolution as is commonly assumed.

Lake Baikal

Spatial extent of a hydrothermal system at Kilauea Volcano, Hawaii, determined from array analyses of shallow long-period seismicity 1. Method

We present a probabilistic method to locate the source of seismic events using seismic antennas. The method is based on a comparison of the event azimuths and slownesses derived from frequency-slowness analyses of array data, with a slowness vector model. Several slowness vector models are considered including both homogeneous and horizontally layered half-spaces and also a more complex medium representing the actual topography and three-dimensional velocity structure of the region under study. In this latter model the slowness vector is obtained from frequency-slowness analyses of synthetic signals. These signals are generated using the finite difference method and include the effects of topography and velocity structure to reproduce as closely as possible the behavior of the observed wave fields. A comparison of these results with those obtained with a homogeneous half-space demonstrates the importance of structural and topographic effects, which, if ignored, lead to a bias in the source location. We use synthetic seismograms to test the accuracy and stability of the method and to investigate the effect of our choice of probability distributions. We conclude that this location method can provide the source position of shallow events within a complex volcanic structure such as Kilauea Volcano with an error of ±200 m.

Hawai'i

Vertical tectonics in northern Escanaba Trough as recorded by thick late Quaternary turbidites

Escanaba Trough, the southernmost segment of the Gorda Ridge, is filled by as much as 500 m of late Quaternary turbidite and hemipelagic sediment. Coring at Deep Sea Drilling Project Site 35 and Ocean Drilling Program (ODP) Sites 1037 and 1038 together with 4.5-kHz deep-tow and 3.5-kHz surface-ship seismic reflection profiles enable a distinct pattern of reflections to be mapped throughout Escanaba Trough in the upper part of this sediment fill. The uppermost 80 m of turbidite sediment, which includes at least 11 turbidity current events, were deposited in <7 kyr. Nine of these turbidites are found throughout Escanaba Trough at water depths >3200 m. The turbidity currents were trapped upon entering Escanaba Trough, resulting in all of the sediment in suspension in the flows being deposited. The thickness of the turbidite layers reflects both the flow thickness and the vertical grain concentration within the flow that deposited the layer. Variations in the turbidite thickness with respect to water depth can be used to estimate the degree of relative vertical movement within the floor of Escanaba Trough. In the area of hydrothermal activity near ODP Site 1038, uplift of as much as 140 m has occurred over the past 8 kyr.

California

Viscoelastic shear zone model of a strike-slip earthquake cycle

I examine the behavior of a two-dimensional (2-D) strike-slip fault system embedded in a 1-D elastic layer (schizosphere) overlying a uniform viscoelastic half-space (plastosphere) and within the boundaries of a finite width shear zone. The viscoelastic coupling model of Savage and Prescott [1978] considers the viscoelastic response of this system, in the absence of the shear zone boundaries, to an earthquake occurring within the upper elastic layer, steady slip beneath a prescribed depth, and the superposition of the responses of multiple earthquakes with characteristic slip occurring at regular intervals. So formulated, the viscoelastic coupling model predicts that sufficiently long after initiation of the system, (1) average fault-parallel velocity at any point is the average slip rate of that side of the fault and (2) far-field velocities equal the same constant rate. Because of the sensitivity to the mechanical properties of the schizosphere-plastosphere system (i.e., elastic layer thickness, plastosphere viscosity), this model has been used to infer such properties from measurements of interseismic velocity. Such inferences exploit the predicted behavior at a known time within the earthquake cycle. By modifying the viscoelastic coupling model to satisfy the additional constraint that the absolute velocity at prescribed shear zone boundaries is constant, I find that even though the time-averaged behavior remains the same, the spatiotemporal pattern of surface deformation (particularly its temporal variation within an earthquake cycle) is markedly different from that predicted by the conventional viscoelastic coupling model. These differences are magnified as plastosphere viscosity is reduced or as the recurrence interval of periodic earthquakes is lengthened. Application to the interseismic velocity field along the Mojave section of the San Andreas fault suggests that the region behaves mechanically like a ???600-km-wide shear zone accommodating 50 mm/yr fault-parallel motion distributed between the San Andreas fault system and Eastern California Shear Zone. Copyright 2001 by the American Geophysical Union.

California

Lower crustal deformation beneath the central Transverse Ranges, southern California: Results from the Los Angeles Region Seismic Experiment

We present a P wave velocity model derived from active source seismic data collected during the 1994 Los Angeles Region Seismic Experiment. Our model extends previously published upper crustal velocity models to mantle depths. Our model was developed by both ray tracing through a layered model and calculating travel times through a gridded model. It includes an 8-km-thick crustal root centered beneath the surface trace of the San Andreas fault, north of the highest topography in the San Gabriel Mountains. A simple mass balance calculation suggests that ∼36 km of north-south shortening across the San Andreas fault in the central Transverse Ranges could have formed this root. If north-south compression began when the “Big Bend” in the San Andreas fault formed at ∼5 Ma, 36 km of shortening implies a north-south contraction rate of ∼7.1 mm/yr across the central Transverse Ranges. If, instead, north-south compression began when the Transverse Ranges formed at 3.4–3.9 Ma, 36 km of shortening implies a contraction rate of 9.2–10.6 mm/yr. North of the San Andreas fault, the Mojave Desert crust has a low-velocity (6.3 km/s) mid and lower crust and a 28-km-deep Moho. South of the San Andreas fault, beneath the Los Angeles and San Gabriel Valley basins, there is a fast (6.6–6.8 km/s), thick (10–12 km) lower crust with a 27-km-deep Moho. Farther south still, the lower crust of the Continental Borderland is fast (6.6–6.8 km/s) and thin (5 km) with a shallow (22 km deep) Moho.

California