Geology Reports⌕ Search

SEARCH · Geology Reports

Results for “Holocene”

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 685 records · Page 38Linked to original sources

The Quaternary Silver Creek Fault Beneath the Santa Clara Valley, California

The northwest-trending Silver Creek Fault is a 40-km-long strike-slip fault in the eastern Santa Clara Valley, California, that has exhibited different behaviors within a changing San Andreas Fault system over the past 10-15 Ma. Quaternary alluvium several hundred meters thick that buries the northern half of the Silver Creek Fault, and that has been sampled by drilling and imaged in a detailed seismic reflection profile, provides a record of the Quaternary history of the fault. We assemble evidence from areal geology, stratigraphy, paleomagnetics, ground-water hydrology, potential-field geophysics, and reflection and earthquake seismology to determine the long history of the fault in order to evaluate its current behavior. The fault formed in the Miocene more than 100 km to the southeast, as the southwestern fault in a 5-km-wide right step to the Hayward Fault, within which the 40-km-long Evergreen pull-apart basin formed. Later, this basin was obliquely cut by the newly recognized Mt. Misery Fault to form a more direct connection to the Hayward Fault, although continued growth of the basin was sufficient to accommodate at least some late Pliocene alluvium. Large offset along the San Andreas-Calaveras-Mt Misery-Hayward Faults carried the basin northwestward almost to its present position when, about 2 Ma, the fault system was reorganized. This led to near abandonment of the faults bounding the pull-apart basin in favor of right slip extending the Calaveras Fault farther north before stepping west to the Hayward Fault, as it does today. Despite these changes, the Silver Creek Fault experienced a further 200 m of dip slip in the early Quaternary, from which we infer an associated 1.6 km or so of right slip, based on the ratio of the 40-km length of the strike-slip fault to a 5-km depth of the Evergreen Basin. This dip slip ends at a mid-Quaternary unconformity, above which the upper 300 m of alluvial cover exhibits a structural sag at the fault that we interpret as a negative flower structure. This structure implies some continuing strike slip on the Silver Creek Fault in the late Quaternary as well, with a transtensional component but no dip slip. Our only basis for estimating the rate of this later Quaternary strike slip on the Silver Creek Fault is to assume continuation of the inferred early Quaternary rate of less than 2 mm/yr. Faulting evident in a detailed seismic reflection profile across the Silver Creek Fault extends up to the limit of data at a depth of 50 m and age of about 140 ka, and the course of Coyote Creek suggests Holocene capture in a structural depression along the fault. No surface trace is evident on the alluvial plain, however, and convincing evidence of Holocene offset is lacking. Few instrumentally recorded earthquakes are located near the fault, and those that are near its southern end represent cross-fault shortening, not strike slip. The fault might have been responsible, however, for two poorly located moderate earthquakes that occurred in the area in 1903. Its southeastern end does mark an abrupt change in the pattern of abundant instrumentally recorded earthquakes along the Calaveras Fault-in both its strike and in the depth distribution of hypocenters-that could indicate continuing influence by the Silver Creek Fault. In the absence of convincing evidence to the contrary, and as a conservative estimate, we presume that the Silver Creek Fault has continued its strike-slip movement through the Holocene, but at a very slow rate. Such a slow rate would, at most, yield very infrequent damaging earthquakes. If the 1903 earthquakes did, in fact, occur on the Silver Creek Fault, they would have greatly reduced the short-term future potential for large earthquakes on the fault.

Open-File Report↗

Apalachicola Bay interpreted seismic horizons and updated IRIS chirp seismic-reflection data

Apalachicola Bay and St. George Sound contain the largest oyster fishery in Florida, and the growth and distribution of the numerous oyster reefs here are the combined product of modern estuarine conditions and the late Holocene evolution of the bay. A suite of geophysical data and cores were collected during a cooperative study by the U.S. Geological Survey, the National Oceanic and Atmospheric Administration Coastal Services Center, and the Apalachicola National Estuarine Research Reserve to refine the geology of the bay floor as well as the bay's Holocene stratigraphy. Sidescan-sonar imagery, bathymetry, high-resolution seismic profiles, and cores show that oyster reefs occupy the crests of sandy shoals that range from 1 to 7 kilometers in length, while most of the remainder of the bay floor is covered by mud. The sandy shoals are the surficial expression of broader sand deposits associated with deltas that advanced southward into the bay between 6,400 and 4,400 years before present. The seismic and core data indicate that the extent of oyster reefs was greatest between 2,400 and 1,200 years before present and has decreased since then due to the continued input of mud to the bay by the Apalachicola River. The association of oyster reefs with the middle to late Holocene sandy delta deposits indicates that the present distribution of oyster beds is controlled in part by the geologic evolution of the estuary.

Florida↗

Regional and other general factors bearing on evaluation of earthquake and other geologic hazards to coastal communities of southeastern Alaska

The great Alaska earthquake of March 27, 1964, brought into sharp focus the need for engineering geologic studies in seismically active regions. As a result, nine communities in southeastern Alaska were selected for reconnaissance investigations as an integral part of an overall program to evaluate earthquake and other geologic hazards in most of the larger Alaska coastal communities. This report gives background information on the regional and other general factors that bear on these evaluations. Southeastern Alaska, about 525 miles long and averaging about 125 miles in width, consists of a narrow mainland strip and numerous islands. For the most part, it is a region of rugged relief with numerous glaciers capping many of the higher mountainous areas and with long linear fiords forming the inland waterways. A maritime climate prevails with mild winters and cool summers. The southeastern part of the region receives the highest precipitation in the continental United States. Ketchikan, with a population of 6,994 in 1970, is the largest city. Geology and structure of the area are complex. Igneous, metamorphic, and sedimentary rocks crop out and range in age from Paleozoic to Tertiary. Surficial deposits of Pleistocene and Holocene age mantle many areas. All of southeastern Alaska, except probably the highest peaks, was covered by glacier ice advances of late Pleistocene age. Major deglaciation was well advanced by 10,000 years ago--a time which approximately marks the end of the Pleistocene and the beginning of the Holocene. There followed a period of warm climate called the Hypsithermal, which in southeastern Alaska began 7,000-8,000 years ago and ended about 4,800-3,500 years ago. Glaciers in most places receded back of their present positions. The Hypsithermal was followed by an interval (termed Neoglaciation) of cooler climate and resurgence of glacier ice which continues to the present, although most glaciers are now rapidly receding. During the past 10,000 years worldwide sea level has risen about 100 feet, but during the past 4,000 years it has risen only about 10 feet or about 0.03 inch per year. With sea level used as a datum, the amount of sea-level rise must be added to the apparent uplift of land for the time under consideration to determine the actual amount of land uplift. The widespread presence of emergent marine deposits, several hundred feet above sea level, demonstrates that the land in southeastern Alaska has been uplifted since the last major deglaciation. The greatest known has been uplifted since the last major deglaciation. The greatest known uplift is in the vicinity of Juneau where glaciomarine deposits are present 750 feet above present sea level. Part of southeastern Alaska is presently undergoing one of the most rapid rates of uplift of any place in the world. The fastest emergence is occurring in the Glacier Bay area where the land is being uplifted relative to sea level approximately 3.9 cm per year. Most or all of the uplift appears to be due to rebound as a result of deglaciation. Southeastern Alaska lies within the circum-Pacific earthquake belt, one of the world's greatest zones of seismic activity. During historic time, there have been five earthquakes in the region with magnitudes of 8 or greater, three with magnitudes of 7 to 8, eight with magnitudes of 6 to 7, more than 15 with magnitudes of 5 to 6, and about 140 recorded earthquakes with magnitudes smaller than 5 or of unassigned magnitudes. All of the earthquakes with magnitudes 8 or greater, and a large proportion of the others, appear to be related to the active Fairweather- Queen Charlotte Islands fault system or its western extension, the Chugach-St. Elias fault. Earthquake epicenters on the Denali fault system, the other major fault system in southeastern Alaska, are few in comparison. However, because high microearthquake activity has been recorded recently on this system and earthquakes of moderate size have occurred on some of its segments, the Denali fault system probably should not be dismissed as a relict fault system of no current tectonic importance. There are numerous other known faults, as well as lineaments that may be faults of varying degrees of tectonic activity in southeastern Alaska, adjacent Canada, and eastern Alaska. One of these elements is the Totschunda fault system, which connects with the Denali fault system in eastern Alaska; it has been very active during Holocene time but few historical earthquake epicenters appear to be related to it. Both historical seismicity and geologic conditions, such as frequency and recency of faulting, must be considered together to permit an assessment of the future earthquake probability of an area. Data are too few for both factors for an accurate evaluation to be made of earthquake probability in southeastern Alaska. However, information compiled in the form of strain-release and seismic-zone maps permit some generalizations. Thus, it is tentatively concluded that most, if not all, of southeastern Alaska should be placed in seismic zone 3, a zone in which earthquakes of magnitude greater than 6 will occur from time to time and where there may be major damage to manmade structures. Inferred effects from future earthquakes in southeastern Alaska include: (1) surface displacement along faults and other tectonic land-level changes, (2) ground shaking, (3) compaction, (4) liquefaction in cohesionless materials, (5) reaction of sensitive and quick clays, (6) water-sediment ejection and associated subsidence and ground fracturing, (7) earthquake-induced sub aerial slides and slumps, (8) earthquake induced subaqueous slides, (9) effects on glaciers and related features, (10) effects on ground water and stream flow, and (11) tsunamis, seiches, and other abnormal water waves. Because of the reconnaissance nature of our studies in the coastal communities and the sparsity of laboratory data on physical properties of geologic units in each area studied, the inferred effects must be largely empirical and generalized. Therefore, the inferences are based in large part upon the effects of past major earthquakes in Alaska and elsewhere, particularly upon the well-documented effects of the Alaska earthquake of March 27, 1964. Buildings, highways, bridges, tunnels, harbor facilities, pipelines, canals, and other manmade structures may be severely damaged or destroyed by fault displacement or related tectonic land-level changes in southeastern Alaska. Direct damage from fault rupture would be restricted virtually to structures built directly athwart the fault. In California and Nevada, fault rupture almost always accompanies shocks of magnitude 6.5 or greater. The Alaska earthquake of March 27, 1964, and the Chilean earthquake of May 22, 1960, dramatically illustrated the severe adverse effects that can result from uplift or subsidence over a wide area. The variable most responsible for the degree of shaking at any epicentral distance is the type of ground. Generally, shaking is considerably greater in poorly consolidated deposits than in hard bedrock, particularly if the deposits are water saturated. Severe shaking of alluvial deposits and manmade fill, with resultant heavy damage, is well documented from the records of many past earthquakes. Damage commonly has been heavy as a result of ground settlement caused by compaction of loose sediments by shaking during an earthquake. This has been especially true where compaction was accompanied by tectonic downdrop of land, such as occurred during the Chilean earthquake of 1960 and the Alaska earthquake of 1964. Loosely emplaced manmade fill, deltaic deposits, beach deposits, and alluvial deposits may be susceptible to compaction in southeastern Alaska during a severe earthquake. Liquefaction of sand and silt is a fairly common effect of large earthquakes. It was well illustrated at Niigata, Japan, during the earthquake of June 16, 1964, and resulted in extensive damage. When part of a sloping soil mass liquefies, the entire mass can undergo catastrophic failure and can flow as a high-density liquid. In southeastern Alaska, deltaic deposits probably would be most susceptible to liquefaction. Sensitive and quick clays, which lose a considerable part of their strength when shaken, commonly fail during an earthquake and become rapid earthflows. Extensive studies were made of the sensitivity of the Bootlegger Cove Clay at Anchorage because of the marked loss of shear strength and dramatic failures of the deposits during the Alaska earthquake of 1964. If similar sensitive clays are present in some places in southeastern Alaska, they most likely are in some of the emergent fine-grained marine deposits; supporting data to confirm their presence, however, are largely lacking. Records of some 50 major earthquakes show that in at least half of the instances water and sediment have been ejected from surficial deposits Water-sediment ejection and associated subsidence and ground fracturing commonly cause extensive damage to the works of man. Ejecta may fill basements and other low-lying parts of buildings. Agricultural land can be covered with a blanket of infertile soils, and small ponds can be filled or made shallow. In southeastern Alaska these phenomena are most likely to occur on valley floors, deltas, tidal flats, alluvial fans, swamps, and lakeshores. Earthquake-induced sliding on land generally is confined to steep slopes but may take place in fine-grained deposits on moderately to nearly flat surfaces if the deposits are subject to liquefaction. A large rockslide triggered by the Lituya Bay, Alaska, earthquake of July 10, 1958, generated a wave that surged up the opposite wall of the inlet to a record height of 1,740 feet. During the Hebgen Lake, Montana, earthquake of August 17, 1959, a spectacular rockslide plunged into the Madison River canyon, buried 28 people, dammed the river, and created a large lake. Earthquake-records are replete with accounts of sliding of surficial deposits during moderate to large earthquakes. Most or all of the general factors that favor subaerial landsliding are present in southeastern Alaska. Earthquake-induced subaqueous slides can produce adverse effects both nearshore and some distance offshore. Nearshore sliding may progress shoreward and destroy harbor facilities and other structures, commonly with substantial loss of life. Disastrous large submarine slides occurred along the fronts of deltas in Seward and Valdez during the Alaska earthquake of 1964. In similar fashion, the largest submarine slides in southeastern Alaska likely will be triggered along the larger delta fronts. Sliding farther offshore can constitute a threat to navigation because of changes in water depths. Also underwater sliding can break communication cables. Glaciers were not greatly affected by the Alaska earthquake of 1964 despite the fact that about 20 percent of the area that underwent strong shaking is covered by ice. In contrast, the cataclysmic avalanche of ice and rock that fell from a high glacier-covered peak in Peru during the earthquake of May 31, 1970, produced devastating effects downvalley on man and his works in the form of mudflows. Most towns in southeastern Alaska are sufficiently distant from glaciers so as not be to directly affected. Both the Alaska earthquake of 1964 and the Hebgen Lake, Montana, earthquake of 1959 significantly affected ground- and surface-water regimens. Water levels in some wells declined whereas in others flow increased. Some springs discharged at a rate three times as much as normal; flow of others decreased or stopped. Discharge of many streams increased markedly. Most or all of the effects described above could occur in parts of southeastern Alaska during future large earthquakes. Tsunamis, seiches, and other abnormal water waves associated with large earthquakes commonly cause vast property damage and heavy loss of life. Tsunami effects can be devastating to coastal areas as far as many thousands of miles from their generation source. Seiche effects generally are confined to inland bodies of water or to relatively enclosed coastal bodies of water. Abnormal waves generated by submarine sliding or by subaerial sliding into water generally produce only local effects but may be highly devastating. Tsunami waves resulting from the Chilean earthquake of 1960 inflicted extensive damage and loss of life on coastal communities throughout a large part of southern Chile, and significant runups and damage were recorded in many places throughout the Pacific Ocean area. The tsunami waves generated by the Alaska earthquake of 1964 struck with devastating force along a broad stretch of the Alaska coast and produced heavy property damage and loss of life as far away as Crescent City, Calif. Seiche waves generated by that earthquake reached runup heights of 20-30 feet on some lakes in Alaska, and water-level fluctuations were recorded on streams, reservoirs, lakes, and swimming pools in States bordering the Gulf of Mexico. Waves generated by submarine sliding struck violently at a number of places during or immediately after the quake and were the major cause of loss of life and damage to property. Slide-generated waves probably would have a higher destructive potential in southeastern Alaska than either tsunami waves or seiche waves because of their possibly higher local runups and because they can hit the shores almost without warning during or immediately after an earthquake. Nonearthquake-related geologic hazards, although generally far less dramatic than those related to earthquakes, tend to occur so much more frequently or persistently that their aggregate effects can be significant. Three kinds of geologic hazards of this type are discussed: (1) nonearthquake-induced landsliding and subaqueous sliding, (2) flooding, and (3) land uplift. The potential for nonearthquake-triggered landsliding in southeastern Alaska ranges widely from place to place. Past sliding generally furnishes the clue in the prediction of where and in what materials future sliding will occur. Fast-moving rockslides, debris slides, and mudflows can be expected to occur from time to time on steep slopes and be highly destructive to highways, power plants, pipelines, buildings, and other facilities located on a slope or at its base. Present slow downslope movement of talus can be expected to continue at the same general rate unless conditions are changed by man or there are climatic changes. Snow and debris avalanches can be especially hazardous during winter months. Long-inactive landslides may be triggered into renewed activity or new slides may be created by man-induced modifications. Accelerated slope erosion and debris flows may follow large-scale clearing and cutting of timber. Subaqueous sliding can be expected to occur periodically along fronts of deltas and on other oversteepened underwater slopes. Floods have been common in parts of southeastern Alaska because of heavy precipitation and rapid runoff from steep slopes with resulting heavy damage to roads and other facilities. Continued damage can be expected in the future unless more remedial measures are taken. Current uplift of land in southeastern Alaska, although probably not affecting man significantly in a short period of time, may have some adverse long-term effects. These long-term effects should be borne in mind when facilities such as docks and boat harbors are constructed on or near the shore, where there is a critical relation between height of land and water.

Alaska↗

Preliminary geologic map of the Simi 7.5' quadrangle, southern California: A digital database

The Simi Quadrangle covers an area of about 62 square miles in southern Ventura County. The Santa Clara River Valley occupies the northwestern corner of the quadrangle. Mountainous terrain of South Mountain and Oak Ridge characterizes the northern and central area. Elevation within the quadrangle ranges from about 250 feet along the arroyo bottoms to over 2200 feet. Steep, highly dissected slopes form much of the boundary of the area. In the southeast, Little Simi Valley, drained by Arroyo Simi/Arroyo Las Posas, separates the southern flank of Oak Ridge from the Las Posas Hills. The Las Posas upland area, a broad elevated region that slopes gently to the south, separates the South Mountain-Oak Ridge highlands from the Las Posas-Camarillo Hills between Little Simi Valley on the east and the Oxnard Plain on the west. This relatively low-lying area is also referred to as the Las Posas Valley. Numerous north-south-trending drainages cut South Mountain and Oak Ridge creating steep narrow canyons on north-facing slopes and wide flat-bottomed canyons with incised streams on south-facing slopes. A network of residential streets and ranch and oilfield roads that traverse the area from U.S. Highway 101 and State Highways 118, 23, and 126 provides access to the area. Current land use includes citrus and avocado orchards, oil well drilling and production, sand and gravel quarries, decorative-rock quarries, cattle grazing, suburban residential development, and golf courses. The oldest geologic unit mapped in the Simi Quadrangle is the upper Eocene to lower Miocene Sespe Formation. The Sespe Formation consists of alluvial fan and floodplain deposits of interbedded pebble-cobble conglomerate, massive to thick-bedded sandstone, and thin-bedded siltstone and claystone. In the northern part of the map area, Sespe Formation is overlain by and interfingers with the upper Oligocene to lower Miocene Vaqueros Formation that is composed of transitional and marine sandstone, siltstone, and claystone with local sandy coquina beds. In the Las Posas Hills, Sespe Formation is unconformably overlain by marine sandstones of the middle Miocene Topanga Group that are interlayered with and intruded by basalt flows, breccia, and diabase dikes of the Conejo Volcanics. Deep-marine strata of the upper Miocene Modelo Formation cover the Vaqueros Formation and Topanga Group along the crests and southern flanks of South Mountain and Oak Ridge. They also occur as isolated outcrops in the Las Posas Hills. Locally, Modelo Formation consists of interbedded diatomaceous shale, claystone, mudstone, and siltstone with minor sandstone, limestone, chert, and tuff beds. The most widely exposed rock units in the area are the Plio-Pleistocene marine and non-marine Pico and Saugus Formations that crop out on the southern flank of South Mountain-Oak Ridge. Locally, the Pico Formation consists of marine siltstone and silty shale with minor sandstone and pebbly sandstone. The Saugus Formation overlies and interfingers with the Pico Formation and is composed of interbedded shallow-marine to brackish water sandstone, siltstone, pebble-to-cobble conglomerate, and coquina beds that grade laterally and vertically into non-marine sandstone, siltstone, and conglomerate. A local member of the Saugus Formation is exposed in the southwest corner of the map area. It is predominantly a volcanic breccia conglomerate that resembles the Conejo Volcanics breccia, but is believed to represent remnants of landslide debris shed from the Conejo Volcanics into a local trough during Saugus time. Quaternary surficial deposits cover the floor and margins of the Little Simi Valley, Santa Clara River Valley in the north, and Arroyo Las Posas in the south, and extend up into the larger canyons that drain South Mountain and Oak Ridge. Extensive surficial deposits are also present in the Las Posas upland area in the southwest. These upper Pleistocene to Holocene sediments consist of older and younger alluvial fan and valley deposits, colluvium, active alluvial fans, and active stream deposits. Pleistocene- to Holocene-age landslide deposits are widespread throughout the Simi Quadrangle, especially in the finer grained Tertiary sedimentary units where bedding planes are dip slopes. In addition, massive slumps are present in the Sespe and Vaqueros Formations on anti-dip slopes. Seismic and well data from the San Fernando Valley (SFV) document evolution of that region from mid-Miocene rifting to north-south contraction. Formations in the western SFV subsurface (Cretaceous to Paleogene strata, and Miocene Topanga and Modelo Formations) trace southward to outcrops in the Santa Monica Mountains that constrain faulting along the valley's south basin edge. Cretaceous strata in the Simi Uplift to the west are over 2 km higher than equivalent strata beneath the western SFV across a boundary marked by the Chatsworth Reservoir fault, and Neogene thinning and offlap. The Simi fault, located at the eastern end of the Simi-Santa Rosa fault system, bounds the northern margins of the Simi and Tierra Rejada Valleys. West of Simi Valley, the Simi fault has placed Miocene Conejo Volcanics over Plio-Pleistocene Saugus Formation rocks. The 15.5 ± 0.8 m.y.a. base of the Conejo Volcanics, identified in oil well logs, is inferred to have a dip-slip separation of about 425 to 550 m, suggesting a low long-term slip rate of about 0.03 mm/yr. However, substantial late Quaternary offset is suggested by the presence of more than 150 m of Pleistocene and younger alluvium that fills the east-west trending, down-dropped bedrock trough beneath western Simi Valley. In addition, trenching within faulted colluvial deposits in Tierra Rejada Valley has revealed evidence of multiple shears within Holocene (?) deposits.

California↗

Catalog of the historically active volcanoes of Alaska

Alaska hosts within its borders over 80 major volcanic centers that have erupted during Holocene time (< 10,000 years). At least 29 of these volcanic centers (table 1) had historical eruptions and 12 additional volcanic centers may have had historical eruptions. Historical in Alaska generally means the period since 1760 when explorers, travelers, and inhabitants kept written records. These 41 volcanic centers have been the source for >265 eruptions reported from Alaska volcanoes. With the exception of Wrangell volcano, all the centers are in, or near, the Aleutian volcanic arc, which extends 2500 km from Hayes volcano 145 km west of Anchorage in the Alaska-Aleutian Range to Buldir Island in the western Aleutian Islands (fig. 1). The volcanic arc, a subduction-related feature associated with underthrusting of the Pacific plate beneath the North American plate is divided between oceanic island arc and continental margin segments, the boundary occurring at about 165° W longitude (fig. 1). An additional 7 volcanic centers in the Aleutian arc (table 2; fig. 1 A) have active fumarole fields but no reported historical eruptions. This report discusses the location, physiography and structure, eruptive history, and geology of those volcanoes in Alaska that have experienced one or more eruptions that have been recorded in the written history (i.e., in historical time). It is part of the group of catalogs entitled Catalogue of Active Volcanoes of the World published beginning in 1951 under the auspices of the International Association of Volcanology and Chemistry of the Earth's Interior (IAVCEI). A knowledge of the information contained in such catalogs aids in understanding the type and scale of activity that might be expected during a particular eruption, the hazards the eruption may pose, and even the prediction of eruptions. The catalog will thus be of value not only to the inhabitants of Alaska but to government agencies concerned with emergency response, air traffic operations, and weather, as well as to industry and scientists. The combination of the hazard posed by volcanic ash to jet aircraft and the heavy use of international air routes located parallel to, and on either side of, the Aleutian volcanic arc means that even remote volcanoes in Alaska now pose significant hazards to life and property. Although this report is concerned with historical eruptions from Alaskan volcanoes, other volcanoes in Alaska have erupted in the past 10,000 years and might therefore be expected to erupt again. Several Holocene volcanic centers in the Aleutian arc have no reported historical activity. Elsewhere in Alaska the Bering Sea basalt fields cover large areas of the Yukon Delta, Seward Peninsula, and several of the islands of the Bering Sea. Holocene centers also occur in the Wrangell Mountains and in isolated occurrences in the interior and southeastern Alaska. Eruptions from these centers have occurred within the past several hundred years but none were transcribed in the written record. Moodie and others (1992), however, report oral traditions among the Northern Athapaskan Indians of the southwestern Yukon Territory that may record the second and younger deposition of the White River Ash circa A.D. 720. This lobe of the White River Ash was deposited during the paroxysmal eruption of Churchill volcano in the Wrangell Mountains of eastcentral Alaska (McGimsey and others, 1992; Richter and others, 1995).

Alaska↗

Seismic stability of the Duwamish River Delta, Seattle, Washington

The delta front of the Duwamish River valley near Elliott Bay and Harbor Island is founded on young Holocene deposits shaped by sea-level rise, episodic volcanism, and seismicity. These river-mouth deposits are highly susceptible to seismic soil liquefaction and are potentially prone to submarine landsliding and disintegrative flow failure. A highly developed commercial-industrial corridor, extending from the City of Kent to the Elliott Bay/Harbor Island marine terminal facilities, is founded on the young Holocene deposits of the Duwamish River valley. The deposits of this Holocene delta have been shaped not only by relative sea-level rise but also by episodic volcanism and seismicity. Ground-penetrating radar (GPR), cores, in situ testing, and outcrops are being used to examine the delta stratigraphy and to infer how these deposits will respond to future volcanic eruptions and earthquakes in the region. A geotechnical investigation of these river-mouth deposits indicates high initial liquefaction susceptibility during earthquakes, and possibly the potential for unlimited-strain disintegrative flow failure of the delta front.

Washington↗

Field and laboratory data From an earthquake history study of scarps of the Lake Creek-Boundary Creek fault between the Elwha River and Siebert Creek, Clallam County, Washington

Fault scarps recently discovered on Airborne Laser Swath Mapping (ALSM; also known as LiDAR) imagery show Holocene movement on the Lake Creek–Boundary Creek fault on the north flank of the Olympic Mountains of northwestern Washington State. Such recent movement suggests the fault is a potential source of large earthquakes. As part of the effort to assess seismic hazard in the Puget Sound region, we map scarps on ALSM imagery and show primary field and laboratory data from backhoe trenches across scarps that are being used to develop a latest Pleistocene and Holocene history of large earthquakes on the fault. Although some scarp segments 0.5–2 km long along the fault are remarkably straight and distinct on shaded ASLM imagery, most scarps displace the ground surface <1 m, and, therefore, are difficult to locate in dense brush and forest. We are confident of a surface-faulting or folding origin and a latest Pleistocene to Holocene age only for scarps between Lake Aldwell and the easternmost fork of Siebert Creek, a distance of 22 km. Stratigraphy in five trenches at four sites help determine the history of surface-deforming earthquakes since glacier recession and alluvial deposition 11–17 ka. Although the trend and plunge of indicators of fault slip were measured only in the weathered basalt exposed in one trench, upward-splaying fault patterns and inconsistent displacement of successive beds along faults in three of the five trenches suggest significant lateral as well as vertical slip during the surface-faulting or folding earthquakes that produced the scarps. Radiocarbon ages on fragments of wood charcoal from two wedges of scarp-derived colluvium in a graben-fault trench suggest two surface-faulting earthquakes between 2,000 and 700 years ago. The three youngest of nine radiocarbon ages on charcoal fragments from probable scarp-derived colluvum in a fold-scarp trench 1.2 km to the west suggest a possible earlier surface-faulting earthquake less than 5,000 years ago.

Washington↗

Geologic map of the Craters of the Moon 30' x 60' Quadrangle, Idaho

The Craters of the Moon 30 x 60 minute quadrangle shows the geology of the northern two-thirds of the Craters of the Moon (COM) lava field and volcanic structures of the northern and central parts of the Great Rift volcanic rift zone. The COM lava field is the largest, predominantly Holocene lava field in the conterminous United States. The northwest corner of the map shows older sedimentary, intrusive, and volcanic rocks that range in age from Ordovician to Miocene. These rocks provide evidence of compressional fold and thrust events of the Antler and Sevier orogenies. Compression was followed by voluminous volcanism represented by the Challis Volcanic Group. Basin-and-Range faulting followed in Neogene time. The COM lava field covers about 1,600 square kilometers and contains about 30 cubic kilometers of lava flows and associated vent deposits. Stratigraphic relationships, paleomagnetic studies, and radiocarbon ages indicate that the field formed during eight eruptive periods designated as H, the oldest, to A, the youngest. Each eruptive period was several hundred years or less in duration and separated from other eruptive periods by non-eruptive recurrence intervals of several hundred to about 3,000 years. The first eruptive period began about 15,000 carbon-14 years ago and the latest one ended about 2,100 carbon-14 years ago. All available field, paleomagnetic, radiocarbon, and argon-40/argon-39 data are incorporated in this map and they quantitatively refine the volcanic and paleomagnetic history of the pre-Holocene lava fields and the COM lava field. In a sense, these data determine the 'pulse rate' for Pleistocene and Holocene basaltic volcanism in the area of this map. Twenty-three new argon-40/argon-39 geochronologic data reveal a fairly complete and continuous record of basaltic volcanism in the Craters of the Moon 30 x 60 minute quadrangle for the last 500 ka. The ages cluster into age groupings at ~30 ka, 50-70 ka, 100-125 ka, 260-290 ka, 320-340 ka, and 475 ka. There are apparent periods of ~30 to 60 ka duration when little or no volcanic activity took place between groups. Magnetic polarity and remanent inclination and declination directions for most lava flows in the quadrangle have normal magnetic polarity; they were emplaced during the Brunhes Normal Polarity Chron and are younger than 780,000 years. Directions of remanent magnetization and the new argon-40/argon-39 ages were used to correlate and approximately date lava flows and lava fields for this map.

Scientific Investigations Map↗

Geologic map of the north half of the Lake Walcott 30'×60' quadrangle, Idaho

The geologic map of the northern half of the Lake Walcott 30ʹ×60ʹ quadrangle shows the volcanic geology of the southern part of the Craters of the Moon lava field, the complex geologic features of the Holocene Kings Bowl and Wapi lava fields, and the southern part of the Great Rift volcanic rift zone. The long extent and distribution of skylights in lava-tube systems of the Horse Butte and Wapi Park lava fields are depicted on this map. 40 Ar/ 39 Ar and K/Ar age determinations give detail to the Holocene, late Pleistocene, and late middle Pleistocene volcanic lava fields in this quadrangle. Most of the younger basalt eruptions (less than 150 thousand years [ka]) have occurred along the Great Rift volcanic rift zone, but two of the younger lava fields are located in the western part of the quadrangle. Kimama Butte, a shield volcano, is 87±11 ka, and Shale Butte is dated at 11±6 ka. Paleomagnetic studies have shown that the Horse Butte-Inferno Chasm eruptive fissure system has at least five paleomagnetic-correlative lava fields, the Claasen vent complex consists of at least seven correlative lava fields, and the Streifling-Flat Top vent complex includes at least four correlative lava fields. This map provides geologic, geochronologic, and paleomagnetic data for Holocene lava fields along the southern part of the Great Rift, and for late Pleistocene and late middle Pleistocene lava fields in the central and western parts of the quadrangle. These data can contribute to wise management and preservation of the Craters of the Moon National Monument and for broad-scale understanding of the basaltic-volcanic evolution of the eastern Snake River Plain.

Idaho↗

Geologic map of the northern Harrat Rahat volcanic field, Kingdom of Saudi Arabia

Harrat Rahat, in the west-central part of the Kingdom of Saudi Arabia, is the largest of 15 Cenozoic harrats (Arabic for “volcanic field”) distributed on the Arabian plate. It extends more than 300 km north-south and 50 to 75 km east-west, and it covers an area of approximately 20,000 km 2 , has a volume of approximately 2,000 km 3 , and encompasses more than 900 observable vents. Volcanism commenced around 10 Ma and has continued into historic time, the most recent eruption occurring in 1256 C.E. Volcanic products are dominated by alkali basalt and hawaiite lava flows, with subordinate mugearite lava flows, as well as benmoreite and trachyte lava flows, domes, and pyroclastic flows. This geologic map distinguishes 239 eruptive units that cover an area of 3,340 km 2 in northern Harrat Rahat and the adjacent city of Al-Madinah. Results are presented as a geologic map of the study area at 1:75,000 scale and of smaller regions of particular interest at 1:25,000 scale, along with interpretive text. Most units are basaltic lava flows that erupted from the broadly north-northwest-trending main vent axis that constructed the topographic crest of the volcanic field. This 300- to 400-m-high vent axis, which has a width of 6 to 10 km, lies in the eastern one-third of northern Harrat Rahat. Basalt and hawaiite lava flows can extend as far as 27 km from their vents, but most are 10 to 15 km long. Evolved products such as mugearites, benmoreites, and trachytes are less extensive; the trachytic pyroclastic flows extend as far as 9 km from their source vents, although most only reach 4 to 6 km. Vents of the evolved products are restricted to the main vent axis or its flanks. No volcanic rocks older than 1.2 Ma are exposed in the map area, and about 90 percent of the exposed volcanic rocks erupted during the past 570 thousand years. As depicted on the geologic maps, eruption ages and field relations define 12 eruptive stages for northern Harrat Rahat for the past 1.2 million years. Other important geochronological findings include (1) several late Pleistocene lava flows near Al-Madinah, which previously were interpreted as Holocene from archeological evidence; (2) the eruption age of a cluster of cinder cones and small lava flows in the western outskirts of Al-Madinah (previously ascribed to an eruption in 641 C.E.) is actually 13.3±1.9 ka, close to the Pleistocene-Holocene boundary; and (3) only two Holocene eruptions have been identified in the map area, those of the historically described basalt of Al Labah in 1256 C.E. and the dome and pyroclastic flows of the trachyte of Um Rgaibah at 4.2±5.2 ka.

Scientific Investigations Map↗

Surficial geology and Quaternary tectonics of the Madison Valley and fault zone, Madison, Gallatin, and Beaverhead Counties, southwest Montana

The north-northwest-striking Madison fault is approximately 95 kilometers in length, lying at the confluence of the northeastern Basin and Range province and the Yellowstone tectonic parabola. The fault zone consists primarily of west-dipping normal faults that have east-dipping antithetic faults, which create the Madison Valley graben and several northeast-trending intrabasin faults. The Madison fault and associated sections discussed herein refer to the main west-dipping, range-bounding fault along the eastern side of the valley. Detailed geologic mapping (1:12,000 scale) of the entire fault zone and fault scarp profiling (total of 102 profiles) of the Madison fault reveal greater late Quaternary paleoseismic activity towards the south, including at least three paleoevents along the southern part of the fault that postdate Pinedale glaciation. Early to middle Holocene alluvial fans have vertical surface offsets that average between 2.0 and 3.0 meters and define the characteristic single-event surface offset. Pinedale lateral moraines have vertical surface offsets as great as 12.0 meters. Late Pleistocene to Holocene multiple-event fault scarps show little evidence of beveling, suggesting short seismic recurrence intervals and potential late Pleistocene and Holocene temporal clustering. Long-term average tectonic activity rates indicate slip rates ranging from 0.18–0.6 millimeters per year. Based on a comparison of fault-scarp height versus maximum slope angle of known regression lines developed from other paleoseismic investigations, the most recent event ranges from 5–1 ka. The northern section of the fault zone is defined by multiple normal faults, which detached the hanging walls of Laramide thrust faults within the Paleozoic and Mesozoic strata. This resulted in the partitioning of extension along multiple preexisting structures and less displacement along individual normal fault strands. Structural controls on lateral propagation of individual paleoevents involve the position of lateral ramps along preexisting Laramide contractional faults. This resulted in greater displacement within the larger basement-cored structures along the southern section of the fault zone, where extension is accommodated by one inferred principal basement-involved normal fault. Inferred east-northeast trending, intrabasin, normal faults within the southern half of the fault zone have no late Pleistocene displacement.

Idaho, Montana, Wyoming↗

Map of topographic lineaments interpreted as recent surface ruptures along the Bennett Valley and Southern Maacama Fault Zones, Sonoma County, California

This study documents lidar-illuminated topographic lineaments interpreted as evidence of recent surface fault ruptures and surface ruptures related to distributed deformation along the Bennett Valley Fault Zone and the southernmost Maacama Fault Zone in the northern San Francisco Bay area (fig. 1, on map sheet). Together, these fault zones form a structural connection across a right stepover between the main Maacama Fault Zone and the Rodgers Creek Fault, overlapping principal strands of the San Andreas plate boundary system north of San Francisco Bay (figs. 1, 2, on map sheet) that accommodate about a quarter of the ~40 millimeters per year of regional dextral tectonic slip (Parsons and others, 2013). Although much of the Bennett Valley Fault Zone is seismically active (McLaughlin and others, 2012; Sowers and others, 2010), only the north end of the fault zone (the “Spring Valley strand,” indicated on southern map panel of map sheet and labeled “D” in figure 2, on map sheet) was previously known to have youthful geomorphic expression and to displace Holocene deposits (McLaughlin and others, 2008; Sowers and others, 2010; Sowers and others, 2016). Holocene activity along the length of the Bennett Valley Fault Zone and southernmost Maacama Fault Zone, with possible implications for rupture propagation, continuity, and slip transfer, had not been identified. However, fault splays that project eastward from the Rodgers Creek Fault have been mapped and hypothesized to accommodate slip transfer to the Bennett Valley Fault Zone (Hecker and Randolph Loar, 2018). The generally subtle and distributed nature of surface-rupture evidence along the Bennett Valley and southernmost Maacama Fault Zones, and extensive vegetation cover, had left recent faulting previously unmapped along most of the zone. The map presented here represents a new compilation of inferred surface-rupture features detected using high-resolution topographic lidar data from an airborne lidar survey of Sonoma County, California (OpenTopography, 2014). These data, which enable subtle topographic features to be discerned, indicate that recent (likely Holocene) surface ruptures extend throughout the Bennett Valley and southernmost Maacama Fault Zones. The rupture-lineament map was created using a GIS (geographic information system) framework and is included herein as an image map at a scale of 1:36,000 and as digital datasets (included as supplemental information to this report). The mapping is intended to lay the groundwork for future studies designed to better characterize how plate-boundary slip is accommodated on this important and complex system of faults.

California↗

An index of ecological integrity for the Mississippi alluvial plain ecoregion: index development and relations to selected landscape variables

Macroinvertebrate community, fish community, water-quality, and habitat data collected from 36 sites in the Mississippi Alluvial Plain Ecoregion during 1996-98 by the U.S. Geological Survey were considered for a multimetric test of ecological integrity. Test metrics were correlated to site scores of a Detrended Correspondence Analysis of the fish community (the biological community that was the most statistically significant for indicating ecological conditions in the ecoregion) and six metrics--four fish metrics, one chemical metric (total ammonia plus organic nitrogen) and one physical metric (turbidity)--having the highest correlations were selected for the index. Index results indicate that sites in the northern half of the study unit (in Arkansas and Missouri) were less degraded than sites in the southern half of the study unit (in Louisiana and Mississippi). Of 148 landscape variables evaluated, the percentage of Holocene deposits and cotton insecticide rates had the highest correlations to index of ecological integrity results. sites having the highest (best) index scores had the lowest percentages of Holocene deposits and the lowest cotton insecticide use rates, indicating that factors relating to the amount of Holocene deposits and cotton insecticide use rates partially explain differences in ecological conditions throughout the Mississippi Alluvial Plain Ecoregion.

Arkansas, Kentucky, Louisiana, Mississippi, Missou↗

Geology and quality of water in the Modesto-Merced area, San Joaquin Valley, California, with a brief section on hydrology

The Modesto-Merced area includes about 1,800 square miles on the northeast side of the San Joaquin Valley. The physiographic units in the area are (1) Sierra Nevada, (2) dissected uplands, (3) low alluvial plains and fans, (4) river flood plains and channels, and (5) overflow lands. Geologic units consist of consolidated rocks and unconsolidated deposits. The dominant structure of the geologic units is that of a homocline, which reflects the southwestward-tilted fault block of the Sierra Nevada. The consolidated rocks include: (1) basement complex (pre-Tertiary), (2) marine sandstone and shale (Cretaceous), (3) Ione Formation and other sedimentary rocks (Eocene), (4) Valley Springs Formation (late? Miocene), and (5) the Mehrten Formation (Miocene and Pliocene). In the eastern part of the area, the consolidated rocks generally yield small quantities of water to wells except for the Mehrten Formation which is an important aquifer. The unconsolidated deposits include: (1) continental deposits (Pliocene and Pleistocene?), (2) lacustrine and marsh deposits (Pleistocene), (3) older alluvium (Pleistocene and Holocene?), (4) younger alluvium (Holocene), and (5) flood-basin deposits (Holocene). The continental deposits and older alluvium are the main water-yielding units in the unconsolidated deposits. The lacustrine and marsh deposits (E-clay) and the flood-basin deposits yield little water to wells, and the younger alluvium in most places probably yields only moderate quantities of water to wells. There are three ground-water bodies in the Modesto-Merced area: (1) the unconfined water body, (2) the confined water body, and (3) the water body in consolidated rocks. The unconfined water body occurs in the unconsolidated deposits above and east of the E-clay, except in the western and southern parts of the area where clay lenses occur and semiconfined conditions exist. The confined water body occurs in the unconsolidated deposits below the E-clay and extends downward to the base of fresh water. The water body in consolidated rocks occurs under both perched and confined conditions. Ground-water movement in the unconfined and the confined water bodies is generally westward toward the valley trough. In the unconfined water body, the water also moves toward the major rivers, and toward pumping depressions at Modesto and near El Nido. Because of the higher head in the unconfined water body, water slowly moves from it through the E-clay to the underlying confined water body. Surface water is used extensively for irrigation in most of the study area; consequently, shallow water in the unconfined water body has to be controlled by pumping. Nevertheless, water levels near Modesto declined about 6 feet between 1958 and 1962, for the most part during the dry years 1959-61. Near El Nido, water levels declined about 70 feet between 1942 and 1967. In the confined water body, water levels were high in the winter and spring and low in the summer and fall, reflecting irrigation practices. In the southwestern part of the area, water levels in the confined water body declined about 14 feet from 1962 to 1968 and rose slightly after 1968, whereas in the northwestern part water levels remained fairly constant. Water from the upper reaches of the Stanislaus, Tuolumne, Merced, and Chowchilla Rivers is calcium bicarbonate in chemical type and is of excellent quality. As those rivers cross the valley floor, their water quality is generally degraded by return flows from irrigated land, and in the Tuolumne River by saline water discharged from abandoned gas wells. Average dissolvedsolids content in the major rivers, as indicated by chemical analyses, do not exceed 400 mg/1 (milligrams per liter) except in the San Joaquin River where average dissolved-solids content has not exceeded 1,050 mg/l. Water from minor streams in the area is bicarbonate in chemical type with dissolved-solids content ranging from 56 mg/1 to about 350 mg/l. Although chloride-type fresh ground water occurs in the unconfined and confined water bodies and in the water body in consolidated rocks, most of the fresh ground water is a bicarbonate type that has a dissolved-solids content of less than 500 mg/l. Water having dissolved solids in excess of about 2,000 mg/1 is considered to be saline. Saline water extends below the base of fresh water to the basement complex and, except in the extreme eastern part, underlies most of the study area. Saline water also occurs as lenses above the base of fresh water in the unconfined and confined water bodies.

California↗

Geohydrologic reconnaissance of the Crofton Unit, northeastern Nebraska

This report describes the results of a reconnaissance of the sources, location, quantity, and quality of ground water in the proposed irrigation project area called the "Crofton Unit," an area of approximately 350 square miles in northeastern Nebraska. Precambrian Sioux Quartzite underlies the entire study area. Throughout most of the area, the Sioux is overlain by undifferentiated Paleozoic rocks. The Sioux and, where present, the Paleozoic rocks are overlain unconformably by Cretaceous formations, including both the Dakota and Niobrara Formations. Surface strata are Pleistocene and Holocene deposits of sand and gravel and Pleistocene deposits of glacial till. The Dakota and Niobrara Formations are each capable of yielding sufficient water for irrigation in parts of the area. Pleistocene and Holocene deposits yield adequate water for irrigation in the southern part of the area where saturated thicknesses range from 25 to 176 feet and along the flood plain of the Missouri River where saturated thicknesses range from 40 to 101 feet. Reported yields of 13 Pleistocene irrigation wells sampled for water-quality analyses averaged 857 gallons per minute. Water from wells developed in the Dakota Formation is a calcium sulfate type. Use of this water for successful long-term irrigation would necessitate burdensome salinity-control measures because of the salinity of the water. Water from wells developed in the Niobrara Formation is a calcium bicarbonate type and is suitable for irrigation of most crops. Water from wells developed in both Niobrara and Pleistocene deposits generally is a calcium sulfate type and may be used for irrigation of salt-tolerant crops with special measures for salinity control. The quality of water from Pleistocene and Holocene deposits is highly variable. In both the southern part of the area and along the Missouri River, water from these deposits is a calcium bicarbonate type and generally suitable for irrigation.

Nebraska↗

Hydrogeology and predevelopment flow in the Texas Gulf Coast aquifer systems

A multilayered ground-water flow system exists in the Coastal Plain sediments of Texas. The Tertiary and Quaternary clastic deposits have an area! extent of 128,000 square miles onshore and in the Gulf of Mexico. Two distinct aquifer systems are recognized for the sediments, which range in thickness from a few feet to more than 12,000 feet. The older system the Texas coastal uplands aquifer system consists of four aquifers and two confining units in the Wilcox and Claiborne Groups. It is bounded from below by the practically impermeable Midway confining unit or by the top of the geopressured zone. It is bounded from above by the poorly permeable Vicksburg-Jackson confining unit, which separates it from the younger coastal lowlands aquifer system. The coastal lowlands aquifer system consists of five permeable zones and two confining units that range in age from Oligocene to Holocene. The hydrogeologic units of both systems are exposed in bands that parallel the coastline. The units dip and thicken toward the Gulf. Quality of water in the aquifer systems varies greatly, with dissolved solids ranging from a few hundred to more than 200,000 milligrams per liter. A three-dimensional, variable-density digital model was developed to simulate predevelopment flow in the aquifer systems, for which steady-state conditions were assumed. Horizontal hydraulic conductivities of the aquifers and permeable zones in the calibrated model range from 15 feet per day for the middle Wilcox aquifer, to 170 feet per day for the Holocene-upper Pleistocene aquifer. Vertical hydraulic conductivities range from 1 x 10 -5 foot per day for the Vicksburg-Jackson confining unit, to 1 x 10 -2 foot per day for four of the aquifers and permeable zones. The simulated values of transmissivity and leakance are functions of the percent of sand that is present in each model grid block. There is a large range in precipitation across the study area, from about 21 inches per year in the west to about 56 inches per year in the east. Eastward from a line through Corpus Christi and San Antonio, average annual precipitation ranges from about 30 to about 56 inches. A few inches per year reaches the saturated zone in topographically high areas and is discharged in low areas as evapotranspiration, seepage, springflow, and stream base flow. A smaller amount of water flows through the aquifers and permeable zones downdip from the outcrop areas. This flow results in upward or downward leakage into adjacent hydrogeologic units, but is generally upward into overlying units. Westward from the line through Corpus Christi and San Antonio, average annual precipitation ranges from about 30 to about 21 inches. The general pattern of flow in the aquifers and permeable zones is similar to that in the east, but rates of flow are somewhat smaller. In contrast to the east, ground-water discharge in the west is generally not visible. Evapotranspiration is the main mechanism for ground-water discharge, with most ground water being discharged through evapotranspiration by phreatophytes. Simulated discharge and recharge rates in the combined outcrop areas of all units do not exceed 6 inches per year. The large rates occur in small, local topographically low and high areas. The average discharge rate simulated in the outcrops of the units is 0.45 inch per year. The recharge area is considerably smaller than the discharge area, and the average recharge rate over this smaller area is 0.74 inch per year. Total simulated recharge in the outcrop areas is 269 million cubic feet per day, which is offset by an equal amount of discharge in the outcrop areas. The smallest rates of leakage are across the Vicksburg-Jackson confining unit, with downward and upward rates of less than one million cubic feet per day. The greatest rate of leakage is 47 million cubic feet per day upward into the Holocene-upper Pleistocene permeable zone.

Texas↗

Late-Glacial paleoecology of the Middle Susitna Valley, Alaska: Environmental context for human dispersal

We present here the results of multi-proxy analyses (sediment geochemistry, diatoms, and pollen) from sediment cores collected at four lakes in the middle Susitna Valley, Alaska. These lakes form a transect from the tundra to the boreal forest. The retrieved cores span from ~12,000 cal yr BP to the present, with age control provided by radiometric dates and tephra deposits, some of which are newly identified. Results indicate that deglaciation occurred before 12,000 cal yr BP and that by that time, the lakes were deep, productive, and surrounded by shrub tundra. The lake with the highest sampling resolution indicates a brief climatic reversal ~11,500 cal yr BP with decreased diatom-inferred lake level and lowered lake productivity, and reduced shrub presence. During the early to middle Holocene, all of the sedimentary records provide evidence of climatic amelioration with tree expansion and productive lakes. A middle to late Holocene climatic deterioration with reduced trees and a shallower, less productive lake is also indicated. In addition, the prominent Watana tephra at ~4,000 cal yr BP likely reduced lake productivity and affected the vegetation. Even though the region was relatively productive soon after deglaciation, people did not occupy the region until ~11,000 cal yr BP, about 1000 years later, and then only sparsely. By the middle and late Holocene, the region was more densely populated and this shift in human occupancy presumably reflects changes in resource abundance, especially caribou. Whether the Watana ashfall influenced caribou abundance and thus people, is still under investigation, but given the tephra’s effect on vegetation and lake productivity, it seems likely.

Alaska↗

A statistical framework for integrating nonparametric proxy distributions into geological reconstructions of relative sea level

Robust, proxy-based reconstructions of relative sea-level (RSL) change are critical to distinguishing the processes that drive spatial and temporal sea-level variability. The relationships between individual proxies and RSL can be complex and are often poorly represented by traditional methods that assume Gaussian likelihood distributions. We develop a new statistical framework to estimate past RSL change based on nonparametric, empirical modern distributions of proxies in relation to RSL, applying the framework to corals and mangroves as an illustrative example. We validate our model by comparing its skill in reconstructing RSL and rates of change to two previous RSL models using synthetic time-series datasets based on Holocene sea-level data from South Florida. The new framework results in lower bias, better model fit, and greater accuracy and precision than the two previous RSL models. We also perform sensitivity tests using sea-level scenarios based on two periods of interest – meltwater pulses (MWPs) and the Holocene – to analyze the sensitivity of the statistical reconstructions to the quantity and precision of proxy data; we define high-precision indicators, such as mangroves and the reef-crest coral Acropora palmata , with 2 σ vertical uncertainties within ± 3 m and lower-precision indicators, such as Orbicella spp., with 2 σ vertical uncertainties within ± 10 m. For reconstructing rapid rates of change in RSL of up to ∼ 40 m kyr −1 , such as those that may have characterized MWPs during deglacial periods, we find that employing the nonparametric model with 5 to 10 high-precision data points per kiloyear enables us to constrain rates to within ± 3 m kyr −1 (1 σ ). For reconstructing RSL with rates of up to ∼ 15 m kyr −1 , as observed during the Holocene, we conclude that employing the model with 5 to 10 high-precision (or a combination of high- and low-precision) data points per kiloyear enables precise estimates of RSL within &#xB1; &#x223C; "> 2 m (2 σ ) and accurate RSL reconstructions with errors ≲ 0.7 m. Employing the nonparametric model with only lower-precision indicators also produces fairly accurate estimates of RSL with errors ≲ 1.50 m, although with less precision, only constraining RSL to &#xB1; &#x223C; "> 3–4 m (2 σ ). Although the model performs better than previous models in terms of bias, model fit, accuracy, and precision, it is computationally expensive to run because it requires inverting large matrices for every sample. The new model also provides minimal gains over similar models when a large quantity of high-precision data are available. Therefore, we recommend incorporating the nonparametric likelihood distributions when no other information (e.g., reef facies or epibionts indicative of shallow-water environments to refine coral elevational uncertainties) or no high-precision data are available at a location or during a given time period of interest.

Advances in Statistical Climatology, Meteorology a↗