Bridging From Monitoring to Solutions-Based Thinking: Lessons From CalCOFI for Understanding and Adapting to Marine Climate Change Impacts

Multidisciplinary, integrated ocean observing programs provide critical data for monitoring the effects of climate change on marine ecosystems. California Cooperative Oceanic Fisheries Investigations (CalCOFI) samples along the US West Coast and is one of the world’s longest-running and most comprehensive time series, with hydrographic and biological data collected since 1949. The pairing of ecological and physical measurements across this long time series informs our understanding of how the California Current marine ecosystem responds to climate variability. By providing a baseline to monitor change, the CalCOFI time series serves as a Keeling Curve for the California Current. However, challenges remain in connecting the data collected from long-term monitoring programs with the needs of stakeholders concerned with climate change adaptation (i.e. resource managers, policy makers, and the public), including for the fisheries and aquaculture sectors. We use the CalCOFI program as a case study to ask: how can long-term ocean observing programs inform ecosystem based management efforts and create data flows that meet the needs of stakeholders working on climate change adaptation? Addressing this question and identifying solutions requires working across sectors and recognizing stakeholder needs. Lessons learned from CalCOFI can inform other regional monitoring programs around the world, including those done at a smaller scale in developing countries.


INTRODUCTION
The US West Coast is one of the best studied marine regions of the world due to a high concentration of academic institutions, investment in fisheries, and political support. Consequently, the data-rich environment provides an excellent opportunity to explore how oceanographic monitoring programs can support ecosystem management strategies that are adaptive to climate change. Ideally, more comprehensive oceanographic monitoring should lead to more climate-resilient management (Skern-Mauritzen et al., 2016;Davidson et al., 2019), however, there are unique challenges in linking physical and biological oceanographic data to management needs (Field and Francis, 2006). These include, for example, difficulties in linking meaningful climate indices to fisheries trends to inform management decisions (Myers, 1998;Smith et al., 2007;Pitcher et al., 2009), and in creating data products that support the needs of a diverse community of stakeholders (Rayner et al., 2019).
One reason the US West Coast is so well studied is due to the long-running California Cooperative Oceanic Fisheries Investigations (CalCOFI) Program. CalCOFI began in the late 1940s in response to the collapse of the economically important sardine fishery on the US West Coast (Hewitt, 1988;Scheiber, 1990). Competing hypotheses argued for overfishing versus environmental change as the dominant driver responsible for the sardine decline, and the CalCOFI Program was founded to investigate the relationship between oceanographic conditions and sardine abundance (Hewitt, 1988). From the beginning, CalCOFI represented a unique partnership between academic, federal, and state partners (Scheiber, 1990), which now consists of the Scripps Institution of Oceanography at UCSD, federal scientists at NOAA's National Marine Fisheries Service (NMFS), and scientists working for the California Department of Fish and Wildlife (Ohman and Venrick, 2003;McClatchie, 2014). Since 2004, CalCOFI also contributes to the California Current Ecosystem Long Term Ecological Research (CCE LTER) Network (Ohman and Hobbie, 2008).
While initially focused on sardines, the scientists who founded CalCOFI had the vision to create an ecosystem monitoring program at a time when ecosystem research was in its infancy (Scheiber, 1997), thus beginning one of the most comprehensive ecosystem monitoring programs that extends the full length of the Keeling Curve (Keeling and Keeling, 2017) -a time series that is fundamental to the scientific understanding of increasing greenhouse gases and anthropogenic climate change. Presently, the overarching objective of CalCOFI is to understand the effects of longterm changes in the California Current Ecosystem (CCE) (Ohman and Venrick, 2003).
The goal of this study is to identify how the CalCOFI Program can best support stakeholder needs for understanding climate change impacts and developing climate change adaptation strategies. To do this, we first review the data that are collected by the CalCOFI Program (1949-current) and synthesize how these data have informed our understanding of the dynamics and climate sensitivity of the CCE. In doing so, we also examine the CalCOFI record for climate relevant variables and compare changes observed over the last 70 years to changes predicted by climate models for the future. We then consider how CalCOFI data are currently used to support ecosystem and fisheries management decision-making, explore these within the context of adaptation to climate change, and consider how CalCOFI data can be used by a diversity of stakeholders to support climate change adaptation needs in the future.

What Data Does the CalCOFI Program Collect?
At present, CalCOFI conducts quarterly ship-based surveys in a grid design at fixed stations in the Southern California Current Ecosystem (SCCE) and collects both hydrographic and biological samples. The earliest samples were collected in 1949, and the sampling protocols were standardized by the early 1950s, allowing for comparisons across the 70 years time series. Although the survey region has expanded and contracted, the "core" sampling grid, consistently maintained throughout the time series, is composed of 66 stations along 6 grid lines (Figure 1). This area ranges from the US-Mexico Border to slightly north of Point Conception, California, and extends over 500 km offshore. In 2004, nine additional coastal stations were added to the "core" CalCOFI grid, as part of the Southern California Coastal Ocean Observing System (SCCOOS), to extend sampling further inshore. During the winter and spring cruises, the survey grid typically expands further north to the San Francisco Bay region (Figure 1) to capture most of the sardine spawning spatial domain. The extended CalCOFI grid includes 113 stations, 16 of which fall in California state waters, 89 in U.S. federal waters, and 8 are located in international waters outside of the U.S. exclusive economic zone (Figure 1). Due to fluctuating funding between 1969 and 1983, for this period quarterly cruises were only conducted on a triennial basis. Line 90 represents the best sampled CalCOFI line, and spring and summer samples are more complete within the CalCOFI time-series than winter or fall (Peabody et al., 2018).
A suite of hydrographic and biological data are collected at CalCOFI stations and during transit between stations. Hydrographic bottle samples for temperature, salinity, oxygen, and phosphate were collected since 1949 at specific depths (up to 500 m). Nutrient analyses expanded to include silicate, nitrate and nitrite for bottle samples starting in 1961, and then chlorophyll in 1973, dissolved inorganic carbon (DIC) in 1983, and primary production using 14 C incubations in 1984. CTD profiling began in 1993, and at each station an SBE 911plus CTD is deployed to measure water column temperature, salinity, dissolved oxygen, chlorophyll a, and photosynthetically active radiation up to a depth of ∼515 m. The CTD is an oceanographic sensor package that is surrounded by a rosette of Niskin bottles for discreet water sampling and is deployed vertically through the water column to determine the hydrographic profiles of elements of interest (Williams, 2009). An SBE 18 pH sensor was added to the CTD sensor package starting in 2009.
Biological sampling at stations is conducted predominantly with nets. Vertical, oblique and surface tows are routinely collected. Oblique tows are a consistent feature of CalCOFI, and a tow is conducted at each station to sample the plankton. In 1978, bongo nets replaced ring nets to improve sampling, and an adjustment is applied for abundance estimates for some FIGURE 1 | Map of 113 CalCOFI stations as relates to U.S. jurisdictional boundaries. In CalCOFI stations are shown based on if they occur in state (orange), federal (blue), or international (black) waters. The US exclusive economic zone is shown with a black outline and California state waters are shown in gray surrounding land (tan). Circles represent the 66 station "core grid," squares represent the 38 extended grid stations, and the 9 nearshore stations are represented with asterisks. CalCOFI lines (horizontal) and stations (vertical) are indicated by the labels on the left of the grid and at the top and bottom of the grid, respectively. Not all station labels are shown. A scale bar is provided for spatial reference.
fish and zooplankton taxa to correct for changes in catchability (Ohman and Lavaniegos, 2002;Thompson et al., 2017b). A flow meter is attached to the net allowing for a measurement of sampling volume as the bongo net is towed obliquely from 212 m to the surface at 1-2 knots. Zooplankton displacement volume is measured for each sample, and all ichthyoplankton are removed. Ichthyoplankton are then identified to lowest possible taxon and enumerated. The ichthyoplankton time series is one of the best studied data products from the CalCOFI Program and includes over 400 species of fish. In recent years, certain invertebrate species of commercial importance have also been sorted and time series are available since 1997 for Market Squid (Doryteuthis opalescens) paralarvae and since 2007 for Spiny Lobster (Palinuridae spp.) phyllosoma. Ichthyoplankton and zooplankton from each tow are preserved in 5% formalin or 95% ethanol and kept in collections. Genetic barcoding has been conducted to identify species such as some rockfishes from ethanol-preserved samples that are difficult or impossible to identify based solely on morphology (Thompson et al., , 2017a. Additionally, zooplankton are sampled using vertically towed PairoVET nets deployed to 70 m depth, and manta nets (1977-present), which are towed at the surface. Information about net sampling procedures are presented in Kramer et al. (1972) and Smith and Richardson (1977). All hydrographic 1 and ichthyoplankton data 2 are quality-controlled and publicly available; invertebrate zooplankton data are partially available 3 .
Larval abundances from the CalCOFI ichthyoplankton timeseries are used to estimate or index adult spawning stock biomass. Egg and larval abundance data have been used routinely in conjunction with trawl surveys to estimate the spawning stock of Northern Anchovy (Engraulis mordax) and Pacific Sardine (Sardinops sagax) via the daily egg-production method (Lasker, 1985). In the absence of additional trawl data, larval abundance data alone have been used to provide a coarser index of adult spawning stock biomass for Northern Anchovy and Pacific Sardine , rockfish (Sebastes spp.; Moser et al., 2000), California Halibut (Paralichthys californicus; Moser and Watson, 1990), and for other rocky nearshore taxa (Moser et al., 2000). The fact that larval fish are predominantly sampled during the preflexion stage when net-sampler avoidance is minimal and post-spawning losses to mortality are still low relative to later larval stages makes them an effective index of spawning stock biomass . However, ichthyoplankton counts likely do not correlate with recruitment success due to high and variable mortality at the larval stage. CalCOFI samples are well suited to document dynamics of hundreds of species of fishes because although adult populations reside in different habitats (e.g., pelagic, mesopelagic, benthic), larvae of almost every species are found in the upper 200 m of the water column (Moser, 1996;Moser and Watson, 2006). Thus, the program focuses its efforts on pelagic sampling and no benthic samples are collected.
Since 2014, CalCOFI has incorporated genomic sampling through collaboration with scientists at the J. Craig Venter Institute, La Jolla, CA, United States. Seawater samples at specific stations along the CalCOFI grid are collected and sequenced for RNA and DNA to study the microbial and eukaryotic community composition and gene expression in the water column. This is called the NOAA CalCOFI Genomics Project.
A series of underway samples are also collected on CalCOFI cruises which include ADCP measurements, meteorological data, and marine mammal and seabird observations. Since 1996, the Continuous Underway Fish Egg Sampler (CUFES) (Checkley et al., 1997(Checkley et al., , 2000 has been used to collect continuous underway data from which Northern Anchovy, Pacific Sardine, Jack Mackerel (Trachurus symmetricus), Pacific Hake (Merluccius productus), and squid (Decapodiformes spp.) eggs are identified. Additional autonomous tools have also expanded the scope of CalCOFI sampling. Since 2005, underwater spray gliders collect hydrographic data, including temperature and salinity, along CalCOFI lines, which have included lines 66.7, 80, 90, and 93 (Davis et al., 2008). Dissolved oxygen sensors are also now being deployed on spray gliders. Additional details on the CalCOFI sampling methodology that are outside of the scope of this study can be found in McClatchie (2014).

Looking Back at the Past 70 Years and Forward to the Future
To examine the CalCOFI record in relation to future changes predicted by climate models, we first compare the length and trends of the CalCOFI time series for climate relevant variables to that of the Keeling Curve (Figure 2). The Keeling Curve is a record of atmospheric carbon dioxide measurements taken at the Mauna Loa Observatory in Hawaii since 1958, and has been fundamental to tracking anthropogenic greenhouse gas forcing and the subsequent relationship to climate change (Keeling and Keeling, 2017). For the CalCOFI time series, we examine trends in winter sea surface temperature (SST) ( Figure 2B) and winter dissolved oxygen concentration at 200 m ( Figure 2D) across the core 66-station CalCOFI grid. For ocean acidification-relevant measurements, we synthesize a timeline of existing measurements ( Figure 2C). As an example of the ichthyoplankton dataset, trends in the annual mean densities of Panama Lightfish, Vinciguerria lucetia, which is a warm-affinity mesopelagic species, are calculated across the core 66-station CalCOFI grid ( Figure 2E).
We then compare climate trends observed from the CalCOFI dataset with expected warming predicted from a climate model. We use the Large Ensemble Community Project (LENS) (Kay et al., 2015) to look at predicted changes in annual mean SST, pH, and oxygen concentration for the CCE under the "business as usual" emissions scenario (RCP 8.5). LENS includes a 40member ensemble of fully-coupled Community Earth System Model version 1 (CESM1) simulations for the period 1920-2100. From 1920 to 2005, the simulations use historic emissions data and after 2005 use RCP 8.5 projected emissions (Kay et al., 2015). We selected 1950-2100 as the period of interest to capture the period which corresponds to the CalCOFI time series.

Identifying Management and Stakeholder Use of CalCOFI Data
To evaluate how CalCOFI data currently informs fisheries and ecosystem management, we examine the fisheries and ecosystem management structures for the US West Coast, including federal and state policies, and identify instances where CalCOFI data are used for informing management. We also examine the spatial overlap between the CalCOFI grid and relevant management jurisdictions including state, federal, and international waters, and different marine conservation areas (California Marine Protected Areas, National Marine Sanctuaries, and rockfish conservation areas). We then identify stakeholder user groups of CalCOFI data, by considering both known users (i.e., those who have published CalCOFI data or have actively used it for their organization) and potential users (i.e., new organizations or groups that have not used CalCOFI data but whose activities overlap with the CalCOFI grid). Based on their known activities, we then consider the specific needs of different stakeholder groups related to questions of climate change impacts and adaptation, and where possible identify opportunities for better aligning stakeholder needs with CalCOFI monitoring.
By evaluating the climate change-relevant data collected by an ocean observing program and considering it within the landscape of ocean management policies and spatial overlap with regional stakeholders, other oceanographic observing programs can follow this process to better link their basic monitoring data to stakeholder needs relating to climate change adaptation. shows the increase in atmospheric CO 2 at Mauna Loa through the twentieth and twenty-first century and is an essential record for understanding climate forcing. The CalCOFI record (B-E) extends a similar length and is an essential record of how the ocean responds to climate forcing through hydrographic and biological changes. Unlike the atmospheric record, the oceanic record is noisy and shows pronounced multi-decadal natural variability. Trends in winter sea surface temperature (SST) across the core 66-station CalCOFI grid are shown in (B), the extent of all carbonate chemistry measurements taken on CalCOFI cruises is shown in (C), trends in winter dissolved oxygen at 200 m across the core 66-station CalCOFI grid are shown in (D), and trends in the annual mean density of Panama Lightfish, a warm-affinity mesopelagic species, across the core 66-station CalCOFI grid is shown in (E). Blue trend lines show smoothed conditional means.
One of the dominant modes of low-frequency variability (∼30 years) that affects the CCE is the Pacific Decadal Oscillation (PDO) which is based on variability of sea surface temperatures in the North Pacific (Mantua et al., 1997). Negative phases of the PDO are associated with colder, more productive overall conditions in the CCE, and positive phases are associated with warmer, less productive conditions. The North Pacific Gyre Oscillation (NPGO) is also a dominant mode of lowfrequency variability and is significantly correlated with multiple oceanographic conditions in the CalCOFI region, including upwelling strength, nutrient fluxes, and sea surface salinity anomalies . Higher frequency variability (2-7 years) is driven by the El Niño Southern Oscillation, which was first recognized as a Pacific-wide phenomenon at a CalCOFI Symposium (Sette and Isaacs, 1960). El Niño conditions are characterized by warm, less productive conditions in the CCE, while La Niña conditions are characterized by colder, more productive conditions. Stronger and more frequent El Niño conditions occur during positive PDO phases, while La Niña conditions are more pronounced and frequent during negative PDO phases. The CalCOFI record captures a negative PDO phase which lasted from 1950 to 1976 and a positive PDO phase that commenced after 1976. Although it was speculated that the PDO returned to a cool phase following a strong La Niña in 1999 (Zwolinski and Demer, 2012), the string of mostly negative PDO values between 1998 and 2003 now appear to have been more of a transient event rather than an actual regime shift (Bond et al., 2003). Consistent with this, there was no evidence of a prolonged change in the CalCOFI ichthyoplankton assemblage after 1998 (Peabody et al., 2018) and record high water temperatures occurred between 2014 and 2016 . Therefore, many studies consider the CalCOFI record in terms of the "cold" (1950-1976), "warm" (1977-1998), and present (1999-present) period.
Anthropogenic climate change is a directional change on which these natural modes of variability are layered, and the high amount of natural variability in the system makes it difficult to discern the anthropogenic signal (Figure 2). However, projections of the climate in the North Pacific for the twentyfirst century predict an increasing importance of anthropogenic warming over internal modes of variability (Bond et al., 2003;Di Lorenzo et al., 2008;Bonfils and Santer, 2011). Here, we briefly review how the CalCOFI record has informed our understanding of how the SCCE responds to different timescales of climate variability.

Warming
Despite strong inter-annual variability, temperatures within the CalCOFI region have shown a secular warming trend from 1949 to 2000, with stronger warming observed inshore than offshore (McClatchie, 2014). From 1949 to 2000, sea surface temperatures (SST) warmed by 1.3 • C or 0.025 • C per year on average (Di Lorenzo et al., 2005). The upper 100 m of the water column warmed 0.8 • C during the period 1950-1992, representing an average warming of 0.019 • C per year (Roemmich, 1992). Comparing the CalCOFI record with the temperature time series maintained at the Scripps pier, showed a slightly stronger warming trend inshore at the Scripps pier (Roemmich, 1992). From 1950 to 1993, the heat content integrated over the upper 200 m of the water column increased by 6.2-9.1% inshore and 2.4-7.1% offshore in the central and southern CCE (Palacios et al., 2004). From 1951 to 2001, heat content in the upper water column (20-200 m) increased with time roughly ten times more than the heat content of the deeper waters (300-500 m) (Hsieh et al., 2009).

Ocean Acidification
In the CCE, changes in carbonate chemistry can be driven both by changing water mass properties and circulation patterns (Meinvielle and Johnson, 2013), as well as the addition of anthropogenic CO 2 from the atmosphere into the water (Bednarsek et al., 2014). McClatchie et al. (2016) examined changes in carbonate system parameters using the empirical equations derived for the California Current region (Alin et al., 2012). Near the coast, more acidic conditions were found at shallower depths than at CalCOFI stations offshore; the aragonite saturation horizon ( = 1) depth ranged from 84 to 267 m offshore and 50 to 250 m inshore (McClatchie et al., 2016). From 1985 to 2011, an acidification trend was detected; the average aragonite saturation horizon ( = 1) depth shoaled by 18 m, and an 11% increase in acidity was observed at nearshore stations for depths between 15 and 100 m (McClatchie et al., 2016). This analysis likely underestimated acidification because it did not incorporate changes due to the rising atmospheric CO 2 concentration. For the California Current region, Bednarsek et al. (2014) estimated a potential anthropogenic contribution of 1.19 µmol kg −1 yr −1 to source-water DIC that is upwelled off California, based on an increasing trend in North Pacific surface water DIC.
While a continuous time-series of inorganic carbon system measurements is not available from CalCOFI to study ocean acidification (OA) trends, a multi-decadal record of measurements exists ( Figure 2C). From 1983 to 1994, DIC measurements from ∼10 m depth were collected by Charles David Keeling from a number of CalCOFI stations. Since 2008, DIC measurements have been collected at certain CalCOFI stations at depths between 0 and 200 m and analyzed by Andrew G. Dickson's lab at Scripps Institution of Oceanography, La Jolla, CA. Due to historical factors, Line 90 in the CalCOFI grid has the best coverage for DIC measurements. Alkalinity measurements of CalCOFI bottle samples have also been made during the periods 1984-1993 and 2007-present, and more recently pH measurements have been made from 2015-present. Additionally, underway pCO 2 measurements were conducted in collaboration with MBARI scientists from the mid 1990s to late 2000s on CalCOFI cruises, and more recently are being collected by Todd Martz's group at Scripps Institution of Oceanography, La Jolla, CA, since 2015. Since 2009, a pH sensor has also been added to the CTD package. Integrating these various carbonate system measurements across the CalCOFI time series to look at spatial and temporal OA trends is an area of current research.

Ocean Deoxygenation
Oxygen levels in the ocean are declining (Schmidtko et al., 2017;Breitburg et al., 2018), driven by changes in ocean ventilation, transport, and oxygen consumption linked to climate change (Levin, 2018;Oschlies et al., 2018). CalCOFI oxygen measurements are essential to understanding and resolving the oxygen dynamics in the SCCE. From 1984 to 2006, strong linear declines in oxygen were observed; declines were strongest nearshore, the largest relative declines occurred at 200-300 m depth, while the largest absolute declines occurred in the upper water column with oxygen losses of 1-2 µmol O 2 kg −1 year −1 (Bograd et al., 2008). At a nearshore CalCOFI station off San Diego, the depth of the hypoxic boundary (O 2 < 60 µmol kg −1 ) shoaled by up to 90 m during this 23 years period, and shoaled by 41 m on average across the CalCOFI grid (Bograd et al., 2008).
Compared to oxygen loss observed across the world's oceans (Schmidtko et al., 2017), the SCCE has one of the most rapid observed oxygen declines. However, the unique length of the CalCOFI time-series allowed these changes to be contextualized within longer term regional trends. McClatchie et al. (2010) analyzed the extended CalCOFI time-series  and found that the recent linear decrease in oxygen was not representative of the whole time-series. Instead, a clear multidecadal signal emerged, with the period between 1950 and 1957 characterized by low oxygen conditions, the period between 1980 and 1987 by higher oxygen conditions, and 2000-2007 by low oxygen conditions ( Figure 2D, McClatchie et al., 2010).
The non-monotonic nature of changing oxygen conditions in the CCE is related to natural low-frequency variability within the system (Deutsch et al., 2005). Observed declines in oxygen are thought to be associated with increased volume of advected Pacific Equatorial Water which is transported northward by the California Undercurrent (Meinvielle and Johnson, 2013;Bograd et al., 2015;Ren et al., 2018) and has a low oxygen signature. At shallower depths (<150 m), increases in respiration may also contribute to oxygen decline (Booth et al., 2014;Bograd et al., 2015;Ren et al., 2018). At shorter timescales, ENSO events affect the oxygen dynamics in the SCCE by altering the thermocline depth, with La Niña conditions characterized by lower oxygen and pH conditions (Nam et al., 2011). In general, any declines in oxygen are also coupled with decreases in pH; thus, deoxygenation and acidification will be simultaneous stressors for marine organisms in the CCE (Levin, 2018). Due to the high level of natural variability of oxygen in the CCE, forced declines in oxygen due to secular climate change are only predicted to emerge by ∼2030 (Long et al., 2016). The CalCOFI dataset is essential to ongoing monitoring of deoxygenation trends in the SCCE.

Ecosystem Responses to Climate Forcing
Marine communities are not static, but change through time in response to climate forcing. Changes can manifest through shifting species composition, changes in abundance, diversity, distribution, and trophic interactions (Poloczanska et al., 2016;Griffiths et al., 2017). Since CalCOFI collects both hydrographic and biological samples, it offers insight into how the SCCE responds to climate forcing at different timescales.
The CalCOFI ichthyoplankton time series has revealed that the fish community exhibits changes in distribution, abundance, and diversity in relation to climate forcing (Hsieh et al., 2005(Hsieh et al., , 2008(Hsieh et al., , 2009Koslow et al., , 2013Koslow et al., , 2015Koslow et al., , 2017Koslow et al., , 2018McClatchie et al., 2016McClatchie et al., , 2018Peabody et al., 2018;Siegelman-Charbit et al., 2018). Hsieh et al. (2008Hsieh et al. ( , 2009 found that from 1951 to 2002, oceanic and coastal fishes in the SCCE shifted their distribution poleward in response to warming, with dielvertical migrating species showing higher responsiveness to environmental change than non-migrating species (Hsieh et al., 2009). Abundance of oceanic species was significantly positively correlated with temperature, irrespective of if fish were warm water, cold water, or broadly distributed (Hsieh et al., 2009). Hsieh et al. (2009) posit that intensified stratification observed in the SCCE during the warm period created a more suitable habitat for oceanic species, thus leading to a spatial distribution shift shoreward during this period, and higher abundances within the CalCOFI grid. Similar trends were found when analyzing relationships across interannual and multidecadal scales, suggesting that these relationships may also manifest with secular climate change.
Trends in the CalCOFI ichthyoplankton time-series have also been compared to an independent time-series of ichthyoplankton collected from California power plant cooling water intakes in Southern California. Siegelman-Charbit et al. (2018) found that the abundances of coastal cold-water affinity fish taxa in the SCCE declined dramatically from 1970 to the 2000s, and this decline was related to ocean warming, increased coastal upwelling, reduced offshore upwelling and advection, as well as a decline in productivity and zooplankton prey availability. Koslow et al. (2015) also reported declining larval fish abundances across the same period, driven by declines in Pacific hake, Northern anchovy, and Pacific sardines. These results contrast with findings by Hsieh et al. (2009) for oceanic fish species, and indicate that oceanic and coastal species respond differently to climate forcing in the SCCE.
In general, warm periods in the CalCOFI record are associated with increases in fish species richness due to an influx of tropical or warm-water affinity species from the south (McClatchie et al., 2016(McClatchie et al., , 2018Koslow et al., 2017). This trend has been seen for both the 1997/1998 El Niño as well as during periods associated with increased intrusion of warm southern water (1985( -1996( and 1999( -2014( ) (McClatchie et al., 2016. Species richness increased despite the more acidic signature of southern water, suggesting that fish community species richness is not negatively impacted by current levels of acidification (McClatchie et al., 2016). Across the CalCOFI time-series (1969-2011) fish species richness is positively correlated with periods of warmer temperatures, higher midwater oxygen concentration, and the warm phase of the PDO . It is worth noting that conditions that give rise to higher ichthyoplankton diversity (i.e., species richness and community evenness) are associated with poorer primary productivity conditions, and thus correspond to lower fisheries productivity .
Analyses of the CalCOFI ichthyoplankton time-series have also informed our understanding of warm and coldassociated fish species (Hsieh et al., 2005;McClatchie et al., 2016McClatchie et al., , 2018. Mesopelagic cool-assemblage species include the Northern Lampfish (Stenobrachius leucopsarus) and California Smoothtongue (Leuroglossus stilbius), while mesopelagic warmassemblage species include the Mexican Lampfish (Triphoturus mexicanus) and the Panama Lightfish (Vinciguerria lucetia) (McClatchie et al., 2016). Warm assemblage species tend to increase in abundance during warm events (McClatchie et al., 2016;Koslow et al., 2018) and these species have shown a secular trend of increasing relative abundance from 1951 to 2016(McClatchie et al., 2018. From 2011 to 2015, mesopelagic fish with tropical or subtropical distributions became increasingly dominant in the SCCE . Commercially important species also show relationships with climate forcing. Bocaccio Rockfish (Sebastes paucispinis) ichthyoplankton abundance is negatively correlated with CalCOFI SST, while Pacific Sardine and Pacific Mackerel (Scomber japonicus) ichthyoplankton are positively correlated with CalCOFI SST (Hsieh et al., 2005). To examine changes in fishing pressure, an annual review of California fisheries (CDFW, 2015) can be examined.
Oxygen dynamics in the SCCE also correlate with trends in the CalCOFI ichthyoplankton time-series. Mesopelagic fish abundance showed a 63% decline between periods of high and low oxygen in the CalCOFI record, suggesting that deoxygenation may result in declining mesopelagic fish abundances in the future . Lindegren et al. (2017) also found that mesopelagic fish abundance was positively correlated with oxygen concentration, and exhibited a threshold relationship. In contrast, oxygen was not found to be a significant variable in explaining trends in the abundance of zooplankton, euphausiids, pelagic fish, or Pacific hake in the SCCE (Lindegren et al., 2017). The vertical distributions of adult fish may also be affected by oxygen changes but this cannot be seen with the CalCOFI dataset, which only samples the early life history stages.
Key findings have also been made with the CalCOFI zooplankton dataset. Between 1951 and 1993, zooplankton displacement volume decreased significantly, correlated with a warming of surface waters (Roemmich and McGowan, 1995). From 1951 to 2005, the mean zooplankton displacement volume decreased by 64% off southern California and 68% off central California, which was attributed to a decline in high volumelow carbon zooplankton (i.e., "cool phase" salps) in the latter half of the time-series (Lavaniegos and Ohman, 2007). Zooplankton biomass in terms of carbon showed no long-term trend, highlighting that the zooplankton displacement volume is not directly correlated with zooplankton biomass in terms of carbon (Lavaniegos and Ohman, 2007). Planktonic copepods, followed by euphausiids, dominate the carbon biomass in the zooplankton samples (Lavaniegos and Ohman, 2007). The zooplankton data have also shown that certain krill species (euphausiids) have warm (e.g., Nyctiphanes simplex and Euphausia eximia) and cold (e.g., Euphausia pacifica and Thysanoessa spinifera) affinities, and that krill distributions respond to natural climate forcing through ENSO and PDO cycles (Brinton and Townsend, 2003).
Since CalCOFI samples capture species that occupy multiple trophic levels as adults, the data have been used to study how climate forcing affects food web interactions and ecosystem resilience (Lindegren et al., 2016(Lindegren et al., , 2017. Despite the presence of pronounced multidecadal variability and associated changes in the abundance of certain species and groups in the SCCE, ecosystem functioning has largely remained unchanged, due to complementarity of species within trophic levels (Lindegren et al., 2016). In other words, functionally similar species have opposite responses to climate drivers, resulting in overall ecosystem stability. Bottom-up control is the primary mode of ecosystem regulation in the SCCE (Lindegren et al., 2017), meaning that changes to primary production are carried up the food chain and directly affect top predator abundances. However, under El Niño-like conditions (i.e., weak upwelling, low nutrient concentrations, and low primary production), Lindegren et al. (2017) found that bottom-up forcing is reduced in the SCCE and the system exhibits "wasp-waist" dynamics with mid-trophic levels exerting moderate top-down forcing (Cury et al., 2000). Nitrate concentration, wind stress curl-driven upwelling, and temperature were found to be the best single predictors of ecosystem state across multiple trophic levels (Lindegren et al., 2016), suggesting that climate impacts on these variables will have the greatest impact on the overall ecosystem state of the SCCE.
Regime shifts, meaning rapid, abrupt and persistent changes in ecosystem structure and function (Mollmann et al., 2015), have also been detected in the SCCE using the CalCOFI dataset. Based on the ichthyoplankton time-series, regime changes were detected in 1965 and 1976 for both the summer and spring larval fish assemblage, and additionally in 1983 and 1990 for the spring assemblage (Peabody et al., 2018). The dominant pattern since the strong 1976 regime shift has been a long-term trend toward a more southern and offshore assemblage composition in the CalCOFI region (McClatchie et al., 2018;Peabody et al., 2018). When looking at trends across multiple trophic levels in the CalCOFI dataset, significant break points were identified in 1965, 1976, 1987, and 1999(Lindegren et al., 2016. Major ecosystem regime shifts in the CCE often occur when the PDO and NPGO show strong, simultaneous, and opposite sign reversals . Therefore, climate change impacts in the SCCE will depend on the phase relationship between fluctuations of the PDO and NPGO modes  along with the forcing of secular climate warming.

Future Climate Change Projections for the California Current
While CalCOFI has provided many insights into how the SCCE has responded to natural and anthropogenic climate forcing over the last 70 years, the twenty-first century will likely be defined by more rapid change and a larger contribution of anthropogenic forcing. Several valuable modeling products are available that allow for projecting future conditions in the California Current under different climate scenarios. However, developing accurate models at timescales that matter to resource managers is still an active area of research. Some of the key challenges involve: properly reproducing biogeochemical dynamics, which are essential for projecting changes in nutrients, oxygen, and pH; improving model resolution, which is important for resolving processes such as upwelling at smaller spatial scales; and accurately reproducing coastal dynamics, which are key for nearshore species such as Dungeness Crab (Metacarcinus magister) and California Spiny Lobster (Panulirus interruptus) (Holt et al., 2009(Holt et al., , 2014Stock et al., 2011). For managing some living marine resources, climate projections on the 3-10 years time scale are the most useful, however, these timescales are not well represented within global climate and ocean models (Busch et al., 2016).
From 1950 to 2019, the Large Ensemble Community Project (LENS) model (Kay et al., 2015) predicted an average increase in SST by 0.5 • C, a decrease in surface pH by 0.1, and a decrease in surface oxygen concentration of 2.56 mmol m −3 for the SCCE. By 2050, an additional 1.25 • C warming, a 0.1 surface pH decrease, and a 5.75 mmol m −3 decrease in surface oxygen concentration is predicted under the RCP 8.5 "business as usual" emissions scenario (Figure 3). By end of century (2100), conditions in the CalCOFI region are projected to be 3.2 • C warmer (SST), surface pH will be 0.3 units lower, and surface oxygen concentrations will have decreased by 13.83 mmol m −3 compared to 2019 conditions under the RCP 8.5 emissions scenario (Figure 3).
The CalCOFI time series shows a stronger warming trend than that projected by the LENS ensemble from 1950 to 2019. While the LENS ensemble showed an average increase in SST by 0.5 • C, CalCOFI observations show an average SST warming of 1.45 • C for the core CalCOFI grid. For 1950For -1952, the average SST across all CalCOFI core stations was 15.16 • C (ranging from 10.11 to 21.14 • C), while in 2016-2018, average SST in the same region was 16.61 • C (ranging from 11.20 to 23.23 • C). Thus, observations show a warming trend between ∼1 and 2 • C for the mean, minimum and maximum SST values observed across the time series. Thus, the LENS ensemble predictions for warming appear to be conservative compared to CalCOFI observations. Gruber et al. (2012) and Hauri et al. (2013) used a regional ocean modeling system (ROMS) model at eddy-resolving resolution for the CCE to look at how OA may progress from 1995 to 2050. Under the high emission scenario, which represents 541 ppm P CO2 by 2050, the mean annual pH in the CCE decreased by 0.2 pH units by 2050, compared to pre-industrial conditions, and the aragonite saturation horizon shoaled by ∼250-300 m (Gruber et al., 2012). During the summer upwelling season, low pH conditions were prominent, and by 2050 under the high emission scenario, large stretches of the nearshore central CCE became undersaturated with respect to aragonite (Gruber et al., 2012). These findings suggest that OA will progress rapidly in the CCE, and that nearshore areas (<10 km) are more susceptible to low pH and undersaturated conditions (Gruber et al., 2012). Under current conditions in the CCE, there is a strong positive correlation between the proportion of pteropods with severe shell dissolution and the percentage of water in the upper 100 m that is undersaturated with respect to aragonite (Bednarsek et al., 2014). Therefore, habitat suitability for animals that are sensitive to undersaturated conditions is predicted to decline in the CCE under climate change scenarios (Bednarsek et al., 2014).
Based on observed trends of rapid oxygen loss in the SCCE (Bograd et al., 2008, predicting future oxygen changes is an area of management interest. A recent study found an intriguing statistical relationship between subsurface anomalies in the core of the North Pacific Current and decadal changes in oxygen in the CCE (Pozo Buil and Di Lorenzo, 2017). Based on ocean observations and an ocean reanalysis product for the period 1950-2010, they found that anomalies in the North Pacific can be used to accurately predict oxygen trends on the US West Coast over a time-scale of 9-12 years and predicted a new period of strong oxygen decline by 2020 (Pozo Buil and Di Lorenzo, 2017). If accurate, predictions at these time-scales could be helpful for guiding management decisions if there exist known relationships between the abundance of fisheries species and environmental oxygen conditions. Trends in salinity showed high correlation with trends in oxygen, suggesting that salinity can be used as a passive tracer of circulation-driven changes in the oxygen dynamics (Pozo Buil and Di Lorenzo, 2017), in contrast with changes due to warming or biological use, which may become more pronounced with climate change in the CCE by year ∼2030 (Long et al., 2016).
Reduction in primary production is an additional climate change stressor predicted to impact marine ecosystems (Bopp et al., 2013). However, there are large regional differences and high levels of uncertainty relating to projected changes in primary production across Earth system models (Steinacher et al., 2010). In eastern boundary upwelling systems (EBUS) such as the CCE, changes in upwelling-favorable winds may drive geographic and temporal changes in primary productivity. A meta-analysis found that upwelling-producing winds have strengthened over the past ∼60 years in three EBUS including the California Current (Sydeman et al., 2014). Coastal upwelling-favorable winds have intensified and will continue to do so in poleward portions of EBUS . A climate model ensemble further predicts seasonal variability in CCE upwelling trends with intensification occurring in the spring but weakening in the summer (Brady et al., 2017). Changes in nutrient supply could also modulate productivity: Rykaczewski and Dunne (2010) explored this question using a biogeochemistrycoupled earth system model, showing increases in nitrate supply and productivity for the CCE under climate change, despite increases in stratification and potential dampening of vertical transport. While increased primary production has positive, bottom-up effects for fisheries production, these changes could be accompanied by decreases in oxygen and pH (Rykaczewski and Dunne, 2010).
If upwelling patterns follow current trends and predictions, there may be spatially heterogeneous variability in patterns of ocean warming and primary production in the future, though the CCE is particularly uncertain (Wang et al., 2015). Increased warming and stratification offshore would result in decreased primary production, while the influx of cold, nutrient-rich water due to increased coastal upwelling could mitigate warming and actually increase primary production near the coast (Garcia-Reyes et al., 2015). A pattern of warm offshore waters coupled with cool, productive coastal waters was observed in the central CCE in 2015 (Santora et al., 2017). This unusual oceanographic state resulted in a unique assemblage where subtropical mesopelagic taxa as well as taxa such as rockfishes that thrive under cooler, productive conditions were both very abundant in the same region, and overall species richness was the highest recorded since sampling began in 1990 (Santora et al., 2017).

Fisheries Management on the US West Coast and CalCOFI
To effectively explore how monitoring data can best be incorporated into management decisions from the perspective of climate change adaptation, a basic understanding of the structure of the policy and management framework is necessary. For the US West Coast, several jurisdictional management bodies and policies guide marine resource management. This section briefly discusses how the data collected by CalCOFI fits within the management priorities for the US West Coast. CalCOFI already contributes directly to US West Coast fisheries population assessments and management (McClatchie, 2014) by providing the PFMC with information about fish population abundances, oceanographic conditions, and ecosystem status. CalCOFI larval abundance time-series are directly used as indices of spawning stock biomass in the Bocaccio (He and Field, 2017) and Cowcod (Sebastes levis; Dick and MacCall, 2014) rockfish stock assessments, and were recently used to monitor population trends of Blue Rockfish (Sebastes mystinus; Dick et al., 2017). CalCOFI data are also used to quantify Daily Egg Production values that contribute to Pacific Sardine (Hill et al., 2017) and Northern Anchovy (Fissel et al., 2011;Dorval et al., 2018) stock assessments. Additionally, CalCOFI larval data track the status of harvested species such as Market Squid (Koslow and Allen, 2011) and Jack Mackerel , and taxa such as mesopelagic fishes that are considered by the PFMC to be important forage species for higher trophic levels (McClatchie et al., 2018). This information is used by the PFMC to set harvest limits at levels that support sustainable fisheries.
In addition to evaluating the status of individual species, fisheries management in the United States also emphasizes holistic monitoring of ecosystem dynamics (Levin et al., 2008), and NOAA Fisheries is moving toward a next generation stock assessment approach which incorporates ecosystem links . The inclusion of ecosystem factors such as oceanographic variables into single-species stock assessments, termed an Ecosystem Approach to Fisheries Management (EAFM), can allow for better stock-focused management decisions (Link, 2016). Currently, most stock assessments do not take environmental conditions into account due to the difficulty of elucidating predictable recruitment-environment relationships that explain enough recruitment variability to be useful for management (Myers, 1998;Jacobsen and McClatchie, 2013;Lindegren et al., 2013;Tolimieri et al., 2018). A global review of 1200 marine fish stocks revealed that only 24 stocks were managed in a way that incorporated ecosystem drivers such as climate (Skern-Mauritzen et al., 2016). The Pacific Sardine is one example in which oceanographic data directly informs an aspect of single-species stock assessment. Pacific Sardine stock productivity is higher under warmer conditions, so Pacific Sardine management uses SST averaged across the core CalCOFI grid to calculate the over fishing limit and the acceptable biological catch (Hill et al., 2017), with fishing effort increasing under higher temperatures. However, recent efforts to re-evaluate the temperature dependence of Pacific Sardine stock productivity concluded that the earlier positive correlation with SST may no longer be valid (Zwolinski and Demer, 2019).
Beyond a single-species focus, Fisheries Ecosystem Plans (FEPs) have been developed to guide fisheries management toward implementing ecological principles (Marshall et al., 2018). FEPs are encouraged but not mandated by federal fisheries policy, and five of eight regional fisheries management councils, including the PFMC (Pacific Fishery Management Council [PFMC], 2013), have developed FEPs. The PFMC also currently includes an Ecosystem Consideration appendix to their groundfish, coastal pelagic species, and salmon FMPs as an effort to inject monitoring data more directly into fisheries management decisions, which includes some CalCOFI data.
The CalCOFI dataset is well suited to not only inform the status of individual species but also to provide insight into the ecosystem state of the SCCE and how it responds to both short and long-term environmental forcing (McClatchie, 2014). Integrated ecosystem assessments (IEAs) have been developed by NOAA to support Ecosystem Based Management (EBM) (Levin et al., 2009), and the California Current was selected as the first region for the development of a full IEA, which began in 2009 and is called the California Current Integrated Ecosystem Assessment (CCIEA). CalCOFI contributes several of the indicators used by the CCIEA (Figure 4), including the abundances of larval Northern Anchovies, Market Squid, Pacific Hake, Jack Mackerel, Sanddabs, Pacific Sardine, Shortbelly Rockfish (Sebastes jordani), and a category of warmwater and coolwater mesopelagic fish larvae. Since 2013, the IEA team has presented an annual report to the PFMC on the state of the CCE including indicator trends using CalCOFI data (Harvey et al., 2018). Organizations that use the CCIEA to inform decision-making include the PFMC, the West Coast National Marine Sanctuaries, the West Coast Governors Alliance, and the North Pacific Marine Science Organization (Samhouri et al., 2013).

How Can CalCOFI Support Stakeholder Needs for Addressing Climate Change Adaptation?
While the CalCOFI dataset already supports efforts at the state and federal level to consider management approaches that are adaptive to climate change, progress can be made in linking monitoring data with actionable strategies for climate change adaptation. One approach to this is through recognizing and meeting the needs of stakeholders concerned with climate change adaptation. Here we consider several different stakeholder user groups of CalCOFI data products (Figure 5).

Academic Community
One of the most obvious stakeholders is the science community, which includes stock assessors, the climate modeling community, as well as scientists studying the ecological impacts of climate change. For these communities, time series and sampling consistency, attention to accuracy of measurements, rapid data availability following cruises, and easy data access are key. Given that the CalCOFI time series currently extends almost 70 years, it allows analysis of climate related interactions at multidecadal, interannual, and seasonal time scales in addition to secular (longterm) climate change. Consequently, it supports two key research needs identified by the 2009 National Climate Assessment report: improving understanding of the climate system and drivers, and improving understanding of climate impacts and vulnerability. Survey data are essential for validating regional ocean models, Global Circulation Models, and Global Climate Models and CalCOFI data are integrated into regional modeling efforts in terms of model calibration, data assimilation, hindcasting and reanalysis model products. Data-assimilative ROMS which incorporate CalCOFI data have been developed and are used by research groups at the University of California, Santa Cruz and the University of California, Los Angeles.

Fisheries and Aquaculture Community
The management community is also a key stakeholder of CalCOFI data products. Improving the paths by which monitoring data are used by management bodies will improve the ability for managers to make decisions that are more responsive to climate change. Answering how climate-related effects can be incorporated into adaptive marine resource management processes is a key question for NMFS (Busch et al., 2016), and the PFMC is currently considering this as part of its Climate and Fishing Communities Initiative.
Similarly, state management agencies such as CDFW are stakeholders of CalCOFI data products. CalCOFI provides fishery-independent time series that can be incorporated into management for key state-managed fisheries, such as California Spiny Lobster and Market Squid (Koslow and Allen, 2011;Koslow et al., 2012), and the nearshore sampling stations provide a look at more nearshore changes that are relevant to the state. There are no fishery-independent time series for many recreational fisheries, and CalCOFI data are well suited to inform stock status for taxa such as Saltwater Basses (Paralabrax spp.) that may be overfished (Jarvis et al., 2014). Early life history stages of these forage species are captured in the CalCOFI dataset, providing an existing time series index of spawning stock biomass for reference.
Climate change can also lead to fish stocks shifting their distribution, which can have important management implications regarding allocation. When a stock is shared across national boundaries, this can add additional complexity. The southern end of the CalCOFI domain samples along the US-Mexican border meaning that it can provide information on shifting distributions across international boundaries. Data from the Mexican exclusive economic zone are more limited, although data from IMECOCAL (Baumgartner et al., 2008) and some catch data from Mexico are incorporated into the Pacific Sardine stock assessment. In contrast to shared US-Canadian stocks, no explicit fisheries stock sharing agreements exist between the US and Mexico. One potential emerging fishery of interest is the Humboldt Squid, which is fished in Mexican waters but is only periodically caught off California. Humboldt Squid appears to expand its range northward during warm, low oxygen periods (Gilly et al., 2013) and could emerge as a valuable fishery in a warmer future. Since commercial fishermen need specific permits to allow for fishing new species, early indicators of a shift in distribution are key to supporting responsive management actions.
The fisheries industry is also a stakeholder of CalCOFI data. The CalCOFI region captures an area of great interest to the recreational fisheries industry, which is an engine of the economy in Southern California. Many recreational fisheries species are not federally managed and therefore population trends of these species are not as well known. Since CalCOFI samples the early life history stages of over 400 fish species, it provides a valuable record of fisheries species trends that receive less management attention, but may be impacted by climate change. Sport fishing organizations and recreational fishing party boats can all be thought of as potential stakeholders of CalCOFI data. CalCOFI is also important for the commercial fishing industry. The California Wetfish Producers Association pays close attention to CalCOFI data since most wetfish (i.e., sardines, anchovies, squid, and mackerels) in the state come from southern and central California. This group is actively involved with the PFMC and coastal pelagic species management, and industry-funded researchers for the California Wetfish Producers Association use CalCOFI paralarvae counts of Market Squid in their monitoring studies.
Beyond impacts on fisheries species, there is also a need to consider adaptive capacity of the fisheries industry, which depends on the types of fisheries that are present. For example, while stock assessments focus on single-species, fishermen can fish multiple species and switch effort based on both biological and economic factors. Fuller (2018) shows that fisheries in the CCE are strongly connected by human participation, FIGURE 5 | Linking basic monitoring to decision-making for climate change adaptation requires integrating different data streams with the needs of stakeholders. CalCOFI exists alongside a wealth of additional sampling programs that, in concert, allow for comparing changes seen in the CalCOFI record across regions (upper, left), monitoring frameworks within a region (upper, middle), and sampling programs on different fish life stages (upper, right). Trends observed from these monitoring programs can support different stakeholder groups (bottom) including the scientific, management, fishing, and conservation communities on climate change adaptation strategies and actions. Arrows show that information can flow both ways.
with participants in the Dungeness Crab pot fishery showing high connectedness to many different fisheries. Incorporating considerations about effort-switching may be critical in thinking about climate change adaptation for fisheries. Since CalCOFI samples multiple species, analyses could be catered to reflect fishing practices by grouping species by connectivity of fisheries efforts, as opposed to focusing on single-species analyses or ecological groupings.
While there is currently little offshore aquaculture along the southern U.S. West Coast, there are ongoing efforts to increase its development. Brander (2010) points out that aquaculture resembles terrestrial animal husbandry more than it does wild capture fisheries, and therefore has greater scope for decisionmaking that is adaptive to climate change. In the future, the aquaculture industry could become an additional stakeholder of CalCOFI data. The Catalina Sea Ranch is a mussel aquaculture farm, which operates in federal waters offshore of Los Angeles, CA. It is the first aquaculture facility permitted for commercial use in federal waters in the US, and the company plans to expand its 100 acre underwater shellfish farm to 1,000 acres. Shellfish aquaculture may be especially concerned about OA impacts. There have also been proposals for finfish aquaculture in federal waters off California, such as a recent proposal by Pacific Ocean Aquafarms (previously the Rose Canyon Fisheries Farm) to raise California Yellowtail (Seriola lalandi) and White Seabass in offshore net pens off San Diego. While it is unclear if this venture will be approved, it is worth noting that just south of the border in Mexico, tuna aquaculture is a widespread practice, with several companies, including Baja Aqua Farms, raising tuna in offshore pens. While aquaculture groups will likely invest in their own ocean monitoring platforms, CalCOFI can provide information on larger-scale trends (i.e., changes in pH or oxygen) that are occurring independent of aquaculture practices.

Ecosystem Management Community
Planners and managers of marine conservation areas can also be considered stakeholders of CalCOFI data. Rockfish conservation areas, marine protected areas (MPAs), and national marine sanctuaries all occur within the CalCOFI region (Figure 6). While these conservation areas have different mandates and roles, CalCOFI can provide regional context for interpreting local trends that may be related to specific management or marine use actions. Nine CalCOFI stations occur within National Marine Sanctuaries (NMS) (seven within the Monterey Bay NMS, one within the Greater Farallones NMS, and one within the Channel Islands NMS) ( Figure 6A). The Channel Islands NMS makes use of CalCOFI data in its formal status reports but there are clear opportunities for CalCOFI data to further inform the Monterey Bay and Greater Farallones NMS. CalCOFI stations also occur within three State Marine Reserves (South Point, Begg Rock, and Point Conception) and two State Marine Conservation Areas (Point Dume and Dana Point) ( Figure 6B). Six CalCOFI stations also fall within the Cowcod Conservation Area (CCA) (Figure 6C). Data from CalCOFI stations within and outside of the CCA have been used to show trends in recovery of overfished rockfish stocks (Thompson et al., 2017a). Understanding vulnerability of marine conservation areas to climate change and considering adaptation measures are active areas of interest for planners and managers of these areas; the CalCOFI dataset can provide important regional insight for climate change studies and opportunities to evaluate management action.
For the northern and central CCS, Sydeman et al. (2014) developed a Multivariate Ocean Climate Indicator (MOCI) to aid in the evaluation of ecosystem management decisions. The MOCI synthesizes 14 basin-and regional-scale atmospheric and oceanographic indicators from a variety of time series (Sydeman et al., 2014). If this multivariate index proves useful in guiding management decisions, a similar index could be developed for Southern California by incorporating CalCOFI data. Currently, CalCOFI scientists produce an annual State of the California Current report, which provides a review of environmental conditions in the CCE (e.g., Thompson et al., 2018). Species that are responsive to climate forcing and are dominant in abundance in the community can be good bioindicator species and several taxa from the CalCOFI time series are used as bio-indicators for the CCIEA.
Certain non-governmental organizations such as EcoAdapt and The Nature Conservancy Climate Fisheries group are also actively involved in developing climate change vulnerability assessments and strategies for climate change adaptation for fisheries and marine ecosystem management. These groups may use CalCOFI data or products such as the State of the California Current report. Additionally, the CalCOFI Program now exists alongside a wealth of additional sampling programs that have developed in the region through time. These include: NOAA's Northwest Fisheries Science Center Groundfish Trawl Survey and the Southwest Fisheries Science Center Pelagic Juvenile Rockfish Recruitment and Ecosystem Assessment Survey and summer Coastal Pelagic Survey, the 40 + years ichthyoplankton time-series from the Southern California coastal power plant intakes (Miller and McGowan, 2013;Siegelman-Charbit et al., 2018), the CCE LTER Program, commercial and recreational fishing data (Perry et al., 2010), the IMECOCAL oceanographic survey south of the border (Baumgartner et al., 2008), the southern (SCCOOS) and central California (CeNCOOS) coastal ocean observing programs, as well as others. Leveraging knowledge, collaborations, and investment across these different groups will allow for better scientific guidance for climate change adaptation for fisheries and ecosystem management.

Climate Change Impacts on Fisheries
CalCOFI was originally developed to inform Pacific Sardine management, and continues to inform sustainable fisheries management on the US West Coast. As climate change progresses, CalCOFI can be used to monitor climate change impacts on fisheries species. The effects of climate change on fisheries production is a concern for fisheries management around the world (Barange et al., 2018;Cheung, 2018).
Climate change can affect fisheries species through direct effects of hydrographic changes (e.g., temperature, oxygen, pH) on physiology and the indirect effects of ecological changes mediated through the food web and multi-species interactions (Brander, 2010). Overall, trends indicate increased catch potential in high latitudes and decreased catch potential in lower latitudes (i.e., tropics) with climate change (Cheung et al., 2010). On the US West Coast, fisheries for Albacore Tuna (Thunnus alalunga) may increase with climate change and new fisheries, such as on the Humboldt Squid (Dosidicus gigas), may develop (Chavez et al., 2017).
Managing fisheries effectively under climate change will require the ability to identify and track trends, develop indices of change, and develop early indicators of potential shifts in biological communities, which are all areas that CalCOFI directly contributes to. Management bodies need synthesized data products, including relevant indices, to make decisions, as well as reliable climate forecasting abilities over 1-10 years time frames that guide decision-making. To do this, linkages between environmental indices and fisheries need to be effectively identified so that forecasting can be used to inform management decisions.
NOAA Fisheries has developed a Climate Science Strategy (CSS) for including climate-related information in managing fisheries and protected species, which includes conducting climate vulnerability analyses, establishing ecosystem indicators, and developing capacity for conducting management strategy evaluations (MSEs) of climate-related impacts on management targets, priorities, and goals (Busch et al., 2016). MSEs are a modeling tool which allow scientists and managers to compare the performance of different management strategies in meeting specific management objectives (i.e., fishery resource status) (Punt et al., 2016) and can be used to manage climate change risks (Plaganyi et al., 2013). To implement the CSS, regional NMFS offices develop Regional Action Plans, which are developed jointly by scientists, managers, and stakeholders. For the US West Coast, a Western Regional Action Plan has been developed (Noaa Nw/Sw Fisheries Science Centers, 2016) that lays out the NOAA Fisheries climate change strategy and the Western Regional Implementation Plan provides the means to execute the strategy. Incorporating environmental data into stock assessments to improve fisheries management under climate change is also a recognized priority for NMFS (Karp et al., 2018;Lynch et al., 2018).
While the work on climate-ready fisheries management is relatively new, several adaptation approaches have been proposed including: addressing cumulative impacts to reduce stressors for fisheries species, preparing for emerging fisheries, accounting for climate effects in stock assessments, adjusting spatial boundaries for management areas, coordinating across static boundaries such as national borders, developing regional climate change projections to guide management decisions, expanding opportunities for real-time responses to climate, and promoting social-ecological resilience (Pinsky and Mantua, 2014). Vulnerability assessments, modeling tools and MSEs, and the use of adaptation frameworks that incorporate socio-ecological considerations are all tools that can be used to address fisheries adaptation to climate change (Lindegren and Brander, 2018).
Based on CalCOFI data, in the SCCE, secular climate change will likely lead to increased fish community diversity (McClatchie et al., 2016(McClatchie et al., , 2018Koslow et al., 2017), as warmaffinity species expand their ranges, but may be associated with lower fisheries productivity for current community dominants like Northern Anchovy and Pacific Hake . Changes in community composition are also expected. If conditions continue to warm, a change in dominance is expected from California Smoothtongue to Mexican Lampfish inshore and Northern Lampfish to Panama Lightfish offshore (McClatchie et al., 2016). A comparison of the 15 most abundant ichthyoplankton taxa by sampling region between the CalCOFI region and the region to the south off Baja California from the IMECOCAL Program, suggest that warming conditions may lead to increased relative abundances of Sanddabs (Citharichthys spp.), Pacific Sardine, and Pacific Mackerel, and decreased abundances of Pacific Hake, Rockfishes, and Croaker (Scianidae spp.) . Additionally, decreases in primary productivity could give rise to greater influence of mid-trophic levels on ecosystem dynamics in the SCCE (Lindegren et al., 2017).
States are also increasingly interested in understanding how climate change may affect their fisheries. In California, a report on readying California's fisheries for climate change was recently developed for consideration by the CDFW to help inform the state's process to amend the Marine Life Management Act (MLMA) Master Plan (Chavez et al., 2017). CDFW has also been developing vulnerability assessments of its primary fisheries species, and the California Ocean Science Trust, Ocean Protection Council, and the Ocean and Coastal Climate Advisory Team are actively engaged in considering climateready management strategies in California. A challenge of considering climate change in fisheries management is that across "short" time-scales relevant to decision makers (<10 years), the anthropogenic component of climate change only adds a small increment of change compared to the natural variability experienced by fisheries species in the region (Brander, 2010). However, extreme climate events can lead to rapid effects on a fisheries species, so changes to the frequency or timing of extreme events is also a climate change vulnerability (Brander, 2010).
The CalCOFI dataset captures the early life history stages of many of these species, allowing for monitoring of larval relative abundance trends through time as climate change progresses. Analysis of the CalCOFI time series suggests that commercially exploited fish species show clearer distributions shifts in response to climate forcing than unexploited species (Hsieh et al., 2008), and thus may be more sensitive to climate change. The CalCOFI record has also shown that the abundance of several groundfish, such as English Sole (Parophrys vetulus), Cabezon (Scorpaenichthys marmoratus), Aurora Rockfish (Sebastes aurora), Bocaccio Rockfish, Shortbelly Rockfish, and other rockfishes is especially low after El Niño events (1958,1983,1997), likely due to the low reproductive output of the adult population for those years (Lindegren et al., 2017). Following El Niño events, both dominance and assemblage spatial structure of the ichthyoplankton community tended to return to the long-term climatological average, suggesting resilience to temporary perturbations (McClatchie et al., 2016(McClatchie et al., , 2018. Given the considerable changes in warming, pH, and oxygen projected for the SCCE for the twenty-first century (Figure 3), as well as the uncertainty surrounding changes in primary production, the role of CalCOFI in monitoring climate change impacts in the SCCE will become even more critical.
Both "bottom up" and "top down" approaches exist to climate adaptation. One example of a "bottom up" approach is Australia's research program on Marine Climate Impacts and Adaptation (Hobday and Cvitanovic, 2017) that provides local fishermen and enterprises with knowledge and climate forecast products such as seasonal forecasts at small spatial scales. This directly provides fishermen with tools to guide their decisions. In contrast, the European Maritime and Fisheries Fund uses a "top down" approach that emphasizes strategic planning at large spatial and organizational scales (Lindegren and Brander, 2018). The US strategy for climate adaptation is more like the European approach, however, there may be opportunities to incorporate more "bottom up" elements that may increase participation from local stakeholders. Wilson et al. (2018) proposes that developing climate-ready strategies for fisheries will require increasing the level of participation in the management process by bringing together fishermen, non-governmental organizations, and traditional scientific and management bodies within an adaptive co-management framework. Using a "bottom up" approach to climate adaptation by hosting workshops or roundtables could increase participation, knowledge, and investment from the local community, and connect CalCOFI more directly with the fisheries industry stakeholders.

Limitations and Areas for Improvement
One limitation of using past trends to predict future changes is that ENSO, PDO, and NPGO climate forcing may not be representative of changes predicted with secular climate change. Koslow et al. (2018) show an intriguing recent breakdown in the observed relationship in oceanic forcing on the response of mesopelagic fish across the CalCOFI time series. While mesopelagic fish abundance was positively correlated with midwater oxygen concentration for much of the CalCOFI record, this relationship appears to be breaking down in the twentyfirst century as more hypoxia-tolerant warm-water affinity mesopelagic fishes such as Panama Lightfish and Mexican Lampfish move into the region and increase in abundance, despite declining oxygen levels . It is important to keep in mind that the ranges of most species extend beyond the CalCOFI grid, and therefore species abundance trends within the core CalCOFI grid do not necessarily equate to trends in overall population size (Weber and McClatchie, 2012).
The concept that environment ∼ abundance trends based on correlative data over the past 7 decades may not hold under climate change was also evident in 2018. Northern Anchovy have historically been categorized as a species that does well during relatively cold conditions while Pacific Sardine have been thought to thrive under warm conditions (Chavez et al., 2003). Biomass of adult and larval Northern Anchovy, however, approached record highs throughout the CCE from 2017 to 2018 . These spawning adults had to have been born and successfully recruited during the warmest sustained waters on record between 2014 and 2016 , thus challenging the notion that Northern Anchovy thrive under cold conditions. By contrast, Pacific Sardine remained at low abundances over the past decade even though ocean conditions superficially supported resurgence. These findings underscore the difficulty in predicting how species will respond to climate change and emphasize the need to understand the precise mechanisms affecting fisheries population dynamics . Even with 7 decades of observations, the patterns observed provide correlations in oceanographic processes that support (but do not necessarily prove) mechanistic drivers, such as external climate forcing. For this reason the CCE LTER is built around the goal of demonstrating the mechanisms driving patterns observed in CalCOFI's observations via experimental oceanography (Ohman and Hobbie, 2008).
In this study, we identify two additional areas where CalCOFI can better support stakeholder needs. The first is developing a more robust framework for OA monitoring within the CalCOFI Program, which aligns with state needs vocalized by the Ocean Protection Council. While many carbonate system measurements have been collected on CalCOFI cruises over the years, there have been no comprehensive efforts to integrate these datasets and couple them with the biological samples. However, these samples can provide a rich baseline for studying OA trends which are predicted to manifest rapidly in the CCS (Gruber et al., 2012;Hauri et al., 2013). Knowledge of baseline carbonate system conditions may also be important in informing the development of offshore shellfish aquaculture in the region and considering future OA vulnerability for this industry. Therefore, supporting continuous OA monitoring along CalCOFI lines should be a high priority.
Lastly, one new area which CalCOFI can contribute to through the NOAA CalCOFI Genomic Project is in monitoring changes in the microbial community. There is a concern that marine infectious diseases will increase under climate change (Burge et al., 2014). Monitoring pathogen abundances and understanding conditions that give rise to outbreaks can help manage this vulnerability. HABs may also increase in frequency and severity under climate change (Wells et al., 2015) with economic implications for the fisheries, aquaculture, and coastal sectors. Since novel high-throughput sequencing approaches can provide a detailed analysis of the composition as well as activity (i.e., gene expression) of the microbial community, these novel methodologies could help inform these emerging climate change adaptation issues. Additionally, these data may be useful for environmental consulting companies that monitor sewage and powerplant outfalls, by providing a regional baseline for their samples. Systematically incorporating this sampling into the CalCOFI cruises would further the ecosystem perspective gained from the CalCOFI Program. Frusher et al. (2014) outlined four preconditioning factors that are necessary for developing an effective climate adaptation process for fisheries management. These include: (1) early observations of rapid ocean change, coinciding with (2) observations of biological changes that provide focus for action, (3) within a region that has a strong history of marine focus and management, and (4) the presence of well-developed networks involving multiple marine resource users. Based on these, the US West Coast is well positioned to be able to respond to climate change effectively with an adaptive management framework for fisheries, aquaculture, and marine ecosystem use. The CalCOFI Program contributes in an essential way to all four of these factors, by collecting paired oceanographic and biological data on the marine ecosystem, fostering the long record and history of fisheries oceanography and studies on the effects of climate variability on marine resources for the region, and the presence of many partnerships (other scientific surveys, academic contributions, governmental scientists, and NGOs). For CalCOFI information to be most useful, managers need a fuller understanding of the linkages between biological and oceanographic changes in managed species and the whole ecosystem.

CONCLUSION
The time-series length as well as the spatial scale of CalCOFI sampling has shown that, overall, trends for warming, acidification, and deoxygenation appear to be more pronounced inshore than offshore, which suggests that coastal and nearshore ecosystems may be more vulnerable to rapid climate change. Thus, climate change effects may be felt more rapidly within state waters than federal waters, and nearshore marine protected areas may also experience more rapid change than offshore areas. While the CalCOFI record provides many examples of how marine communities in the SCCE respond to climate forcing, Koslow et al. (2018) caution that the impacts of secular climate change are unlikely to manifest as simply an extension of the warm phase of the ENSO or PDO cycles, and therefore past trends may not always be reliable guides for future ecosystem responses, which emphasizes the need for continued monitoring.
Monitoring programs will be essential to track the progression of anthropogenic climate change and understand associated ecosystem impacts. We suggest that existing oceanographic monitoring programs conduct an internal review to assess how the data collected from these programs can better contribute to stakeholder needs related to climate change adaptation. For countries working on developing new oceanographic monitoring programs to address climate issues, the following elements of CalCOFI have contributed in a valuable way to understanding the effects of climate change on the marine ecosystem: sampling inshore to offshore has revealed that nearshore trends are more pronounced, sampling below the thermocline in the upper 500 m of the water column has shown different manifestations of climate forcing, the collection of paired biological and hydrographic data allows for analysis of trends between environmental change and ecological responses, the positioning of the CalCOFI grid in a biogeographic boundary zone means that ecological communities are more responsive to climate change, and sampling across multiple trophic levels allows for an understanding of ecosystem functioning and food web impacts of climate change. The biological sampling strategy of CalCOFI is powerful in that early life history stages of many species are captured, allowing for a comparison of community trends through time for both fished and unfished species. However, the methodology also provides challenges in inferring adult habitat distributions, adult fitness, and recruitment success, which require complementary additional surveys to address.
Given the wealth of research and data along the US West Coast, this region is well primed to develop and test innovative strategies for climate change adaptation for marine ecosystem management and the CalCOFI Program is a key part of this process.

AUTHOR CONTRIBUTIONS
NG and BS conceptualized the manuscript. NG performed the analysis and wrote the manuscript. ED contributed to the climate projection analysis. All authors contributed to the manuscript revision, and read and approved the submitted version.

FUNDING
Funding for this study was provided by a NOAA QUEST grant to BS and NG and an NSF OCE1600283 grant supported ED.