Geology Reports⌕ Search

Geology topics

Jonathan T. Hagstrum

Publications and source records attributed to Jonathan T. Hagstrum.

At least 19 recordsLinked to original sources

Avian navigation: Comparing the olfactory navigational “map” and the infrasound direction-finding hypotheses to aeronautics

Animal navigation has long been a fascinating but bewildering subject. Humans and animals might well share similar navigational strategies because they developed within the same physical environments. A “map-and-compass” model has been proposed to explain the two-step avian navigational process, but the “map” step has remained elusive. Although scalar values from bicoordinate geomagnetic or atmospheric olfactory gradients have been considered foundational to the avian map, neither has proved convincing engendering decades of controversy. The olfactory map, and an alternative infrasound direction-finding (IDF) hypothesis, are discussed in this review. The olfactory map hypothesis currently requires extensive stable gradients of trace-odor ratios, but such gradients are highly unlikely within a turbulent and rapidly mixed lower atmosphere. The IDF hypothesis, on the other hand, postulates a two-step navigational model analogous to the maritime and aeronautical radio direction-finding technique. This review was also written to encourage further investigation, and direct testing, of the acoustic navigational process. The IDF hypothesis, at present, appears the better explanation of observed avian navigational behavior and accuracy within the atmosphere’s physical environment.

Journal of Comparative Physiology A↗

New paleomagnetic constraints on the eruption timing, stratigraphy, and post-emplacement deformation of the Picture Gorge Basalt within the Columbia River Basalt Group

The Picture Gorge Basalt (PGB) is part of the Miocene Columbia River Basalt Group (CRBG). The PGB, which outcrops in eastern Oregon, is considered coincident in time with the voluminous Grande Ronde Basalt. New radiometric ages have expanded the age‐range of the PGB, suggesting it began erupting prior to the Steens Basalt to its south and continued until after cessation of the Grande Ronde Basalt eruptions, an interval of 1.5 Ma. However, the existing paleomagnetism of the PGB implies this eruption timeline is an overestimate. To reconcile the radiometric and paleomagnetic timescales for the PGB, we conducted a paleomagnetic study on sections of the PGB to construct a detailed, high‐quality magnetostratigraphy. Our data indicate the stratigraphically lowest lava flows in the PGB are of reversed polarity, revealing a new paleomagnetic transition with the PGB and a reversed (R)–normal (N)–reversed (R) sequence. This suggests one of two timeline possibilities for PGB volcanism: (a) eruptions began and during through CRBG polarity chrons R0–N0–R1, penecontemporaneous with Steens Basalt, or (b) eruptions began and persisted during CRBG polarity chrons R1–N1–R2. Our work supports a longer interval of PGB volcanism than was suggested by previous paleomagnetic data but is at odds with the suggestion that PGB eruptions lasted through the entire main CRBG. We favor a scenario wherein PGB eruptions begin with R0 and continue into the R1 paleomagnetic interval. The paleomagnetic results also record a ∼18° vertical‐axis rotation of east‐central Oregon after ∼16 Ma with respect to stable North America.

Oregon↗

A reinterpretation of “Homing pigeons’ flight over and under low stratus” based on atmospheric propagation modeling of infrasonic navigational cues

Pigeons flying above temperature inversion and related low-stratus layers appear to lack important navigational cues, and a reinterpretation of Wagner’s 1978 study suggests that these cues are low-frequency acoustic signals (infrasound). Wagner released homing pigeons above opaque stratus over the Swiss Plateau to determine whether they could locate their loft beneath it. Birds above the clouds appeared lost, while those that descended beneath them returned home directly. Atmospheric propagation modeling of infrasonic waves virtually transmitted from the loft area shows that these signals would have been ducted beneath the inversion layer, and would not have reached the release sites above it. The absence of homeward infrasonic cues above temperature inversions could explain the disorientation of Wagner’s birds, especially if such signals are the predominant cues used by pigeons to home. The possible generation of infrasonic navigational signals in the loft area and recent queries concerning the infrasound navigational “map” hypothesis are also discussed.

Journal of Comparative Physiology A↗

Impact-related microspherules in Late Pleistocene Alaskan and Yukon “muck” deposits signify recurrent episodes of catastrophic emplacement

Large quantities of impact-related microspherules have been found in fine-grained sediments retained within seven out of nine, radiocarbon-dated, Late Pleistocene mammoth ( Mammuthus primigenius ) and bison ( Bison priscus ) skull fragments. The well-preserved fossils were recovered from frozen “muck” deposits (organic-rich silt) exposed within the Fairbanks and Klondike mining districts of Alaska, USA, and the Yukon Territory, Canada. In addition, elevated platinum abundances were found in sediment analysed from three out of four fossil skulls. In view of this new evidence, the mucks and their well-preserved but highly disrupted and damaged vertebrate and botanical remains are reinterpreted in part as blast deposits that resulted from several episodes of airbursts and ground/ice impacts within the northern hemisphere during Late Pleistocene time (~46–11 ka B.P.). Such a scenario might be explained by encounters with cometary debris in Earth-crossing orbits (Taurid Complex) that was generated by fragmentation of a large short-period comet within the inner Solar System.

Alaska, Yukon Territoriy↗

Paleomagnetism and 40Ar/39Ar geochronology of the Plio-Pleistocene Boring Volcanic Field: Implications for the geomagnetic polarity time scale and paleosecular variation

Paleomagnetic directions and 40 Ar/ 39 Ar ages have been determined for samples of lava flows from the same outcrops, where possible, for 84 eruptive units ranging in age from 3200 ka to 60 ka within the Boring Volcanic Field (BVF) of the Pacific Northwest, USA. This study expands upon our previous results for the BVF, and compares the combined results with the current geomagnetic polarity time scale (GPTS). Lava flows with transitional directions were found within the BVF at the Matuyama-Brunhes and Jaramillo-Matuyama polarity boundaries, and replicate ages corresponding to these and other boundaries have been newly ascertained. Although the BVF data generally agree with GPTS chronozone boundaries, they indicate that onset of the Gauss-Matuyama transition and Olduvai subchron occurred significantly earlier than given in the current time scale calibration. Additional comparisons show that the BVF results are consistent with recent statistical models of geomagnetic paleosecular variation.

Physics of the Earth and Planetary Interiors↗

Seasonal changes in atmospheric noise levels and the annual variation in pigeon homing performance

Repeated releases of experienced homing pigeons from single sites were conducted between 1972 and 1974 near Cornell University in upstate New York and between 1982 and 1983 near the University of Pittsburgh in western Pennsylvania, USA. No annual variation in homing performance was observed at these sites in eastern North America, in contrast to results from a number of similar experiments in Europe. Assuming pigeons home using low-frequency infrasonic signals (~0.1–0.3 Hz), as has been previously proposed, the annual and geographic variability in homing performance within the northern hemisphere might be explained, to a first order, by seasonal changes in low-frequency atmospheric background noise levels related to storm activity in the North Atlantic Ocean, and by acoustic waveguides formed between the surface and seasonally reversing stratospheric winds. In addition, increased dispersion among departure bearings of test birds on some North American release days was possibly caused by infrasonic noise from severe weather events during tornado and Atlantic hurricane seasons.

Journal of Comparative Physiology A↗

40 Ar/ 39 Ar geochronology, paleomagnetism, and evolution of the Boring volcanic field, Oregon and Washington, USA

The 40 Ar/ 39 Ar investigations of a large suite of fine-grained basaltic rocks of the Boring volcanic field (BVF), Oregon and Washington (USA), yielded two primary results. (1) Using age control from paleomagnetic polarity, stratigraphy, and available plateau ages, 40 Ar/ 39 Ar recoil model ages are defined that provide reliable age results in the absence of an age plateau, even in cases of significant Ar redistribution. (2) Grouping of eruptive ages either by period of activity or by composition defines a broadly northward progression of BVF volcanism during latest Pliocene and Pleistocene time that reflects rates consistent with regional plate movements. Based on the frequency distribution of measured ages, periods of greatest volcanic activity within the BVF occurred 2.7–2.2 Ma, 1.7–0.5 Ma, and 350–50 ka. Grouped by eruptive episode, geographic distributions of samples define a series of northeast-southwest–trending strips whose centers migrate from south-southeast to north-northwest at an average rate of 9.3 ± 1.6 mm/yr. Volcanic activity in the western part of the BVF migrated more rapidly than that to the east, causing trends of eruptive episodes to progress in an irregular, clockwise sense. The K 2 O and CaO values of dated samples exhibit well-defined temporal trends, decreasing and increasing, respectively, with age of eruption. Divided into two groups by K 2 O, the centers of these two distributions define a northward migration rate similar to that determined from eruptive age groups. This age and compositional migration rate of Boring volcanism is similar to the clockwise rotation rate of the Oregon Coast Range with respect to North America, and might reflect localized extension on the trailing edge of that rotating crustal block.

Oregon, Washington↗

Atmospheric propagation modeling indicates homing pigeons use loft-specific infrasonic ‘map’ cues

Results from an acoustic ray-tracing program using daily meteorological profiles are presented to explain ‘release-site biases’ for homing pigeons at three experimental sites in upstate New York where W. T. Keeton and his co-workers at Cornell University conducted extensive releases between 1968 and 1987 in their investigations of the avian navigational ‘map’. The sites are the Jersey Hill and Castor Hill fire towers, and another near Weedsport, where control pigeons from the Cornell loft vanished in random directions, in directions consistently >50 deg clockwise and in directions ∼15 deg clockwise from the homeward bearing, respectively. Because Cornell pigeons were disoriented at Jersey Hill whereas birds from other lofts were not, it is inferred that Jersey Hill lies within an acoustic ‘shadow’ zone relative to infrasonic signals originating from the Cornell loft’s vicinity. Such signals could arise from ground-to-air coupling of near-continuous microseisms, or from scattering of direct microbaroms off terrain features, both of which are initially generated by wave–wave interactions in the deep ocean. HARPA runs show that little or no infrasound from the loft area arrived at Jersey Hill on days when Cornell pigeons were disoriented there, and that homeward infrasonic signals could have arrived at all three sites from directions consistent with pigeon departure bearings, especially on days when these bearings were unusual. The general stability of release-site biases might be due to influences of terrain on transmission of the homeward signals under prevailing weather patterns, whereas short-term changes in biases might be caused by rapid shifts in atmospheric conditions.

New York↗

Paleomagnetism of three upper Jurassic ash-flow sheets in southeastern Arizona: Implications for regional deformation

A sequence of three petrologically distinct Upper Jurassic ash-flow tuffs has been sampled for paleomagnetic analysis in the Huachuca Mountains, Canelo Hills, and Mustang Mountains of southeastern Arizona. Site-mean paleomagnetic directions for these units indicate significant vertical-axis rotations between sampling localities (≃15°–40°), and a comparison of these and previously published data with reference directions from the Colorado Plateau implies that the Jurassic rocks in southeastern Arizona have been rotated in a clockwise sense with respect to stable North America. These new data are consistent with paleomagnetic results showing clockwise rotation of Upper Cretaceous rocks in south central Arizona, and further indicate that studies of North American apparent polar wander based on rocks from southern Arizona could be subject to error. The clockwise rotations in southern Arizona are likely related to Late Cretaceous and early Tertiary strike-slip movement on northwest-trending high-angle faults in the region, and this movement may be linked to an increase in oblique Farallon-North American plate convergence at that time.

Arizona↗

Paleomagnetic and 40Ar/39Ar evidence for remagnetization of Mesozoic oceanic rocks on the Vizcaino Peninsula, Baja California Sur, Mexico

Previously published paleomagnetic data for Upper Jurassic pillow lavas of the Vizcaino Peninsula indicate that they were deposited near a paleolatitude of 14°N or S. Whether or not this result implies northward transport with respect to the continental interior has been controversial due to the lack of reliable Jurassic reference poles for the North American plate. Available paleomagnetic data for nearby Upper Triassic pillow basalts and overlying pelagic sediments at Punta San Hipólito, however, fail a fold test indicating that these rocks were remagnetized post-folding. Indistinguishable in-situ paleomagnetic directions and perturbed 40 Ar/ 39 Ar age spectra for the Triassic and Jurassic pillow lavas are consistent with resetting of their magnetic and isotopic systems in the middle Cretaceous, probably during burial by the overlap Valle Formation (>10 km thick). Resetting apparently occurred post-deformation so the paleolatitude of remagnetization is unknown. High-coercivity directions from a few samples of the Triassic rocks might represent an older magnetization acquired during deposition at paleolatitudes near 6°N or S.

Geophysical Research Letters↗

Faulting within the Mount St. Helens conduit and implications for volcanic earthquakes

The 2004–2008 eruption of Mount St. Helens produced seven dacite spines mantled by cataclastic fault rocks, comprising an outer fault core and an inner damage zone. These fault rocks provide remarkable insights into the mechanical processes that accompany extrusion of degassed magma, insights that are useful in forecasting dome-forming eruptions. The outermost part of the fault core consists of finely comminuted fault gouge that is host to 1- to 3-mm-thick layers of extremely fine-grained slickenside-bearing ultracataclasite. Interior to the fault core, there is an ∼2-m-thick damage zone composed of cataclastic breccia and sheared dacite, and interior to the damage zone, there is massive to flow-banded dacite lava of the spine interior. Structures and microtextures indicate entirely brittle deformation, including rock breakage, tensional dilation, shearing, grain flow, and microfaulting, as well as gas and fluid migration through intergranular pores and fractures in the damage zone. Slickenside lineations and consistent orientations of Riedel shears indicate upward shear of the extruding spines against adjacent conduit wall rocks. Paleomagnetic directions, demagnetization paths, oxide mineralogy, and petrology indicate that cataclasis took place within dacite in a solidified steeply dipping volcanic conduit at temperatures above 500 °C. Low water content of matrix glass is consistent with brittle behavior at these relatively high temperatures, and the presence of tridymite indicates solidification depths of <1 km. Cataclasis was coincident with the eruption’s seismogenic zone at <1.5 km. More than a million small and low-frequency “drumbeat” earthquakes with coda magnitudes (M d ) <2.0 and frequencies <5 Hz occurred during the 2004–2008 eruption. Our field data provide a means with which to estimate slip-patch dimensions for shear planes and to compare these with estimates of slip patches based on seismic moments and shear moduli for dacite rock and granular fault gouge. Based on these comparisons, we find that aseismic creep is achieved by micron-scale displacements on Riedel shears and by granular flow, whereas the drumbeat earthquakes require millimeter to centimeter displacements on relatively large (e.g., ∼1000 m 2 ) slip patches, possibly along observed extensive principal shear zones within the fault core but probably not along the smaller Riedel shears. Although our field and structural data are compatible with stick-slip models, they do not rule out seismic and infrasound models that call on resonance of steam-filled fractures to generate the drumbeat earthquakes. We suggest that stick-slip and gas release processes may be coupled, and that regardless of the source mechanism, the distinctive drumbeat earthquakes are proving to be an effective precursor for dome-forming eruptions. Our data document a continuous cycle of deformation along the conduit margins beginning with episodes of fracture in the damage zone and followed by transfer of motion to the fault core. We illustrate the cycle of deformation using a hypothetical cross section of the Mount St. Helens conduit, extending from the surface to the depth of magmatic solidification.

Washington↗

North American Jurassic APW: The current dilemma

Geologists are in a quandary over the correct interpretation of paleomagnetic data for Jurassic rocks of the North American plate. Conflicting reference paleopoles and alternate methods of constructing apparent polar wander (APW) paths have led to a controversy regarding the configuration of Jurassic APW for cratonic North America. These differences have been recently disputed at meetings (see Eos , Spring Meeting Supplement, April 7, 1992, p. 94) and in an exchange of letters in the Journal of Geophysical Research-Solid Earth [ Butler et al ., 1992; Van Fossen and Kent , 1992a]. At stake is important information concerning the Jurassic paleogeography of North America and the whole of Pangea, as well as the nature and driving mechanism of plate movements. In addition, cratonic paleopoles provide a reference frame for measuring relative displacements of tectonostratigraphic terranes.

Eos Science News↗

Remagnetization and northward coastwise transport of Franciscan Complex rocks, northern California: A reinterpretation of the paleomagnetic data

Previous paleomagnetic studies of accreted oceanic rocks within the Franciscan Complex of northern California have concluded that these rocks originated far to the south of their present positions with respect to the North American continent. Based on positive “fold” tests, the characteristic remanent magnetizations were inferred to predate accretion-related deformation and metamorphism, and to have been acquired during or soon after deposition of these rocks. Thus, the paleomagnetic data were thought to provide direct information on ancient oceanic plate motions. However, the plate motions implied by some of these paleomagnetic data are problematic (e.g., exceptional plate velocities), and uniform-polarity magnetizations in almost all of these rocks indicate the possibility of remagnetization. Recent work on oceanic rocks in similar subduction complexes of Japan and Mexico have shown that they were most likely chemically remagnetized during accretion prior to disruption of the original stratigraphic sequences. Modern analogs indicate that the oceanic rocks in Mexico were probably remagnetized while still part of a shallow-dipping subducting slab (<10°) at the base of an accretionary prism. Assuming these rocks were near horizontal at the time of remagnetization, paleolatitudes at which these rocks were subducted and subsequent arc-parallel displacements along the western margin of North America can be inferred. In this paper, Franciscan rocks in northern California are reinterpreted as also having been remagnetized prior to accretion-related deformation. This scenario satisfies both geologic and paleomagnetic constraints for these rocks, and resolves conflicts between data indicating both remagnetization and tectonic displacement. Transport of the Laytonville Limestone from the southern hemisphere is not required. Paleolatitudes of subduction and remagnetization in the northern hemisphere (12° to 33°) appear to be inversely proportional to age of accretion (middle Cretaceous to Oligocene) for the Franciscan rocks. Subsequent northward diplacements (800 to 3700 km) and clockwise rotations (56° and 154°) of these rocks inferred from the paleomagnetic data are consistent with potential displacements along the western margin of North America during late Mesozoic and Cenozoic time calculated using examples of modern subduction zones and current plate reconstruction models.

California↗

Archeomagnetic dating in western North America: An updated reference curve based on paleomagnetic and archeomagnetic data sets

[1] A robust database of paleomagnetic (PM) and archeomagnetic (AM) data for western North America covering the past 4 millennia is assembled from three data sets: (1) published PM data for dated lava flows, (2) recently published AM data generated by R. L. DuBois, and (3) published AM data generated by J. L. Eighmy. Chronologic information for the AM data sets has been reassessed, and the data points have been assigned to age ranges depending on their archeological contexts. We have selected the computer program ArcheoCourbe, available from M. Le Goff (http://www.ipgp.fr/∼legoff) or the auxiliary material (see also http://EarthRef.org), to construct AM reference curves and to quantitatively estimate AM dates for archeological and geological samples. The updated curve provides a more detailed record of prehistoric secular variation for western North America than was previously available and shows a short interval of rapid secular variation (SV) in the 10th or 11th centuries at near-excursional rates. Similar and coeval features in the Hawaiian and, to a lesser extent, Arkansan SV data sets possibly indicate a rapid amplitudinal shift in the regional nondipole field at that time.

Arizona, California, Colorado, Idaho, New Mexico, ↗

The Boring Volcanic Field of the Portland-Vancouver area, Oregon and Washington: Tectonically anomalous forearc volcanism in an urban setting

More than 80 small volcanoes are scattered throughout the Portland-Vancouver metropolitan area of northwestern Oregon and southwestern Washington. These volcanoes constitute the Boring Volcanic Field, which is centered in the Neogene Portland Basin and merges to the east with coeval volcanic centers of the High Cascade volcanic arc. Although the character of volcanic activity is typical of many monogenetic volcanic fields, its tectonic setting is not, being located in the forearc of the Cascadia subduction system well trenchward of the volcanic-arc axis. The history and petrology of this anomalous volcanic field have been elucidated by a comprehensive program of geologic mapping, geochemistry, 40Ar/39Ar geochronology, and paleomag-netic studies. Volcanism began at 2.6 Ma with eruption of low-K tholeiite and related lavas in the southern part of the Portland Basin. At 1.6 Ma, following a hiatus of ~0.8 m.y., similar lavas erupted a few kilometers to the north, after which volcanism became widely dispersed, compositionally variable, and more or less continuous, with an average recurrence interval of 15,000 yr. The youngest centers, 50&ndash;130 ka, are found in the northern part of the field. Boring centers are generally monogenetic and mafic but a few larger edifices, ranging from basalt to low-SiO2 andesite, were also constructed. Low-K to high-K calc-alkaline compositions similar to those of the nearby volcanic arc dominate the field, but many centers erupted magmas that exhibit little influence of fluids derived from the subducting slab. The timing and compositional characteristics of Boring volcanism suggest a genetic relationship with late Neogene intra-arc rifting.

Portland, Vancouver↗

The Columbia River Basalt Group: from the gorge to the sea

Miocene flood basalts of the Columbia River Basalt Group inundated eastern Washington, Oregon, and adjacent Idaho between 17 and 6 Ma. Some of the more voluminous flows followed the ancestral Columbia River across the Cascade arc, Puget-Willamette trough, and the Coast Range to the Pacific Ocean. We have used field mapping, chemistry, and paleomagnetic directions to trace individual flows and flow packages from the Columbia River Gorge westward into the Astoria Basin, where they form pillow palagonite complexes and mega-invasive bodies into older marine sedimentary rocks. Flows of the Grande Ronde, Wanapum, and Saddle Mountains Basalts all made it to the ocean; at least 33 flows are recognized in the western Columbia River Gorge, 50 in the Willamette Valley, 16 in the lower Columbia River Valley, and at least 12 on the Oregon side of the Astoria Basin. In the Astoria Basin, the basalt flows loaded and invaded the wet marine sediments, producing peperite breccias, soft sediment deformation, and complex invasive relations. Mega-invasive sills up to 500 m thick were emplaced into strata as old as Eocene, and invasive dikes up to 90 m thick can be traced continuously for 25 km near the basin margin. Mega-pillow complexes up to a kilometer thick are interpreted as the remains of lava deltas that prograded onto the shelf and a filled submarine canyon southeast of Astoria, possibly providing the hydraulic head for injection of invasive sills and dikes at depth.

Oregon;Washington↗

Distribution, 14 C chronology, and paleomagnetism of latest Pleistocene and Holocene lava flows at Haleakala volcano, Island of Maui, Hawai'i: A revision of lava flow hazard zones

New mapping and 60 new radiocarbon ages define the age and distribution of latest Pleistocene and Holocene (past 13,000 years) lava flows at Haleakalā volcano, Island of Maui. Paleomagnetic directions were determined for 118 sites, of which 89 are in lava flows younger than 13,000 years. The paleomagnetic data, in conjunction with a reference paleosecular variation (PSV) curve for the Hawaiian Islands, are combined with our knowledge of age limitations based on stratigraphic control to refine age estimates for some of the undated lava flows. The resulting volumetric rate calculations indicate that within analytical error, the extrusion rate has remained nearly constant during the past 13,000 years, in the range 0.05&ndash;0.15 km 3 /kyr, only about half the long-term rate required to produce the postshield strata emplaced in the past &sim;1 Myr. Haleakalā's eruptive frequency is similar to that of Hualālai volcano on the Island of Hawai&lsquo;i, but its lava flows cover substantially less area per unit time. The reduced rates of lava coverage indicate a lower volcanic hazard than in similar zones at Hualālai.

Hawaii↗

Antipodal hotspots and bipolar catastrophes: Were oceanic large-body impacts the cause?

One aspect of the hotspot distribution that has received little attention is its antipodal character. Of 45 ‘primary’ hotspots found in most hotspot compilations 22 (49%) form antipodal pairs within observed hotspot drift limits (≤ 20 mm/yr). In addition, the available ages, or possible age ranges, for both hotspots of an antipodal pair tend to be similar (≤ 10 Myr difference) or overlap. Monte Carlo simulations indicate that the antipodal primary hotspots' locations and ages are not due to chance at the > 99% confidence level ( p < 0.01). All hotspot pairs include at least one oceanic hotspot, and these are consistently opposite those hotspots related to large igneous provinces (LIPs) and continental volcanism. A mechanism of formation is considered in which minor hotspot volcanism is induced at, and flood basalt volcanism is triggered by seismic energy focused antipodal to, oceanic large-body impact sites. Because continental impacts are expected to have lower seismic efficiencies, continents possibly acted as shields to the formation of antipodal hotspot pairs. Published numerical models indicate that large oceanic impacts (10-km-diameter bolide) generate megatsunami capable of altering coastal depositional environments on a global scale. Past impact-generated megatsunami, consequently, could have left widespread stratigraphic records, possibly misinterpreted as indicating large rapid changes in eustatic sea level, and widely disrupted continental and marine sediment reservoirs responsible for abrupt changes in the isotopic composition of seawater. Phanerozoic mass extinction events, therefore, might have resulted primarily from catastrophic megatsunami in a dominantly oceanic hemisphere and the near contemporaneous effusion of vast quantities of noxious gases from flood basalt eruptions in a dominantly continental one.

Earth and Planetary Science Letters↗