Generalized ocean color inversion model for retrieving marine inherent optical properties

Ocean color measured from satellites provides daily, global estimates of marine inherent optical properties (IOPs). Semi-analytical algorithms (SAAs) provide one mechanism for inverting the color of the water observed by the satellite into IOPs. While numerous SAAs exist, most are similarly constructed and few are appropriately parameterized for all water masses for all seasons. To initiate community-wide discussion of these limitations, NASA organized two workshops that deconstructed SAAs to identify similarities and uniqueness and to progress toward consensus on a unified SAA. This effort resulted in the development of the generalized IOP (GIOP) model software that allows for the construction of different SAAs at runtime by selection from an assortment of model parameterizations. As such, GIOP permits isolation and evaluation of specific modeling assumptions, construction of SAAs, development of regionally tuned SAAs, and execution of ensemble inversion modeling. Working groups associated with the workshops proposed a preliminary default configuration for GIOP (GIOP-DC), with alternative model parameterizations and features defined for subsequent evaluation. In this paper, we: (1) describe the theoretical basis of GIOP; (2) present GIOP-DC and verify its comparable performance to other popular SAAs using both in situ and synthetic data sets; and, (3) quantify the sensitivities of their output to their parameterization. We use the latter to develop a hierarchical sensitivity of SAAs to various model parameterizations, to identify components of SAAs that merit focus in future research, and to provide material for discussion on algorithm uncertainties and future emsemble applications.

The Ocean Colour Climate Change Initiative: III. A round-robin comparison on in-water bio-optical algorithms

Satellite-derived remote-sensing reflectance (Rrs) can be used for mapping biogeochemically relevant variables, such as the chlorophyll concentration and the Inherent Optical Properties (IOPs) of the water, at global scale for use in climate-change studies. Prior to generating such products, suitable algorithms have to be selected that are appropriate for the purpose. Algorithm selection needs to account for both qualitative and quantitative requirements. In this paper we develop an objective methodology designed to rank the quantitative performance of a suite of bio-optical models. The objective classification is applied using the NASA bio-Optical Marine Algorithm Dataset (NOMAD). Using in situRrs as input to the models, the performance of eleven semi-analytical models, as well as five empirical chlorophyll algorithms and an empirical diffuse attenuation coefficient algorithm, is ranked for spectrally-resolved IOPs, chlorophyll concentration and the diffuse attenuation coefficient at 489 nm. The sensitivity of the objective classification and the uncertainty in the ranking are tested using a Monte-Carlo approach (bootstrapping). Results indicate that the performance of the semi-analytical models varies depending on the product and wavelength of interest. For chlorophyll retrieval, empirical algorithms perform better than semi-analytical models, in general. The performance of these empirical models reflects either their immunity to scale errors or instrument noise in Rrs data, or simply that the data used for model parameterisation were not independent of NOMAD. Nonetheless, uncertainty in the classification suggests that the performance of some semi-analytical algorithms at retrieving chlorophyll is comparable with the empirical algorithms. For phytoplankton absorption at 443 nm, some semi-analytical models also perform with similar accuracy to an empirical model. We discuss the potential biases, limitations and uncertainty in the approach, as well as additional qualitative considerations for algorithm selection for climate-change studies. Our classification has the potential to be routinely implemented, such that the performance of emerging algorithms can be compared with existing algorithms as they become available. In the long-term, such an approach will further aid algorithm development for ocean-colour studies.

Analytical phytoplankton carbon measurements spanning diverse ecosystems

The measurement of phytoplankton carbon (Cphyto) in the field has been a long-sought but elusive goal in oceanography. Proxy measurements of Cphyto have been employed in the past, but are subject to many confounding influences that undermine their accuracy. Here we report the first directly measured Cphyto values from the open ocean. The Cphyto samples were collected from a diversity of environments, ranging from Pacific and Atlantic oligotrophic gyres to equatorial upwelling systems to temperate spring conditions. When compared to earlier proxies, direct measurements of Cphyto exhibit the strongest relationship with particulate backscattering coefficients (bbp) (R2=0.69). Chlorophyll concentration and total particulate organic carbon (POC) concentration accounted for ~20% less variability in Cphyto than bbp. Ratios of Cphyto to Chl a span an order of magnitude moving across and within distinct ecosystems. Similarly, Cphyto:POC ratios were variable with the lowest values coming from productive temperate waters and the highest from oligotrophic gyres. A strong relationship between Cphyto and bbp is particularly significant because bbp is a property retrievable from satellite ocean color measurements. Our results, therefore, are highly encouraging for the global monitoring of phytoplankton biomass from space. The continued application of our Cphyto measurement approach will enable validation of satellite retrievals and contribute to an improved understanding of environmental controls on phytoplankton biomass and physiology.

A method for semi-automated objective quantification of linear bedforms from multi-scale Digital Elevation Models

The increasing availability of large, detailed digital representations of the Earth’s surface demands the application of objective and quantitative analyses. Given recent advances in the understanding of the mechanisms of formation of linear bedform features from a range of environments, objective measurement of their wavelength, orientation, crest and trough positions, height and asymmetry is highly desirable. These parameters are also of use when determining observation-based parameters for use in many applications such as numerical modelling, surface classification and sediment transport pathway analysis. Here, we (i) adapt and extend extant techniques to provide a suite of semi-automatic tools which calculate crest orientation, wavelength, height, asymmetry direction and asymmetry ratios of bedforms, and then (ii) undertake sensitivity tests on synthetic data, increasingly complex seabeds and a very large-scale (39 000km2) aeolian dune system. The automated results are compared with traditional, manually derived,measurements at each stage. This new approach successfully analyses different types of topographic data (from aeolian and marine environments) from a range of sources, with tens of millions of data points being processed in a semi-automated and objective manner within minutes rather than hours or days. The results from these analyses show there is significant variability in all measurable parameters in what might otherwise be considered uniform bedform fields. For example, the dunes of the Rub’ al Khali on the Arabian peninsula are shown to exhibit deviations in dimensions from global trends. Morphological and dune asymmetry analysis of the Rub’ al Khali suggests parts of the sand sea may be adjusting to a changed wind regime from that during their formation 100 to 10 ka BP.

Coastal upwelling and downwelling forcing of circulation in a semi-enclosed bay: Ria de Vigo

Semi-enclosed bays in upwelling regions are exposed to forcing related to winds, currents and buoyancy over the shelf. The influence of this external forcing is moderated by factors such as connectivity to the open ocean, shelter by surrounding topography, dimensions of the bay, and freshwater outflows. Such bays, preferred locations for ports, mariculture, marine industry, recreational activities and coastal settlement, present a range of characteristics, understanding of which is necessary to their rational management. Observations in such a semi-enclosed bay, the Ria de Vigo in Spain, are used to characterize the influence of upwelling and downwelling pulses on its circulation. In this location, near the northern limit of the Iberian upwelling system, upwelling events dominate during a short summer season and downwelling events the rest of the year. The ria response to the external forcing is central to nutrient supply and resultant plankton productivity that supports its high level of cultured mussel production. Intensive field studies in September 2006 and June 2007 captured a downwelling event and an upwelling event, respectively. Data from eight current profiler moorings and boat-based MiniBat/ADCP surveys provided an unprecedented quasi-synoptic view of the distribution of water masses and circulation patterns in any ria. In the outer ria, circulation was dominated by the introduction of wind-driven alongshore flow from the external continental shelf through the ria entrances and its interaction with the topography. In the middle ria, circulation was primarily related to the upwelling/downwelling cycle, with a cool, salty and dense lower layer penetrating to the inner ria during upwelling over the shelf. A warmer, lower salinity and less dense surface layer of coastal waters flowed inward during downwelling. Without external forcing, the inner ria responded primarily to tides and buoyancy changes related to land runoff. Under both upwelling and downwelling conditions, the flushing of the ria involved shelf responses to wind pulses. Their persistence for a few days was sufficient to allow waters from the continental shelf to penetrate the innermost ria. Longer term observations supported by numerical modeling are required to confirm the generality of such flushing events in the ria and determine their typical frequency, while comparative studies should explore how these scenarios fit into the range of conditions experienced in other semi-enclosed bays.