Search

Enhanced Antarctic ice sheet mass loss yields ocean surface freshening, cooling and sea ice expansion, which result in changes in the atmospheric conditions. Using the Southern Ocean Freshwater Input from Antarctica (SOFIA) multi-model ensemble, we study the atmospheric response to a 100-year idealized freshwater release of 0.1 Sv. All models simulate a surface-intensified tropospheric cooling and lower-stratospheric warming south of 35°S. Tropospheric cooling is attributed to sea ice expansion and the associated albedo enhancement in winter and a colder sea surface in summer. This cooling yields a downward displacement of the tropopause, reduced stratospheric water vapor content and ultimately warming around 200 hPa. An enhanced southward eddy heat flux explains warming at 10?100 hPa during austral winter. Despite a temporally (and spatially) uniform prescribed freshwater flux, a prominent sea ice seasonal cycle and atmosphere dynamics result in a distinct seasonal pattern in the occurrence and magnitude of the temperature responses.

The morphology and molecular study of the penguin brain are crucial to define its survival in the extreme conditions of Antarctica. The present study focusses on extracting different optical parameters of the penguin brain using label-free optical imaging and spectroscopic techniques. In label-free optical imaging, we have used quantitative phase imaging, which provides morphological information about the neurons in brain tissue, giving the quantitative phase value of 5 to 20 radians corresponding to the 8 µm tissue section. In label-free spectroscopic techniques, we have used autofluorescence and Raman spectroscopy. Autofluorescence spectroscopy provides molecular information about nicotinamide dinucleotide, flavins, lipofuscins, and porphyrins in the brain’s spectral range of 420 nm to 700 nm. Raman spectroscopy provides multiple peaks associated with different molecules in the brain; among them, few signals are observed at approximately 1305 cm−1, 1448 cm−1, and 1661 cm−1, which correspond to vibrational modes indicative of vibrational features within lipids and protein structures, as well as the presence of amide groups within brain tissue constituents. All these techniques provide the microscopic and molecular fingerprint of the penguin brain, which can be useful for understanding penguin’s anatomical, physiological, and social behavior.

Active subglacial lakes beneath the Antarctic Ice Sheet provide insights into the dynamic subglacial environment, with implications for ice-sheet dynamics and mass balance. Most previously identified lakes have been found upstream (>100 km) of fast-flowing glaciers in West Antarctica, and none have been found in the coastal region of Dronning Maud Land (DML) in East Antarctica. The regional distribution and extent of lakes as well as their timescales and mechanisms of filling–draining activity remain poorly understood. We present local ice surface elevation changes in the coastal DML region that we interpret as unique evidence of seven active subglacial lakes located under slowly moving ice near the grounding line margin. Laser altimetry data from ICESat-2 and ICESat (Ice, Cloud, and Land Elevation Satellites) combined with multi-temporal Reference Digital Elevation Model of Antarctica (REMA) strips reveal that these lakes actively fill and drain over periods of several years. Stochastic analyses of subglacial water routing together with visible surface lineations on ice shelves indicate that these lakes discharge meltwater across the grounding line. Two lakes are within 15 km of the grounding line, while another three are within 54 km. Ice flows 17–172 m a−1 near these lakes, much slower than the mean ice flow speed near other active lakes within 100 km of the grounding line (303 m a−1). Our results improve knowledge of subglacial meltwater dynamics and evolution in this region of East Antarctica and provide new observational data to refine subglacial hydrological models.

Euphausia superba is a well-known Antarctic crustacean of great economic and ecological importance, whose management requires accurate and precise abundance and distribution estimates. Such estimates are difficult to achieve given the remoteness, extension, and large spatio-temporal variability of its geographic distribution. Acoustic data collected on board krill fishing vessels during normal fishing operation has a great potential to enhance such abundance and distribution estimates. In the present work we test the hypothesis that design-free hydroacoustic data collected during regular fishing operations can be used to produce abundance and distribution estimates with similar accuracy and precision than design-based scientific surveys. Thus, we produced and compared distribution and abundance estimates produced using either design-free hydroacoustic data collected during regular fishing operations or design-based data from scientific surveys conducted off the South Orkney Islands during summer 2017 and 2019. Following a Bayesian geostatistical approach that considered and fitted simultaneously the spatial and temporal correlation of the data, we tested different auto-correlation structures and selected the most informative models. The comparison included the means and coefficients of variation (CV) of the probability of presence (p), conditional density (d) and relative abundance index (RAI) estimates. In addition, we also simulated scenarios of parallel and orthogonal transects and obtained RAI estimates from each scenario to compare with design-based and design-free estimates for each year. In 2017, the mean RAI estimated using design-free data (94 421 m2; CV: 14 %) was ∼ 50 % higher than the one estimated with design-based data (60 232 m2; CV: 42 %), both within the fishing area. In 2019, the mean RAI estimated using design-free data (509 413 m2 CV: 6 %) was ∼ 5-fold higher than the one obtained using design-based data (113 654 m2; CV: 33 %) in the same area. Design-free RAI estimates were highly sensitive to extrapolating the inference area from fishing to the high-density sub-area. On the other hand, changing from an hourly-resolved spatio-temporal model to a purely spatial model resulted in neglectable changes. Despite observed differences in mean estimates, both methods identified similar areas of high presence and density of Antarctic krill north and north-west of the South Orkney Islands. The 2017 estimate from design-free data was probably affected by a larger dispersion of krill, and a less observed effective area during regular fishing operations. Our results show that despite using state-of-the-art methods for processing and analyzing design-free, acoustic data collected by the fishing fleet, it still yielded unreliable RAI estimates. The bias and uncertainty related to design-free data were reduced when parallel or orthogonal transects were applied although orthogonal transects yielded results with increased accuracy as they were only 21 % lower and 0.02 % higher than the true value in 2017 and 2019, respectively. Other possible approach to minimize bias would be integrating hydroacoustic information from multiple vessels.

Plastic particles are present in biotic and abiotic matrices; hence, plastic pollution is a global issue involving terrestrial and marine fauna and poses a threat to humans. Ocean circulation is a crucial vector of microplastics worldwide. Plastic pollution is among the significant threats to the ocean ecosystem. Studies and papers on plastic pollution in the oceans worldwide have been reported. However, the distribution, characterization, and abundance of micro- and nano plastics in the global ocean still need to be carefully investigated. Once plastics are present in the environment, they denature, degrade, and are more prone to fragmentation. It is well established that large plastic objects and macroplastics fragment into mesoplastics and large microplastics through photodegradation and weathering. Hence microplastics easily break up into fragments <100 µm (small microplastics, SMPs) or even into sub-micrometric particles, the nanoplastics. The small size of these SMPs and nanoplastics allows them to be ingested by different organisms according to their mouthparts’ size. Besides, this fragmentation will enable additives and plasticizers to be released into the environment, where they may pose a threat to biota throughout the trophic web in various ecosystems, e.g., from oceans and soils to glaciers. Micro- and nanoplastics (MNPs) can be transported over long distances, together with the other airborne particles. As a result of long-range transport and short-range transport, airborne MNPs can be carried from worldwide to mountain glaciers; from mid-latitudes, they can reach the very high and very low latitudes, i.e., the Arctic and Antarctica. Due to global climate change, warm ocean streams heavily affect the sea circulation in polar areas, carrying regulated and emerging pollutants, microplastics being among them. In this scenario, polar environments may be significantly enriched by MNPs carried by warmer ocean currents intruding into the polar oceans and those in atmospheric aerosol. MNPs may threaten the sea ice formation and enhance the melting of glaciers. The melting and disappearance of glaciers and the intrusion of warm currents into polar areas are also compounded by the thawing of permafrost, which can release pollutants, including MNPs. This Research Topic aims to study the interconnected pathways of MNPs that are paramount to understanding the global microplastic cycle and how climate change alters polar environments and the rest of the world. Furthermore, we aim to identify bioindicators in marine species, populations, and ecosystems, while acknowledging the interconnectedness of freshwater, terrestrial and atmospheric environments to the polar environment. Research on world glaciers will provide a comprehensive evaluation of the impacts of plastic pollution on the marine polar environment and biota, including impacts on humans.

Abstract Antarctic sea ice is one of the largest biomes on Earth providing a critical habitat for ice algae. Measurements of primary production in Antarctic sea ice remain scarce and an observation-based estimate of primary production has not been revisited in over 30 years. We fill this knowledge gap by presenting a newly compiled circumpolar data set of particulate and dissolved organic carbon from 362 ice cores, sampled between 1989 and 2019, to estimate sea-ice net community production using a carbon biomass accumulation approach. Our estimate of 26.8?32.9 Tg C yr?1 accounts for at least 15%?18% of the total primary production in the Antarctic sea-ice zone, less than a previous observation-based estimate (63?70 Tg C yr?1) and consistent with recent modeled estimates. The results underpin the ecological significance of sea-ice algae as an early season resource for pelagic food webs.

Observations of water stable isotopes in Antarctic surface snow, precipitation and water vapor are key for improving our understanding of the atmospheric water cycle and past climate reconstructions from ice cores. In this study, we use isotopic observations in Antarctica to assess the skill of the isotope-enabled atmospheric general circulation model LMDZ6, nudged to ERA5 above the boundary layer (1980?2023 period). The model has no significant bias for time-mean temperature and snow accumulation over the ice sheet. Sensitivity test on parameterized supersaturation strength highlights its opposite effect on precipitation ${\delta }^{18}$O and d-excess. Selecting an intermediate supersaturation strength resulted in a minimal bias for surface snow ${\delta }^{18}$O across the continent, with a reduced but systematic positive bias in surface snow d-excess ( ${\sim} $5?). We then assessed seasonal and diurnal isotope variability with daily precipitation and continuous vapor isotopes at Dumont d?Urville (DDU, coastal station) and Concordia (inland station). On a seasonal scale, LMDZ6iso accurately reproduces the seasonal cycle of precipitation ${\delta }^{18}$O and d-excess at both stations. Moving from statistical evaluation to physical analysis, we use the individual process contributions to boundary-layer water vapor isotopes to identify the main drivers controlling the clear-sky isotopic daily cycles. At Concordia, daily isotope variations are mainly driven by surface sublimation, whereas at DDU they are driven by surface sublimation and advection by the katabatic flow. Our results suggest that to further improve water isotopes in LMDZ6iso, fractionation during surface sublimation should be included and fractionation at condensation for low temperature should be better constrained.

The management strategy for the Antarctic krill (Euphausia superba) fishery is being revised. A key aim is to spatially and temporally allocate catches in a manner that minimizes impacts to both the krill stock and dependent predators. This process requires spatial information on the distribution and abundance of krill, yet gaps exist for an important fishing area surrounding the South Orkney Islands in the south Scotia Sea. To fill this need, we create a dynamic distribution model for krill in this region. We used data from a spatially and temporally consistent acoustic survey (2011-2020) and year-specific environmental covariates within a two-part hurdle model. The model successfully captured observed spatial and temporal patterns in krill density. The covariates found to be most important included distance from shelf break, distance from summer sea ice extent, and salinity. The northern and eastern shelf edges of the South Orkney Islands were areas of consistently high krill density and displayed strong spatial overlap between intense fishing activity and foraging chinstrap penguins. High mean krill density was also linked to oceanographic features located within the Weddell Sea. Our data suggest that years in which these features were closer to the South Orkney shelf were also years of positive Southern Annular Mode and higher observed krill densities. Our findings highlight existing fishery?predator?prey overlap in the region and support the hypothesis that Weddell Sea oceanography may play a role in transporting krill into this region. These results will feed into the next phase of krill fisheries management assessment.

Over the last decade, the Southern Ocean has experienced episodes of severe sea ice area decline. Abrupt events of sea ice loss are challenging to predict, in part due to incomplete understanding of processes occurring at the scale of individual ice floes. Here, we use high-resolution altimetry (ICESat-2) to quantify the seasonal life cycle of floes in the perennial sea ice pack of the Weddell Sea. The evolution of the floe chord distribution (FCD) shows an increase in the proportion of smaller floes between November and February, which coincides with the asymmetric melt–freeze cycle of the pack. The freeboard ice thickness distribution (fITD) suggests mirrored seasonality between the western and southern sections of the Weddell Sea ice cover, with an increasing proportion of thicker floes between October and March in the south and the opposite in the west. Throughout the seasonal cycle, there is a positive correlation between the mean chord length of floes and their average freeboard thickness. Composited floe profiles reveal that smaller floes are more vertically round than larger floes and that the mean roundness of floes increases during the melt season. These results show that regional differences in ice concentration and type at larger scales occur in conjunction with different behaviors at the small scale. We therefore suggest that floe-derived metrics obtained from altimetry could provide useful diagnostics for floe-aware models and improve our understanding of sea ice processes across scales.

Antarctic sea ice has changed significantly over the past four decades; yet limited understanding of fundamental processes, including its seasonal cycle, hinders our ability to interpret these changes. Here, we examine the processes determining the moment when sea ice locally disappears each spring, defined as the retreat date, using satellite observations over 1994?2020. We find that climatological retreat date is driven by sea ice melt in most of the seasonal ice zone and strongly constrained by the seasonal maximum ice thickness. Ice removal due to drifting ice export predominantly drives retreat only in coastal polynyas. At interannual timescales, retreat date anomalies are also preconditioned by prior maximum ice thickness, which affects melt-driven spring ice loss through the ice-albedo feedback, though this effect appears limited to specific regions. Winds emerge as a primary driver of interannual variability in the retreat date, influencing both drift- and melt-related spring ice removal processes.

During the Quaternary, ice sheets experienced several retreat–advance cycles, strongly influencing climate patterns. In order to properly simulate these phenomena, it is preferable to use physics-based models instead of parameterizations to estimate the surface mass balance (SMB), which strongly influences the evolution of the ice sheet. To further investigate the potential of these SMB models, this work evaluates the BErgen Snow SImulator (BESSI), a multi-layer snow model with high computational efficiency, as an alternative to providing the SMB for the Earth system model iLOVECLIM for multi-millennial simulations as in paleostudies. We compare the behaviors of BESSI and insolation temperature melt (ITM), an existing SMB scheme of iLOVECLIM during the Last Interglacial (LIG). Firstly, we validate the two SMB models using the regional climate model Mod- èle Atmosphérique Régional (MAR) as forcing and reference for the present-day climate over the Greenland and Antarctic ice sheets. The evolution of the SMB over the LIG (130–116 ka) is computed by forcing BESSI and ITM with transient climate forcing obtained from iLOVECLIM for both ice sheets. For present-day climate conditions, both BESSI and ITM exhibit good performance compared to MAR despite a much simpler model setup. While BESSI performs well for both Antarctica and Greenland for the same set of parame- ters, the ITM parameters need to be adapted specifically for each ice sheet. This suggests that the physics embedded in BESSI allows better capture of SMB changes across varying climate conditions, while ITM displays a much stronger sen- sitivity to its tunable parameters. The findings suggest that BESSI can provide more reliable SMB estimations for the iLOVECLIM framework to improve the model simulations of the ice sheet evolution and interactions with climate for multi-millennial simulations.

Snowfall is an important component of the mass balance of ice sheets and glaciers in Antarctica. In coastal Victoria Land (VL), changes to snowfall can impact ice masses, landscapes, and coastal ecosystems. Coastal VL is characterized by strong gradients in snowfall rates between the polar desert of the McMurdo Dry Valleys and the high accumulation in northern VL. Extreme precipitation events significantly contribute to total precipitation, with the largest contribution in the Terra Nova Bay area. We present a comprehensive analysis of snowfall dynamics in this region, using a Lagrangian moisture source diagnostic to study moisture sources and Self-Organizing Maps (SOM) to link these to different synoptic weather types. The moisture for snowfall in VL originates from the Southern Ocean, with more local sources in the Ross Sea embayment in summer when sea ice is reduced. We show a strong division in moisture sources between northern and southern VL, with the north receiving precipitation from moisture sources to the west and southern VL from the east. Precipitation in northern VL results from meridional transport of marine air from lower latitudes, while precipitation in southern VL is related to cyclonic disturbances in the Ross Sea that bring moisture from the east. Extreme precipitation in northern VL occurs during blocking highs that intensify meridional transport. Such intrusions of marine air, sometimes in the form of atmospheric rivers, do not impact the more isolated western Ross Ice Shelf and southern VL further in the Ross Sea embayment.

In the subauroral zone at the boundary of the auroral oval in the evening and night hours during geomagnetic disturbances, a narrow (about 1°–2°) and extended structure (several hours in longitude) is formed. It is known as a polarization jet (PJ) or the subauroral ion drift (SAID). The PJ/SAID is a fast westward ion drift and is one of the main signatures of a geomagnetic disturbance in the subauroral ionosphere at the altitudes of the F-layer, when the geomagnetic AE index reaches more than 500 nT. Plasma speed in the PJ/SAID can reach several kilometres per second, and the size of plasma irregularities inside it can reach scales from tens of meters to several hundred meters. Such high velocities and structured plasma can affect trans-ionospheric radio waves and lead to scintillations in the received signal. We show that at the moment of auroral activity intensification, an increase in the magnitude of phase scintillation index (σϕ) as well as loss of satellite signals lock were observed in the region of the PJ/SAID equatorward of the auroral oval over Dronning Maud Land (Queen Maud Land) in Antarctica. We find that fluctuations inside the PJ/SAID can lead to serious deterioration of radio communication or navigational services. We emphasize the importance of considering the geometry of the beam passing from the GNSS satellite to the receiver on the ground. We highlight the mutual contribution of the PJ/SAID and the diffuse aurora boundary, which are almost impossible to separate in practice. Our results demonstrate the importance of considering the subauroral zone, where very dynamic plasma formations can occur with a strong flow and various-scale irregularities inside that lead to serious interference in satellite communications.

Understanding the connection between maturity stages and morphology in relation to size selectivity in trawls is essential for assessing the impact of various fishing gear on the population structures of harvested species, their fishing mortality rates, and the efficiency of the gear used. The Antarctic krill (Euphausia superba) fishery is the largest in the Southern Ocean by volume, and there is increasing interest in expanding the industry. The krill fishery employs different trawl designs and is not currently subject to technical regulations specifying the types of fishing gear and mesh sizes that can legally be used. There is a need to establish a robust model predicting size selectivity that includes the morphological variation in the population of krill. Male and female Antarctic krill are described with 12 maturity stages, from juveniles to sexually mature adults, each with distinct morphological features. The current study established a morphological description of each individual krill maturity stage to identify and parameterize what determines size selectivity using the FISHSELECT framework. This framework is used to predict size selectivity for each of the different stages in various mesh sizes and openings relevant to the krill fishery, in both actual and virtual populations. The results can be used to assess size selectivity for specific fishing gears and population structures, facilitating more accurate understanding and modeling of the fishery’s impact on the demographic composition of the krill stock.

Abstract Basal melting of ice shelves is fundamental to Antarctic ice sheet mass loss, yet direct observations remain sparse. We present the first year-round melt record (2017-2021) from a phase-sensitive radar on Fimbulisen, one of the fastest flowing ice shelves in Dronning Maud Land, East Antarctica. The observed long-term mean ablation rate at 350 m depth below the central ice shelf was 1.0 ± 0.5 m yr?1, marked by substantial sub-weekly variability ranging from 0.4 to 3.5 m yr?1. 36-h filtered basal melt rate fluctuations closely align with ocean velocity. On seasonal time scales, melt rates peak during austral spring to autumn (September-March), driven by both elevated ocean velocities and thermal driving near the base. The combined effect of thermal driving and current speed explains the majority of the melt rate variability (r = 0.84), highlighting the dominant role of shear-driven turbulence. This relationship enables parameterization of melt rates for the decade-long ocean record (2010?2021), although deviations appear under low and high forcing conditions. Both observed and parameterized melt rates show similar yearly mean magnitudes compared to satellite-derived melt rates but with a tenfold lower seasonal amplitude and a 3-month delay in seasonality. These detailed concurrent ice?ocean observations provide essential validation data for remote sensing and numerical models that aim to quantify and project ice-shelf response to a change in ocean forcing. In situ measurements and continued monitoring are crucial for accurately assessing and modeling future basal melt rates, and for understanding the complex dynamics driving ice-shelf stability and sea-level change.

The Southern Ocean (south of 30°S) contributes significantly to global ocean carbon uptake through the solubility, physical and biological pumps. Many studies have estimated carbon export to the deep ocean, but very few have attempted a basin-scale perspective, or accounted for the sea-ice zone (SIZ). In this study, we use an extensive array of BGC-Argo floats to improve previous estimates of carbon export across basins and frontal zones, specifically including the SIZ. Using a new method involving changes in particulate organic carbon and dissolved oxygen along the mesopelagic layer, we find that the total Southern Ocean carbon export from 2014 to 2022 is 2.69 ± 1.23 PgC y?1. The polar Antarctic zone contributes the most (41%) with 1.09 ± 0.46 PgC y?1. Conversely, the SIZ contributes the least (8%) with 0.21 ± 0.09 PgC y?1 and displays a strong shallow respiration in the upper 200 m. However, the SIZ contribution can increase up to 14% depending on the depth range investigated. We also consider vertical turbulent fluxes, which can be neglected at depth but are important near the surface. Our work provides a complementary approach to previous studies and is relevant for work that focuses on evaluating the biogeochemical impacts of changes in Antarctic sea-ice extent. Refining estimates of carbon export and understanding its drivers ultimately impacts our comprehension of climate variability at the global ocean scale.

The unique challenges of polar ecosystems, coupled with the necessity for high-precision data, make Unmanned Aerial Vehicles (UAVs) an ideal tool for vegetation monitoring and conservation studies in Antarctica. This review draws on existing studies on Antarctic UAV vegetation mapping, focusing on their methodologies, including surveyed locations, flight guidelines, UAV specifications, sensor technologies, data processing techniques, and the use of vegetation indices. Despite the potential of established Machine-Learning (ML) classifiers such as Random Forest, K Nearest Neighbour, and Support Vector Machine, and gradient boosting in the semantic segmentation of UAV-captured images, there is a notable scarcity of research employing Deep Learning (DL) models in these extreme environments. While initial studies suggest that DL models could match or surpass the performance of established classifiers, even on small datasets, the integration of these advanced models into real-time navigation systems on UAVs remains underexplored. This paper evaluates the feasibility of deploying UAVs equipped with adaptive path-planning and real-time semantic segmentation capabilities, which could significantly enhance the efficiency and safety of mapping missions in Antarctica. This review discusses the technological and logistical constraints observed in previous studies and proposes directions for future research to optimise autonomous drone operations in harsh polar conditions.

Ice-sheet mass loss is one of the clearest manifestations of climate change, with Antarctica discharging mass into the ocean via melting or through calving. The latter produces icebergs that can modify ocean water properties, often at great distances from source. This affects upper-ocean physics and primary productivity, with implications for atmospheric carbon drawdown. A detailed understanding of iceberg modification of ocean waters has hitherto been hindered by a lack of proximal measurements. Here unique measurements of a giant iceberg from an underwater glider enable quantification of meltwater effects on the physical and biological processes in the upper layers of the Southern Ocean, a region disproportionately important for global heat and carbon sequestration. Iceberg basal melting erodes seasonally produced winter water layer stratification, normally forming a strong potential energy barrier to vertical exchange of surface and deep waters, while freshwater run-off increases and shoals near-surface stratification. Nutrient-rich deeper waters, incorporating meltwater loaded with terrigenous material, are ventilated to below this stratification maxima, providing a potential mechanism for alleviating critical phytoplankton-limiting components. Regional historical hydrographic data demonstrate similar stratification changes during the passage of another large iceberg, suggesting that they may be an important pathway of aseasonal winter water modification.

Seabirds can disperse widely when searching for prey, particularly during nonbreeding periods. Conservation measures predominately focus on protecting breeding colonies, while spatial protection at sea is often based on knowledge of the distribution of breeding adults, despite accumulating evidence that marine habitats used by immature birds sometimes differ from those of adults. Juvenile emperor penguins from Atka Bay, west Dronning Maud Land, Antarctica, tracked immediately after fledging performed long migrations to the northern extents of the Convention for the Conservation of Antarctic Marine Living Resources subareas 48.4 and 48.6. Individuals did not remain long at their northern positions, before commencing a rapid southerly movement to within a few hundred km of the marginal ice zone (MIZ). The initial migratory movement was broadly synchronous across individuals. The southward movement and subsequent change to area-restricted searching were consistent with the MIZ representing a potentially important feeding habitat for juvenile emperor penguins. Spatio-temporal management mechanisms may be beneficial in reducing threats to these young penguins.

This study compares CL51 ceilometer observations made at Scott Base, Antarctica, with statistics from the ERA5, JRA55, and MERRA2 reanalyses. To enhance the comparison we use a lidar instrument simulator to derive cloud statistics from the reanalyses which account for instrumental factors. The cloud occurrence in the three reanalyses is slightly overestimated above 3 km, but displays a larger underestimation below 3 km relative to observations. Unlike previous studies, we see no relationship between relative humidity and cloud occurrence biases, suggesting that the cloud biases do not result from the representation of moisture. We also show that the seasonal variation of cloud occurrence and cloud fraction, defined as the vertically integrated cloud occurrence, are small in both the observations and the reanalyses. We also examine the quality of the cloud representation for a set of weather states derived from ERA5 surface winds. The variability associated with grouping cloud occurrence based on weather state is much larger than the seasonal variation, highlighting weather state is a strong control of cloud occurrence. All the reanalyses continue to display underestimates below 3 km and overestimates above 3 km for each weather state. But the variability in ERA5 statistics matches the changes in the observations better than the other reanalyses. We also use a machine learning scheme to estimate the quantity of supercooled liquid water cloud from the ceilometer observations. Ceilometer low-level supercooled liquid water cloud occurrences are considerably larger than values derived from the reanalyses, further highlighting the poor representation of low-level clouds in the reanalyses.