Search USGS⌕ Search

SEARCH · Search USGS

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 127 records · Page 7Linked to original sources
USGS2015

Projected wave conditions in the Eastern North Pacific under the influence of two CMIP5 climate scenarios

Hindcast and 21st century winds, simulated by General Circulation Models (GCMs), were used to drive global- and regional-scale spectral wind-wave generation models in the Pacific Ocean Basin to assess future wave conditions along the margins of the North American west coast and Hawaiian Islands. Three-hourly winds simulated by four separate GCMs were used to generate an ensemble of wave conditions for a recent historical time-period (1976–2005) and projections for the mid and latter parts of the 21st century under two radiative forcing scenarios (RCP 4.5 and RCP 8.5), as defined by the fifth phase of the Coupled Model Inter-comparison Project (CMIP5) experiments. Comparisons of results from historical simulations with wave buoy and ERA-Interim wave reanalysis data indicate acceptable model performance of wave heights, periods, and directions, giving credence to generating projections. Mean and extreme wave heights are projected to decrease along much of the North American west coast. Extreme wave heights are projected to decrease south of ∼50°N and increase to the north, whereas extreme wave periods are projected to mostly increase. Incident wave directions associated with extreme wave heights are projected to rotate clockwise at the eastern end of the Aleutian Islands and counterclockwise offshore of Southern California. Local spatial patterns of the changing wave climate are similar under the RCP 4.5 and RCP 8.5 scenarios, but stronger magnitudes of change are projected under RCP 8.5. Findings of this study are similar to previous work using CMIP3 GCMs that indicates decreasing mean and extreme wave conditions in the Eastern North Pacific, but differ from other studies with respect to magnitude and local patterns of change. This study contributes toward a larger ensemble of global and regional climate projections needed to better assess uncertainty of potential future wave climate change, and provides model boundary conditions for assessing the impacts of climate change on coastal systems.

Ocean Modelling↗
USGS1987

Computer simulations of large asteroid impacts into oceanic and continental sites--Preliminary results on atmospheric, cratering and ejecta dynamics

Computer simulations have been completed that describe passage of a 10-km-diameter asteroid through the Earth's atmosphere and the subsequent cratering and ejecta dynamics caused by impact of the asteroid into both oceanic and continental sites. The asteroid was modeled as a spherical body moving vertically at 20 km/s with a kinetic energy of 2.6 × 10 30 ergs (6.2 × 10 7 Mt ). Detailed material modeling of the asteroid, ocean, crustal units, sedimentary unit, and mantle included effects of strength and fracturing, generic asteroid and rock properties, porosity, saturation, lithostatic stresses, and geothermal contributions, each selected to simulate impact and geologic conditions that were as realistic as possible. Calculation of the passage of the asteroid through a U.S. Standard Atmosphere showed development of a strong bow shock wave followed by a highly shock compressed and heated air mass. Rapid expansion of this shocked air created a large low-density region that also expanded away from the impact area. Shock temperatures in air reached ∼20, 000K near the surface of the uplifting crater rim and were as high as ∼2000K at more than 30 km range and 10 km altitude. Calculations to 30 s showed that the shock fronts in the air and in most of the expanding shocked air mass preceded the formation of the crater, ejecta, and rim uplift and did not interact with them. As cratering developed, uplifted rim and target material were ejected into the very low density, shock-heated air immediately above the forming crater, and complex interactions could be expected. Calculations of the impact events showed equally dramatic effects on the oceanic and continental targets through an interval of 120 s. Despite geologic differences in the targets, both cratering events developed comparable dynamic flow fields and by ∼29s had formed similar-sized transient craters ∼39km deep and ∼62km across. Transient-rim uplift of ocean and crust reached a maximum altitude of nearly 40 km at ∼30s and began to decay at velocities of 500 m/s to develop large-tsunami conditions. After ∼30s, strong gravitational rebound drove both craters toward broad flat-floored shapes. At 120 s, transient crater diameters were ∼80km (continental) and ∼105km (oceanic) and transient depths were ∼27km; crater floors consisting of melted and fragmented hot rock were rebounding rapidly upward. By 60 s, the continental crater had ejected ∼2 × 10 14 t, about twice the mass ejected from the oceanic crater. By 120 s, ∼70, 000km 3 (continental) and ∼90, 000km 3 (oceanic) target material were excavated (no mantle) and massive ejecta blankets were formed around the craters. We estimate that in excess of ∼70% of the ejecta would finally lie within ∼3 crater diameters of the impact, and the remaining ejecta (∼10 13 t), including the vaporized asteroid, would be ejected into the atmosphere to altitudes as high as the ionosphere. Effects of secondary volcanism and return of the ocean over hot oceanic crater floor could also be expected to contribute substantial material to the atmosphere.

International Journal of Impact Engineering↗
USGS2009

Bivergent thrust wedges surrounding oceanic island arcs: Insight from observations and sandbox models of the northeastern caribbean plate

At several localities around the world, thrust belts have developed on both sides of oceanic island arcs (e.g., Java-Timor, Panama, Vanuatu, and the northeastern Caribbean). In these localities, the overall vergence of the backarc thrust belt is opposite to that of the forearc thrust belt. For example, in the northeastern Caribbean, a north-verging accretionary prism lies to the north of the Eastern Greater Antilles arc (Hispaniola and Puerto Rico), whereas a south-verging thrust belt called the Muertos thrust belt lies to the south. Researchers have attributed such bivergent geometry to several processes, including: reversal of subduction polarity; subduction-driven mantle flow; stress transmission across the arc; gravitational spreading of the arc; and magmatic inflation within the arc. New observations of deformational features in the Muertos thrust belt and of fault geometries produced in sandbox kinematic models, along with examination of published studies of island arcs, lead to the conclusion that the bivergence of thrusting in island arcs can develop without reversal of subduction polarity, without subarc mantle flow, and without magmatic inflation. We suggest that the Eastern Greater Antilles arc and comparable arcs are simply crustalscale bivergent (or "doubly vergent") thrust wedges formed during unidirectional subduction. Sandbox kinematic modeling suggests, in addition, that a broad retrowedge containing an imbricate fan of thrusts develops only where the arc behaves relatively rigidly. In such cases, the arc acts as a backstop that transmits compressive stress into the backarc region. Further, modeling shows that when arcs behave as rigid blocks, the strike-slip component of oblique convergence is accommodated entirely within the prowedge and the arc-the retrowedge hosts only dip-slip faulting ("frontal thrusting"). The existence of large retrowedges and the distribution of faulting in an island arc may, therefore, be evidence that the arc is relatively rigid. The rigidity of an island arc may arise from its mafi c composition and has implications for seismic-hazard analysis. ?? 2009 Geological Society of America.

Geological Society of America Bulletin↗
USGS2016

Modeling suitable habitat of invasive red lionfish Pterois volitans (Linnaeus, 1758) in North and South America’s coastal waters

We used two common correlative species-distribution models to predict suitable habitat of invasive red lionfish Pterois volitans (Linnaeus, 1758) in the western Atlantic and eastern Pacific Oceans. The Generalized Linear Model (GLM) and the Maximum Entropy (Maxent) model were applied using the Software for Assisted Habitat Modeling. We compared models developed using native occurrences, using non-native occurrences, and using both native and non-native occurrences. Models were trained using occurrence data collected before 2010 and evaluated with occurrence data collected from the invaded range during or after 2010. We considered a total of 22 marine environmental variables. Models built with non-native only or both native and non-native occurrence data outperformed those that used only native occurrences. Evaluation metrics based on the independent test data were highest for models that used both native and non-native occurrences. Bathymetry was the strongest environmental predictor for all models and showed increasing suitability as ocean floor depth decreased, with salinity ranking the second strongest predictor for models that used native and both native and non-native occurrences, indicating low habitat suitability for salinities <30. Our model results also suggest that red lionfish could continue to invade southern latitudes in the western Atlantic Ocean and may establish localized populations in the eastern Pacific Ocean. We reiterate the importance in the choice of the training data source (native, non-native, or native/non-native) used to develop correlative species distribution models for invasive species.

Aquatic Invasions↗
USGS2015

Modeled changes in extreme wave climates of the tropical Pacific over the 21st century: Implications for U.S. and U.S.-Affiliated atoll islands

Wave heights, periods, and directions were forecast for 2081–2100 using output from four coupled atmosphere–ocean global climate models for representative concentration pathway scenarios RCP4.5 and RCP8.5. Global climate model wind fields were used to drive the global WAVEWATCH-III wave model to generate hourly time-series of bulk wave parameters for 25 islands in the mid to western tropical Pacific. December–February 95th percentile extreme significant wave heights under both climate scenarios decreased by 2100 compared to 1976–2010 historical values. Trends under both scenarios were similar, with the higher-emission RCP8.5 scenario displaying a greater decrease in extreme significant wave heights than where emissions are reduced in the RCP4.5 scenario. Central equatorial Pacific Islands displayed the greatest departure from historical values; significant wave heights decreased there by as much as 0.32 m during December–February and associated wave directions rotated approximately 30° clockwise during June–August compared to hindcast data.

Pacific atoll islands↗
USGS2011

Wave constraints for Titan's Jingpo Lacus and Kraken Mare from VIMS specular reflection lightcurves

Stephan et al. (Stephan, K. et al. [2010]. Geophys. Res. Lett. 37, 7104-+.) first saw the glint of sunlight specularly reflected off of Titan's lakes. We develop a quantitative model for analyzing the photometric lightcurve generated during a flyby in which the specularly reflected light flux depends on the fraction of the solar specular footprint that is covered by liquid. We allow for surface waves that spread out the geographic specular intensity distribution. Applying the model to the VIMS T58 observations shows that the waves on Jingpo Lacus must have slopes of no greater than 0.15??, two orders of magnitude flatter than waves on Earth's oceans. Combining the model with theoretical estimates of the intensity of the specular reflection allows a tighter constraint on the waves: <0.05?? Residual specular signal while the specular point lies on land implies that either the land is wetted, the wave slope distribution is non-Gaussian, or that 5% of the land off the southwest edge of Jingpo Lacus is covered in puddles. Another specular sequence off of Kraken Mare acquired during Cassini's T59 flyby shows rapid flux changes that the static model cannot reproduce. Points just 1. min apart vary in flux by more than a factor of two. The present dataset does not uniquely determine the mechanism causing these rapid changes. We suggest that changing wind conditions, kilometer-wavelength waves, or moving clouds could account for the variability. Future specular observations should be designed with a fast cadence, at least 6 points per minute, in order to differentiate between these hypotheses. Such new data will further constrain the nature of Titan's lakes and their interactions with Titan's atmosphere. ?? 2010 Elsevier Inc.

Icarus↗
USGS2017

Integrated geophysical characteristics of the 2015 Illapel, Chile, earthquake

On 16 September 2015, a Mw 8.3 earthquake ruptured the subduction zone offshore of Illapel, Chile, generating an aftershock sequence with 14 Mw 6.0–7.0 events. A double source W phase moment tensor inversion consists of a Mw 7.2 subevent and the main Mw 8.2 phase. We determine two slip models for the mainshock, one using teleseismic broadband waveforms and the other using static GPS and InSAR surface displacements, which indicate high slip north of the epicenter and west-northwest of the epicenter near the oceanic trench. These models and slip distributions published in other studies suggest spatial slip uncertainties of ~25 km and have peak slip values that vary by a factor of 2. We relocate aftershock hypocenters using a Bayesian multiple-event relocation algorithm, revealing a cluster of aftershocks under the Chilean coast associated with deep (20–45 km depth) mainshock slip. Less vigorous aftershock activity also occurred near the trench and along strike of the main aftershock region. Most aftershocks are thrust-faulting events, except for normal-faulting events near the trench. Coulomb failure stress change amplitudes and signs are uncertain for aftershocks collocated with deeper mainshock slip; other aftershocks are more clearly associated with loading from the mainshock. These observations reveal a frictionally heterogeneous interface that ruptured in patches at seismogenic depths (associated with many aftershocks) and with homogeneous slip (and few aftershocks) up to the trench. This event likely triggered seismicity separate from the main slip region, including along-strike events on the megathrust and intraplate extensional events.

Illapel↗
USGS2009

A Tidally Averaged Sediment-Transport Model for San Francisco Bay, California

A tidally averaged sediment-transport model of San Francisco Bay was incorporated into a tidally averaged salinity box model previously developed and calibrated using salinity, a conservative tracer (Uncles and Peterson, 1995; Knowles, 1996). The Bay is represented in the model by 50 segments composed of two layers: one representing the channel (>5-meter depth) and the other the shallows (0- to 5-meter depth). Calculations are made using a daily time step and simulations can be made on the decadal time scale. The sediment-transport model includes an erosion-deposition algorithm, a bed-sediment algorithm, and sediment boundary conditions. Erosion and deposition of bed sediments are calculated explicitly, and suspended sediment is transported by implicitly solving the advection-dispersion equation. The bed-sediment model simulates the increase in bed strength with depth, owing to consolidation of fine sediments that make up San Francisco Bay mud. The model is calibrated to either net sedimentation calculated from bathymetric-change data or measured suspended-sediment concentration. Specified boundary conditions are the tributary fluxes of suspended sediment and suspended-sediment concentration in the Pacific Ocean. Results of model calibration and validation show that the model simulates the trends in suspended-sediment concentration associated with tidal fluctuations, residual velocity, and wind stress well, although the spring neap tidal suspended-sediment concentration variability was consistently underestimated. Model validation also showed poor simulation of seasonal sediment pulses from the Sacramento-San Joaquin River Delta at Point San Pablo because the pulses enter the Bay over only a few days and the fate of the pulses is determined by intra-tidal deposition and resuspension that are not included in this tidally averaged model. The model was calibrated to net-basin sedimentation to calculate budgets of sediment and sediment-associated contaminants. While simulated net sedimentation in the four basins that comprise San Francisco Bay was correct, the simulations incorrectly eroded shallows while channels deposited because model surface-layer boxes span both shallows and channels, and neglect lateral variability of suspended-sediment concentration. Validation with recent (1983-2005) net sedimentation in South San Francisco Bay was poor, perhaps owing to poorly quantified sediment supply, and to invasive species that altered erosion and deposition processes. This demonstrates that deterministically predicting future sedimentation is difficult in this or any estuary for which boundary conditions are not stationary. The model would best be used as a tool for developing past and present sediment budgets, and for creating scenarios of future sedimentation that are compared to one another rather than considered a deterministic prediction.

Scientific Investigations Report↗
USGS1999

Response of North American freshwater lakes to simulated future climates

We apply a physically based lake model to assess the response of North American lakes to future climate conditions as portrayed by the transient trace-gas simulations conducted with the Max Planck Institute (ECHAM4) and the Canadian Climate Center (CGCM1) atmosphere-ocean general circulation models (A/OGCMs). To quantify spatial patterns of lake responses (temperature, mixing, ice cover, evaporation) we ran the lake model for theoretical lakes of specified area, depth, and transparency over a uniformly spaced (50 km) grid. The simulations were conducted for two 10-year periods that represent present climatic conditions and those around the time of CO2 doubling. Although the climate model output produces simulated lake responses that differ in specific regional details, there is broad agreement with regard to the direction and area of change. In particular, lake temperatures are generally warmer in the future as a result of warmer climatic conditions and a substantial loss (> 100 days/yr) of winter ice cover. Simulated summer lake temperatures are higher than 30??C ever the Midwest and south, suggesting the potential for future disturbance of existing aquatic ecosystems. Overall increases in lake evaporation combine with disparate changes in A/OGCM precipitation to produce future changes in net moisture (precipitation minus evaporation) that are of less fidelity than those of lake temperature.

Conference Paper↗
USGS2013

Integrated carbon budget models for the Everglades terrestrial-coastal-oceanic gradient: Current status and needs for inter-site comparisons

Recent studies suggest that coastal ecosystems can bury significantly more C than tropical forests, indicating that continued coastal development and exposure to sea level rise and storms will have global biogeochemical consequences. The Florida Coastal Everglades Long Term Ecological Research (FCE LTER) site provides an excellent subtropical system for examining carbon (C) balance because of its exposure to historical changes in freshwater distribution and sea level rise and its history of significant long-term carbon-cycling studies. FCE LTER scientists used net ecosystem C balance and net ecosystem exchange data to estimate C budgets for riverine mangrove, freshwater marsh, and seagrass meadows, providing insights into the magnitude of C accumulation and lateral aquatic C transport. Rates of net C production in the riverine mangrove forest exceeded those reported for many tropical systems, including terrestrial forests, but there are considerable uncertainties around those estimates due to the high potential for gain and loss of C through aquatic fluxes. C production was approximately balanced between gain and loss in Everglades marshes; however, the contribution of periphyton increases uncertainty in these estimates. Moreover, while the approaches used for these initial estimates were informative, a resolved approach for addressing areas of uncertainty is critically needed for coastal wetland ecosystems. Once resolved, these C balance estimates, in conjunction with an understanding of drivers and key ecosystem feedbacks, can inform cross-system studies of ecosystem response to long-term changes in climate, hydrologic management, and other land use along coastlines

Florida↗
USGS2009

Earth system sensitivity inferred from Pliocene modelling and data

Quantifying the equilibrium response of global temperatures to an increase in atmospheric carbon dioxide concentrations is one of the cornerstones of climate research. Components of the Earths climate system that vary over long timescales, such as ice sheets and vegetation, could have an important effect on this temperature sensitivity, but have often been neglected. Here we use a coupled atmosphere-ocean general circulation model to simulate the climate of the mid-Pliocene warm period (about three million years ago), and analyse the forcings and feedbacks that contributed to the relatively warm temperatures. Furthermore, we compare our simulation with proxy records of mid-Pliocene sea surface temperature. Taking these lines of evidence together, we estimate that the response of the Earth system to elevated atmospheric carbon dioxide concentrations is 30-50% greater than the response based on those fast-adjusting components of the climate system that are used traditionally to estimate climate sensitivity. We conclude that targets for the long-term stabilization of atmospheric greenhouse-gas concentrations aimed at preventing a dangerous human interference with the climate system should take into account this higher sensitivity of the Earth system. ?? 2010 Macmillan Publishers Limited. All rights reserved.

Nature Geoscience↗
USGS2000

Modeling the effects of snowpack on heterotrophic respiration across northern temperate and high latitude regions: Comparison with measurements of atmospheric carbon dioxide in high latitudes

Simulations by global terrestrial biogeochemical models (TBMs) consistently underestimate the concentration of atmospheric carbon dioxide (CO2) at high latitude monitoring stations during the nongrowing season. We hypothesized that heterotrophic respiration is underestimated during the nongrowing season primarily because TBMs do not generally consider the insulative effects of snowpack on soil temperature. To evaluate this hypothesis, we compared the performance of baseline and modified versions of three TBMs in simulating the seasonal cycle of atmospheric CO2 at high latitude CO2 monitoring stations; the modified version maintained soil temperature at 0 ??C when modeled snowpack was present. The three TBMs include the Carnegie-Ames-Stanford Approach (CASA), Century, and the Terrestrial Ecosystem Model (TEM). In comparison with the baseline simulation of each model, the snowpack simulations caused higher releases of CO2 between November and March and greater uptake of CO2 between June and August for latitudes north of 30??N. We coupled the monthly estimates of CO2 exchange, the seasonal carbon dioxide flux fields generated by the HAMOCC3 seasonal ocean carbon cycle model, and fossil fuel source fields derived from standard sources to the three-dimensional atmospheric transport model TM2 forced by observed winds to simulate the seasonal cycle of atmospheric CO2 at each of seven high latitude monitoring stations, in comparison to the CO2 concentrations simulated with the baseline fluxes of each TBM, concentrations simulated using the snowpack fluxes are generally in better agreement with observed concentrations between August and March at each of the monitoring stations. Thus, representation of the insulative effects of snowpack in TBMs generally improves simulation of atmospheric CO2 concentrations in high latitudes during both the late growing season and nongrowing season. These simulations highlight the global importance of biogeochemical processes during the nongrowing season in estimating carbon balance of ecosystems in northern high and temperate latitudes.

Biogeochemistry↗
USGS2002

An upwelling model for the Phosphoria sea: A Permian, ocean-margin sea in the northwest United States

The Permian Phosphoria Formation, a petroleum source rock and world-class phosphate deposit, was deposited in an epicratonic successor basin on the western margin of North America. We calculate the seawater circulation in the basin during deposition of the lower ore zone in the Meade Peak Member from the accumulation rates of carbonate fluorapatite and trace elements. The model gives the exchange rate of water between the Phosphoria sea and the open ocean to the west in terms of an upwelling rate (84 m yr-1) and residence time (4.2 yr) of seawater in the basin. These hydrographic properties supported a mean rate of primary productivity of 0.87 g m-2 d-1 of carbon in the uppermost few tens of meters of the water column (the photic zone) and denitrifying redox conditions in the bottom water (below approximately 150 m depth). High rain rates, onto the sea floor, of the organic matter that hosted the phosphate and several trace elements contributed to the accumulation of phosphorite, chert, and black shales and mudstones. Evaporation in the Goose Egg basin to the east of the Phosphoria basin ensured the import of surface seawater from the Phosphoria sea. Budgets of water, salt, phosphate, and oxygen, plus the minor accumulation of the biomarker gammacerane, show that exchange of water between the two basins was limited, possibly by the shallow carbonate platform that separated the two basins.

Idaho, Montana, Utah, Wyoming↗
USGS2011

Sensitivity of Pliocene ice sheets to orbital forcing

The stability of the Earth's major ice sheets is a critical uncertainty in predictions of future climate and sea level change. One method of investigating the behaviour of the Greenland and the Antarctic ice sheets in a warmer-than-modern climate is to look back at past warm periods of Earth history, for example the Pliocene. This paper presents climate and ice sheet modelling results for the mid-Pliocene warm period (mPWP; 3.3 to 3.0 million years ago), which has been identified as a key interval for understanding warmer-than-modern climates (Jansen et al., 2007). Using boundary conditions supplied by the United States Geological Survey PRISM Group (Pliocene Research, Interpretation and Synoptic Mapping), the Hadley Centre coupled ocean–atmosphere climate model (HadCM3) and the British Antarctic Survey Ice Sheet Model (BASISM), we show large reductions in the Greenland and East Antarctic Ice Sheets (GrIS and EAIS) compared to modern in standard mPWP experiments. We also present the first results illustrating the variability of the ice sheets due to realistic orbital forcing during the mid-Pliocene. While GrIS volumes are lower than modern under even the most extreme (cold) mid-Pliocene orbit (losing at least 35% of its ice mass), the EAIS can both grow and shrink, losing up to 20% or gaining up to 10% of its present-day volume. The changes in ice sheet volume incurred by altering orbital forcing alone means that global sea level can vary by more than 25 m during the mid-Pliocene. However, we have also shown that the response of the ice sheets to mPWP orbital hemispheric forcing can be in anti-phase, whereby the greatest reductions in EAIS volume are concurrent with the smallest reductions of the GrIS. If this anti-phase relationship is in operation throughout the mPWP, then the total eustatic sea level response would be dampened compared to the ice sheet fluctuations that are theoretically possible. This suggests that maximum eustatic sea level rise does not correspond to orbital maxima, but occurs at times where the anti-phasing of Northern and Southern Hemisphere ice sheet retreat is minimised.

Palaeogeography, Palaeoclimatology, Palaeoecology↗
USGS2017

Millennial-scale variability in the local radiocarbon reservoir age of south Florida during the Holocene

A growing body of research suggests that the marine environments of south Florida provide a critical link between the tropical and high-latitude Atlantic. Changes in the characteristics of water masses off south Florida may therefore have important implications for our understanding of climatic and oceanographic variability over a broad spatial scale; however, the sources of variability within this oceanic corridor remain poorly understood. Measurements of ΔR, the local offset of the radiocarbon reservoir age, from shallow-water marine environments can serve as a powerful tracer of water-mass sources that can be used to reconstruct variability in local-to regional-scale oceanography and hydrology. We combined radiocarbon and U-series measurements of Holocene-aged corals from the shallow-water environments of the Florida Keys reef tract (FKRT) with robust statistical modeling to quantify the millennial-scale variability in ΔR at locations with (“nearshore”) and without (“open ocean”) substantial terrestrial influence. Our reconstructions demonstrate that there was significant spatial and temporal variability in ΔR on the FKRT during the Holocene. Whereas ΔR was similar throughout the region after ∼4000 years ago, nearshore ΔR was significantly higher than in the open ocean during the middle Holocene. We suggest that the elevated nearshore ΔR from ∼8000 to 5000 years ago was most likely the result of greater groundwater influence associated with lower sea level at this time. In the open ocean, which would have been isolated from the influence of groundwater, ΔR was lowest ∼7000 years ago, and was highest ∼3000 years ago. We evaluated our open-ocean model of ΔR variability against records of local-to regional-scale oceanography and conclude that local upwelling was not a significant driver of open-ocean radiocarbon variability in this region. Instead, the millennial-scale trends in open-ocean ΔR were more likely a result of broader-scale changes in western Atlantic circulation associated with an increase in the supply of equatorial South Atlantic water to the Caribbean and shifts in the character of South Atlantic waters resulting from variation in the intensity of upwelling off the southwest coast of Africa. Because accurate estimates of ΔR are critical to precise calibrations of radiocarbon dates from marine samples, we also developed models of nearshore and open-ocean ΔR versus conventional 14 C ages that can be used for regional radiocarbon calibrations for the Holocene. Our study provides new insights into the patterns and drivers of oceanographic and hydrologic variability in the Straits of Florida and highlights the value of the paleoceanographic records from south Florida to our understanding of Holocene changes in climate and ocean circulation throughout the Atlantic.

Florida↗
USGS2008

Attribution of declining Western U.S. Snowpack to human effects

Observations show snowpack has declined across much of the western United States over the period 1950-99. This reduction has important social and economic implications, as water retained in the snowpack from winter storms forms an important part of the hydrological cycle and water supply in the region. A formal model-based detection and attribution (D-A) study of these reductions is performed. The detection variable is the ratio of 1 April snow water equivalent (SWE) to water-year-to-date precipitation (P), chosen to reduce the effect of P variability on the results. Estimates of natural internal climate variability are obtained from 1600 years of two control simulations performed with fully coupled ocean-atmosphere climate models. Estimates of the SWE/P response to anthropogenic greenhouse gases, ozone, and some aerosols are taken from multiple-member ensembles of perturbation experiments run with two models. The D-A shows the observations and anthropogenically forced models have greater SWE/P reductions than can be explained by natural internal climate variability alone. Model-estimated effects of changes in solar and volcanic forcing likewise do not explain the SWE/P reductions. The mean model estimate is that about half of the SWE/P reductions observed in the west from 1950 to 1999 are the result of climate changes forced by anthropogenic greenhouse gases, ozone, and aerosols. ?? 2008 American Meteorological Society.

Journal of Climate↗
USGS2015

Projected future vegetation changes for the northwest United States and southwest Canada at a fine spatial resolution using a dynamic global vegetation model.

Future climate change may significantly alter the distributions of many plant taxa. The effects of climate change may be particularly large in mountainous regions where climate can vary significantly with elevation. Understanding potential future vegetation changes in these regions requires methods that can resolve vegetation responses to climate change at fine spatial resolutions. We used LPJ, a dynamic global vegetation model, to assess potential future vegetation changes for a large topographically complex area of the northwest United States and southwest Canada (38.0&ndash;58.0&deg;N latitude by 136.6&ndash;103.0&deg;W longitude). LPJ is a process-based vegetation model that mechanistically simulates the effect of changing climate and atmospheric CO 2 concentrations on vegetation. It was developed and has been mostly applied at spatial resolutions of 10-minutes or coarser. In this study, we used LPJ at a 30-second (~1-km) spatial resolution to simulate potential vegetation changes for 2070&ndash;2099. LPJ was run using downscaled future climate simulations from five coupled atmosphere-ocean general circulation models (CCSM3, CGCM3.1(T47), GISS-ER, MIROC3.2(medres), UKMO-HadCM3) produced using the A2 greenhouse gases emissions scenario. Under projected future climate and atmospheric CO 2 concentrations, the simulated vegetation changes result in the contraction of alpine, shrub-steppe, and xeric shrub vegetation across the study area and the expansion of woodland and forest vegetation. Large areas of maritime cool forest and cold forest are simulated to persist under projected future conditions. The fine spatial-scale vegetation simulations resolve patterns of vegetation change that are not visible at coarser resolutions and these fine-scale patterns are particularly important for understanding potential future vegetation changes in topographically complex areas.

PLoS ONE↗