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.

SourceResearch connectionFrom yearThrough year

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

At least 181 records · Page 10Linked to original sources

USGS2022

Modeling protected species distributions and habitats to inform siting and management of pioneering ocean industries: A case study for Gulf of Mexico aquaculture

Marine Spatial Planning (MSP) provides a process that uses spatial data and models to evaluate environmental, social, economic, cultural, and management trade-offs when siting (i.e., strategically locating) ocean industries. Aquaculture is the fastest-growing food sector in the world. The United States (U.S.) has substantial opportunity for offshore aquaculture development given the size of its exclusive economic zone, habitat diversity, and variety of candidate species for cultivation. However, promising aquaculture areas overlap many protected species habitats. Aquaculture siting surveys, construction, operations, and decommissioning can alter protected species habitat and behavior. Additionally, aquaculture-associated vessel activity, underwater noise, and physical interactions between protected species and farms can increase the risk of injury and mortality. In 2020, the U.S. Gulf of Mexico was identified as one of the first regions to be evaluated for offshore aquaculture opportunities as directed by a Presidential Executive Order. We developed a transparent and repeatable method to identify aquaculture opportunity areas (AOAs) with the least conflict with protected species. First, we developed a generalized scoring approach for protected species that captures their vulnerability to adverse effects from anthropogenic activities using conservation status and demographic information. Next, we applied this approach to data layers for eight species listed under the Endangered Species Act, including five species of sea turtles, Rice’s whale, smalltooth sawfish, and giant manta ray. Next, we evaluated four methods for mathematically combining scores (i.e., Arithmetic mean, Geometric mean, Product, Lowest Scoring layer) to generate a combined protected species data layer. The Product approach provided the most logical ordering of, and the greatest contrast in, site suitability scores. Finally, we integrated the combined protected species data layer into a multi-criteria decision-making modeling framework for MSP. This process identified AOAs with reduced potential for protected species conflict. These modeling methods are transferable to other regions, to other sensitive or protected species, and for spatial planning for other ocean-uses.

PLoS ONE↗

USGS2006

Sensitivity of Last Glacial Maximum climate to uncertainties in tropical and subtropical ocean temperatures

The faunal and floral gradients that underlie the CLIMAP (1981) sea-surface temperature (SST) reconstructions for the Last Glacial Maximum (LGM) reflect ocean temperature gradients and frontal positions. The transfer functions used to reconstruct SSTs from biologic gradients are biased, however, because at the warmest sites they display inherently low sensitivity in translating fauna to SST and they underestimate SST within the euphotic zones where the pycnocline is strong. Here we assemble available data and apply a statistical approach to adjust for hypothetical biases in the faunal-based SST estimates of LGM temperature. The largest bias adjustments are distributed in the tropics (to address low sensitivity) and subtropics (to address underestimation in the euphotic zones). The resulting SSTs are generally in better agreement than CLIMAP with recent geochemical estimates of glacial-interglacial temperature changes. We conducted a series of model experiments using the GENESIS general atmospheric circulation model to assess the sensitivity of the climate system to our bias-adjusted SSTs. Globally, the new SST field results in a modeled LGM surface-air cooling relative to present of 6.4 ??C (1.9 ??C cooler than that of CLIMAP). Relative to the simulation with CLIMAP SSTs, modeled precipitation over the oceans is reduced by 0.4 mm d-1 (an anomaly -0.4 versus 0.0 mm d-1 for CLIMAP) and increased over land (an anomaly -0.2 versus -0.5 mm d-1 for CLIMAP). Regionally strong responses are induced by changes in SST gradients. Data-model comparisons indicate improvement in agreement relative to CLIMAP, but differences among terrestrial data inferences and simulated moisture and temperature remain. Our SSTs result in positive mass balance over the northern hemisphere ice sheets (primarily through reduced summer ablation), supporting the hypothesis that tropical and subtropical ocean temperatures may have played a role in triggering glacial changes at higher latitudes.

Quaternary Science Reviews↗

USGS2026

Multi-proxy thermal history of basin heating during Cordilleran orogenesis in the Magallanes-Austral retroarc foreland basin, Patagonian Andes

Resolving thermal histories in sedimentary basins is crucial for interpreting orogenic growth, basin burial, and tectonic processes during Cordilleran orogenesis. In the Magallanes–Austral Basin, Patagonian Andes, we integrate new (U-Th)/He thermochronology, vitrinite reflectance (%R o ), calcite-cement clumped isotope data and thermal history modelling to resolve the origin of the regionally extensive Paleogene unconformity (51°S–50°S). Thermal history modelling results require post-depositional heating of Palaeocene (Danian–Selandian) strata below the unconformity and suggest maximum burial temperatures of 87°C–101°C (55–52 Ma) and 89°C–92°C (18–16 Ma). For lower Eocene strata above the unconformity, Miocene burial temperatures (89°C–92°C) are consistent with calcite cement formation temperatures (~62°C–92°C) from carbonate clumped isotopes. Our results indicate that basin burial and heating between ca. 60 and 52 Ma were likely driven by shallowing of the subducting Farallon plate and enhanced plate coupling preceding arrival of the Farallon–Phoenix mid-ocean ridge. Subsequent basin inversion and cooling from ca. 52 to 44 Ma correspond with subduction of this mid-ocean ridge. Refined thermal models, constrained by expanded thermochronometric and organic maturation datasets, indicate that up to ~1.7–2.0 km of proximal foreland basin strata were removed during uplift and erosion across the Paleogene basin margin. A return to basin subsidence beginning ca. 44 Ma may reflect dynamic subsidence after passage of the mid-ocean ridge and renewed coupling between the fold-thrust belt and foreland basin system. Neogene thermal histories document continued subsidence, localized hot orogenic fluid flow along stratigraphic boundaries, followed by a final phase of basin inversion and cooling at ca. 18–16 Ma, which we attribute to regional uplift associated with Chile ridge subduction. Altogether, this study demonstrates that multiple thermal indices when analysed and modelled can provide clarity for tectonic and stratigraphic events that affect foreland basins.

Patagonian Andes↗

USGS1972

Structural profile of the northwestern Caribbean

A seismic reflection and gravity profile across the continental margin of the Yucatan Peninsula, Yucatan Basin, Cayman Ridge, and Cayman Trough suggests that sediments in the Yucatan Basin consist of a thick succession of beds dominated by turbidites that overlie a thick but irregular sequence of beds, probably dominated by pelagic deposits. The so-called “Carib beds”, present elsewhere in the Caribbean, are not evident in the part of the basin crossed by this profile. The sedimentary section rests on a acoustic basement that probably represents the top of oceanic layer 2. A gravity model indicates that the crust beneath the Yucatan Basin is thin and therefore probably is oceanic in character. The crust thickens southward under the Cayman Ridge but thins again beneath the Cayman Trough. This local thickening is consistent with the suggestion that the Cayman Ridge is a rifted part of the Nicaraguan Rise.

Earth and Planetary Science Letters↗

USGS2022

Phenomenological and mechanistic modeling of recreational angling behavior using creel data

Improved understanding of angling behavior in recreational fisheries can help managers account for partial controllability in systems in which angling effort is not directly regulated. Relevant aspects of angling behavior include fishing participation, site choice, and catch-and-release actions and can be considered at scales of both aggregate fishing effort and individual-decision-making. Using creel survey data collected across eight fishery seasons, we predicted aggregate angling effort (e.g., daily number of boating trips engaged in fishing) and characterized influences of individual trip-based site choice for recreational fisheries acting on fall-run Chinook salmon and coho salmon near the mouth of the Columbia River, United States of America. We applied an overdispersed Poisson likelihood-based generalized linear model with an autoregressive structure to aggregate boating effort and a multinomial logit model to trip-based site choice decisions among separate ocean and estuary fishing zones. Predictive models explained 71% and 79% of out-of-sample and in-sample variability in aggregate effort, respectively, and included the covariates weekend status, Chinook salmon catch rate, tidal range, and pre-season expectations of fish abundance. In addition to reinforcing the importance of Chinook salmon catch rate and tidal range, site choice model selection revealed influences of weather, fishery restrictions, expected fishery season lengths, and the individual-specific characteristics of boat length and guide status on decision-making. Model results provide predictive models for short-term fishery planning, highlight heterogeneity in individual decision-making, and illustrate the value of standardized creel data for evaluating angling behavior.

Oregon↗

USGS2019

Hydrodynamic and morphologic response of a back-barrier estuary to an extratropical storm

We investigated the hydrodynamic and morphologic response of Barnegat Bay-Little Egg Harbor, New Jersey, USA to Hurricane Sandy. We implemented a three-dimensional, coupled ocean-wave-sediment transport model of the estuary and explored the role of offshore water levels, offshore waves, local winds and waves by systematically removing forcings from a series of simulations. Offshore water levels had the largest impact on water levels in the bay while waves and local wind forcing created substantial spatial variation along the longitudinal axis of the bay. The shape of the bay and its orientation relative to the storm track influenced the response to winds and restricted the maximum water levels in the northern bay and reduced the maximum volume of surge. Basin-average hydrodynamic residence time was reduced by 40%, though its typical spatial distribution remained during the storm. Wave and current-induced bed shear stress resuspended fine sediment resulting in net erosion from the shoals with ensuing net deposition over fringing low-lying land. The net sediment exchange between the bay and the ocean was several times smaller than the exchange at the peak of the storm resulting in negligible net change in the bay volume. Overall, our results suggest that water level responses are highly sensitive to the specific orientation of storm winds relative to the estuary, thereby limiting the utility of simple inundation models. The sediment transport patterns indicate that storms are an important mechanism for redistributing sediment from shoals to fringing wetlands, while net change to sediment budget can be negligible.

New Jersey↗

USGS2003

Fault slip and seismic moment of the 1700 Cascadia earthquake inferred from Japanese tsunami descriptions

The 1700 Cascadia earthquake attained moment magnitude 9 according to new estimates based on effects of its tsunami in Japan, computed coseismic seafloor deformation for hypothetical ruptures in Cascadia, and tsunami modeling in the Pacific Ocean. Reports of damage and flooding show that the 1700 Casscadia tsunami reached 1-5 m heights at seven shoreline sites in Japan. Three sets of estimated heights express uncertainty about location and depth of reported flooding, landward decline in tsunami heights from shorelines, and post-1700 land-level changes. We compare each set with tsunami heights computed from six Cascadia sources. Each source is vertical seafloor displacement calculated with a three-dimensional elastic dislocation model, for three sources the rupture extends the 1100 km length of the subduction zone and differs in width and shallow dip; for the other sources, ruptures of ordinary width extend 360-670 km. To compute tsunami waveforms, we use a linear long-wave approximation with a finite difference method, and we employ modern bathymetry with nearshore grid spacing as small as 0.4 km. The various combinations of Japanese tsunami heights and Cascadia sources give seismic moment of 1-9 ?? 1022 N m, equivalent to moment magnitude 8.7-9.2. This range excludes several unquantified uncertainties. The most likely earthquake, of moment magnitude 9.0, has 19 m of coseismic slip on an offshore, full-slip zone 1100 km long with linearly decreasing slip on a downdip partial-slip zone. The shorter rupture models require up to 40 m offshore slip and predict land-level changes inconsistent with coastal paleoseismological evidence. Copyright 2003 by the American Geophysical Union.

Journal of Geophysical Research B: Solid Earth↗

USGS2021

Multiple climate change-driven tipping points for coastal systems

As the climate evolves over the next century, the interaction of accelerating sea level rise (SLR) and storms, combined with confining development and infrastructure, will place greater stresses on physical, ecological, and human systems along the ocean-land margin. Many of these valued coastal systems could reach “tipping points,” at which hazard exposure substantially increases and threatens the present-day form, function, and viability of communities, infrastructure, and ecosystems. Determining the timing and nature of these tipping points is essential for effective climate adaptation planning. Here we present a multidisciplinary case study from Santa Barbara, California (USA), to identify potential climate change-related tipping points for various coastal systems. This study integrates numerical and statistical models of the climate, ocean water levels, beach and cliff evolution, and two soft sediment ecosystems, sandy beaches and tidal wetlands. We find that tipping points for beaches and wetlands could be reached with just 0.25 m or less of SLR (~ 2050), with > 50% subsequent habitat loss that would degrade overall biodiversity and ecosystem function. In contrast, the largest projected changes in socioeconomic exposure to flooding for five communities in this region are not anticipated until SLR exceeds 0.75 m for daily flooding and 1.5 m for storm-driven flooding (~ 2100 or later). These changes are less acute relative to community totals and do not qualify as tipping points given the adaptive capacity of communities. Nonetheless, the natural and human built systems are interconnected such that the loss of natural system function could negatively impact the quality of life of residents and disrupt the local economy, resulting in indirect socioeconomic impacts long before built infrastructure is directly impacted by flooding.

California↗

USGS2010

Stochastic model for the long-term transport of stored sediment in a river channel

We develop a stochastic model for the transport of stored sediment down a river channel. The model is based on probabilities of transition of particles among four different sediment storage reservoirs, called active (often mobilized), semiactive, inactive, and stable (hardly ever mobilized). The probabilities are derived from computed sediment residence times. Two aspects of sediment storage are investigated: flushing times of sediment out of a storage reservoir and changes in the quantity of sediment stored in different reservoirs due to seasonal sediment transport into, and out of, a reach. We apply the model to Redwood Creek, a gravel bed river in northern California. Although the Redwood Creek data set is incomplete, the application serves as an example of the sorts of analyses that can be done with the method. The application also provides insights into the sediment storage process. Sediment flushing times are highly dependent on the degree of interaction of the stable reservoir with the more mobile sediment reservoirs. The most infrequent and highest intensity storm events, which mobilize the stable reservoir, are responsible for the long-term shifts in sediment storage. Turnover times of channel sediment in all but the stable reservoir are on the order of 750 years, suggesting this is all the time needed for thorough interchange between these sediment compartments and cycling of most sediment particles from the initial reservoir to the ocean. Finally, the Markov model has adequately characterized sediment storage changes in Redwood Creek for 1947–1982, especially for the active reservoir. The model replicates field observation of the passage of a slug of sediment through the active reservoir of the middle reach of Redwood Creek in the 18 years following a major storm in 1964 that introduced large quantities of landslide debris to the channel.

California↗

USGS2012

Common coastal foraging areas for loggerheads in the Gulf of Mexico: Opportunities for marine conservation

Designing conservation strategies that protect wide-ranging marine species is a significant challenge, but integrating regional telemetry datasets and synthesizing modeled movements and behavior offer promise for uncovering distinct at-sea areas that are important habitats for imperiled marine species. Movement paths of 10 satellite-tracked female loggerheads ( Caretta caretta ) from three separate subpopulations in the Gulf of Mexico, USA, revealed migration to discrete foraging sites in two common areas at-sea in 2008, 2009, and 2010. Foraging sites were 102–904 km away from nesting and tagging sites, and located off southwest Florida and the northern Yucatan Peninsula, Mexico. Within 3–35 days, turtles migrated to foraging sites where they all displayed high site fidelity over time. Core-use foraging areas were 13.0–335.2 km 2 in size, in water <50 m deep, within a mean distance to nearest coastline of 58.5 km, and in areas of relatively high net primary productivity. The existence of shared regional foraging sites highlights an opportunity for marine conservation strategies to protect important at-sea habitats for these imperiled marine turtles, in both USA and international waters. Until now, knowledge of important at-sea foraging areas for adult loggerheads in the Gulf of Mexico has been limited. To better understand the spatial distribution of marine turtles that have complex life-histories, we propose further integration of disparate tracking data-sets at the oceanic scale along with modeling of movements to identify critical at-sea foraging habitats where individuals may be resident during non-nesting periods.

Florida↗

USGS2019

Relevance of wind stress and wave-dependent ocean surface roughness on the generation of winter meteotsunamis in Northern Gulf of Mexico

Meteotsunamis associated with passing squall lines are often observed ahead of cold fronts during winter seasons in Northern Gulf of Mexico. These types of meteotsunamis occur simultaneously with wind speed variations (~5-20 m/s) and sea-level atmospheric pressure oscillations (~1-6 hPa) with periods between 2 hours to several minutes. In order to enhance understanding of meteotsunami generation and propagation mechanisms, a Coupled-Ocean-Atmosphere-Wave-Sediment Transport (COAWST) modeling system is applied to one of the most intense winter meteotsunamis measured in Northern Gulf of Mexico in the last decade (2009-2018). The model verification with sea level and atmospheric observations show that the fully-coupled model is able to reproduce the timing and intensity of the 10-m wind and sea level atmospheric pressure fluctuations. The mean bias between observed and measured wind speeds and atmospheric pressure are 1.73 m/s and 0.63 hPa respectively. The maximum meteotsunami elevation and its timing are successfully captured by modeled (with a 7% underestimation of the maximum elevation). The relative effect of atmospheric pressure and wind stress divergence on meteotsunami generation is assessed with different numerical simulations. Results indicate that both wind stress and atmospheric pressure oscillations contributed to the generation of the meteotsunami. Wind stress was the dominant force in shallow waters (<15 m in this application), while the effects of atmospheric pressure disturbances dominated over areas with Froude number close to one (~40 m in this application). During the passage of the squall line, the sea surface became rougher in a sea state characterized by young and steep local ocean waves. Compared to a purely wind-speed-dependent roughness scheme, the application of a wave-dependent roughness parameterization improved in 37% modeled meteotsunami maximum elevation.

Northern Gulf of Mexico↗

USGS1979

The gravity field of the U.S. Atlantic continental margin

Approximately 39,000 km of marine gravity data collected during 1975 and 1976 have been integrated with U.S. Navy and other available data over the U.S. Atlantic continental margin between Florida and Maine to obtain a 10 mgal contour free-air gravity anomaly map. A maximum typically ranging from 0 to +70 mgal occurs along the edge of the shelf and Blake Plateau, while a minimum typically ranging from −20 to −80 mgal occurs along the base of the continental slope, except for a −140 mgal minimum at the base of the Blake Escarpment. Although the maximum and minimum free-air gravity values are strongly influenced by continental slope topography and by the abrupt change in crustal thickness across the margin, the peaks and troughs in the anomalies terminate abruptly at discrete transverse zones along the margin. These zones appear to mark major NW—SE fractures in the subsided continental margin and adjacent deep ocean basin, which separate the margin into a series of segmented basins and platforms. Rapid differential subsidence of crustal blocks on either side of these fractures during the early stages after separation of North America and Africa (Jurassic and Early Cretaceous) is inferred to be the cause of most of the gravity transitions along the length of margin. The major transverse zones are southeast of Charleston, east of Cape Hatteras, near Norfolk Canyon, off Delaware Bay, just south of Hudson Canyon and south of Cape Cod. Local Airy isostatic anomaly profiles (two-dimensional, without sediment corrections) were computed along eight multichannel seismic profiles. The isostatic anomaly values over major basins beneath the shelf and rise are generally between −10 and −30 mgal while those over the platform areas are typically 0 to +20 mgal. While a few isostatic anomaly profiles show local 10–20 mgal increases seaward of the East Coast Magnetic Anomaly (ECMA: inferred to mark the ocean-continent boundary), the lack of a consistent correlation indicates that the relationship of isostatic gravity anomalies to the magnetic anomalies and the ocean—continent transition is variable. Two-dimensional gravity models have been computed for two profiles off Cape Cod, Massachusetts and Cape May, New Jersey, where excellent reflection, refraction and magnetic control appear to define 10 and 12 km deep sedimentary basins beneath the shelf, respectively and 10 km deep basins beneath the rise. The basins are separated by a 6–8 km deep basement ridge which underlies the ECMA and appears to mark the landward edge of oceanic crust. The gravity models suggest that the oceanic crust is between 11 and 18 km thick beneath the ECMA, but decreases to a thickness of less than 8 km within the first 20–90 km to the southeast. In both profiles, the derived crustal thickness variations support the interpretation that the ECMA occurs over the ocean-continent boundary. The crust underlying the sedimentary cover appears to be 12 to 15 km thick on the landward side of the ECMA and gradually thickens to normal continental values of greater than 25 km within the first 60 to 110 km to the northwest. Multichannel seismic profiles across platform areas, such as Cape Hatteras and Cape Cod, indicate the ocean-continent transition zones there are much narrower than profiles across major sedimentary basins, such as the one off New Jersey.

Tectonophysics↗

USGS2008

Calculating wave-generated bottom orbital velocities from surface-wave parameters

Near-bed wave orbital velocities and shear stresses are important parameters in many sediment-transport and hydrodynamic models of the coastal ocean, estuaries, and lakes. Simple methods for estimating bottom orbital velocities from surface-wave statistics such as significant wave height and peak period often are inaccurate except in very shallow water. This paper briefly reviews approaches for estimating wave-generated bottom orbital velocities from near-bed velocity data, surface-wave spectra, and surface-wave parameters; MATLAB code for each approach is provided. Aspects of this problem have been discussed elsewhere. We add to this work by providing a method for using a general form of the parametric surface-wave spectrum to estimate bottom orbital velocity from significant wave height and peak period, investigating effects of spectral shape on bottom orbital velocity, comparing methods for calculating bottom orbital velocity against values determined from near-bed velocity measurements at two sites on the US east and west coasts, and considering the optimal representation of bottom orbital velocity for calculations of near-bed processes. Bottom orbital velocities calculated using near-bed velocity data, measured wave spectra, and parametric spectra for a site on the northern California shelf and one in the mid-Atlantic Bight compare quite well and are relatively insensitive to spectral shape except when bimodal waves are present with maximum energy at the higher-frequency peak. These conditions, which are most likely to occur at times when bottom orbital velocities are small, can be identified with our method as cases where the measured wave statistics are inconsistent with Donelan's modified form of the Joint North Sea Wave Project (JONSWAP) spectrum. We define the 'effective' forcing for wave-driven, near-bed processes as the product of the magnitude of forcing times its probability of occurrence, and conclude that different bottom orbital velocity statistics may be appropriate for different problems. ?? 2008 Elsevier Ltd.

Computers & Geosciences↗

USGS2021

Simulated global coastal ecosystem responses to a half-century increase in river nitrogen loads

Coastal ecosystems are increasingly threatened by anthropogenic stressors such as harmful algal blooms and hypoxia projected to intensify through the combined effects of eutrophication and warming. As a major terrestrial nitrogen (N) source to the ocean, rivers play a critical role in shaping both coastal and global biogeochemical cycling. Combining an enhanced-resolution (1/4°), global ocean physical-biogeochemical model with dynamic river inputs, we estimate that elevated river nitrogen loads alone resulted in an increase of 16.6 Tg (+5.8%) in the global coastal nitrogen inventory (CNI) over the half century between 1961 and 2010. This change was accompanied by increases in coastal net primary productivity (NPP, +4.6%) and benthic detrital flux (BDF, +7.3%), the latter of which is indicative of an overall higher oxygen demand in coastal sediments. After normalization by area, the ecosystems most sensitive to added river nitrogen (g N m-2 yr-1) were those with long residence times and strong nitrogen limitation. While even enhanced-resolution global models remain limited in their capacity to resolve near-shore responses, these basic sensitivity factors provide two relevant axes for frameworks assessing the comparative susceptibility of globally distributed coastal ecosystems to enhanced nitrogen loading, and the effectiveness of mitigation strategies.

Geophysical Research Letters↗

USGS2020

Meteotsunamis triggered by tropical cyclones

Tropical cyclones are one of the most destructive natural hazards and much of the damage and casualties they cause are flood-related. Accurate characterization and prediction of total water levels during extreme storms is necessary to minimize coastal impacts. While meteotsunamis are known to influence water levels and to produce severe consequences, they have been disregarded during tropical cyclones. This study demonstrates that meteotsunami waves commonly occur during tropical cyclones, and that they can significantly contribute to total water levels. We have discovered that the most extreme meteotsunami events were triggered by inherent features of the structure of tropical cyclones: inner and outer spiral rainbands. While outer distant spiral rainbands produced single-peak meteotsunami waves, inner spiral rainbands triggered longer lasting (~12 hours) wave trains on the front side of the tropical cyclones. We use an idealized coupled ocean-atmosphere-wave numerical model to analyze TC meteotsunami generation and propagation mechanisms.

Nature Communications↗

USGS2005

Recent water temperature trends in the lower Klamath River, California

Elevated water temperatures have been implicated as a factor limiting the recovery of anadromous salmonids in the Klamath River basin. This article reviews evidence of a multi-decade trend of increasing temperatures in the lower main-stem Klamath River above the ocean and, based on model simulations, finds a high probability that water temperature has been increasing by approximately 0.5??C/decade (95% confidence interval [Cl] = 0.42-0.60??C/decade) since the early 1960s. The season of high temperatures that are potentially stressful to salmonids has lengthened by about 1 month over the period studied, and the average length of main-stem river with cool summer temperatures has declined by about 8.2 km/decade. Water temperature trends seem unrelated to any change in main-stem water availability but are consistent with measured basinwide air temperature increases. Main-stem warming may be related to the cyclic Pacific Decadal Oscillation, but if this trend continues it might jeopardize the recovery of anadromous salmonids in the Klamath River basin. ?? Copyright by the American Fisheries Society 2005.

North American Journal of Fisheries Management↗

USGS2024

Assessing habitat use and population dynamics of fisheries resources at Kaloko Fishpond

Throughout Hawai'i, fishponds are considered by their local communities as important cultural touchstones, a source of local, sustainably produced food, and an important component to the development of community-based management for nearshore fisheries. Within Kaloko-Honokōhau National Historical Park, the restoration of Kaloko Fishpond for traditional aquaculture management is a goal of both the National Park Service (NPS) and Hui Kaloko-Honokōhau, a community-based group of kia'i, i.e., caretakers and native Hawaiian cultural practitioners. However, existing data on the demographics and condition of the fish populations within the pond, and the fish-habitat quality are poor to non-existent. Therefore, the objectives of this study were to: catalog fish species composition and distribution in the pond; estimate the abundance of focal species/taxonomic groups; and evaluate the occupancy patterns of the invasive algae Acanthophora spicifera and Upside-down Jellyfish Cassiopea andromeda. As part of these objectives, a survey protocol and analysis framework were designed and evaluated to ensure that the NPS and community group would be able to refine and implement them to continue their monitoring efforts. We conducted dual-observer shore-based visual surveys multiple times per week during September-October 2020 and April-September 2022. A total of 41 species/taxonomic groups were recorded over the course of the surveys. The largest number of species/taxonomic groups were observed at survey stations located on or near the kuapā, or wall separating the fishpond from the ocean. N-mixture models fitted to the data estimated a total population of 353 – 392 mullets, 134 – 192 flagtails (āholehole), and 189 – 277 Milkfish (Awa) Chanos chanos occurring within the 1.2-ha portion of Kaloko Fishpond that could be surveyed visually from the shoreline. Multi-season occupancy models fitted to the surveyed presence of A. spicifera and Upside-down Jellyfish indicted sites throughout most of the pond exhibited moderate and consistent occupancy (ψ = 0.30 – 0.40) throughout much of the pond, except for the northeast corner of the pond (Kaloko Iki) where colonization rates were lower and extinction rates higher than other areas within Kaloko. The visual survey method developed for this study provides a low-cost and effective starting point for the development of methodology that can be used both by NPS personnel and volunteers from the community group. However, we were only able to estimate fish populations for approximately 24% of the area of Kaloko Fishpond with this method. Given that the deeper areas of Kaloko Fishpond are completely inaccessible to the visual survey method used, generating population estimates for the entire pond based on the parameters estimated in the current study is not recommended without further investigation into fish movement and habitat use. Various means to refine this protocol to better meet the needs and abilities of the NPS and community group are proposed.

Hawaii↗

USGS2016

A study of the 2015 Mw 8.3 Illapel earthquake and tsunami: Numerical and analytical approaches

The September 16, 2015 Illapel, Chile earthquake triggered a large tsunami, causing both economic losses and fatalities. To study the coastal effects of this earthquake, and to understand how such hazards might be accurately modeled in the future, different finite fault models of the Illapel rupture are used to define the initial condition for tsunami simulation. The numerical code Non-hydrostatic Evolution of Ocean WAVEs (NEOWAVE) is employed to model the tsunami evolution through the Pacific Ocean. Because only a short time is available for emergency response, and since the earthquake and tsunami sources are close to the coast, gaining a rapid understanding of the near-field run-up behavior is highly relevant to Chile. Therefore, an analytical solution of the 2 ? 1 D shallow water wave equations is considered. With this solution, we show that we can quickly estimate the run-up distribution along the coastline, to first order. After the earthquake and tsunami, field observations were measured in the surrounded coastal region, where the tsunami resulted in significant run-up. First, we compare the analytical and numerical solutions to test the accuracy of the analytical approach and the field observations, implying the analytic approach can accurately model tsunami run-up after an earthquake, without sacrificing the time necessary for a full numerical inversion. Then, we compare both with field run-up measurements. We observe the consistency between the two approaches. To complete the analysis, a tsunami source inversion is performed using run-up field measurements only. These inversion results are compared with seismic models, and are shown to capture the broad-scale details of those models, without the necessity of the detailed data sets they invert.

Illapel↗