USGS Science⌕ Search

SEARCH · USGS Science

Results for “Ocean Modelling”

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 829 records · Page 46Linked to original sources

Sediment transport measurements

Sediment erosion, transport, and deposition in fluvial systems are complex processes that are treated in detail in other sections of this book. Development of methods suitable for the collection of data that contribute to understanding these processes is a still-evolving science. Sediment and ancillary data are fundamental requirements for the proper management of river systems, including the design of structures, the determination of aspects of stream behavior, ascertaining the probable effect of removing an existing structure, estimation of bulk erosion, transport, and sediment delivery to the oceans, ascertaining the long-term usefulness of reservoirs and other public works, tracking movement of solid-phase contaminants, restoration of degraded or otherwise modified streams, and assistance in the calibration and validation of numerical models. This chapter presents techniques for measuring bed-material properties and suspended and bed-load discharges. Well-established and relatively recent, yet adequately tested, sampling equipment and methodologies, with designs that are guided by sound physical and statistical principles, are described. Where appropriate, the theory behind the development of the equipment and guidelines for its use are presented.

Book chapter↗

Vertical spatial coherence model for a transient signal forward-scattered from the sea surface

The treatment of acoustic energy forward scattered from the sea surface, which is modeled as a random communications scatter channel, is the basis for developing an expression for the time-dependent coherence function across a vertical receiving array. The derivation of this model uses linear filter theory applied to the Fresnel-corrected Kirchhoff approximation in obtaining an equation for the covariance function for the forward-scattered problem. The resulting formulation is used to study the dependence of the covariance on experimental and environmental factors. The modeled coherence functions are then formed for various geometrical and environmental parameters and compared to experimental data.

IEEE Journal of Oceanic Engineering↗

Deep magmatic structures of Hawaiian volcanoes, imaged by three-dimensional gravity models

A simplified three-dimensional model for the island of Hawai'i, based on 3300 gravity measurements, provides new insights on magma pathways within the basaltic volcanoes. Gravity anomalies define dense cumulates and intrusions beneath the summits and known rift zones of every volcano. Linear gravity anomalies project southeast from Kohala and Mauna Kea summits and south from Hualālai and Mauna Loa; these presumably express dense cores of previously unrecognized rift zones lacking surface expression. The gravity-modeled dense cores probably define tholeiitic shield–stage structures of the older volcanoes that are now veneered by late alkalic lavas. The three-dimensional gravity method is valuable for characterizing the magmatic systems of basaltic oceanic volcanoes and for defining structures related to landslide and seismic hazards.

Hawaii↗

Development of a seamless multisource topographic/bathymetric elevation model of Tampa Bay

Many applications of geospatial data in coastal environments require knowledge of the nearshore topography and bathymetry. However, because existing topographic and bathymetric data have been collected independently for different purposes, it has been difficult to use them together at the land/water interface owing to differences in format, projection, resolution, accuracy, and datums. As a first step toward solving the problems of integrating diverse coastal datasets, the U.S. Geological Survey (USGS) and the National Oceanic and Atmospheric Administration (NOAA) are collaborating on a joint demonstration project to merge their data for the Tampa Bay region of Florida. The best available topographic and bathymetric data were extracted from the USGS National Elevation Dataset and the NOAA hydrographic survey database, respectively. Before being merged, the topographic and bathymetric datasets were processed with standard geographic information system tools to place them in a common horizontal reference frame. Also, a key part of the preprocessing was transformation to a common vertical reference through the use of VDatum, a new tool created by NOAA's National Geodetic Survey for vertical datum conversions. The final merged product is a seamless topographic/bathymetric model covering the Tampa Bay region at a grid spacing of 1 arc-second. Topographic LIDAR data were processed and merged with the bathymetry to demonstrate the incorporation of recent third party data sources for several test areas. A primary application of a merged topographic/bathymetric elevation model is for user-defined shoreline delineation, in which the user decides on the tidal condition (for example, low or high water) to be superimposed on the elevation data to determine the spatial position of the water line. Such a use of merged topographic/bathymetric data could lead to the development of a shoreline zone, which could reduce redundant mapping efforts by federal, state, and local agencies by allowing them to customize their portrayals of the shoreline using a standard baseline elevation dataset.

Marine Technology Society Journal↗

USGS science for the Nation's changing coasts; shoreline change assessment

The coastline of the United States features some of the most popular tourist and recreational destinations in the world and is the site of intense residential, commercial, and industrial development. The coastal zone also has extensive and pristine natural areas, with diverse ecosystems providing essential habitat and resources that support wildlife, fish, and human use. Coastal erosion is a widespread process along most open-ocean shores of the United States that affects both developed and natural coastlines. As the coast changes, there are a wide range of ways that change can affect coastal communities, habitats, and the physical characteristics of the coast?including beach erosion, shoreline retreat, land loss, and damage to infrastructure. Global climate change will likely increase the rate of coastal change. A recent study of the U.S. Mid-Atlantic coast, for example, found that it is virtually certain that sandy beaches will erode faster in the future as sea level rises because of climate change. The U.S. Geological Survey (USGS) is responsible for conducting research on coastal change hazards, understanding the processes that cause coastal change, and developing models to predict future change. To understand and adapt to shoreline change, accurate information regarding the past and present configurations of the shoreline is essential. A comprehensive, nationally consistent analysis of shoreline movement is needed. To meet this national need, the USGS is conducting an analysis of historical shoreline changes along open-ocean coasts of the conterminous United States and parts of Alaska and Hawaii, as well as the coasts of the Great Lakes.

Fact Sheet↗

Influence of long term climate change on net infiltration at Yucca Mountain, Nevada

Net infiltration and recharge at Yucca Mountain, Nevada, a potential site for a high level nuclear waste repository, are determined both by the rock properties and past and future changes in climate. A 1-dimensional model was constructed to represent a borehole being drilled through the unsaturated zone. The rock properties were matched to the lithologies expected to be encountered in the borehole. As current paleoclimate theory assumes that 18O increases with wetter and cooler global climates, a past climate scenario, built on depletion of 18O from ocean sediments was used as a basis for climate change over the past 700,000 years. The climate change was simulated by assigning net infiltration values as a linear function of 8O. Assuming the rock properties, lithologies and climate scenarios are correct, simulations indicated that Yucca Mountain is not in steady state equilibrium at the surface (<75 meters) when compared to measured data, but that the system could be at steady state conditions at depths of >250 meters. Based on the cyclic climate inputs, the near surface is currently in a long term drying trend (for the last 3,000 years) yet recharge into the water table is continuing to occur at an average rate equivalent to the average input rate of the climate model, indicating that conditions at depth are damped out over very long time periods. The Paintbrush Tuff nonwelded units, positioned between the Tiva Canyon and Topopah Spring welded Tuff Members, do not appear to act as capillary barrier and therefore would not perch water. The low porosity vitric caprock and basal vitrophyre of the Topopah Spring Member, however, act as restrictive layers. The higher porosity rock directly above the caprock reduces the potential for the caprock to perch water leaving the basal vitrophyre as the most likely location for perched water to develop.

Conference Paper↗

Obliquity-paced Pliocene West Antarctic ice sheet oscillations

Thirty years after oxygen isotope records from microfossils deposited in ocean sediments confirmed the hypothesis that variations in the Earth's orbital geometry control the ice ages 1 , fundamental questions remain over the response of the Antarctic ice sheets to orbital cycles 2 . Furthermore, an understanding of the behaviour of the marine-based West Antarctic ice sheet (WAIS) during the 'warmer-than-present' early-Pliocene epoch ( 5–3 Myr ago) is needed to better constrain the possible range of ice-sheet behaviour in the context of future global warming 3 . Here we present a marine glacial record from the upper 600 m of the AND-1B sediment core recovered from beneath the northwest part of the Ross ice shelf by the ANDRILL programme and demonstrate well-dated, 40-kyr cyclic variations in ice-sheet extent linked to cycles in insolation influenced by changes in the Earth's axial tilt (obliquity) during the Pliocene. Our data provide direct evidence for orbitally induced oscillations in the WAIS, which periodically collapsed, resulting in a switch from grounded ice, or ice shelves, to open waters in the Ross embayment when planetary temperatures were up to 3 °C warmer than today 4 and atmospheric CO 2 concentration was as high as 400 p.p.m.v. (refs 5 , 6 ). The evidence is consistent with a new ice-sheet/ice-shelf model 7 that simulates fluctuations in Antarctic ice volume of up to +7 m in equivalent sea level associated with the loss of the WAIS and up to +3 m in equivalent sea level from the East Antarctic ice sheet, in response to ocean-induced melting paced by obliquity. During interglacial times, diatomaceous sediments indicate high surface-water productivity, minimal summer sea ice and air temperatures above freezing, suggesting an additional influence of surface melt 8 under conditions of elevated CO 2 .

Nature↗

Examining the role of unusually warm Indo‐Pacific sea‐surface temperatures in recent African droughts

Southern Africa (SA) and eastern Africa (EA) experienced a sequence of severe droughts in December–February (SA DJF) 2015–2016, October–December (EA OND) 2016 and March–April–May 2017 (EA MAM). This sequence contributed to severe food insecurity. While climate variability in these regions is very complex, the goal of this study is to analyse the role played by unusually warm Indo–Pacific SSTs, where unusual is defined as a 1‐in‐6 year event. We use observed sea‐surface temperatures (SST) and satellite–gauge rainfall observations, a 20‐member ensemble of Community Atmospheric Model version 5.1 simulations (CAM5), and a 40‐member ensemble of climate change simulations from the Community Earth Systems Model version 1 (CESM1) Large Ensemble Community Project (LENS) to explore climate conditions associated with warm events identified based on eastern and western Pacific SSTs. Our analysis suggests that strong El Niño's may be followed by warm western Pacific SST conditions, which can lead to conditions conducive to successive and potentially predictable droughts in SA DJF, EA OND and EA MAM. We show that different regions of warm SST are related to recent droughts—SA DJF: Niño 3.4; EA OND: western equatorial Pacific (WEP); and EA MAM: western North Pacific (WNP). For DJF and MAM, respectively, the CAM5 model driven with observed SST and the same model driven within a climate change experiment indicate that warmer El Niño's and WNP events produce more intense atmospheric responses, potentially associated with more severe droughts. OND climate seems to be strongly influenced by the Indian Ocean Dipole, which corresponds with some WEP events. Given global warming, we suggest that the extreme Niño 3.4 and west Pacific SST events responsible for 2015–2017 droughts are likely to reoccur, thus humanitarian agencies should prepare to predict and respond to multi‐year drought and substantial food insecurity in SA and EA.

Journal of the Royal Meteorological Society↗

Declining global per capita agricultural production and warming oceans threaten food security

Despite accelerating globalization, most people still eat food that is grown locally. Developing countries with weak purchasing power tend to import as little food as possible from global markets, suffering consumption deficits during times of high prices or production declines. Local agricultural production, therefore, is critical to both food security and economic development among the rural poor. The level of local agricultural production, in turn, will be determined by the amount and quality of arable land, the amount and quality of agricultural inputs (fertilizer, seeds, pesticides, etc.), as well as farm-related technology, practices and policies. This paper discusses several emerging threats to global and regional food security, including declining yield gains that are failing to keep up with population increases, and warming in the tropical Indian Ocean and its impact on rainfall. If yields continue to grow more slowly than per capita harvested area, parts of Africa, Asia and Central and Southern America will experience substantial declines in per capita cereal production. Global per capita cereal production will potentially decline by 14% between 2008 and 2030. Climate change is likely to further affect food production, particularly in regions that have very low yields due to lack of technology. Drought, caused by anthropogenic warming in the Indian and Pacific Oceans, may also reduce 21st century food availability in some countries by disrupting moisture transports and bringing down dry air over crop growing areas. The impacts of these circulation changes over Asia remain uncertain. For Africa, however, Indian Ocean warming appears to have already reduced rainfall during the main growing season along the eastern edge of tropical Africa, from southern Somalia to northern parts of the Republic of South Africa. Through a combination of quantitative modeling of food balances and an examination of climate change, this study presents an analysis of emerging threats to global food security.

Food Security↗

The crustal structure of the axis of the Great Valley, California, from seismic refraction measurements

In 1982 the U.S. Geological Survey collected six seismic refraction profiles in the Great Valley of California: three axial profiles with a maximum shot-to-receiver offset of 160 km, and three shorter profiles perpendicular to the valley axis. This paper presents the results of two-dimensional raytracing and synthetic seismogram modeling of the central axial profile. The crust of the central Great Valley is laterally heterogeneous along its axis, but generally consists of a sedimentary section overlying distinct upper, middle, and lower crustal units. The sedimentary rocks are 3-5 km thick along the profile, with velocities increasing with depth from 1.6 to 4.0 km/s. The basement (upper crust) consists of four units: 1. (1) a 1.0-1.5 km thick layer of velocity 5.4-5.8 km/s, 2. (2) a 3-4 km thick layer of velocity 6.0-6.3 km/s, 3. (3) a 1.5-3.0 km thick layer of velocity 6.5-6.6 km/s, and 4. (4) a laterally discontinuous, 1.5 km thick layer of velocity 6.8-7.0 km/s. The mid-crust lies at 11-14 km depth, is 5-8 km thick, and has a velocity of 6.6-6.7 km/s. On the northwest side of our profile the mid-crust is a low-velocity zone beneath the 6.8-7.0 km/s lid. The lower crust lies at 16-19 km depth, is 7-13 km thick, and has a velocity of 6.9-7.2 km/s. Crustal thickness increases from 26 to 29 km from NW to SE in the model. Although an unequivocal determination of crustal composition is not possible from P-wave velocities alone, our model has several geological and tectonic implications. We interpret the upper 7 km of basement on the northwest side of the profile as an ophiolitic fragment, since its thickness and velocity structure are consistent with that of oceanic crust. This fragment, which is not present 10-15 km to the west of the refraction profile, is probably at least partially responsible for the Great Valley gravity and magnetic anomalies, whose peaks lie about 10 km east of our profile. The middle and lower crust are probably gabbroic and the product of magmatic or tectonic underplating, or both. The crustal structure of the Great Valley is dissimilar to that of the adjacent Diablo Range, suggesting the existence of a fault or suture zone throughout the crust between these provinces.

California↗

Heat flow in the Arctic

Defines heat flow as the flux at the earth's solid surface of heat conducted from the interior; the heat-flow-unit (hfu) is on the order of 1-millionth calorie through each sq cm of the surface/sec, which is enough to melt a 4-mm layer of ice over the earth's surface/yr. Earth heat originates from radioactive decay of U, Th and K in the crust and mantle. Although land heat-flow measurements in the Arctic are too few for regional interpretation, those from Cape Thompson, Barrow and Cape Simpson, Northern Alaska are discussed and figured to show what they contribute to understanding of permafrost, climatic change and shoreline movements. Measuring thermal conductivity and gradient is much simpler in ocean basins than on land. Locations of such measurements are mapped, the results for the Alaskan quadrant in more detail. The sharp change in heat flow at the edge of the Alpha Cordillera, shown in a geothermal model, suggests that this feature is a huge accumulation of basalt, rather than mantle material or remnant of a foundering continent as previously postulated. Future Arctic heat flow studies are discussed.

Arctic↗

Watersheds and drainage networks

This topic is an overview of basic concepts about how the distribution of water on the Earth, with specific regard to watersheds, stream and river networks, and waterbodies are represented by geographic data. The flowing and non-flowing bodies of water on the earth’s surface vary in extent largely due to seasonal and annual changes in climate and precipitation. Consequently, modeling the detailed representation of surface water using geographic information is important. The area of land that collects surface runoff and other flowing water and drains to a common outlet location defines a watershed. Terrain and surface features can be naturally divided into watersheds of various sizes. Drainage networks are important data structures for modeling the distribution and movement of surface water over the terrain. Numerous tools and methods exist to extract drainage networks and watersheds from digital elevation models (DEMs). The cartographic representations of surface water are referred to as hydrographic features and consist of a snapshot at a specific time. Hydrographic features can be assigned general feature types, such as lake, pond, river, and ocean. Hydrographic features can be stored, maintained, and distributed for use through vector geospatial databases, such as the National Hydrography Dataset (NHD) for the United States.

Book chapter↗

Deglacial sea level history of the East Siberian Sea and Chukchi Sea margins

Deglacial (12.8–10.7 ka) sea level history on the East Siberian continental shelf and upper continental slope was reconstructed using new geophysical records and sediment cores taken during Leg 2 of the 2014 SWERUS-C3 expedition. The focus of this study is two cores from Herald Canyon, piston core SWERUS-L2-4-PC1 (4-PC1) and multicore SWERUS-L2-4-MC1 (4-MC1), and a gravity core from an East Siberian Sea transect, SWERUS-L2-20-GC1 (20-GC1). Cores 4-PC1 and 20-GC were taken at 120 and 115 m of modern water depth, respectively, only a few meters above the global last glacial maximum (LGM; ∼ 24 kiloannum or ka) minimum sea level of ∼ 125–130 meters below sea level (m b.s.l.). Using calibrated radiocarbon ages mainly on molluscs for chronology and the ecology of benthic foraminifera and ostracode species to estimate paleodepths, the data reveal a dominance of river-proximal species during the early part of the Younger Dryas event (YD, Greenland Stadial GS-1) followed by a rise in river-intermediate species in the late Younger Dryas or the early Holocene (Preboreal) period. A rapid relative sea level rise beginning at roughly 11.4 to 10.8 ka ( ∼ 400 cm of core depth) is indicated by a sharp faunal change and unconformity or condensed zone of sedimentation. Regional sea level at this time was about 108 m b.s.l. at the 4-PC1 site and 102 m b.s.l. at 20-GC1. Regional sea level near the end of the YD was up to 42–47 m lower than predicted by geophysical models corrected for glacio-isostatic adjustment. This discrepancy could be explained by delayed isostatic adjustment caused by a greater volume and/or geographical extent of glacial-age land ice and/or ice shelves in the western Arctic Ocean and adjacent Siberian land areas.

Chukchi Sea, East Siberian Sea↗

Olivine friction at the base of oceanic seismogenic zones

We investigate the strength and frictional behavior of olivine aggregates at temperatures and effective confining pressures similar to those at the base of the seismogenic zone on a typical ridge transform fault. Triaxial compression tests were conducted on dry olivine powder (grain size ???60 ??m) at effective confining pressures between 50 and 300 MPa (using Argon as a pore fluid), temperatures between 600??C and 1000??C, and axial displacement rates from 0.06 to 60 ??m/s (axial strain rates from 3 ?? 10-6 to 3 ?? 10-3 s-1). Yielding shows a negative pressure dependence, consistent with predictions for shear enhanced compaction and with the observation that samples exhibit compaction during the initial stages of the experiments. A combination of mechanical data and microstructural observations demonstrate that deformation was accommodated by frictional processes. Sample strengths were pressure-dependent and nearly independent of temperature. Localized shear zones formed in initially homogeneous aggregates early in the experiments. The frictional response to changes in loading rate is well described by rate and state constitutive laws, with a transition from velocity-weakening to velocity-strengthening at 1000??C. Microstructural observations and physical models indicate that plastic yielding of asperities at high temperatures and low axial strain rates stabilizes frictional sliding. Extrapolation of our experimental data to geologic strain rates indicates that a transition from velocity weakening to velocity strengthening occurs at approximately 600??C, consistent with the focal depths of earthquakes in the oceanic lithosphere. Copyright 2007 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗

Hydrogeology and simulation of groundwater flow in fractured rock in the Newark basin, Rockland County, New York

Groundwater in the Newark basin aquifer flows primarily through discrete water-bearing zones parallel to the strike and dip of bedding, whereas flow perpendicular to the strike is restricted, thereby imparting anisotropy to the groundwater flow field. The finite-element model SUTRA was used to represent bedrock structure in the aquifer by spatially varying the orientation of the hydraulic conductivity tensor to reflect variations in the strike and dip of the bedding. Directions of maximum and medium hydraulic conductivity were oriented parallel to the bedding, and the direction of minimum hydraulic conductivity was oriented perpendicular to the bedding. Groundwater flow models were prepared to simulate local flow in the vicinity of the Spring Valley well field and regional flow through the Newark basin aquifer. The Newark basin contains sedimentary rocks deposited as alluvium during the Late Triassic and is one of a series of basins that developed when Mesozoic rifting of the super continent Pangea created the Atlantic Ocean. The westward-dipping basin is filled with interbedded facies of coarse-grained to fine-grained rocks that were intruded by diabase associated with Jurassic volcanism. The Newark basin aquifer is bounded to the north and east by the Palisades sill and to the west by the Ramapo Fault. Although the general dip of bedding is toward the fault, mapping of conglomerate beds indicates the rocks are folded into broad anticlines and synclines. An alternative, more uniform pattern of regional structure, based on interpolated strike and dip measurements from a number of sources, has also been proposed. Two groundwater flow models (A for the former type of bedrock structure and B for the latter type) were developed to represent these alternative depictions of bedrock structure. Transient simulations were calibrated to reproduce measured water-level recoveries in a 9.3 mi&sup2 area surrounding the Spring Valley well field during a 5-day aquifer test in 1992. The models represented a 330-ft thick rock mass divided vertically into 10 equally spaced layers and were calibrated through nonlinear regression. Results of model B best matched the observed water-level recoveries with an estimated hydraulic conductivity of 9.5 ft/day, specific storage of 7.6 x 10 -6 ft -1, and K max : K min anisotropy ratio (hydraulic conductivity parallel to bedding: perpendicular to bedding) of 72:1. Model error was 50 percent greater in model A because the assumed structure did not match the actual strike of bedding in this area. Steady-state simulations of regional flow through the 85.4-mi2 modeled extent of the Newark basin aquifer represented both the alluvial aquifer beneath the Mawah River and the fractured bedrock. The rock mass was divided into two aquifer units: an upper 500-ft thick unit divided into 10 equally spaced layers through which most groundwater is assumed to flow and a lower unit divided into 7 layers with increasing thickness. Models were calibrated through nonlinear regression to average water levels measured in 140 wells from August 2005 through April 2007. Water levels simulated using the two models were similar and generally matched those observed, and the average recharge rate estimated using both models was 19 inches/year for the simulated period. Estimated transmissivity parallel to the strike of bedding (1,100 ft&sup2/d) was uniform in two transmissivity (T) zones in model A, but in model B the transmissivity of a high T zone (1,600 ft&sup2/d), delineated on the basis of aquifer test data, was slightly greater than in a low T zone (1,300 ft&sup2/d). The K max : K min anisotropy was estimated to be 58:1 in model A and 410:1 in model B, so the proportion of flow perpendicular to bedding is less in model B than in model A. Distributions of groundwater age simulated with models A and B are similar and indicate that most shallow ground-water (225 ft below the bedrock surface) is 5 t

Scientific Investigations Report↗

Evolution of CAM and C 4 carbon-concentrating mechanisms

Mechanisms for concentrating carbon around the Rubisco enzyme, which drives the carbon-reducing steps in photosynthesis, are widespread in plants; in vascular plants they are known as crassulacean acid metabolism (CAM) and C 4 photosynthesis. CAM is common in desert succulents, tropical epiphytes, and aquatic plants and is characterized by nighttime fixation of CO2. The proximal selective factor driving the evolution of this CO2-concentrating pathway is low daytime CO2, which results from the unusual reverse stomatal behavior of terrestrial CAM species or from patterns of ambient CO2 availability for aquatic CAM species. In terrestrials the ultimate selective factor is water stress that has selected for increased water use efficiency. In aquatics the ultimate selective factor is diel fluctuations in CO2 availability for palustrine species and extreme oligotrophic conditions for lacustrine species. C4 photosynthesis is based on similar biochemistry but carboxylation steps are spatially separated in the leaf rather than temporally as in CAM. This biochemical pathway is most commonly associated with a specialized leaf anatomy known as Kranz anatomy; however, there are exceptions. The ultimate selective factor driving the evolution of this pathway is excessively high photorespiration that inhibits normal C3 photosynthesis under high light and high temperature in both terrestrial and aquatic habitats. CAM is an ancient pathway that likely has been present since the Paleozoic era in aquatic species from shallow-water palustrine habitats. While atmospheric CO2 levels have undoubtedly affected the evolution of terrestrial plant carbon-concentrating mechanisms, there is reason to believe that past atmospheric changes have not played as important a selective role in the aquatic milieu since palustrine habitats today are not generally carbon sinks, and the selective factors driving aquatic CAM are autogenic. Terrestrial CAM, in contrast, is of increasing selective value under extreme water deficits, and undoubtedly, high Mesozoic CO2 levels reduced the amount of landscape perceived by plants as water limited. Late Tertiary and Quaternary reductions in atmospheric CO2, coupled with increasing seasonality, were probably times of substantial species radiation and ecological expansion for CAM plants. C4 photosynthesis occurs in only about half as many families as CAM, and three-fourths of C4 species are either grasses or sedges. Molecular phylogenies indicate C4 is a more recent innovation than CAM and that it originated in the mid-Tertiary, 20–30 Ma, although some data support an earlier origin. While the timing of the origin of C4 remains controversial, the nearly explosive increase in C4 species is clearly documented in the late Miocene, 4–7 Ma. Increasing seasonality has been widely suggested as an important climatic stimulus for this C4 expansion. Alternatively, based on models of photosynthetic quantum yield at different temperatures and CO2 concentration, it has been hypothesized that the late Miocene C4 expansion resulted from declining atmospheric CO2 levels. This model is most appropriate for explaining the transition from C3 grasslands to C4 grasslands but by itself may not be sufficient to explain the more likely scenario of a late Miocene transition from C3 woodland/ savanna to C4 grasslands. A largely unexplored hypothesis is that climatic changes in late Miocene altered disturbance regimes, in particular the incidence of fires, which today are often associated with maintenance of C4 grasslands. Oceanic charcoal sediments that appear to represent Aeolian deposits from continental wildfires follow a strikingly similar pattern of explosive increase in late Miocene. Climate, CO2, and disturbance are not mutually exclusive explanations and probably all acted in concert to promote the expansion of C4 grasslands. More recently, late Quaternary changes in CO2 may have been responsible for driving major changes in the landscape distribution of C4 species. The theory is sound; however, many of the studies cited in support of this model are open to alternative interpretations, and none has eliminated climatic factors as important selective agents. CAM and C4 evolution required coupling of biochemical pathways with structural changes in photosynthetic tissues, succulence in CAM and Kranz in C4. This was apparently accomplished by piecemeal evolution beginning with mechanisms for recapturing respiratory CO2, although this need not have been so in aquatic CAM species. It has been proposed that the extreme rarity of both pathways in the same plant results from biochemical and structural incompatibilities (Sage 2002). Equally important is the fact that the selective environments are quite different, with CAM evolution thriving on stressful sites inhospitable to C3 species whereas C4 evolution has selected for rapid growth capable of outcompeting associated C3 plants.

International Journal of Plant Sciences↗

Geophysical evidence for the evolution of the California Inner Continental Borderland as a metamorphic core complex

We use new seismic and gravity data collected during the 1994 Los Angeles Region Seismic Experiment (LARSE) to discuss the origin of the California Inner Continental Borderland (ICB) as an extended terrain possibly in a metamorphic core complex mode. The data provide detailed crustal structure of the Borderland and its transition to mainland southern California. Using tomographic inversion as well as traditional forward ray tracing to model the wide-angle seismic data, we find little or no sediments, low (≤6.6 km/s) P wave velocity extending down to the crust-mantle boundary, and a thin crust (19 to 23 km thick). Coincident multichannel seismic reflection data show a reflective lower crust under Catalina Ridge. Contrary to other parts of coastal California, we do not find evidence for an underplated fossil oceanic layer at the base of the crust. Coincident gravity data suggest an abrupt increase in crustal thickness under the shelf edge, which represents the transition to the western Transverse Ranges. On the shelf the Palos Verdes Fault merges downward into a landward dipping surface which separates "basement" from low-velocity sediments, but interpretation of this surface as a detachment fault is inconclusive. The seismic velocity structure is interpreted to represent Catalina Schist rocks extending from top to bottom of the crust. This interpretation is compatible with a model for the origin of the ICB as an autochthonous formerly hot highly extended region that was filled with the exhumed metamorphic rocks. The basin and ridge topography and the protracted volcanism probably represent continued extension as a wide rift until ∼13 m.y. ago. Subduction of the young and hot Monterey and Arguello microplates under the Continental Borderland, followed by rotation and translation of the western Transverse Ranges, may have provided the necessary thermomechanical conditions for this extension and crustal inflow.

California↗

Updates to the regional groundwater-flow model of the New Jersey Coastal Plain, 1980–2013

A 21-layer three-dimensional transient groundwater-flow model of the New Jersey Coastal Plain was developed and calibrated by the U.S. Geological Survey (USGS) in cooperation with the New Jersey Department of Environmental Protection to simulate groundwater-flow conditions during 1980–2013, incorporating average annual groundwater withdrawals and average annual groundwater recharge. This model is the third version of the New Jersey Coastal Plain regional groundwater-flow model that was initially developed as part of the USGS Regional Aquifer System Analysis (RASA) program. The model simulates groundwater flow in 11 aquifers and 10 intervening confining units of the New Jersey Coastal Plain to provide a regional overview of groundwater conditions. Averaged groundwater withdrawal data for 1980 to 2013 were used in the model. The 11 aquifers in New Jersey are, from shallowest to deepest, the Holly Beach water-bearing zone and the confined Cohansey aquifer in Cape May County; the Rio Grande water-bearing zone; the Atlantic City 800-foot sand; the Piney Point, Vincentown, and Wenonah-Mount Laurel aquifers; the Englishtown aquifer system; and the upper, middle, and lower aquifers of the Potomac-Raritan-Magothy (PRM) aquifer system. The model was developed with the MODFLOW–2005 numerical code and the UCODE parameter estimation technique and calibrated using water-level and base-flow observations. A total of 3,453 water-level observations from 392 wells in New Jersey and 48 wells in Delaware from 1983 to 2013 were used in model calibration, which includes historical water-level trends for 29 wells in New Jersey during 1980–2013 presented in time-series hydrographs. In addition, derived observations also were included by calculating the vertical gradient at 33 pairs of nested observation wells in New Jersey, for a total of 210 observations. Changes in water levels over time were calculated for 134 wells in New Jersey and four wells in Delaware where water levels had varied substantially (approximately 10 ft) over the 30-year span of synoptic water-level measurements, for a total of 767 observations. A total of 1,485 base-flow observations in 47 surface-water basins in New Jersey from 1980 to 2013 were used in model calibration. Updates to the groundwater-flow model include the conversion to a fully three-dimensional model from the previous quasi-three-dimensional model. The new model will allow for potential future uses such as particle tracking or simulation of variable-density groundwater flow that could not be accomplished with earlier versions of the model. Spatially and temporally variable recharge estimated by using a soil-water balance model resulted in a spatially and temporally finer discretization. The Rio Grande water-bearing zone was added to the model as an aquifer layer to refine estimates of simulated flow in Atlantic and Cape May Counties, New Jersey. Hydrogeologic parameters were updated to include the confining units in New Jersey and corresponding hydrogeologic units in Delaware and eastern Maryland. The simulated water levels for the New Jersey Coastal Plain aquifers were compared to water-level measurements made during 1980–2013. The average residual for 4,243 water-level observations for New Jersey (simulated water levels minus measured water levels) is 1.5 feet. The simulated water-level contours for the confined aquifers for 2013 were compared to potentiometric surfaces produced from water levels measured during 2013. Simulated water levels generally matched the 2013 potentiometric surfaces of the confined aquifers in the areas of large withdrawals. Hydrographs of wells in the confined Coastal Plain aquifers of New Jersey show that simulated water levels generally match the magnitude and seasonal variation of the observed water levels. Hydrographs of base flow for the 47 streamgaging stations in New Jersey indicate that most of the simulated and estimated data match reasonably well. Groundwater withdrawals are an important resource for water supply, agricultural, industrial, and commercial needs in the New Jersey Coastal Plain. Groundwater withdrawals from the New Jersey Coastal Plain aquifers have resulted in persistent, regionally extensive cones of depression in the Englishtown aquifer system and Wenonah-Mount Laurel aquifer in Ocean and Monmouth Counties; Wenonah-Mount Laurel and upper, middle, and lower PRM aquifers in Camden County; and Atlantic City 800-foot sand in Atlantic County. Because hydrologic stresses and water-management needs change with time, periodic updates to the groundwater-flow model are required to provide current information about hydrologic conditions in the New Jersey Coastal Plain and to maintain its usefulness as a tool to manage water resources and develop water-resource strategies. The current updates will support the continued application of this model as a tool for evaluating the regional effects of changes in groundwater withdrawals and of current and potential future water-management strategies on groundwater levels in the New Jersey Coastal Plain.

New Jersey Coastal Plain↗