Geology Reports⌕ Search

SEARCH · Geology Reports

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 1,153 records · Page 64Linked to original sources

Numerical modeling of late Glacial Laurentide advance of ice across Hudson Strait: Insights into terrestrial and marine geology, mass balance, and calving flux

A time-dependent finite element model was used to reconstruct the advance of ice from a late Glacial dome on northern Quebec/Labrador across Hudson Strait to Meta Incognita Peninsula (Baffin Island) and subsequently to the 9.9–9.6 ka 14 C Gold Cove position on Hall Peninsula. Terrestrial geological and geophysical information from Quebec and Labrador was used to constrain initial and boundary conditions, and the model results are compared with terrestrial geological information from Baffin Island and considered in the context of the marine event DC-0 and the Younger Dryas cooling. We conclude that advance across Hudson Strait from Ungava Bay to Baffin Island is possible using realistic glacier physics under a variety of reasonable boundary conditions. Production of ice flux from a dome centered on northeastern Quebec and Labrador sufficient to deliver geologically inferred ice thickness at Gold Cove (Hall Peninsula) appears to require extensive penetration of sliding south from Ungava Bay. The discharge of ice into the ocean associated with advance and retreat across Hudson Strait does not peak at a time coincident with the start of the Younger Dryas and is less than minimum values proposed to influence North Atlantic thermohaline circulation; nevertheless, a significant fraction of freshwater input to the North Atlantic may have been provided abruptly and at a critical time by this event.

Hudson Strait↗

Petrology and tectonics of Phanerozoic continent formation: From island arcs to accretion and continental arc magmatism

Mesozoic continental arcs in the North American Cordillera were examined here to establish a baseline model for Phanerozoic continent formation. We combine new trace-element data on lower crustal xenoliths from the Mesozoic Sierra Nevada Batholith with an extensive grid-based geochemical map of the Peninsular Ranges Batholith, the southern equivalent of the Sierras. Collectively, these observations give a three-dimensional view of the crust, which permits the petrogenesis and tectonics of Phanerozoic crust formation to be linked in space and time. Subduction of the Farallon plate beneath North America during the Triassic to early Cretaceous was characterized by trench retreat and slab rollback because old and cold oceanic lithosphere was being subducted. This generated an extensional subduction zone, which created fringing island arcs just off the Paleozoic continental margin. However, as the age of the Farallon plate at the time of subduction decreased, the extensional environment waned, allowing the fringing island arc to accrete onto the continental margin. With continued subduction, a continental arc was born and a progressively more compressional environment developed as the age of subducting slab continued to young. Refinement into a felsic crust occurred after accretion, that is, during the continental arc stage, wherein a thickened crustal and lithospheric column permitted a longer differentiation column. New basaltic arc magmas underplate and intrude the accreted terrane, suture, and former continental margin. Interaction of these basaltic magmas with pre-existing crust and lithospheric mantle created garnet pyroxenitic mafic cumulates by fractional crystallization at depth as well as gabbroic and garnet pyroxenitic restites at shallower levels by melting of pre-existing lower crust. The complementary felsic plutons formed by these deep-seated differentiation processes rose into the upper crust, stitching together the accreted terrane, suture and former continental margin. The mafic cumulates and restites, owing to their high densities, eventually foundered into the mantle, leaving behind a more felsic crust. Our grid-based sampling allows us to estimate an unbiased average upper crustal composition for the Peninsular Ranges Batholith. Major and trace-element compositions are very similar to global continental crust averaged over space and time, but in detail, the Peninsular Ranges are slightly lower in compatible to mildly incompatible elements, MgO, Mg#, V, Sc, Co, and Cr. The compositional similarities suggest a strong arc component in global continental crust, but the slight discrepancies suggest that additional crust formation processes are also important in continent formation as a whole. Finally, the delaminated Sierran garnet pyroxenites have some of the lowest U/Pb ratios ever measured for silicate rocks. Such material, if recycled and stored in the deep mantle, would generate a reservoir with very unradiogenic Pb, providing one solution to the global Pb isotope paradox.

Earth and Planetary Science Letters↗

Protracted construction of gabbroic crust at a slow spreading ridge: Constraints from 206 Pb/ 238 U zircon ages from Atlantis Massif and IODP Hole U1309D (30°N, MAR)

Sensitive high-resolution ion microprobe (SHRIMP) U-Pb zircon ages of 24 samples from oceanic crust recovered in Integrated Ocean Drilling Program (IODP) Hole U1309D and from the surface of Atlantis Massif, Mid-Atlantic Ridge (MAR) (30°N) document a protracted history of accretion in the footwall to an oceanic detachment fault. Ages for 18 samples of evolved Fe-Ti oxide gabbro and felsic dikes collected 40–1415 m below seafloor in U1309D yield a weighted mean of 1.20 ± 0.03 Ma (mean square of weighted deviates = 7.1). However, the ages range from 1.08 ± 0.07 Ma and 1.28 ± 0.05 Ma indicating crustal construction occurred over a minimum of 100–200 ka. The zircon ages, along with petrologic observations, indicate at least 2 major periods of intrusive activity with age peaks separated by 70 ka. The oldest ages are observed below 600 mbsf, an observation inconsistent with models requiring constant depth melt intrusion beneath a detachment fault. The data are most consistent with a “multiple sill” model whereby sills intrude at random depths below the ridge axis over a length scale greater than 1.4 km. Zircon ages from broadly spaced samples collected along the southern ridge of Atlantis Massif yield a detachment fault slip rate of 28.7 ± 6.7 mm/a and imply significant asymmetric plate spreading (up to 100% on the North American plate) for at least 200 ka during core complex formation.

Geochemistry, Geophysics, Geosystems↗

Thin, low‐velocity crust beneath the southern Yukon‐Tanana Terrane, east central Alaska: Results from Trans‐Alaska crustal transect refraction/wide‐angle reflection data

A seismic refraction/wide‐angle reflection survey for the Trans‐Alaska Crustal Transect program reveals a thin, reflective crust beneath the southern Yukon‐Tanana terrane (YTT) in east central Alaska. These data are the first detailed refraction survey of the southern YTT and compose a 130‐km‐long reversed profile along the Alaska and Richardson highways. Results from this study indicate that low‐velocity (≤ 6.4 km/s) rocks extend to approximately 27 km in depth. Based on these low velocities and an average Poisson's ratio of 0.23 determined for depths of ≤27 km, an overall silicic composition is interpreted for this portion of the crust beneath the Yukon‐Tanana terrane. From approximately 8 to 27 km depth the crust exhibits an increase in reflectivity. This middle to lower crustal reflectivity is modeled as alternating high‐ and low‐velocity lamellae with an average velocity of 6.1 km/s at 10 km depth to an average velocity of 6.4 km/s at 27 km depth. Beneath these reflective, low‐velocity rocks a 3‐ to 5‐km‐thick, 7.0 km/s basal crustal layer produces a prominent reflection that extends to offsets of up to 280 km. The crust‐mantle boundary, modeled at an average depth of 30 km, produces a variable PmP reflection, which may indicate lateral heterogeneity of this boundary, and a weak and emergent Pn refraction with a velocity of 8.2 km/s. We interpret the crustal section as follows: the low‐velocity rocks of the southern YTT extend from the surface to depths of approximately 10 km; underthrust Mesozoic flysch of the Kahiltna terrane, rocks of the Gravina arc, and basement of the Wrangellia(?) terrane extend from 10 to 27 km depth; a 3‐ to 5‐km‐thick layer of mantle‐derived mafic rocks, relic oceanic crust, or Wrangellia(?) terrane lower crust extends from 27 to approximately 30 km depth; a tectonically young Moho beneath the southern YTT is found at an average depth of 30 km; and it is underlain by a mantle that may be relatively cool and/or olivine rich. In this interpretation, the Yukon‐Tanana terrane is a thin‐skinned terrane. Our results indicate that tectonic, and possibly magmatic, underplating has played a significant role in crustal growth for central Alaska.

Alaska↗

Landscape dominance of introduced herpetofauna on an oceanic island

Habitat loss and fragmentation can negatively impact native wildlife and facilitate establishment of introduced species. On islands, introduced species are a primary cause of extinction and can alter community membership through predation or competition for resources. Consequently, elucidating the distribution of introduced and native species can improve understanding of the potential synergistic effects of land use and introduced species on native island species. The island of Saipan in the Commonwealth of the Northern Mariana Islands has a long history of herpetofaunal introductions. However, the distribution of native and introduced herpetofauna on Saipan is poorly understood, and the effects of land use and introduced species on the native herpetofauna are largely unknown. In 2018, we randomly selected sites on Saipan from established island-wide transects and used multiple detection methods, repeated surveys, and single-season occupancy models to account for imperfect detection and construct species distribution models. We investigated the role of biotic and abiotic factors in explaining occupancy and detection of Saipan’s native and introduced herpetofauna. We recorded 2568 observations representing 1 amphibian and 11 reptile species. Introduced species were encountered 3.8 times more frequently and were 33% more broadly distributed than native species. Occupancy for the native forest dwelling species Emoia caeruleocauda and Perochirus ateles was positively associated with elevation and enhanced vegetation index, suggesting that human habitat modification at lower elevations may be reducing their suitable habitat.We found evidence of a range contraction for the imperiled P. ateles and suggest potential avenues for implementing targeted conservation of native herpetofauna.

Global Ecology and Conservation↗

Characterization of calcium isotopes in natural and synthetic barite

The mineral barite (BaSO4) accommodates calcium in its crystal lattice, providing an archive of Ca-isotopes in the highly stable sulfate mineral. Holocene marine (pelagic) barite samples from the major ocean basins are isotopically indistinguishable from each other (??44/40Ca = -2.01 ?? 0.15???) but are different from hydrothermal and cold seep barite samples (??44/40Ca = -4.13 to -2.72???). Laboratory precipitated (synthetic) barite samples are more depleted in the heavy Ca-isotopes than pelagic marine barite and span a range of Ca-isotope compositions, ??44/40Ca = -3.42 to -2.40???. Temperature, saturation state, a Ba2 + / a SO42 -, and aCa2+/aBa2+ each influence the fractionation of Ca-isotopes in synthetic barite; however, the fractionation in marine barite samples is not strongly related to any measured environmental parameter. First-principles lattice dynamical modeling predicts that at equilibrium Ca-substituted barite will have much lower 44Ca/40Ca than calcite, by -9??? at 0 ??C and -8??? at 25 ??C. Based on this model, none of the measured barite samples appear to be in isotopic equilibrium with their parent solutions, although as predicted they do record lower ??44/40Ca values than seawater and calcite. Kinetic fractionation processes therefore most likely control the extent of isotopic fractionation exhibited in barite. Potential fractionation mechanisms include factors influencing Ca2+ substitution for Ba2+ in barite (e.g. ionic strength and trace element concentration of the solution, competing complexation reactions, precipitation or growth rate, temperature, pressure, and saturation state) as well as nucleation and crystal growth rates. These factors should be considered when investigating controls on isotopic fractionation of Ca2+ and other elements in inorganic and biogenic minerals. ?? 2008 Elsevier Ltd.

Geochimica et Cosmochimica Acta↗

Paleozoic ophiolitic assemblages within the southern New England orogen of eastern Australia: Implications for growth of the Gondwana margin

Several ophiolitic assemblages occur in the southern New England orogen. The development of these rocks and their relations to the rest of the orogen have major implications for the tectonic evolution of eastern Gondwana. A major, narrow but elongate belt of Early Cambrian suprasubduction zone ophiolite crops out along and near the PeelManning Fault System and is juxtaposed against younger arc and subduction complex terranes. No pre-Permian links with the rest of the New England orogen have been established for this terrane. It may represent portions of Lachlan Fold Belt basement which underlies younger, westward overthrust New England terranes, and has been diapirically emplaced at higher crustal levels as serpentinite-matrix melange. Middle to Late Devonian ophiolitic rocks in the Yarras Complex comprise basement to the Birpai subterrane and represent a crustal cross section through a rifted island arc. Correlatives of this terrane also occur within the more extensive Gamilaroi terrane to the west of which deeper crustal levels are not exposed. The various components of serpentinite-matrix melange in the Ngamba terrane at Port Macquarie superficially appear to represent a dismembered ophiolite association. However, the various components of the melange exhibit a wide range of ages, metamorphism, and tectonic affinities, rendering a genetically related origin unlikely. This terrane includes fragments of ocean floor accreted into a Late Devonian subduction complex, which was later affected by Early Carboniferous forearc serpentinite diapirism and high Mg series magmatism. Zircon inheritance in Triassic dikes, which intrude the melange attest to the development or later emplacement of this forearc region over an older Lachlan Fold Belt basement. Ultramafic rocks of the Bundjalung terrane in the east of the New England region probably formed at deep levels in an intraoceanic island arc and are intruded by boninitic dikes. The tectonic development of the NEO was significantly more complicated than has been suggested in earlier published models. Periodic accretion of island arc systems, some of which are now represented by suprasubduction zone ophiolites, to the eastern margin of Gondwana suggests multiple phases of subduction with the possibility of polarity reversals throughout the history of accretion. Lateral accretion was not the only means by which Gondwana continental crust grew, and there was considerable postaccretion continentward overthrusting of younger terranes.

New South Wales↗

Petrology and age of volcanic-arc rocks from the continental margin of the Bering Sea: Implications for Early Eocene relocation of plate boundaries

Eocene volcanic flow and dike rocks from the Beringian margin have arc characteristics, implying a convergent history for this region during the early Tertiary. The extrusive rocks are basalt, basaltic andesite, andesite, and minor dacite and rhyolite. The intrusive sample is from a quartz diorite dike intruding serpentinized peridotite. Major-element oxide contents, particularly FeO*/MgO versus SiO 2 , identify both tholeiitic and calc-alkalic basalt; more silicic lavas have calc-alkalic affinities. Consistent with volcanic-arc compositions, spidergrams show pronounced Nb–Ta depletion and alkali enrichment relative to light-rare-earth-element ( LREE ) abundance. Chondrite-normalized REE plots show relatively flat patterns, with only slight LREE enrichment for tholeiitic compositions and greater LREE enrichment and lower heavy-rare-earth-element ( HREE ) abundance for calc-alkalic compositions. The samples, particularly those with calc-alkalic compositions, are rich in plagioclase that is strongly zoned; the more silicic samples contain orthopyroxene, clinopyroxene, and primary amphibole. The quartz diorite dike contains iron-rich almandine phenocrysts that appear to be magmatic, suggesting emplacement at great depth near the base of the crust or upper mantle.Chemical and mineralogical compositions are similar to those of modern Aleutian-arc lavas. They also resemble volcanic-arc compositions from western mainland Alaska, although greater chemical diversity and a stronger continental influence are observed in the Alaskan mainland rocks.Early Eocene ages of 54.4–50.2 Ma for the Beringian samples are well constrained by conventional K–Ar ages of nine plagioclase separates and by concordant 40 Ar/ 39 Ar incremental heating and total-fusion experiments. A concordant U–Pb zircon age of 53 Ma for the quartz-diorite dike is in good agreement with the K–Ar data.Plate motion studies of the North Pacific Ocean indicate more northerly directed subduction prior to the Tertiary and a continuous belt of arc-type volcanism extending from Siberia, along the Beringian margin, into mainland Alaska. Around 56 Ma (chron 25–24), subduction changed to a more westerly direction and subduction-related volcanism ceased for most of mainland Alaska. The increasingly oblique angle of convergence should have ended subduction along the Beringian margin as well. However, consistent ages of 54–50 Ma indicate a final pulse in arc-type magmatism during this period of plate adjustment, which may be explained by three different models: (1) The northern and central part of the Beringian margin maintained a higher angle of convergence, allowing a final pulse of arc-type magmatism. (2) The rocks erupted in an early, or proto, Aleutian arc and were rafted against the continental margin along transform faults. (3) The rocks erupted along a leaky transform fault, analogous to calc-alkalic volcanism in the southern California borderland.

Alaska↗

Geochemistry of southern Pagan Island lavas, Mariana arc: The role of subduction zone processes

New major and trace element abundances, and Pb, Sr, and Nd isotopic ratios of Quaternary lavas from two adjacent volcanoes (South Pagan and the Central Volcanic Region, or CVR) located on Pagan Island allow us to investigate the mantle source (i.e., slab components) and melting dynamics within the Mariana intra-oceanic arc. Geologic mapping reveals a pre-caldera (780-9.4ka) and post-caldera (<9.4ka) eruptive stage for South Pagan, whereas the eruptive history of the older CVR is poorly constrained. Crystal fractionation and magma mixing were important crustal processes for lavas from both volcanoes. Geochemical and isotopic variations indicate that South Pagan and CVR lavas, and lavas from the northern volcano on the island, Mt. Pagan, originated from compositionally distinct parental magmas due to variations in slab contributions (sediment and aqueous fluid) to the mantle wedge and the extent of mantle partial melting. A mixing model based on Pb and Nd isotopic ratios suggests that the average amount of sediment in the source of CVR (~2.1%) and South Pagan (~1.8%) lavas is slightly higher than Mt. Pagan (~1.4%) lavas. These estimates span the range of sediment-poor Guguan (~1.3%) and sediment-rich Agrigan (~2.0%) lavas for the Mariana arc. Melt modeling demonstrates that the saucer-shaped normalized rare earth element (REE) patterns observed in Pagan lavas can arise from partial melting of a mixed source of depleted mantle and enriched sediment, and do not require amphibole interaction or fractionation to depress the middle REE abundances of the lavas. The modeled degree of mantle partial melting for Agrigan (2-5%), Pagan (3-7%), and Guguan (9-15%) lavas correlates with indicators of fluid addition (e.g., Ba/Th). This relationship suggests that the fluid flux to the mantle wedge is the dominant control on the extent of partial melting beneath Mariana arc volcanoes. A decrease in the amount of fluid addition (lower Ba/Th) and extent of melting (higher Sm/Yb), and an increase in the sediment contribution (higher Th/Nb, La/Sm, and Pb isotopic ratios) from Mt. Pagan to South Pagan could reflect systematic cross-arc or irregular along-arc melting variations. These observations indicate that the length scale of compositional heterogeneity in the mantle wedge beneath Mariana arc volcanoes is small (~10km).

Contributions to Mineralogy and Petrology↗

Generating a reference flow network with improved connectivity to support durable data integration and reproducibility in the coterminous US

This report presents a reference flow network for the conterminous United States that is built from the best available information from the U.S. Geological Survey, the National Oceanic and Atmospheric Administration National Weather Service, and the U.S. Environmental Protection Agency. The work is intended to support durable data integration and reproducibility. Originating from the National Hydrography Dataset Plus (NHDPlus) V2.1, the reference flow network incorporates network connectivity enhancements from federal agency efforts. After incorporating these network improvements, many original NHDPlus attributes were regenerated to enable network navigation and related operations. After introducing the motivation and background for this work, this report describes the attribute generation workflow and data quality checks that were performed in preparation of the dataset. The reference flow network follows the NHDPlus data model and is described using terms defined in the Mainstem and Drainage Basin logical model and WaterML2 Part3: Surface Hydrology Features conceptual model.

Environmental Modelling and Software↗

Geologic interpretation and multibeam bathymetry of the sea floor in southeastern Long Island Sound

Digital terrain models (DTMs) produced from multibeam echosounder (MBES) bathymetric data provide valuable base maps for marine geological interpretations (e.g. Todd and others, 1999; Mosher and Thomson, 2002; ten Brink and others, 2004; Poppe and others, 2006a,b). These maps help define the geological variability of the sea floor (one of the primary controls of benthic habitat diversity); improve our understanding of the processes that control the distribution and transport of bottom sediments, the distribution of benthic habitats and associated infaunal community structures; and provide a detailed framework for future research, monitoring, and management activities. The bathymetric survey interpreted herein (National Oceanic and Atmospheric Administration (NOAA) survey H11255) covers roughly 95 km? of sea floor in southeastern Long Island Sound (fig. 1). This bathymetry has been examined in relation to seismic reflection data collected concurrently, as well as archived seismic profiles acquired as part of a long-standing geologic mapping partnership between the State of Connecticut and the U.S. Geological Survey (USGS). The objective of this work was to use these geophysical data sets to interpret geomorphological attributes of the sea floor in terms of the Quaternary geologic history and modern sedimentary processes within Long Island Sound.

Open-File Report↗

Geologic interpretation and multibeam bathymetry of the sea floor in the vicinity of the Race, eastern Long Island Sound

Digital terrain models (DTMs) produced from multibeam bathymetric data provide valuable base maps for marine geological interpretations (Todd and others, 1999; Mosher and Thomson, 2002; ten Brink and others, 2004; Poppe and others, 2006a, b, c, d). These maps help define the geological variability of the sea floor (one of the primary controls of benthic habitat diversity), improve our understanding of the processes that control the distribution and transport of bottom sediments and the distribution of benthic habitats and associated infaunal community structures, and provide a detailed framework for future research, monitoring, and management activities. The bathymetric survey interpreted herein (National Oceanic and Atmospheric Administration (NOAA) survey H11250) covers roughly 94 km² of sea floor in an area where a depression along the Orient Point-Fishers Island segment of the Harbor Hill-Roanoke Point-Charlestown Moraine forms the Race, the eastern opening to Long Island Sound. The Race also divides easternmost Long Island Sound from northwestern Block Island Sound (fig. 1). This bathymetry has been examined in relation to seismic reflection data collected concurrently, as well as archived seismic profiles acquired as part of a long-standing geologic mapping partnership between the State of Connecticut and the U.S. Geological Survey (USGS). The objective of this work was to use these acoustic data sets to interpret geomorphological attributes of the sea floor, and to use these interpretations to better understand the Quaternary geologic history and modern sedimentary processes.

Connecticut;New York↗

Long-term apparent survival of a cold-stunned subpopulation of juveniles green turtles

Understanding the effects of extreme weather on animal populations is fundamental to ecological and conservation sciences and species management. Climate change has resulted in both warm and cold temperature extremes, including an increased frequency of severe cold snaps at middle latitudes in North America. These unusually cold air masses cause rapid declines in nearshore ocean temperatures in coastal areas, with detrimental effects on marine organisms. Acute cold-stun events (hereafter cold stuns) occur when hundreds to thousands of resident juvenile sea turtles fail to escape shallow water during cold snaps. Human intervention through rescue and recovery largely mitigates direct juvenile sea turtle mortality, but delayed effects of cold stuns on rescued individuals are not well understood. Our objective was to examine long-term juvenile green turtle ( Chelonia mydas ) survival across four cold stuns of varying severity in St. Joseph Bay, Florida, between 2010 and 2018. We used the classic Cormack–Jolly–Seber model in a hierarchical Bayesian framework to estimate apparent survival (i.e., emigration and mortality) of rescued turtles at different time intervals. Our results indicated about half of a cohort rescued during a severe cold stun in January 2010 likely remained in the population 1 year later, with 10%–20% remaining 4 years later, and as few as 5% by 2018. The results also suggested higher apparent survival for cohorts rescued during two subsequent milder cold stuns. Emigration was a more plausible ecological explanation for low apparent survival than delayed mortality. Potential ecological mechanisms underlying emigration include a reduction in food availability and a behavioral response to either the severe weather event or handling during rescue (or both). However, the typical annual turnover of juvenile green turtles, though assumed low, is not well known in St. Joseph Bay. Thus, our apparent survival estimates may be reflective of higher-than-expected emigration in the broader population. Our study provides important baseline information about long-term juvenile sea turtle survival after cold stuns in temperate regions. We also highlight the importance of strategic monitoring between cold stuns to examine additional ecological questions.

Florida↗

Global climate changes recorded in coastal wetland sediments: Empirical observations linked to theoretical predictions

Whether coastal areas are experiencing, and responding to, an accelerated rate of global sea-level rise (GSLR) is critically important for the ∼2 billion people living near Earth's oceans. Accretion rates from a suite of physiographically diverse coastal wetlands surrounding Long Island, NY accelerated during the 20th century at 2.3 ± 0.2 × 10 −2 mm yr −2 , which is comparable to reported rates of GSLR acceleration and global temperature changes. Wetlands varied in tidal range, salinity and geomorphic setting, and were located in embayments with limited human impacts in a region with limited and constant rates of subsidence. From geochronologies with temporal resolutions of 2–5 yr, we constructed new composite histories of sediment accretion and mineral deposition. Wetland dynamics are consistent with predictions from sedimentology and a numerical model of ecogeomorphic response, suggesting that these systems, and likely others worldwide, are responding to accelerated GSLR and related climatic changes.

New York↗

An Antarctic ecosystem value index to quantify ecological value across trophic levels and over time

The Southern Ocean around Antarctica is one of the fastest changing regions on the planet and an emerging resource frontier for fisheries. Here, we present the Antarctic Ecosystem Value Index created by merging ecosystem information across food web trophic levels, from phytoplankton to fish and penguins, to quantify the ecological value of marine areas around the Antarctic continent. We find that coastal polynyas - areas of reduced sea-ice - have Index values 31–72% higher than surrounding areas, suggesting that these areas are biologically valuable hot spots for a number of ice-dependent Antarctic Species. Using output from an Earth system model to generate future projections of the Index, we find that high-value locations, often within polynyas, are likely to continue to be valuable throughout the 21st century despite environmental changes. The Antarctic Ecosystem Value Index indicates that penguins lose importance as their habitat becomes increasingly unsuitable, so protecting high-value habitat areas may be critical for these species. This study also shows that while many high-value Index areas are within existing or proposed Marine Protected Areas, there are several opportunities for adopting additional protection, particularly in East Antarctica and the Amundsen Sea.

Nature Communications↗

Potential effects of sea-level rise on the depth to saturated sediments of the Sagamore and Monomoy flow lenses on Cape Cod, Massachusetts

In 2014, the U.S. Geological Survey, in cooperation with the Association to Preserve Cape Cod, the Cape Cod Commission, and the Massachusetts Environmental Trust, began an evaluation of the potential effects of sea-level rise on water table altitudes and depths to water on central and western Cape Cod, Massachusetts. Increases in atmospheric and oceanic temperatures arising, in part, from the release of greenhouse gases likely will result in higher sea levels globally. Increasing water table altitudes in shallow, unconfined coastal aquifer systems could adversely affect infrastructure—roads, utilities, basements, and septic systems—particularly in low-lying urbanized areas. The Sagamore and Monomoy flow lenses on Cape Cod are the largest and most populous of the six flow lenses that comprise the region’s aquifer system, the Cape Cod glacial aquifer. The potential effects of sea-level rise on water table altitude and depths to water were evaluated by use of numerical models of the region. The Sagamore and Monomoy flow lenses have a number of large surface water drainages that receive a substantial amount of groundwater discharge, 47 and 29 percent of the total, respectively. The median increase in the simulated water table altitude following a 6-foot sea-level rise across both flow lenses was 2.11 feet, or 35 percent when expressed as a percentage of the total sea-level rise. The response is nearly the same as the sea-level rise (6 feet) in some coastal areas and less than 0.1 foot near some large inland streams. Median water table responses differ substantially between the Sagamore and Monomoy flow lenses—at 29 and 49 percent, respectively—because larger surface water discharge on the Sagamore flow lens results in increased dampening of the water table response than in the Monomoy flow lens. Surface waters dampen water table altitude increases because streams are fixed-altitude boundaries that cause hydraulic gradients and streamflow to increase as sea-level rises, partially fixing the local water table altitude. The region has a generally thick vadose zone with a mean of about 38 feet; areas with depths to water of 5 feet or less, as estimated from light detection and ranging (lidar) data from 2011 and simulated water table altitudes, currently [2011] occur over about 24.9 square miles, or about 8.4 percent of the total land area of the Sagamore and Monomoy flow lenses, generally in low-lying coastal areas and inland near ponds and streams. Excluding potentially submerged areas, an additional 4.5, 9.8, and 15.9 square miles would have shallow depths to water (5 feet or less) for projected sea-level rises of 2, 4, and 6 feet above levels in 2011. The additional areas with shallow depths to water generally occur in the same areas as the areas with current [2011] depths to water of 5 feet or less: low-lying coastal areas and near inland surface water features. Additional areas with shallow depths to water for the largest sea-level rise prediction (6 feet) account for about 5.7 percent of the total land area, excluding areas likely to be inundated by seawater. The numerous surface water drainages will dampen the response of the water table to sea-level rise. This dampening, combined with the region’s thick vadose zone, likely will mitigate the potential for groundwater inundation in most areas. The potential does exist for groundwater inundation in some areas, but the effects of sea-level rise on depths to water and infrastructure likely will not be substantial on a regional level.

Massachusetts↗

Use of vertical temperature gradients for prediction of tidal flat sediment characteristics

Sediment characteristics largely govern tidal flat morphologic evolution; however, conventional methods of investigating spatial variability in lithology on tidal flats are difficult to employ in these highly dynamic regions. In response, a series of laboratory experiments was designed to investigate the use of temperature diffusion toward sediment characterization. A vertical thermistor array was used to quantify temperature gradients in simulated tidal flat sediments of varying compositions. Thermal conductivity estimates derived from these arrays were similar to measurements from a standard heated needle probe, which substantiates the thermistor methodology. While the thermal diffusivities of dry homogeneous sediments were similar, diffusivities for saturated homogeneous sediments ranged approximately one order of magnitude. The thermal diffusivity of saturated sand was five times the thermal diffusivity of saturated kaolin and more than eight times the thermal diffusivity of saturated bentonite. This suggests that vertical temperature gradients can be used for distinguishing homogeneous saturated sands from homogeneous saturated clays and perhaps even between homogeneous saturated clay types. However, experiments with more realistic tidal flat mixtures were less discriminating. Relationships between thermal diffusivity and percent fines for saturated mixtures varied depending upon clay composition, indicating that clay hydration and/or water content controls thermal gradients. Furthermore, existing models for the bulk conductivity of sediment mixtures were improved only through the use of calibrated estimates of homogeneous end-member conductivity and water content values. Our findings suggest that remotely sensed observations of water content and thermal diffusivity could only be used to qualitatively estimate tidal flat sediment characteristics.

Journal of Geophysical Research C: Oceans↗

Comparing simulated carbon budget of a Lei bamboo forest with flux tower data

Bamboo forest ecosystem is the part of the forest ecosystem. The distribution area of bamboo forest is limited, but in somewhere, like south China, it has been cultivate for a long time with human management. As the climate change has been take great effect on forest carbon budget, many researchers pay attention to the carbon budget in bamboo forest. Moreover cultivative management had a significant impact on the bamboo forest carbon budget. In this study, we modified a terrestrial ecosystem model named Integrated Biosphere Simulator (IBIS) according the management of Lei bamboo forest. Some management, like fertilization, shoots harvesting and organic mulching in winter, had been incorporated into model. Then we had compared model results with the observation data from a Lei bamboo flux tower. The simulated and observed results had achieved good consistency. Our simulated Lei bamboo forest yearly net ecosystem productivity (NEP) was 0.41 kgC a -1 of carbon, which is very close to the observation data 0.45 kgC a -1 of carbon. And the monthly simulated results can take the change of carbon budget in each month, similar to the data we got from flux tower. It reflects that the modified IBIS model can characterize the growth of bamboo forest and perform the simulation well. And then two groups of simulations were set to evaluate effects of cultivative managements on Lei bamboo forests carbon budget. And results showed that both fertilization and organic mulching had taken positive effects on Lei bamboo forests carbon sequestration.

Terrestrial, Atmospheric and Oceanic Sciences↗