Abstract
The Antarctic Peninsula (West Antarctica) marine ecosystem has undergone substantial changes due to climate-induced shifts in atmospheric and oceanic temperatures since the 1950s. Using 25 years of satellite data (1998-2022), this study presents evidence that phytoplankton biomass and bloom phenology in the West Antarctic Peninsula are significantly changing as a response to anthropogenic climate change. Enhanced phytoplankton biomass was observed along the West Antarctic Peninsula, particularly in the early austral autumn, resulting in longer blooms. Long-term sea ice decline was identified as the main driver enabling phytoplankton growth in early spring and autumn, in parallel with a recent intensification of the Southern Annular Mode (2010-ongoing), which was observed to influence regional variability. Our findings contribute to the understanding of the complex interplay between environmental changes and phytoplankton responses in this climatically key region of the Southern Ocean and raise important questions regarding the far-reaching consequences that these ecological changes may have on global carbon sequestration and Antarctic food webs in the future.
Similar content being viewed by others
Introduction
Marine ecosystems along the coast of the Antarctic Peninsula (West Antarctica; Fig. 1) have been undergoing significant change over the past few decades1,2,3. Among the most important is the rise of atmospheric and oceanic temperatures due to global warming, particularly during the winter1,4, which has led to the generalised retreat of glaciers, ice shelves and to a decrease in the extent, thickness, and volume of sea ice5,6,7,8,9. All marine trophic levels have been impacted by these changes2,3,10,11. Concerning primary producers, alterations in phytoplankton biomass, composition, and cell size have been reported, with potential consequences for the Antarctic food web2,6,12,13,14,15,16,17. Nevertheless, phytoplankton communities appear to be affected differently along the Western Antarctic Peninsula (WAP)2,12,13,14,15. In the southern-mid WAP (south of Anvers Island; 64°33â²S), studies have shown increases in phytoplankton biomass and cell size2,14,15,17, while in the northern Peninsula, the current consensus is that phytoplankton biomass and cell size are decreasing2,12,13,18. However, it should be noted that the phytoplankton dynamics in the northern part of the Peninsula have been less studied overall and are currently much less understood than the southern-mid WAP12.
While the continuous in-situ observation of marine phytoplankton in the WAP has been undertaken since the second half of the 20th century19, the latter was mostly restricted to coastal sites along the WAP in the proximity of active research stations and to oceanographic expeditions during the austral summer. This raises difficulties in understanding whether rapid biological changes are already taking place in the region. Ocean colour remote sensing provides a complementary tool to in-situ sampling, which has already proven to crucially contribute to the unravelling of short- and long-term phytoplankton patterns in the ocean at both regional and global scales20. Nonetheless, applications of satellite data to Antarctic coastal waters have been hampered by limited data availability due to persistent cloud cover, sea ice, and low solar angles, inter alia21,22. Additionally, global satellite algorithms have consistently underestimated chlorophyll-a (Chl-a; a proxy for phytoplankton biomass) in the Antarctic Peninsula by up to a factor of 222,23,24,25.
In this study, we implement the new OC4-SO regional algorithm for Chl-a, which was specifically calibrated for the Antarctic Peninsula using the largest high-quality in-situ regional dataset22 to advance existing knowledge of phytoplankton dynamics and the impact of climate change on primary production in this critical Southern Ocean region. We use 25 years of continuous multi-sensor remote sensing data (ESA OC-CCI26) to disentangle patterns of phytoplankton biomass and bloom phenology throughout the Antarctic Peninsula marine ecosystem in relation to regional atmospheric and oceanographic changes related to ongoing climate change. We observe enhanced phytoplankton biomass along the WAP, particularly in the early austral autumn, resulting in longer blooms. Long-term sea ice decline appears to be the main driver enabling growth in early spring and autumn, in parallel with a recent intensification of the Southern Annular Mode (SAM; 2010-ongoing).
Results and discussion
Phytoplankton biomass dynamics along the Antarctic Peninsula
Five marine subregions with distinctive seasonal phytoplankton patterns along the Antarctic Peninsula were identified after applying a hierarchical clustering analysis to key phytoplankton biomass and bloom phenology metrics (Fig. 2). From north to south: (i) DRA, which includes waters in the southern Drake Passage, north of the Southern Antarctic Circumpolar Current Front (SACCF; Fig. 2a), (ii) BRS, which extends from Elephant Island to the offshore waters south of Anvers Island, including the Bransfield Strait and the South Shetland Islands (Fig. 2b), (iii) WEDN, which is situated between the Bransfield Strait and the outer NW Weddell Sea, off the Northern Antarctic Peninsula (Fig. 2c), (iv) GES, which contains the coastal waters stretching south from the Gerlache Strait (Fig. 2d), and (v) WEDS, the NW Weddell coastal waters east of the Antarctic Peninsula, within the Larsen A and B embayments (Fig. 2e).
Phytoplankton blooms off the Antarctic Peninsula generally began in late October or November, with peaks in productivity occurring between December and February and ending in late March or April, coastal locations showing later blooms and higher Chl-a peaks (Fig. 2; Supplementary Table 1). DRA, the northernmost and most open ocean subregion, was where phytoplankton started to bloom earlier, right after austral winter. The low biomass levels and early short-lived biomass peak observed in this subregion suggest that micronutrient limitation (potentially iron) may occur shortly after December, corroborating previous in-situ studies that have reported lower biomass in the southern Drake Passage associated with iron limitation27,28, a characteristic of the more oceanic waters off the WAP during summer29. The nature of phytoplankton dynamics among subregions is also related to their physical environment. Regions within the WAP exhibited generally warmer waters, lower sea ice coverage, and stronger winds than the regions in the eastern sector of the Antarctic Peninsula (Supplementary Fig. 1; Supplementary Table 2). Sea ice concentration was the primary factor influencing the timing of bloom initiation around the Peninsula (Supplementary Fig. 2; p valueâ<â0.05), as regions with more sea ice during spring were often associated with later blooms (e.g., WEDN, GES, and WEDS). The importance of environmental conditions in shaping phytoplankton dynamics was also evidenced in BRS, as seen by its high variability in both bloom phenology and biomass. Despite its relatively low sea ice concentration, the Bransfield Strait, at the centre of BRS, exhibits a complex current system that includes cyclonic and anticyclonic eddies, surface and subsurface thermal fronts, as well as intrusions of warmer waters from the west (the Circumpolar Deep Water; CDW) and colder waters from the northeast (from the NW Weddell Sea)30,31,32,33, an overall complexity that translates to high spatial and temporal variability.
Enhanced biomass and longer blooms in the WAP
Results from a linear trend analysis indicated that the mean phytoplankton biomass from September to April (austral spring to autumn) has significantly increased between 1998-2022 in the WAP, particularly in the Bransfield Strait (p valueâ<â0.05; Fig. 3a, Supplementary Fig. 3). At the same time, it is not clear from our results whether biomass also changed in the eastern sector of the Antarctic Peninsula, which may be a consequence of the low number of satellite observations in this sector and lower accuracy caused by the high sea ice coverage34. The phytoplankton biomass increase observed in the WAP seems to have occurred mainly in the early spring and early autumn (Table 1). This is evidenced by the increases of +0.003âmg mâ3 Chl-a yearâ1 in the DRA and of +0.009âmg mâ3 Chl-a yearâ1 in the BRS during September (early spring). This was even more the case during early autumn (March and April), when biomass trends were nearly three times higher compared with the early spring: +0.009âmg·mâ3 Chl-a yearâ1 in DRA, +0.026âmg mâ3 Chl-a yearâ1 in BRS, and +0.061âmg mâ3 Chl-a yearâ1 in GES. Therefore, most of the biomass increase in the WAP occurred in early autumn, which corroborates the higher average autumn biomass observed in 2011â2020 compared to 2001â2010 (Fig. 2a, b, d).
The increase in biomass seen in the GES partially corroborates the findings of previous in-situ and satellite studies, although for different temporal periods2,14,15,16,17,35. For instance, Moreau et al17. observed an increase in the annual integrated primary production from 1997 to 2010 using exclusively satellite data17 but found no clear trend in satellite Chl-a. These authors attributed this increase in primary production to the long-term trend towards earlier sea ice retreats and its effects on phytoplankton growth during the spring17. Therefore, they did not consider the large phytoplankton growth in autumn that we observed here, which cannot be related to the earlier sea ice retreat. However, Turner et al. using satellite Chl-a from 1997 to 2022, also observed an increase in phytoplankton biomass in the autumn in northern WAP16. Previous studies using in-situ data, such as Montes-Hugo et al2. (from 1978 to 2006) and Brown et al.15 (from 1993 to 2017), also observed a similar decadal enhancement of phytoplankton biomass in this region, yet, unlike our satellite-detected trends, they only observed it for the summer (the only period for which the current availability of in situ data allows for more accurate long-term trend assessments). These studies suggested that this biomass increase was a consequence of a more stable upper mixed layer during the summer due to the meltwater resulting from climate-induced sea ice decline and glacier retreat2,15. While our region-wise satellite observations do not show a clear increasing trend of biomass during summer in the mid/southern WAP (GES), the average summer biomass between 2011 and 2020 does appear to have increased in comparison to the previous decade (Fig. 2d). Moreover, a slight increasing biomass trend was also observed for December (although only at a p valueâ=â0.09; Table 1) and pixel-wise trend analyses for the summer do suggest that biomass might be increasing in coastal waters south of the Gerlache Strait, although the significant increase is locally observed (Supplementary Fig. 3). A reanalysis of in-situ Chl-a data collected along the WAP (from Palmer LTER and the Brazilian High Latitude Oceanography Group) also suggests that summer biomass has increased in the BRS (â+â0.06âm ·mâ3 yearâ1) and GES (â+â0.18âmg mâ3 yearâ1) (Fig. 3b; Supplementary Fig. 4). Therefore, there may have been, in truth, an increase in summer biomass that satellite data did not accurately capture.
Our findings also help to shed some light on the current understanding of the northern WAP (DRA and BRA regions). Our results contradict the hypothesis proposed by Montes-Hugo et al2. in 2009, who suggested that regional phytoplankton biomass is decreasing due to the deepening of the upper mixed layer caused by less sea ice concentration and enhanced wind stress2, although it is important to note that this study compared the period between 1978â1986 and 1998â2006. Other phytoplankton ecology studies have since echoed this hypothesis, although no further study has confirmed it using long-term datasets. Here, however, we demonstrate that phytoplankton biomass is increasing since 1998, an increase concurrent with wind stress (Supplementary Fig. 5), which suggests that phytoplankton growth has not been hindered by a long-term deepening of the mixed layer resulting from stronger winds35. Our observations is supported by recent studies reporting localized long-term increases in summer biomass in Potter Cove (King George Island; 62.2382°S, 58.6673°W) and the region between the South Shetland Islands and the Elephant and Clarence islands, although neither of these works report increases in autumn biomass36,37. Contrarily, Turner et al. observed a recent increase in autumn biomass in the northernmost waters of the WAP, yet did not report long-term increase in annual Chl-a16.
The changes in phytoplankton biomass have led to changes in bloom phenology (Table 2), with blooms ending later in two of the three studied sub-regions of the WAP (â+â0.4 weeks yearâ1 in GES and +0.1 weeks yearâ1 in DRA). As a result, blooms have become longer in both regions, thus resulting in higher total biomass production, particularly at the DRA (p valueâ<â0.05). While the timing of the bloom termination changed, there were no apparent shifts in the timings of bloom initiation or bloom peak. Surprisingly, no trends were observed for BRS despite also being a region where biomass increased significantly in the autumn. This was likely due to the high interannual variability associated with bloom timing in this region, which might prevent the identification of linear trends within the 25 years of our dataset. Previous long-term studies on phytoplankton bloom phenology in the WAP are rare due to the scarcity of data that characterizes the region. Using +20-year data (1992-2016) from three different stations along the WAP, Kim et al37. showed that the northernmost station (Carlini Station; BRS region) exhibited a longer and less pronounced (i.e., with a lower biomass peak) bloom than the other stations (Palmer and Rothera)37, in line with our results. However, these authors did not assess long-term changes in bloom timing. Recently, Thomalla et al38. suggested that there has been an overall tendency between 1998 and 2022 for blooms in the more heavily ice-influenced regions of the Southern Ocean to initiate earlier and have longer durations38, which contradicts our results for the WAP. It is important to mention that Thomalla et al. focused on the entire Southern Ocean, thus being less detailed at a regional scale, unlike the present study for the WAP. Another recently published study, focused on the WAP, also reported that phytoplankton blooms in the marginal ice zone and the continental shelf have shifted later16, corroborating our results and reinforcing the importance of regional studies for the understanding of phytoplankton phenology in complex regions.
Sea ice decline is enabling phytoplankton growth
The changes observed in phytoplankton biomass and bloom dynamics in spring and autumn can be mainly attributed to the long-term decline in sea ice coverage seen along the WAP over the past four decades, particularly during the late autumn â early winter transition5,6,7,8,9,17. A decrease in the mean sea ice extent was seen in both spring and autumn periods from 1982â1997 to 1998â2022 for each region of the WAP (DRA, BRS, and GES; Fig. 4a; Supplementary Fig. 6). A key consequence of this decline is that sea ice retreat (in spring) and advance (in autumn) are occurring increasingly earlier and later, respectively, leading to an increase in the number of sea ice-free days (Fig. 4a). This allows for enhanced phytoplankton growth during spring and autumn in areas previously occupied by sea ice since the additional space and the improved light conditions are now more favourable for phytoplankton accumulation. The observed sea ice decline across the entire WAP likely explains why the phytoplankton growth in the autumn was widespread along the region.
On average, the autumn phytoplankton biomass in the WAP increased from 29% to 42% of the total biomass from September to April, while spring and summer biomass contribution slightly decreased (Fig. 4b). Several reasons can explain this difference. Firstly, the rate at which sea ice advance is becoming later is higher than the rate at which sea ice retreat is becoming earlier5, i.e., the amount of the new ice-free days caused by the long-term warming is higher in the autumn than in the spring. Secondly, the advance of sea ice in the WAP is delayed by persistent northerly winds that occur during the autumn5, meaning that sea ice coverage tends to be lower in the autumn than in the early spring. As such, there is more new space for phytoplankton to grow and accumulate biomass in the autumn than in the spring. Moreover, the lower sea ice coverage in the autumn may allow for higher light penetration in the upper layers of the water column12, promoting phytoplankton growth in waters where biomass could not previously accumulate under the sea ice. Finally, phytoplankton biomass typically reaches much higher concentrations in the late summer compared to the early spring (Fig. 3a, b, d). Apart from the higher temperatures and lower sea ice coverage, a key reason for this is that the phytoplankton biomass already accumulated in the late summer can be maintained throughout the early autumn if enough nutrients are available, while phytoplankton in the early spring must grow from the typically very low winter biomass.
These findings put into perspective the future trends for sea ice extent in Antarctica. In the long term, sea ice is expected to decrease as warming continues39, similar to what is also projected for the Arctic40. The resulting decrease in the duration of the sea ice season is likely to be followed by enhanced phytoplankton growth at the limits of the sea ice season (spring and autumn). Nevertheless, even if all required nutrients for growth are available, the duration of phytoplankton blooms will ultimately be limited by the low light availability that characterizes the austral winter in Antarctica. In the Arctic, Ardyna et al.41 showed that recent sea ice loss have triggered fall blooms analogous to what we observe in the WAP, with wind-driven vertical mixing playing an important role in supporting phytoplankton growth. More high-resolution in-situ studies focused on the austral spring and, particularly, autumn are essential for a better understanding of how phytoplankton dynamics will change. In addition, complementary sampling techniques to in-situ sampling expeditions such as multiparametric (biogeochemical) floats and underwater gliders are becoming more frequent and are expected to contribute significantly to the acquisition of in-situ observations from other seasons besides summer42,43,44.
SAM as a key driver of regional variability
While the overall declining trend in sea ice explains most of the findings observed in our work, it does not explain some of the regional differences observed along the WAP. For instance, it does not explain why the more statistically significant pixel-wise biomass trends were found in the northernmost and more offshore areas of the WAP (Fig. 2a). This can be mainly attributed to an intensification in the Southern Annular Mode, a key climate mode modulating the marine environment of the WAP, particularly its northernmost sector. Under positive SAM conditions, circumpolar westerlies strengthen, bringing warmer air into the WAP and increasing cloud coverage45. In the marine ecosystems of the WAP, positive SAM conditions are typically associated with warmer surface temperatures and reduced sea ice36. Since the 1950s, the SAM index has shown a generally upward trend linked to anthropogenic climate change, although a plateau was observed between the late 1990s and the 2000s36,45,46. Since 2010, however, the SAM index has intensified again and is becoming increasingly positive at a fast rate, as seen by the large increase in the cumulative SAM index (Fig. 5a).
We hypothesize that this recent intensification of SAM may have contributed to drive the observed region increase in phytoplankton biomass as well as the reason why this was more prevalent in the northwestern Peninsula. We observed a positive and significant pixel-wise correlation (considering the average between September-April period) between the average Chl-a concentration and SAM index in the waters along most of the Western Antarctic Peninsula, particularly for the southern Drake Passage and the northern tip of the Antarctic Peninsula (Fig. 5b). As previously mentioned, wind intensification associated with positive SAM results in deeper mixed layers along the Antarctic Peninsula, especially in offshore waters38. The increased vertical mixing caused by stronger winds could be increasing micronutrient availability, particularly iron (Fe), for the DRA and the non-sheltered areas of the BRS, which are less affected by sea ice conditions and have a stronger oceanic influence. This would help prevent nutrient limitation in these areas, especially in the autumn, when a greater amount of nutrients has already been consumed27,28, coinciding with the period when most of the biomass enhancement was observed in our study (Table 1). Accordingly, Tagliabue et al.47 observed that Fe supply from winter mixing is crucial to offshore primary production in the Southern Ocean. Furthermore, the stronger westerlies associated with positive spring SAM can also enhance warm, nutrient-rich CDW intrusions over the shelf during summer36,48,49. Such intrusions propagate along the shelf, eventually reaching the northern Antarctic Peninsula and its more offshore areas36. Several studies have already shown high phytoplankton biomass linked to CDW intrusions48,49, suggesting that these ocean circulation processes coupled with positive SAM trends could be another contributing factor to the increased biomass in these regions.
Contrastingly, we do not observe a positive correlation between SAM and Chl-a in the southernmost coastal WAP (GES; Fig. 5b). Enhanced CDW intrusions have been reported to force and increase glacier retreat due to its warmer temperatures, leading to more frequent and stabilized water column structures during spring and summer months15. While such changes in the upper ocean physical compartments may have initially led to an increase in summer biomass15, the recent intensification of SAM could be driving shifts in phytoplankton community composition and size structure and, consequently, on biomass accumulation. Recent studies along the WAP have reported frequent shifts from large centric diatoms to smaller-sized cryptophytes in regions under localized glacier meltwater input during the austral summer13,18,50,51,52. As a result, this shift could decrease the total summer phytoplankton biomass accumulation under positive SAM periods, explaining the slightly negative correlation between SAM and Chl-a for the coastal regions observed here. Moreover, it will also be essential to consider the influence the El-Niño Southern Oscillation can have on the interannual dynamics of phytoplankton, as summers following extreme El-Niño events have already been associated with large diatom blooms53,54.
A preferential feeder on large diatoms55,56,57, Antarctic krill (Euphasia superba) is a key species in the Southern Ocean ecosystem and an essential link between phytoplankton and higher trophic levels such as whales, seals, and penguins4,11,12. Therefore, localized shifts from diatoms to cryptophytes could likely contribute to shaping the phytoplankton dynamics along the coastal WAP by alleviating the grazing pressure from the krill on phytoplankton biomass accumulation, while also potentially impacting the krill stocks. Nevertheless, it is still unclear what impacts such changes in phytoplankton community composition and size structure may have on the WAP marine ecosystems. Further in situ studies are required to understand if a phytoplankton community shift to cryptophytes will significantly weaken biological carbon uptake or lead to structural changes in the WAP coastal marine food web.
Methods
A summary of all datasets used, including their temporal coverage, spatial and temporal resolution, units, product name and relevant references is available in the Supplementary Table 3. The following subsections present further detail on each dataset and presents all methodologies used. All analysis were performed in Python 3.8.8.
Satellite data
Daily satellite remote sensing reflectance data with 4âkm spatial resolution was extracted from the ESA Ocean Colour Climate Change Initiative (OC-CCI) v6.0 product26 for 1998-2022. Chlorophyll-a (mg mâ3) was estimated by applying the OC4-SO algorithm22, a regional Chlorophyll-a algorithm specifically devised for the Antarctic Peninsula region to increase the accuracy of Chlorophyll-a measurements22. Satellite sea surface temperature (°C) and sea ice concentration (%) fields with daily temporal resolution and ~5âkm spatial resolution were acquired from the Operational Sea Surface Temperature and Ice Analysis (OSTIA)58, available at the Copernicus Marine Environment Monitoring Service (https://marine.copernicus.eu/). The sea ice concentration was used to calculate daily sea ice extent (i.e., the area occupied by sea ice) for the clusters DRA, BRS and GES by multiplying the area of each pixel (~25âkm2) by the number of pixels covered by sea ice (using a minimum threshold of 15% of sea ice concentration per pixel5). For each cluster, the annual sea ice duration, day of sea ice advance, and day of sea ice retreat were calculated following Stammerjohn et al5. Following this definition, the day of sea ice advance corresponded to the first day in which sea ice concentration is higher than 15% for at least 5 consecutive days. The day of sea ice retreat is the inverse, i.e., the first day after the winter when the sea ice concentration drops below 15% (a 5 consecutive day window is also used)5.
Daily, 4âkm resolution photosynthetically active radiation (PAR; Einstein mâ2 dâ1) data were retrieved from the ESA GlobColour project59. Monthly Southern Annular Mode (SAM) index values were retrieved from the NOAA Centre for Weather and Climate Prediction Centre (https://www.cpc.ncep.noaa.gov/).
In-situ data
HPLC-derived in-situ surface measurements of Chl-a collected across the Western Antarctic Peninsula were used. The full dataset (Nâ=â4322) included published data available in Valente et al.60 and & Palmer LTER61, as well as unpublished data collected by the Brazilian High-Latitude Oceanography Group (GOAL-FURG).
The GOAL-FURG in-situ Chl-a measurements were collected during 12 austral summer expeditions to the WAP aboard the Brazilian vessels NP Almirante Maximiano and NP Ary Rongel (Nâ=â533). Seawater samples (ranging from 0.5 to 2.5âL) were filtered under low vacuum using GF/F filters (25âmm diameter, 0.7 μm pore size) and subsequently stored at â80â°C. Prior to analysis, the filters were placed in screw-cap centrifuge tubes with 3âmL of 95% methanol (2% ammonium acetate) containing 0.05 mg Lâ1 trans-β-apo-8â-carotenal as an internal standard. Samples were sonicated for 5âmin, stored at â20â°C for one hour and then centrifuged at 1100âÃâg for 5âmin at 3â°C. The supernatants were filtered through PTFE membrane filters (0.2âµm pore size) to ensure no residues were inserted in the HPLC system. For pigment analysis, a monomeric C8 column and a mobile phase containing pyridine was used. Before injection, 1âmL of sample was mixed with 0.4âmL of Milli-Q water in 2.0âmL glass sample vials and placed in the HPLC cooling rack of a Shimadzu Prominence LC-20A Modular HPLC System. The methodâs detection and quantification limits have been calculated and discussed by Mendes et al62. Pigment standards from DHI (Institute for Water and Environment, Denmark) were used for calibration of the HPLC system. Pigments were identified based on their absorbance spectra and retention times from the signals in the photodiode array detector (SPD-M20A; 190â800ânm; 1ânm wavelength accuracy) or fluorescence detector (RF-10AXL; Ex. 430ânm/Em.670ânm). Peaks were integrated using the LC-Solution software, but all peak integrations were manually verified and corrected if necessary. A quality assurance threshold procedure, through the application of limits of quantification and detection, was applied to the pigment data to reduce the uncertainty of pigments found in low concentrations63. Pigment concentrations were normalized to the internal standard to account for losses and volume changes.
Combined, the three databases result in a comprehensive dataset of in situ Chl-a measurements spanning from 1998 to 2020. In-situ atmospheric wind speed measured at Palmer Station were extracted from the âPalmer Station WeatherâDaily Averages (#28)â dataset, available in EDI64. In-situ atmospheric wind speed measured at Marambio and Great Wall stations were extracted from the Global Historical Climatology Network daily (GHCNd)65.
Clustering areas with coherent phytoplankton seasonal patterns
Due to the extent and complexity of the marine environments off the Antarctic Peninsula, an agglomerative hierarchical clustering analysis was first used to cluster areas within the Antarctic Peninsula with coherent phytoplankton seasonal patterns using remote sensing data from 1998 to 2022. First, for each pixel, the following metrics were calculated: (i) the first month when phytoplankton growth (i.e., satellite Chl-a) is, on average, detected after the winter, (ii) the last month when phytoplankton growth is, on average, detected after the summer, (iii) the average month when the peak satellite Chl-a is detected, and (iv) the numerical integration of mean Chl-a concentrations during the austral summer (November-February), calculated using Simpsonâs rule. All metrics were normalized. Finally, the spatial resolution of the pixels was reduced to 10âkm to avoid performance issues when running the clustering analysis. The distance and linkage metrics used for clustering were the Manhattan distance and average linkage, respectively. Four clusters were initially identified. After closer inspection, one cluster was divided in two due to a partial influence of sea ice, which could affect phytoplankton bloom phenology, thus bringing the total number of clusters to 5. All pixels that were, on average, covered with sea ice during the entire year (i.e., sea ice concentration above 15%) were discarded. Subsequently, pixels indicative of noise (e.g., pixels of a cluster located inside another cluster without a logical reason) were manually removed. Finally, the spatial resolution of the resulting map of clusters was upscaled to 4âkm. All analyses were performed using Scipy and Sklearn modules in Python 3.8.8.
Phytoplankton bloom phenology
For each of the five identified clusters, phytoplankton bloom phenology metrics were derived from satellite Chl-a. First, daily satellite Chl-a was temporally averaged to an 8-day week frequency and spatially averaged within each cluster. Subsequently, the phytoplankton austral spring-summer bloom was identified as the event, no shorter than 2 weeks, in which Chl-a surpasses a threshold of 5% of the annual Chl-a median66,67,68. Five bloom phenology metrics were then calculated. The bloom initiation and bloom termination dates were defined as the first and last week of the bloom, respectively (i.e., when Chl-a rises above and subsequently falls below the 5% threshold). Bloom duration corresponds to the number of weeks between bloom initiation and bloom termination. The bloom peak date is the week within the blooming period when the maximum Chl-a value is found.
Analysis of phytoplankton trends
To understand if phytoplankton biomass and bloom phenology metrics in the Western Antarctic Peninsula have changed since 1998, the time series for each previously defined cluster were assessed for statistically significant trends or environmental shifts. For each year, the mean Chl-a during each month between September and April (full austral spring-summer period) was calculated, as well as the mean Chl-a during this entire period (September-April) and during only the peak of the summer (December-February). The regions within the eastern sector of the Antarctic Peninsula (WEDi and WEDo) were excluded from these analyses since satellite observations are very limited due to the persistence of sea ice that characterizes this sector. Since low data availability and sea ice adjacency is likely to affect the validity of trend analyses. Linear trends were estimated by fitting linear regressions to each time series and extracting the slope of the resulting line. Outliers were excluded by removing all Chl-a values over (under) the mean Chl-aâ+â(â) 3 standard deviations. Pixel-wise linear regressions were also performed for the mean Chl-a during September-April in order to produce maps showcasing the linear trend for Chl-a. In order to understand if a similar trend was observed in in-situ data, the mean Chl-a during September-April was also calculated using only in-situ data from the GES and BRS regions (the only regions of the WAP where data was abundant enough to accurately run trend analyses). To understand how the trends in phytoplankton biomass may have changed between decades in the WAP, the average Chl-a concentration for each season (spring, summer, and autumn) was calculated for 2001â2010 and for 2011â2020. Data from the three regions of the WAP (DRA, BRS and GES) were used for this comparison. The proportion of the Chl-a attributed to each season to the total Chl-a (the sum of the average spring, summer, and autumn mean Chl-a) was then calculated. Finally, to better understand the relationship between Chl-a and SAM, a pixel-wise linear regression was performed for the full WAP region.
Reporting summary
Further information on research design is available in the Nature Portfolio Reporting Summary linked to this article.
Data availability
The previously unpublished in-situ chlorophyll-a dataset from GOAL utilized in this study is now accessible at https://doi.org/10.5281/zenodo.12580624. All other data utilized in this study have been previously published and are available online. The OC-CCI satellite chlorophyll-a product is available at https://www.oceancolour.org/. The in-situ chlorophyll-a data acquired from Valente et al60. are available at https://doi.org/10.5194/essd-14-5737-2022. The in-situ chlorophyll-a data acquired from Palmer-LTER are available https://doi.org/10.6073/pasta/9d0c1561ca8c5540227df6efa37c61b5. The OSTIA product, containing sea surface temperature and sea ice concentration data is available at https://doi.org/10.48670/moi-00168. The GlobColour PAR data is available at https://hermes.acri.fr/. The wind data from ERA5 Reanalysis are available at https://doi.org/10.24381/cds.adbb2d47. The in-situ wind data from Palmer Station are available at https://doi.org/10.6073/pasta/3eefb45dbfb784c3cabe3690ea46fe9e. The in-situ wind datasets from Marambio Station and Great Wall Station are available at https://www.ncei.noaa.gov/products/land-based-station/global-historical-climatology-network-daily. The Southern Annular Mode (SAM) data is available at https://www.cpc.ncep.noaa.gov/. Source data are provided with this paper.
Code availability
All custom code developed within this work is available at: https://github.com/afonsomferreira/antarcticpeninsula-trends.
References
Meredith, M. P. & King, J. C. Rapid climate change in the ocean west of the Antarctic Peninsula during the second half of the 20th century. Geophys. Res. Lett. 32, L19604 (2005).
Montes-Hugo, M. et al. Recent changes in phytoplankton communities associated with rapid regional climate change along the Western Antarctic Peninsula. Science 323, 1470â1473 (2009).
Henley, S. F. et al. Changing biogeochemistry of the Southern Ocean and its ecosystem implications. Front. Marine Sci. 7, 581 (2020).
Turner, J. et al. Antarctic climate change during the last 50 years. Int. J. Climatol. 25, 279â294 (2005).
Stammerjohn, S. E., Martinson, D. G., Smith, R. C. & Iannuzzi, R. A. Sea ice in the western Antarctic Peninsula region: Spatio-temporal variability from ecological and climate change perspectives. Deep Sea Res. Part II: Topical Stud. Oceanogr. 55, 2041â2058 (2008).
Lin, Y. et al. Decline in plankton diversity and carbon flux with reduced sea ice extent along the Western Antarctic Peninsula. Nat. Commun. 12, 4948 (2021).
Mar Flexas, M. et al. Antarctic Peninsula warming triggers enhanced basal melt rates throughout West Antarctica. Sci. Adv. 8, eabj9134 (2022).
Andreasen, J. R. et al. Change in Antarctic ice shelf area from 2009 to 2019. Cryosphere 17, 2059â2072 (2023).
Davison, B. J. et al. Widespread increase in discharge from West Antarctic Peninsula glaciers since 2018, EGUsphere [preprint] (2024).
Ducklow, H. W. et al. Marine pelagic ecosystems: the West Antarctic Peninsula. Philos. Trans. R. Soc. B: Biol. Sci. 362, 67â94 (2006).
Atkinson, A. et al. Krill (Euphausia superba) distribution contracts southward during rapid regional warming. Nat. Clim. Change 9, 142â147 (2019).
Ferreira, A. et al. Changes in phytoplankton communities Along the Northern Antarctic Peninsula: Causes, impacts and research priorities. Front. Marine Sci. 7, 576254 (2020).
Mendes, C. R. B. et al. Shifts in the dominance between diatoms and cryptophytes during three late summers in the Bransfield Strait (Antarctic Peninsula). Polar Biol. 36, 537â547 (2013).
Schofield, O. et al. Decadal variability in coastal phytoplankton community composition in a changing West Antarctic Peninsula. Deep Sea Res. Part I: Oceanographic Res. Pap. 124, 42â54 (2017).
Brown, M. S. et al. Enhanced oceanic CO2 uptake along the rapidly changing West Antarctic Peninsula. Nat. Clim. Chang 9, 678â683 (2019).
Turner, J. S. et al. Changing phytoplankton phenology in the marginal ice zone west of the Antarctic Peninsula. Mar. Ecol. Prog. Ser. 734, 1â21 (2024).
Moreau, S. et al. Climate change enhances primary production in the western Antarctic Peninsula. Glob. Change Biol. 21, 2191â2205 (2015).
Mendes, C. R. B. et al. Cryptophytes: An emerging algal group in the rapidly changing Antarctic Peninsula marine environments. Glob. Change Biol. 29, 1791â1808 (2023).
Boyce, D. G., Lewis, M. R. & Worm, B. Global phytoplankton decline over the past century. Nature 466, 591â596 (2010).
Groom, S. et al. Satellite ocean colour: Current status and future perspective. Front. Marine Sci. 6, (2019).
Pope, A. et al. Community review of Southern Ocean satellite data needs. Antarct. Sci. 29, 97â138 (2017).
Ferreira, A. et al. OC4-SO: A new chlorophyll-a algorithm for the Western Antarctic peninsula using multi-sensor satellite data. Remote Sens. 14, 1052 (2022).
Dierssen, H. M. & Smith, R. C. Bio-optical properties and remote sensing ocean color algorithms for Antarctic Peninsula waters. J. Geophys. Res.: Oceans 105, 26301â26312 (2000).
Szeto, M., Werdell, P. J., Moore, T. S. & Campbell, J. W. Are the worldâs oceans optically different? J. Geophys. Res.: Oceans 116, C00H04 (2011).
Johnson, R., Strutton, P. G., Wright, S. W., McMinn, A. & Meiners, K. M. Three improved satellite chlorophyll algorithms for the Southern Ocean. J. Geophys. Res.: Oceans 118, 3694â3703 (2013).
Sathyendranath, S. et al. An ocean-colour time series for use in climate studies: The experience of the ocean-colour climate change initiative (OC-CCI). Sensors 19, 4285 (2019).
Trimborn, S., Hoppe, C. J. M., Taylor, B. B., Bracher, A. & Hassler, C. Physiological characteristics of open ocean and coastal phytoplankton communities of Western Antarctic Peninsula and Drake Passage waters. Deep Sea Res. Part I: Oceanogr. Res. Pap. 98, 115â124 (2015).
Balaguer, J., Koch, F., Hassler, C. & Trimborn, S. Iron and manganese co-limit the growth of two phytoplankton groups dominant at two locations of the Drake Passage. Commun. Biol. 5, 1â12 (2022).
Annett, A. L. et al. Controls on dissolved and particulate iron distributions in surface waters of the Western Antarctic Peninsula shelf. Mar. Chem. 196, 81â97 (2017).
Sangrà , P. et al. The Bransfield current system. Deep Sea Res. Part I: Oceanogr. Res. Pap. 58, 390â402 (2011).
Sangrà , P. et al. The Bransfield gravity current. Deep Sea Res. Part I: Oceanogr. Res. Pap. 119, 1â15 (2017).
Dotto, T. S., Kerr, R., Mata, M. M. & Garcia, C. A. E. Multidecadal freshening and lightening in the deep waters of the Bransfield Strait, Antarctica. J. Geophys. Res.: Oceans 121, 3741â3756 (2016).
Yuri Damini, B. et al. Antarctica Slope Front bifurcation eddy: A stationary feature influencing CO2 dynamics in the northern Antarctic Peninsula. Prog. Oceanogr. 212, 102985 (2023).
Bélanger, S., Ehn, J. K. & Babin, M. Impact of sea ice on the retrieval of water-leaving reflectance, chlorophyll a concentration and inherent optical properties from satellite ocean color data. Remote Sens. Environ. 111, 51â68 (2007).
Schofield, O. et al. Changes in the upper ocean mixed layer and phytoplankton productivity along the West Antarctic Peninsula. Philos. Trans. R. Soc. A: Math., Phys. Eng. Sci. 376, 20170173 (2018).
Zhang, Z. et al. Linkage of the physical environments in the northern Antarctic Peninsula region to the Southern Annular Mode and the implications for the phytoplankton production. Prog. Oceanogr. 188, 102416 (2020).
Kim, H. et al. Inter-decadal variability of phytoplankton biomass along the coastal West Antarctic Peninsula. Philos. Trans. R. Soc. A: Math., Phys. Eng. Sci. 376, 20170174 (2018).
Thomalla, S. J. et al. Widespread changes in Southern Ocean phytoplankton blooms linked to climate drivers. Nat. Clim. Chang. 13, 975â984 (2023).
Siegert, M. et al. The Antarctic Peninsula Under a 1.5â°C Global Warming Scenario. Front. Environ. Sci. 7, 102 (2019).
Parkinson, C. L. & DiGirolamo, N. E. Sea ice extents continue to set new records: Arctic, Antarctic, and global results. Remote Sens. Environ. 267, 112753 (2021).
Ardyna, M. et al. Recent Arctic Ocean sea ice loss triggers novel fall phytoplankton blooms. Geophys. Res. Lett. 41, 6207â6212 (2014).
Thomalla, S. J., Ogunkoya, A. G., Vichi, M. & Swart, S. Using optical sensors on gliders to estimate phytoplankton carbon concentrations and chlorophyll-to-carbon ratios in the Southern Ocean. Front. Mar. Sci. 4, 34 (2017).
Moreau, S., Boyd, P. W. & Strutton, P. G. Remote assessment of the fate of phytoplankton in the Southern Ocean sea-ice zone. Nat. Commun. 11, 3108 (2020).
Newman, L. et al. Delivering sustained, coordinated, and integrated observations of the southern ocean for global impact. Front. Mar. Sci. 6, 433 (2019).
Fogt, R. L. & Marshall, G. J. The Southern Annular Mode: Variability, trends, and climate impacts across the Southern Hemisphere. WIREs Clim. Change 11, e652 (2020).
Turner, J. et al. Absence of 21st century warming on Antarctic Peninsula consistent with natural variability. Nature 535, 411â415 (2016).
Tagliabue, A. et al. Surface-water iron supplies in the Southern Ocean sustained by deep winter mixing. Nat. Geosci. 7, 314â320 (2014).
Prezelin, B. B., Hofmann, E. E., Mengelt, C. & Klinck, J. M. The linkage between Upper Circumpolar Deep Water (UCDW) and phytoplankton assemblages on the west Antarctic Peninsula continental shelf. J. Mar. Res. 58, 165â202 (2000).
Clem, K. R., Renwick, J. A., McGregor, J. & Fogt, R. L. The relative influence of ENSO and SAM on Antarctic Peninsula climate. J. Geophys. Res.: Atmospheres 121, 9324â9341 (2016).
Mascioni, M., Almandoz, G. O., Ekern, L. & Pan, B. J. & Vernet. Microplanktonic diatom assemblages dominated the primary production but not the biomass in an Antarctic fjord. J. Mar. Syst. 224, 103624 (2021).
Mendes, C. R. B. et al. New insights on the dominance of cryptophytes in Antarctic coastal waters: A case study in Gerlache Strait. Deep Sea Res. Part II: Topical Stud. Oceanogr. 149, 161â170 (2018).
Costa, R. R. et al. Physical-biological drivers modulating phytoplankton seasonal succession along the Northern Antarctic Peninsula. Environ. Res. 231, 116273 (2023).
Costa, R. R. et al. Dynamics of an intense diatom bloom in the Northern Antarctic Peninsula, February 2016. Limnol. Oceanogr. 65, 2056â2075 (2020).
Costa, R. R. et al. Large diatom bloom off the Antarctic Peninsula during cool conditions associated with the 2015/2016 El Niño. Commun. Earth Environ. 2, 1â11 (2021).
Santa Cruz, F., Ernst, B., Arata, J. A. & Parada, C. Spatial and temporal dynamics of the Antarctic krill fishery in fishing hotspots in the Bransfield Strait and South Shetland Islands. Fish. Res. 208, 157â166 (2018).
Quetin, L. B. & Ross, R. M. Feeding by Antarctic krill, Euphausia superba: Does size matter? In Antarctic nutrient cycles and food webs 372-377. Berlin, Heidelberg: Springer Berlin Heidelberg (1985).
Haberman, K. L., Ross, R. M. & Quetin, L. B. Diet of the Antarctic krill (Euphausia superba Dana): II. Selective grazing in mixed phytoplankton assemblages. J. Exp. Mar. Biol. Ecol. 283, 97â113 (2003).
Good, S. et al. The current configuration of the OSTIA system for operational production of foundation sea surface temperature and ice concentration analyses. Remote Sens. 12, 720 (2020).
Frouin, R. J., Franz, B. A. & Werdell, J. The SeaWiFS PAR product. In Algorithm Updates for the Fourth SeaWiFS Data Reprocessing, (eds Hooker, S. B. & Firestone, E. R.). Vol 22, 46-50 CC NASA/TM-2003-206892, (2003).
Valente, A. et al. A compilation of global bio-optical in situ data for ocean-colour satellite applications â version three. Earth Syst. Sci. Data 14, 5737â5770 (2022).
Palmer Station Antarctica LTERO. Schofield, M. Vernet, and B. Prezelin. 2022. Photosynthetic pigments of water column samples and analyzed with High Performance Liquid Chromatography (HPLC), collected aboard Palmer LTER annual cruises off the coast of the Western Antarctica Peninsula, 1991 â 2020. ver 9. Environmental Data Initiative. https://doi.org/10.6073/pasta/9d0c1561ca8c5540227df6efa37c61b5 (Accessed 2024-04-23).
Mendes, C. R., Cartaxana, P. & Brotas, V. HPLC determination of phytoplankton and microphytobenthos pigments: comparing resolution and sensitivity of a C18 and a C8 method. Limnol. Oceanogr. Methods 5, 363â370 (2007).
Hooker, S. B. et al. The Fifth SeaWiFS HPLC Analysis Round-Robin Experiment (SeaHARRE-5). NASA Technical Memorandum 2012-217503, NASA Goddard Space Flight Center, Greenbelt, Maryland (2012).
Palmer Station Antarctica LTER. Daily weather averages for Palmer Station, Antarctica (1989-2023) ver 9. Environmental Data Initiative. https://doi.org/10.6073/pasta/3eefb45dbfb784c3cabe3690ea46fe9e (Accessed 2023-11-14).
Menne, M. J., Durre, I., Vose, R. S., Gleason, B. E. & Houston, T. G. An overview of the global historical climatology network-daily database. J. Atmos. Ocean. Technol. 29, 897â910 (2012).
Brody, S. R., Lozier, M. S. & Dunne, J. P. A comparison of methods to determine phytoplankton bloom initiation. J. Geophys. Res. Oceans 119, 2345â2357 (2013).
Cole, H., Henson, S., Martin, A. & Yool, A. Mind the gap: the impact of missing data on the calculation of phytoplankton phenology metrics. J. Geophys. Res. Oceans 117, C08030 (2012).
Ferreira, A., Brotas, V., Palma, C., Borges, C. & Brito, A. C. Assessing phytoplankton bloom phenology in upwelling-influenced regions using ocean color remote sensing. Remote Sens. 13, 675 (2021).
Acknowledgements
A. Ferreira received a Ph.D. grant (SFRH/BD/144586/2019; https://doi.org/10.54499/SFRH/BD/144586/2019) from Fundação para a Ciência e a Tecnologia (FCT). A.C. Brito (CEECIND/00095/2017) and C.V. Guerreiro (CEECIND/00752/2018/C, doi: 10.54499/CEECIND/00752/2018/CP1534/CT0011) were supported by the Scientific Employment Stimulus Programme, also funded by FCT. C.V. Guerreiro also benefited from a Marie Sklodowska-Curie European Fellowship supported by the EU H2020-MSCA-IF-2017 (grant no. 796802) linked to project DUSTCO. C.R.B. Mendes (PQ 312569/2021-1) and E.R. Secchi (PQ 310597/2018-8) were both granted with a researcher fellowship from National Council for Research and Development (CNPq). A MSc fellowship from CAPES was granted to R.R. Costa. This work was also funded by the European Unionâs Horizon 2020 Research and Innovation Programme under grant agreement N 810139: Project Portugal Twinning for Innovation and Excellence in Marine Science and Earth ObservationâPORTWIMS. This is a multidisciplinary study as part of the Brazilian High Latitude Oceanography Group (GOAL) activities in the Brazilian Antarctic Programme (PROANTAR). Financial support was also provided by CNPq and Coordination for the Improvement of Higher Education Personnel (CAPES). This study was conducted within the activities of the PROVOCCAR and ECOPELAGOS projects (CNPq grant nos. 442628/2018-8 and 442637/2018-7, respectively), and is within the scope of two Projects of the Institutional Internationalization Programme (CAPES PrInt-FURG â Call no. 41/2017). CAPES also provided free access to many relevant journals through the portal âPeriódicos CAPESâ. This study had the support of FCT through the strategic projects UIDP/04292/2020 (https://doi.org/10.54499/UIDP/04292/2020) and UIDB/04292/2020 (10.54499/UIDB/04292/2020) awarded to MARE and through project LA/P/0069/2020 (https://doi.org/10.54499/LA/P/0069/2020) granted to the Associate Laboratory ARNET. It also received support from FCT through the Portuguese Polar Programme (PROPOLAR) calls in 2018/2019, 2019/2020, and 2020/2021. This research is framed within the College on Polar and Extreme Environments (Polar2E) of the University of Lisbon. Authors are also indebted to the ESA Ocean ColourâClimate Change Initiative project (https://esa-oceancolour-cci.org/), the Copernicus Marine Environment Monitoring Service (https://marine.copernicus.eu/), the GlobColour project (https://www.globcolour.info/), the NOAA Centre for Weather and Climate Prediction Centre (https://www.cpc.ncep.noaa.gov/), the NOAA Physical Sciences Laboratory (https://psl.noaa.gov/) and the Palmer Station Antarctica Long-Term Ecological Research (https://pallter.marine.rutgers.edu/) for providing free access to the datasets used in this work.
Author information
Authors and Affiliations
Contributions
A.F., C.R.B.M., A.C.B., and V.B. conceptualized the study. C.R.B.M., A.C.B., and V.B supervised the project. A.F., C.R.B.M., A.C.B., R.R.C designed the methodology. A.F., C.R.B.M, R.R.C., V.M.T., and C.V.G. were responsible for the in-situ data collected within this work. A.F. performed the remote sensing and statistical analyses. A.F. wrote the original draft. C.R.B.M, R.R.C., V.B., V.M.T., C.V.G., E.S., and A.C.B. revised and edited the manuscript. C.R.B.M., A.C.B., V.B, and E.S. were responsible for funding acquisition. All authors have read and agreed to the published version of the manuscript.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Peer review
Peer review information
Nature Communications thanks Sebastien Moreau and the other, anonymous, reviewer(s) for their contribution to the peer review of this work. A peer review file is available.
Additional information
Publisherâs note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
Source data
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the articleâs Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the articleâs Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Ferreira, A., Mendes, C.R.B., Costa, R.R. et al. Climate change is associated with higher phytoplankton biomass and longer blooms in the West Antarctic Peninsula. Nat Commun 15, 6536 (2024). https://doi.org/10.1038/s41467-024-50381-2
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41467-024-50381-2