Modeling individual growth reveals decreasing gray whale body length and correlations with ocean climate indices at multiple scales.

Changes in body size have been documented across taxa in response to human activities and climate change. Body size influences many aspects of an individual's physiology, behavior, and ecology, ultimately affecting life history performance and resilience to stressors. In this study, we developed an analytical approach to model individual growth patterns using aerial imagery collected via drones, which can be used to investigate shifts in body size in a population and the associated drivers. We applied the method to a large morphological dataset of gray whales (Eschrichtius robustus) using a distinct foraging ground along the NE Pacific coast, and found that the asymptotic length of these whales has declined since around the year 2000 at an average rate of 0.05-0.12 m/y. The decline has been stronger in females, which are estimated to be now comparable in size to males, minimizing sexual dimorphism. We show that the decline in asymptotic length is correlated with two oceanographic metrics acting as proxies of habitat quality at different scales: the mean Pacific Decadal Oscillation index, and the mean ratio between upwelling intensity in a season and the number of relaxation events. These results suggest that the decline in gray whale body size may represent a plastic response to changing environmental conditions. Decreasing body size could have cascading effects on the population's demography, ability to adjust to environmental changes, and ecological influence on the structure of their community. This finding adds to the mounting evidence that body size is shrinking in several marine populations in association with climate change and other anthropogenic stressors. Our modeling approach is broadly applicable across multiple systems where morphological data on megafauna are collected using drones.

A fungi hotspot deep in the ocean: explaining the presence of Gjaerumia minor in equatorial Pacific bathypelagic waters.

A plant parasite associated with the white haze disease in apples, the Basidiomycota Gjaerumia minor, has been found in most samples of the global bathypelagic ocean. An analysis of environmental 18S rDNA sequences on 12 vertical profiles of the Malaspina 2010 expedition shows that the relative abundance of this cultured species increases with depth while its distribution is remarkably different between the deep waters of the Pacific and Atlantic oceans, being present in higher concentrations in the former. This is evident from sequence analysis and a microscopic survey with a species-specific newly designed TSA-FISH probe. Several hints point to the hypothesis that G. minor is transported to the deep ocean attached to particles, and the absence of G. minor in bathypelagic Atlantic waters could then be explained by the absence of this organism in surface waters of the equatorial Atlantic. The good correlation of G. minor biomass with Apparent Oxygen Utilization, recalcitrant carbon and free-living prokaryotic biomass in South Pacific waters, together with the identification of the observed cells as yeasts and not as resting spores (teliospores), point to the possibility that once arrived at deep layers this species keeps on growing and thriving.

Fungal diversity in sediments of the eastern tropical Pacific oxygen minimum zone revealed by metabarcoding.

Oxygen minimum zones (OMZ) represent ~8% of the ocean, with the Pacific as the largest and top expanding area. These regions influence marine ecosystems, promoting anaerobic microbial communities. Nevertheless, only a fraction of microbial diversity has been studied, with fungi being the less explored component. So, herein we analyzed fungal diversity patterns in surface and subsurface sediments along a bathymetric transect using metabarcoding of the ITS1 region in the OMZ of the Mexican Pacific off Mazatlán. We identified 353 amplicon sequence variants (ASV), within the Ascomycota, Basidiomycota, and Rozellomycota. Spatial patterns evidenced higher alpha diversity in nearshore and subsurface subsamples, probably due to temporal fluctuations in organic matter inputs. Small-scale heterogeneity characterized the community with the majority of ASV (269 ASV) occurring in a single subsample, hinting at the influence of local biogeochemical conditions. This baseline data evidenced a remarkable fungal diversity presenting high variation along a bathymetric and vertical transects.

Production of sounds by squirrelfish during symbiotic relationships with cleaner wrasses.

Examples of symbiotic relationships often include cleaning mutualisms, typically involving interactions between cleaner fish and other fish, called the clients. While these cleaners can cooperate by removing ectoparasites from their clients, they can also deceive by feeding on client mucus, a behavior usually referred to as "cheating behavior" that often leads to a discernible jolt from the client fish. Despite extensive studies of these interactions, most research has focused on the visual aspects of the communication. In this study, we aimed to explore the role of acoustic communication in the mutualistic relationship between cleaner fishes and nine holocentrid client species across four regions of the Indo-Pacific Ocean: French Polynesia, Guam, Seychelles, and the Philippines. Video cameras coupled with hydrophones were positioned at various locations on reefs housing Holocentridae fish to observe their acoustic behaviors during interactions. Our results indicate that all nine species of holocentrids can use acoustic signals to communicate to cleaner fish their refusal of the symbiotic interaction or their desire to terminate the cooperation. These sounds were predominantly observed during agonistic behavior and seem to support visual cues from the client. This study provides a novel example of acoustic communication during a symbiotic relationship in teleosts. Interestingly, these vocalizations often lacked a distinct pattern or structure. This contrasts with numerous other interspecific communication systems where clear and distinguishable signals are essential. This absence of a clear acoustic pattern may be because they are used in interspecific interactions to support visual behavior with no selective pressure for developing specific calls required in conspecific recognition. The different sound types produced could also be correlated with the severity of the client response. There is a need for further research into the effects of acoustic behaviors on the quality and dynamics of these mutualistic interactions.

Annual phenology and migration routes to breeding grounds in western-central North Pacific sei whales.

The sei whale (Balaenoptera borealis) is an important species among baleen whales in the North Pacific and plays a significant role in the ecosystem. Despite the importance of this species, information regarding its migration patterns and breeding locations remains limited. To enhance the understanding of the phenology of North Pacific sei whales, we deployed satellite-monitored tags on these whales in the western and central North Pacific from 2017 to 2023. We fitted 55 sei whale tracks to a state-space model to describe the whales' seasonal movements at feeding grounds and their migratory behavior. The whales typically leave their feeding grounds between November and December, with migration pathways extending from off Japan to the west of the Hawaiian Islands. These southward transits converge in the waters of the Marshall Islands and north of Micronesia between 20° N and 7° N, which appear to be breeding grounds. After a brief stay at these breeding grounds, the whales migrate northward from January to February, reaching their feeding grounds around 30°N by March. To the best of our knowledge, this is the first study to present the phenology of feeding and breeding seasons and the migration pattern of North Pacific sei whales.

Polyethylene terephthalate (PET)-degrading bacteria in the pelagic deep-sea sediments of the Pacific Ocean.

Polyethylene terephthalate (PET) plastic pollution is widely found in deep-sea sediments. Despite being an international environmental issue, it remains unclear whether PET can be degraded through bioremediation in the deep sea. Pelagic sediments obtained from 19 sites across a wide geographic range in the Pacific Ocean were used to screen for bacteria with PET degrading potential. Bacterial consortia that could grow on PET as the sole carbon and energy source were found in 10 of the 19 sites. These bacterial consortia showed PET removal rate of 1.8%-16.2% within two months, which was further confirmed by the decrease of carbonyl and aliphatic hydrocarbon groups using attenuated total reflectance-Fourier-transform infrared analysis (ATR-FTIR). Analysis of microbial diversity revealed that Alcanivorax and Pseudomonas were predominant in all 10 PET degrading consortia. Meanwhile, Thalassospira, Nitratireductor, Nocardioides, Muricauda, and Owenweeksia were also found to possess PET degradation potential. Metabolomic analysis showed that Alcanivorax sp. A02-7 and Pseudomonas sp. A09-2 could turn PET into mono-(2-hydroxyethyl) terephthalate (MHET) even in situ stimulation (40 MPa, 10 °C) conditions. These findings widen the currently knowledge of deep-sea PET biodegrading process with bacteria isolates and degrading mechanisms, and indicating that the marine environment is a source of biotechnologically promising bacterial isolates and enzymes.

When species don't move together: Non-concurrent range shifts in Eastern Pacific kelp forest communities.

Species range shifts due to changing ocean conditions are occurring around the world. As species move, they build new interaction networks as they shift from or into new ecological communities. Typically, species ranges are modeled individually, but biotic interactions have been shown to be important to creating more realistic modeling outputs for species. To understand the importance of consumer interactions in Eastern Pacific kelp forest species distributions, we used a Maxent framework to model a key foundation species, giant kelp (Macrocystis pyrifera), and a dominant herbivore, purple sea urchins (Strongylocentrotus purpuratus). With neither species having previously been modeled in the Eastern Pacific, we found evidence for M. pyrifera expansion in the northern section of its range, with no projected contraction at the southern range edge. Despite its known co-occurrence with M. pyrifera, models of S. purpuratus showed a non-concurrent southern range contraction and a co-occurring northern range expansion. While the co-occurring shifts may lead to increased spatial competition for suitable substrate, this non-concurrent contraction could result in community wide impacts such as herbivore release, tropicalization, or ecosystem restructuring.

Bacteria on the foundational kelp in kelp forest ecosystems: Insights from culturing, whole genome sequencing and metabolic assays.

In coastal marine ecosystems, kelp forests serve as a vital habitat for numerous species and significantly influence local nutrient cycles. Bull kelp, or Nereocystis luetkeana, is a foundational species in the iconic kelp forests of the northeast Pacific Ocean and harbours a complex microbial community with potential implications for kelp health. Here, we report the isolation and functional characterisation of 16 Nereocystis-associated bacterial species, comprising 13 Gammaproteobacteria, 2 Flavobacteriia and 1 Actinomycetia. Genome analyses of these isolates highlight metabolisms potentially beneficial to the host, such as B vitamin synthesis and nitrogen retention. Assays revealed that kelp-associated bacteria thrive on amino acids found in high concentrations in the ocean and in the kelp (glutamine and asparagine), generating ammonium that may facilitate host nitrogen acquisition. Multiple isolates have genes indicative of interactions with key elemental cycles in the ocean, including carbon, nitrogen and sulphur. We thus report a collection of kelp-associated microbial isolates that provide functional insight for the future study of kelp-microbe interactions.

Molecular evidence for a new endemic species of Acartia (Copepoda, Calanoida) from the Southeast Pacific coast.

The loss of biodiversity in marine populations is one of the consequences of the increased events of extreme environmental conditions in the oceans, which can condition the persistence of populations to future scenarios of climate change. Therefore, it is extremely necessary to explore and monitor the genetic diversity of natural populations. In the Southeast Pacific Ocean (SEPO), specifically on the coast of Chile, the presence of the copepod Acartia tonsa has been indicated solely using morphological evidence, due to the absence of genetic information. In the present work, the genetic diversity, population structure and phylogenetic position within the genus Acartia, of populations identified morphologically as A. tonsa, was evaluated by amplification of the mitochondrial cytochrome c oxidase subunit I and nuclear marker 18 s. Our results showed that the populations identified as A. tonsa correspond to a new monophyletic group endemic to SEPO (GMYC = 1.00; PTP = 0.95). The populations showed moderate to high genetic diversity with an incipient structuring between populations and biogeographic zones. Our results suggest that despite the homogenizing effect of the Humboldt Current, isolation by distance and contrasting environmental conditions at different geographic scales have an important influence on the genetic diversity of zooplankton in the SEPO region.

Resistance of rocky intertidal communities to oceanic climate fluctuations.

A powerful way to predict how ecological communities will respond to future climate change is to test how they have responded to the climate of the past. We used climate oscillations including the Pacific Decadal Oscillation (PDO), North Pacific Gyre Oscillation, and El Niño Southern Oscillation (ENSO) and variation in upwelling, air temperature, and sea temperatures to test the sensitivity of nearshore rocky intertidal communities to climate variability. Prior research shows that multiple ecological processes of key taxa (growth, recruitment, and physiology) were sensitive to environmental variation during this time frame. We also investigated the effect of the concurrent sea star wasting disease outbreak in 2013-2014. We surveyed nearly 150 taxa from 11 rocky intertidal sites in Oregon and northern California annually for up to 14-years (2006-2020) to test if community structure (i.e., the abundance of functional groups) and diversity were sensitive to past environmental variation. We found little to no evidence that these communities were sensitive to annual variation in any of the environmental measures, and that each metric was associated with < 8.6% of yearly variation in community structure. Only the years elapsed since the outbreak of sea star wasting disease had a substantial effect on community structure, but in the mid-zone only where spatially dominant mussels are a main prey of the keystone predator sea star, Pisaster ochraceus. We conclude that the established sensitivity of multiple ecological processes to annual fluctuations in climate has not yet scaled up to influence community structure. Hence, the rocky intertidal system along this coastline appears resistant to the range of oceanic climate fluctuations that occurred during the study. However, given ongoing intensification of climate change and increasing frequencies of extreme events, future responses to climate change seem likely.

Interannual variability in marine debris accumulation on Hawaiian shores: The role of North Pacific Ocean basin-scale dynamics.

Community-based marine debris removal efforts on the Hawaiian Islands of Kaua'i and Hawai'i, spanning 2013-2022, provided large datasets and documented remarkable variations in annual amounts of debris, mainly from abandoned, lost and derelict fishing gear. To test the hypothesis that the influx of marine debris on Hawaiian shores is determined by the proximity of the North Pacific garbage patch, whose pattern changes under the control of large-scale ocean dynamics, we compared these observational data with the output of an oceanographic drift model. The high correlations between the total mass of debris collected and the model, ranging between r = 0.81 and r = 0.84, validate the attribution of the strong interannual signal to significant migrations of the garbage patch reproduced in the model experiments. Synchronous variations in marine debris fluxes on the two islands, separated by >500 km, confirm the large scale of the interannual changes in the North Pacific marine debris system.

Macrolitter and microplastics along the East Pacific coasts - A homemade problem needing local solutions.

The East Pacific (EP) region, especially the central and southern EP, has been fairly less studied than other world's regions with respect to marine litter pollution. This comprehensive literature review (257 peer-reviewed publications) showed that both macrolitter (mostly plastics) and microplastics tend to accumulate on EP shorelines. Moreover, they were also reported in all the other compartments investigated: sea surface, water column, seafloor and 'others'. Mostly local, land-based sources (e.g., tourism, poor waste management) were identified across the region, especially at continental sites from low and mid latitudes. Some sea-based sources (e.g., fisheries, long-distance drifting) were also identified at high latitudes and on oceanic islands, likely enhanced by the oceanographic dynamics of the EP that affect transport of floating litter. Our results suggest that effective solutions to the problem require local and preventive strategies to significantly reduce the levels of litter along the EP coasts.

MOLECULAR PHYLOGENY OF THE LEECH GENUS PONTOBDELLA (HIRUDINIDA: PISCICOLIDAE) WITH NOTES ON PONTOBDELLA CALIFORNIANA AND PONTOBDELLA MACROTHELA.

Leech specimens of the genus Pontobdella (Hirudinida: Piscicolidae) were found off the coast of the state of Oaxaca (Pacific) as well as in Veracruz and Tabasco (Gulf of Mexico), Mexico. Based on the specimens collected in Oaxaca, a redescription of Pontobdella californiana is provided, with emphasis on the differences in the reproductive organs with the original description of the species. In addition, leech cocoons assigned to P. californiana were found attached to items hauled by gillnets and studied using scanning electron microscopy and molecular approaches. Samples of Pontobdella macrothela were found in both Pacific and Atlantic oceans, representing new geographic records. The phylogenetic position of P. californiana is investigated for the first time, and with the addition of Mexican samples of both species, the phylogenetic relationships within Pontobdella are reinvestigated. Parsimony and maximum-likelihood phylogenetic analysis were based on mitochondrial (cytochrome oxidase subunit I [COI] and 12S rRNA) and nuclear (18S rRNA and 28S rRNA) DNA sequences. Based on our results, we confirm the monophyly of Pontobdella and the pantropical distribution of P. macrothela with a new record in the Tropical Eastern Pacific.

PREVALENCE OF OPECOELID TREMATODES (TREMATODA: DIGENEA: OPECOELIDAE) IN THE WOOLLY SCULPIN CLINOCOTTUS ANALIS FROM THE PACIFIC COAST OF THE BAJA CALIFORNIA PENINSULA.

Prevalence values for 3 digenean species of the family Opecoelidae were calculated during a half-year period from 20 individuals per month of the fish species Clinocottus analis, collected from an intertidal environment of Baja California, Mexico. Trematode species recovered were Opecoelus adsphaericus, Opecoelus cameroni, and Opecoelus pacificus. Of these, only O. adsphaericus was present throughout the study, whereas O. pacificus and O. cameroni were recorded for 1 or 2 mo, respectively, exhibiting relatively low prevalence values. The decrease in prevalence of O. adsphaericus coincides with the appearance of O. pacificus and O. cameroni; these last 2 species were found only in the largest hosts, and their presence represents new records for Mexican marine fauna.

Genomic characteristics of nine Nitrospirota metagenome-assembled genomes in deep-sea sediments from East Pacific polymetallic nodules zone.

Previously studies have reported that MAGs (Metagenome-assembled genomes) belong to "Candidatus Manganitrophaceae" of phylum Nitrospirota with chemolithoautotrophic manganese oxidation potential exist in freshwater and hydrothermal environments. However, Nitrospirota members with chemolithoautotrophic manganese oxidation potential have not been reported in other marine environments. Through metagenomic sequencing, assembly and binning, nine metagenome-assembled genomes belonging to Nitrospirota are recovered from sediment of different depths in the polymetallic nodule area. Through the key functional genes annotation results, we find that these Nitrospirota have limited potential to oxidize organic carbon because of incomplete tricarboxylic acid cycle and most of them (6/9) have carbon dioxide fixation potential through different pathway (rTCA, WL or CBB). One MAG belongs to order Nitrospirales has the potential to use manganese oxidation to obtain energy for carbon fixation. In addition to manganese ions, the oxidation of inorganic nitrogen, sulfur, hydrogen and carbon monoxide may also provide energy for the growth of these Nitrospirota. In addition, different metal ion transport systems can help those Nitrospirota to resist heavy metal in sediment. Our work expands the understanding of the metabolic potential of Nitrospirota in sediment of polymetallic nodule region and may contributes to promoting the study of chemolithoautotrophic manganese oxidation.

Seasonal variability in energetic value of Crangon alaskensis and effects of marine heatwaves in the Northeast Pacific Ocean.

Research cruises were conducted to sample the invertebrate community along the shelf off the central coast of Oregon from 2010 to 2018. A large marine heatwave (MHW) hit the northeast Pacific in fall 2014 and persisted locally through 2015. Here, we assessed the caloric content changes of Crangon alaskensis (a common sandy shrimp) before, during, and after the 2014-2015 MHW. We found significant reductions in the caloric density of shelf populations of C. alaskensis during summer 2015. Oceanographic indices like the Biologically Effective Upwelling Transport Index (BEUTI) and the Pacific Decadal Oscillation (PDO) had greater predictive power for caloric density and biomass than in situ conditions, although bottom temperature and dissolved oxygen were also significantly correlated with caloric density. Caloric density of C. alaskensis was highest in 2018, indicating favorable conditions after the intense MHW of 2014-2015 allowed the caloric density to rebound.

Variability in the relationship between light scattering and chlorophyll a concentration in oligotrophic tropical regions of the Western Pacific Ocean.

It is important to determine the relationship between the concentration of chlorophyll a (Chla) and the inherent optical properties (IOPs) of ocean water to develop optical models and algorithms that characterize the biogeochemical properties and estimate biological pumping and carbon flux in this environment. However, previous studies reported relatively large variations in the particulate backscattering coefficient (bbp(λ)) and Chla from more eutrophic high-latitude waters to clear oligotrophic waters, especially in oligotrophic oceanic areas where these two variables have little covariation. In this study, we examined the variability of bbp(λ) and Chla in the euphotic layer in oligotrophic areas of the tropical Western Pacific Ocean and determined the sources of these variations by reassessment of in-situ measurements and the biogeochemical-argo (BGC-Argo) database. Our findings identified covariation of bbp(λ) and Chla in the water column below the deep Chla maximum (DCM) layer, and indicated that there was no significant correlation relationship between bbp(λ) and Chla in the upper layer of the DCM. Particles smaller than 3.2 µm that were in the water column above the DCM layer had a large effect on the bbp(λ) in the vertical profile, but particles larger than 3.2 µm and smaller than 10 µm had the largest effect on the bbp(λ) in the water column below the DCM layer. The contribution of non-algal particles (NAPs) to backscattering is up to 50%, which occurs in the water depth of 50 m and not consistent with the distribution of Chla. Phytoplankton and NAPs were modeled as coated spheres and homogeneous spherical particles to simulate the bbp(λ) of the vertical profile by Aden-Kerker method and Mie theory, and the results also indicated that the backscattering caused by particles less than 20 µm were closer to the measured data when they were below and above the DCM layer, respectively. This relationship also reflects the bbp(λ) of particles in the upper water was significantly affected particle size, but bbp(λ) in the lower water was significantly affected by Chla concentration. This effect may have relationship with phytoplankton photoacclimation and the relationship of a phytoplankton biomass maximum with particle size distribution in the water column according to the previous relevant studies. These characteristics also had spatial and seasonal variations due to changes of Chla concentration at the surface and at different depths. There was mostly a linear relationship between Chla and bbp(700) during winter. During other seasons, the relationship between these two variables was better characterized by a power function (or a logarithmic function) in the lower layer of the DCM. The spatial and vertical relationships between the bbp(λ) and Chla and the corresponding variations in the types of particles described in this study provide parameters that can be used for accurate estimation of regional geochemical processes.

Flagellimonas baculiformis sp. nov. and Flagellimonas crocea sp. nov., isolated from surface seawater of the Pacific Ocean.

Two Gram-negative, non-spore-forming, non-motile, non-flagellated bacteria, designated strains D6T and DH64T, were isolated from surface water of the Pacific Ocean. For strain D6T, growth occurred at 10-40 °C, pH 5.5-9.0 and in the presence of 0-8.0 % NaCl (w/v). For strain DH64T, growth occurred at 10-40 °C, pH 5.5-8.5 and in the presence of 0.5-8.0 % NaCl (w/v). Phylogenetic analysis based on 16S rRNA gene sequences indicated that strains D6T and DH64T both belonged to the genera Flagellimonas, with the highest sequence identities to Flagellimonas taeanensis JCM 17757T (98.2 %) and Flagellimonas marinaquae JCM 11811T (98.6 %), respectively. The 16S rRNA gene sequence identity between strains D6T and DH64T was 95.9 %. The average amino acid identity and digital DNA-DNA hybridization values between the two strains and the nearest phylogenetic neighbours were 66.7-93.3 % and 16.1-38.5 %, respectively. The major respiratory quinone of both strains was menaquinone-6. The major polar lipid was phosphatidylethanolamine. The major fatty acids were identified similarly as iso-C15 : 1 G, iso-C15 : 0 and iso-C17 : 0 3-OH. The genomic G+C contents of strains D6T and DH64T were determined to be 45.5 and 42.6 mol%, respectively. The combined genotypic and phenotypic data show that the strains represent two novel species within genera Flagellimonas, for which the names Flagellimonas baculiformis sp. nov. and Flagellimonas crocea sp. nov. are proposed, with type strains D6T (=MCCC M28982T=KCTC 92604T) and DH64T (=MCCC M28986T=KCTC 92975T).

Habitat use of loggerhead (Caretta caretta) and green (Chelonia mydas) turtles at the northern limit of their distribution range of the Northwest Pacific Ocean.

Verifying habitats, including the foraging and nesting areas for sea turtles, enables an understanding of their spatial ecology and successful planning of their conservation and management strategies. Recently, the observation frequency and bycatch of loggerhead (Caretta caretta) and green (Chelonia mydas) turtles have increased in the northern limit of their distribution range, in the northern part of the East China Sea and East (Japan) Sea. We conducted satellite tracking to investigate the habitat use of seven loggerhead and eight green turtles from June 2016 to August 2022 in this area, where little is known about their spatial ecology. We applied a 50 percent volume contour method to determine their main foraging areas and analyzed 6 environmental variables to characterize their habitats. Loggerhead turtles mainly stayed in and used the East China Sea as a foraging area during the tracking period, while two individuals among them also used the East Sea as a seasonal foraging area. Most green turtles also used the East China Sea as a foraging area, near South Korea and Japan, with one individual among them using the lower area of the East Sea as a seasonal foraging area. Notably, one green turtle traveled to Hainan Island in the South China Sea, a historical nesting area. Our results showed that the two sea turtle species included the East Sea as a seasonal foraging area, possibly owing to the abundance of food sources available, despite its relatively lower sea temperature. Considering that loggerhead and green sea turtles were observed using the northern part of the East China Sea and East Sea more frequently than previously known and that the sea temperature gradually increases due to climate change, conservation and management activities are required for sea turtles in these areas.

Three decades of nearshore surveys reveal long-term patterns in gray whale habitat use, distribution, and abundance in the Northern California Current.

The nearshore waters of the Northern California Current support an important seasonal foraging ground for Pacific Coast Feeding Group (PCFG) gray whales. We examine gray whale distribution, habitat use, and abundance over 31 years (1992-2022) using standardized nearshore (< 5 km from shore) surveys spanning a large swath of the PCFG foraging range. Specifically, we generated density surface models, which incorporate detection probability into generalized additive models to assess environmental correlates of gray whale distribution and predict abundance over time. We illustrate the importance of coastal upwelling dynamics, whereby increased upwelling only yields higher gray whale density if interspersed with relaxation events, likely because this combination optimizes influx and retention of nutrients to support recruitment and aggregation of gray whale prey. Several habitat features influence gray whale distribution, including substrate, shelf width, prominent capes, and river estuaries. However, the influence of these features differs between regions, revealing heterogeneity in habitat preferences throughout the PCFG foraging range. Predicted gray whale abundance fluctuated throughout our study period, but without clear directional trends, unlike previous abundance estimates based on mark-recapture models. This study highlights the value of long-term monitoring, shedding light on the impacts of variable environmental conditions on an iconic nearshore marine predator.