Mean dynamic topography and oceanographic parameters estimated from an inverse model and satellite geodesy, with link to model result in one single NetCDF file (392 MB), including the inverse data error covariance

Geostrophic surface velocities can be derived from the gradients of the mean dynamic topography-the difference between the mean sea surface and the geoid. Therefore, independently observed mean dynamic topography data are valuable input parameters and constraints for ocean circulation models. For a successful fit to observational dynamic topography data, not only the mean dynamic topography on the particular ocean model grid is required, but also information about its inverse covariance matrix. The calculation of the mean dynamic topography from satellite-based gravity field models and altimetric sea surface height measurements, however, is not straightforward. For this purpose, we previously developed an integrated approach to combining these two different observation groups in a consistent way without using the common filter approaches (Becker et al. in J Geodyn 59(60):99-110, 2012, doi:10.1016/j.jog.2011.07.0069; Becker in Konsistente Kombination von Schwerefeld, Altimetrie und hydrographischen Daten zur Modellierung der dynamischen Ozeantopographie, 2012, http://nbn-resolving.de/nbn:de:hbz:5n-29199). Within this combination method, the full spectral range of the observations is considered. Further, it allows the direct determination of the normal equations (i.e., the inverse of the error covariance matrix) of the mean dynamic topography on arbitrary grids, which is one of the requirements for ocean data assimilation. In this paper, we report progress through selection and improved processing of altimetric data sets. We focus on the preprocessing steps of along-track altimetry data from Jason-1 and Envisat to obtain a mean sea surface profile. During this procedure, a rigorous variance propagation is accomplished, so that, for the first time, the full covariance matrix of the mean sea surface is available. The combination of the mean profile and a combined GRACE/GOCE gravity field model yields a mean dynamic topography model for the North Atlantic Ocean that is characterized by a defined set of assumptions. We show that including the geodetically derived mean dynamic topography with the full error structure in a 3D stationary inverse ocean model improves modeled oceanographic features over previous estimates.

Geodetic and satellite derived ice surface velocities of Vestfonna ice cap

During 2007 we launched a geodetic campaign on the Svalbard ice cap Vestfonna in order to estimate the velocity field of the ice cap. This was done within the frame of the IPY project KINNVIKA. We present here the velocity measurements derived from our campaigns 2007-2010 and compare the geodetic measurements against InSAR velocity fields from satellite platforms from 1995/96 and 2008. We find the spatial distribution of ice speeds from the InSAR is in good agreement within the uncertainty limits with our geodetic measurements. We observe no clear indication of seasonal ice speed differences, but we find a speed-up of the outlet glacier Franklinbreen between the InSAR campaigns, and speculate the outlet is having a surge phase.

Aircraft and satellite derived arctic BrO profiles of the tropospheric column

We derive tropospheric column BrO during the ARCTAS and ARCPAC field campaigns in spring 2008 using retrievals of total column BrO from the satellite UV nadir sensors OMI and GOME-2 using a radiative transfer model and stratospheric column BrO from a photochemical simulation. We conduct a comprehensive comparison of satellite-derived tropospheric BrO column to aircraft in-situ observations of BrO and related species. The aircraft profiles reveal that tropospheric BrO, when present during April 2008, was distributed over a broad range of altitudes rather than being confined to the planetary boundary layer (PBL). Perturbations to the total column resulting from tropospheric BrO are the same magnitude as perturbations due to longitudinal variations in the stratospheric component, so proper accounting of the stratospheric signal is essential for accurate determination of satellite-derived tropospheric BrO. We find reasonably good agreement between satellite-derived tropospheric BrO and columns found using aircraft in-situ BrO profiles, particularly when satellite radiances were obtained over bright surfaces (albedo >0.7), for solar zenith angle <80° and clear sky conditions. The rapid activation of BrO due to surface processes (the bromine explosion) is apparent in both the OMI and GOME-2 based tropospheric columns. The wide orbital swath of OMI allows examination of the evolution of tropospheric BrO on about hourly time intervals near the pole. Low surface pressure, strong wind, and high PBL height are associated with an observed BrO activation event, supporting the notion of bromine activation by high winds over snow.

(Table 1) Date, deployment location, tag longevity, and satellite data of humpback whales (Megaptera novaeangliae) from the Cook Islands

Humpback whales (Megaptera novaeangliae) undertake extensive seasonal migrations from summer feeding areas in high latitudes to winter mating and calving grounds in tropical waters (Clapham and Mead 1999, http://www.jstor.org/stable/3504352). In the Southern Hemisphere, seven populations are recognized by the International Whaling Commission (IWC). In this study, we report the movements of seven whales satellite-tagged in the Cook Islands, including the first documented migration to an antarctic feeding ground. In September 2006 and 2007 we attached Argos satellite-monitored tags to eight humpback whales of various sex and behavioral classes. All whales were tagged in the nearshore waters of Rarotonga (the largest island in the Cooks group).

(Table 1) Date, deployment location, tag longevity, and satellite data of humpback whales (Megaptera novaeangliae) off New Caledonia

Knowledge of the local and migratory movements of humpback whales (Megaptera novaeangliae) from New Caledonia is very limited. To investigate this topic, we attached satellite-monitored tags to 12 whales off southern New Caledonia. Tag longevity ranged from 1 to 52 days (X = 22.5 days). Tagged whales generally moved to the south or southeast, with several spending time in a previously unknown seamount habitat named Antigonia before resuming movement, generally toward Norfolk Island or New Zealand. However, 1 female with a calf traveled the entire length of the western coast of New Caledonia (~450 km) and then west in the direction of the Chesterfield Reefs, a 19th century American (''Yankee'') whaling ground. None of the New Caledonia whales traveled to or toward eastern Australia, which is broadly consistent with the low rate of interchange observed from photo-identification comparisons between these 2 areas. The connections between New Caledonia and New Zealand, together with the relatively low numbers of whales seen in these places generally, support the idea that whales from these 2 areas constitute a single population that remains small and unrecovered.

Planktic foraminifera analyses of surface sediments and satellite data in the Canary Islands region

The Canary Islands region occupies a key position with respect to biogeochemical cycles, with the zonal transition from oligotrophic to nutrient-rich waters and the contribution of Saharan dust to the particle flux. We present the distribution of geochemical proxies (TOC, carbonate, d15N, d13Corg, C/N-ratio) and micropaleontological parameters (diatoms, dinoflagellates, foraminifera, pteropods), in 80 surface-sediment samples in order to characterise the influence of coastally upwelled water on the domain of the subtropical gyre. Results of the surface-sediment analyses confirmed the high biomass gradient from the coast to the open ocean inferred from satellite data of surface chlorophyll or SST. The distribution of total dinoflagellate cysts, the planktic foraminifera species Globigerina bulloides, the diatom resting spore Chaetoceros spp., and TOC concentration coincided well with the areas of strong filament production off Cape Ghir and Cape Yubi. The warm-water planktic foraminifera Globigerinoides ruber (white), the diatom Nitzschia spp., and the d15N-values showed the opposite trend with high values in the open ocean. Factor analyses on the planktic foraminifera species distribution indicated three major assemblages in the Canary Islands region that represent the present surface-water conditions from the upwelling influenced region via a mixing area towards the subtropical gyre.

Atmospheric, oceanic and hydrological polar motion excitation mechanisms based on a combination of geometric and gravimetric space observations

The goal of our study is to determine accurate time series of geophysical Earth rotation excitations to learn more about global dynamic processes in the Earth system. For this purpose, we developed an adjustment model which allows to combine precise observations from space geodetic observation systems, such as Satellite Laser Ranging (SLR), Global Navigation Satellite Systems (GNSS), Very Long Baseline Interferometry (VLBI), Doppler Orbit determination and Radiopositioning Integrated on Satellite (DORIS), satellite altimetry and satellite gravimetry in order to separate geophysical excitation mechanisms of Earth rotation. Three polar motion time series are applied to derive the polar motion excitation functions (integral effect). Furthermore we use five time variable gravity field solutions from Gravity Recovery and Climate Experiment (GRACE) to determine not only the integral mass effect but also the oceanic and hydrological mass effects by applying suitable filter techniques and a land-ocean mask. For comparison the integral mass effect is also derived from degree 2 potential coefficients that are estimated from SLR observations. The oceanic mass effect is also determined from sea level anomalies observed by satellite altimetry by reducing the steric sea level anomalies derived from temperature and salinity fields of the oceans. Due to the combination of all geodetic estimated excitations the weaknesses of the individual processing strategies can be reduced and the technique-specific strengths can be accounted for. The formal errors of the adjusted geodetic solutions are smaller than the RMS differences of the geophysical model solutions. The improved excitation time series can be used to improve the geophysical modeling.