A nonlinear optimization algorithm for WindSat wind vector retrievals

WindSat is a space-based polarimetric microwave radiometer designed to demonstrate the capability to measure the ocean surface wind vector using a radiometer. We describe a nonlinear iterative algorithm for simultaneous retrieval of sea surface temperature, columnar water vapor, columnar cloud liquid water, and the ocean surface wind vector from WindSat measurements. The algorithm uses a physically based forward model function for the WindSat brightness temperatures. Empirical corrections to the physically based model are discussed. We present evaluations of initial retrieval performance using a six-month dataset of WindSat measurements and collocated data from other satellites and a numerical weather model. We focus primarily on the application to wind vector retrievals.     

An emissivity-based wind vector retrieval algorithm for the WindSat polarimetric radiometer

The Naval Research Laboratory WindSat polarimetric radiometer was launched on January 6, 2003 and is the first fully polarimetric radiometer to be flown in space. WindSat has three fully polarimetric channels at 10.7, 18.7, and 37.0 GHz and vertically and horizontally polarized channels at 6.8 and 23.8 GHz. A first-generation wind vector retrieval algorithm for the WindSat polarimetric radiometer is developed in this study. An atmospheric clearing algorithm is used to estimate the surface emissivity from the measured WindSat brightness temperature at each channel. A specular correction factor is introduced in the radiative transfer equation to account for excess reflected atmospheric brightness, compared to the specular assumption, as a function wind speed. An empirical geophysical model function relating the surface emissivity to the wind vector is derived using coincident QuikSCAT scatterometer wind vector measurements. The confidence in the derived harmonics for the polarimetric channels is high and should be considered suitable to validate analytical surface scattering models for polarized ocean surface emission. The performance of the retrieval algorithm is assessed with comparisons to Global Data Assimilation System (GDAS) wind vector outputs. The root mean square (RMS) uncertainty of the closest wind direction ambiguity is less than 20/spl deg/ for wind speeds greater than 6 m/s and less than 15/spl deg/ at 10 m/s and greater. The retrieval skill, the percentage of retrievals in which the first-rank solution is the closest to the GDAS reference, is 75% at 7 m/s and 85% or higher above 10 m/s. The wind speed is retrieved with an RMS uncertainty of 1.5 m/s.     

High-resolution passive polarimetric microwave mapping of oceansurface wind vector fields

The retrieval of ocean surface wind fields in both one and two dimensions is demonstrated using passive polarimetric microwave imagery obtained from a conical-scanning airborne polarimeter. The retrieval method is based on an empirical geophysical model function (GMF) for ocean surface thermal emission and an adaptive maximum likelihood (ML) wind vector estimator. Data for the GMF were obtained using the polarimetric scanning radiometer/digital (PSR/D) on the NASA P-3 aircraft during the Labrador Sea Deep Convection Experiment in 1997. To develop the GMF, a number of buoy overflights and GPS dropsondes were used, out of which a GMF of 10.7, 18.7, and 37.0 GHz azimuthal harmonics for the first three Stokes parameters was constructed for the SSM/I incident angle of 53.1°. The data show repeatable azimuthal harmonic coefficient amplitudes of ~2-3 K peak-to-peak, with a 100% increase in harmonic amplitudes as the frequency is increased from 10.7 to 37 GHz. The GMF is consistent with and extends the results of two independent studies of SSM/I data and also provides a model for the third Stokes parameter over wind speeds up to 20 m/s. The aircraft data show that the polarimetric channels are much less susceptible to geophysical noise associated with maritime convection than the first two Stokes parameters. The polarimetric measurement technique used in the PSR/D also demonstrates the viability of digital correlation radiometry for aircraft or satellite measurements of the full Stokes vector. The ML retrieval algorithm incorporates the additional information on wind direction available from multiple looks and polarimetric channels in a straightforward manner and accommodates the reduced SNRs of the first two Stokes parameters in the presence of convection by weighting these channels by their inverse SNR     

Calibration of WindSat polarimetric channels with a vicarious cold reference

Absolute calibration of WindSat's third and fourth Stokes brightness temperatures (T/sub 3/ and T/sub 4/) is needed at the tenth of Kelvin level in order to adequately resolve their dependence on wind direction. Previous aircraft based fully polarimetric microwave radiometers have generally relied on "circle flights", during which a single area of the ocean is observed at all azimuth angles, to estimate residual biases in the calibration of its polarimetric channels. WindSat, the first spaceborne fully polarimetric microwave radiometer, operates in low Earth orbit and thus cannot execute this traditional calibration technique. A new method is presented to estimate the residual biases that are present in WindSat's T/sub 3/ and T/sub 4/ estimates. The method uses a vicarious cold reference brightness temperature applied to measurements made by WindSat at /spl plusmn/45/spl deg/ slant linear (T/sub P/ and T/sub M/) and left- and right-hand circular (T/sub L/ and T/sub R/) polarization. WindSat derives the third and fourth Stokes brightness temperatures by the differences T/sub P/-T/sub M/ and T/sub L/-T/sub R/, respectively. The method is demonstrated by applying it to the 10.7-GHz WindSat observations. Calibration biases of 0.2-0.6 K are determined with a precision of 0.04 K.     

Polarimetric measurements of sea surface brightness temperaturesusing an aircraft K-band radiometer

Presents the first experimental evidence that the polarimetric brightness temperatures of sea surfaces are sensitive to ocean wind direction in the incidence angle range of 30 to 50°. The experimental data were collected by a K-band (19.35 GHz) polarimetric wind radiometer (WINDRAD) mounted on the NASA DC-8 aircraft. A set of aircraft radiometer flights was successfully completed in November 1993. The authors performed circle flights over National Data Buoy Center (NDBC) moored buoys deployed off the northern California coast, which provided ocean wind measurements. The results indicate that passive polarimetric radiometry has a strong potential for global ocean wind speed and direction measurements from space     

Effects of foam on ocean surface microwave emission inferred from radiometric observations of reproducible breaking waves

WindSat, the first satellite polarimetric microwave radiometer, and the NPOESS Conical Microwave Imager/Sounder both have as a key objective the retrieval of the ocean surface wind vector from radiometric brightness temperatures. Available observations and models to date show that the wind direction signal is only 1-3 K peak-to-peak at 19 and 37 GHz, much smaller than the wind speed signal. In order to obtain sufficient accuracy for reliable wind direction retrieval, uncertainties in geophysical modeling of the sea surface emission on the order of 0.2 K need to be removed. The surface roughness spectrum has been addressed by many studies, but the azimuthal signature of the microwave emission from breaking waves and foam has not been adequately addressed. Recently, a number of experiments have been conducted to quantify the increase in sea surface microwave emission due to foam. Measurements from the Floating Instrumentation Platform indicated that the increase in ocean surface emission due to breaking waves may depend on the incidence and azimuth angles of observation. The need to quantify this dependence motivated systematic measurement of the microwave emission from reproducible breaking waves as a function of incidence and azimuth angles. A number of empirical parameterizations of whitecap coverage with wind speed were used to estimate the increase in brightness temperatures measured by a satellite microwave radiometer due to wave breaking in the field of view. These results provide the first empirically based parameterization with wind speed of the effect of breaking waves and foam on satellite brightness temperatures at 10.8, 19, and 37 GHz.     

A Comparison of the Impact of QuikScat and WindSat Wind Vector Products on Met Office Analyses and Forecasts

Several studies have demonstrated that retrievals of wind vectors from the WindSat polarimetric radiometer are of sufficient quality to be considered for assimilation in operational numerical weather prediction models. In this paper, WindSat data are used in a state-of-the-art global meteorological analysis and forecasting system. Each wind vector contains a directional ambiguity and so is assimilated in a similar way to that of scatterometer data. The forecast impact of using analyses containing information from WindSat data was investigated for a period during August and September of 2005, when a large number of tropical cyclones were present. Forecast errors were reduced in the surface pressure fields, and the average improvement across the forecast range was found to be 1.0%. This is comparable to the improvement of 1.1% found in the same fields when winds were assimilated from the QuikScat scatterometer. The impact on tropical cyclone tracks in the forecasts was also studied. The scatterometer improved (reduced) the track errors markedly by 25% in the analyses. When impacts across the forecast range out to five days were also included, the improvement was found to be 8%. In contrast, the assimilation of WindSat data improved the analysis track errors by 7%, although this figure was found to be 10% across the complete forecast range.      

Deep-space calibration of the WindSat radiometer

The WindSat microwave polarimetric radiometer consists of 22 channels of polarized brightness temperatures operating at five frequencies: 6.8, 10.7, 18.7, 23.8, and 37.0 GHz. The 10.7-, 18.7-, and 37.0-GHz channels are fully polarimetric (vertical/horizontal, /spl plusmn/45/spl deg/ and left-hand and right-hand circularly polarized) to measure the four Stokes radiometric parameters. The principal objective of this Naval Research Laboratory experiment, which flys on the USAF Coriolis satellite, is to provide the proof of concept of the first passive measurement of ocean surface wind vector from space. This paper presents details of the on-orbit absolute radiometric calibration procedure, which was performed during of a series of satellite pitch maneuvers. During these special tests, the satellite pitch was slowly ramped to +45/spl deg/ (and -45/spl deg/), which caused the WindSat conical spinning antenna to view deep space during the forward (or aft portion) of the azimuth scan. When viewing the homogeneous and isotropic brightness of space (2.73 K) through both the main reflector and the cold-load calibration reflector, it is possible to determine the absolute calibration of the individual channels and the relative calibration bias between polarimetric channels. Results demonstrate consistent and stable channel calibrations (with very small brightness biases) that exceed the mission radiometric calibration requirements.     

Polarimetric brightness temperatures of sea surfaces measured withaircraft K- and Ka-band radiometers

Dual-frequency (19 and 37 GHz), multi-incidence measurements of the Stokes parameters of sea surface microwave emission are reported. A series of aircraft polarimetric radiometer flights were carried out over the National Data Buoy Center (NDBC) moored buoys deployed off the northern California coast in July and August 1994. Measured radiometric temperatures showed a few Kelvin azimuth modulations in all Stokes parameters with respect to the wind direction. Wind directional signals observed in the 37-GHz channel were similar to those in the 19-GHz channel. This indicates that the wind direction signals in sea surface brightness temperatures have a weak frequency dependence in the range of 19-37 GHz. Harmonic coefficients of the wind direction signals were derived from experimental data versus incidence angle. It was found that the first harmonic coefficients, which are caused by the up and downwind asymmetric surface features, had a small increasing trend with the incidence angle. In contrast, the second harmonic coefficients, caused by the up and crosswind asymmetry, showed significant variations in T _v and U data, with a sign change when the incidence angle increased from 45° to 65°. Besides the first three Stokes parameters, the fourth Stokes parameter, V, which had never been measured before for sea surfaces, was measured using our 19-GHz channel. The Stokes parameter V. Has an odd symmetry just like that of the third Stokes parameter U, and increases with increasing incidence angles. In summary, sea surface features created by surface winds are anisotropic in azimuth direction and modulate all Stokes parameters of sea surface microwave brightness temperatures by as large as a few Kelvin in the range of incidence angles from 45° to 65° applicable to spaceborne observations     

The WindSat spaceborne polarimetric microwave radiometer: sensor description and early orbit performance

The global ocean surface wind vector is a key parameter for short-term weather forecasting, the issuing of timely weather warnings, and the gathering of general climatological data. In addition, it affects a broad range of naval missions, including strategic ship movement and positioning, aircraft carrier operations, aircraft deployment, effective weapons use, underway replenishment, and littoral operations. WindSat is a satellite-based multifrequency polarimetric microwave radiometer developed by the Naval Research Laboratory for the U.S. Navy and the National Polar-orbiting Operational Environmental Satellite System Integrated Program Office. It is designed to demonstrate the capability of polarimetric microwave radiometry to measure the ocean surface wind vector from space. The sensor provides risk reduction for the development of the Conical Microwave Imager Sounder, which is planned to provide wind vector data operationally starting in 2010. WindSat is the primary payload on the Department of Defense Coriolis satellite, which was launched on January 6, 2003. It is in an 840-km circular sun-synchronous orbit. The WindSat payload is performing well and is currently undergoing rigorous calibration and validation to verify mission success.     

Comparison of WindSat brightness temperatures with two-scale model predictions

Predictions of the polarized microwave brightness temperatures over the ocean are made using a two-scale surface bidirectional reflectance model combined with an atmospheric radiative transfer model. The reflected atmospheric radiation is found to contribute significantly to the magnitude and directional dependence of the brightness temperatures. The predicted brightness temperatures are also sensitive to the form of the shortwave spectrum. Calculations are made using a new physically based model for the wave spectrum, and preliminary comparisons are made with WindSat observations at 10.7, 18.7, and 37 GHz, for wind speeds ranging from 0-20 m/s and for vertically integrated atmospheric water vapor concentrations from 0-70 mm. Predictions of the mean (azimuthally averaged) brightness temperatures for vertical and horizontal polarization agree quite well with WindSat observations over this range of wind speeds and water vapor concentrations. The predicted azimuthal variations of the third and fourth Stokes parameters also agree fairly well with the observations, except for the fourth Stokes parameter at 37 GHz. Further adjustments of the wave spectrum are expected to improve the agreement.     

WindSat Passive Microwave Polarimetric Signatures of the Greenland Ice Sheet

WindSat has systematically collected the first global fully polarimetric passive microwave data over both land and ocean. As the first spaceborne polarimetric microwave radiometer, it was designed to measure ocean surface wind speed and direction by including the third and fourth Stokes parameters, which are mostly related to the asymmetric structures of the ocean surface roughness. Although designed for wind vector retrieval, WindSat data are also collected over land and ice, and this new data has revealed, for the first time, significant land signals in the third and fourth Stokes parameter channels, particularly over Greenland and the Antarctic ice sheets. The third and fourth Stokes parameters show well-defined large azimuth modulations that appear to be correlated with geophysical variations, particularly snow structure, melting, and metamorphism, and have distinct seasonal variation. The polarimetric signatures are relatively weak in the summer and are strongest around spring. This corresponds well with the formation and erosion of the sastrugi in the dry snow zone and snowmelt in the soaked zone. In this paper, we present the full polarimetric signatures obtained from WindSat over Greenland, and use a simple empirical observation model to quantify the azimuthal variations of the signatures in space and time.      

Wind direction over the ocean determined by an airborne, imaging,polarimetric radiometer system

The speed and direction of winds over the ocean can be determined by polarimetric radiometers. This has been established by theoretical work and demonstrated experimentally using airborne radiometers carrying out circle flights and thus measuring the full 360° azimuthal response from the sea surface. An airborne experiment, with the aim of measuring wind direction over the ocean using an imaging polarimetric radiometer, is described. A polarimetric radiometer system of the correlation type, measuring all four Stokes brightness parameters, is used. Imaging is achieved using a 1-m aperture conically scanning antenna. The polarimetric azimuthal signature of the ocean is known from modeling and circle flight experiments. Combining the signature with the measured brightness data from just a single flight track enables the wind direction to be determined on a pixel-by-pixel basis in the radiometer imagery     

HUT fully polarimetric calibration standard for microwave radiometry

This paper describes the Helsinki University of Technology's Fully Polarimetric Calibration Standard (FPCS). The developed standard generates a complete Stokes reference vector and it is applied for the end-to-end absolute calibration of a fully polarimetric microwave radiometer at 36.5 GHz. The FPCS is based on the function principle of a Gasiewski-Kunkee linearly polarized (tripolarimetric) standard, with an additional phase retardation plate to generate the fourth Stokes parameter. Design considerations and operational aspects of the standard are discussed in this paper. An advanced calibration procedure, which takes advantage of both the tripolarimetric and fully polarimetric calibration scenes to suppress calibration uncertainties, is introduced. The feasibility of the standard has been verified and the generated brightness temperatures in a sample calibration are presented. An extensive set of tests has been performed to evaluate the characteristics and performance of the calibration standard. Furthermore, the use of the advanced calibration procedure to measure the characteristics of the phase retardation plate has been successfully demonstrated. The achievable calibration accuracy is analyzed and discussed relative to requirements for maritime wind vector measurements; the results indicate that the pixel-to-pixel retrieval of the wind speed is possible with high accuracy and the retrieval of the wind direction with at least moderate accuracy. In addition to calibration of a fully polarimetric radiometer, other potential applications, e.g., linearity measurements, are discussed.     

An updated analysis of the ocean surface wind direction signal in passive microwave brightness temperatures

We analyze the wind direction signal for vertically (v) and horizontally (h) polarized microwave radiation at 37 GHz, 19 GHz, and 11 GHz; and an Earth incidence angle of 53/spl deg/. We use brightness temperatures from SSM/I and TMI and wind vectors from buoys and the QUIKSCAT scatterometer. The wind vectors are space and time collocated with the radiometer measurements. Water vapor, cloud water and sea surface temperature are obtained from independent measurements and are uncorrelated with the wind direction. We find a wind direction signal that is noticeably smaller at low and moderate wind speeds than a previous analysis had indicated. We attribute the discrepancy to errors in the atmospheric parameters that were present in the data set of the earlier study. We show that the polarization combination 2v-h is almost insensitive to atmospheric changes and agrees with the earlier results. The strength of our new signals agrees well with JPL aircraft radiometer measurements. It is significantly smaller than the prediction of the two-scale sea surface emission model for low and intermediate wind speeds.     

Evaluation of hurricane ocean vector winds from WindSat

The ability to accurately measure ocean surface wind vectors from space in all weather conditions is important in many scientific and operational usages. One highly desirable application of satellite-based wind vector retrievals is to provide realistic estimates of tropical cyclone intensity for hurricane monitoring. Historically, the extreme environmental conditions in tropical cyclones (TCs) have been a challenge to traditional space-based wind vector sensing provided by microwave scatterometers. With the advent of passive microwave polarimetry, an alternate tool for estimating surface wind conditions in the TC has become available. This paper evaluates the WindSat polarimetric radiometer's ability to accurately sense winds within TCs. Three anecdotal cases studies are presented from the 2003 Atlantic Hurricane season. Independent surface wind estimates from aircraft flights and other platforms are used to provide surface wind fields for comparison to WindSat retrievals. Results of a subjective comparison of wind flow patterns are presented as well as quantitative statistics for point location comparisons of wind speed and direction.     

Modeling of wind direction signals in polarimetric sea surfacebrightness temperatures

There has been an increasing interest in the applications of polarimetric microwave radiometers for ocean wind remote sensing. Aircraft and spaceborne radiometers have found a few Kelvins wind direction signals in sea surface brightness temperatures, in addition to their sensitivities to wind speeds. However, it was not clear what physical scattering mechanisms produced the observed brightness dependence on wind direction. To this end, polarimetric microwave emissions from wind-generated sea surfaces are investigated with a polarimetric two-scale scattering model, which relates the directional wind-wave spectrum to passive microwave signatures of sea surfaces. Theoretical azimuthal modulations are found to agree well with experimental observations for all Stokes parameters from near nadir to 65° incidence angles. The upwind and downwind asymmetries of brightness temperatures were interpreted using the hydrodynamic modulation. The contributions of Bragg scattering by short waves, geometric optics scattering by long waves and sea foam are examined. The geometric optics scattering mechanism underestimates the directional signals in the first three Stokes parameters, and predicts no signals in the fourth Stokes parameter (V). In contrast, the Bragg scattering was found to dominate the wind direction signals from the two-scale model and correctly predicted the phase changes of the upwind and crosswind asymmetries in T_υ and U from middle to high incidence angles. The phase changes predicted by the Bragg scattering theory for radiometric emission from water ripples is corroborated by the numerical Monte Carlo simulation of rough surface scattering. This theoretical interpretation indicates the potential use of polarimetric brightness temperatures for retrieving the directional wave spectrum of short gravity and capillary waves     

Evaluation of WindSat wind vector performance with respect to QuikSCAT estimates

The WindSat instrument was launched on January 6, 2003 as part of a risk reduction effort to assess the potential of using spaceborne fully polarimetric radiometry to measure the marine wind vector. Microwave radiometry on the Special Sensor Microwave/Imager onboard the Defense Meteorological Satellite Program satellites has long provided wind speed measurements. Fully polarimetric radiometry offers the additional possibility of obtaining wind direction as well. By contrast, the QuikSCAT satellite uses active microwave measurements to estimate the wind vector from space. It represents the most comprehensive satellite dataset against which to compare WindSat measurements. In this paper, we systematically compare temporally and spatially coincident WindSat and QuikSCAT wind vector measurements against the design goals of the WindSat instrument, taking into consideration expected differences related to instrument precision and the spatial and temporal variability of the wind field.     

Efficient Model-Based Estimation of Atmospheric Transmittance and Ocean Wind Vectors From WindSat Data

A new method for estimating the atmospheric transmittance and wind speed over the ocean from WindSat data is derived using a simplified model for the ocean surface reflectivity. The simplified reflectivity model is used to calculate both the surface emissivity and the reflection of downwelling atmospheric radiation. The wind-speed dependence of the surface reflectivity is parameterized using simple rational functions with coefficients determined from the WindSat data. Because the vertically polarized brightness temperature depends primarily on the atmospheric state, it is used to obtain an initial estimate of the atmospheric transmittance at each spatial location. These estimates are then combined with the horizontally polarized brightness temperature to estimate the wind speed at each location. The first wind-speed estimate is used to refine the estimate of the transmittance, and the process is repeated until the estimates converge, resulting in a simultaneous solution for the atmospheric transmittance and the wind speed. The results are illustrated for two WindSat data sets collected on September 12 and 14, 2003. We have also investigated two methods of estimating wind direction using WindSat measurements of the third and fourth Stokes parameters. The first method involves an algebraic solution for the wind direction from simultaneous measurements of the third and fourth Stokes parameters. The second method involves measurements of the third Stokes parameter from two look directions (fore and aft scan angles), made possible by the conical scanning geometry of WindSat. A comparison and evaluation of these methods is made using the same data sets.      

The accuracy of preliminary WindSat vector wind measurements: comparisons with NDBC buoys and QuikSCAT

Two preliminary, six-month long global WindSat vector wind datasets are validated using buoys and QuikSCAT measurements. Buoy comparisons yield speed and direction root mean square accuracies of 1.4 m/s and 25/spl deg/ for the "NESDIS0" product and 1.3 m/s and 23/spl deg/ for the more recently produced "B1" product from the Naval Research Laboratory. WindSat along- and across-wind random component errors of 0.7-1.0 and 2.6-2.8 m/s (respectively) are larger than those calculated for QuikSCAT in the same period. Global WindSat-QuikSCAT comparisons generally confirmed the buoy analyses. While simple rain flags based directly on WindSat brightness temperature measurements alone are shown to overflag for rain systematically, the advanced "Environmental Data Record" rain flag in the B1 product matches well with Special Sensor Microwave/Imager rain detection frequency and preserves the accuracy of the unflagged vector wind measurements.