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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 739 records · Page 41Linked to original sources

USGS1988

Tectonic evolution of Gorda Ridge inferred from sidescan sonar images

Gorda Ridge is the southern segment of the Juan de Fuca Ridge complex, in the north-east Pacific. Along-strike spreading-rate variation on Gorda Ridge and deformation of Gorda Plate are evidence for compression between the Pacific and Gorda Plates. GLORIA sidescan sonographs allow the spreading fabric associated with Gorda Ridge to be mapped in detail. Between 5 and 2 Ma, a pair of propagating rifts re-orientated the northern segment of Gorda Ridge by about 10?? clockwise, accommodating a clockwise shift in Pacific-Juan de Fuca plate motion that occurred around 5 Ma. Deformation of Gorda Plate, associated with southward decreasing spreading rates along southern Gorda Ridge, is accommodated by a combination of clockwise rotation of Gorda Plate crust, coupled with left-lateral motion on the original normal faults of the ocean crust. Segments of Gorda Plate which have rotated by different amounts are separated by narrow deformation zones across which sharp changes in ocean fabric trend are seen. Although minor lateral movement may occur on these NW to WNW structures, no major right-lateral movement, as predicted by previous models, is observed. ?? 1988 Kluwer Academic Publishers.

Marine Geophysical Research↗

USGS2003

Future coral reef habitat marginality: Temporal and spatial effects of climate change in the Pacific basin

Marginal reef habitats are regarded as regions where coral reefs and coral communities reflect the effects of steady-state or long-term average environmental limitations. We used classifications based on this concept with predicted time-variant conditions of future climate to develop a scenario for the evolution of future marginality. Model results based on a conservative scenario of atmospheric CO2 increase were used to examine changes in sea surface temperature and aragonite saturation state over the Pacific Ocean basin until 2069. Results of the projections indicated that essentially all reef locations are likely to become marginal with respect to aragonite saturation state. Significant areas, including some with the highest biodiversity, are expected to experience high-temperature regimes that may be marginal, and additional areas will enter the borderline high temperature range that have experienced significant ENSO-related bleaching in the recent past. The positive effects of warming in areas that are presently marginal in terms of low temperature were limited. Conditions of the late 21st century do not lie outside the ranges in which present-day marginal reef systems occur. Adaptive and acclimative capabilities of organisms and communities will be critical in determining the future of coral reef ecosystems.

Conference Paper↗

USGS2006

Interpreting the spatio-temporal patterns of sea turtle strandings: Going with the flow

Knowledge of the spatial and temporal distribution of specific mortality sources is crucial for management of species that are vulnerable to human interactions. Beachcast carcasses represent an unknown fraction of at-sea mortalities. While a variety of physical (e.g., water temperature) and biological (e.g., decomposition) factors as well as the distribution of animals and their mortality sources likely affect the probability of carcass stranding, physical oceanography plays a major role in where and when carcasses strand. Here, we evaluate the influence of nearshore physical oceanographic and wind regimes on sea turtle strandings to decipher seasonal trends and make qualitative predictions about stranding patterns along oceanfront beaches. We use results from oceanic drift-bottle experiments to check our predictions and provide an upper limit on stranding proportions. We compare predicted current regimes from a 3D physical oceanographic model to spatial and temporal locations of both sea turtle carcass strandings and drift bottle landfalls. Drift bottle return rates suggest an upper limit for the proportion of sea turtle carcasses that strand (about 20%). In the South Atlantic Bight, seasonal development of along-shelf flow coincides with increased numbers of strandings of both turtles and drift bottles in late spring and early summer. The model also predicts net offshore flow of surface waters during winter - the season with the fewest relative strandings. The drift bottle data provide a reasonable upper bound on how likely carcasses are to reach land from points offshore and bound the general timeframe for stranding post-mortem (< two weeks). Our findings suggest that marine turtle strandings follow a seasonal regime predictable from physical oceanography and mimicked by drift bottle experiments. Managers can use these findings to reevaluate incidental strandings limits and fishery takes for both nearshore and offshore mortality sources. ?? 2005 Elsevier Ltd. All rights reserved.

Biological Conservation↗

USGS2020

Mysterious tsunami in the Caribbean Sea following the 2010 Haiti earthquake possibly generated by dynamically triggered early aftershocks

Dynamically triggered offshore aftershocks, caused by passing seismic waves from main shocks located on land, are currently not considered in tsunami warnings. The M7.0 2010 Haiti earthquake epicenter was located on land 27 km north of the Caribbean Sea and its focal mechanism was oblique strike-slip. Nevertheless, a tsunami recorded on a Caribbean Deep-Ocean Assessment and Reporting of Tsunami (DART) buoy and a tide gauge produced runup heights of 1–3 m along Haiti southeast coast. Earthquake finite-fault model inversions of the DART waveform suggest that a reverse fault doublet with magnitudes of M6.8 and M6.5 located 85 km southwest of the epicenter may have excited the tsunami. This doublet collocates with dynamically triggered aftershocks, derived from back-projection analysis, that occurred 20-60 s after the main shock of the Haiti earthquake. The aftershocks are within a region of maximum dynamic strain predicted by the main shock, on a possibly tectonically active submarine ridge southwest of Haiti's Southern Peninsula. The agreement between the tsunami finite-fault source models and the seismic and tectonic evidence suggests that earthquakes on land, even strike-slip faults, can generate tsunamis by dynamically triggering offshore aftershocks.

Earth and Planetary Science Letters↗

USGS2007

Generation and propagation of nonlinear internal waves in Massachusetts Bay

During the summer, nonlinear internal waves (NLIWs) are commonly observed propagating in Massachusetts Bay. The topography of the area is unique in the sense that the generation area (over Stellwagen Bank) is only 25 km away from the shoaling area, and thus it represents an excellent natural laboratory to study the life cycle of NLIWs. To assist in the interpretation of the data collected during the 1998 Massachusetts Bay Internal Wave Experiment (MBIWE98), a fully nonlinear and nonhydrostatic model covering the generation/shoaling region was developed, to investigate the response of the system to the range of background and driving conditions observed. Simplified models were also used to elucidate the role of nonlinearity and dispersion in shaping the NLIW field. This paper concentrates on the generation process and the subsequent evolution in the basin. The model was found to reproduce well the range of propagation characteristics observed (arrival time, propagation speed, amplitude), and provided a coherent framework to interpret the observations. Comparison with a fully nonlinear hydrostatic model shows that during the generation and initial evolution of the waves as they move away from Stellwagen Bank, dispersive effects play a negligible role. Thus the problem can be well understood considering the geometry of the characteristics along which the Riemann invariants of the hydrostatic problem propagate. Dispersion plays a role only during the evolution of the undular bore in the middle of Stellwagen Basin. The consequences for modeling NLIWs within hydrostatic models are briefly discussed at the end.

Massachusetts↗

USGS2004

One-dimensional wave bottom boundary layer model comparison: Specific eddy viscosity and turbulence closure models

Six one-dimensional-vertical wave bottom boundary layer models are analyzed based on different methods for estimating the turbulent eddy viscosity: Laminar, linear, parabolic, k —one equation turbulence closure, k−ε —two equation turbulence closure, and k−ω —two equation turbulence closure. Resultant velocity profiles, bed shear stresses, and turbulent kinetic energy are compared to laboratory data of oscillatory flow over smooth and rough beds. Bed shear stress estimates for the smooth bed case were most closely predicted by the k−ω model. Normalized errors between model predictions and measurements of velocity profiles over the entire computational domain collected at 15° intervals for one-half a wave cycle show that overall the linear model was most accurate. The least accurate were the laminar and k−ε models. Normalized errors between model predictions and turbulence kinetic energy profiles showed that the k−ω model was most accurate. Based on these findings, when the smallest overall velocity profile prediction error is required, the processing requirements and error analysis suggest that the linear eddy viscosity model is adequate. However, if accurate estimates of bed shear stress and TKE are required then, of the models tested, the k−ω model should be used.

Journal of Waterway, Port, Coastal and Ocean Engin↗

USGS2004

Evidence for increased latent heat transport during the Cretaceous (Albian) greenhouse warming

Quantitative estimates of increased heat transfer by atmospheric H 2O vapor during the Albian greenhouse warming suggest that the intensified hydrologic cycle played a greater role in warming high latitudes than at present and thus represents a viable alternative to oceanic heat transport. Sphaerosiderite ??18O values in paleosols of the North American Cretaceous Western Interior Basin are a proxy for meteoric ??18O values, and mass-balance modeling results suggest that Albian precipitation rates exceeded modern rates at both mid and high latitudes. Comparison of modeled Albian and modern precipitation minus evaporation values suggests amplification of the Albian moisture deficit in the tropics and moisture surplus in the mid to high latitudes. The tropical moisture deficit represents an average heat loss of ???75 W/m2 at 10??N paleolatitude (at present, 21 W/m2). The increased precipitation at higher latitudes implies an average heat gain of ???83 W/m2 at 45??N (at present, 23 W/m2) and of 19 W/m2 at 75??N (at present, 4 W/m2). These estimates of increased poleward heat transfer by H2O vapor during the Albian may help to explain the reduced equator-to-pole temperature gradients. ?? 2004 Geological Society of America.

Geology↗

USGS2025

Simulated hydrologic responses to proposed wastewater-returnflow scenarios in Falmouth, Massachusetts

The Cape Cod aquifer is the sole source of drinking water for communities on Cape Cod, Massachusetts, including the Town of Falmouth, where the aquifer is currently threatened by contamination from septic-system-derived nitrogen. To address this problem, the Town is proposing to sewer areas of Falmouth, treat the wastewater at the Town’s Main Wastewater Treatment Facility (a nitrogen removing/tertiary treatment facility), and discharge the treated wastewater to an ocean outfall pipe in Nantucket Sound. The U.S. Geological Survey, in cooperation with the Town of Falmouth, updated a three-dimensional steady-state groundwater flow model to represent current (defined as 2019–23) average hydrologic conditions and to simulate the long-term average freshwater hydrologic response to two wastewater-return-flow scenarios. Scenario 1 involves the sewering of all properties south of Route 28 in Falmouth, which approximates the Town’s possible sewer expansion over the next 20–30 years. Scenario 2 involves sewering of all properties in Falmouth to demonstrate the maximum potential effect of sewering on the aquifer. Overall, the simulated hydrologic response of water-table altitudes and streamflow in both scenarios was relatively small compared to fluctuations from natural recharge. In scenario 1, the water-table altitude decreased by about 0.1 feet south of Route 28, where the conversion to municipal sewers removed wastewater-return flow from onsite septic systems. The water-table altitude decreased by about 0.1–0.2 feet over a larger area in Falmouth under town-wide sewering in scenario 2. The greatest decrease in water-table altitude in both scenarios occurred near the Main Wastewater Treatment Facility, with a decrease of about 1.1 feet in scenario 1 and about 1.3 feet in scenario 2. Simulated decreases in streamflow also were estimated for six selected streams in Falmouth and Mashpee. In both scenarios, the largest simulated decreases in streamflow were at the Coonamessett River, which is the closest stream to the Main Wastewater Treatment Facility. In scenario 1, the average annual decrease in flow at the Coonamessett River was 0.1 cubic feet per second, a 1.1 percent decrease from current (2019–23) conditions. In scenario 2, streamflow at the Coonamessett River decreased by 0.6 cubic feet per second, a 5.4 percent decrease from current (2019–23) conditions.

Massachusetts↗

USGS2003

A delayed action oscillator shared by biennial, interannual, and decadal signals in the Pacific Basin

Biennial, interannual, and decadal signals in the Pacific basin are observed to share patterns and evolution in covarying sea surface temperature (SST), 18??C isotherm depth (Z18), zonal surface wind (ZSW), and wind stress curl (WSC) anomalies from 1955 to 1999. Each signal has warm SST anomalies propagating slowly eastward along the equator, generating westerly ZSW anomalies in their wake. These westerly ZSW anomalies produce cyclonic WSC anomalies off the equator which pump baroclinic Rossby waves in the western/central tropical North Pacific Ocean. These Rossby waves propagate westward, taking ???6, ???12, and ???36 months to reach the western boundary near ???7??N, ???12??N, and ???18??N on biennial, interannual, and decadal period scales, respectively. There, they reflect as equatorial coupled waves, propagating slowly eastward in covarying SST, Z18, and ZSW anomalies, taking ???6, ???12, and ???24 months to reach the central/eastern equatorial ocean. These equatorial coupled waves produce a delayed-negative feedback to the warm SST anomalies there. The decrease in Rossby wave phase speed with latitude, the increase in meridional scale of equatorial SST anomalies with period scale, and the associated increase in latitude of Rossby wave forcing are consistent with the delayed action oscillator (DAO) model used to explain El Nin??o. However, this is not true of the western-boundary reflection of Rossby waves into slow equatorial coupled waves. This requires modification of the extant DAO model. We construct a modified DAO model, demonstrating how the various mechanisms and the size and sources of their delays yield the resulting frequency of each signal.

Journal of Geophysical Research C: Oceans↗

USGS2016

Simulated effects of groundwater withdrawals from aquifers in Ocean County and vicinity, New Jersey

Rapid population growth since the 1930s in Ocean County and vicinity, New Jersey, has placed increasing demands upon the area’s freshwater resources. To examine effects of groundwater withdrawals, a three-dimensional groundwater-flow model was developed to simulate the groundwater-flow systems of five area aquifers: the unconfined Kirkwood-Cohansey aquifer system and Vincentown aquifer, and three confined aquifers— the Rio Grande water-bearing zone, the Atlantic City 800-foot sand, and the Piney Point aquifer. The influence of withdrawals is evaluated by using transient groundwater-flow model simulations that incorporate three withdrawal schemes. These are (1) no-withdrawal conditions; (2) 2000–03 withdrawal conditions, using reported monthly withdrawals at all production wells from January 2000 through December 2003; and (3) maximum-allocation withdrawal conditions using the maximum withdrawal allowed by New Jersey Department of Environmental Protection permits at each well. Particle tracking analysis, using results from model simulations, delineated particle flow paths from production wells to the point of recharge, and estimated particle travel times. Compared with no-withdrawal conditions, 2000–03 withdrawal conditions reduced the amount of groundwater flow out of the Kirkwood-Cohansey aquifer system into streams, increased the net flow of water into other layers, reduced net flow into or out of storage, and reduced flow from the Kirkwood-Cohansey aquifer system to constant head cells. Freshwater discharging to the Barnegat Bay-Little Egg Harbor estuary from streams and groundwater is essential to maintaining the ecology of the bay. Examination of selected stress periods indicates that simulated base flow in streams flowing into the Barnegat Bay-Little Egg Harbor estuary is reduced by as much as 49 cubic feet per second for 2000 to 2003 withdrawal conditions when compared with no-withdrawal conditions. In the three confined aquifers, water levels during periods of low recharge and high withdrawals, and high recharge and low withdrawals, were examined to determine seasonal effects on the confined flow systems. The simulated potentiometric surface of the Rio Grande water-bearing zone and the Atlantic City 800-foot sand during selected stress periods indicates substantial declines from no-withdrawal conditions to 2000–03 conditions as a result of groundwater withdrawals. Cones of depression in Toms River Township, Seaside Heights and Seaside Park Boroughs, and Barnegat Light Borough developed in the potentiometric surface of the Piney Point aquifer in response to withdrawals. Maximum-allocation withdrawals decreased flow out of the Kirkwood-Cohansey aquifer system to constant head cells, increased flow out of the aquifer system to adjacent and lower layers, and reduced groundwater discharge to streams when compared with 2000–03 withdrawal conditions. Increases in withdrawals from the Rio Grande water-bearing zone, the Atlantic City 800-foot sand, and the Piney Point aquifer result in an increase in simulated net groundwater flow into these aquifers. Base-flow reduction from 2000–03 conditions to maximum-allocation conditions of 25 to 29 cubic feet per second in all streams draining to the Barnegat Bay-Little Egg Harbor also is indicated. Potentiometric surfaces of the Rio Grande water-bearing zone, Atlantic City 800-foot sand, and the Piney Point aquifer during two stress periods of simulated maximum-allocation withdrawal conditions indicated the expansion of several cones of depression developed during 2000–03 withdrawals. Simulation of average 2000–03 withdrawal conditions indicated the extent to which the groundwater-flow system is susceptible to potential saltwater intrusion into near-shore wells. Travel time from recharge to discharge location ranged from 11 to approximately 50,700 years in near-shore Kirkwood-Cohansey aquifer system wells. Those in Seaside Heights Borough, in Island Beach State Park (Berkeley Township), and in Ship Bottom Borough have particle travel times from 140 to 12,000 years and flow paths that originated under Barnegat Bay or the Atlantic Ocean from the simulation of average maximum-allocation withdrawal conditions. Travel time along flow paths to wells screened in the Rio Grande water-bearing zone and the Atlantic City 800-foot sand from recharge to discharge point ranged from nearly 530 years to greater than 3.73 million years from the simulation of average 2000–03 withdrawal conditions. Particle tracking indicated that most wells screened in these aquifers derived a large part of their recharge from the Oswego River Basin, with a small portion of flow originating either beneath Barnegat Bay or to the east beneath the Atlantic Ocean. Travel time along flow paths that start beneath either Barnegat Bay or the Atlantic Ocean ranged from 2,300 to approximately 134,000 years from the simulation of average maximum-allocation withdrawal conditions."

New Jersey↗

USGS2014

Geoelectrical structure of the central zone of Piton de la Fournaise volcano (Reunion)

A study of the geoelectrical structure of the central part of Piton de la Fournaise volcano (Reunion, Indian Ocean) was made using direct current electrical (DC) and transient electromagnetic soundings (TEM). Piton de la Fournaise is a highly active oceanic basaltic shield and has been active for more than half a million years. Joint interpretation of the DC and TEM data allows us to obtain reliable 1D models of the resistivity distribution. The depth of investigation is of the order of 1.5 km but varies with the resistivity pattern encountered at each sounding. Two-dimensional resistivity cross sections were constructed by interpolation between the soundings of the 1D interpreted models. Conductors with resistivities less than 100 ohm-m are present at depth beneath all of the soundings and are located high in the volcanic edifice at elevations between 2000 and 1200 m. The deepest conductor has a resistivity less than 20 ohm-m for soundings located inside the Enclos and less than 60-100 ohm-m for soundings outside the Enclos. From the resistivity distributions, two zones are distinguished: (a) the central zone of the Enclos; and (b) the outer zone beyond the Enclos. Beneath the highly active summit area, the conductor rises to within a few hundred meters of the surface. This bulge coincides with a 2000-mV self-potential anomaly. Low-resistivity zones are inferred to show the presence of a hydrothermal system where alteration by steam and hot water has lowered the resistivity of the rocks. Farther from the summit, but inside the Enclos the depth to the conductive layers increases to approximately 1 km and is inferred to be a deepening of the hydrothermally altered zone. Outside of the Enclos, the nature of the deep, conductive layers is not established. The observed resistivities suggest the presence of hydrated minerals, which could be found in landslide breccias, in hydrothermally altered zones, or in thick pyroclastic layers. Such formations often create perched water tables. The known occurrence of large eastward-moving landslides in the evolution of Piton de la Fournaise strongly suggests that large volumes of breccias should exist in the interior of the volcano; however, extensive breccia deposits are not observed at the bottom of the deep valleys that incise the volcano to elevations lower than those determined for the top of the conductors. The presence of the center of Piton de la Fournaise beneath the Plaine des Sables area during earlier volcanic stages (ca. 0.5 to 0.150 Ma) may have resulted in broad hydrothermal alteration of this zone. However, this interpretation cannot account for the low resistivities in peripheral zones. It is not presently possible to discriminate between these general interpretations. In addition, the nature of the deep conductors may be different in each zone. Whatever the geologic nature of these conductive layers, their presence indicates a major change of lithology at depth, unexpected for a shield volcano such as Piton de la Fournaise.

Bulletin of Volcanology↗

USGS2006

Airborne microorganisms in the African desert dust corridor over the mid-Atlantic ridge, Ocean Drilling Program, Leg 209

The objective of this study was to enhance our understanding of the fate and trans-Atlantic transport of dustborne microorganisms from Northern Africa to the Caribbean and Americas, and more specifically to determine if culturable populations could be detected at a mid-ocean site, closer to the source of dust relative to land-based Caribbean sites, during the early summer months of May and June. Between the dates of 22 May and 30 June 2003, daily air samples were collected and evaluated for the presence of culturable bacterial and fungal colony-forming units (CFU). Here we report a statistically significant correlation between daily atmospheric CFU counts at a mid-ocean research site (???15??N, 45??W) and daily desert dust concentrations as determined by the U.S. Navy's Naval Aerosol Analysis and Prediction System (NAAPS) Global Aerosol Model (Honrath et al. (2004). Journal of Geophysical Research, 109; Johnson et al. (2003). Global Biogeochemical Cycles, 17, 1063; Reid et al. (2004). Geophysical Research Letters, 31; Schollaert, Yoder, Westphal, & O'Reilly (2003). Journal of Geophysical Research, 108, 3191). ?? Springer Science+Business Media B.V. 2006.

Aerobiologia↗

USGS2025

Reconciliation of geochronology and paleozoogeography for Quaternary marine terraces, San Luis Obispo Bay area, California, USA

In the San Luis Obispo Bay area of central California, interpretations of marine terrace ages have been hampered by inconsistent results from geochronological indicators (U-series ages of corals and correlations using amino acid racemization of mollusks) and seemingly contradictory paleozoogeographic aspects of fossil faunal assemblages. New U-series ages of corals, amino acid data from mollusks, and detailed analyses of fossil mollusk assemblages are presented that reconcile many of the apparently discordant results from previous studies. The two lowest-elevation terraces are Q1 (lower) and Q2 (upper); both are thought to date from high-sea stands of marine isotope stage (MIS) 5, the last interglacial complex. A combination of U-series dating and amino acid racemization results indicates that the Q1 terrace probably dates to the ∼80 ka high-sea stand of MIS 5a. U-series analyses of corals from the Q2 terrace show open-system histories, but consideration of two alternative open-system histories indicates that this terrace likely hosts corals dating to the high-sea stands of MIS 5c (∼105–95 ka) and MIS 5e (∼130–115 ka). Amino acid ratios in the fossil bivalve Leukoma staminea (Conrad, 1837) support the age differences between the two terraces and the open-system models suggested by the U-series data. Alkenone studies along much of the coast of California have shown that sea surface temperatures (SSTs) of the eastern Pacific Ocean were warmer than present during MIS 5e. During MIS 5c, on the other hand, SSTs were either similar to the present or cooler, but in either case, much cooler than during MIS 5e. The paleozoogeographic aspects of the molluscan faunas from the Q2 terrace support an interpretation that these deposits contain a mixture of fossils from both the MIS 5c and 5e high-sea stands. Fossils from the Q2 terrace include southern or southward-ranging species (likely dating to MIS 5e) and northern or northward-ranging species (likely dating to MIS 5c). A mixture of MIS 5e and MIS 5c fossils on the Q2 terrace is consistent with similar data reported for other localities in California. This phenomenon is explained by glacial isostatic adjustment (GIA) processes, which brought about a higher-than-present sea level during MIS 5c. GIA processes, combined with a low uplift rate, allowed mixing of MIS 5e fossils during the MIS 5c high-sea stand through capture of an MIS 5e terrace that had experienced minimal uplift in ∼20,000 years. Terrace reoccupation and fossil mixing such as that in the San Luis Obispo Bay area can be expected along any part of the Pacific Coast of North America that was affected by GIA processes and where uplift rates are low.

California↗

USGS2004

The effects of mixotrophy on the stability and dynamics of a simple planktonic food web

Recognition of the microbial loop as an important part of aquatic ecosystems disrupted the notion of simple linear food chains. However, current research suggests that even the microbial loop paradigm is a gross simplification of microbial interactions due to the presence of mixotrophs—organisms that both photosynthesize and graze. We present a simple food web model with four trophic species, three of them arranged in a food chain (nutrients–autotrophs–herbivores) and the fourth as a mixotroph with links to both the nutrients and the autotrophs. This model is used to study the general implications of inclusion of the mixotrophic link in microbial food webs and the specific predictions for a parameterization that describes open ocean mixed layer plankton dynamics. The analysis indicates that the system parameters reside in a region of the parameter space where the dynamics converge to a stable equilibrium rather than displaying periodic or chaotic solutions. However, convergence requires weeks to months, suggesting that the system would never reach equilibrium in the ocean due to alteration of the physical forcing regime. Most importantly, the mixotrophic grazing link seems to stabilize the system in this region of the parameter space, particularly when nutrient recycling feedback loops are included.

Theoretical Population Biology↗

USGS2014

Termini of calving glaciers as self-organized critical systems

Over the next century, one of the largest contributions to sea level rise will come from ice sheets and glaciers calving ice into the ocean 1 . Factors controlling the rapid and nonlinear variations in calving fluxes are poorly understood, and therefore difficult to include in prognostic climate-forced land-ice models. Here we analyse globally distributed calving data sets from Svalbard, Alaska (USA), Greenland and Antarctica in combination with simulations from a first-principles, particle-based numerical calving model to investigate the size and inter-event time of calving events. We find that calving events triggered by the brittle fracture of glacier ice are governed by the same power-law distributions as avalanches in the canonical Abelian sandpile model 2 . This similarity suggests that calving termini behave as self-organized critical systems that readily flip between states of sub-critical advance and super-critical retreat in response to changes in climate and geometric conditions. Observations of sudden ice-shelf collapse and tidewater glacier retreat in response to gradual warming of their environment 3 are consistent with a system fluctuating around its critical point in response to changing external forcing. We propose that self-organized criticality provides a yet unexplored framework for investigations into calving and projections of sea level rise.

Nature Geoscience↗

USGS2010

Thermochronology and tectonics of the Leeward Antilles: Evolution of the southern Caribbean Plate boundary zone

Tectonic reconstructions of the Caribbean Plate are severely hampered by a paucity of geochronologic and exhumation constraints from anastomosed basement blocks along its southern margin. New U/Pb, 40 Ar/ 39 Ar, apatite fission track, and apatite (U-Th)/He data constrain quantitative thermal and exhumation histories, which have been used to propose a model for the tectonic evolution of the emergent parts of the Bonaire Block and the southern Caribbean Plate boundary zone. An east facing arc system intruded through an oceanic plateau during ~90 to ~87 Ma and crops out on Aruba. Subsequent structural displacements resulted in >80°C of cooling on Aruba during 70–60 Ma. In contrast, exhumation of the island arc sequence exposed on Bonaire occurred at 85–80 Ma and 55–45 Ma. Santonian exhumation on Bonaire occurred immediately subsequent to burial metamorphism and may have been driven by the collision of a west facing island arc with the Caribbean Plate. Island arc rocks intruded oceanic plateau rocks on Gran Roque at ~65 Ma and exhumed rapidly at 55–45 Ma. We attribute Maastrichtian-Danian exhumation on Aruba and early Eocene exhumation on Bonaire and Gran Roque to sequential diachronous accretion of their basement units to the South American Plate. Widespread unconformities indicate late Eocene subaerial exposure. Late Oligocene–early Miocene dextral transtension within the Bonaire Block drove subsidence and burial of crystalline basement rocks of the Leeward Antilles to ≤1 km. Late Miocene–recent transpression caused inversion and ≤1 km of exhumation, possibly as a result of the northward escape of the Maracaibo Block.

Tectonics↗

USGS2011

Natural radium and radon tracers to quantify water exchange and movement in reservoirs

Radon and radium isotopes are routinely used to quantify exchange rates between different hydrologic reservoirs. Since their recognition as oceanic tracers in the 1960s, both radon and radium have been used to examine processes such as air-sea exchange, deep oceanic mixing, benthic inputs, and many others. Recently, the application of radon-222 and the radium-quartet (223,224,226,228Ra) as coastal tracers has seen a revelation with the growing interest in coastal groundwater dynamics. The enrichment of these isotopes in benthic fluids including groundwater makes both radium and radon ideal tracers of coastal benthic processes (e.g. submarine groundwater discharge). In this chapter we review traditional and recent advances in the application of radon and radium isotopes to understand mixing and exchange between various hydrologic reservoirs, specifically: (1) atmosphere and ocean, (2) deep and shallow oceanic water masses, (3) coastal groundwater/benthic pore waters and surface ocean, and (4) aquifer-lakes. While the isotopes themselves and their distribution in the environment provide qualitative information about the exchange processes, it is mixing/exchange and transport models for these isotopes that provide specific quantitative information about these processes. Brief introductions of these models and mixing parameters are provided for both historical and more recent studies.

Book chapter↗

USGS2006

Photometric properties of Titan's surface from Cassini VIMS: Relevance to titan's hemispherical albedo dichotomy and surface stability

The Visual and Infrared Mapping Spectrometer (VIMS) instrument on the Cassini Saturn Orbiter returned spectral imaging data as the spacecraft undertook six close encounters with Titan beginning 7 July, 2004. Three of these flybys each produced overlapping coverage of two distinct regions of Titan's surface. Twenty-four points were selected on approximately opposite hemispheres to serve as photometric controls. Six points were selected in each of four reflectance classes. On one hemisphere each control point was observed at three distinct phase angles. From the derived phase coefficients, preliminary normal reflectances were derived for each reflectance class. The normal reflectance of Titan's surface units at 2.0178 ??m ranged from 0.079 to 0.185 for the most absorbing to the most reflective units assuming no contribution from absorbing haze. When a modest haze contribution of ??=0.1 is considered these numbers increase to 0.089-0.215. We find that the lowest three reflectance classes have comparable normal reflectance on either hemisphere. However, for the highest brightness class the normal reflectance is higher on the hemisphere encompassing longitude 14-65?? compared to the same high brightness class for the hemisphere encompassing 122-156?? longitude. We conclude that an albedo dichotomy observed in continental sized units on Titan is due not only to one unit having more areal coverage of reflective material than the other but the material on the brighter unit is intrinsically more reflective than the most reflective material on the other unit. This suggests that surface renewal processes are more widespread on Titan's more reflective units than on its less reflective units. We note that one of our photometric control points has increased in reflectance by 12% relative to the surrounding terrain from July of 2004 to April and May of 2005. Possible causes of this effect include atmospheric processes such as ground fog or orographic clouds; the suggestion of active volcanism cannot be ruled out. Several interesting circular features which resembled impact craters were identified on Titan's surface at the time of the initial Titan flyby in July of 2004. We traced photometric profiles through two of these candidate craters and attempted to fit these profiles to the photometric properties expected from model depressions. We find that the best-fit attempt to model these features as craters requires that they be unrealistically deep, approximately 70 km deep. We conclude that despite their appearance, these circular features are not craters, however, the possibility that they are palimpsests cannot be ruled out. We used two methods to test for the presence of vast expanses of liquids on Titan's surface that had been suggested to resemble oceans. Specular reflection of sunlight would be indicative of widespread liquids on the surface; we found no evidence of this. A large liquid body should also show uniformity in photometric profile; we found the profiles to be highly variable. The lack of specular reflection and the high photometric variability in the profiles across candidate oceans is inconsistent with the presence of vast expanses of flat-lying liquids on Titan's surface. While liquid accumulation may be present as small, sub-pixel-sized bodies, or in areas of the surface which still remain to be observed by VIMS, the presence of large ocean-sized accumulations of liquids can be ruled out. The Cassini orbital tour offers the opportunity for VIMS to image the same parts of Titan's surface repeatedly at many different illumination and observation geometries. This creates the possibility of understanding the properties of Titan's atmosphere and haze by iteratively adapting models to create a best fit to the surface reflectance properties. ?? 2006 Elsevier Ltd. All rights reserved.

Planetary and Space Science↗