Search

Elongated open-water areas in sea ice (leads) release sea spray particles to the atmosphere. However, there is limited knowledge on the amount, properties and drivers of sea spray emitted from leads, and no existing parameterization of this process is available for use in models. In this work, we use measurements of aerosol fluxes from Nilsson et al. (2001) to produce an estimate of the location, timing and amount of sea spray emissions from leads at the scale of the Arctic Ocean for 1 year. Lead fractions are derived using sea ice data sets from numerical models and satellite detection. The proposed parameterization estimates that leads account for 0.3 %–9.8 % of the annual sea salt aerosol number emissions in the Arctic Ocean regions where sea ice concentration is greater than 80 %. Assuming similar size distributions to those from emissions from the open ocean, leads account for 30 %–85 % of mass emissions in sea ice regions. The total annual mass of sea salt emitted from leads, 0.1–2.1 Tg yr−1, is comparable to the mass of sea salt aerosol transported above sea ice from the open ocean, according to the MERRA-2 reanalysis. In addition to providing the first estimates of possible upper and lower bounds of sea spray emissions from leads, the conceptual model developed in this work is implemented and tested in the regional atmospheric chemistry model WRF-Chem. Given the estimates obtained in this work, the impact of sea spray from leads on Arctic clouds and radiative budget needs to be further explored.

We have implemented the brittle Bingham–Maxwell sea ice rheology (BBM) into SI3, the sea ice component of NEMO. After discussing the numerical aspects and requirements that are specific to the implementation of a brittle rheology in the Eulerian, finite-difference, Arakawa C-grid framework, we detail the approach we have used. This approach relies on the introduction of an additional set of prognostic stress tensor components, sea ice damage, and sea ice velocity vector, following a grid point arrangement that expands the C-grid into the Arakawa E-grid. The newly implemented BBM rheology is first assessed by means of a set of idealized SI3 simulations at different spatial resolutions. Then, sea ice deformation rates obtained from simulations of the Arctic at a 1/4° spatial resolution, performed with the coupled ocean–sea ice setup of NEMO, are assessed against satellite observations. For all these simulations, results obtained with the default current workhorse setup of SI3 are provided to serve as a reference. Our results show that using a brittle type of rheology, such as BBM, allows SI3 to simulate the highly localized deformation pattern of sea ice, as well as its scaling properties, from the scale of the model's computational grid up to the basin scale.

In this study, we investigate the fully multivariate state and parameter estimation through idealised simulations of a dynamics-only model that uses the novel Maxwell elasto-brittle (MEB) sea-ice rheology and in which we estimate not only the sea-ice concentration, thickness and velocity, but also its level of damage, internal stress and cohesion. Specifically, we estimate the air drag coefficient and the so-called damage parameter of the MEB model. Mimicking the realistic observation network with different combinations of observations, we demonstrate that various issues can potentially arise in a complex sea-ice model, especially in instances for which the external forcing dominates the model forecast error growth. Even though further investigation will be needed using an operational (a coupled dynamics–thermodynamics) sea-ice model, we show that, with the current observation network, it is possible to improve both the observed and the unobserved model state forecast and parameter accuracy.

Short-term sea ice predictability is challenging despite recent advancements in sea ice modelling and new observations of sea ice deformation that capture small-scale features (open leads and ridges) at the kilometre scale. A new method for assimilation of satellite-derived sea ice deformation into numerical sea ice models is presented. Ice deformation provided by the Copernicus Marine Service is computed from sea ice drift derived from synthetic aperture radar at a high spatio-temporal resolution. We show that high values of ice deformation can be interpreted as reduced ice concentration or increased ice damage – i.e. scalar variables responsible for ice strength in brittle or visco-plastic sea ice dynamical models. This method is tested as a proof of concept with the neXt-generation Sea Ice Model (neXtSIM), where the assimilation scheme uses a data insertion approach and forecasting with one member. We obtain statistics of assimilation impact over a long test period with many realisations starting from different initial times. Assimilation and forecasting experiments are run on synthetic and real observations in January 2021 and show increased accuracy of deformation prediction for the first 3–4 d. Similar conclusions are obtained using both brittle and visco-plastic rheologies implemented in neXtSIM. Thus, the forecasts improve due to the update of sea ice mechanical properties rather than the exact rheological formulation. It is demonstrated that the assimilated information can be extrapolated in space – gaps in spatially discontinuous satellite observations of deformation are filled with a realistic pattern of ice cracks, confirmed by later satellite observations. The limitations and usefulness of the proposed assimilation approach are discussed in a context of ensemble forecasts. Pathways to estimate intrinsic predictability of sea ice deformation are proposed.

Multiyear sea ice (MYI) cover in the Arctic has been monitored for decades using increasingly sophisticated remote sensing techniques, and these have documented a significant decline in MYI over time. However, such techniques are unable to differentiate between the processes affecting the evolution of the MYI. Further, estimating the thickness and thus the volume of MYI remains challenging. In this study we employ a sea ice–ocean model to investigate the changes to MYI over the period 2000–2018. We exploit the Lagrangian framework of the sea ice model to introduce a new method of tracking MYI area and volume which is based on identifying MYI during freeze onset each autumn. The model is found to successfully reproduce the spatial distribution and evolution of observed MYI extent. We discuss the balance of the processes (melt, ridging, export, and replenishment) linked to the general decline in MYI cover. The model suggests that rather than one process dominating the losses, there is an episodic imbalance between the different sources and sinks of MYI. We identify those key to the significant observed declines in 2007 and 2012; while melt and replenishment are important in 2012, sea ice dynamics play a significant role in 2007. Notably, the model suggests that in years such as 2007, convergence of the ice, through ridging, can result in large reductions in MYI area without a corresponding loss of MYI volume. This highlights the benefit of using models alongside satellite observations to aid interpretation of the observed MYI evolution in the Arctic.

Advanced data assimilation (DA) methods, widely used in geophysical and climate studies to merge observations with numerical models, can improve state estimates and consequent forecasts. We interface the deterministic ensemble Kalman filter (DEnKF) to the Lagrangian neXt generation Sea Ice Model, neXtSIM. The ensemble is generated by perturbing the atmospheric and oceanic forcing throughout the simulations and randomly initialized ice cohesion. Our ensemble–DA system assimilates sea ice concentration (SIC) from the Ocean and Sea Ice Satellite Application Facility (OSI-SAF) and sea ice thickness (SIT) from the merged CryoSat-2 and SMOS datasets (CS2SMOS). Because neXtSIM is computationally solved on a time-dependent evolving mesh, it is a challenging application for ensemble–DA. As a solution, we perform the DEnKF analysis on a fixed and regular reference mesh, on which model variables are interpolated before the DA and then back to each member's mesh after the DA. We evaluate the impact of assimilating different types of sea ice observations on the model's forecast skills of the Arctic sea ice by comparing satellite observations and a free-run ensemble in an Arctic winter period, 2019–2020. Significant improvements in modeled SIT indicate the importance of assimilating weekly CS2SMOS SIT, while the improvements of SIC and ice extent are moderate but benefit from daily ingestion of the OSI-SAF SIC. For most of the winter, the correlation between SIT and SIC is weaker, which results in little cross-inference between the two variables in the assimilation step. Overall, the ensemble–DA system based on the stand-alone sea ice model demonstrates the feasibility of winter Arctic sea ice prediction with good computational efficiency. These results open the path toward operational implementation and the extension to multi-year assimilation.

Sea ice is a key component of the Earth's climate system as it modulates the energy exchanges and associated feedback processes at the air–sea interface in polar regions. These exchanges have been suggested to strongly depend on openings in the sea ice cover, which are associated with fine-scale sea ice deformations, but the importance of these processes remains poorly understood as most numerical models struggle to represent these deformations without using very costly horizontal resolutions (≃ 5 km). In this study, we present results from a 12 km resolution ocean–sea ice coupled model, the first that uses a brittle rheology to represent the mechanical behaviour of sea ice. This rheology has been shown to reproduce observed characteristics and complexity of fine-scale sea ice deformations at relatively coarse resolutions. We evaluate and discuss the Arctic sea ice mass balance of this coupled model for the period 2000–2018. We first assess sea ice quantities relevant for climate (volume, extent, and drift) and find that they are consistent with satellite observations. We evaluate components of the mass balance for which observations are available, i.e. sea ice volume export through Fram Strait and winter mass balance in the Arctic marginal seas for the period 2003–2018. Model values show a good match with observations, remaining within the estimated uncertainty, and the interannual variability of the dynamic contribution to the winter mass balance is generally well captured. We discuss the relative contributions of dynamics and thermodynamics to the sea ice mass balance in the Arctic Basin for 2000–2018. Using the ability of the model to represent divergence motions at different scales, we investigate the role of leads and polynyas in ice production. We suggest a way to estimate the contribution of leads and polynyas to ice growth in winter, and we estimate this contribution to add up to 25 %–35 % of the total ice growth in pack ice from January to March. This contribution shows a significant increase over 2000–2018. This coupled framework opens up new opportunities to understand and quantify the interplay between small-scale sea ice dynamics and ocean properties.

This article describes the observation techniques and suggests processing methods to estimate dynamical sea-ice parameters from data of the Earth Explorer 10 candidate Harmony. The two Harmony satellites will fly in a reconfigurable formation with Sentinel-1D. Both will be equipped with a multi-angle thermal infrared sensor and a passive radar receiver, which receives the reflected Sentinel-1D signals using two antennas. During the lifetime of the mission, two different formations will be flown. In the stereo formation, the Harmony satellites will fly approximately 300 km in front and behind Sentinel-1, which allows for the estimation of instantaneous sea-ice drift vectors. We demonstrate that the addition of instantaneous sea-ice drift estimates on top of the daily integrated values from feature tracking have benefits in terms of interpretation, sampling and resolution. The wide-swath instantaneous drift observations of Harmony also help to put high-temporal-resolution instantaneous buoy observations into a spatial context. Additionally, it allows for the extraction of deformation parameters, such as shear and divergence. As a result, Harmony's data will help to improve sea-ice statistics and parametrizations to constrain sea-ice models. In the cross-track interferometry (XTI) mode, Harmony's satellites will fly in close formation with an XTI baseline to be able to estimate surface elevations. This will allow for improved estimates of sea-ice volume and also enables the retrieval of full, two-dimensional swell-wave spectra in sea-ice-covered regions without any gaps. In stereo formation, the line-of-sight diversity allows the inference of swell properties in both directions using traditional velocity bunching approaches. In XTI mode, Harmony's phase differences are only sensitive to the ground-range direction swell. To fully recover two-dimensional swell-wave spectra, a synergy between XTI height spectra and intensity spectra is required. If selected, the Harmony mission will be launched in 2028.

The neXtSIM-F (neXtSIM forecast) forecasting system consists of a stand-alone sea ice model, neXtSIM (neXt-generation Sea Ice Model), forced by the TOPAZ ocean forecast and the ECMWF atmospheric forecast, combined with daily data assimilation of sea ice concentration. It uses the novel brittle Bingham–Maxwell (BBM) sea ice rheology, making it the first forecast based on a continuum model not to use the viscous–plastic (VP) rheology. It was tested in the Arctic for the time period November 2018–June 2020 and was found to perform well, although there are some shortcomings. Despite drift not being assimilated in our system, the sea ice drift is good throughout the year, being relatively unbiased, even for longer lead times like 5 d. The RMSE in speed and the total RMSE are also good for the first 3 or so days, although they both increase steadily with lead time. The thickness distribution is relatively good, although there are some regions that experience excessive thickening with negative implications for the summertime sea ice extent, particularly in the Greenland Sea. The neXtSIM-F forecasting system assimilates OSI SAF sea ice concentration products (both SSMIS and AMSR2) by modifying the initial conditions daily and adding a compensating heat flux to prevent removed ice growing back too quickly. The assimilation greatly improves the sea ice extent for the forecast duration.

We explore several statistical properties of the observed and simulated Arctic sea-ice lead fraction, as well as the statistics of simulated Arctic ocean–atmosphere heat fluxes. First we show that the observed lead fraction in the Central Arctic has a monofractal spatial scaling, which we relate to the multifractal spatial scaling present in sea-ice deformation rates. We then show that the relevant statistics of the observed lead fraction in the Central Arctic are well represented by our model, neXtSIM. Given that the heat flux through leads may be up to 2 orders of magnitude larger than that through unbroken ice, we then explore the statistical properties (probability distribution function – PDF – and spatial scaling) of the heat fluxes simulated by neXtSIM. We demonstrate that the modelled heat fluxes present a multifractal scaling in the Central Arctic, where heat fluxes through leads dominate the high-flux tail of the PDF. This multifractal character relates to the multi- and monofractal character of deformation rates and the lead fraction. In the wider Arctic, the high-flux tail of the PDF is dominated by an exponential decay, which we attribute to the presence of coastal polynyas. Finally, we show that the scaling of the simulated lead fraction and heat fluxes depends weakly on the model resolution and discuss the role sub-grid-scale parameterisations of the ice heterogeneity may have in improving this result.

As sea ice extent decreases in the Arctic, surface ocean waves have more time and space to develop and grow, exposing the marginal ice zone (MIZ) to more frequent and more energetic wave events. Waves can fragment the ice cover over tens of kilometres, and the prospect of increasing wave activity has sparked recent interest in the interactions between wave-induced sea ice fragmentation and lateral melting. The impact of this fragmentation on sea ice dynamics, however, remains mostly unknown, although it is thought that fragmented sea ice experiences less resistance to deformation than pack ice. Here, we introduce a new coupled framework involving the spectral wave model WAVEWATCH III and the sea ice model neXtSIM, which includes a Maxwell elasto-brittle rheology. This rheological framework enables the model to efficiently track and keep a “memory” of the level of sea ice damage. We propose that the level of sea ice damage increases when wave-induced fragmentation occurs. We used this coupled modelling system to investigate the potential impact of such a local mechanism on sea ice kinematics. Focusing on the Barents Sea, we found that the internal stress decrease of sea ice resulting from its fragmentation by waves resulted in a more dynamical MIZ, particularly in areas where sea ice is compact. Sea ice drift is enhanced for both on-ice and off-ice wind conditions. Our results stress the importance of considering wave–sea-ice interactions for forecast applications. They also suggest that waves likely modulate the area of sea ice that is advected away from the pack by the ocean, potentially contributing to the observed past, current and future sea ice cover decline in the Arctic.

In this paper, we evaluate the neXtSIM sea ice model with respect to the observed scaling invariance properties of sea ice deformation in the spatial and temporal domains. Using an Arctic setup with realistic initial conditions, state-of-the-art atmospheric reanalysis forcing and geostrophic currents retrieved from satellite data, we show that the model is able to reproduce the observed properties of this scaling in both the spatial and temporal domains over a wide range of scales, as well as their multi-fractality. The variability of these properties during the winter season is also captured by the model. We also show that the simulated scaling exhibits a space–time coupling, a suggested property of brittle deformation at geophysical scales. The ability to reproduce the multi-fractality of this scaling is crucial in the context of downscaling model simulation outputs to infer sea ice variables at the sub-grid scale and also has implications for modeling the statistical properties of deformation-related quantities, such as lead fractions and heat and salt fluxes.

Numerical models solved on adaptive moving meshes have become increasingly prevalent in recent years. Motivating problems include the study of fluids in a Lagrangian frame and the presence of highly localized structures such as shock waves or interfaces. In the former case, Lagrangian solvers move the nodes of the mesh with the dynamical flow; in the latter, mesh resolution is increased in the proximity of the localized structure. Mesh adaptation can include remeshing, a procedure that adds or removes mesh nodes according to specific rules reflecting constraints in the numerical solver. In this case, the number of mesh nodes will change during the integration and, as a result, the dimension of the model's state vector will not be conserved. This work presents a novel approach to the formulation of ensemble data assimilation (DA) for models with this underlying computational structure. The challenge lies in the fact that remeshing entails a different state space dimension across members of the ensemble, thus impeding the usual computation of consistent ensemble-based statistics. Our methodology adds one forward and one backward mapping step before and after the ensemble Kalman filter (EnKF) analysis, respectively. This mapping takes all the ensemble members onto a fixed, uniform reference mesh where the EnKF analysis can be performed. We consider a high-resolution (HR) and a low-resolution (LR) fixed uniform reference mesh, whose resolutions are determined by the remeshing tolerances. This way the reference meshes embed the model numerical constraints and are also upper and lower uniform meshes bounding the resolutions of the individual ensemble meshes. Numerical experiments are carried out using 1-D prototypical models: Burgers and Kuramoto–Sivashinsky equations and both Eulerian and Lagrangian synthetic observations. While the HR strategy generally outperforms that of LR, their skill difference can be reduced substantially by an optimal tuning of the data assimilation parameters. The LR case is appealing in high dimensions because of its lower computational burden. Lagrangian observations are shown to be very effective in that fewer of them are able to keep the analysis error at a level comparable to the more numerous observers for the Eulerian case. This study is motivated by the development of suitable EnKF strategies for 2-D models of the sea ice that are numerically solved on a Lagrangian mesh with remeshing.

A new algorithm for estimating sea ice age (SIA) distribution based on the Eulerian advection scheme is presented. The advection scheme accounts for the observed divergence or convergence and freezing or melting of sea ice and predicts consequent generation or loss of new ice. The algorithm uses daily gridded sea ice drift and sea ice concentration products from the Ocean and Sea Ice Satellite Application Facility. The major advantage of the new algorithm is the ability to generate individual ice age fractions in each pixel of the output product or, in other words, to provide a frequency distribution of the ice age allowing to apply mean, median, weighted average or other statistical measures. Comparison with the National Snow and Ice Data Center SIA product revealed several improvements of the new SIA maps and time series. First, the application of the Eulerian scheme provides smooth distribution of the ice age parameters and prevents product undersampling which may occur when a Lagrangian tracking approach is used. Second, utilization of the new sea ice drift product void of artifacts from EUMETSAT OSI SAF resulted in more accurate and reliable spatial distribution of ice age fractions. Third, constraining SIA computations by the observed sea ice concentration expectedly led to considerable reduction of multi-year ice (MYI) fractions. MYI concentration is computed as a sum of all MYI fractions and compares well to the MYI products based on passive and active microwave and SAR products.

We propose a satellite mission that uses a near-nadir Ka-band Doppler radar to measure surface currents, ice drift and ocean waves at spatial scales of 40 km and more, with snapshots at least every day for latitudes 75 to 82°, and every few days for other latitudes. The use of incidence angles of 6 and 12° allows for measurement of the directional wave spectrum, which yields accurate corrections of the wave-induced bias in the current measurements. The instrument's design, an algorithm for current vector retrieval and the expected mission performance are presented here. The instrument proposed can reveal features of tropical ocean and marginal ice zone (MIZ) dynamics that are inaccessible to other measurement systems, and providing global monitoring of the ocean mesoscale that surpasses the capability of today's nadir altimeters. Measuring ocean wave properties has many applications, including examining wave–current interactions, air–sea fluxes, the transport and convergence of marine plastic debris and assessment of marine and coastal hazards.

We present a sensitivity analysis and discuss the probabilistic forecast capabilities of the novel sea ice model neXtSIM used in hindcast mode. The study pertains to the response of the model to the uncertainty on winds using probabilistic forecasts of ice trajectories. neXtSIM is a continuous Lagrangian numerical model that uses an elasto-brittle rheology to simulate the ice response to external forces. The sensitivity analysis is based on a Monte Carlo sampling of 12 members. The response of the model to the uncertainties is evaluated in terms of simulated ice drift distances from their initial positions, and from the mean position of the ensemble, over the mid-term forecast horizon of 10 days. The simulated ice drift is decomposed into advective and diffusive parts that are characterised separately both spatially and temporally and compared to what is obtained with a free-drift model, that is, when the ice rheology does not play any role in the modelled physics of the ice. The seasonal variability of the model sensitivity is presented and shows the role of the ice compactness and rheology in the ice drift response at both local and regional scales in the Arctic. Indeed, the ice drift simulated by neXtSIM in summer is close to the one obtained with the free-drift model, while the more compact and solid ice pack shows a significantly different mechanical and drift behaviour in winter. For the winter period analysed in this study, we also show that, in contrast to the free-drift model, neXtSIM reproduces the sea ice Lagrangian diffusion regimes as found from observed trajectories. The forecast capability of neXtSIM is also evaluated using a large set of real buoy's trajectories and compared to the capability of the free-drift model. We found that neXtSIM performs significantly better in simulating sea ice drift, both in terms of forecast error and as a tool to assist search and rescue operations, although the sources of uncertainties assumed for the present experiment are not sufficient for complete coverage of the observed IABP positions.

Ruxolitinib reduces spleen volume, improves symptoms, and increases survival in patients with intermediate- or high-risk myelofibrosis. However, suboptimal response may occur, potentially because of signaling via the phosphoinositide 3-kinase (PI3K)/protein kinase B pathway. This phase 2 study evaluated dosing, efficacy, and safety of add-on PI3Kδ inhibitor parsaclisib for patients with primary or secondary myelofibrosis with suboptimal response to ruxolitinib. Eligible patients remained on a stable ruxolitinib dose and received add-on parsaclisib 10 or 20 mg, once daily for 8 weeks, and once weekly thereafter (daily-to-weekly dosing; n = 32); or parsaclisib 5 or 20 mg, once daily for 8 weeks, then 5 mg once daily thereafter (all-daily dosing; n = 42). Proportion of patients achieving a ≥10% decrease in spleen volume at 12 weeks was 28% for daily-to-weekly dosing and 59.5% for all-daily dosing. Proportions of patients achieving ≥50% decrease at week 12 in Myelofibrosis Symptom Assessment Form and Myeloproliferative Neoplasms Symptom Assessment Form symptom scores were 14% and 18% for daily-to-weekly dosing, and 28% and 32% for all-daily dosing, respectively. Most common nonhematologic treatment-emergent adverse events were nausea (23%), diarrhea (22%), abdominal pain and fatigue (each 19%), and cough and dyspnea (each 18%). New-onset grade 3 and 4 thrombocytopenia were observed in 19% of patients, each dosed daily-to-weekly, and in 26% and 7% of patients dosed all-daily, respectively, managed with dose interruptions. Hemoglobin levels remained steady. The addition of parsaclisib to stable-dose ruxolitinib can reduce splenomegaly and improve symptoms, with manageable toxicity in patients with myelofibrosis with suboptimal response to ruxolitinib. This trial was registered at www.clinicaltrials.gov as #NCT02718300.

In this paper we describe a waves-in-ice model (WIM), which calculates ice breakage and the wave radiation stress (WRS). This WIM is then coupled to the new sea-ice model neXtSIM, which is based on the elasto-brittle (EB) rheology. We highlight some numerical issues involved in the coupling and investigate the impact of the WRS, and of modifying the EB rheology to lower the stiffness of the ice in the area where the ice has broken up (the marginal ice zone or MIZ). In experiments in the absence of wind, we find that wind waves can produce noticeable movement of the ice edge in loose ice (concentration around 70 %) – up to 36 km, depending on the material parameters of the ice that are used and the dynamical model used for the broken ice. The ice edge position is unaffected by the WRS if the initial concentration is higher (≳ 0.9). Swell waves (monochromatic waves with low frequency) do not affect the ice edge location (even for loose ice), as they are attenuated much less than the higher-frequency components of a wind wave spectrum, and so consequently produce a much lower WRS (by about an order of magnitude at least).In the presence of wind, we find that the wind stress dominates the WRS, which, while large near the ice edge, decays exponentially away from it. This is in contrast to the wind stress, which is applied over a much larger ice area. In this case (when wind is present) the dynamical model for the MIZ has more impact than the WRS, although that effect too is relatively modest. When the stiffness in the MIZ is lowered due to ice breakage, we find that on-ice winds produce more compression in the MIZ than in the pack, while off-ice winds can cause the MIZ to be separated from the pack ice.

Sea ice freeboard derived from satellite altimetry is the basis for the estimation of sea ice thickness using the assumption of hydrostatic equilibrium. High accuracy of altimeter measurements and freeboard retrieval procedure are, therefore, required. As of today, two approaches for estimating the freeboard using laser altimeter measurements from Ice, Cloud, and land Elevation Satellite (ICESat), referred to as tie points (TP) and lowest-level elevation (LLE) methods, have been developed and applied in different studies. We reproduced these methods for the ICESat observation periods (2003–2008) in order to assess and analyse the sources of differences found in the retrieved freeboard and corresponding thickness estimates of the Arctic sea ice as produced by the Jet Propulsion Laboratory (JPL) and Goddard Space Flight Center (GSFC). Three main factors are found to affect the freeboard differences when applying these methods: (a) the approach used for calculation of the local sea surface references in leads (TP or LLE methods), (b) the along-track averaging scales used for this calculation, and (c) the corrections for lead width relative to the ICESat footprint and for snow depth accumulated in refrozen leads. The LLE method with 100 km averaging scale, as used to produce the GSFC data set, and the LLE method with a shorter averaging scale of 25 km both give larger freeboard estimates comparing to those derived by applying the TP method with 25 km averaging scale as used for the JPL product. Two factors, (a) and (b), contribute to the freeboard differences in approximately equal proportions, and their combined effect is, on average, about 6–7 cm. The effect of using different methods varies spatially: the LLE method tends to give lower freeboards (by up to 15 cm) over the thick multiyear ice and higher freeboards (by up to 10 cm) over first-year ice and the thin part of multiyear ice; the higher freeboards dominate. We show that the freeboard underestimation over most of these thinner parts of sea ice can be reduced to less than 2 cm when using the improved TP method proposed in this paper. The corrections for snow depth in leads and lead width, (c), are applied only for the JPL product and increase the freeboard estimates by about 7 cm on average. Thus, different approaches to calculating sea surface references and different along-track averaging scales from one side and the freeboard corrections as applied when producing the JPL data set from the other side roughly compensate each other with respect to freeboard estimation. Therefore, one may conclude that the difference in the mean sea ice thickness between the JPL and GSFC data sets reported in previous studies should be attributed mostly to different parameters used in the freeboard-to-thickness conversion.

We characterize sea-ice drift by applying a Lagrangian diffusion analysis to buoy trajectories from the International Arctic Buoy Programme (IABP) dataset and from two different models: the standalone Lagrangian sea-ice model neXtSIM and the Eulerian coupled ice–ocean model used for the TOPAZ reanalysis. By applying the diffusion analysis to the IABP buoy trajectories over the period 1979–2011, we confirm that sea-ice diffusion follows two distinct regimes (ballistic and Brownian) and we provide accurate values for the diffusivity and integral timescale that could be used in Eulerian or Lagrangian passive tracers models to simulate the transport and diffusion of particles moving with the ice. We discuss how these values are linked to the evolution of the fluctuating displacements variance and how this information could be used to define the size of the search area around the position predicted by the mean drift. By comparing observed and simulated sea-ice trajectories for three consecutive winter seasons (2007–2011), we show how the characteristics of the simulated motion may differ from or agree well with observations. This comparison illustrates the usefulness of first applying a diffusion analysis to evaluate the output of modeling systems that include a sea-ice model before using these in, e.g., oil spill trajectory models or, more generally, to simulate the transport of passive tracers in sea ice.

The Arctic sea ice cover has changed drastically over the last decades. Associated with these changes is a shift in dynamical regime seen by an increase of extreme fracturing events and an acceleration of sea ice drift. The highly non-linear dynamical response of sea ice to external forcing makes modelling these changes and the future evolution of Arctic sea ice a challenge for current models. It is, however, increasingly important that this challenge be better met, both because of the important role of sea ice in the climate system and because of the steady increase of industrial operations in the Arctic. In this paper we present a new dynamical/thermodynamical sea ice model called neXtSIM that is designed to address this challenge. neXtSIM is a continuous and fully Lagrangian model, whose momentum equation is discretised with the finite-element method. In this model, sea ice physics are driven by the combination of two core components: a model for sea ice dynamics built on a mechanical framework using an elasto-brittle rheology, and a model for sea ice thermodynamics providing damage healing for the mechanical framework. The evaluation of the model performance for the Arctic is presented for the period September 2007 to October 2008 and shows that observed multi-scale statistical properties of sea ice drift and deformation are well captured as well as the seasonal cycles of ice volume, area, and extent. These results show that neXtSIM is an appropriate tool for simulating sea ice over a wide range of spatial and temporal scales.

Leads within consolidated sea ice control heat exchange between the ocean and the atmosphere during winter, thus constituting an important climate parameter. These narrow elongated features occur when sea ice is fracturing under the action of wind and currents, reducing the local mechanical strength of the ice cover, which in turn impact the sea ice drift pattern. This creates a high demand for a high-quality lead fraction (LF) data set for sea ice model evaluation, initialization, and for the assimilation of such data in regional models. In this context, an available LF data set retrieved from satellite passive microwave observations (Advanced Microwave Scanning Radiometer – Earth Observing System, AMSR-E) is of great value, which has been providing pan-Arctic light- and cloud-independent daily coverage since 2002. In this study errors in this data set are quantified using accurate LF estimates retrieved from Synthetic Aperture Radar (SAR) images employing a threshold technique. A consistent overestimation of LF by a factor of 2–4 is found in the AMSR-E LF product. It is shown that a simple adjustment of the upper tie point used in the method to estimate the LF can reduce the pixel-wise error by a factor of 2 on average. Applying such an adjustment to the full data set may thus significantly increase the quality and value of the original data set.

We propose a method to reduce the error generated when computing sea ice deformation fields from synthetic aperture radar (SAR)-derived sea ice motion. The method is based on two steps. The first step consists of using a triangulation of the positions taken from the sea ice trajectories to define a mesh on which a first estimate of sea ice deformation is computed. The second step consists of applying a specific smoother to the deformation field to reduce the artificial noise that arises along discontinuities in the sea ice motion field. This method is here applied to RADARSAT Geophysical Processor System (RGPS) sea ice trajectories having a temporal and spatial resolution of about 3 days and 10 km, respectively. From the comparison between unfiltered and filtered fields, we estimate that the artificial noise causes an overestimation of about 60% of opening and closing. The artificial noise also has a strong impact on the statistical distribution of the deformation and on the scaling exponents estimated with multifractal analysis. We also show that a similar noise is present in the deformation fields provided in the widely used four-point deformation RGPS data set. These findings may have serious implications for previous studies as the constant overestimation of the opening and closing could lead to a large overestimation of freezing in leads, salt rejection and sea ice ridging.

We leveraged the power of ChatGPT and Bayesian optimization in the development of a multi-AI-driven system, backed by seven large language model-based assistants and equipped with machine learning algorithms, that seamlessly orchestrates a multitude of research aspects in a chemistry laboratory (termed the ChatGPT Research Group). Our approach accelerated the discovery of optimal microwave synthesis conditions, enhancing the crystallinity of MOF-321, MOF-322, and COF-323 and achieving the desired porosity and water capacity. In this system, human researchers gained assistance from these diverse AI collaborators, each with a unique role within the laboratory environment, spanning strategy planning, literature search, coding, robotic operation, labware design, safety inspection, and data analysis. Such a comprehensive approach enables a single researcher working in concert with AI to achieve productivity levels analogous to those of an entire traditional scientific team. Furthermore, by reducing human biases in screening experimental conditions and deftly balancing the exploration and exploitation of synthesis parameters, our Bayesian search approach precisely zeroed in on optimal synthesis conditions from a pool of 6 million within a significantly shortened time scale. This work serves as a compelling proof of concept for an AI-driven revolution in the chemistry laboratory, painting a future where AI becomes an efficient collaborator, liberating us from routine tasks to focus on pushing the boundaries of innovation.

Sea ice volume has decreased in the last decades, evoked by changes in sea ice area and thickness. Estimates of sea ice area and thickness rely on a number of geophysical parameters which introduce large uncertainties. To quantify these uncertainties we use freeboard retrievals from ICESat and investigate different assumptions about snow depth, sea ice density and area. We find that uncertainties in ice area are of minor importance for the estimates of sea ice volume during the cold season in the Arctic basin. The choice of mean ice density used when converting sea ice freeboard into thickness mainly influences the resulting mean sea ice thickness, while snow depth on top of the ice is the main driver for the year-to-year variability, particularly in late winter. The absolute uncertainty in the mean sea ice thickness is 0.28 m in February/March and 0.21 m in October/November. The uncertainty in snow depth contributes up to 70% of the total uncertainty and the ice density 30–35%, with higher values in October/November. We find large uncertainties in the total sea ice volume and trend. The mean total sea ice volume is 10 120 ± 1280 km3 in October/November and 13 250 ± 1860 km3 in February/March for the time period 2005–2007. Based on these uncertainties we obtain trends in sea ice volume of −1450 ± 530 km3 a−1 in October/November and −880 ± 260 km3 a−1 in February/March over the ICESat period (2003–2008). Our results indicate that, taking into account the uncertainties, the decline in sea ice volume in the Arctic between the ICESat (2003–2008) and CryoSat-2 (2010–2012) periods may have been less dramatic than reported in previous studies. However, more work and validation is required to quantify these changes and analyse possible unresolved biases in the freeboard retrievals.