Data and Code archive for the interpolation of surface ocean carbon dioxide

We have developed a statistical gap-filling method adapted to the specific coverage and properties of observed fugacity of surface ocean CO2 (fCO2). We have used this method to interpolate the Surface Ocean CO2 Atlas (SOCAT) v2 database on a 2.5°×2.5° global grid (south of 70°N) for 1985-2011 at monthly resolution. The method combines a spatial interpolation based on a 'radius of influence' to determine nearby similar fCO2 values with temporal harmonic and cubic spline curve-fitting, and also fits long term trends and seasonal cycles. Interannual variability is established using deviations of observations from the fitted trends and seasonal cycles. An uncertainty is computed for all interpolated values based on the spatial and temporal range of the interpolation. Tests of the method using model data show that it performs as well as or better than previous regional interpolation methods, but in addition it provides a near-global and interannual coverage.

Continuous 25-years aerosol records at coastal Antarctica: Part I. Variability and climate implications of ionic compounds

The aerosol climatology at the coastal Antarctic Neumayer Station (NM) was investigated based on continuous, 25-yr long observations of biogenic sulphur components (methanesulfonate and non-sea salt sulphate), sea salt and nitrate. Although significant long-term trends could only be detected for nitrate (-3.6 ± 2.5% per year between 1983 and 1993 and +4.0 ± 3.2% per year from 1993-2007), non-harmonic periodicities between 2 and 5 yr were typical for all species. Dedicated time series analyses revealed that relations to sea ice extent and various circulation indices are weak at best or not significant. In particular, no consistent link between sea ice extent and sea salt loadings was evident suggesting only a rather local relevance of the NM sea salt record. Nevertheless, a higher Southern Annular Mode index tended to entail a lower biogenic sulphur signal. In examining the spatial uniformity of the NM findings we contrasted them to respective 17 yr records from the coastal Dumont d'Urville Station. We found similar long-term trends for nitrate, indicating an Antarctic-wide but not identifiable atmospheric signal, although any significant impact of solar activity or pollution could be ruled out. No inter-site variability on the multiannual scale was evident for the other ionic compounds.

Gravity potential spherical harmonic series

Time variable gravity fields, reflecting variations of mass distribution in the system Earth is one of the key parameters to understand the changing Earth. Mass variations are caused either by redistribution of mass in, on or above the Earth's surface or by geophysical processes in the Earth's interior. The first set of observations of monthly variations of the Earth gravity field was provided by the US/German GRACE satellite mission beginning in 2002. This mission is still providing valuable information to the science community. However, as GRACE has outlived its expected lifetime, the geoscience community is currently seeking successor missions in order to maintain the long time series of climate change that was begun by GRACE. Several studies on science requirements and technical feasibility have been conducted in the recent years. These studies required a realistic model of the time variable gravity field in order to perform simulation studies on sensitivity of satellites and their instrumentation. This was the primary reason for the European Space Agency (ESA) to initiate a study on ''Monitoring and Modelling individual Sources of Mass Distribution and Transport in the Earth System by Means of Satellites''. The goal of this interdisciplinary study was to create as realistic as possible simulated time variable gravity fields based on coupled geophysical models, which could be used in the simulation processes in a controlled environment. For this purpose global atmosphere, ocean, continental hydrology and ice models were used. The coupling was performed by using consistent forcing throughout the models and by including water flow between the different domains of the Earth system. In addition gravity field changes due to solid Earth processes like continuous glacial isostatic adjustment (GIA) and a sudden earthquake with co-seismic and post-seismic signals were modelled. All individual model results were combined and converted to gravity field spherical harmonic series, which is the quantity commonly used to describe the Earth's global gravity field. The result of this study is a twelve-year time-series of 6-hourly time variable gravity field spherical harmonics up to degree and order 180 corresponding to a global spatial resolution of 1 degree in latitude and longitude. In this paper, we outline the input data sets and the process of combining these data sets into a coherent model of temporal gravity field changes. The resulting time series was used in some follow-on studies and is available to anybody interested.

Magnetic susceptibility, XRF and color reflectance data of DSDP Hole 72-516F, cores 113 and 114

Deep marine successions of early Campanian age from DSDP site 516F drilled at low paleolatitudes in the South Atlantic reveal distinct sub-Milankovitch variability in addition to precession and eccentricity related variations. Elemental abundance ratios point to a similar 5 climatic origin for these variations and exclude a quadripartite structure - as observed in the Mediterranean Neogene - of the precession related cycles as an explanation for the inferred semi-precession cyclicity in MS. However, the semi-precession cycle itself is likely an artifact, reflecting the first harmonic of the precession signal. The sub-Milankovitch variability is best approximated by a ~ 7 kyr cycle as shown by 10 spectral analysis and bandpass filtering. The presence of sub-Milankovitch cycles with a period similar to that of Heinrich events of the last glacial cycle is consistent with linking the latter to low-latitude climate change caused by a non-linear response to precession induced variations in insolation between the tropics.

Stable carbon and oxygen isotope ratios of foraminifera in surface sediments of the South China Sea

The relationship between planktonic and benthic foraminiferal stable-isotope values and oceanographic conditions and factors controlling isotopic variations are discussed on the basis of oxygen and carbon isotopic analyses of 192 modern surface and Last Glacial Maximum (LGM) samples from the South China Sea (SCS). The harmonic variation of benthic delta18O in surface sediments with water depth and temperature implies that the temperature is the main factor influencing benthic delta18O variations. Planktonic delta18O fluctuates with sea surface temperature (SST) and salinity (SSS). The N-S temperature gradient results in planktonic delta18O decreasing from the northeast to the south. Cool, saline waters driven by the winter monsoon are interpreted to have been responsible for the high delta18O values in the northeast SCS. The East Asian monsoons not only bring nutrients into the South China Sea and maintain high nutrient concentration levels at the southwestern and northeastern ends, which cause depleted delta13C both in planktonic (surface) and benthic (bottom) samples but also reduce planktonic/benthic delta18O differences. The distribution of delta18O and delta13C in the surface and LGM samples are strikingly similar, indicating that the impact of SST and SSS has been maintained, and nutrient inputs, mainly from the northeastern and southwestern ends, have been controlled by monsoons since the LGM. Comparisons of the modern and LGM delta18O indicate a difference of about 3.6 °C in bottom-water temperature and a large surface-to-bottom temperature gradient during the LGM as compared to today.

Code and readme of a MATLAB-based graphical user interface program for computing functionals of the geopotential

We present a novel graphical user interface program GrafLab (GRAvity Field LABoratory) for spherical harmonic synthesis (SHS) created in MATLAB®. This program allows to comfortably compute 38 various functionals of the geopotential up to ultra-high degrees and orders of spherical harmonic expansion. For the most difficult part of the SHS, namely the evaluation of the fully normalized associated Legendre functions (fnALFs), we used three different approaches according to required maximum degree: (i) the standard forward column method (up to maximum degree 1800, in some cases up to degree 2190); (ii) the modified forward column method combined with Horner's scheme (up to maximum degree 2700); (iii) the extended-range arithmetic (up to an arbitrary maximum degree). For the maximum degree 2190, the SHS with fnALFs evaluated using the extended-range arithmetic approach takes only approximately 2-3 times longer than its standard arithmetic counterpart, i.e. the standard forward column method. In the GrafLab, the functionals of the geopotential can be evaluated on a regular grid or point-wise, while the input coordinates can either be read from a data file or entered manually. For the computation on a regular grid we decided to apply the lumped coefficients approach due to significant time-efficiency of this method. Furthermore, if a full variance-covariance matrix of spherical harmonic coefficients is available, it is possible to compute the commission errors of the functionals. When computing on a regular grid, the output functionals or their commission errors may be depicted on a map using automatically selected cartographic projection.

Strontium measurements and sea surface temperatures of corals from Xisha Island, South China Sea

A Porites coral collected from Xisha Island, South China Sea, represents a skeleton secreted in the period from 1906 to 1994. The Sr contents of the coral vary linearly with the instrument-measured sea-surface temperature (SST), giving a Sr thermometer: SST = -1.9658 x Sr + 193.26. The reconstructed SST data show that the late 20th century was warmer (about 1°C) than the early 20th century and that two cooling (1915/1916 and 1947/1948) and three warming (1935/1936, 1960/1961, and 1976/1977) shifts occurred in the century. The temperature shifts are more pronounced for winters, implying a close effect of the west Pacific warm pool and Asian monsoon and suggesting that the former is a primary force controlling the climatic system of the region. Results of this study and previously published data indicate a close link of temperature shifts between the boreal summer and the austral winter or the boreal winter and the austral summer. The annual SST anomalies in the South China Sea and the South Pacific reveal the existence of harmonic but opposite SST variations between the two regions. On the decadal scale the comparative annual SST anomalies for the South China Sea and for the equatorial west Pacific show a similarity in temperature variations, implying that the South China Sea climate is coherent with climatic regime of the tropical west Pacific.

Temperature, precipitation and Normalized Difference Vegetation Index (NDVI) from 1982 to 2009 in the Koshi River Basin of the Middle Himalayas

The aim of this paper is to find out if there is a significant difference in using NDVI dataset processed by harmonic analysis method to evaluate its dynamic and response to climate change, compared with the original data.

Stable carbon and oxygen isotope record od IODP Hole 306-U1313

The Mid-Pleistocene transition (MPT) was the time when quasi-periodic (? 100 kyr), high-amplitude glacial variability developed in the absence of any significant change in the character of orbital forcing, leading to the establishment of the characteristic pattern of late Pleistocene climate variability. It has long been known that the interval around 900 ka stands out as a critical point of the MPT, when major glaciations started occurring most notably in the northern hemisphere. Here we examine the record of climatic conditions during this significant interval, using high-resolution stable isotope records from benthic and planktonic foraminifera from a sediment core in the North Atlantic (Integrated Ocean Drilling Program Expedition 306, Site U1313). We have considered the time interval from late in Marine Isotope Stage (MIS) 23 to MIS 20 (910 to 790 ka). Our data indicate that interglacial MIS 21 was a climatically unstable period and was broken into four interstadial periods, which have been identified and correlated across the North Atlantic region. These extra peaks tend to contradict previous studies that interpreted the MIS 21 variability as consisting essentially of a linear response to cyclical changes in orbital parameters. Cooling events in the surface record during MIS 21 were associated with low benthic carbon isotope excursions, suggesting a coupling between surface temperature changes and the strength of the Atlantic meridional overturning circulation. Time series analysis performed on the whole interval indicates that benthic and planktonic oxygen isotopes have significant concentrations of spectral power centered on periods of 10.7 kyr and 6 kyr, which is in agreement with the second and forth harmonic of precession. The excellent correspondence between the foraminifera d18O records and insolation variations at the Equator in March and September suggests that a mechanism related to low-latitude precession variations, advected to the high latitudes by tropical convective processes, might have generated such a response. This scenario accounts for the presence of oscillations at frequencies equal to precession harmonics at Site U1313, as well as the occurrence of higher amplitude oscillations between the MIS22/21 transition and most of MIS 21, times of enhanced insolation variability.

(Table 1) Annual mean concentrations of atmospheric a- and g-Hexachlorocyclohexane at stations Alert, Canada and Zeppelin, Norway

Twelve year datasets of weekly atmospheric concentrations of alpha- and gamma-HCH were compared between the two Arctic monitoring stations of Alert, Nunavut, Canada, and Zeppelin Mountain, Svalbard, Norway. Time-series analysis was conducted with the use of dynamic harmonic regression (DHR), which provided a very good model fit, to examine both the seasonal behaviour in these isomers and the longer-term, underlying trends. Strong spatial differences were not apparent between the two sites, although subtle differences in seasonal behaviour and composition were identified. For example, the composition of gamma-HCH to total HCH (alpha + gamma) was greater at Zeppelin compared to Alert, probably reflecting this site's proximity to major use regions of lindane. Pronounced seasonality in air concentrations for gamma-HCH was marked by a 'spring maximum event' (SME), confirming earlier studies. For alpha-HCH, the SME was much weaker and only evident at Alert, whereas at Zeppelin, seasonal fluctuations for alpha-HCH were marked by elevated concentrations in summer and lower concentrations during winter, with this pattern most apparent for the years after 2000. We attribute this difference in spatial and temporal patterns to the Arctic oscillation. A similar climatic pattern was not evident at either site in the gamma-HCH data. Seasonally adjusted, long-term trends revealed declining concentrations at both sites for alpha- and gamma-HCH over the entire time-series. Recent legislation affecting lindane use appear to account for this decline in gamma-HCH, with little evidence of a delay or 'lag' between the banning of lindane in Europe (a main source region) or Canada, and a decline in air concentrations observed at both Arctic sites.