Assimilation of sea-ice concentration in a global climate model — physical and statistical aspects

We investigate the initialisation of Northern Hemisphere sea ice in the global climate model ECHAM5/MPI-OM by assimilating sea-ice concentration data. The analysis updates for concentration are given by Newtonian relaxation, and we discuss different ways of specifying the analysis updates for mean thickness. Because the conservation of mean ice thickness or actual ice thickness in the analysis updates leads to poor assimilation performance, we introduce a proportional dependence between concentration and mean thickness analysis updates. Assimilation with these proportional mean-thickness analysis updates leads to good assimilation performance for sea-ice concentration and thickness, both in identical-twin experiments and when assimilating sea-ice observations. The simulation of other Arctic surface fields in the coupled model is, however, not significantly improved by the assimilation. To understand the physical aspects of assimilation errors, we construct a simple prognostic model of the sea-ice thermodynamics, and analyse its response to the assimilation. We find that an adjustment of mean ice thickness in the analysis update is essential to arrive at plausible state estimates. To understand the statistical aspects of assimilation errors, we study the model background error covariance between ice concentration and ice thickness. We find that the spatial structure of covariances is best represented by the proportional mean-thickness analysis updates. Both physical and statistical evidence supports the experimental finding that assimilation with proportional mean-thickness updates outperforms the other two methods considered. The method described here is very simple to implement, and gives results that are sufficiently good to be used for initialising sea ice in a global climate model for seasonal to decadal predictions

Satellite passive microwave sea-ice concentration data set inter-comparison for Arctic summer conditions

We report on results of a systematic inter-comparison of 10 global sea-ice concentration (SIC) data products at 12.5 to 50.0 km grid resolution from satellite passive microwave (PMW) observations for the Arctic during summer. The products are compared against SIC and net ice surface fraction (ISF) - SIC minus the per-grid-cell melt pond fraction (MPF) on sea ice - as derived from MODerate resolution Imaging Spectroradiometer (MODIS) satellite observations and observed from ice-going vessels. Like in Kern et al. (2019), we group the 10 products based on the concept of the SIC retrieval used. Group I consists of products of the European Organisation for the Exploitation of Meteorological Satellites (EUMETSAT) Ocean and Sea Ice Satellite Application Facility (OSI SAF) and European Space Agency (ESA) Climate Change Initiative (CCI) algorithms. Group II consists of products derived with the Comiso bootstrap algorithm and the National Oceanographic and Atmospheric Administration (NOAA) National Snow and Ice Data Center (NSIDC) SIC climate data record (CDR). Group III consists of Arctic Radiation and Turbulence Interaction Study (ARTIST) Sea Ice (ASI) and National Aeronautics and Space Administration (NASA) Team (NT) algorithm products, and group IV consists of products of the enhanced NASA Team algorithm (NT2). We find widespread positive and negative differences between PMW and MODIS SIC with magnitudes frequently reaching up to 20 %-25 % for groups I and III and up to 30 %-35 % for groups II and IV. On a pan-Arctic scale these differences may cancel out: Arctic average SIC from group I products agrees with MODIS within 2 %-5 % accuracy during the entire melt period from May through September. Group II and IV products overestimate MODIS Arctic average SIC by 5 %-10 %. Out of group III, ASI is similar to group I products while NT SIC underestimates MODIS Arctic average SIC by 5 %-10 %. These differences, when translated into the impact computing Arctic sea-ice area (SIA), match well with the differences in SIA between the four groups reported for the summer months by Kern et al. (2019). MODIS ISF is systematically overestimated by all products; NT provides the smallest overestimations (up to 25 %) and group II and IV products the largest overestimations (up to 45 %). The spatial distribution of the observed overestimation of MODIS ISF agrees reasonably well with the spatial distribution of the MODIS MPF and we find a robust linear relationship between PMW SIC and MODIS ISF for group I and III products during peak melt, i.e. July and August. We discuss different cases taking into account the expected influence of ice surface properties other than melt ponds, i.e. wet snow and coarse-grained snow/refrozen surface, on brightness temperatures and their ratios used as input to the SIC retrieval algorithms. Based on this discussion we identify the mismatch between the actually observed surface properties and those represented by the ice tie points as the most likely reason for (i) the observed differences between PMW SIC and MODIS ISF and for (ii) the often surprisingly small difference between PMW and MODIS SIC in areas of high melt pond fraction. We conclude that all 10 SIC products are highly inaccurate during summer melt. We hypothesize that the unknown number of melt pond signatures likely included in the ice tie points plays an important role - particularly for groups I and II - and recommend conducting further research in this field

Zwei Jahre pro-Q – Bilanz eines biokonformen Eutergesundheitsprogrammes

The pro-q-project is an udder health programme especially for organic farmers under Swiss conditions. The aim is to minimize the use of s antibiotics and to assure a good udder health state on the farms. At the end of the third project year of 100 project farms 16 participate for more than two years. The udder health situation of these 16 farms improved moderately in the second project year (37% of the cows with a so-matic cell count > 100’000/ml) compared to the year before project start (40% of the cows with a somatic cell count > 100’000/ml). Furthermore, a nearly linear increase of the arithmetic mean of herd lactation number (3.27 year before project start, 3.33 first project year, 3.44 second project year, resp.) could be found. The number of antibiotic treatments per cow and year decreased from 0.37 in the year before project start to 0.24 in the second project year

Management factors affecting udder health and effects of a one year extension program in organic dairy herds

The first part of this study was a cross-sectional analysis of the impact of 29 management factors on udder health in organic dairy farms in Switzerland. All 77 farms joined the extension program ‘pro-Q'. As a measure of udder health the theoretical bulk milk somatic cell count (TBMSCC) calculated by the monthly cow composite somatic cell count over a time period of 1 year was chosen. The basic udder health of the farms was determined by TBMSCC during the year prior to the start of the project (mean for all farms = 176 460 cells/ml). In the multivariable analysis, the five factors ‘swiss brown breed', ‘alpine summer pasturing', ‘calf feeding with milk from mastitis diseased cows', ‘hard bedding' and ‘no post-milking' remained as significant risk factors on udder health. In the second part of the study, the development of management factors and the udder health situation affected by an extension program after 1 year was investigated. A partial improvement of the management factors on the farms but no overall improvement on udder health and no association between management changes and udder health changes were found. Improvement of udder health was more likely in farms with higher basic TBMSCC than in those farms with less udder health problems at the beginning of the projec

Trends in sea-ice variability on the way to an ice-free Arctic

This is the author accepted manuscript. The final version is available from the publisher via the DOI in this record. The final author version was published under the title: Statistical indicators of Arctic sea-ice stability-prospects and limitations and is available in ORE via https://ore.exeter.ac.uk/repository/handle/10871/23493It has been widely debated whether Arctic sea-ice loss can reach a tipping point beyond which a large sea-ice area disappears abruptly. The theory of dynamical systems predicts a slowing down when a system destabilises towards a tipping point. In simple stochastic systems this can result in increasing variance and autocorrelation, potentially yielding an early warning of an abrupt change. Here we aim to establish whether the loss of Arctic sea ice would follow these conceptual predictions, and which trends in sea ice variability can be expected from pre-industrial conditions toward an Arctic that is ice-free during the whole year. To this end, we apply a model hierarchy consisting of two box models and one comprehensive Earth system model. We find a consistent and robust decrease of the ice volume's annual relaxation time before summer ice is lost because thinner ice can adjust more quickly to perturbations. Thereafter, the relaxation time increases, mainly because the system becomes dominated by the ocean water's large heat capacity when the ice-free season becomes longer. Both trends carry over to the autocorrelation of sea ice thickness in time series. These changes are robust to the nature and origin of climate variability in the models and hardly depend on the balance of feedbacks. Therefore, characteristic trends can be expected in the future. As these trends are not specific to the existence of abrupt ice loss, the prospects for early warnings seem very limited. This result also has implications for statistical indicators in other systems whose effective mass changes over time, affecting the trend of their relaxation time. However, the robust relation between state and variability would allow an estimate of sea-ice variability from only short observations. This could help one to estimate the likelihood and persistence of extreme events in the future.This work was carried out under the programme of the Netherlands Earth System Science Centre (NESSC), financially supported by the Ministry of Education, Culture and Science (OCW). We also acknowledge the World Climate Research Programme’s Working Group on Coupled Modelling, which is responsible for CMIP, and we thank the climate modelling groups for producing and making available their model output. We thank Vasilis Dakos for helping to apply his early warnings R package and Chao Li for making available the MPI-ESM model output. S. B. gratefully acknowledges Arie Staal for his fruitful and revealing approaches to savour scientific achievements. We are also indebted to Till Wagner and Ian Eisenman for their valuable comments and their very amiable and cooperative spirit

Statistical indicators of Arctic sea-ice stability-prospects and limitations

This is the final version of the article. Available from the European Geosciences Union via the DOI in this record.We examine the relationship between the mean and the variability of Arctic sea-ice coverage and volume in a large range of climates from globally ice-covered to globally ice-free conditions. Using a hierarchy of two column models and several comprehensive Earth system models, we consolidate the results of earlier studies and show that mechanisms found in simple models also dominate the interannual variability of Arctic sea ice in complex models. In contrast to predictions based on very idealised dynamical systems, we find a consistent and robust decrease of variance and autocorrelation of sea-ice volume before summer sea ice is lost. We attribute this to the fact that thinner ice can adjust more quickly to perturbations. Thereafter, the autocorrelation increases, mainly because it becomes dominated by the ocean water's large heat capacity when the ice-free season becomes longer. We show that these changes are robust to the nature and origin of climate variability in the models and do not depend on whether Arctic sea-ice loss occurs abruptly or irreversibly. We also show that our climate is changing too rapidly to detect reliable changes in autocorrelation of annual time series. Based on these results, the prospects of detecting statistical early warning signals before an abrupt sea-ice loss at a "tipping point" seem very limited. However, the robust relation between state and variability can be useful to build simple stochastic climate models and to make inferences about past and future sea-ice variability from only short observations or reconstructions.This work was carried out under the programme of the Netherlands Earth System Science Centre (NESSC), financially supported by the Ministry of Education, Culture and Science (OCW). We also acknowledge the World Climate Research Programme’s Working Group on Coupled Modelling, which is responsible for CMIP, and we thank the climate modelling groups for producing and making available their model output. We thank Vasilis Dakos for helping to apply his early warnings R package and Chao Li for making available the MPI-ESM model output. S. B. gratefully acknowledges Arie Staal for his fruitful and revealing approaches to savour scientific achievements. We are also indebted to Till Wagner and Ian Eisenman for their valuable comments and their very amiable and cooperative spirit. Finally, we acknowledge two anonymous reviewers who helped us to improve the manuscript

ICON-O: The Ocean Component of the ICON Earth System Model - Global simulation characteristics and local telescoping capability

Abstract We describe the ocean general circulation model ICON-O of the Max Planck Institute for Meteorology, which forms the ocean-sea ice component of the Earth system model ICON-ESM. ICON-O relies on innovative structure-preserving finite volume numerics. We demonstrate the fundamental ability of ICON-O to simulate key features of global ocean dynamics at both uniform and non-uniform resolution. Two experiments are analyzed and compared with observations, one with a nearly uniform and eddy-rich resolution of ?10?km and another with a telescoping configuration whose resolution varies smoothly from globally ?80?km to ?10?km in a focal region in the North Atlantic. Our results show first, that ICON-O on the nearly uniform grid simulates an ocean circulation that compares well with observations and second, that ICON-O in its telescope configuration is capable of reproducing the dynamics in the focal region over decadal time scales at a fraction of the computational cost of the uniform-grid simulation. The telescopic technique offers an alternative to the established regionalization approaches. It can be used either to resolve local circulation more accurately or to represent local scales that cannot be simulated globally while remaining within a global modeling framework

Satellite passive microwave sea-ice concentration data set intercomparison: closed ice and ship-based observations

We report on results of a systematic intercomparison of 10 global sea-ice concentration (SIC) data products at 12.5 to 50.0 km grid resolution for both the Arctic and the Antarctic. The products are compared with each other with respect to differences in SIC, sea-ice area (SIA), and sea-ice extent (SIE), and they are compared against a global wintertime near-100 % reference SIC data set for closed pack ice conditions and against global year-round ship-based visual observations of the sea-ice cover. We can group the products based on the concept of their SIC retrieval algorithms. Group I consists of data sets using the self-optimizing EUMETSAT OSI SAF and ESA CCI algorithms. Group II includes data using the Comiso bootstrap algorithm and the NOAA NSIDC sea-ice concentration climate data record (CDR). The standard NASA Team and the ARTIST Sea Ice (ASI) algorithms are put into group III, and NASA Team 2 is the only element of group IV. The three CDRs of group I (SICCI-25km, SICCI-50km, and OSI-450) are biased low compared to a 100 % reference SIC data set with biases of - 0.4 % to -1.0 % (Arctic) and -0.3 % to -1.1 % (Antarctic). Products of group II appear to be mostly biased high in the Arctic by between +1.0 % and +3.5 %, while their biases in the Antarctic range from -0.2 % to +0.9 %. Group III product biases are different for the Arctic, +0.9 % (NASA Team) and -3.7 % (ASI), but similar for the Antarctic, -5.4 % and -5.6 %, respectively. The standard deviation is smaller in the Arctic for the quoted group I products (1.9 % to 2.9 %) and Antarctic (2.5 % to 3.1 %) than for group II and III products: 3.6 % to 5.0 % for the Arctic and 4.0 % to 6.5 % for the Antarctic. We refer to the paper to understand why we could not give values for group IV here. We discuss the impact of truncating the SIC distribution, as naturally retrieved by the algorithms around the 100 % sea-ice concentration end. We show that evaluation studies of such truncated SIC products can result in misleading statistics and favour data sets that systematically overestimate SIC. We describe a method to reconstruct the non-truncated distribution of SIC before the evaluation is performed. On the basis of this evaluation, we open a discussion about the overestimation of SIC in data products, with far-reaching consequences for surface heat flux estimations in winter. We also document inconsistencies in the behaviour of the weather filters used in products of group II, and we suggest advancing studies about the influence of these weather filters on SIA and SIE time series and their trends