Geochemistry of Cenezoic volcanic rocks, Baja California, Mexico: Implications for the petrogenesis of post-subduction magmas

Late Cenozoic volcanism in Baja California records the effects of cessation of subduction at a previously convergent, plate margin. Prior to 12.5 m.y., when subduction along the margin of Baja ceased, the predominant volcanic activity had a calc-alkaline signature, ranging in composition from basalt to rhyolite. Acidic pyroclastic activity was common, and possibly represented the westermost, distal edge of the Sierra Madre Occidental province. After 12.5 m.y., however, the style and composition of the magmatic products changed dramatically. The dominant rock type within the Jaraguay and San Borja volcanic fields is a magnesian andesite, with up to 8% MgO at 57% SiO₂, low Fe/Mg ratios, and high Na/K ratios. These rocks have unusual trace-element characteristics, with high abundances of Sr (up to 3000 ppm), low contents of Rb; K/Rb ratios are very high (usually over 1000, and up to 2500), and Rb/Sr ratios are low (less than 0.01). Furthermore, La_n/Yb_n ratios are high, consistent with derivation from a mantle source with fractionated REE patterns. ⁸⁷Sr/⁸⁶Sr ratios are less than 0.7048, and usually less than 0.7040, whereas the pre-12.5 m.y. lavas have ⁸⁷Sr/⁸⁶Sr ratios between 0.7038 and 0.7063. We have previously termed these rocks bajaites, in order to distinguish them from other magnesian andesites. Bajaites also occur in southernmost Chile and the Aleutian Islands, areas which also have histories of attempted or successful ridge subduction.It is proposed that the bajaite series is produced during the unusual physico-chemical conditions operating during the subduction of young oceanic lithosphere, or subduction of a spreading centre. During normal subduction, the oceanic crust dehydrates, releasing volatiles (water, Rb and other large-ion lithophile elements) into the overlying wedge. Subduction of younger crust will result in a progressive decrease, and eventual cessation of the transfer of volatiles when subduction stops. Thermal rebound of the mantle may cause the slab to melt, perhaps under eclogitestable conditions. The resulting melt will be heavy-REE-depleted, perhaps dacitic, but will otherwise inherit MORB-like Rb/Sr and K/Rb ratios. The ascending melt will react with the mantle to form the source of the bajaitic rocks. Furthermore, any amphibole in the mantle, stabilised during the higher PH₂O conditions of earlier subduction, will break down and contribute a high-K/Rb ratio component.The implications of this study are that firstly, the subducted slab does not contribute a highly fractionated REE component in most modern arcs (i.e. the slab does not melt); secondly, Rb has a very short residence time in the mantle, and its abundance in arc rocks is a direct reflection of the input from the dehydrating slab; and thirdly, bajaitelike rocks may provide recognition of attempted or successful ridge subduction in the geologic past.     

Dunes of the Gran Desierto Sand-Sea,Sonora, Mexico

The Gran Desierto Sand-Sea contains dunes of crescentic and star form in simple, compound, and complex varieties. The dunes have developed in bimodal to complex wind regimes of intermediate energy. Transitions from simple through compound to complex crescentic dunes are associated with regional changes in wind regimes. Growth of large star dunes takes place by merging of smaller crescentic and reversing dunes from southerly directions and reworking of sand by northerly and westerly winds. Although wind regimes appear to be the major control of dune morphology in this sand-sea, there is a close correlation between the spacing of simple crescentic dunes and the grain size of the coarse 20th percentile.     

Shallow Miocene basaltic magma reservoirs in the Bahia de Los Angeles basin, Baja California, Mexico

The basement in the Bahía de Los Angeles basin consists of Paleozoic metamorphic rocks and Cretaceous granitoids. The Neogene stratigraphy overlying the basement is formed, from the base to the top, by andesitic lava flows and plugs, sandstone and conglomeratic horizons, and Miocene pyroclastic flow units and basaltic flows. Basaltic dikes also intrude the whole section. To further define its structure, a detailed gravimetric survey was conducted across the basin about 1 km north of the Sierra Las Flores. In spite of the rough and lineal topography along the foothills of the Sierra La Libertad, we found no evidence for large-scale faulting. Gravity data indicates that the basin has a maximum depth of 120 m in the Valle Las Tinajas and averages 75 m along the gravimetric profile. High density bodies below the northern part of the Sierra Las Flores and Valle Las Tinajas are interpreted to be part of basaltic dikes. The intrusive body located north of the Sierra Las Flores is 2.5 km wide and its top is about 500 m deep. The lava flows of the top of the Sierra Las Flores, together with the distribution of basaltic activity north of this sierra, suggests that this intrusive body continues for 20 km along a NNW-trending strike. Between the sierras Las Flores and Las Animas, a 0.5-km-wide, 300-m-thick intrusive body is interpreted at a depth of about 100 m. This dike could be part of the basaltic activity of the Cerro Las Tinajas and the small mounds along the foothills of western Sierra Las Animas. The observed local normal faulting in the basin is inferred to be mostly associated with the emplacement of the shallow magma reservoirs below Las Flores and Las Tinajas.     

Anomalous Seismicity and Accelerating Moment Release Preceding the 2001 and 2002 Earthquakes in Northern Baja California, Mexico

— An algorithm recently developed by RUNDLE et al. (2002) to find regions of anomalous seismic activity associated with large earthquakes identified the location of an M _w = 5.6 earthquake near Calexico, Mexico. In this paper we analyze the regional seismicity before this event, and a nearby M _w = 5.7 event, using time-to-failure algorithms developed by BOWMAN et al. (1998) and BOWMAN and KING (2001a,b). The former finds the radius of a circular region surrounding the epicenter that optimizes the time-to-failure acceleration of seismic release. The latter optimizes acceleration based on the expected stress accumulation pattern for a dislocation source. Both methods found a period of accelerating seismicity in an optimal region, the size of which agrees with previously proposed scaling relations. This positive result suggests that the Rundle algorithm may provide a useful technique to identify regions of accelerating seismicity, which can then be analyzed using signal optimization time-to-failure techniques.     

The 2010 M w 7.2 El Mayor-Cucapah Earthquake Sequence, Baja California, Mexico and Southernmost California, USA: Active Seismotectonics along the Mexican Pacific Margin

The El Mayor-Cucapah earthquake sequence started with a few foreshocks in March 2010, and a second sequence of 15 foreshocks of M?>?2 (up to M4.4) that occurred during the 24?h preceding the mainshock. The foreshocks occurred along a north?Csouth trend near the mainshock epicenter. The M _w 7.2 mainshock on April 4 exhibited complex faulting, possibly starting with a ~M6 normal faulting event, followed ~15?s later by the main event, which included simultaneous normal and right-lateral strike-slip faulting. The aftershock zone extends for 120?km from the south end of the Elsinore fault zone north of the US?CMexico border almost to the northern tip of the Gulf of California. The waveform-relocated aftershocks form two abutting clusters, each about 50?km long, as well as a 10?km north?Csouth aftershock zone just north of the epicenter of the mainshock. Even though the Baja California data are included, the magnitude of completeness and the hypocentral errors increase gradually with distance south of the international border. The spatial distribution of large aftershocks is asymmetric with five M5+ aftershocks located to the south of the mainshock, and only one M5.7 aftershock, but numerous smaller aftershocks to the north. Further, the northwest aftershock cluster exhibits complex faulting on both northwest and northeast planes. Thus, the aftershocks also express a complex pattern of stress release along strike. The overall rate of decay of the aftershocks is similar to the rate of decay of a generic California aftershock sequence. In addition, some triggered seismicity was recorded along the Elsinore and San Jacinto faults to the north, but significant northward migration of aftershocks has not occurred. The synthesis of the El Mayor-Cucapah sequence reveals transtensional regional tectonics, including the westward growth of the Mexicali Valley and the transfer of Pacific?CNorth America plate motion from the Gulf of California in the south into the southernmost San Andreas fault system to the north. We propose that the location of the 2010 El Mayor-Cucapah, as well as the 1992 Landers and 1999 Hector Mine earthquakes, may have been controlled by the bends in the plate boundary.     

Reconnaissance geophysics and geology of the Pinacate craters,Sonora, Mexico

The Pinacate volcanic field of northwest Sonora, Mexico contains eleven tuff-rimmed explosion craters as well as hundreds of cinder cones. A multiphase eruptive and collapse history for each of the craters is evidenced by: cinder cones exposed in cross-section in crater rims, rim-breaking cinder cones, major eccentric vents, erosion hiatuses between tuff deposition, and coalescing craters exhibiting differeng degrees of erosion. Basaltic intrusions occurred along ring fractures and in the center of at least one crater. Magnetic surveys indicate that explosion craters have tuffaceous feeder vents whereas cinder cones have basaltic, highly magnetic plugs. Computer models illustrate that magnetic terrain effects can contribute to anomalies observed associated with craters, and need to be considered where magnetic surveys include irregularities such as rims or peaks.     

Structural Geometry of a Sector of the Colorado River Delta,Baja California,Mexico, Based on Seismic Reflections

A structural study in the SW section of the Colorado River delta using seismic reflection data is presented. The study area is located along the Cerro Prieto transform fault, which extends from the northern Gulf of California through the Mexicali Valley and is an active fault within the Pacific-North American plate boundary zone. The research was supported by a database of five seismic profiles with a total length of 215 km, collected in the early 80’s by Petróleos Mexicanos. The results show a high density of faults, most of which are buried by sediments. Within the Cerro Prieto fault zone, several faults were identified, such as: Palmas, Mesa, and Pangas Viejas, until now unknown. In addition, even though the Indiviso fault was investigated and superficially identify prior to this work, herein mapped at depth. West of the Cerro Prieto fault zone lies the Las Tinajas basin, bound by the Dunas and Saldaña faults and by the Montague basin to the southeast. The deformation zone along the plate boundary is 18-km-wide, stretching from the Cerro Prieto fault in the east to the Pangas Viejas fault in the west. The orientations of the faults are NW–SE, and if projected from the southern side of the Sierra Cucapah southward, the faults tend to join the Cerro Prieto fault. In the Las Tinajas basin, the acoustic basement is deeper than 5,000 m. Some of the largest vertical displacements generated by the 2010 7.2-Mw El Mayor-Cucapah earthquake occurred southeast of the epicenter and coincided with the location of the Pangas Viejas Fault, which is buried by sediments. Before this event, seismic activity was very low, and no structures were known in the area. In this paper, we demonstrate that there are at least seven major faults that may now pose a high seismic hazard.     

P- and S-Wave Site Response of the Seismic Network RESNOM Determined from Earthquakes of Northern Baja California, Mexico

—We determined the response to P- and S-wave incidence of the permanent stations of the seismic network of Baja California (RESNOM) using two independent methods. We selected 65 events with magnitudes between 2.2 and 4.8 and hypocentral distances ranging between 5 and 330 km. The site response of the ten stations analyzed was first estimated using average spectral ratios between the horizontal and the vertical components of motion (H/V ratios). As a second approach we performed a simultaneous inversion for source and site. In order to invert the spectral records to determine the site response, we made an independent estimate of the attenuation for two different source-station path regions. Then we corrected the spectral records for the attenuation effect before we made the inversion. Although the average H/V ratio of many sites is inside the error bars of the site response estimated with the spectral inversion, the spectral inversion tends to give higher values. For the S wave some sites show similar frequency of predominant peak when comparing the responses obtained with both methods. In contrast, for the P waves the H/V ratios disagree with the results of the inversion. In general, the site response of the stations is strongly frequency dependent for both P and S waves. We also found that the natural frequency of resonance of the sites is near 0.5 Hz for P and near 0.8 Hz for the S waves.     

Flow direction and magnetic mineralogy of lava flows from the central parts of the Peninsula of Baja California,Mexico

Magnetic properties of two apparently unrelated lava suites (one tholeiitic and the other alkalic) coexisting in the central parts of the Peninsula of Baja California, Mexico, were measured in this study. Macroscopic indicators and measurements of anisotropy of magnetic susceptibility (AMS) were combined to infer the flow direction on those lavas. These measurements were used to investigate the likelihood of the existence of an abnormally long tholeiitic lava flow in this part of the Peninsula. The obtained results indicate that the tholeiitic flows in the region constitute regular length flows, and are likely to have been issued through vents located within the area of study. Additionally, comparison of the magnetic signature of tholeiitic and alkalic lavas reveals a systematic difference in some of the magnetic parameters (bulk susceptibility and Curie temperatures). In particular, it is proposed that the narrow range of Curie temperatures characteristic of each lava type can be used as a reliable proxy for the identification of each lava type in the region. All of these findings show that the alkalic and tholeiitic lavas of the region are more closely related to each other than previously suspected.Editorial responsibility: R. Cioni     

Spatial Distribution,Scaling and Self-similar Behavior of Fracture Arrays in the Los Planes Fault,Baja California Sur,Mexico

We present a case of detailed analysis of fracture arrays spanning four orders of magnitude in length; all of them measured at a single natural site by acquiring images at progressively larger scales. There is a high dispersion of cumulative-length exponents, box dimensions and fracture densities. However, the fractal analysis supports the fractal nature of fracture arrays. Our data indicate the existence of an upper limit for the density parameters, as similarly reported by other authors. We prove that box dimension is in inverse relation with fracture concentration and in direct relation with fracture density. These relations are also observed in our data and additionally there is an upper limit for the box dimensions. We interpret the dispersion in our results as more fundamental than methodological problems. It represents a truncation in the complete evolution of the fracture systems because in natural cases strain initiates overprinting of previous fracture arrays. Considering that larger fractures accommodate strain more efficiently than small fractures, the generation of small fractures is inhibited in the presence of pre-existing larger fractures. Maximum values of fracture density prevent accommodating an excess of strain in a single or restricted range of scales; we claim this condition produces migration of fracturing to larger scales originating fracture scaling.     

Geology, ore deposits and hydrothermal venting in Bahía Concepción, Baja California Sur, Mexico

Abstract Bahía Concepción is located in the eastern coast of the Baja California peninsula and it is shaped by northwestern–southeastern normal faults. These are associated with a 12–6 Ma rifting episode, although some have been reactivated since the Pliocene. The most abundant rocks correspond to the arc related Comondú Group, Oligocene to Miocene, which forms a mainly calc‐alkaline volcanic and volcaniclastic sequence. There are less extensive outcrops of sedimentary rocks, lava flows, domes and pyroclastic rocks of Pliocene to Quaternary ages. The Neogene volcanism in the area indicates a shift from a subduction regime to an intraplate volcanism related to continental extension and the opening of an oceanic basin. The Bahía Concepción area contains numerous Mn ore deposits, being the biggest at El Gavilán and Guadalupe. The Mn deposits occur as veins, breccias and stockworks, and are composed by Mn oxides (pyrolusite, coronadite, romanechite), dolomite, quartz and barite. The deposits are hosted in volcanic rocks of the Comondú Group and, locally, in Pliocene sedimentary rocks. Thus, the Mn deposits formed between the Middle Miocene and the Pliocene. The mineralized structures are associated with Miocene northwestern–southeastern fault systems, which are analogous to those associated with the Cu‐Co‐Zn‐Mn deposits of El Boleo. The Bahía Concepción area also bears subaerial and submarine hot springs, which are associated with the same fault systems and host rocks. The submarine and subaerial geothermal manifestations south of the bay are possibly related with recent volcanism. The geothermal manifestations within the bay are intertidal hot springs and shallow submarine diffuse venting areas. Around the submarine vents (5–15 m deep, 87°C), Fe‐oxyhydroxide crusts with pyrite and cinnabar precipitate. In the intertidal vents (62°C), aggregates of opal, calcite, barite and Ba‐rich Mn oxides occur covered by silica‐carbonate stromatolitic sinters. Some 10–30 cm thick crustiform veins formed by chalcedony, calcite and barite are also found close to the vents. The hydrothermal fluids exhibit mixed isotopic compositions between δ¹⁸O‐enriched meteoric and local marine water. The precipitation of Ba‐rich Mn oxides around the vent sites could be an active analog for the processes that produced Miocene to Pliocene hydrothermal Mn‐deposits.     

Benthic foraminifera show some resilience to ocean acidification in the northern Gulf of California,Mexico

Extensive CO₂ vents have been discovered in the Wagner Basin, northern Gulf of California, where they create large areas with lowered seawater pH. Such areas are suitable for investigations of long-term biological effects of ocean acidification and effects of CO₂ leakage from subsea carbon capture storage. Here, we show responses of benthic foraminifera to seawater pH gradients at 74–207 m water depth. Living (rose Bengal stained) benthic foraminifera included Nonionella basispinata, Epistominella bradyana and Bulimina marginata. Studies on foraminifera at CO₂ vents in the Mediterranean and off Papua New Guinea have shown dramatic long-term effects of acidified seawater. We found living calcareous benthic foraminifera in low pH conditions in the northern Gulf of California, although there was an impoverished species assemblage and evidence of post-mortem test dissolution.     

A constrained 2D gravity model of the Sebastián Vizcaíno Basin, Baja California Sur, Mexico

The subsurface geometry of the Sebastián Vizcaíno Basin is obtained from the 2D inversion of gravity data, constrained by a density-versus-depth relationship derived from an oil exploration deep hole. The basin accumulated a thick pile of marine sediments that evolved in the fore-arc region of the compressive margin prevalent along western North America during Mesozoic and Tertiary times. Our interpretation indicates that the sedimentary infill in the Sebastián Vizcaíno Basin reaches a maximum thickness of about 4 km at the centre of a relatively symmetric basin. At the location of the Suaro-1 hole, the depth to the basement derived from this work agrees with the drilled interface between calcareous and volcaniclastic members of the Alisitos Formation. A sensitivity analysis strongly suggests that the assumed density function leads to a nearly unique solution of the inverse problem.     

Upper crust attenuation in the Basin and Range Province of Sonora,Mexico

Journal of Seismology - We estimated the S-wave quality factor Q(f) of the upper crust in northeastern Sonora, Mexico, using a two-layer model. The first layer of the model has a thickness of 1 km,...     

Faulting processes of the 1956 San Miguel,Baja California,earthquake sequence

Body waveform modeling is used to determine the source processes of three large earthquakes (magnitude 6.8, 6.4, 6.3) occurring between February 9 and 15, 1956 along the San Miguel fault in northern Baja California, Mexico. Results of the modeling suggest that the mainshock on February 9 was responsible for the 20 km of surface faulting observed during the sequence. Although previous researchers have suggested a complex rupture history for the mainshock, uncertainty estimates of source-time function shape indicate single or double source models fit the observed waveforms equally well. The February 15 aftershock, however, appears to have consisted of two events. Locations and focal mechanisms obtained for the three events suggest that the rupture process may have been controlled by cross faults to the main trace of the San Miguel fault. The good correlation between source parameter information and the surficial geology/geometry of the San Miguel fault zone demonstrates the usefulness of waveform modeling studies in unravelling the complexities of historic multi-event earthquake sequences.