Comparisons of observed ice properties and Envisat Advanced Synthetic Aperture Radar (ASAR) data during Nathanial B. Palmer cruise NBP09-01 and Oden cruise OSO08/09

Envisat Advanced Synthetic Aperture Radar (ASAR) Wide Swath Mode (WSM) images are used to derive C-band HH-polarization normalized radar cross sections (NRCS). These are compared with ice-core analysis and visual ship-based observations of snow and ice properties observed according to the Antarctic Sea Ice Processes and Climate (ASPeCt) protocol during two International Polar Year summer cruises (Oden 2008 and Palmer 2009) in West Antarctica. Thick first-year (TFY) and multi-year (MY) ice were the dominant ice types. The NRCS value ranges between -16.3 ± 1.1 and -7.6 ± 1.0 dB for TFY ice, and is -12.6 ± 1.3 dB for MY ice; for TFY ice, NRCS values increase from ~-15 dB to -9 dB from December/January to mid-February. In situ and ASPeCt observations are not, however, detailed enough to interpret the observed NRCS change over time. Co-located Advanced Microwave Scanning Radiometer-Earth Observing System (AMSR-E) vertically polarized 37 GHz brightness temperatures (TB37V), 7 day and 1 day averages as well as the TB37V difference between ascending and descending AMSR-E overpasses suggest the low NRCS values (-15 dB) are associated with snowmelt being still in progress, while the change towards higher NRCS values (-9dB) is caused by commencement of melt-refreeze cycles after about mid-January.

(Table 1) Ice thickness at lake Nam Co, Tibet, during 2008-2009

Lake ice change is one of the sensitive indicators of regional and global climate change. Different sources of data are used in monitoring lake ice phenology nowadays. Visible and Near Infrared bands of imagery (VNIR) are well suited for the observation of freshwater ice change, for example data from AVHRR and MODIS. Active and passive microwave data are also used for the observation of lake ice, e.g., from satellite altimetry and radiometry, backscattering coefficient from QuickSCAT, brightness temperature (Tb) from SSM/I, SMMR, and AMSR-E. Most of the studies are about lake ice cover phenology, while few studies focus on lake ice thickness. For example, Hall et al. using 5 GHz (6 cm) radiometer data showed a good relationship between Tb and ice thickness. Kang et al. found the seasonal evolution of Tb at 10.65 GHz and 18.7 GHz from AMSR-E to be strongly influenced by ice thickness. Many studies on lake ice phenology have been carried out since the 1970s in cold regions, especially in Canada, the USA, Europe, the Arctic, and Antarctica. However, on the Tibetan Plateau, very little research has focused on lake ice-cover change; only a small number of published papers on Qinghai Lake ice observations. The main goal of this study is to investigate the change in lake ice phenology at Nam Co on the Tibetan Plateau using MODIS and AMSR-E data (monitoring the date of freeze onset, the formation of stable ice cover, first appearance of water, and the complete disappearance of ice) during the period 2000-2009.

(Table 1) Mean temperatures in january and july, and total annual snowfall at Great Bear Lake and Great Slave Lake, Canada

The sensitivity of brightness temperature (T(B)) at 6.9, 10.7, and 18.7 GHz from Advanced Microwave Scanning Radiometer-Earth Observing System (AMSR-E) observations is investigated over five winter seasons (2002-2007) on Great Bear Lake and Great Slave Lake, Northwest Territories, Canada. The T(B) measurements are compared to ice thicknesses obtained with a previously validated thermodynamic lake ice model. Lake ice thickness is found to explain much of the increase of T(B) at 10.7 and 18.7 GHz. T(B) acquired at 18.7 GHz (V-pol) and 10.7 GHz (H-pol) shows the strongest relation with simulated lake ice thickness over the period of study (R**2 > 0.90). A comparison of the seasonal evolution of T(B) for a cold winter (2003-2004) and a warm winter (2005-2006) reveals that the relationship between T(B) and ice growth is stronger in the cold winter (2003-2004). Overall, this letter shows the high sensitivity of T(B) to ice growth and, thus, the potential of AMSR-E mid-frequency channels to estimate ice thickness on large northern lakes.

Near-real-time service provision during effusive crises at Etna and Stromboli: basis and implementation of satellite-based IR operations

Using the NEODAAS-Dundee AVHRR receiving station (Scotland), NEODAAS-Plymouth can provide calibrated brightness temperature data to end users or interim users in near-real time. Between 2000 and 2009 these data were used to undertake volcano hot spot detection, reporting and time-average discharge rate dissemination during effusive crises at Mount Etna and Stromboli (Italy). Data were passed via FTP, within an hour of image generation, to the hot spot detection system maintained at Hawaii Institute of Geophysics and Planetology (HIGP, University of Hawaii at Manoa, Honolulu, USA). Final product generation and quality control were completed manually at HIGP once a day, so as to provide information to onsite monitoring agencies for their incorporation into daily reporting duties to Italian Civil Protection. We here describe the processing and dissemination chain, which was designed so as to provide timely, useable, quality-controlled and relevant information for ‘one voice’ reporting by the responsible monitoring agencies.

Near-real-time service provision during effusive crises at Etna and Stromboli: basis and implementation of satellite-based IR operations

Using the NEODAAS-Dundee AVHRR receiving station (Scotland), NEODAAS-Plymouth can provide calibrated brightness temperature data to end users or interim users in near-real time. Between 2000 and 2009 these data were used to undertake volcano hot spot detection, reporting and time-average discharge rate dissemination during effusive crises at Mount Etna and Stromboli (Italy). Data were passed via FTP, within an hour of image generation, to the hot spot detection system maintained at Hawaii Institute of Geophysics and Planetology (HIGP, University of Hawaii at Manoa, Honolulu, USA). Final product generation and quality control were completed manually at HIGP once a day, so as to provide information to onsite monitoring agencies for their incorporation into daily reporting duties to Italian Civil Protection. We here describe the processing and dissemination chain, which was designed so as to provide timely, useable, quality-controlled and relevant information for ‘one voice’ reporting by the responsible monitoring agencies.

Burned Area Mapping in the North American Boreal Forest Using Terra-MODIS LTDR (2001–2011): A Comparison with the MCD45A1, MCD64A1 and BA GEOLAND-2 Products

An algorithm based on a Bayesian network classifier was adapted to produce 10-day burned area (BA) maps from the Long Term Data Record Version 3 (LTDR) at a spatial resolution of 0.05° (~5 km) for the North American boreal region from 2001 to 2011. The modified algorithm used the Brightness Temperature channel from the Moderate Resolution Imaging Spectroradiometer (MODIS) band 31 T31 (11.03 μm) instead of the Advanced Very High Resolution Radiometer (AVHRR) band T3 (3.75 μm). The accuracy of the BA-LTDR, the Collection 5.1 MODIS Burned Area (MCD45A1), the MODIS Collection 5.1 Direct Broadcast Monthly Burned Area (MCD64A1) and the Burned Area GEOLAND-2 (BA GEOLAND-2) products was assessed using reference data from the Alaska Fire Service (AFS) and the Canadian Forest Service National Fire Database (CFSNFD). The linear regression analysis of the burned area percentages of the MCD64A1 product using 40 km × 40 km grids versus the reference data for the years from 2001 to 2011 showed an agreement of R2 = 0.84 and a slope = 0.76, while the BA-LTDR showed an agreement of R2 = 0.75 and a slope = 0.69. These results represent an improvement over the MCD45A1 product, which showed an agreement of R2 = 0.67 and a slope = 0.42. The MCD64A1, BA-LTDR and MCD45A1 products underestimated the total burned area in the study region, whereas the BA GEOLAND-2 product overestimated it by approximately five-fold, with an agreement of R2 = 0.05. Despite MCD64A1 showing the best overall results, the BA-LTDR product proved to be an alternative for mapping burned areas in the North American boreal forest region compared with the other global BA products, even those with higher spatial/spectral resolution

New Possibilities for Geophysical Parameter Retrievals Opened by GCOM-W1 AMSR2

A new approach to retrieve sea surface wind speed (SWS) in tropical cyclones (TCs) from the Advanced Microwave Scanning Radiometer 2 (AMSR2) data is presented. Analysis of all six AMSR2 C- and X-band channel measurements over TCs is shown to efficiently help to separate the rain contribution. Corrected measurements at 6.9 and 10.65 GHz are then used to retrieve the SWS. Spatial and temporal collocation of AMSR2 and tropical rain measurement mission (TRMM) microwave instrument (TMI) data is then further used to empirically relate TMI rain rate (RR) product to RR estimates from AMSR2 in hurricanes. SWS estimates are validated with measurements from the stepped frequency microwave radiometer (SFMR). As further tested, more than 100 North Atlantic and North Pacific TCs are analyzed for the 2012–2014 period. Despite few particular cases, most SWS fields are in a very good agreement with TC center data on maximum wind speeds, radii of storm, and hurricane winds. As also compared, very high consistency between AMSR2 and L-band SMOS wind speed estimates are obtained, especially for the super typhoon Haiyan, to prove the high potential of AMSR2 measurements in TCs.

Radio-Frequency Interference Identification Over Oceans for C- and X-Band AMSR2 Channels

A new method for radio-frequency interference (RFI) contamination identification over open oceans for the two C-subbands and X-band of Advanced Microwave Scanning Radiometer 2 (AMSR2) channel measurements is suggested. The method is based both on the AMSR2 brightness temperature (T-B) modeling and on the analysis of AMSR2 measurements over oceans. The joint analysis of T-B spectral differences allowed to identify the relations between them and the limits of their variability, which are ensured by the changes in the environmental conditions. It was found that the constraints, based on the ratio of spectral differences, are more regionally and seasonally independent than the spectral differences themselves. Although not all possible RFI combinations are considered, the developed simple criteria can be used to detect most RFI-contaminated pixels over the World Ocean for AMSR2 measurements in two C-subbands and the X-band.

Temperature-dependent internal quantum efficiency of blue high-brightness light-emitting diodes

Internal quantum efficiency (IQE) of a blue high-brightness InGaN/GaN light-emitting diode (LED) was evaluated from the external quantum efficiency measured as a function of current at various temperatures ranged between 13 and 440 K. Processing the data with a novel evaluation procedure based on the ABC-model, we have determined the temperature-dependent IQE of the LED structure and light extraction efficiency of the LED chip. Separate evaluation of these parameters is helpful for further optimization of the heterostructure and chip designs. The data obtained enable making a guess on the temperature dependence of the radiative and Auger recombination coefficients, which may be important for identification of dominant mechanisms responsible for the efficiency droop in III-nitride LEDs. Thermal degradation of the LED performance in terms of the emission efficiency is also considered.

The role of bound chlorine in the brightness reversion of bleached hardwood kraft pulp

Our previous paper showed fragmentary evidence that pulp brightness reversion may be negatively affected by its organically bound chlorine (OX) content. A thorough investigation on eucalyptus kraft pulp led to the conclusion that OX increases reversion of certain pulps but this trend is not universal. Alkaline bleaching stages decrease reversion regardless of pulp OX content. Pulps bleached with high temperature chlorine dioxide revert less than those bleached with conventional chlorine dioxide in sequences ending with a chlorine dioxide stage but similarly in sequences ending with a final peroxide stage. The use of secondary condensate for pulp washing decreases reversion.

Extended optimal estimation techniques for sea surface temperature from the Spinning Enhanced Visible and Infra-Red Imager (SEVIRI)

Sea surface temperature (SST) can be estimated from day and night observations of the Spinning Enhanced Visible and Infra-Red Imager (SEVIRI) by optimal estimation (OE). We show that exploiting the 8.7 μm channel, in addition to the “traditional” wavelengths of 10.8 and 12.0 μm, improves OE SST retrieval statistics in validation. However, the main benefit is an improvement in the sensitivity of the SST estimate to variability in true SST. In a fair, single-pixel comparison, the 3-channel OE gives better results than the SST estimation technique presently operational within the Ocean and Sea Ice Satellite Application Facility. This operational technique is to use SST retrieval coefficients, followed by a bias-correction step informed by radiative transfer simulation. However, the operational technique has an additional “atmospheric correction smoothing”, which improves its noise performance, and hitherto had no analogue within the OE framework. Here, we propose an analogue to atmospheric correction smoothing, based on the expectation that atmospheric total column water vapour has a longer spatial correlation length scale than SST features. The approach extends the observations input to the OE to include the averaged brightness temperatures (BTs) of nearby clear-sky pixels, in addition to the BTs of the pixel for which SST is being retrieved. The retrieved quantities are then the single-pixel SST and the clear-sky total column water vapour averaged over the vicinity of the pixel. This reduces the noise in the retrieved SST significantly. The robust standard deviation of the new OE SST compared to matched drifting buoys becomes 0.39 K for all data. The smoothed OE gives SST sensitivity of 98% on average. This means that diurnal temperature variability and ocean frontal gradients are more faithfully estimated, and that the influence of the prior SST used is minimal (2%). This benefit is not available using traditional atmospheric correction smoothing.

Estimation of sea surface temperature from the Spinning Enhanced Visible and Infra Red Imager, improved using numerical weather prediction

Most of the operational Sea Surface Temperature (SST) products derived from satellite infrared radiometry use multi-spectral algorithms. They show, in general, reasonable performances with root mean square (RMS) residuals around 0.5 K when validated against buoy measurements, but have limitations, particularly a component of the retrieval error that relates to such algorithms' limited ability to cope with the full variability of atmospheric absorption and emission. We propose to use forecast atmospheric profiles and a radiative transfer model to simulate the algorithmic errors of multi-spectral algorithms. In the practical case of SST derived from the Spinning Enhanced Visible and Infrared Imager (SEVIRI) onboard Meteosat Second Generation (MSG), we demonstrate that simulated algorithmic errors do explain a significant component of the actual errors observed for the non linear (NL) split window algorithm in operational use at the Centre de Météorologie Spatiale (CMS). The simulated errors, used as correction terms, reduce significantly the regional biases of the NL algorithm as well as the standard deviation of the differences with drifting buoy measurements. The availability of atmospheric profiles associated with observed satellite-buoy differences allows us to analyze the origins of the main algorithmic errors observed in the SEVIRI field of view: a negative bias in the inter-tropical zone, and a mid-latitude positive bias. We demonstrate how these errors are explained by the sensitivity of observed brightness temperatures to the vertical distribution of water vapour, propagated through the SST retrieval algorithm.

Sea surface temperature from a geostationary satellite by optimal estimation

Optimal estimation (OE) is applied as a technique for retrieving sea surface temperature (SST) from thermal imagery obtained by the Spinning Enhanced Visible and Infra-Red Imager (SEVIRI) on Meteosat 9. OE requires simulation of observations as part of the retrieval process, and this is done here using numerical weather prediction fields and a fast radiative transfer model. Bias correction of the simulated brightness temperatures (BTs) is found to be a necessary step before retrieval, and is achieved by filtered averaging of simulations minus observations over a time period of 20 days and spatial scale of 2.5° in latitude and longitude. Throughout this study, BT observations are clear-sky averages over cells of size 0.5° in latitude and longitude. Results for the OE SST are compared to results using a traditional non-linear retrieval algorithm (“NLSST”), both validated against a set of 30108 night-time matches with drifting buoy observations. For the OE SST the mean difference with respect to drifter SSTs is − 0.01 K and the standard deviation is 0.47 K, compared to − 0.38 K and 0.70 K respectively for the NLSST algorithm. Perhaps more importantly, systematic biases in NLSST with respect to geographical location, atmospheric water vapour and satellite zenith angle are greatly reduced for the OE SST. However, the OE SST is calculated to have a lower sensitivity of retrieved SST to true SST variations than the NLSST. This feature would be a disadvantage for observing SST fronts and diurnal variability, and raises questions as to how best to exploit OE techniques at SEVIRI's full spatial resolution.

Optimal estimation of sea surface temperature from split-window observations

Optimal estimation (OE) improves sea surface temperature (SST) estimated from satellite infrared imagery in the “split-window”, in comparison to SST retrieved using the usual multi-channel (MCSST) or non-linear (NLSST) estimators. This is demonstrated using three months of observations of the Advanced Very High Resolution Radiometer (AVHRR) on the first Meteorological Operational satellite (Metop-A), matched in time and space to drifter SSTs collected on the global telecommunications system. There are 32,175 matches. The prior for the OE is forecast atmospheric fields from the Météo-France global numerical weather prediction system (ARPEGE), the forward model is RTTOV8.7, and a reduced state vector comprising SST and total column water vapour (TCWV) is used. Operational NLSST coefficients give mean and standard deviation (SD) of the difference between satellite and drifter SSTs of 0.00 and 0.72 K. The “best possible” NLSST and MCSST coefficients, empirically regressed on the data themselves, give zero mean difference and SDs of 0.66 K and 0.73 K respectively. Significant contributions to the global SD arise from regional systematic errors (biases) of several tenths of kelvin in the NLSST. With no bias corrections to either prior fields or forward model, the SSTs retrieved by OE minus drifter SSTs have mean and SD of − 0.16 and 0.49 K respectively. The reduction in SD below the “best possible” regression results shows that OE deals with structural limitations of the NLSST and MCSST algorithms. Using simple empirical bias corrections to improve the OE, retrieved minus drifter SSTs are obtained with mean and SD of − 0.06 and 0.44 K respectively. Regional biases are greatly reduced, such that the absolute bias is less than 0.1 K in 61% of 10°-latitude by 30°-longitude cells. OE also allows a statistic of the agreement between modelled and measured brightness temperatures to be calculated. We show that this measure is more efficient than the current system of confidence levels at identifying reliable retrievals, and that the best 75% of satellite SSTs by this measure have negligible bias and retrieval error of order 0.25 K.

Retrievals of sea surface temperature from infrared imagery: origin and form of systematic errors

We show that retrievals of sea surface temperature from satellite infrared imagery are prone to two forms of systematic error: prior error (familiar from the theory of atmospheric sounding) and error arising from nonlinearity. These errors have different complex geographical variations, related to the differing geographical distributions of the main geophysical variables that determine clear-sky brightness-temperatures over the oceans. We show that such errors arise as an intrinsic consequence of the form of the retrieval (rather than as a consequence of sub-optimally specified retrieval coefficients, as is often assumed) and that the pattern of observed errors can be simulated in detail using radiative-transfer modelling. The prior error has the linear form familiar from atmospheric sounding. A quadratic equation for nonlinearity error is derived, and it is verified that the nonlinearity error exhibits predominantly quadratic behaviour in this case.