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 91 records · Page 5Linked to original sources
USGS2011

Bathymetric controls on Pliocene North Atlantic and Arctic sea surface temperature and deepwater production

The mid-Pliocene warm period (MPWP; ~ 3.3 to 3.0 Ma) is the most recent interval in Earth's history in which global temperatures reached and remained at levels similar to those projected for the near future. The distribution of global warmth, however, was different than today in that the high latitudes warmed more than the tropics. Multiple temperature proxies indicate significant sea surface warming in the North Atlantic and Arctic Oceans during the MPWP, but predictions from a fully coupled ocean–atmosphere model (HadCM3) have so far been unable to fully predict the large scale of sea surface warming in the high latitudes. If climate proxies accurately represent Pliocene conditions, and if no weakness exists in the physics of the model, then model boundary conditions may be in error. Here we alter a single boundary condition (bathymetry) to examine if Pliocene high latitude warming was aided by an increase in poleward heat transport due to changes in the subsidence of North Atlantic Ocean ridges. We find an increase in both Arctic sea surface temperature and deepwater production in model experiments that incorporate a deepened Greenland–Scotland Ridge. These results offer both a mechanism for the warming in the North Atlantic and Arctic Oceans indicated by numerous proxies and an explanation for the apparent disparity between proxy data and model simulations of Pliocene northern North Atlantic and Arctic Ocean conditions. Determining the causes of Pliocene warmth remains critical to fully understanding comparisons of the Pliocene warm period to possible future climate change scenarios.

Palaeogeography, Palaeoclimatology, Palaeoecology↗
USGS2011

Genetic discontinuity among regional populations of Lophelia perfusa in the North Atlantic Ocean

Knowledge of the degree to which populations are connected through larval dispersal is imperative to effective management, yet little is known about larval dispersal ability or population connectivity in Lophelia pertusa, the dominant framework-forming coral on the continental slope in the North Atlantic Ocean. Using nine microsatellite DNA markers, we assessed the spatial scale and pattern of genetic connectivity across a large portion of the range of L. pertusa in the North Atlantic Ocean. A Bayesian modeling approach found four distinct genetic groupings corresponding to ocean regions: Gulf of Mexico, coastal southeastern U.S., New England Seamounts, and eastern North Atlantic Ocean. An isolation-by-distance pattern was supported across the study area. Estimates of pairwise population differentiation were greatest with the deepest populations, the New England Seamounts (average FST = 0.156). Differentiation was intermediate with the eastern North Atlantic populations (FST = 0.085), and smallest between southeastern U.S. and Gulf of Mexico populations (FST = 0.019), with evidence of admixture off the southeastern Florida peninsula. Connectivity across larger geographic distances within regions suggests that some larvae are broadly dispersed. Heterozygote deficiencies were detected within the majority of localities suggesting deviation from random mating. Gene flow between ocean regions appears restricted, thus, the most effective management scheme for L. pertusa involves regional reserve networks

Conservation Genetics↗
USGS2011

Building transparent data access for ocean observatories: Coordination of U.S. IOOS DMAC with NSF's OOI Cyberinfrastructure

The NOAA-led U.S. Integrated Ocean Observing System (IOOS) and the National Science Foundation's Ocean Observatories Initiative (OOI) have been collaborating since 2007 on advanced tools and technologies that ensure open access to ocean observations and models. Initial collaboration focused on serving ocean data via cloud computing-a key component of the OOI cyberinfrastructure (CI) architecture. As the OOI transitioned from planning to execution in the Fall of 2009, an OOI/IOOS team developed a customer-based "use case" to align more closely with the emerging objectives of OOI-CI team's first software release scheduled for Summer 2011 and provide a quantitative capacity for stress-testing these tools and protocols. A requirements process was initiated with coastal modelers, focusing on improved workflows to deliver ocean observation data. Accomplishments to date include the documentation and assessment of scientific workflows for two "early adopter" modeling teams from IOOS Regional partners (Rutgers-the State University of New Jersey and University of Hawaii's School of Ocean and Earth Science and Technology) to enable full understanding of data sources and needs; generation of all-inclusive lists of the data sets required and those obtainable through IOOS; a more complete understanding of areas where IOOS can expand data access capabilities to better serve the needs of the modeling community; and development of "data set agents" (software) to facilitate data acquisition from numerous data providers and conversions of the data format to the OOI-CI canonical form. ?? 2011 MTS.

Conference Paper↗
USGS2012

The Geoid: Effect of compensated topography and uncompensated oceanic trenches

The geoid is becoming increasingly important in interpretation of global tectonics. Most of the topography of the earth is isostatically compensated, so removal of its effect from the geoid is appropriate before tectonic modeling. The oceanic trenches, however, are dynamically depressed features and cannot be isostatically compensated in the classical way. Continental topography compensated at 35 km gives intracontinental geoidal undulations of up to 15 m over mountain ranges in a spherical harmonic expansion to order and degree 22. Oceanic topography compensated at 40 km, reasonable for the thermally supported long wavelengths, matches the +10 m difference between old continents and old oceans in a detailed NASA/GSFC geoid. Removing the assumed compensation for the oceanic trenches leaves negative anomalies of up to 9 m amplitude caused by their uncompensated mass deficit. This mass deficit acts as a partial "regional compensation" for the excess mass of the subducting slabs, and partly explains why geoidal (and gravity) anomalies over the cold slabs are less than thermal models predict.

Geophysical Research Letters↗
USGS2005

Movement of the saltwater interface in the surficial aquifer system in response to hydrologic stresses and water-management practices, Broward County, Florida

A study was conducted to evaluate the relation between water-level fluctuations and saltwater intrusion in Broward County, Florida. The objective was achieved through data collection at selected wells in Broward County and through the development of a variable-density ground-water flow model. The numerical model is representative of many locations in Broward County that contain a well field, control structure, canal, the Intracoastal Waterway, and the Atlantic Ocean. The model was used to simulate short-term movement (from tidal fluctuations to monthly changes) and long-term movement (greater than 10 years) of the saltwater interface resulting from changes in rainfall, well-field withdrawals, sea-level rise, and upstream canal stage. The SEAWAT code, which is a combined version of the computer codes, MODFLOW and MT3D, was used to simulate the complex variable-density flow patterns. Model results indicated that the canal, control structure, and sea level have major effects on ground-water flow. For periods greater than 10 years, the upstream canal stage controls the movement and location of the saltwater interface. If upstream canal stage is decreased by 1 foot (0.3048 meter), the saltwater interface takes 50 years to move inland and stabilize. If the upstream canal stage is then increased by 1 foot (0.3048 meter), the saltwater interface takes 90 years to move seaward and stabilize. If sea level rises about 48 centimeters over the next 100 year as predicted, then inland movement of the saltwater interface may cause well-field contamination. For periods less than 10 years, simulation results indicated that a 3-year drought with increased well-field withdrawals probably will not have long-term effects on the position of the saltwater interface in the Biscayne aquifer. The saltwater interface returns to its original position in less than 10 years. Model results, however, indicated that the interface location in the lower part of the surficial aquifer system takes longer than 10 years to recover from a drought. Additionally, rainfall seems to have the greatest effect on saltwater interface movement in areas some distance from canals, but the upstream canal stage has the greatest effect on the movement of the saltwater interface near canals. Field data indicated that saltwater interface movement includes short-term fluctuations caused by tidal fluctuations and long-term seasonal fluctuations. Statistical analyses of daily-averaged data indicated that the saltwater interface moves in response to pumpage, rainfall, and upstream canal stage. In areas near the canal, the saltwater interface is most affected by canal stage because water-management structures control the stage in the upstream part of the canal and allow movement of the saltwater interface. In areas away from the canal, the saltwater interface is most affected by pumpage and rainfall, depending on the location of well fields. Data analyses also revealed that rainfall changes the vertical flow direction in the Biscayne aquifer. Results from the study indicated that upstream canal stage substantially affects the long-term position of the saltwater interface in the surficial aquifer system. The saltwater interface moves faster inland than seaward because of changes in upstream canal stage. For short-term problems, such as drought, the threat of saltwater intrusion in the Biscayne aquifer does not appear to be severe if the well-field withdrawal is increased; however, this conclusion is based on the assumption that well-field withdrawals will decrease once the drought is over. Sea-level rise may be a potential threat to the water supply in Broward County as the saltwater interface moves inland toward well fields.

Florida↗
USGS2022

Impacts of the ocean-atmosphere coupling into the very short range prediction system during the impact of Hurricane Matthew on Cuba

The main goal of this investigation is analyzing the impact of insert the ocean-atmosphere coupling into the very short range prediction system of Cuba. The ocean-atmosphere coupled components of the Coupled Ocean-Atmosphere-Wave-Sediment Transport Modeling System are used for this purpose and the hurricane Matthew is selected as study case. Two experiments are performed: first, using a dynamic sea surface temperature, updated daily in the atmospheric model WRF; and second using adynamic coupling between the atmospheric and an oceanic models. For the simulated track, the best results are obtained with the coupled system. The impact of coupling on the maximum wind velocities and minimum central pressure is studied. In the coupled system the sea surface temperature has more influence in the surface latent heat fluxes. Also, with this methodology the dry footprint and the behavior of the precipitation field in the presence of a hurricane are studied. This analysis shows that the hurricane acts like an open and self-sustained system in the numerical experiments. The highest differences in the precipitation simulations are in the significant convective area inside the hurricane.

Ciência e Natura↗
USGS2000

Streamflow changes in the Sierra Nevada, California, simulated using a statistically downscaled general circulation model scenario of climate change

Simulations of future climate using general circulation models (GCMs) suggest that rising concentrations of greenhouse gases may have significant consequences for the global climate. Of less certainty is the extent to which regional scale (i.e., sub-GCM grid) environmental processes will be affected. In this chapter, a range of downscaling techniques are critiqued. Then a relatively simple (yet robust) statistical downscaling technique and its use in the modelling of future runoff scenarios for three river basins in the Sierra Nevada, California, is described. This region was selected because GCM experiments driven by combined greenhouse-gas and sulphate-aerosol forcings consistently show major changes in the hydro-climate of the southwest United States by the end of the 21st century. The regression-based downscaling method was used to simulate daily rainfall and temperature series for streamflow modelling in three Californian river basins under current-and future-climate conditions. The downscaling involved just three predictor variables (specific humidity, zonal velocity component of airflow, and 500 hPa geopotential heights) supplied by the U.K. Meteorological Office couple ocean-atmosphere model (HadCM2) for the grid point nearest the target basins. When evaluated using independent data, the model showed reasonable skill at reproducing observed area-average precipitation, temperature, and concomitant streamflow variations. Overall, the downscaled data resulted in slight underestimates of mean annual streamflow due to underestimates of precipitation in spring and positive temperature biases in winter. Differences in the skill of simulated streamflows amongst the three basins were attributed to the smoothing effects of snowpack on streamflow responses to climate forcing. The Merced and American River basins drain the western, windward slope of the Sierra Nevada and are snowmelt dominated, whereas the Carson River drains the eastern, leeward slope and is a mix of rainfall runoff and snowmelt runoff. Simulated streamflow in the American River responds rapidly and sensitively to daily-scale temperature and precipitation fluctuations and errors; in the Merced and Carson Rivers, the response to the same short-term influences is much less. Consequently, the skill of simulated flows was significantly lower in the American River model than in the Carson and Merced. The physiography of the three basins also accounts for differences in their sensitivities to future climate change. Increases in winter precipitation exceeding +100% coupled with mean temperature rises greater than +2°C result in increased winter streamflows in all three basins. In the Merced and Carson basins, these streamflow increases reflect large changes in winter snowpack, whereas the streamflow changes in the lower elevation American basin are driven primarily by rainfall runoff. Furthermore, reductions in winter snowpack in the American River basin, owing to less precipitation falling as snow and earlier melting of snow at middle elevations, lead to less spring and summer streamflow. Taken collectively, the downscaling results suggest significant changes to both the timing and magnitude of streamflows in the Sierra Nevada by the end of the 21st Century. In the higher elevation basins, the HadCM2 scenario implies more annual streamflow and more streamflow during the spring and summer months that are critical for water-resources management in California. Depending on the relative significance of rainfall runoff and snowmelt, each basin responds in its own way to regional climate forcing. Generally, then, climate scenarios need to be specified — by whatever means — with sufficient temporal and spatial resolution to capture subtle orographic influences if projections of climate-change responses are to be useful and reproducible.

California↗
USGS1998

CRUST 5.1: A global crustal model at 5° x 5°

We present a new global model for the Earth's crust based on seismic refraction data published in the period 1948-1995 and a detailed compilation of ice and sediment thickness. An extensive compilation of seismic refraction measurements has been used to determine the crustal structure on continents and their margins. Oceanic crust is modeled with both a standard model for normal oceanic crust, and variants for nonstandard regions, such as oceanic plateaus. Our model (CRUST 5.1) consists of 2592 5° × 5° tiles in which the crust and uppermost mantle are described by eight layers: (1) ice, (2) water, (3) soft sediments, (4) hard sediments, (5) crystalline upper, (6) middle, (7) lower crust, and (8) uppermost mantle. Topography and bathymetry are adopted from a Mandard database (ETOPO-5). Compressional wave velocity in each layer is based on field measurements, and shear wave velocity and density are estimated using recently published empirical V p -V s and V p -density relationships. The crustal model differs from previous models in that (1) the thickness and seismic/density structure of sedimentary basins is accounted for more completely, (2) the velocity structure of unmeasured regions is estimated using statistical averages that are based on a significantly larger database of crustal structure, (3) the compressional wave, shear wave, and density structure have been explicitly specified using newly available constraints from field and laboratory studies. Thus - this global crustal model is based on substantially more data than previous models and differs from them in many important respects. A new map of the thickness of the Earth's crust is presented, and we illustrate the application of this model by using it to provide the crustal correction for surface wave phase velocity maps. Love waves at 40 s are dominantly sensitive to crustal structure, and there is a very close correspondence between observed phase velocities at this period and those predicted by CRUST 5.1. We find that the application of crustal corrections to long-period (167 s) Rayleigh waves significantly increases the variance in the phase velocity maps and strengthens the upper mantle , velocity anomalies beneath stable continental regions. A simple calculation of crustal isostacy indicates significant lateral variations in upper mantle density. The model CRUST 5.1 provides a complete description of the physical properties of the Earth's crust at a scale of 5° × 5° and can be used for a wide range of seismological and nonseismological problems.

Journal of Geophysical Research B: Solid Earth↗
USGS2003

Crustal magnetization and accretion at the Southwest Indian Ridge near the Atlantis II fracture zone, 0-25 Ma

We analyze geophysical data that extend from 0 to 25-Myr-old seafloor on both flanks of the Southwest Indian Ridge (SWIR). Lineated marine magnetic anomalies are consistent and identifiable within the study area, even over seafloor lacking a basaltic upper crust. The full spreading rate of 14 km/Myr has remained nearly constant since at least 20 Ma, but crustal accretion has been highly asymmetric, with half rates of 8.5 and 5.5 km/Myr on the Antarctic and African flanks, respectively. This asymmetry may be unique to a ???400 km wide corridor between large-offset fracture zones of the SWIR. In contrast to the Mid-Atlantic Ridge, crustal magnetization amplitudes correlate directly with seafloor topography along the present-day rift valleys. This pattern appears to be primarily a function of along-axis variations in crustal thickness, rather than magnetic mineralogy. Off-axis, magnetization amplitudes at paleo-segment ends are more positive than at paleo-segment midpoints, suggesting the presence of an induced component of magnetization within the lower crust or serpentinized upper mantle. Alteration of the magnetic source layer at paleo-segment midpoints reduces magnetization amplitudes by 70-80% within 20 Myr of accretion. Magnetic and Ocean Drilling Program (ODP) Hole 735B data suggest that the lower crust cooled quickly enough to lock in a primary thermoremanent magnetization that is in phase with that of the overlying upper crust. Thus magnetic polarity boundaries within the intrusive lower crust may be steeper than envisioned in prior models of ocean crustal magnetization. As the crust ages, the lower crust becomes increasingly important in preserving marine magnetic stripes.

Journal of Geophysical Research B: Solid Earth↗
USGS2002

Nearly frictionless faulting by unclamping in long-term interaction models

In defiance of direct rock-friction observations, some transform faults appear to slide with little resistance. In this paper finite element models are used to show how strain energy is minimized by interacting faults that can cause long-term reduction in fault-normal stresses (unclamping). A model fault contained within a sheared elastic medium concentrates stress at its end points with increasing slip. If accommodating structures free up the ends, then the fault responds by rotating, lengthening, and unclamping. This concept is illustrated by a comparison between simple strike-slip faulting and a mid-ocean-ridge model with the same total transform length; calculations show that the more complex system unclapms the transforms and operates at lower energy. In another example, the overlapping San Andreas fault system in the San Francisco Bay region is modeled; this system is complicated by junctions and stepovers. A finite element model indicates that the normal stress along parts of the faults could be reduced to hydrostatic levels after ???60-100 k.y. of system-wide slip. If this process occurs in the earth, then parts of major transform fault zones could appear nearly frictionless.

Geology↗
USGS2005

Performance of four turbulence closure models implemented using a generic length scale method

A two-equation turbulence model (one equation for turbulence kinetic energy and a second for a generic turbulence length-scale quantity) proposed by Umlauf and Burchard [J. Marine Research 61 (2003) 235] is implemented in a three-dimensional oceanographic model (Regional Oceanographic Modeling System; ROMS v2.0). These two equations, along with several stability functions, can represent many popular turbulence closures, including the k-kl (Mellor-Yamada Level 2.5), k-??, and k-?? schemes. The implementation adds flexibility to the model by providing an unprecedented range of turbulence closure selections in a single 3D oceanographic model and allows comparison and evaluation of turbulence models in an otherwise identical numerical environment. This also allows evaluation of the effect of turbulence models on other processes such as suspended-sediment distribution or ecological processes. Performance of the turbulence models and sediment-transport schemes is investigated with three test cases for (1) steady barotropic flow in a rectangular channel, (2) wind-induced surface mixed-layer deepening in a stratified fluid, and (3) oscillatory stratified pressure-gradient driven flow (estuarine circulation) in a rectangular channel. Results from k-??, k-??, and gen (a new closure proposed by Umlauf and Burchard [J. Marine Research 61 (2003) 235]) are very similar for these cases, but the k-kl closure results depend on a wall-proximity function that must be chosen to suit the flow. Greater variations appear in simulations of suspended-sediment concentrations than in salinity simulations because the transport of suspended-sediment amplifies minor variations in the methods. The amplification is caused by the added physics of a vertical settling rate, bottom stress dependent resuspension, and diffusive transport of sediment in regions of well mixed salt and temperature. Despite the amplified sensitivity of sediment to turbulence models in the estuary test case, the four closures investigated here all generated estuarine turbidity maxima that were similar in their shape, location, and concentrations.

Ocean Modelling↗
USGS2003

Mid-Holocene climates of the Americas: a dynamical response to changed seasonality

Simulations of the climatic response to mid-Holocene (6 ka BP) orbital forcing with two coupled ocean–atmosphere models (FOAM and CSM) show enhancement of monsoonal precipitation in parts of the American Southwest, Central America and northernmost South America during Northern Hemisphere summer. The enhanced onshore flow that brings precipitation into Central America is caused by a northward displacement of the inter-tropical convergence zone, driven by cooling of the equatorial and warming of the northern subtropical and mid-latitude ocean. Ocean feedbacks also enhance precipitation over the American Southwest, although the increase in monsoon precipitation there is largely driven by increases in land-surface temperature. The northward shift in the equatorial precipitation band that causes enhanced precipitation in Central America and the American Southwest has a negative feedback effect on monsoonal precipitation in northern South America. The simulations demonstrate that mid-Holocene aridity in the mid-continent of North America is dynamically linked to the orbitally induced enhancement of the summer monsoon in the American Southwest, with a spatial structure (wet in the Southwest and dry in the mid-continent) similar to that found in strong monsoon years today. Changes in winter precipitation along the west coast of North America, in Central America and along the Gulf Coast, caused by southward-displacement of the westerly storm tracks, indicate that changes in the Northern Hemisphere winter monsoon also play a role in regional climate changes during the mid-Holocene. Although the simulations with FOAM and CSM differ in detail, the general mechanisms and patterns are common to both. The model results thus provide a coherent dynamical explanation for regional patterns of increased or decreased aridity shown by vegetation, lake status and aeolian data from the Americas.

Climate Dynamics↗
USGS2020

Modeling soil porewater salinity response to drought in tidal freshwater forested wetlands

There is a growing concern about the adverse effects of saltwater intrusion via tidal rivers, streams and creeks into tidal freshwater forested wetlands (TFFW) due to sea‐level rise (SLR) and intense and extended drought events. However, the magnitude and duration of porewater salinity in exceedance of plant salinity stress threshold (2 practical salinity units, psu) and the controlling factors remain unclear. In this study, we developed a TFFW soil porewater salinity model, in which the feedback mechanisms between soil salinity and evapotranspiration and hydraulic conductivity were incorporated. We selected sites (upper, middle, lower tidal freshwater forest sites and oligohaline marsh site) along the coastal floodplains of two rivers, the Waccamaw River (SC, USA) and the Savannah River (GA and SC, USA), that represent landscape salinity gradients from tidal influence of the Atlantic Ocean. The model results agreed well with field measurements and revealed that with drought‐induced saltwater intrusion, the mean annual soil porewater salinity and duration of elevated soil porewater salinity (> 2 psu) increased significantly compared to the normal (non‐drought) condition, posing a threat to the health and ecosystem services of TFFW even in the absence of SLR. Model results also showed more severe salinity stress under drought for the lower forest sites along the two rivers, where soil salinity values have already been at or in exceedance of the 2 psu threshold.

Georgia, South Carolina↗
USGS1999

Stasis and extinction of Silurian (Llandovery-Wenlock) trilobite associations related to oceanic cyclicity

Silurian trilobites of the central United States belong to a series of temporally-successive associations which appeared abruptly, maintained taxonomic stasis for a time, and then disappeared abruptly. Their disappearance resulted from global perturbations of short-term duration and moderate magnitude, which caused substantial taxonomic replacement but no reorganization of major ecosystems. The most significant extinction and replacement in Silurian trilobite associations in the study area occurs near the Llandovery-Wenlock boundary. This turnover in trilobite associations appears to correspond to Jeppsson's Ireviken Event in his model of oceanic and climatic cyclicity. Major sea-level changes earlier in the Llandovery did not have a similar impact on trilobite associations.

Journal of Paleontology↗
USGS2021

An efficient method to calculate depth-integrated, phase-averaged momentum balances in non-hydrostatic models

Analysis of the mean (wave-averaged) momentum balance is a common approach used to explain the physical forcing driving wave set-up and mean currents in the nearshore zone. Traditionally this approach has been applied to phase-averaged models but has more recently been applied to phase-resolving models using post-processing, whereby model output is used to calculate each of the momentum terms. While phase-resolving models have the advantage of capturing the nonlinear properties of waves propagating in the nearshore (making them advantageous to enhance understanding of nearshore processes), the post-processing calculation of the momentum terms does not guarantee that the momentum balance closes. We show that this is largely due to the difficulty (or impossibility) of being consistent with the numerical approach. If the residual is of a similar magnitude as any of the relevant momentum terms (which is common with post-processing methods as we show), the analysis is largely compromised. Here we present a new method to internally calculate and extract the depth-integrated, mean momentum terms in the phase-resolving non-hydrostatic wave-flow model SWASH in a manner that is consistent with the numerical implementation. Further, we demonstrate the utility of the new method with two existing physical model studies. By being consistent with the numerical framework, the internal method calculates the momentum terms with a much lower residual at computer precision, combined with greatly reduced calculation time and output storage requirements compared to post-processing techniques. The method developed here allows the accurate evaluation of the depth-integrated, mean momentum terms of wave-driven flows while taking advantage of the more complete representation of the wave dynamics offered by phase-resolving models. Furthermore, it provides an opportunity for advances in the understanding of nearshore processes particularly at more complex sites where wave nonlinearity and energy transfers are important.

Ocean Modelling↗
USGS2016

Dynamic reusable workflows for ocean science

Digital catalogs of ocean data have been available for decades, but advances in standardized services and software for catalog search and data access make it now possible to create catalog-driven workflows that automate — end-to-end — data search, analysis and visualization of data from multiple distributed sources. Further, these workflows may be shared, reused and adapted with ease. Here we describe a workflow developed within the US Integrated Ocean Observing System (IOOS) which automates the skill-assessment of water temperature forecasts from multiple ocean forecast models, allowing improved forecast products to be delivered for an open water swim event. A series of Jupyter Notebooks are used to capture and document the end-to-end workflow using a collection of Python tools that facilitate working with standardized catalog and data services. The workflow first searches a catalog of metadata using the Open Geospatial Consortium (OGC) Catalog Service for the Web (CSW), then accesses data service endpoints found in the metadata records using the OGC Sensor Observation Service (SOS) for in situ sensor data and OPeNDAP services for remotely-sensed and model data. Skill metrics are computed and time series comparisons of forecast model and observed data are displayed interactively, leveraging the capabilities of modern web browsers. The resulting workflow not only solves a challenging specific problem, but highlights the benefits of dynamic, reusable workflows in general. These workflows adapt as new data enters the data system, facilitate reproducible science, provide templates from which new scientific workflows can be developed, and encourage data providers to use standardized services. As applied to the ocean swim event, the workflow exposed problems with two of the ocean forecast products which led to improved regional forecasts once errors were corrected. While the example is specific, the approach is general, and we hope to see increased use of dynamic notebooks across the geoscience domains.

Journal of Marine Science and Engineering↗
USGS1983

Crustal structure of southwestern Saudi Arabia

The southwestern Arabian Shield is composed of uplifted Proterozoic metamorphic and plutonic rocks. The Shield is bordered on the southwest by Cenozoic sedimentary and igneous rocks of the Red Sea paar and on the east by the Arabian Platform, an area of basin sedimentation throughout Phanerozoic time. The Shield appears to have been formed by successive episodes of island arc volcanism and sea-floor spreading, followed by several cycles of compressive tectonism and metamorphism. An interpretation and synthesis of a deep-refraction seismic profile from the Riyadh area to the Farasan Islands, and regional gravity, aeromagnetic, heat flow, and surface geologic data have yielded a self-consistent regional-scale model of the crust and upper mantle for this area. The model consists of two 20 km-thick layers of crust with an average compressional wave velocity in the upper crust of about 6.3 km/s and an average velocity in the lower. crust of about 7.0 km/s. This crust thins abruptly to less than 20 km near the southwestern end of the profile where Precambrian outcrops abut the Cenozoic rocks and to 8 km beneath the Farasan Islands. The data over the coastal plain and Red Sea shelf areas are fit satisfactorily by an oceanic crustal model. A major lateral velocity inhomogeneity in the crust is inferred about 25 km northeast of Sabhah and is supported by surface geologic evidence. The major velocity discontinuities occur at about the same depth across the entire Shield and are interpreted to indicate horizontal metamorphic stratification of the Precambrian crust. Several lateral inhomogenities in both the upper and lower .crust of the . Shield are interpreted, to indicate bulk compositional variations. The subcrustal portion of the model is composed of a hot, low-density lithosphere beneath the Red Sea which is systematically cooler and denser to the northeast. This model provides a mechanism which explains the observed topographic uplift, regional gravity pattern, heat flow, and mantle compressional wave velocities. Such a lithosphere could be produced by upwelling of hot asthenosphere beneath the Red Sea which then flows laterally beneath the lithosphere of the Arabian Plate.

Open-File Report↗
USGS2004

DEM generation and tidal deformation detection for Sulzberger Ice Shelf, West Antarctica using SAR interferometry

In this study we generated a relative Digital Elevation Model (DEM) over the Sulzberger Ice Shelf, West Antarctica using ERS1/2 synthetic aperture radar (SAR) interferometry data. Four repeat pass differential interferograms are used to find the grounding zone and to classify the study area. An interferometrically derived DEM is compared with laser altimetry profile from ICESat. Standard deviation of the relative height difference is 5.12 m and 1.34 m in total length of the profile and at the center of the profile respectively. The magnitude and the direction of tidal changes estimated from interferogram are compared with those predicted tidal differences from four ocean tide models. Tidal deformation measured in InSAR is -16.7 cm and it agrees well within 3 cm with predicted ones from tide models.

Conference Paper↗