Raw and processed acoustic data were collected in East Antarctica from the RSV Aurora Australis during two surveys: the Krill Availability, Community Trophodynamics and AMISOR Surveys (KACTAS) and the Krill Acoustics and Oceanography Survey (KAOS) in the East Antarctic. Seabed alias example file used to make Figure 3. 38_false_bottom.sv.csv (file size: 10.7 MB) This file is an example of seabed aliasing (false bottom echo) that occurred during the KAOS survey. The data in this file are acoustic (mean volume backscattering strength, Sv) sample-by-sample and in a CSV format. Ping_index - ping number Distance_GPS - along track distance from the vessel’s GPS (nautical miles). Distance_vl - along track distance from the vessel’s log (nautical miles). Not used here, so is populated by Echoview’s ‘don’t care’ value (-9.90E+37) to keep the file format consistent. Ping_date - format yyyy-mm-dd Ping_time - format hh:mm:ss Ping_milliseconds - format (integer; ms) Latitude - position from the vessel’s GPS (degrees) Longitude position from the vessel’s GPS (degrees) Depth_start - start depth of vessel echosounder logging range (m) Depth_stop - stop depth of vessel echosounder logging range (m) Range_start - start range of vessel echosounder logging range (m) Range_stop - stop range of vessel echosounder logging range (m) Sample_count number of samples in a ping. Acoustic Sv samples follow in column-wise vector (dB re 1 m-1) Transect metadata all_transects.csv (file size: 14.8 K) This is the transect metadata for both the KACTAS and KAOS surveys: Transect - transect number startDate - start date of transect dd/mm/yyyy startTime - start time of transect hh:mm endDate - end date of transect dd/mm/yyyy endTime - end time of transect dd/mm/yyyy Ping_subset - a ping subset specified between two timestamps to isolate acoustic data that occurred on transect , i.e start timestamp to stop timestamp (yyyy-mm-dd hh:mm = yyyy-mm-dd hh:mm). Direction - Direction traveled along the transect (N - north or S- south). Light - day or night when transect was observed Survey - Either KACTAS or KAOS Leg - Krill box 1 or 2 (there were two surveys only during the KAOS voyage). Pass - Sampling bout for a transect in a given direction. Example R-code FigureAndDataprocessingExample.R (12KB) This R-code provides examples of scripting acoustic data processing using EchoviewR, specifically, using ping-subsets to isolate acoustic data along a transect of interest, detecting schools and exporting echo integrations for 38, 120 and 200 kHz. The R-code for making the figures in the paper is also given. GPS vessel positions for both the KACTAS and KAOS surveys KACTASandKAOS_GPS.csv (3 MB) This CSV file gives the vessel track for both the KACTAS and KAOS surveys. GPS_date - dd/mm/yyyy GPS_time - HH:MM:SS GPS_milliseconds - integer Latitude - position from the vessel’s GPS (degrees) Longitude - position from the vessel’s GPS (degrees) x - relative grid position (x) used for plotting Figure 1 y - relative grid position (y) used for plotting Figure 1 Survey - either: KACTAS_Krillbox, KAOS_Krillbox1, or KACTAS_Krillbox2 Krill swarms for both the KACTAS and KAOS surveys KACTASandKAOSswarms.csv (2.1 MB) Krill swarms descriptors for the KACTAS and KAOS surveys in CSV format (see Table 2 for description of the data fields). Echoview file for the KACTAS survey KACTAS-survey.EV (21.8 MB) An Echoview file (version 12.0) for the KACTAS acoustic data analysis KACTAS EK60 scientific echosounder calibration values KACTAS_EK500_calibration.ecs 3,623 16/05/2022 21:21 -a-- An Echoview format calibration file for the KACTAS survey (see Table 3 for calibration values and Demer et al.20 for a description of the calibration parameters). Echoview file for the KAOS survey KAOS-survey.EV (16.4MB) An Echoview file (version 12.0) for the KAOS acoustic data analysis KAOS EK60 scientific echosounder calibration values KAOS_EK60_calibration.ecs 5,711 01/11/2021 04:36 -a– An Echoview format calibration

This is data describing acoustically observed krill swarms that was used in the Bestley et al. (2017) paper 'Predicting krill swarm characteristics important for marine predators foraging off East Antarctica' (http://onlinelibrary.wiley.com/doi/10.1111/ecog.03080/full). Abstract of the paper presented here: Open ocean predator-prey interactions are often difficult to interpret because of a lack of information on prey fields at scales relevant to predator behaviour. Hence, there is strong interest in identifying the biological and physical factors influencing the distribution and abundance of prey species, which may be of broad predictive use for conservation planning and evaluating effects of environmental change. This study focuses on a key Southern Ocean prey species, Antarctic krill Euphausia superba, using acoustic observations of individual swarms (aggregations) from a large-scale survey off East Antarctica. We developed two sets of statistical models describing swarm characteristics, one set using underway survey data for the explanatory variables, and the other using their satellite remotely sensed analogues. While survey data are in situ and contemporaneous with the swarm data, remotely sensed data are all that is available for prediction and inference about prey distribution in other areas or at other times. The fitted models showed that the primary biophysical influences on krill swarm characteristics included daylight (solar elevation/radiation) and proximity to the Antarctic continental slope, but there were also complex relationships with current velocities and gradients. Overall model performance was similar regardless of whether underway or remotely sensed predictors were used. We applied the latter models to generate regional-scale spatial predictions using a 10-yr remotely-sensed time series. This retrospective modelling identified areas off east Antarctica where relatively dense krill swarms were consistently predicted during austral mid-summers, which may underpin key foraging areas for marine predators. Spatiotemporal predictions along Antarctic predator satellite tracks, from independent studies, illustrate the potential for uptake into further quantitative modelling of predator movements and foraging. The approach is widely applicable to other krill-dependent ecosystems, and our findings are relevant to similar efforts examining biophysical linkages elsewhere in the Southern Ocean and beyond. This comma separated variable (CSV) file contains the krill swarm data used in: Bestley, S., Raymond, B., Gales, N.J., Harcourt, R.G., Hindell, M.A., Jonsen, I.D., Nicol, S., Peron, C., Sumner, M.D., Weimerskirch, H. and Wotherspoon, S.J., Cox, M.J. (2017). Predicting krill swarm characteristics important for marine predators foraging off East Antarctica. Ecography. The column descriptions are: Depth_mean_m = (units m) mean depth of a krill swarm Date = (YYYYMMDD) observation date (UTC) Time = (HH:mm:ss.ss) observation time (UTC) Lat = (dd.ddddd) latitude Lon = (ddd.ddddd) longitude transect = BROKE West transect number 7 to 11 (see Fig. 1, Bestley et al. 2017) denVolgm3 = (units g wet mass m-3) internal krill swarm density in gram wet mass per cubic metre.

This dataset contains hydroacoustic results from the Antarctic Division Biomass Experiment II (ADBEX II) cruise of the Nella Dan. This cruise is the third in a series of six cruises, performing a long term survey of krill and other zooplankton distribution and abundance. Australia was to have participated in the Second International Biomass Experiment I (SIBEX I), but withdrew due to resupply problems. ADBEX II is a reduced sampling program of what was to have been sampled during SIBEX I. Three transects were made off Antarctica in the Mawson region of the Australian sector, in January to March 1984, covering a survey area of 70,000 square kilometers. Quantitative and geographic krill distribution, abundance, mean and variance of the krill weight density, and total krill biomass were obtained. Biomass estimates for ADBEX II are given as 3.5 million tonnes, obtained by extrapolating over the survey area used on the SIBEX II cruise (1.28x10^6 square kilometers). Temperature, nutrient and salinty data were also obtained, as well as trawl results. Summary results are listed in the documentation. The fields in this dataset are: pressure temperature salinity volume geopotential samples deviation

These data represent the results of the first study to use Earth System Model (ESM) outputs of SST and chlorophyll-a to simulate circumpolar krill growth potential for the recent past (1960-1989) and future climate change scenarios (2070-2099). Growth potential is obtained using an empirically-derived krill growth model (Atkinson et al. 2006, Limnol. Oceanogr.), where growth is modeled as a function of SST and chlorophyll-a. It serves as an approximation of habitat quality, as areas that support high growth rates are assumed to be good habitat (see Murphy et al., 2017, Sci Rep). To increase confidence in the future projections, ESMs were selected and weighted for each season based on their skill at reproducing observation-based krill growth potential for the recent past. First, eleven ESMs which provided SST and chlorophyll-a outputs were obtained from the Coupled Model Inter-comparison Project 5 archive. These included: CanESM2, CMCC-CESM, CNRM-CM5, GFL-ESM2G, GFDL-ESM2M, GISS-E2-H-CC, HadGEM2-CC, IPSL-CM5A-LR, MPI-ESM-MR, MRI-ESM1 and NorESM1-ME. For each ESM, seasonal surface averages of SST and chlorophyll-a were used to calculate growth potential for the historical scenario (1960-1989), which was then bilinearly interpolated on to the same 1°x1° grid. Satellite observation-based datasets for SST and chlorophyll-a were used to calculate observation-based growth potential for the recent past (1997-2010). These comprised seasonal surface averages of SST (from the OISST v2 daily dataset, 1/4⁰ horizontal resolution) and chlorophyll-a (the mean of the SeaWiFS and Johnson et al. (2013) corrected estimate of SeaWiFS daily datasets, 1/12⁰ horizontal resolution). Observation-based growth potential was then bilinearly interpolated onto the same grid as the ESMs. ESM skill for each season was subsequently assessed against observation-based growth potential using a Taylor Diagram. The ESMs were selected and weighted according to their performance to produce a weighted subset (see "ESM_weighting_method.pdf" file). Of the netcdfs provided, "hist_mean_ensemble.nc" represents the unweighted mean of seasonal growth potential, calculated from the initial ensemble of eleven ESMs for the historical scenario. The "hist_mean_subset.nc" file represents the analogous output of the weighted subset. Future projections of seasonal growth potential for Representative Concentration Pathways (RCPs) 4.5 and 8.5 were obtained using the weighted subset for the period of 2070-2099. These projected seasonal surface averages are provided in the "rcp45_mean_subset.nc" and "rcp85_mean_subset.nc" files. RCPs represent standard climate change scenarios developed by the Intergovernmental Panel on Climate Change, with 4.5 reflecting some mitigation of carbon emissions, and 8.5 being the "business as usual" scenario. Analogous netcdfs for the weighted subset outputs of chlorophyll-a (chl) and SST (tos) for the historical and RCP scenarios are also provided in the "chl_tos_netcdfs.zip" file so that the driving environmental variables underlying growth potential can be examined.

Echosounder data were collected on a multidisciplinary research voyage conducted from the RV Tangaroa, operated by New Zealand’s National Institute of Water and Atmospheric Research Limited (NIWA). The voyage lasted 42 days, departing from Wellington, New Zealand on January 29th , 2015 and returning to the same port on 11th March 2015. Active acoustic data were obtained continuously using a calibrated scientific echosounder (Simrad EK60, Horten, Norway). The echosounder operated at 38 and 120 kHz for the duration of the voyage with a pulse duration of 1.024 ms, a pulse repetition rate of one ping per second and a 7° beam width. The echosounder data here are a subset of that collected throughout the voyage and include only data from south of 65°S. This subset of data focuses on research questions pertaining to Antarctic blue whales and krill.

This dataset contains hydroacoustic results from the First International Biomass Experiment (FIBEX) cruise of the Nella Dan. This cruise is the first in a series of six cruises, performing a long term survey of krill and other zooplankton distribution and abundance. FIBEX was an international project, done in co-operation with Japan, South Africa and France. Fourteen transects were made by Australia off Antarctica in the Australian sector (Mawson to Davis region) and Prydz Bay in January to March 1981, covering a survey area of 1.09x10^6 square kilometers. Quantitative and geographic krill distribution, abundance, mean and variance of the krill weight density, and total krill biomass were obtained. The standing stock of krill for the Prydz Bay region was calculated at the second workshop to be approximately 1.3 million tonnes, or 1.6 million tonnes if extrapolated over the survey area of the SIBEX II cruise. Temperature, nutrient and salinty data were also obtained, as well as trawl results. Summary results are listed in the documentation. A scanned copy of the acoustics log book from the voyage is available for download from the provided URL.

This data were collected on the sixth Nella Dan voyage of a long term field survey project being conducted by the Australian Antarctic Division, to collect distribution, abundance and population structure data for the krill Euphausia superba in the Prydz Bay region, Antarctica. This voyage, the Australian Antarctic Marine Biological Ecosystem Research 1986/1987 (AAMBER 86/87) cruise, operated between February and April 1987. During March, a survey of the krill population and zooplankton community was conducted, to determine the late summer distribution and abundance of krill, especially the larvae. This was done as a follow up to SIBEX II in mid-summer (Janurary) 1985. The major species investigated were Euphausia superba, Euphausia frigidia, Euphausia crystallorophias and Thysanoessa marcuria. Phytoplankton pigment analysis was also conducted at each CTD station site.

Regular Trawl At each regular trawl station a quantitative standard double oblique tow was conducted from the surface down to 200 m (or to within 10 m of the bottom at stations shallower than 200 m). Such a depth range is considered to be the best compromise between the time available for sampling and the likely vertical depth range of krill. During the hauls, ship speed was maintained at a constant 2.5 plus or minus 0.5 knots. Wire speed of 0.7 to 0.8 m/s during paying out and of 0.3 m/sec during hauling (approx. 0.5 m/s and 0.2 m/s respectively at vertical depth change rate). The net mouth angle is remarkably constant during hauling within the speed ranges given above. When the net reaches maximum depth, the winch was stopped for about 30 seconds to allow the net to stabilise before starting retrieval. When hauling, propeller thrust was turned off when the net reached a depth of 15 to 20 m; this was to minimise the effects of the propeller action on the net operation and avoids damage of the samples. Target Trawl Whenever interesting targets were seen on the echo-sounder, or large amounts of krill were required for any purpose, target trawls were performed. Once the position of the target was marked, the ship was turned and navigated to run over the target from direction required within navigation capacity. The ship speed was lowered down to below 2.0 knots before hitting the target, so that the net could be lowered down to the desired depth whenever the net reached the target. Fine adjustments were made throughout the trawl by monitoring the echo-sounder in the aft control room. For live krill target trawl, ship speed was kept as slow as possible to avoid any damage to krill. Sample processing for all regular trawl stations: RMT-8 1.Measure the total sample volume (Drain water, then measure using water replacement; mandatory only for the regular hauls) 2.Sort out all Antarctic krill and count their number. If the sample mainly consists of krill and the volume is more than ~1L, a known portion of the whole sample was sub-sampled for the further processing. 3.Stage (TL, Carapace Length, Maturity) of all krill (or subsample), up to 50 to 150 individuals, and digestive gland size (the longest axis) of up to 50 individuals were measured using digital calipers. 4.Other zooplankton groups were immediately sorted out from the catch and their numbers were recorded. Preservation of RMT-8 samples Krill (including those used for onboard demography measurements) were fixed in 10% formalin for their further analysis. Whenever excess amount of krill catch were made, they were sampled and frozen for POP (persistent organic pollutant) measurements, preserved in 80% ethanol for genetic analysis, and frozen under -80C/ liquid nitrogen for chemical analysis. Fish were preserved in formalin, EtOH, or frozen. Squids were preserved in ethanol. RMT-1 1.The whole sample was fixed with 10 % formalin. 2.If the sample volume was too large, then a known proportion of catch was randomly sub-sampled and fixed. This work was completed as part of ASAC projects 2655 and 2679 (ASAC_2655, ASAC_2679).

Taken from the "Supporting Information" for the main paper. See the referenced papers for more information. Our results are based on numerical simulation of Southern Ocean sea ice, conducted using the Los Alamos numerical sea-ice model CICE version 4.0 [CICE4; Bailey et al., 2010] configured in stand-alone mode on a 0.25 degree x 0.25 degree grid, extending to 45 degrees S, with 3-hourly output [Stevens, 2013]. The atmospheric forcing for CICE4 came from the hemispheric forecasting model Polar Limited Area Prediction Systems [Polar- LAPS; Adams, 2006] and ocean forcing from the global ocean general circulation model Australian Climate Ocean Model [AusCOM; Bi and Marsland, 2010]. The model is well-constrained in its representation of processes of sea ice formation and melt, and comparison with observed areal ice extent shows minimal deviations over the 1998-2003 period, particularly during winter [Stevens 2013]. Stevens [2013] evaluates the sensitivity of the model to the number of ice thickness categories. Sea ice thickness sensitivities in the CICE model are considered in detail in Hunke [2010, 2014]. For the warm climate scenario, changes were implemented that are consistent with the A1B scenario from the Fourth Assessment from the IPCC [Meehl et al., 2007]. This is a mid-range scenario that assumes rapid economic growth before introduction of new and more efficient technologies mid century. Specifically, the following changes were applied uniformly to the current climate forcing field for a single year: a 2 degrees C increase in air temperature, a 0.2 mm/day increase in rain, a 1.5% increase in cloud fraction, a -2.3 hPa change in surface air pressure, a 25% increase in wind, a 12 Wm-2 increase in long wave downward radiation and a 20% increase in humidity. Outputs and forcings from CICE4 that are relevant for consideration of under-ice habitats for larval krill include: snow depth, ice thickness, ice concentration, movement, ridging rate, day length (dependent on day-of-year and latitude), radiation above the ice (influenced by cloud cover), and radiation below the ice (influenced by ice and snow depth). Table 1 in the main text describes how these were used in the following two filters and one overlay for evaluating the location and suitability of potential larval krill habitat during winter. Taken from the abstract of the main paper: Over-wintering of larvae underneath Antarctic pack ice is a critical stage in the life cycle of Antarctic krill. However, there are no circumpolar assessments of available habitat for larval krill, making it difficult to evaluate how climate change may impact this life stage. We use outputs from a circumpolar sea-ice model, together with a set of simple assumptions regarding key habitat features, to identify possible regions of larval krill habitat around Antarctica during winter. In particular we assume that the location and suitability of habitat is determined by both food availability and three dimensional complexity of the sea ice. We then compare the combined area of these regions under current conditions to that under a warm climate scenario. Results indicate that, while total areal sea-ice extent decreases, there is a consistently larger area of potential larval krill habitat under warm conditions. These findings highlight that decreases in sea-ice extent may not necessarily be detrimental for krill populations and underline the complexity of predicting future trajectories for this key species in the Antarctic ecosystem.

---- Public Summary from Project ---- This project is designed to provide an understanding of the interactions between krill, other zooplankton, the physical environment and the predators dependent on krill. This will directly address a number of pressing problems facing CCAMLR (the Commission for the Conservation of Antarctic Marine Living Resources) in its attempts to manage the krill fishery using an 'ecosystem approach'. Expected outcomes: As a result of logistic operations (i.e. diversion to Casey) the 29 days on site allocated to this work was reduced to 10 days. Hence only a fraction of the intended program of work was conducted. Acoustics: Acoustics data (for 38, 120, 200kHz) was collected for the top 250m of the water column for nine and a half of the planned 13 transects in our 60 x 60 nautical mile survey region.