Showing papers in "Earth System Science Data in 2012"

Database of diazotrophs in global ocean: abundance, biomass and nitrogen fixation rates

Abstract: . Marine N2 fixing microorganisms, termed diazotrophs, are a key functional group in marine pelagic ecosystems. The biological fixation of dinitrogen (N2) to bioavailable nitrogen provides an important new source of nitrogen for pelagic marine ecosystems and influences primary productivity and organic matter export to the deep ocean. As one of a series of efforts to collect biomass and rates specific to different phytoplankton functional groups, we have constructed a database on diazotrophic organisms in the global pelagic upper ocean by compiling about 12 000 direct field measurements of cyanobacterial diazotroph abundances (based on microscopic cell counts or qPCR assays targeting the nifH genes) and N2 fixation rates. Biomass conversion factors are estimated based on cell sizes to convert abundance data to diazotrophic biomass. The database is limited spatially, lacking large regions of the ocean especially in the Indian Ocean. The data are approximately log-normal distributed, and large variances exist in most sub-databases with non-zero values differing 5 to 8 orders of magnitude. Reporting the geometric mean and the range of one geometric standard error below and above the geometric mean, the pelagic N2 fixation rate in the global ocean is estimated to be 62 (52–73) Tg N yr−1 and the pelagic diazotrophic biomass in the global ocean is estimated to be 2.1 (1.4–3.1) Tg C from cell counts and to 89 (43–150) Tg C from nifH-based abundances. Reporting the arithmetic mean and one standard error instead, these three global estimates are 140 p 9.2 Tg N yr−1, 18 p 1.8 Tg C and 590 p 70 Tg C, respectively. Uncertainties related to biomass conversion factors can change the estimate of geometric mean pelagic diazotrophic biomass in the global ocean by about p70%. It was recently established that the most commonly applied method used to measure N2 fixation has underestimated the true rates. As a result, one can expect that future rate measurements will shift the mean N2 fixation rate upward and may result in significantly higher estimates for the global N2 fixation. The evolving database can nevertheless be used to study spatial and temporal distributions and variations of marine N2 fixation, to validate geochemical estimates and to parameterize and validate biogeochemical models, keeping in mind that future rate measurements may rise in the future. The database is stored in PANGAEA ( doi:10.1594/PANGAEA.774851 ).

Picophytoplankton biomass distribution in the global ocean

Abstract: . The smallest marine phytoplankton, collectively termed picophytoplankton, have been routinely enumerated by flow cytometry since the late 1980s during cruises throughout most of the world ocean. We compiled a database of 40 946 data points, with separate abundance entries for Prochlorococcus, Synechococcus and picoeukaryotes. We use average conversion factors for each of the three groups to convert the abundance data to carbon biomass. After gridding with 1° spacing, the database covers 2.4% of the ocean surface area, with the best data coverage in the North Atlantic, the South Pacific and North Indian basins, and at least some data in all other basins. The average picophytoplankton biomass is 12 p 22 μg C l−1 or 1.9 g C m−2. We estimate a total global picophytoplankton biomass of 0.53–1.32 Pg C (17–39% Prochlorococcus, 12–15% Synechococcus and 49–69% picoeukaryotes), with an intermediate/best estimate of 0.74 Pg C. Future efforts in this area of research should focus on reporting calibrated cell size and collecting data in undersampled regions. http://doi.pangaea.de/10.1594/PANGAEA.777385

A global diatom database – abundance, biovolume and biomass in the world ocean

Abstract: . Phytoplankton identification and abundance data are now commonly feeding plankton distribution databases worldwide. This study is a first attempt to compile the largest possible body of data available from different databases as well as from individual published or unpublished datasets regarding diatom distribution in the world ocean. The data obtained originate from time series studies as well as spatial studies. This effort is supported by the Marine Ecosystem Model Inter-Comparison Project (MAREMIP), which aims at building consistent datasets for the main plankton functional types (PFTs) in order to help validate biogeochemical ocean models by using carbon (C) biomass derived from abundance data. In this study we collected over 293 000 individual geo-referenced data points with diatom abundances from bottle and net sampling. Sampling site distribution was not homogeneous, with 58% of data in the Atlantic, 20% in the Arctic, 12% in the Pacific, 8% in the Indian and 1% in the Southern Ocean. A total of 136 different genera and 607 different species were identified after spell checking and name correction. Only a small fraction of these data were also documented for biovolumes and an even smaller fraction was converted to C biomass. As it is virtually impossible to reconstruct everyone's method for biovolume calculation, which is usually not indicated in the datasets, we decided to undertake the effort to document, for every distinct species, the minimum and maximum cell dimensions, and to convert all the available abundance data into biovolumes and C biomass using a single standardized method. Statistical correction of the database was also adopted to exclude potential outliers and suspicious data points. The final database contains 90 648 data points with converted C biomass. Diatom C biomass calculated from cell sizes spans over eight orders of magnitude. The mean diatom biomass for individual locations, dates and depths is 141.19 μg C l−1, while the median value is 11.16 μg C l−1. Regarding biomass distribution, 19% of data are in the range 0–1 μg C l−1, 29% in the range 1–10 μg C l−1, 31% in the range 10–100 μg C l−1, 18% in the range 100–1000 μg C l−1, and only 3% > 1000 μg C l−1. Interestingly, less than 50 species contributed to > 90% of global biomass, among which centric species were dominant. Thus, placing significant efforts on cell size measurements, process studies and C quota calculations of these species should considerably improve biomass estimates in the upcoming years. A first-order estimate of the diatom biomass for the global ocean ranges from 444 to 582 Tg C, which converts to 3 to 4 Tmol Si and to an average Si biomass turnover rate of 0.15 to 0.19 d−1. Link to the dataset: doi:10.1594/PANGAEA.777384 .

An 18-yr long (1993–2011) snow and meteorological dataset from a mid-altitude mountain site (Col de Porte, France, 1325 m alt.) for driving and evaluating snowpack models

Abstract: . A quality-controlled snow and meteorological dataset spanning the period 1 August 1993–31 July 2011 is presented, originating from the experimental station Col de Porte (1325 m altitude, Chartreuse range, France). Emphasis is placed on meteorological data relevant to the observation and modelling of the seasonal snowpack. In-situ driving data, at the hourly resolution, consist of measurements of air temperature, relative humidity, windspeed, incoming short-wave and long-wave radiation, precipitation rate partitioned between snow- and rainfall, with a focus on the snow-dominated season. Meteorological data for the three summer months (generally from 10 June to 20 September), when the continuity of the field record is not warranted, are taken from a local meteorological reanalysis (SAFRAN), in order to provide a continuous and consistent gap-free record. Data relevant to snowpack properties are provided at the daily (snow depth, snow water equivalent, runoff and albedo) and hourly (snow depth, albedo, runoff, surface temperature, soil temperature) time resolution. Internal snowpack information is provided from weekly manual snowpit observations (mostly consisting in penetration resistance, snow type, snow temperature and density profiles) and from a hourly record of temperature and height of vertically free ''settling'' disks. This dataset has been partially used in the past to assist in developing snowpack models and is presented here comprehensively for the purpose of multi-year model performance assessment. The data is placed on the PANGAEA repository ( http://dx.doi.org/10.1594/PANGAEA.774249 ) as well as on the public ftp server ftp://ftp-cnrm.meteo.fr/pub-cencdp/ .

A new 100-m Digital Elevation Model of the Antarctic Peninsula derived from ASTER Global DEM: methods and accuracy assessment

Abstract: . A high resolution surface topography Digital Elevation Model (DEM) is required to underpin studies of the complex glacier system on the Antarctic Peninsula. A complete DEM with better than 200 m pixel size and high positional and vertical accuracy would enable mapping of all significant glacial basins and provide a dataset for glacier morphology analyses. No currently available DEM meets these specifications. We present a new 100-m DEM of the Antarctic Peninsula (63–70° S), based on ASTER Global Digital Elevation Model (GDEM) data. The raw GDEM products are of high-quality on the rugged terrain and coastal-regions of the Antarctic Peninsula and have good geospatial accuracy, but they also contain large errors on ice-covered terrain and we seek to minimise these artefacts. Conventional data correction techniques do not work so we have developed a method that significantly improves the dataset, smoothing the erroneous regions and hence creating a DEM with a pixel size of 100 m that will be suitable for many glaciological applications. We evaluate the new DEM using ICESat-derived elevations, and perform horizontal and vertical accuracy assessments based on GPS positions, SPOT-5 DEMs and the Landsat Image Mosaic of Antarctica (LIMA) imagery. The new DEM has a mean elevation difference of −4 m (p 25 m RMSE) from ICESat (compared to −13 m mean and p97 m RMSE for the original ASTER GDEM), and a horizontal error of less than 2 pixels, although elevation accuracies are lower on mountain peaks and steep-sided slopes. The correction method significantly reduces errors on low relief slopes and therefore the DEM can be regarded as suitable for topographical studies such as measuring the geometry and ice flow properties of glaciers on the Antarctic Peninsula. The DEM is available for download from the NSIDC website: http://nsidc.org/data/nsidc-0516.html ( doi:10.5060/D47P8W9D ).

The global distribution of pteropods and their contribution to carbonate and carbon biomass in the modern ocean

Abstract: . Pteropods are a group of holoplanktonic gastropods for which global biomass distribution patterns remain poorly described. The aim of this study was to collect and synthesise existing pteropod (Gymnosomata, Thecosomata and Pseudothecosomata) abundance and biomass data, in order to evaluate the global distribution of pteropod carbon biomass, with a particular emphasis on temporal and spatial patterns. We collected 25 939 data points from several online databases and 41 scientific articles. These data points corresponded to observations from 15 134 stations, where 93% of observations were of shelled pteropods (Thecosomata) and 7% of non-shelled pteropods (Gymnosomata). The biomass data has been gridded onto a 360 × 180° grid, with a vertical resolution of 33 depth levels. Both the raw data file and the gridded data in NetCDF format can be downloaded from PANGAEA, doi:10.1594/PANGAEA.777387 . Data were collected between 1950–2010, with sampling depths ranging from 0–2000 m. Pteropod biomass data was either extracted directly or derived through converting abundance to biomass with pteropod-specific length to carbon biomass conversion algorithms. In the Northern Hemisphere (NH), the data were distributed quite evenly throughout the year, whereas sampling in the Southern Hemisphere (SH) was biased towards winter and summer values. 86% of all biomass values were located in the NH, most (37%) within the latitudinal band of 30–60° N. The range of global biomass values spanned over four orders of magnitude, with mean and median (non-zero) biomass values of 4.6 mg C m−3 (SD = 62.5) and 0.015 mg C m−3, respectively. The highest mean biomass was located in the SH within the 70–80° S latitudinal band (39.71 mg C m−3, SD = 93.00), while the highest median biomass was in the NH, between 40–50° S (0.06 mg C m−3, SD = 79.94). Shelled pteropods constituted a mean global carbonate biomass of 23.17 mg CaCO3 m−3 (based on non-zero records). Total biomass values were lowest in the equatorial regions and equally high at both poles. Pteropods were found at least to depths of 1000 m, with the highest biomass values located in the surface layer (0–10 m) and gradually decreasing with depth, with values in excess of 100 mg C m−3 only found above 200 m depth. Tropical species tended to concentrate at greater depths than temperate or high-latitude species. Global biomass levels in the NH were relatively invariant over the seasonal cycle, but more seasonally variable in the SH. The collected database provides a valuable tool for modellers for the study of marine ecosystem processes and global biogeochemical cycles. By extrapolating regional biomass to a global scale, we established global pteropod biomass to add up to 500 Tg C.

Global marine plankton functional type biomass distributions: Phaeocystis spp.

Abstract: . The planktonic haptophyte Phaeocystis has been suggested to play a fundamental role in the global biogeochemical cycling of carbon and sulphur, but little is known about its global biomass distribution. We have collected global microscopy data of the genus Phaeocystis and converted abundance data to carbon biomass using species-specific carbon conversion factors. Microscopic counts of single-celled and colonial Phaeocystis were obtained both through the mining of online databases and by accepting direct submissions (both published and unpublished) from Phaeocystis specialists. We recorded abundance data from a total of 1595 depth-resolved stations sampled between 1955–2009. The quality-controlled dataset includes 5057 counts of individual Phaeocystis cells resolved to species level and information regarding life-stages from 3526 samples. 83% of stations were located in the Northern Hemisphere while 17% were located in the Southern Hemisphere. Most data were located in the latitude range of 50–70° N. While the seasonal distribution of Northern Hemisphere data was well-balanced, Southern Hemisphere data was biased towards summer months. Mean species- and form-specific cell diameters were determined from previously published studies. Cell diameters were used to calculate the cellular biovolume of Phaeocystis cells, assuming spherical geometry. Cell biomass was calculated using a carbon conversion factor for prymnesiophytes. For colonies, the number of cells per colony was derived from the colony volume. Cell numbers were then converted to carbon concentrations. An estimation of colonial mucus carbon was included a posteriori, assuming a mean colony size for each species. Carbon content per cell ranged from 9 pg C cell−1 (single-celled Phaeocystis antarctica) to 29 pg C cell−1 (colonial Phaeocystis globosa). Non-zero Phaeocystis cell biomasses (without mucus carbon) range from 2.9 × 10−5 to 5.4 × 103 μg C l−1, with a mean of 45.7 μg C l−1 and a median of 3.0 μg C l−1. The highest biomasses occur in the Southern Ocean below 70° S (up to 783.9 μg C l−1) and in the North Atlantic around 50° N (up to 5.4 × 103 μg C l−1). The original and gridded data can be downloaded from PANGAEA, doi:10.1594/PANGAEA.779101 .

Two weather radar time series of the altitude of the volcanic plume during the May 2011 eruption of Grímsvötn, Iceland

Abstract: . The eruption of Grimsvotn volcano in Iceland in 2011 lasted for a week, 21–28 May. The eruption was explosive and peaked during the first hours, with the eruption plume reaching 20–25 km altitude. The height of the plume was monitored every 5 min with a C-band weather radar located at Keflavik International Airport and a mobile X-band radar, 257 km and 75 km distance from the volcano respectively. In addition, photographs taken during the first half-hour of the eruption give information regarding the initial rise. Time series of the plume-top altitude were constructed from the radar observations. This paper presents the two independent radar time series. The series have been cross validated and there is a good agreement between them. The echo top radar series of the altitude of the volcanic plume are publicly available from the Pangaea Data Publisher ( doi:10.1594/PANGAEA.778390 ).

Future Flows Climate: an ensemble of 1-km climate change projections for hydrological application in Great Britain

Abstract: . The dataset Future Flows Climate was developed as part of the project ''Future Flows and Groundwater Levels'' to provide a consistent set of climate change projections for the whole of Great Britain at both space and time resolutions appropriate for hydrological applications, and to enable climate change uncertainty and climate variability to be accounted for in the assessment of their possible impacts on the environment. Future Flows Climate is derived from the Hadley Centre's ensemble projection HadRM3-PPE that is part of the basis of UKCP09 and includes projections in available precipitation (water available to hydrological processes after snow and ice storages have been accounted for) and potential evapotranspiration. It corresponds to an 11-member ensemble of transient projections from January 1950 to December 2098, each a single realisation from a different variant of HadRM3. Data are provided on a 1-km grid over the HadRM3 land areas at a daily (available precipitation) and monthly (PE) time step as netCDF files. Because systematic biases in temperature and precipitation were found between HadRM3-PPE and gridded temperature and precipitation observations for the 1962–1991 period, a monthly bias correction procedure was undertaken, based on a linear correction for temperature and a quantile-mapping correction (using the gamma distribution) for precipitation followed by a spatial downscaling. Available precipitation was derived from the bias-corrected precipitation and temperature time series using a simple elevation-dependant snow-melt model. Potential evapotranspiration time series were calculated for each month using the FAO-56 Penman-Monteith equations and bias-corrected temperature, cloud cover, relative humidity and wind speed from HadRM3-PPE along with latitude of the grid and the day of the year. Future Flows Climate is freely available for non-commercial use under certain licensing conditions. It is the dataset used to generate Future Flows Hydrology, an ensemble of transient projections of daily river flow and monthly groundwater time series for representative river basins and boreholes in Great Britain. doi:10.5285/bad1514f-119e-44a4-8e1e-442735bb9797 .

First-order estimate of the planktic foraminifer biomass in the modern ocean

Abstract: . Planktic foraminifera are heterotrophic mesozooplankton of global marine abundance. The position of planktic foraminifers in the marine food web is different compared to other protozoans and ranges above the base of heterotrophic consumers. Being secondary producers with an omnivorous diet, which ranges from algae to small metazoans, planktic foraminifers are not limited to a single food source, and are assumed to occur at a balanced abundance displaying the overall marine biological productivity at a regional scale. With a new non-destructive protocol developed from the bicinchoninic acid (BCA) method and nano-photospectrometry, we have analysed the protein-biomass, along with test size and weight, of 754 individual planktic foraminifers from 21 different species and morphotypes. From additional CHN analysis, it can be assumed that protein-biomass equals carbon-biomass. Accordingly, the average individual planktic foraminifer protein- and carbon-biomass amounts to 0.845 μg. Samples include symbiont bearing and symbiont-barren species from the sea surface down to 2500 m water depth. Conversion factors between individual biomass and assemblage-biomass are calculated for test sizes between 72 and 845 μm (minimum test diameter). Assemblage-biomass data presented here include 1128 sites and water depth intervals. The regional coverage of data includes the North Atlantic, Arabian Sea, Red Sea, and Caribbean as well as literature data from the eastern and western North Pacific, and covers a wide range of oligotrophic to eutrophic waters over six orders of magnitude of planktic-foraminifer assemblage-biomass (PFAB). A first order estimate of the average global planktic foraminifer biomass production (>125 μm) ranges from 8.2–32.7 Tg C yr−1 (i.e. 0.008–0.033 Gt C yr−1), and might be more than three times as high including neanic and juvenile individuals adding up to 25–100 Tg C yr−1. However, this is a first estimate of regional PFAB extrapolated to the global scale, and future estimates based on larger data sets might considerably deviate from the one presented here. This paper is supported by, and a contribution to the Marine Ecosystem Data project (MAREDAT). Data are available from http://www.pangaea.de ( http://doi.pangaea.de/10.1594/PANGAEA.777386 ).

Proglacial river stage, discharge, and temperature datasets from the Akuliarusiarsuup Kuua River northern tributary, Southwest Greenland, 2008–2011

Abstract: . Pressing scientific questions concerning the Greenland ice sheet's climatic sensitivity, hydrology, and contributions to current and future sea level rise require hydrological datasets to resolve. While direct observations of ice sheet meltwater losses can be obtained in terrestrial rivers draining the ice sheet and from lake levels, few such datasets exist. We present a new hydrologic dataset from previously unmonitored sites in the vicinity of Kangerlussuaq, Southwest Greenland. This dataset contains measurements of river stage and discharge for three sites along the Akuliarusiarsuup Kuua (Watson) River's northern tributary, with 30 min temporal resolution between June 2008 and July 2011. Additional data of water temperature, air pressure, and lake stage are also provided. Flow velocity and depth measurements were collected at sites with incised bedrock or structurally reinforced channels to maximize data quality. However, like most proglacial rivers, high turbulence and bedload transport introduce considerable uncertainty to the derived discharge estimates. Eleven propagating error sources were quantified, and reveal that largest uncertainties are associated with flow depth observations. Mean discharge uncertainties (approximately the 68% confidence interval) are two to four times larger (p19% to p43%) than previously published estimates for Greenland rivers. Despite these uncertainties, this dataset offers a rare collection of direct measurements of ice sheet runoff to the global ocean and is freely available for scientific use at http://dx.doi.org/10.1594/PANGAEA.762818 .

Picoheterotroph ( Bacteria and Archaea ) biomass distribution in the global ocean

Abstract: . We compiled a database of 39 766 data points consisting of flow cytometric and microscopical measurements of picoheterotroph abundance, including both Bacteria and Archaea. After gridding with 1° spacing, the database covers 1.3% of the ocean surface. There are data covering all ocean basins and depths except the Southern Hemisphere below 350 m or from April until June. The average picoheterotroph biomass is 3.9 p 3.6 μg C l−1 with a 20-fold decrease between the surface and the deep sea. We estimate a total ocean inventory of about 1.3 × 1029 picoheterotroph cells. Surprisingly, the abundance in the coastal regions is the same as at the same depths in the open ocean. Using an average of published open ocean measurements for the conversion from abundance to carbon biomass of 9.1 fg cell−1, we calculate a picoheterotroph carbon inventory of about 1.2 Pg C. The main source of uncertainty in this inventory is the conversion factor from abundance to biomass. Picoheterotroph biomass is ~2 times higher in the tropics than in the polar oceans. doi:10.1594/PANGAEA.779142

Homogenization of Portuguese long-term temperature data series: Lisbon, Coimbra and Porto

Abstract: . Three long-term temperature data series measured in Portugal were studied to detect and correct non-climatic homogeneity breaks and are now available for future studies of climate variability. Series of monthly minimum (Tmin) and maximum (Tmax) temperatures measured in the three Portuguese meteorological stations of Lisbon (from 1856 to 2008), Coimbra (from 1865 to 2005) and Porto (from 1888 to 2001) were studied to detect and correct non-climatic breaks. These series, together with monthly series of average temperature (Taver) and temperature range (DTR) derived from them, were tested in order to detect breaks, using firstly metadata, secondly a visual analysis, and thirdly four widely used homogeneity tests: von Neumann ratio test, Buishand test, standard normal homogeneity test, and Pettitt test. The homogeneity tests were used in absolute (using temperature series themselves) and relative (using sea-surface temperature anomalies series obtained from HadISST2.0.0.0 close to the Portuguese coast or already corrected temperature series as reference series) modes. We considered the Tmin, Tmax and DTR series as most informative for the detection of breaks due to the fact that Tmin and Tmax could respond differently to changes in position of a thermometer or other changes in the instrument's environment; Taver series have been used mainly as control. The homogeneity tests showed strong inhomogeneity of the original data series, which could have both internal climatic and non-climatic origins. Breaks that were identified by the last three mentioned homogeneity tests were compared with available metadata containing data such as instrument changes, changes in station location and environment, observation procedures, etc. Significant breaks (significance 95% or more) that coincided with known dates of instrumental changes were corrected using standard procedures. It was also noted that some significant breaks, which could not be connected to known dates of any changes in the park of instruments or stations location and environment, were probably caused by large volcanic eruptions. The corrected series were again tested for homogeneity; the corrected series were considered free of non-climatic breaks when the tests of most of monthly series showed no significant (significance 95% or more) breaks that coincide with dates of known instrument changes. Corrected series are now available within the framework of ERA-CLIM FP7 project for future studies of climate variability ( doi:10.1594/PANGAEA.785377 ).

King George Island ice cap geometry updated with airborne GPR measurements

Abstract: . Ice geometry is a mandatory requirement for numerical modelling purposes. In this paper we present a consistent data set for the ice thickness, the bedrock topography and the ice surface topography of the King George Island ice cap (Arctowski icefield and the adjacent central part). The new data set is composed of ground based and airborne ground penetrating radar (GPR) and differential GPS (DGPS) measurements, obtained during several field campaigns. Blindow et al. (2010) already provided a comprehensive overview of the ground based measurements carried out in the safely accessible area of the ice cap. The updated data set incorporates airborne measurements in the heavily crevassed coastal areas. Therefore, in this paper special attention is paid to the airborne measurements by addressing the instrument used, survey procedure, and data processing in more detail. In particular, the inclusion of airborne GPR measurements with the 30 MHz BGR-P30-System developed at the Institute of Geophysics (University of Munster) completes the picture of the ice geometry substantially. The compiled digital elevation model of the bedrock shows a rough, highly variable topography with pronounced valleys, ridges, and troughs. Mean ice thickness is 240 p 6 m, with a maximum value of 422 p 10 m in the surveyed area. Noticeable are bounded areas in the bedrock topography below sea level where marine based ice exists. The provided data set is required as a basis for future monitoring attempts or as input for numerical modelling experiments. The data set is available from the PANGAEA database at http://dx.doi.org/10.1594/PANGAEA.770567 .

A compilation of tropospheric measurements of gas-phase and aerosol chemistry in polar regions

Abstract: . Measurements of atmospheric chemistry in polar regions have been made for more than half a century. Probably the first Antarctic ozone data were recorded in 1958 during the International Geophysical Year. Since then, many measurement campaigns followed, and the results are now spread over many publications in several journals. Here, we have compiled measurements of tropospheric gas-phase and aerosol chemistry made in the Arctic and the Antarctic. It is hoped that this data collection is worth more than the sum of its components and serves as a basis for future analyses of spatial and temporal trends in polar atmospheric chemistry.

Dobson, Brewer, ERA-40 and ERA-Interim original and merged total ozone data sets – evaluation of differences: a case study, Hradec Králové (Czech), 1961–2010

Abstract: . Homogenized data series of total ozone measurements taken by the regularly and well calibrated Dobson and Brewer spectrophotometers at Hradec Kralove (Czech) and the data from the re-analyses ERA-40 and ERA-Interim were merged and compared to investigate differences between the particular data sets originated in Central Europe, the Northern Hemisphere (NH) mid-latitudes. The Dobson-to-Brewer transfer function and the algorithm for approximation of the data from the re-analyses were developed, tested and applied for creation of instrumentally consistent and completed total ozone data series of the 50-yr period 1961–2010 of observations. This correction has reduced the well-known seasonal differences between Dobson and Brewer data below the 1% calibration limit of the spectrophotometers. Incorporation of the ERA-40 and ERA-Interim total ozone data on days with missing measurements significantly improved completeness and reliability of the data series mainly in the first two decades of the period concerned. Consistent behaviour of the original and corrected/merged data sets was found in the pre-ozone-hole period (1961–1985). In the post-Pinatubo (1994–2010) era the data series show seasonal differences that can introduce uncertainty in estimation of ozone recovery mainly in the winter-spring season when the effect of the Montreal Protocol and its Amendments is expected. All the data sets confirm substantial depletion of ozone also in the summer months that gives rise to the question about its origin. The merged and completed data series of total ozone will be further analyzed to quantify chemical ozone losses and contribution of natural atmospheric processes to the ozone depletion over the region. This case study points out the importance of selection and evaluation of the quality and consistency of the input data sets used in estimation of long-term ozone changes including recovery of the ozone layer over the selected areas. Data are available from the PANGAEA database at doi:10.1594/PANGAEA.779819 .