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.

Bathymetric map of Heron Reef, Australia, derived from airborne hyperspectral data at 1 m resolution

A simple method for efficient inversion of arbitrary radiative transfer models for image analysis is presented. The method operates by representing the shape of the function that maps model parameters to spectral reflectance by an adaptive look-up tree (ALUT) that evenly distributes the discretization error of tabulated reflectances in spectral space. A post-processing step organizes the data into a binary space partitioning tree that facilitates an efficient inversion search algorithm. In an example shallow water remote sensing application, the method performs faster than an implementation of previously published methodology and has the same accuracy in bathymetric retrievals. The method has no user configuration parameters requiring expert knowledge and minimizes the number of forward model runs required, making it highly suitable for routine operational implementation of image analysis methods. For the research community, straightforward and robust inversion allows research to focus on improving the radiative transfer models themselves without the added complication of devising an inversion strategy.

Hyperspectral downwelling irradiance during POLARSTERN cruise ANT-XXVI/1

A Monte Carlo based radiative transfer model has been developed for calculating the availability of solar radiation within the top 100 m of the ocean. The model is optimized for simulations of spatial high resolution downwelling irradiance Ed fluctuations that arise from the lensing effect of waves at the water surface. In a first step the accuracy of simulation results has been verified by measurements of the oceanic underwater light field and through intercomparison with an established radiative transfer model. Secondly the potential depth-impact of nonlinear shaped single waves, from capillary to swell waves, is assessed by considering the most favorable conditions for light focusing, i.e. monochromatic light at 490 nm, very clear oceanic water with a low chlorophyll a content of 0.1 mg/m**3 and high sun elevation. Finally light fields below irregular wave profiles accounting for realistic sea states were simulated. Our simulation results suggest that under open ocean conditions light flashes with 50% irradiance enhancements can appear down to 35 m depth, and light variability in the range of ±10% compared to the mean Ed is still possible in 100 m depth.

The CoastColour Round Robin datasets: a database to evaluate the performance of algorithms for the retrieval of water quality parameters in coastal waters

The CoastColour project Round Robin (CCRR) project (http://www.coastcolour.org) funded by the European Space Agency (ESA) was designed to bring together a variety of reference datasets and to use these to test algorithms and assess their accuracy for retrieving water quality parameters. This information was then developed to help end-users of remote sensing products to select the most accurate algorithms for their coastal region. To facilitate this, an inter-comparison of the performance of algorithms for the retrieval of in-water properties over coastal waters was carried out. The comparison used three types of datasets on which ocean colour algorithms were tested. The description and comparison of the three datasets are the focus of this paper, and include the Medium Resolution Imaging Spectrometer (MERIS) Level 2 match-ups, in situ reflectance measurements and data generated by a radiative transfer model (HydroLight). The datasets mainly consisted of 6,484 marine reflectance associated with various geometrical (sensor viewing and solar angles) and sky conditions and water constituents: Total Suspended Matter (TSM) and Chlorophyll-a (CHL) concentrations, and the absorption of Coloured Dissolved Organic Matter (CDOM). Inherent optical properties were also provided in the simulated datasets (5,000 simulations) and from 3,054 match-up locations. The distributions of reflectance at selected MERIS bands and band ratios, CHL and TSM as a function of reflectance, from the three datasets are compared. Match-up and in situ sites where deviations occur are identified. The distribution of the three reflectance datasets are also compared to the simulated and in situ reflectances used previously by the International Ocean Colour Coordinating Group (IOCCG, 2006) for algorithm testing, showing a clear extension of the CCRR data which covers more turbid waters.

Airborne hyperspectral image data of Heron Reef, Australia

This airborne hyperspectral (19 bands) image data of Heron Reef, Great Barrier Reef, Australia is derived from Compact Airborne Spectrographic Imager (CASI) data acquired on 1st and 3rd of July 2002, latitude -23.45, longitude 151.92. Processing and correction to at-surface data was completed by Karen Joyce (Joyce, 2004). Raw imagery consisted several images corresponding to the number of flight paths taken to cover the entire Heron Reef. Spatial resolution is one meter. Radiometric corrections converted the at-sensor digital number values to at surface spectral radiance values using sensor specific calibration coefficients and CSIRO's c-WomBat-c atmospheric correction software. Geometric corrections were done using field collected coordinates of features identified in the image. Projection used was Universal Transverse Mercator Zone 56 South and Datum used was WGS 84. Image data is in TIFF format.

Simulated Bromoform emission and concentration using the ocean biogeochemistry model MPIOM/HAMOCC forced by 6-hourly NCEP data

Bromoform (CHBr3) is one important precursor of atmospheric reactive bromine species that are involved in ozone depletion in the troposphere and stratosphere. In the open ocean bromoform production is linked to phytoplankton that contains the enzyme bromoperoxidase. Coastal sources of bromoform are higher than open ocean sources. However, open ocean emissions are important because the transfer of tracers into higher altitude in the air, i.e. into the ozone layer, strongly depends on the location of emissions. For example, emissions in the tropics are more rapidly transported into the upper atmosphere than emissions from higher latitudes. Global spatio-temporal features of bromoform emissions are poorly constrained. Here, a global three-dimensional ocean biogeochemistry model (MPIOM-HAMOCC) is used to simulate bromoform cycling in the ocean and emissions into the atmosphere using recently published data of global atmospheric concentrations (Ziska et al., 2013) as upper boundary conditions. Our simulated surface concentrations of CHBr3 match the observations well. Simulated global annual emissions based on monthly mean model output are lower than previous estimates, including the estimate by Ziska et al. (2013), because the gas exchange reverses when less bromoform is produced in non-blooming seasons. This is the case for higher latitudes, i.e. the polar regions and northern North Atlantic. Further model experiments show that future model studies may need to distinguish different bromoform-producing phytoplankton species and reveal that the transport of CHBr3 from the coast considerably alters open ocean bromoform concentrations, in particular in the northern sub-polar and polar regions.

(Table 2) Age model of ODP Site 184-1146

Orbital tuning of benthic d18O is a common approach for assigning ages to ocean sediment records. Similar environmental forcing of the northern South China Sea and the southeast Asian cave regions allows for transfer of the speleothem d18O radiometric chronology to the planktonic and benthic d18O records from Ocean Drilling Program Site 1146, yielding a new chronology with 41 radiometrically calibrated datums, spanning the past 350 kyr. This approach also provides for an independent assessment of the accuracy of the orbitally tuned benthic d18O chronology for the last 350 kyr. The largest differences relative to the latest chronology occur in marine isotope stages (MIS) 5.4, 5.5, 6, 7, and 9.3. Prominent suborbital-scale structure believed to be global in nature is identified within MIS 5.4 and MIS 7.2. On the basis of the radiometrically calibrated chronology, the time constant of the ice sheet is found to be 5.4 kyr at the precession band (light d18O lags precession minima by -55.4°) and 10.4 kyr at the obliquity band (light d18O lags obliquity maxima by 57.4°). These values are significantly shorter than the single 17 kyr time constant originally estimated by Imbrie et al. (1984), based primarily on the timing of terminations I and II and the 15 kyr time constant used by Lisiecki and Raymo (2005, doi:10.1029/2004PA001071).