Bayesian retrieval of scatterometer wind fields

The retrieval of wind fields from scatterometer observations has traditionally been separated into two phases; local wind vector retrieval and ambiguity removal. Operationally, a forward model relating wind vector to backscatter is inverted, typically using look up tables, to retrieve up to four local wind vector solutions. A heuristic procedure, using numerical weather prediction forecast wind vectors and, often, some neighbourhood comparison is then used to select the correct solution. In this paper we develop a Bayesian method for wind field retrieval, and show how a direct local inverse model, relating backscatter to wind vector, improves the wind vector retrieval accuracy. We compare these results with the operational U.K. Meteorological Office retrievals, our own CMOD4 retrievals and a neural network based local forward model retrieval. We suggest that the neural network based inverse model, which is extremely fast to use, improves upon current forward models when used in a variational data assimilation scheme.

Non-zero mean Gaussian process prior wind field models

This report outlines the derivation and application of a non-zero mean, polynomial-exponential covariance function based Gaussian process which forms the prior wind field model used in 'autonomous' disambiguation. It is principally used since the non-zero mean permits the computation of realistic local wind vector prior probabilities which are required when applying the scaled-likelihood trick, as the marginals of the full wind field prior. As the full prior is multi-variate normal, these marginals are very simple to compute.

Bayesian inference for wind field retrieval

In many problems in spatial statistics it is necessary to infer a global problem solution by combining local models. A principled approach to this problem is to develop a global probabilistic model for the relationships between local variables and to use this as the prior in a Bayesian inference procedure. We show how a Gaussian process with hyper-parameters estimated from Numerical Weather Prediction Models yields meteorologically convincing wind fields. We use neural networks to make local estimates of wind vector probabilities. The resulting inference problem cannot be solved analytically, but Markov Chain Monte Carlo methods allow us to retrieve accurate wind fields.

Modelling frontal discontinuities in wind fields

A Bayesian procedure for the retrieval of wind vectors over the ocean using satellite borne scatterometers requires realistic prior near-surface wind field models over the oceans. We have implemented carefully chosen vector Gaussian Process models; however in some cases these models are too smooth to reproduce real atmospheric features, such as fronts. At the scale of the scatterometer observations, fronts appear as discontinuities in wind direction. Due to the nature of the retrieval problem a simple discontinuity model is not feasible, and hence we have developed a constrained discontinuity vector Gaussian Process model which ensures realistic fronts. We describe the generative model and show how to compute the data likelihood given the model. We show the results of inference using the model with Markov Chain Monte Carlo methods on both synthetic and real data.

Bayesian inference for wind field retrieval

In many problems in spatial statistics it is necessary to infer a global problem solution by combining local models. A principled approach to this problem is to develop a global probabilistic model for the relationships between local variables and to use this as the prior in a Bayesian inference procedure. We show how a Gaussian process with hyper-parameters estimated from Numerical Weather Prediction Models yields meteorologically convincing wind fields. We use neural networks to make local estimates of wind vector probabilities. The resulting inference problem cannot be solved analytically, but Markov Chain Monte Carlo methods allow us to retrieve accurate wind fields.

Variability of the permetric response in normals and in glaucoma

This study investigated the variability of response associated with various perimetric techniques, with the aim of improving the clinical interpretation of automated static threshold perirnetry. Evaluation of a third generation of perimetric threshold algorithms (SITA) demonstrated a reduction in test duration by approximately 50% both in normal subjects and in glaucoma patients. SITA produced a slightly higher, but clinically insignificant, Mean Sensitivity than with the previous generations of algorithms. This was associated with a decreased between-subject variability in sensitivity and hence, lower confidence intervals for normality. In glaucoma, the SITA algorithms gave rise to more statistically significant visual field defects and a similar between-visit repeatability to the Full Threshold and FASTPAC algorithms. The higher estimated sensitivity observed with SITA compared to Full Threshold and FASTPAC were not attributed to a reduction in the fatigue effect. The investigation of a novel method of maintaining patient fixation, a roving fixation target which paused immediately prior lo the stimulus presentation, revealed a greater degree of fixational instability with the roving fixation target compared to the conventional static fixation target. Previous experience with traditional white-white perimetry did not eradicate the learning effect in short-wavelength automated perimetry (SWAP) in a group of ocular hypertensive patients. The learning effect was smaller in an experienced group of patients compared to a naive group of patients, but was still at a significant level to require that patients should undertake a series of at least three familiarisation tests with SWAP.

A realtime Internet-based quality control application for hurricane surface wind observations

This thesis chronicles the design and implementation of a Internet/Intranet and database based application for the quality control of hurricane surface wind observations. A quality control session consists of selecting desired observation types to be viewed and determining a storm track based time window for viewing the data. All observations of the selected types are then plotted in a storm relative view for the chosen time window and geography is positioned for the storm-center time about which an objective analysis can be performed. Users then make decisions about data validity through visual nearest-neighbor comparison and inspection. The project employed an Object Oriented iterative development method from beginning to end and its implementation primarily features the Java programming language. ^

Spatial, geomorphological, and seasonal variability of CDOM in estuaries of the Florida Coastal Everglades

This paper demonstrates the usefulness of fluorescence techniques for long-term monitoring and assessment of the dynamics (sources, transport and fate) of chromophoric dissolved organic matter (CDOM) in highly compartmentalized estuarine regions with non-point water sources. Water samples were collected monthly from a total of 73 sampling stations in the Florida Coastal Everglades (FCE) estuaries during 2001 and 2002. Spatial and seasonal variability of CDOM characteristics were investigated for geomorphologically distinct sub-regions within Florida Bay (FB), the Ten Thousand Islands (TTI), and Whitewater Bay (WWB). These variations were observed in both quantity and quality of CDOM. TOC concentrations in the FCE estuaries were generally higher during the wet season (June–October), reflecting high freshwater loadings from the Everglades in TTI, and a high primary productivity of marine biomass in FB. Fluorescence parameters suggested that the CDOM in FB is mainly of marine/microbial origin, while for TTI and WWB a terrestrial origin from Everglades marsh plants and mangroves was evident. Variations in CDOM quality seemed mainly controlled by tidal exchange/mixing of Everglades freshwater with Florida Shelf waters, tidally controlled releases of CDOM from fringe mangroves, primary productivity of marine vegetation in FB and diagenetic processes such as photodegradation (particularly for WWB). The source and dynamics of CDOM in these subtropical estuaries is complex and found to be influenced by many factors including hydrology, geomorphology, vegetation cover, landuse and biogeochemical processes. Simple, easy to measure, high sample throughput fluorescence parameters for surface waters can add valuable information on CDOM dynamics to long-term water quality studies which can not be obtained from quantitative determinations alone.

Spatial and Temporal Variability of Dissolved Organic Matter Quantity and Composition in an Oilgotrophic Subtropical Coastal Wetland

Dissolved organic matter (DOM) is an essential component of the carbon cycle and a critical driver in controlling variety of biogeochemical and ecological processes in wetlands. The quality of this DOM as it relates to composition and reactivity is directly related to its sources and may vary on temporal and spatial scales. However, large scale, long-term studies of DOM dynamics in wetlands are still scarce in the literature. Here we present a multi-year DOM characterization study for monthly surface water samples collected at 14 sampling stations along two transects within the greater Everglades, a subtropical, oligotrophic, coastal freshwater wetland-mangrove-estuarine ecosystem. In an attempt to assess quantitative and qualitative variations of DOM on both spatial and temporal scales, we determined dissolved organic carbon (DOC) values and DOM optical properties, respectively. DOM quality was assessed using, excitation emission matrix (EEM) fluorescence coupled with parallel factor analysis (PARAFAC). Variations of the PARAFAC components abundance and composition were clearly observed on spatial and seasonal scales. Dry versus wet season DOC concentrations were affected by dry-down and re-wetting processes in the freshwater marshes, while DOM compositional features were controlled by soil and higher plant versus periphyton sources respectively. Peat-soil based freshwater marsh sites could be clearly differentiated from marl-soil based sites based on EEM–PARAFAC data. Freshwater marsh DOM was enriched in higher plant and soil-derived humic-like compounds, compared to estuarine sites which were more controlled by algae- and microbial-derived inputs. DOM from fringe mangrove sites could be differentiated between tidally influenced sites and sites exposed to long inundation periods. As such coastal estuarine sites were significantly controlled by hydrology, while DOM dynamics in Florida Bay were seasonally driven by both primary productivity and hydrology. This study exemplifies the application of long term optical properties monitoring as an effective technique to investigate DOM dynamics in aquatic ecosystems. The work presented here also serves as a pre-restoration condition dataset for DOM in the context of the Comprehensive Everglades Restoration Plan (CERP).

A realtime internet based quality control application for hurricane surface wind observations

This thesis chronicles the design and implementation of a Intemet/Intranet and database based application for the quality control of hurricane surface wind observations. A quality control session consists of selecting desired observation types to be viewed and determining a storm track based time window for viewing the data. All observations of the selected types are then plotted in a storm relative view for the chosen time window and geography is positioned for the storm-center time about which an objective analysis can be performed. Users then make decisions about data validity through visual nearestneighbor comparison and inspection. The project employed an Object Oriented iterative development method from beginning to end and its implementation primarily features the Java programming language.

Stable isotope ratios on planktonic foraminifer N. pachyderma of surface sediments from the Arctic Ocean (Table 1)

Planktic foraminifers Neogloboquadrina pachyderma (sin.) from 87 eastern and central Arctic Ocean surface sediment samples were analyzed for stable oxygen and carbon isotope composition. Additional results from 52 stations were taken from the literature. The lateral distribution of delta18O (18O/16O) values in the Arctic Ocean reveals a pattern of roughly parallel, W-E stretching zones in the Eurasian Basin, each ~0.5 per mil wide on the delta18O scale. The low horizontal and vertical temperature variability in the Arctic halocline waters (0-100 m) suggests only little influence of temperature on the oxygen isotope distribution of N. pachyderma (sin.). The zone of maximum delta18O values of up to 3.8 per mil is situated in the southern Nansen Basin and relates to the tongue of saline (> 33%.) Atlantic waters entering the Arctic Ocean through the Fram Strait. delta18O values decrease both to the Barents Shelf and to the North Pole, in accordance with the decreasing salinities of the halocline waters. In the Nansen Basin, a strong N-S delta18O gradient is in contrast with a relatively low salinity change and suggests contributions from different freshwater sources, i.e. salinity reduction from sea ice meltwater in the south and from light isotope waters (meteoric precipitation and river-runoff) in the northern part of the basin. North of the Gakkel Ridge, delta18O and salinity gradients are in good accordance and suggest less influence of sea ice melting processes. The delta13C (13C/12C) values of N. pachyderma (sin.) from Arctic Ocean surface sediment samples are generally high (0.75-0.95 per mil). Lower values in the southern Eurasian Basin appear to be related to the intrusion of Atlantic waters. The high delta13C values are evidence for well ventilated surface waters. Because the perennial Arctic sea ice cover largely prevents atmosphere-ocean gas exchange, ventilation on the seasonally open shelves must be of major importance. Lack of delta13C gradients along the main routes of the ice drift from the Siberian shelves to the Fram Strait suggests that primary production (i.e. CO2 consumption) does probably not change the CO2 budget of the Arctic Ocean significantly.

High-Resolution In Situ Oxygen-Argon Studies of Surface Biological and Physical Processes in the Polar Oceans

The Arctic Ocean and Western Antarctic Peninsula (WAP) are the fastest warming regions on the planet and are undergoing rapid climate and ecosystem changes. Until we can fully resolve the coupling between biological and physical processes we cannot predict how warming will influence carbon cycling and ecosystem function and structure in these sensitive and climactically important regions. My dissertation centers on the use of high-resolution measurements of surface dissolved gases, primarily O2 and Ar, as tracers or physical and biological functioning that we measure underway using an optode and Equilibrator Inlet Mass Spectrometry (EIMS). Total O2 measurements are common throughout the historical and autonomous record but are influenced by biological (net metabolic balance) and physical (temperature, salinity, pressure changes, ice melt/freeze, mixing, bubbles and diffusive gas exchange) processes. We use Ar, an inert gas with similar solubility properties to O2, to devolve distinct records of biological (O2/Ar) and physical (Ar) oxygen. These high-resolution measurements that expose intersystem coupling and submesoscale variability were central to studies in the Arctic Ocean, WAP and open Southern Ocean that make up this dissertation.

Key findings of this work include the documentation of under ice and ice-edge blooms and basin scale net sea ice freeze/melt processes in the Arctic Ocean. In the WAP O2 and pCO2 are both biologically driven and net community production (NCP) variability is controlled by Fe and light availability tied to glacial and sea ice meltwater input. Further, we present a feasibility study that shows the ability to use modeled Ar to derive NCP from total O2 records. This approach has the potential to unlock critical carbon flux estimates from historical and autonomous O2 measurements in the global oceans.

Surface wind stress modeled for the Panama Seaway with the Community Climate System Model 3

The early Pliocene warm phase was characterized by high sea surface temperatures and a deep thermocline in the eastern equatorial Pacific. A new hypothesis suggests that the progressive closure of the Panamanian seaway contributed substantially to the termination of this zonally symmetric state in the equatorial Pacific. According to this hypothesis, intensification of the Atlantic meridional overturning circulation (AMOC) - induced by the closure of the gateway - was the principal cause of equatorial Pacific thermocline shoaling during the Pliocene. In this study, twelve Panama seaway sensitivity experiments from eight ocean/climate models of different complexity are analyzed to examine the effect of an open gateway on AMOC strength and thermocline depth. All models show an eastward Panamanian net throughflow, leading to a reduction in AMOC strength compared to the corresponding closed-Panama case. In those models that do not include a dynamic atmosphere, deepening of the equatorial Pacific thermocline appears to scale almost linearly with the throughflow-induced reduction in AMOC strength. Models with dynamic atmosphere do not follow this simple relation. There are indications that in four out of five models equatorial wind-stress anomalies amplify the tropical Pacific thermocline deepening. In summary, the models provide strong support for the hypothesized relationship between Panama closure and equatorial Pacific thermocline shoaling.