Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Evolution and Development”

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 163 records · Page 9Linked to original sources

Advanced Tools for River Science: EAARL and MD_SWMS: Chapter 3

Disruption of flow regimes and sediment supplies, induced by anthropogenic or climatic factors, can produce dramatic alterations in river form, vegetation patterns, and associated habitat conditions. To improve habitat in these fluvial systems, resource managers may choose from a variety of treatments including flow and/or sediment prescriptions, vegetation management, or engineered approaches. Monitoring protocols developed to assess the morphologic response of these treatments require techniques that can measure topographic changes above and below the water surface efficiently, accurately, and in a standardized, cost-effective manner. Similarly, modeling of flow, sediment transport, habitat, and channel evolution requires characterization of river morphology for model input and verification. Recent developments by the U.S. Geological Survey with regard to both remotely sensed methods (the Experimental Advanced Airborne Research LiDAR; EAARL) and computational modeling software (the Multi-Dimensional Surface-Water Modeling System; MD_SWMS) have produced advanced tools for spatially explicit monitoring and modeling in aquatic environments. In this paper, we present a pilot study conducted along the Platte River, Nebraska, that demonstrates the combined use of these river science tools.

Report↗

Mid-Cretaceous oblique rifting of West Antarctica: Emplacement and rapid cooling of the Fosdick Mountains migmatite-cored gneiss dome

In Marie Byrd Land, West Antarctica, the Fosdick Mountains migmatite-cored gneiss dome was exhumed from mid- to lower middle crustal depths during the incipient stage of the West Antarctic Rift system in the mid-Cretaceous. Prior to and during exhumation, major crustal melting and deformation included transfer and emplacement of voluminous granitic material and numerous intrusions of mantle-derived diorite in dikes. A succession of melt- and magma-related structures formed at temperatures in excess of 665 ± 50 °C based on Ti-in-zircon thermometry. These record a transition from wrench to oblique extensional deformation that culminated in the development of the oblique South Fosdick Detachment zone. Solid-state fabrics within the detachment zone and overprinting brittle structures record translation of the detachment zone and dome to shallow levels. To determine the duration of exhumation and cooling, we sampled granite and gneisses at high spatial resolution for U–Pb zircon geochronology and 40 Ar/ 39 Ar hornblende and biotite thermochronology. U–Pb zircon crystallization ages for the youngest granites are 102 Ma. Three hornblende ages are 103 to 100 Ma and 12 biotite ages are 101 to 99 Ma. All overlap within uncertainty. The coincidence of zircon crystallization ages with 40 Ar/ 39 Ar cooling ages indicates cooling rates > 100 °C/m.y. that, when considered together with overprinting structures, indicates rapid exhumation of granite and migmatite from deep to shallow crustal levels within a transcurrent setting. Orientations of structures and age-constrained crosscutting relationships indicate counterclockwise rotation of stretching axes from oblique extension into nearly orthogonal extension with respect to the Marie Byrd Land margin. The rotation may be a result of localized extension arising from unroofing and arching of the Fosdick dome, extensional opening within a pull-apart zone, or changes in plate boundary configuration. The rapid tectonic and temperature evolution of the Fosdick Mountains dome lends support to recently developed numerical models of crustal flow and cooling in orogenic crust undergoing extension/transtension, and accords with numerous studies of migmatite-cored gneiss domes in transcurrent settings.

LITHOS↗

A holistic modelling approach to project the evolution of inlet-interrupted coastlines over the 21st century

Approximately one quarter of the World’s sandy beaches, most of which are interrupted by tidal inlets, are eroding. Understanding the long-term (50-100 year) evolution of inlet-interrupted coasts in a changing climate is therefore of great importance for coastal zone planners and managers. This study therefore focuses on the development and piloting of an innovative model that can simulate the climate-change driven evolution of inlet-interrupted coasts at 50-100 year time scales, while taking into account the contributions from catchment-estuary-coastal systems in a holistic manner. In this new model, the evolution of inlet-interrupted coasts is determined by: (1) computing the variation of total sediment volume exchange between the inlet-estuary system and its adjacent coast, and (2) distributing the computed sediment volume along the inlet-interrupted coast as a spatially and temporally varying quantity. The exchange volume, as computed here, consists of three major components: variation in fluvial sediment supply; basin (or estuarine) infilling due to the sea-level rise-induced increase in accommodation space; and estuarine sediment volume change due to variations in river discharge. To pilot the model, it is here applied to three different catchment-estuary-coastal systems: the Alsea estuary (Oregon, USA), Dyfi estuary (Wales, UK), and Kalutara inlet (Sri Lanka). Results indicate that all three systems will experience sediment deficits by 2100 (i.e. sediment importing estuaries). However, processes and system characteristics governing the total sediment exchange volume, and thus coastline change, vary markedly among the systems due to differences in geomorphic settings and projected climatic conditions. These results underline the importance of accounting for the different governing processes when assessing the future evolution of inlet-interrupted coastlines.

Frontiers in Marine Science↗

The international river interface cooperative: Public domain flow and morphodynamics software for education and applications

This paper describes a new, public-domain interface for modeling flow, sediment transport and morphodynamics in rivers and other geophysical flows. The interface is named after the International River Interface Cooperative (iRIC), the group that constructed the interface and many of the current solvers included in iRIC. The interface is entirely free to any user and currently houses thirteen models ranging from simple one-dimensional models through three-dimensional large-eddy simulation models. Solvers are only loosely coupled to the interface so it is straightforward to modify existing solvers or to introduce other solvers into the system. Six of the most widely-used solvers are described in detail including example calculations to serve as an aid for users choosing what approach might be most appropriate for their own applications. The example calculations range from practical computations of bed evolution in natural rivers to highly detailed predictions of the development of small-scale bedforms on an initially flat bed. The remaining solvers are also briefly described. Although the focus of most solvers is coupled flow and morphodynamics, several of the solvers are also specifically aimed at providing flood inundation predictions over large spatial domains. Potential users can download the application, solvers, manuals, and educational materials including detailed tutorials at www.-i-ric.org . The iRIC development group encourages scientists and engineers to use the tool and to consider adding their own methods to the iRIC suite of tools.

Advances in Water Resources↗

Modeling barrier island response to sea-level rise in the Outer Banks, North Carolina

An 8500-year Holocene simulation developed in GEOMBEST provides a possible scenario to explain the evolution of barrier coast between Rodanthe and Cape Hatteras, NC. Sensitivity analyses suggest that in the Outer Banks, the rate of sea-level rise is the most important factor in determining how barrier islands evolve. The Holocene simulation provides a basis for future simulations, which suggest that if sea level rises up to 0.88 m by AD 2100, as predicted by the highest estimates of the Intergovernmental Panel on Climate Change, the barrier in the study area may migrate on the order of 2.5 times more rapidly than at present. If sea level rises beyond IPCC predictions to reach 1.4–1.9 m above modern sea level by AD 2100, model results suggest that barrier islands in the Outer Banks may become vulnerable to threshold collapse, disintegrating during storm events, by the end of the next century. Consistent with sensitivity analyses, additional simulations indicate that anthropogenic activities, such as increasing the rate of sediment supply through beach nourishment, will only slightly affect barrier island migration rates and barrier island vulnerability to collapse.

North Carolina↗

Social effects of rabies infection in male vampire bats (Desmodus rotundus)

Rabies virus (RABV) transmitted by the common vampire bat ( Desmodus rotundus ) poses a threat to agricultural development and public health throughout the Neotropics. The ecology and evolution of rabies host-pathogen dynamics are influenced by two infection-induced behavioral changes. RABV-infected hosts often exhibit increased aggression which facilitates transmission, and rabies also leads to reduced activity and paralysis prior to death. Although several studies document rabies-induced behavioral changes in rodents and other dead-end hosts, surprisingly few studies have measured these changes in vampire bats, the key natural reservoir throughout Latin America. Here, we take advantage of an experiment designed to test the safety and efficacy of an oral rabies vaccine in captive male vampire bats to quantify for the first time how rabies affects allogrooming and aggressive behaviors in the vampire bat. Compared to non-rabid vampire bats, rabid individuals reduced their allogrooming prior to death, but we did not detect increases in aggression among bats. To put our results in context, we review what is known and what remains unclear about behavioral changes of rabid vampire bats.

Biology Letters↗

Developing numeric chlorophyll a criteria through science, enhanced monitoring, and a collaborative stakeholder partnership

Chlorophyll a water quality criteria provide the means for regulating nutrient pollution in surface water systems without many of the technical challenges associated with nitrogen and phosphorus thresholds. However, the development of defensible chlorophyll a criteria is not without its own challenges, particularly in complex coastal ecosystems. We describe the evolution of water column chlorophyll a criteria specific to the Chesapeake Bay estuary, from a general narrative to numeric thresholds developed to meet multiple regulatory management goals. We then present a unique approach to numeric chlorophyll a criteria derivation, one that integrates continuous and underway water quality monitoring datasets. These datasets were used for the dual purposes of characterizing baseline conditions and determining chlorophyll a concentrations that confer unacceptable levels of risk to a specific water body in both space and time. This novel approach was used by the Commonwealth of Virginia to update the numeric chlorophyll a criteria it had originally established for the tidal James River, a subestuary of the Chesapeake Bay, in 2005. Tidal James River chlorophyll a concentrations were evaluated in terms of their relationships with harmful algal blooms, dissolved oxygen, pH, and water clarity. The importance of spatially and temporally intensive monitoring datasets to chlorophyll a criteria development is highlighted in this case study.

Virginia↗

Mid-crustal flow during Tertiary extension in the Ruby Mountains core complex, Nevada

Structural analysis and geochronologic data indicate a nearly orthogonal, late Eocene–Oligocene flow pattern in migmatitic infrastructure immediately beneath the kilometer-thick, extensional, mylonitic shear zone of the Ruby Mountains metamorphic core complex, Nevada. New U-Pb radiometric dating indicates that the development of a northward-trending lineation in the infrastructure is partly coeval with the development of a pervasive, west-northwest–trending lineation in the mylonitic shear zone. U-Pb monazite data from the leucogranite orthogneiss of Thorpe Creek indicate a crystallization age of ca. 36–39 Ma. Zircon fractions from a biotite monzogranite dike yield an age of ca. 29 Ma. The three dated samples from these units exhibit a penetrative, approximately north-south–trending elongation lineation. This lineation is commonly defined by oriented bundles of sillimanite and/or elongated aggregates of quartz and feldspar, indicating a synmetamorphic and syndeformational origin. The elongation lineation can be interpreted as a slip line in the flow plane of the migmatitic, nonmylonitic infrastructural core of the northern Ruby Mountains. A portion of this midcrustal flow is coeval with the well-documented, west-northwest sense of slip in the structurally overlying kilometer-thick, mid-Tertiary mylonitic shear zone. Lineations in the mylonitic zone are orthogonal to those in the deeper infrastructure, suggesting fundamental plastic decoupling between structural levels in this core complex. Furthermore, the infrastructure is characterized by overlapping, oppositely verging fold nappes, which are rooted to the east and west. One of the nappes may be synkinematic with the intrusion of the late Eocene orthogneiss of Thorpe Creek. In addition, the penetrative, elongation lineation in the infrastructure is subparallel to hinge lines of parasitic folds developed synchronous with the fold nappes, suggesting a kinematically related evolution.

Nevada, Utah↗

Sediment transport, terrace incision, and sediment chemistry in subbasins of differing land use, Stinking Water Creek Basin near Rangely, Colorado

In many areas of western Colorado, Late Cretaceous Mancos Shale is a source of sediment, salinity, and selenium to surface water. Natural weathering processes and human activity can expose and move sediments from the Mancos Shale, leading to increased salinity and selenium concentrations in surface water. The U.S. Geological Survey, in cooperation with the Bureau of Land Management, studied four subbasins of Stinking Water Creek near Rangely, Colorado, to provide information regarding landscapes on Mancos Shale that have undergone land-use change. Rangeland grazing is present in all four subbasins. Subbasins A1 and A2 have land use related to energy development and subbasins B1 and B2 do not. Aerial photographs, sediment chemistry, optically stimulated luminescence (OSL) ages, cross-section surveys, and Watershed Erosion Prediction Project (WEPP) modeling were used to characterize the history of land-use change by estimating stream channel incision rates and sediment, salinity, and selenium loads. All subbasins had upper and lower terraces and a modern streambed channel, and these surfaces were dated using OSL. In subbasins A1, A2, and B1, the streambed age was older than the terraces, indicating the channels are incising and widening into the older alluvium. In subbasin B2, however, the streambed age was much younger than the terraces, indicating aggradation is the main process. Differences in channel morphology were not always related to land use but to the subbasin position (upstream compared to downstream). Channels in subbasins A1 and A2 were generally wider and shallower than channels in subbasins B1 and B2. Landscape disturbances related to energy development in subbasins A1 and A2 could have caused channel evolution to progress from incision to widening processes.

Colorado↗

Repeating earthquakes during multiple phases of unrest and eruption at Mount Agung, Bali, Indonesia, 2017

In 2017, Mount Agung produced a small (VEI 2) eruption that was preceded by an energetic volcano-tectonic (VT) swarm (>800 earthquakes per day up to M4.9) and two months of declining activity. The period of decreased seismic activity complicated forecasting efforts for scientists monitoring the volcano. We examine the time history of earthquake families at Mount Agung in search of additional insight into the temporal changes in the shallow crust prior to eruption. Specifically, we analyze the period of declining seismic activity about five weeks prior to the eruption when forecasting uncertainty was greatest. We use REDPy (Hotovec-Ellis and Jeffries, 2016) to build a catalog of 6,508 earthquakes from 18 October 2017–15 February 2018 and group them into families of repeating earthquakes based on waveform similarity using a cross-correlation coefficient threshold of 0.8. We show that the evolution of earthquake families provides evidence that Mount Agung was progressing toward eruption even though overall earthquake rates and seismic-energy-release declined. We find that earthquake families that dominated seismicity during the beginning of the crisis ceased near the onset of tremor on 12 November 2017. Then, earthquake families took on characteristics commonly observed during effusive phases of eruptions on 15 November—a full six days before the first phreatomagmatic eruption on 21 November 2017 and a full ten days before the actual onset of lava effusion on 25 November 2017. We interpret the transitions in seismicity as the manifestation of a three-phase physical model including an Intrusion Phase, a Transition Phase, and a Eruptive Phase. During the Intrusion Phase, seismicity was dominated by VT earthquakes with a relatively high percentage of repeaters (59%) grouped into numerous (65) simultaneous families. During the Eruptive Phase, seismicity included both VT and low frequency earthquakes that grouped into relatively long-lived families despite a low overall percentage of repeaters (14%). The Transition Phase exhibited characteristics of earthquake families between the Intrusion Phase and Eruptive Phase. We conclude that the time history of earthquake families provides insight into the evolution of the stress distribution in the volcanic edifice, the development of the volcanic conduit, and seismogenesis of magma effusion. Finally, we discuss the role that repeating earthquakes could play in real-time monitoring at restless volcanoes. Our work suggests eruption forecasts can be improved by incorporating automatic processing codes to assist seismologists during sustained periods of high earthquake rates, even at sparsely monitored volcanoes.

Bali, Mount Agung↗

Tick control: Trapping, bio-control, host management and other alternative strategies

Biology of Ticks is the most comprehensive work on tick biology and tick-borne diseases. This second edition is a multi-authored work, featuring the research and analyses of renowned experts across the globe. Spanning two volumes, the book examines the systematics, biology, structure, ecological adaptations, evolution, genomics and the molecular processes that underpin the growth, development and survival of these important disease-transmitting parasites. Also discussed is the remarkable array of diseases transmitted (or caused) by ticks, as well as modern methods for their control. This book should serve as a modern reference for students, scientists, physicians, veterinarians and other specialists. Volume II includes chapters on the ecology of non-nidicolous and nidicolous ticks, genetics and genomics (including the genome of the Lyme disease vector Ixodes scapularis) and immunity, including host immune responses to tick feeding and tick-host interactions, as well as the tick's innate immune system that prevents and/or controls microbial infections. Six chapters cover in depth the many diseases caused by the major tick-borne pathogens, including tick-borne protozoa, viruses, rickettsiae of all types, other types of bacteria (e.g., the Lyme disease agent) and diseases related to tick paralytic agents and toxins. The remaining chapters are devoted to tick control using vaccines, acaricides, repellents, biocontrol, and, finally, techniques for breeding ticks in order to develop tick colonies for scientific study.

Book chapter↗

Mechanics of flow and sediment transport in delta distributary channels

Predicting the planform and dimensions of a channel downstream from a confluence of two smaller channels with known sediment and water supplies is a fundamental, well-studied problem in geomorphology and engineering. An analogous but less well understood problem is found well downstream of such confluences, where large river channels split into two or more distributary channels on a river delta. In this case, both the flow and sediment supplies in the downstream distributaries are set by the dynamics near the bifurcation of the upstream channel and by downstream boundary conditions. Over time, the pattern of erosion and deposition in the distributary channels gives rise to variations in the amount of water and sediment routed into them. In the simplest case, this results in channel switching on deltas, but in a more general sense these dynamics produce a rich suite of interesting morphologic change contributing both to the evolution of the channel distributary network and the overall evolution of the delta. As part of a study to develop a better understanding of these processes, we conducted a field study measuring the detailed morphology of the Hong-Luoc junction on the Red River downstream of Hanoi, Vietnam. This junction was selected for such a study because it has a 1,000-year history, modern observations suggest that it is currently switching (changing the proportion of sediment and streamflow provided to each of the distributary channels), and hydrologic configuration of the junction allows for the study of two bifurcations and one confluence in a single junction complex. In this paper, our morphologic observations are used in computational flow models to show how flow and sediment transport changes as a function of total discharge upstream of the junction. This is a key component of understanding how the junction attains stability over a range of flows or how imbalances in the distribution of flow and sediment transport lead to destabilization of the channel bifurcation.

Conference Paper↗

Hydrous pyrolysis of New Albany Shale: A study examining maturation changes and porosity development

The characterization of nanoscale organic structures has improved our understanding of porosity development within source-rock reservoirs, but research linking organic porosity evolution to thermal maturity has generated conflicting results. To better understand this connection, an immature (0.25% solid bitumen reflectance; BR o ) sample of the New Albany Shale was used in four isothermal hydrous pyrolysis (HP) experiment sequences at 300°, 320°, 340°, and 370°C, with residues collected periodically for a maximum of 103 days. The HP residues, along with the original immature sample and two naturally matured (1.49 and 1.56% BR o ) New Albany Shale samples were analyzed for organic petrology, total organic carbon (TOC) content, and organic porosity evaluation using correlative light and electron microscopy (CLEM). All of the HP series increased in thermal maturity with increasing duration of pyrolysis, though reflectance for each series plateaued within 25 days of maturation. Initially, TOC in the HP residues decreases (from 14.24 wt. %) with increasing thermal maturity until ∼1.0% BR o where TOC remains at ∼9–10 wt. % for all remaining residues. Qualitative CLEM observations within the 50–100 day 300° and 340°C HP sequences (0.95–1.70% BR o ), and the naturally matured samples, develop organic porosity in smaller (<5 μm in diameter), void-filling solid bitumen that occurs in spaces between clays and other fine-grained minerals. The 370°C HP residues developed significant organic porosity, relative to the other HP temperature series in all solid bitumen accumulations regardless of size. Overall, the study indicates that temperature and duration of artificial maturation play an important role in the abundance of pores in the HP residues. This work expands on our understanding of the conditions needed for the generation and development of organic porosity in the New Albany Shale and potentially in other marine source-rock petroleum systems.

Illinois, Indiana↗

Karst of the Mid-Atlantic region in Maryland, West Virginia, and Virginia

The Mid-Atlantic region hosts some of the most mature karst landscapes in North America, developed in highly deformed rocks within the Piedmont and Valley and Ridge physiographic provinces. This guide describes a three-day excursion to examine karst development in various carbonate rocks by following Interstate 70 west from Baltimore across the eastern Piedmont, across the Frederick Valley, and into the Great Valley proper. The localities were chosen in order to examine the structural and lithological controls on karst feature development in marble, limestone, and dolostone rocks with an eye toward the implications for ancient landscape evolution, as well as for modern subsidence hazards. A number of caves will be visited, including two commercial caverns that reveal strikingly different histories of speleogenesis. Links between karst landscape development, hydrologic dynamics, and water resource sustainability will also be emphasized through visits to locally important springs. Recent work on quantitative dye tracing, spring water geochemistry, and groundwater modeling reveal the interaction between shallow and deep circulation of groundwater that has given rise to the modern karst landscape. Geologic and karst feature mapping conducted with the benefit of lidar data help reveal the strong bedrock structural controls on karst feature development, and illustrate the utility of geologic maps for assessment of sinkhole susceptibility.

Book chapter↗

Detailed 40Ar/39Ar geochronology of the Loyalty and Three Kings Ridges clarifies the extent and sequential development of Eocene to Miocene southwest Pacific remnant volcanic arcs

The 2015 VESPA voyage (Volcanic Evolution of South Pacific Arcs) was a seismic and rock dredging expedition to the Loyalty and Three Kings Ridges and South Fiji Basin. In this paper we present 33 40 Ar/ 39 Ar, 22 micropaleontological, and two U/Pb ages for igneous and sedimentary rocks from 33 dredge sites in this little-studied part of the southwest Pacific Ocean. Igneous rocks include basalts, dolerites, basaltic andesites, trachyandesites, and a granite. Successful Ar/Ar dating of altered and/or low-K basalts was achieved through careful sample selection and processing, detailed petrographic and element mapping of groundmass, and incremental heating experiments on both phenocryst and groundmass separates to interpret the complex spectra produced by samples having multiple K reservoirs. The 40 Ar/ 39 Ar ages of most of the sampled lavas, irrespective of composition, are latest Oligocene to earliest Miocene (25–22 Ma); two are Eocene (39–36 Ma). The granite has a U/Pb zircon age of 23.6 ± 0.3 Ma. 40 Ar/ 39 Ar lava ages are corroborated by microfossil ages from associated sedimentary rocks. The VESPA lavas are part of a >3,000 km long disrupted belt of Eocene to Miocene subduction-related volcanic rocks. The belt includes arc rocks in Northland New Zealand, Northland Plateau, Three Kings Ridge, and Loyalty Ridge and, speculatively, D’Entrecasteaux Ridge. This belt is the product of superimposed Eocene and Oligocene-Miocene remnant volcanic arcs that were stranded along and near the edge of Zealandia while still-active arc belts migrated east with the Pacific trench.

Geochemistry, Geophysics, Geosystems↗

Probabilistic application of an integrated catchment-estuary-coastal system model to assess the evolution of inlet-interrupted coasts over the 21st century

Inlet-interrupted sandy coasts are dynamic and complex coastal systems with continuously evolving geomorphological behaviors under the influences of both climate change and human activities. These coastal systems are of great importance to society (e.g., providing habitats, navigation, and recreational activities) and are affected by both oceanic and terrestrial processes. Therefore, the evolution of these inlet-interrupted coasts is better assessed by considering the entirety of the Catchment-Estuary-Coastal (CEC) systems, under plausible future scenarios for climate change and increasing pressures due to population growth and human activities. Such a holistic assessment of the long-term evolution of CEC systems can be achieved via reduced-complexity modeling techniques, which are also ably quantifying the uncertainties associated with the projections due to their lower simulation times. Here, we develop a novel probabilistic modeling framework to quantify the input-driven uncertainties associated with the evolution of CEC systems over the 21 st century. In this new approach, probabilistic assessment of the evolution of inlet-interrupted coasts is achieved by (1) probabilistically computing the exchange sediment volume between the inlet-estuary system and its adjacent coast, and (2) distributing the computed sediment volumes along the inlet-interrupted coast. The model is applied at three case study sites: Alsea estuary (United States), Dyfi estuary (United Kingdom), and Kalutara inlet (Sri Lanka). Model results indicate that there are significant uncertainties in projected volume exchange at all the CEC systems (min-max range of 2.0 million cubic meters in 2100 for RCP 8.5), and the uncertainties in these projected volumes illustrate the need for probabilistic modeling approaches to evaluate the long-term evolution of CEC systems. A comparison of 50 th percentile probabilistic projections with deterministic estimates shows that the deterministic approach overestimates the sediment volume exchange in 2100 by 15–30% at Alsea and Kalutara estuary systems. Projections of coastline change obtained for the case study sites show that accounting for all key processes governing coastline change along inlet-interrupted coasts in computing coastline change results in projections that are between 20 and 134% greater than the projections that would be obtained if only the Bruun effect were taken into account, underlining the inaccuracies associated with using the Bruun rule at inlet-interrupted coasts.

Frontiers in Applied Mathematics and Statistics↗

Magmatic evolution of the Gulf of California rift

Magmatic events in the Baja peninsula-Gulf of California region are closely related to sequential tectonic regimes of: (1) late Tertiary subduction beneath the continental margin until ~16 to 12.5 Ma; and (2) continental to oceanic rifting that began about 13 Ma. Orogenic calcalkaline volcanic rocks formed subparallel belts of rhyolite ignimbrite of Oligocene age (~34 to 27 Ma) east of the Gulf in the Sierra Madre Occidental and andesite of Miocene age (~24 to 11 Ma) along the eastern Baja peninsula. Shutoff of the Miocene andesitic arc broadly corresponds to migration of the Pacific-Farallon (Guadalupe)-North America triple junction along the Baja peninsula. Orogenic andesitic volcanism ended at ~16 Ma in northern Baja, and at ~11 Ma in southern Baja, within ~1 to 2 Ma of cessation of subduction. Waning orogenic magmatism persisted along the southern Baja peninsula during the initial stages of rifting. Paleogeographic implications of the distribution of circum-Gulf volcanic rocks and their inferred easterly source areas suggest that, from about 13 to 8 Ma, the rift consisted of a narrow seaway along the eastern side of the present Gulf. By about 6 Ma, the Gulf had broadened westward to approximately its present width. The period from 13 to 10 Ma was a time of tectonic transition and magmatic diversity. During this interval, medium-K calcalkaline, high-K calcalkaline, alkalic, and tholeiitic magmas erupted from the central to southern part of the Baja peninsula-Gulf region, and from about 14 to 8 Ma, rhyolite ignimbrite erupted in the northern Gulf region. Since 10 Ma, volcanism on the western Baja peninsula has been dominated by alkalic magmas, the Gulf margins by calcalkaline magmas, and the Gulf by tholeiitic magmas. Postsubduction calcalkaline andesite to rhyolite erupted sporadically and locally along the northern Gulf margins from ~9 Ma through Holocene time. Alkalic magmas associated with rifting in the Gulf of California are unlike alkalic magmas from intraplate settings. They comprise nepheline- to quartz-normative, basaltic to andesitic rocks characterized by diverse trace element ratios. They are broadly similar to intraplate alkalic rocks in having high K, P, Ba, Sr and the light REE (rare earth elements), but have some trace element characteristics typical of orogenic magmas (low Nb and Ta relative to K, Ba, Sr, and La) that distinguish these alkalic rocks from intraplate or cratonic rift alkalic rocks. Rb abundances (mostly <25 ppm, as low as 5 ppm) are unusually low in these alkalic rocks, both in contrast to the high K (~1.7-4.1 wt% K 2 O), Ba (~800-2200 ppm), and Sr (~1000-3700 ppm), and in comparison to Rb abundances in mafic lavas in general; many samples have Rb at levels typically found in low-K tholeiites. The high K/Rb ratios, mostly ~1000 to 5000, are among the highest reported for any lavas. Samples ranging in SiO 2 from 48 to 61% have Mg’ ≥ 65, indicating that a broad range of compositions may represent primary melts. The generally high Mg’, and the low Rb, Th, Nb, and Ta in some lavas, suggest that depleted refractory mantle equivalent to a mid-ocean ridge basalt (MORB) source contributed to these melts. The distinctive trace element enrichments are attributed to the stabilization in the source of amphibole and apatite from metasomatic fluids related to prior subduction beneath the peninsula. In discrimination diagrams, these alkalic rocks typically cluster outside previously defined fields for various magma types. The common tectonic setting among these and geochemically similar alkalic suites is the occurrence in a continental margin that is undergoing extension, located along a recently or currently active convergent plate boundary. Tholeiitic basalts have erupted since the earliest stages of continental rifting through development of oceanic spreading centers. Gulf of California tholeiites exhibit an evolution distinguished particularly by differences in rare-earth-element (REE) abundances. Early-rift magmas are similar in trace element geochemistry to some intraplate tholeiites from ocean islands and continental flood basalts, and have convex-upward REE patterns. The most recent tholeiites from the East Pacific Rise in the mouth of the Gulf are equivalent to N-MORB (normal mid-ocean ridge basalt) from other mid-ocean rifts, and are depleted in the light-REE. The transition in tholeiite geochemistry from intraplate tholeiite to MORB is attributed to progressive fusion of more refractory mantle components. Incipient-rift tholeiites are derived from selective fusion of clinopyroxenite veins, and with continued melting, N-MORB are formed by fusion of peridotite. This tholeiite sequence reflects the evolution of mantle source regions in rifts that sustain an ensialic to oceanic transition.

Baja California, Baja California Sur↗

A compartmentalized solute transport model for redox zones in contaminated aquifers: 1. Theory and development

This paper, the first of two parts [see Abrams and Loague , this issue], takes the compartmentalized approach for the geochemical evolution of redox zones presented by Abrams et al . [1998] and embeds it within a solute transport framework. In this paper the compartmentalized approach is generalized to facilitate the description of its incorporation into a solute transport simulator. An equivalent formulation is developed which removes any discontinuities that may occur when switching compartments. Rate‐limited redox reactions are modeled with a modified Monod relationship that allows either the organic substrate or the electron acceptor to be the rate‐limiting reactant. Thermodynamic constraints are used to inhibit lower‐energy redox reactions from occurring under infeasible geochemical conditions without imposing equilibrium on the lower‐energy reactions. The procedure used allows any redox reaction to be simulated as being kinetically limited or thermodynamically limited, depending on local geochemical conditions. Empirical reaction inhibition methods are not needed. The sequential iteration approach (SIA), a technique which allows the number of solute transport equations to be reduced, is adopted to solve the coupled geochemical/solute transport problem. When the compartmentalized approach is embedded within the SIA, with the total analytical concentration of each component as the dependent variable in the transport equation, it is possible to reduce the number of transport equations even further than with the unmodified SIA. A one‐dimensional, coupled geochemical/solute transport simulation is presented in which redox zones evolve dynamically in time and space. The compartmentalized solute transport (COMPTRAN) model described in this paper enables the development of redox zones to be simulated under both kinetic and thermodynamic constraints. The modular design of COMPTRAN facilitates the use of many different, preexisting solute transport and geochemical codes. The companion paper [ Abrams and Loague , this issue] presents examples of the application of COMPTRAN to field‐scale problems.

Water Resources Research↗