Next Article in Journal
ARCTIS — A MATLAB® Toolbox for Archaeological Imaging Spectroscopy
Previous Article in Journal
Land Cover Characterization and Classification of Arctic Tundra Environments by Means of Polarized Synthetic Aperture X- and C-Band Radar (PolSAR) and Landsat 8 Multispectral Imagery — Richards Island, Canada
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Inter-Calibration of Satellite Passive Microwave Land Observations from AMSR-E and AMSR2 Using Overlapping FY3B-MWRI Sensor Measurements

1
Flathead Lake Biological Station, The University of Montana, Polson, MT 59860, USA
2
Numerical Terradynamic Simulation Group, The University of Montana, Missoula, MT 59812, USA
3
State Key Laboratory of Remote Sensing Science, Institute of Remote Sensing and Digital Earth, Chinese Academy of Sciences, Beijing 100101, China
4
National Satellite Meteorological Center, Beijing 100081, China
5
Earth System Science Interdisciplinary Center (ESSIC), University of Maryland, 5825 University Research Court, Suite 4001, College Park, MD 20740-3823, USA
*
Author to whom correspondence should be addressed.
Remote Sens. 2014, 6(9), 8594-8616; https://doi.org/10.3390/rs6098594
Submission received: 22 May 2014 / Revised: 1 September 2014 / Accepted: 10 September 2014 / Published: 16 September 2014

Abstract

: The development and continuity of consistent long-term data records from similar overlapping satellite observations is critical for global monitoring and environmental change assessments. We developed an empirical approach for inter-calibration of satellite microwave brightness temperature (Tb) records over land from the Advanced Microwave Scanning Radiometer for EOS (AMSR-E) and Microwave Scanning Radiometer 2 (AMSR2) using overlapping Tb observations from the Microwave Radiation Imager (MWRI). Double Differencing (DD) calculations revealed significant AMSR2 and MWRI biases relative to AMSR-E. Pixel-wise linear relationships were established from overlapping Tb records and used for calibrating MWRI and AMSR2 records to the AMSR-E baseline. The integrated multi-sensor Tb record was largely consistent over the major global vegetation and climate zones; sensor biases were generally well calibrated, though residual Tb differences inherent to different sensor configurations were still present. Daily surface air temperature estimates from the calibrated AMSR2 Tb inputs also showed favorable accuracy against independent measurements from 142 global weather stations (R2 ≥ 0.75, RMSE ≤ 3.64 °C), but with slightly lower accuracy than the AMSR-E baseline (R2 ≥ 0.78, RMSE ≤ 3.46 °C). The proposed method is promising for generating consistent, uninterrupted global land parameter records spanning the AMSR-E and continuing AMSR2 missions.

Graphical Abstract

1. Introduction

Microwave radiometers onboard Earth observing satellites have been widely used for observing critical land surface parameters supporting environmental monitoring and global change studies [14]. High-quality and consistent geophysical retrievals from well-calibrated, stable radiometer observations and long-term data records are essential for accurately delineating regional patterns and distinguishing natural climate variations from more subtle environmental trends.

The Advanced Microwave Scanning Radiometer for EOS (AMSR-E) was operational on board the NASA Aqua satellite from June 2002 to October 2011 and provided vertically and horizontally polarized, multi-frequency brightness temperature (Tb) observations of the global land surface at 1–3 day intervals [5]. A variety of global land products and geophysical retrievals have been derived from the AMSR-E Tb record, including atmosphere precipitable water vapor, surface air temperature, open water inundation, surface soil moisture, landscape freeze-thaw status, snow water equivalent, and vegetation canopy biomass optical depth [616]. These products have been refined through ongoing evaluations and reprocessings of the AMSR-E record to support global water, energy and carbon cycle studies [1721].

The NASA AMSR-E sensor ceased operations in October 2011 due to a failure of the rotational antenna spin mechanism, while continuing observations from the AMSR2 sensor onboard the JAXA GCOM-W1 (Global Change Observation Mission 1st-Water) satellite began in May 2012 and effectively replaces AMSR-E [22]. In principle, similar geophysical retrieval algorithms from AMSR-E can be applied to AMSR2 since the two sensors have similar configurations including frequency, overpass time, and spatial resolution; these cross-sensor synergies facilitate the development of consistent long-term data records needed by the science community. However, AMSR-E and AMSR2 have differing sensor calibration and design characteristics, and non-overlapping operational records that potentially constrain development of consistent long-term data records from these sensors. Careful inter-comparison and inter-calibration of the AMSR (henceforth referring to AMSR-E and AMSR2 collectively) observations are a pre-requisite for developing consistent land parameter retrievals encompassing both sensor records.

Besides AMSR-E and AMSR2, post-launch radiometer calibration and validation (Cal/Val) activities have been underway for many years involving similar sensors, including the Defense Meteorological Satellite Program (DMSP) Special Sensor Microwave Imager (SSM/I) [23,24], Tropical Rainfall Measuring Mission Microwave Imager (TRMM TMI) [25], and the recent Global Precipitation Measurement (GPM) mission [26]. As summarized in [27], the methodologies involved in sensor inter-calibration activities can be classified into three general groups: (1) the simultaneous nadir overpass (SNO) approach; (2) statistical inter-calibration; and (3) double differencing (DD) methods. The statistical inter-calibration approach compares sensor observations directly, and calibrates the target instrument data using empirical models developed from multiple (e.g., thousands) collocated observations; this approach is best suited for two instruments having collocated and near simultaneous observations. On the other hand, for instruments without acceptable collocation thresholds, the DD method, which utilizes an intermediate reference such as a third instrument or radiative transfer model (RTM), is often applied [27]. The AMSR sensors have non-overlapping operational records, without temporally collocated Tb observations; however, the DD method can be applied using an intermediate reference common to these records to quantify collocated Tb differences for consistent calibration of the AMSR-E and AMSR2 records. An advantage of the DD method is that the calibration is less likely to be affected by sensor configuration differences, including incident angle and central frequency of the Tb retrievals.

A critical issue affecting AMSR inter-calibration is a temporal gap of approximately nine months between the end of AMSR-E operations and subsequent initiation of the AMSR2 data record. Spatially collocated Tb observations from the two radiometers are inappropriate for direct comparisons due to temporally dynamic land and atmosphere parameters. Therefore, current inter-calibration efforts for AMSR-E and AMSR2 are based on the general DD method. Preliminary inter-calibration results for AMSR2 have been reported from the Intersatellite Radiometer Calibration Working Group (XCAL) and JAXA [28,29]. Their approach involves using RTM simulations as an intermediate reference for calibrating AMSR2 Tb retrievals, while selecting TMI or WindSat as the baseline sensor.

We adopted an alternative calibration approach aimed at producing consistent Tb and resulting land surface air temperature estimates extending from the AMSR-E to AMSR2 portions of record. The calibration process utilized the DD method and a third microwave radiometer record with collocated observations to both AMSR-E and AMSR2 to calibrate AMSR2 Tb observations to the earlier AMSR-E baseline record. The reference sensor used for calibration is the MWRI (Microwave Radiation Imager) on-board the Chinese FengYun 3B (FY3B) satellite, which was launched in November 2010 [30]. Similar overlapping MWRI observations to both AMSR-E and AMSR2 provided a unique opportunity for comparing and calibrating the sensor records. The MWRI has a similar instrument configuration and nearly simultaneous satellite overpass and data acquisition times as the AMSR sensors (i.e., 1:40 p.m./a.m. and 1:30 p.m./a.m. ascending/descending orbital nodes for MWRI and the AMSR sensors, respectively). The similar MWRI and AMSR overpass times mitigate potential Tb differences from diurnal changes relative to other potential reference sensors with different overpass times, including WindSat (∼6 AM/PM) and TMI (variable overpass time). The descriptions of the sensors and datasets used in this investigation are presented in Section 2, while the calibration process and validation results are described in Sections 3 and 4.

2. Instruments and Data

2.1. Instrument Similarities and Differences

AMSR-E on NASA’s EOS Aqua satellite was launched 4 May 2002 and operated until 4 October 2011. As the successor to AMSR-E, AMSR2 on-board the sun-synchronous satellite GCOM-W1 was launched on 18 May 2012 by JAXA. GCOM-W1 is the first of three planned satellites of the GCOM-W project, which is a 13-year mission designed for global and continuous observation of Earth’s water and energy cycles. AMSR2 has most of the characteristics of AMSR-E [31]. Both sensors are conical scanning passive microwave radiometers with similar instrument configurations. The AMSR instruments measure microwave emissions from the Earth surface twice daily for descending/ascending orbital equatorial crossings at 1:30 AM/PM local time, with vertically (V) and horizontally (H) polarized Tb retrievals at six frequencies (6.9, 10.7, 18.7, 23.8, 36.5, 89.0 GHz). Major changes in AMSR2, distinctive from AMSR-E, include an additional frequency at 7.3 GHz designed for mitigating Radio Frequency Interference (RFI), and a larger (2.0 m diameter) main reflector for enhanced spatial resolution. The MWRI onboard the FY-3B satellite was launched 5 November 2010. Similar to the channel setting and view geometry of AMSR-E, MWRI conically scans the Earth at five frequencies (10.7, 18.7, 23.8, 36.5, 89.0 GHz) and dual (H, V) polarizations. The MWRI has four major differences from the AMSR instruments, including: (a) no C-band channel relative to the AMSR 10.7 GHz channels due to widespread C-band RFI; (b) coarser spatial resolution and slightly narrower orbital swath for all frequencies relative to the AMSR sensors; (c) an approximate 10-min satellite overpass time difference between MWRI and AMSR-E; and (d) an approximate 53 degree earth incident angle (EIA) of measurement instead of 55 degrees for the AMSR sensors. Detailed sensor configurations for AMSR-E, AMSR2 and MWRI are summarized in Table 1.

2.2. Ancillary Datasets and Study Sites

Four ancillary datasets were employed for calibrating the AMSR2 Tb record to the AMSR-E baseline. The datasets used for calibration include: (a) AMSR-E L2A Version 3 (V003) re-sampled swath Tb data, which is the most recent version with complete recalibration to the Remote Sensing Systems (RSS) Version 7 standard; (b) AMSR2 L1R re-sampled swath Tb records (version 1.1); (c) the MWRI Level 1 swath Tb record; (d) daily minimum and maximum air temperature (Tmn and Tmx) measurements from global WMO weather stations that were used for independent validation of AMSR daily air temperature estimates derived from the calibrated Tb results.

All of the satellite ascending and descending orbital swath Tb data were spatially re-sampled and re-projected to a consistent 25-km resolution global Equal Area Scalable Earth (EASE) Grid (version 1) following previously established methods [32,33]. Based on the gridded data, four-month daily Tb records of spatially and temporally collocated observations for MWRI and AMSR-E (1 June 2011 to 30 September 2011), and MWRI and AMSR2 (1 June 2013 to 30 September 2013) were extracted over all global land areas. The two sets of matched datasets were then used to evaluate sensor Tb differences and develop pixel-wise empirical relationships for developing a consistent calibrated AMSR Tb record. For the sensor inter-calibration evaluation, collocated Tb observations of the three sensors involving both calibrated and uncalibrated results were extracted over five selected relatively homogeneous test regions, each approximately 1° × 1° in size, and representing the major Köppen-Geiger global climate classification zones (Polar, Continental, Tropical, Dry, and Temperate Climate types). Among the regions, the Tropical test area in Amazonia (Latitude: −6.35° to −7.14°, Longitude: −64.82° to −65.85°) was previously used for inter-sensor comparisons between AMSR-E and WindSat [34]. The Tb calibration was also assessed within major global land cover areas defined from a 0.05 degree resolution MODIS (MCD12C1) IGBP (International Geosphere-Biosphere Programme) land cover classification [35].

Surface (∼2 m height) air temperature estimates were also derived from both calibrated and uncalibrated Tb AM and PM orbital crossings over a three-year composited sensor record. The daily temperatures were derived for all global land areas using multi-frequency, ascending and descending orbit Tb inputs from the composited sensor record and temperature estimation algorithms previously developed for AMSR-E [7]. The resulting air temperatures determined from the AMSR-E baseline and calibrated AMSR2 Tb inputs were then compared with coincident in-situ daily air temperature measurements from 142 globally distributed WMO weather stations. Daily air temperatures determined from descending (∼1:30 a.m.) and ascending (∼1:30 p.m.) orbit AMSR Tb inputs were compared against collocated Tmn and Tmx station measurements based on the assumption that the satellite overpasses are approximately coincident with the respective diurnal temperature minima and maxima. These comparisons were used for evaluating potential differences in temperature estimation accuracy between the AMSR-E baseline and calibrated AMSR2 records, and as a relative quality assessment of the AMSR2 Tb calibration.

The composited three-year daily Tb record used for air temperature estimation was derived from one-year of AMSR-E data (1 September 2010–31 August 2011), one-year of calibrated AMSR2 data (1 September 2012–31 August 2013) and one-year of calibrated MWRI data (1 September 2011–31 August 2012). Both the MWRI and AMSR2 Tb records were calibrated to the AMSR-E baseline using the pixel-wise DD approach and collocated MWRI and AMSR Tb observations for all global land areas. The WMO air temperature datasets supporting the three-year comparisons were obtained from the National Climate Data Center (NCDC) Global Summary of the Day (GSOD version 7). Similar to previous work [7], the WMO validation stations were selected within major global land cover classes determined from the MODIS IGBP land cover classification [35], while the number of WMO stations selected within each land cover class was proportional to its global land area representation. The locations of the 142 global WMO temperature validation sites and the five Tb calibration test regions representing the major global climate zones are presented in Figure 1 and the corresponding inset map, respectively.

3. Methods

3.1. Double Differencing Method

As introduced in Section 1, Double Differencing (DD) is a straightforward and effective method for sensor inter-comparison and inter-calibration. By utilizing a third instrument as an intermediate reference, the DD method makes direct comparisons between target and reference sensors while avoiding potential uncertainty introduced by using a RTM, and mitigating negative impacts from systematic inconsistencies between the intermediate reference and inter-calibrated sensor records. The mathematical form of DD approach used in this study is expressed as:

DD ( AMSR 2 , AMSR-E ) = SD ( AMSR 2 , MWRI ) SD ( AMSR-E , MWRI ) SD ( AMSR 2 , MWRI ) = T b _ AMSR 2 T b _ MWRI SD ( AMSR-E , MWRI ) = T b _ AMSR E T b _ MWRI
where SD (Kelvin) represents the single difference of the collocated sensor Tb observations for a similar orbital retrieval, frequency and polarization. In this study, Equation (1) was applied for the two four-month matched sensor Tb datasets with calculations made on a per pixel basis for each channel except C-band, which was not part of the MWRI sensor design. According to the available literature, comparisons between the MWRI, AMSR-E and SSMIS sensors indicate relatively stable and accurate MWRI Tb measurements with calibration error of approximately ±1.5 K for all frequencies [36]. These favorable MWRI characteristics and similar collocated observations spanning the AMSR sensor records enable the observation bias of AMSR2 in comparison to AMSR-E to be estimated from Equation (1). However, despite general similarities between these sensors, there are four major differences between MWRI and AMSR sensors as described in Section 2. In particular, the MWRI has coarser sampling footprints for all channel frequencies and an approximate 10-min difference in sensor acquisition overpass time from the AMSR sensors. This implies that the MWRI observations may not be fully consistent with the AMSR observations, especially for complex terrain and heterogeneous landscapes. Therefore, correlations between the MWRI and AMSR Tb observations were first computed on a per grid cell basis over all global land areas. The subsequent estimation of DD values in Equation (1) were only determined for highly correlated pixels with Tb correlations (R) greater than 0.95 (p-value < 0.05).

3.2. Inter-Calibration of AMSR2 and AMSR-E

After describing observation biases between AMSR-E and AMSR2 using the above DD method, we then employed a linear calibration adjustment of the AMSR2 observations against the AMSR-E baseline. The linear relationships between collocated MWRI and AMSR Tb observations for the well-correlated pixels were described as:

T b _ AMSR E = a 1 + b 1 T b _ MWRI
T b _ AMSR 2 = a 2 + b 2 T b _ MWRI
where a1, b1, a2, b2 are empirical regression coefficients derived from the matched datasets for each pixel. Therefore,
T b _ AMSR E = a + b T b _ AMSR 2
where a = a1a2·b1/b2 and b = b1/b2.

For pixels where the Tb observations between AMSR and MWRI were well correlated (R > 0.95), coefficients, a and b, were calculated directly based on the pixel-wise linear regression results; otherwise an Inverse Distance Weighting (IDW) method was used for spatial interpolation estimation of a and b from nearby pixels having the same land cover type and successfully calculated parameter (a and b) values. The pixel-wise regression and IDW method are intended to account for possible spatial variations of sensor biases, which are illustrated in the following Subsection 4.1.

3.3. Air Temperature Estimation Method

Surface air temperature is an essential climate variable for determining the state of the climate system and monitoring environmental change. It also plays a critical role in regulating environmental processes and land–atmosphere interactions governing surface water, energy and carbon budgets [37,38]. In [7], Tmn and Tmx over the Northern Hemisphere were estimated based on an iterative method using multi-frequency Tb records from AMSR-E. After screening out potential precipitation, Radio Frequency Interference (RFI), snow and frozen conditions based on the procedure by [7], three Tb ratios from the AMSR-E 18.7 and 23.8 GHz channels are first employed iteratively to determine the optical thickness of atmospheric water vapor, land surface open water fraction and vegetation transmissivity. This information is then used to solve for the effective surface temperature, Ts, which was found to be well correlated with surface air temperature. By comparing with daily measurements from WMO weather stations, the resulting Ts retrievals for respective morning (AM) and afternoon (PM) overpasses were empirically related to Tmn and Tmx in the similar empirical regression form as Equations (5) and (6). The algorithm was developed for the Version 2 AMSR-E/Aqua L2A dataset and showed favorable accuracy (R2 ≥ 0.79, RMSE = 3.5 K) against daily temperature measurements from independent weather stations representing a broad range of climate and land cover conditions [7].

In this study, we applied the same algorithm as [7] but with regression coefficients adjusted to reflect reprocessing updates of the AMSR-E Tb record from Remote Sensing Systems (RSS) Version 6 to Version 7 as:

T mn = 2.05 + 0.89 T s + 1.83 T c 4.47 T c 2 0.10 ( abs ( Lat ) )
T mx = 6.21 + 0.86 T s + 19.71 T c + 23.15 T c 2 0.10 ( abs ( Lat ) )
where Tmn and Tmx (Celsius) are the respective daily minimum and maximum surface air temperature estimation of the satellite morning and afternoon overpasses; Tc is the dimensionless frequency dependent vegetation transmissivity; Ts is the effective surface temperature in Celsius, and Lat is the geographic latitude of the temperature estimation.

4. Results and Discussion

4.1. Inter-Calibration between AMSR2, AMSR-E and MWRI

Stable and well-correlated Tb observations from AMSR-E, AMSR2 and MWRI form the basis of the sensor inter-calibrations. The channel specific Tb correlations between the AMSR and MWRI sensors were thus examined first for the two, four-month overlapping data records. The linear correlation coefficient R was calculated from the overlapping sensor records for each 25 km resolution pixel within the global land domain, while the mean Tb sample size used to determine R for each pixel was 50.4 ± 29.6 days for the AMSR-E and MWRI overlapping data record and 49.2 ± 29.5 days for the AMSR2 and MWRI overlapping record. Generally favorable Tb correlations were observed over the global domain from the overlapping MWRI and AMSR records for all analyzed frequencies, orbits and polarizations, though the locations and global extent of significant correlations varied according to the different channel frequencies (e.g., Tables 2 and 3). As an illustrated example, correlation maps of collocated H-Polarized Tb ascending observations for the intermediate 23 GHz channel from the overlapping MWRI and AMSR records are presented in Figure 2. The global patterns of 23 GHz results are similar to the other frequency results, while the 23 GHz channel is also important for the subsequent air temperature estimation. More detailed summaries representing all frequencies are listed in Tables 2 and 3. The 23 GHz results from the MWRI and AMSR sensor records are highly correlated (R > 0.95) over more than half of the global domain. Regions with lower Tb correlations (R < 0.85) were located near large water bodies, including coastal areas, the Great Lakes and Amazon River. This is expected due to the variable spatial sampling and coverage of the collocated elliptical sensor footprints, which reflect satellite orbital differences and discrepancies in sensor spatial resolutions and geolocation accuracy. In particular, the inclusion of different fractional water coverage in the collocated footprints can lead to significant differences in the Tb observations due to the large dielectric contrast between water and land features. To minimize the effect of these uncertainties on the sensor Tb inter-comparison and calibration, only highly correlated pixels (R > 0.95) were selected for further evaluation of sensor observation biases. The overall Tb bias of MWRI against AMSR2 or AMSR-E was directly calculated by the SD method while the bias of AMSR2 against AMSR-E was indirectly estimated by the DD method using Equation (1). The resulting global distribution of estimated AMSR2 bias in relation to the baseline AMSR-E record is illustrated in Figure 3a,b for respective 23 GHz and 18 GHz, H-Polarized ascending Tb records. The AMSR2 ascending Tb record at 23 GHz and 18 GHz is approximately 2.10 K and 2.46 K higher than the AMSR-E record over the global domain. Higher Tb bias was generally found in regions of colder radiance, such as the Polar regions and Sahara Desert, where lower Tb observations are caused by low surface temperature and/or low H-pol emissivity. This finding is consistent with the previous study from JAXA that showed much larger AMSR2 biases toward the cold end of the Tb range [28,29].

The complete Tb inter-comparison results for the three sensor records and all matching frequencies are summarized in Table 2 for ascending orbits and Table 3 for descending orbits. The statistical metrics for the overlapping data records include the mean Tb difference or bias, the standard deviation (STD) of the Tb difference and the percentage of the global land area with R > 0.95. By referring to AMSR-E as the standard or baseline, it can be seen from Tables 2 and 3 that AMSR2 tends to generate higher Tb readings for all channels, while MWRI generally produces lower Tb values. Meanwhile, MWRI and AMSR observations are highly correlated (R > 0.95) over approximately 27%–69% of the global land area, depending on Tb frequency and polarization. As discussed in the previous paragraph, variable spatial sampling and coverage of the collocated elliptical sensor footprints may lead to low correlation of the observations. The cross-sensor Tb correlation is also lowest for the 89 GHz channel as indicated from Tables 2 and 3. As introduced in Subsection 2.1, there is an approximate 10-min overpass time difference between AMSR and MWRI and changing atmosphere conditions may lead to a lower correlation. The calibration correction of AMSR2 Tb to the AMSR-E baseline was conducted individually for each channel frequency, polarization and orbit. Based on Tables 2 and 3, global mean biases between AMSR2 and AMSR-E observations can be estimated for each sensor channel. As is seen in Figure 4, the AMSR2 observations are systematically deviated from the AMSR-E baseline by approximately 0.15 K–3.06 K, depending on the Tb frequency, polarization and ascending/descending orbits.

Based on the observed strong linear correlations between the sensor observations, empirical linear relationships were derived between the AMSR Tb observations using Equations (2)(4) to calibrate the AMSR2 record to the AMSR-E baseline. Based on the MWRI matched datasets, the intercept and slope of the relationship between AMSR2 and AMSR-E were directly calculated in Equations (2)(4) for well correlated land pixels (R > 0.95), and spatially interpolated over the remaining land areas using the IDW method. The calculation was carried out for each channel on a per pixel basis. Example Tb slope and intercept distribution maps describing relations between the AMSR2 and AMSR-E records is given in Figure 5 for the 23 GHz, H-Polarized ascending Tb observations. Assuming the relationships between the sensors are stable over the observation record, the empirical slope and intercept maps provide effective ancillary inputs for calibrating the AMSR2 and MWRI Tb records to the AMSR-E baseline. Using this approach and the overlapping sensor Tb observations, the available AMSR2 (July 2012–September 2013) and MWRI (June 2011–September 2013) records were linearly calibrated to the AMSR-E (June 2002–October 2011) Tb record. The resulting integrated Tb observations from the three sensors spanned a continuous, calibrated daily observation record extending from June 2002–September 2013, and consistent with the AMSR-E baseline.

4.2. Evaluation of the Sensor Inter-Calibration Results

The performance of the sensor Tb inter-calibration results was examined for five, homogeneous 1° × 1° Tb calibration test regions representing major Köppen-Geiger global climate classification zones, including Polar, Continental, Tropical, Dry, and Temperate climate types (Figure 1). Spatial mean Tb values were calculated for each test region. To ensure relatively homogenous conditions suitable for inter-calibration, the test regions were selected on the basis of having standard deviations of vertical polarization Tb values not exceeding 2 K, and standard deviations of horizontal polarization Tb values not exceeding 3 K [39]. The Tb calibration test results representing all frequencies and polarizations are summarized Tables 4 to 5. These results indicate that the Tb inter-calibration method significantly mitigates sensor observation biases for all channels over the broad range of global climate conditions represented by the five test regions.

The selected tropical Tb calibration test area in central Amazonia (Figure 1) was identified as relatively homogenous at microwave frequencies and suitable for sensor inter-comparisons from a previous study [34]. Both un-calibrated and calibrated AMSR2 and MWRI Tb (H-Polarization, and Ascending orbit) values for the 23 GHz, 18 GHz and 10 GHz channels are compared with collocated AMSR-E baseline observations for this region in Figure 6. For the 23 GHz example, there are 59 point pairs evaluated in both Figure 6a,b. These results indicate that the AMSR-E and AMSR2 Tb values (Y-axis) as a whole are well correlated with the MWRI data (R = 0.90) but also show large Tb differences of approximately −3.14 K (Bias) and 3.31 K (RMSE). The resulting empirical calibration of AMSR2 and MWRI observations to the AMSR-E baseline described in Section 3 produced generally consistent Tb observations across the three sensor records (Figure 6b). The resulting calibrated AMSR2 and MWRI Tb values for the reference area correspond favorably with the AMSR-E baseline (R = 0.99, RMSE = 0.42 K, Bias = 0.04 K), and are consistent with the inter-calibration results from the other climate test regions. Similarly, instrument biases were corrected for the other frequencies. As illustrated in Figure 6c–f, R increased from 0.82–0.98 for 18 GHz and 0.76–0.98 for 10 GHz, respectively; bias also decreased from −0.45 K to −0.10 K for 18 GHz and −3.60 K to 0.24 K for 10 GHz, respectively. Note that Tb observations from AMSR sensors are generally higher than the uncalibrated MWRI observations except for 18 GHz channels where slightly lower AMSR-E Tb values were found as shown in Figure 6c and also Tables 2 and 3.

The global Tb calibration results were also summarized for the major global land cover classes according to Tb ascending (Table 6) and descending (Table 7) orbits. The statistical parameters presented include the mean bias, RMSE and correlation coefficient which are averaged for all channels for the uncalibrated and calibrated results within each land cover class. Similar to climate based Tb calibration test, the land cover assessment results also indicate that the inter-calibration process mitigated most observation biases. Further analysis of residual errors from the resulting calibration showed no apparent non-linear bias or other artifacts in the integrated Tb record. However, the RMSE differences were only partially reduced and the correlations between datasets were slightly improved. The Tb differences from AMSR and MWRI are not only caused by different instrument calibrations, but also by other factors including different sensor spatial coverage, temporal sampling and random uncertainty in the observations. With observation biases mitigated, the calibrated AMSR2 and MWRI Tb values are similar to the AMSR-E baseline, though further integration refinement of the sensor records would likely benefit by more explicit accounting of the different sensor spatial and temporal characteristics.

4.3. Air Temperature Estimation Results

The global integrated sensor Tb record was used as a daily ancillary input for estimating surface air temperatures using algorithms Equations (5) and (6) adapted from AMSR-E [7]. Daily Tmx and Tmn estimates were derived from respective ascending (PM) and descending (AM) orbit Tb inputs and used for further assessment of the integrated sensor record. The resulting satellite temperature (Tmx and Tmn) estimates were validated against coincident in-situ air temperature measurements from 142 independent weather stations encompassing a broad range of global climate and land cover conditions. The air temperature estimates were extracted for three, one-year periods representing baseline AMSR-E and calibrated MWRI and AMSR2 portions of the integrated Tb record used for the validation assessment.

As introduced in Section 3.3, the air temperature estimates directly depend on the Tb observations at 18 and 23 GHz channels; the temperature estimates also indirectly depend on other channels used to identify and screen out unfavorable retrieval conditions, including high fractional water inundation, active precipitation or frozen surface conditions [7]. Therefore, the air temperature estimation is sensitive to all Tb channels ranging from 10–89 GHz. We assumed that the overall quality of the in-situ measurements from the 142 WMO validation sites was consistent over the three sensor years, so that the resulting temperature validation reflects the quality and consistency of both the Tb calibration and temperature estimation from the different sensors. The comparisons between measured and estimated Tmx from the AMSR-E baseline and calibrated AMSR2 Tb records are shown in Figure 7 for the 142 validation sites. The global temperature estimates from the AMSR-E baseline show generally good correlation (R2 = 0.86) and low error (RMSE = 3.37 °C) relative to the in situ station measurements (Figure 6a), and consistent with temperature validation results from a previous Northern Hemisphere study [7]. The temperature estimates from the calibrated AMSR2 Tb record were also favorable (R2 = 0.83 and RMSE = 3.64 °C) but represented a small, but significant (p < 0.05) accuracy reduction from the AMSR-E baseline. The overall statistical results of the temperature validation comparisons are summarized in Table 8. Overall, the AMSR temperature estimates are well correlated with the in-situ measurements and show favorable accuracy (R2 ≥ 0.75, RMSE ≤ 3.64 °C). These results also indicate that the calibrated AMSR2 data record can be used for extending the AMSR-E land parameter record with reasonable accuracy and consistency. However, the calibrated MWRI Tb record produced substantially lower temperature estimation accuracy relative to the station measurements (Table 8), with both lower temperature correspondence (R2 = 0.80 for Tmx and R2 = 0.69 for Tmn) and larger RMSE (4.12 °C for Tmx and 4.11 °C for Tmn) differences than the AMSR results. The lower MWRI temperature accuracy was primarily attributed to the coarser spatial resolution of the MWRI footprint relative to the AMSR sensors, and other factors described in the following section.

4.4. Discussion

The relative Tb bias from AMSR2 was detected, calibrated and mitigated against the AMSR-E baseline. However, slightly lower temperature estimation accuracy from the calibrated AMSR2 Tb record was also found. The lower AMSR2 temperature estimation accuracy may reflect a limited Tb calibration period represented by the available four-month matched datasets, which do not fully represent potential seasonal or annual changes in the overlapping sensor relationships. As indicated in Subsections 4.2 and 4.3, a number of factors may cause the observation differences between sensors, including different instrument calibrations, differences in spatial coverage (e.g., spatial resolution, geo-location of the footprint center, and footprint sampling), temporal sampling and random uncertainty in the observations. The DD method is capable of minimizing observation biases and partially mitigating the effects of the different sensor configurations on the Tb retrievals, though the residual impacts of these differences, as well as observation uncertainties inherent to each instrument cannot be completely eliminated using the relatively simple DD inter-calibration method.

The impacts of the inherent differences in sensor observations are also reflected from the lower MWRI performance. Coarser MWRI spatial footprints relative to the other sensors may be the primary factor contributing to lower sensor performance. For example, the 18 GHz channel footprints from AMSR-E and AMSR2 have respective mean spatial resolutions of 21 km and 18 km, whereas the corresponding MWRI resolution is 39 km, or roughly four times the AMSR footprints. Although the orbital swath observations from the sensors were gridded to a consistent 25 km resolution global EASE-grid, the resulting Tb records after bias correction by the DD method still preserve underlying differences in footprint resolution and sampling from the individual channel frequencies. Moreover, the resampling process itself inevitably causes data losses and its related uncertainties can also affect inter-calibration accuracy.

Future data records with longer overlapping periods and more sophisticated calibrations that consider the inherent differences in sensor configurations, especially different sensor spatial resolutions, would likely improve the accuracy of cross-sensor data integration efforts. Moreover, both MWRI and AMSR2 are relatively new sensors, while AMSR-E encompasses a much longer data record and has undergone multiple data reprocessing events, each with relative improvements over earlier versions. With continuing observations and similar data reprocessing and sensor calibration enhancements, both AMSR2 and MWRI are expected to show improved performance leading to more accurate integrated long-term environmental data records.

5. Conclusions

Consistent Tb records from similar overlapping satellite microwave sensors provide the basis for generating long-term global land parameter observation records for documenting and diagnosing environmental change. We developed an empirical approach for inter-calibrating AMSR-E and AMSR2 Tb records over land using similar overlapping satellite passive microwave observations from the FY3B-MWRI sensor. Our global pixel-wise correlation analysis indicates that the MWRI and AMSR sensor observations are highly correlated (R > 0.95) over approximately 27%–69% of the global land area, depending on sensor frequency and polarization. The Tb correlation analysis was developed from limited, four-month overlapping MWRI observations with the AMSR2 and AMSR-E sensor records, but indicated significant MWRI and AMSR2 Tb biases in relation to the AMSR-E baseline observations. The sensor Tb biases vary with channel frequency and polarization, and satellite orbital nodes, though AMSR2 generally produced the highest Tb readings while MWRI had the lowest. As a prerequisite for integrating the individual sensor observations into a continuous and consistent data record, empirical linear relationships were established between the AMSR2 and AMSR-E Tb retrievals for pixels where both AMSR sensor records were highly correlated (R > 0.95) with the overlapping MWRI observations. Pixels not meeting this high correlation requirement generally represented coastal boundaries and other heterogeneous land areas, where the sensor Tb relationships were estimated by spatially interpolating regression coefficients from nearby pixels having the same land cover type, but strong sensor correlations. Similar relationships were established between the MWRI and AMSR-E baseline Tb observations. The AMSR2 and MWRI Tb records were then calibrated to the AMSR-E baseline using the DD approach. The resulting calibrated 23 GHz H-Pol Tb observations from the different sensors were found to be generally consistent with each other for a homogeneous Amazon reference area (R = 0.99, RMSE = 0.42 K, Bias = 0.04 K) and other Tb calibration test regions representing major global climate zones. For the larger global domain, the calibrated AMSR2 and MWRI Tb values are comparable with the AMSR-E baseline, with mean biases for the major land cover areas ranging from −0.81 K to 0.09 K. We also found no evidence of non-linear variations in the sensor Tb relationships, though our analysis was limited by the relatively short, four month overlapping MWRI and AMSR records. The calibrated sensor records were used to obtain daily surface air temperature (Tmx, Tmn) estimates using multi-frequency Tb inputs and a modified iterative RT retrieval algorithm [7]. The resulting temperature estimates were validated against independent temperature measurements from 142 global weather stations, as a means for further assessing the relative quality and performance of the integrated satellite Tb record. The temperature results indicate successful global calibration and integration of the AMSR Tb records, and favorable temperature estimation accuracy (R2 ≥ 0.75, RMSE ≤ 3.64 °C), but with lower temperature accuracy from the calibrated MWRI record (R2 ≥ 0.69, RMSE ≤ 4.12 °C). Despite the lower MWRI performance, similar, albeit coarser resolution, overlapping Tb observations from this operational sensor generally provide an effective means for harmonizing ongoing AMSR2 observations with the earlier AMSR-E baseline record. The overall Tb calibration and temperature estimation accuracy from the integrated satellite record is expected to improve through continuing observations and reprocessing improvements of the MWRI and AMSR2 records.

Multi-frequency passive microwave sensors, including AMSR and MWRI, are less affected by solar illumination and atmosphere constraints relative to optical-infrared sensors and are well-suited for near-daily global monitoring. As indicated in this study, the combined observations from these instruments enable consistent and continuing global environmental monitoring spanning multiple years and satellite data records. Spatially and temporally complementary, multi-frequency Tb observations from these instruments also provide a large degree of sensitivity and potential information content for isolating different contributions to the total microwave emission signal, including atmosphere, soil moisture, vegetation and open water effects, enabling land parameter retrievals with favorable accuracy. Finally, our results indicate that synergistic observations from the MWRI and AMSR sensors provide an effective means for developing continuing, long-term global land parameter records with suitable performance and accuracy for global change assessment and monitoring.

Acknowledgments

The AMSR2 L1R brightness temperature data used for this study were provided courtesy of JAXA. AMSR-E data are produced by Remote Sensing Systems and sponsored by the NASA Earth Science MEaSUREs DISCOVER Project and the AMSR-E Science Team. AMSR-E data were also provided courtesy of the National Snow and Ice Data Center (NSIDC). The WMO weather station data were provided courtesy of the National Climate Data Center (NCDC). This work was conducted at the University of Montana with funding provided by the NASA Terra/Aqua Science and Making Earth System data records for Use in Research Environments (MEaSUREs) programs (NNX11AD46G, NNX14AB20A).

Author Contributions

Jinyang Du performed inter-sensor calibrations and air temperature estimations using the integrated satellite data record; John S. Kimball developed the general scheme of the sensor inter-calibration work and contributed substantially to the writing of the manuscript; Jiancheng Shi performed evaluations on the inter-sensor calibration results; Lucas A. Jones developed AMSR-E air temperature estimation algorithms and evaluated the inter-calibration results; Shengli Wu, Ruijing Sun and Hu Yang provided initial processes on MWRI datasets.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Jones, M.O.; Kimball, J.S.; McDonald, K.C.; Jones, L.A. Utilizing satellite passive microwave remote sensing for monitoring global land surface phenology. Remote Sens. Environ 2011, 115, 1102–1114. [Google Scholar]
  2. Kim, Y.; Kimball, J.S.; Zhang, K.; McDonald, K.C. Satellite detection of increasing northern hemisphere non-frozen seasons from 1979 to 2008: Implications for regional vegetation growth. Remote Sens. Environ 2012, 121, 472–487. [Google Scholar]
  3. Watts, J.; Kimball, J.S.; Jones, L.A.; Schroeder, R.; McDonald, K.C. Satellite microwave remote sensing of contrasting surface water inundation changes within the Arctic-Boreal region. Remote Sens. Environ 2012, 127, 223–236. [Google Scholar]
  4. Chuang, T.W.; Henebry, G.M.; Kimball, J.S.; VanRoekel-Patton, D.L.; Hildreth, M.B.; Wimberly, M.C. Satellite microwave remote sensing for mosquito population dynamics. Remote Sens. Environ 2012, 125, 147–156. [Google Scholar]
  5. Kawanishi, T.J.; Sezai, T.; Ito, Y.; Imaoka, K.; Takeshima, T.; Ishido, Y.; Shibata, A.; Miura, M.; Inahata, H.; Spencer, R.W.; et al. The Advanced Scanning Microwave Radiometer for the Earth Observing System (AMSR-E): NASDA’s contribution to the EOS for global energy and water cycle studies. IEEE Trans. Geosci. Remote Sens 2003, 41, 184–194. [Google Scholar]
  6. Daily Global Land Surface Parameters Derived from AMSR-E, Boulder Colorado USA: National Snow and Ice Data Center. 2010. Available online: http://nsidc.org/data/nsidc-0451.html (accessed on 15 September 2014).
  7. Jones, L.A.; Ferguson, C.R.; Kimball, J.S.; Zhang, K.; Chan, S.T.K.; McDonald, K.C.; Njoku, E.; Wood, E. Satellite microwave remote sensing of daily land surface air temperature minima and maxima from AMSR-E. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens 2010, 3, 111–123. [Google Scholar]
  8. MEaSUREs Global Record of Daily Landscape Freeze/Thaw Status, Version 02 [1979 to 2010]; National Snow and Ice Data Center: Boulder, CO, USA, 2012. Available online: http://nsidc.org/data/nsidc-0477.html (accessed on 15 September 2014).
  9. Jackson, T.J.; Cosh, M.H.; Bindlish, R.; Starks, P.J.; Bosch, D.; Seyfried, M.; Goodrich, D.C.; Moran, M.S.; Du, J. Validation of advanced microwave scanning radiometer soil moisture products. IEEE Trans. Geosci. Remote Sens 2010, 48, 4256–4272. [Google Scholar]
  10. Koike, T.; Nakamura, Y.; Kaihotsu, I.; Davva, G.; Matsuura, N.; Tamagawa, K.; Fujii, H. Development of an Advanced Microwave Scanning Radiometer (AMSR-E) algorithm of soil moisture and vegetation water content. Proc. Hydraul. Eng 2004, 48, 217–222. [Google Scholar]
  11. Mladenova, I.E.; Jackson, T.J.; Njoku, E.; Bindlish, R.; Chan, S.; Cosh, M.H.; Holmes, T.R.H.; de Jeu, R.A.M.; Jones, L.; Kimball, J.; et al. Remote monitoring of soil moisture using passive microwave based techniques—Theoretical basis and overview of selected algorithms for AMSR-E. Remote Sens. Environ 2014, 144, 197–213. [Google Scholar]
  12. Njoku, E.G.; Chan, S.K. Vegetation and surface roughness effects on AMSR-E land observations. Remote Sens. Environ 2006, 100, 190–199. [Google Scholar]
  13. Paloscia, S.; Macelloni, G.; Santi, E. Soil moisture estimates from AMSR-E brightness temperatures by using a dual-frequency algorithm. IEEE Trans. Geosci. Remote Sens 2006, 44, 3135–3144. [Google Scholar]
  14. Kelly, R.E.; Chang, A.T.; Tsang, L.; Foster, J.L. A prototype AMSR-E global snow area and snow depth algorithm. IEEE Trans. Geosci. Remote Sens 2003, 41, 230–242. [Google Scholar]
  15. Santi, E.; Pettinato, S.; Paloscia, S.; Pampaloni, P.; Fontanelli, G.; Crepaz, A.; Valt, M. Monitoring of Alpine snow using satellite radiometers and artificial neural networks. Remote Sens. Environ 2014, 144, 179–186. [Google Scholar]
  16. Tedesco, M.; Narvekar, P.S. Assessment of the NASA AMSR-E SWE product. IEEE J. Sel. Top. Appl. Earth Obs. Remote Sens 2010, 3, 141–159. [Google Scholar]
  17. Kimball, J.S.; Jones, L.A.; Zhang, K.; Heinsch, F.A.; McDonald, K.C.; Oechel, W.C. A satellite approach to estimate land-atmosphere CO2 exchange for Boreal and Arctic biomes using MODIS and AMSR-E. IEEE Trans. Geosci. Remote Sens 2009, 47, 569–587. [Google Scholar]
  18. Meissner, T.; Wentz, F. Intercalibration of AMSR-E and Windsat brightness temperature measurements over land scenes. Proceedings of the IEEE International Geoscience and Remote Sensing Symposium, Honolulu, HI, USA, 25–30 July 2010; pp. 3218–3219.
  19. Mu, Q.; Jones, L.A.; Kimball, J.S.; McDonald, K.C.; Running, S.W. Satellite assessment of land surface evapotranspiration for the pan-Arctic domain. Water Resour. Res 2009, 45. [Google Scholar] [CrossRef]
  20. Jones, M.O.; Kimball, J.S.; Jones, L.A.; McDonald, K.C. Satellite passive microwave detection of North America start of season phenology. Remote Sens. Environ 2012, 123, 324–333. [Google Scholar]
  21. Tuttle, S.E.; Salvucci, G.D. A new approach for validating satellite estimates of soil moisture using large-scale precipitation: Comparing AMSR-E products. Remote Sens. Environ 2014, 142, 207–222. [Google Scholar]
  22. Imaoka, K.; Takashi, M.; Misako, K.; Marehito, K.; Norimasa, I.; Keizo, N. Status of AMSR2 instrument on GCOM-W1, earth observing missions and sensors: Development, implementation, and characterization II. Proc. SPIE 2012, 852815. [Google Scholar] [CrossRef]
  23. Colton, M.; Poe, G. Intersensor calibration of DMSP SSM/I’s: F-8 to F-14, 1987–1997. IEEE Trans. Geosci. Remote Sens 1999, 37, 418–439. [Google Scholar]
  24. Yang, S.; Weng, F.; Yan, B.; Sun, N.; Goldberg, M. Special Sensor Microwave Imager (SSM/I) intersensor calibration using a simultaneous conical overpass technique. J. Appl. Meteorol. Climatol 2011, 50, 77–95. [Google Scholar]
  25. Wentz, F.; Ashcroft, P.; Gentemann, C. Post-launch calibration of the TRMM microwave imager. IEEE Trans. Geosci. Remote Sens 2001, 39, 415–422. [Google Scholar]
  26. Biswas, S.; Farrar, S.; Gopalan, K.; Santos-Gracia, A.; Jones, W.L.; Billanow, S. Intercalibration of microwave radiometer brightness temperatures for the global precipitation measurement mission. IEEE Trans. Geosci. Remote Sens 2013, 51, 1465–1477. [Google Scholar]
  27. Chander, G.; Hewison, T.J.; Fox, N.; Wu, X.; Xiong, X.; Blackwell, W. Overview of inter-calibration of satellite instruments. IEEE Trans. Geosci. Remote Sens 2013, 51, 1056–1080. [Google Scholar]
  28. GPM X-Cal Results for AMSR2. Available online: http://smsc.cnes.fr/MEGHAT/XCAL-Workshop/XCal_Toulouse_AMSR2.pdf (accessed on 1 April 2014).
  29. Results of Intercalibration between AMSR2 and TMI/AMSR-E. Available online: http://suzaku.eorc.jaxa.jp/GCOM_W/materials/product/AMSR2_XcalResults_1-1.pdf (accessed on 1 April 2014).
  30. Yang, H.; Weng, F.; Lv, L.; Lu, N.; Liu, G.; Bai, M.; Qian, Q.; He, J.; Xu, H. The FengYun-3 microwave radiation imager on-orbit verification. IEEE Trans. Geosci. Remote Sens 2011, 49, 4552–4560. [Google Scholar]
  31. Imaoka, K.; Kachi, M.; Kasahara, M.; Ito, N.; Nakagawa, K.; Oki, T. Instrument performance and calibration Of AMSR-E and AMSR2. Proceedings of ISPRS Commission VIII Mid-Term Symposium “Networking the World with Remote Sensing”, Kyoto, Japan, 9–12 August 2010; XXXVIII. Part 8..
  32. Algorithm Theoretical Basis Document, AMSR 2A Algorithm; RSS Technical Report 121 599B-1; Santa Rosa, CA, USA, 1999. Available online: https://nsidc.org/data/amsre/pdfs/amsr_atbd_level2a.pdf (accessed on 1 April 2014).
  33. Armstrong, R.L.; Brodzik, M.J. An earth-gridded SSM/I data set for cryospheric studies and global change monitoring. Adv. Space Res 1995, 16, 155–163. [Google Scholar]
  34. Das, N.N.; Colliander, A.; Chan, S.K.; Njoku, E.G.; Li, L. Intercomparisons of brightness temperature observations over land from AMSR-E and WindSat. IEEE Trans. Geosci. Remote Sens 2014, 52, 452–464. [Google Scholar]
  35. Justice, C.O.; Townshend, J.R.G.; Vermote, E.F.; Masuoka, E.; Wolfe, R.E.; Saleous, N.; Roy, D.P.; Morisette, J.T. An overview of MODIS land data processing and product status. Remote Sens. Environ 2002, 83, 3–15. [Google Scholar]
  36. Yang, H.; Zou, X.; Li, X.; You, R. Environmental data records from FengYun-3B microwave radiation imager. IEEE Trans. Geosci. Remote Sens 2012, 50, 4986–4993. [Google Scholar]
  37. Bristow, K.L.; Campbell, G.S. On the relationship between incoming solar radiation and daily maximum and minimum temperature. Agric. For. Meteorol 1984, 31, 159–166. [Google Scholar]
  38. Dickinson, R.E. Land surface processes and climate modeling. Bull. Am. Meteorol. Soc 1996, 76, 1445–1448. [Google Scholar]
  39. Liang, H. Inter-Satellite Microwave Radiometer Calibration. Ph.D. Dissertation. University of Central Florida, Orlando, FL, USA, 2008. [Google Scholar]
Figure 1. Global distribution of the WMO weather stations (white circles) used for validating the satellite surface air temperature estimates; the validation sites represent the global range of climate and land cover conditions represented by the MODIS IGBP global land cover map; The locations of five homogeneous test regions representing the major climate zones are represented by black squares and labeled as P, C, TR, D and T for Polar, Continental, Tropical, Dry, and Temperate Climate types, respectively. These test regions were used for partial validation of the sensor Tb inter-calibration.
Figure 1. Global distribution of the WMO weather stations (white circles) used for validating the satellite surface air temperature estimates; the validation sites represent the global range of climate and land cover conditions represented by the MODIS IGBP global land cover map; The locations of five homogeneous test regions representing the major climate zones are represented by black squares and labeled as P, C, TR, D and T for Polar, Continental, Tropical, Dry, and Temperate Climate types, respectively. These test regions were used for partial validation of the sensor Tb inter-calibration.
Remotesensing 06 08594f1 1024
Figure 2. Spatial distribution of the correlation coefficients calculated from four-month overlapping (a) AMSR-E and MWRI and (b) AMSR2 and MWRI datasets for selected 23 GHz, H-Polarized ascending orbit Tb observations; lower correlation areas are associated with coastal boundaries having large dielectric differences between land and open water features, and other heterogeneous land areas.
Figure 2. Spatial distribution of the correlation coefficients calculated from four-month overlapping (a) AMSR-E and MWRI and (b) AMSR2 and MWRI datasets for selected 23 GHz, H-Polarized ascending orbit Tb observations; lower correlation areas are associated with coastal boundaries having large dielectric differences between land and open water features, and other heterogeneous land areas.
Remotesensing 06 08594f2 1024
Figure 3. Example global distributions of estimated ascending orbit Tb biases between uncalibrated AMSR2 and AMSR-E baseline observations for the H-Polarized (a) 23 GHz, and (b) 18 GHz channels, respectively (areas with correlation coefficient R < 0.95 are marked in grey).
Figure 3. Example global distributions of estimated ascending orbit Tb biases between uncalibrated AMSR2 and AMSR-E baseline observations for the H-Polarized (a) 23 GHz, and (b) 18 GHz channels, respectively (areas with correlation coefficient R < 0.95 are marked in grey).
Remotesensing 06 08594f3 1024
Figure 4. Global mean biases (Unit: Kelvin) between AMSR2 and AMSR-E estimated from overlapping AMSR and MWRI records for well correlated (R > 0.95) pixels for each frequency and polarization (blue line is for the ascending orbits and red line for descending orbits).
Figure 4. Global mean biases (Unit: Kelvin) between AMSR2 and AMSR-E estimated from overlapping AMSR and MWRI records for well correlated (R > 0.95) pixels for each frequency and polarization (blue line is for the ascending orbits and red line for descending orbits).
Remotesensing 06 08594f4 1024
Figure 5. Example global distributions of the empirical (a) intercept and (b) slope values that define the linear relationship between AMSR2 Tb observations (dependent variable) and the AMSR-E baseline (independent variable) for ascending H-Polarized observations at 23 GHz.
Figure 5. Example global distributions of the empirical (a) intercept and (b) slope values that define the linear relationship between AMSR2 Tb observations (dependent variable) and the AMSR-E baseline (independent variable) for ascending H-Polarized observations at 23 GHz.
Remotesensing 06 08594f5a 1024Remotesensing 06 08594f5b 1024
Figure 6. H-Polarized Tb comparisons from overlapping MWRI and AMSR ascending orbit observations for the selected Amazon tropical reference area; the comparisons represent both (left panel) original uncalibrated and (right panel) calibrated (a,b) 23 GHz (c,d) 18 GHz (e,f) 10 GHz MWRI and AMSR2 Tb results relative to the AMSR-E baseline for 23 GHz, 18 GHz and 10 GHz channels.
Figure 6. H-Polarized Tb comparisons from overlapping MWRI and AMSR ascending orbit observations for the selected Amazon tropical reference area; the comparisons represent both (left panel) original uncalibrated and (right panel) calibrated (a,b) 23 GHz (c,d) 18 GHz (e,f) 10 GHz MWRI and AMSR2 Tb results relative to the AMSR-E baseline for 23 GHz, 18 GHz and 10 GHz channels.
Remotesensing 06 08594f6 1024
Figure 7. Example correspondence between in-situ daily Tmx measurements from 142 WMO validation sites and corresponding temperature estimates from the (a) AMSR-E baseline and (b) calibrated AMSR2 observations.
Figure 7. Example correspondence between in-situ daily Tmx measurements from 142 WMO validation sites and corresponding temperature estimates from the (a) AMSR-E baseline and (b) calibrated AMSR2 observations.
Remotesensing 06 08594f7 1024
Table 1. Satellite microwave sensor configurations used in this study.
Table 1. Satellite microwave sensor configurations used in this study.
Instrument Configurations

SpecificationsAMSR-2AMSR-EMWRI

Satellite PlatformGCOM-W1AQUAFY3B
Altitude700 km705 km836 km
Equator Crossing Time (Local time zone)1:30 p.m. Ascending1:30 p.m. Ascending1:40 p.m. Ascending
1:30 a.m. Descending1:30 a.m. Descending1:40 a.m. Descending
Antenna Size2 m (Diameter)1.6 m (Diameter)0.977 m × 0.897 m
Incident Angle555553

Spatial Resolution (km × km)

Band (GHz)AMSR-2AMSR-EMWRI

6.9362 × 3575 × 43N/A
7.362 × 35N/AN/A
10.6542 × 2451 × 2985 × 51
18.722 × 1427 × 1650 × 30
23.819 × 1132 × 1845 × 27
36.512 × 714 × 830 × 18
89.05 × 36 × 415 × 9
Table 2. Summary of sensor ascending orbit Tb mean biases (Unit: Kelvin) and Standard Deviation (STD) (Unit: Kelvin) values estimated from overlapping AMSR and MWRI records for well correlated (R > 0.95) pixels, and the percentage of the global land domain with R > 0.95.
Table 2. Summary of sensor ascending orbit Tb mean biases (Unit: Kelvin) and Standard Deviation (STD) (Unit: Kelvin) values estimated from overlapping AMSR and MWRI records for well correlated (R > 0.95) pixels, and the percentage of the global land domain with R > 0.95.
Pairwise ComparisonsH_PolV_Pol

BiasSTDPercentageBiasSTDPercentage
10 GHz

MWRI-AMSR-E−2.782.0942.86−3.131.2959.46
MWRI-AMSR2−5.852.2240.51−4.611.5956.43

18 GHz

MWRI-AMSR-E−0.341.6847.220.201.1261.18
MWRI-AMSR2−2.801.5449.50−1.351.3461.79

23 GHz

MWRI-AMSR-E−2.581.3460.83−2.360.9866.01
MWRI-AMSR2−4.671.0365.21−3.750.9469.28

36 GHz

MWRI-AMSR-E−3.061.6949.07−4.011.1860.06
MWRI-AMSR2−4.391.2954.96−4.830.9863.20

89 GHz

MWRI-AMSR-E−1.831.6627.24−1.661.0632.95
MWRI-AMSR2−1.981.5428.94−2.331.0336.15
Table 3. Summary of sensor Descending orbit Tb mean biases (Unit: Kelvin) and Standard Deviation (STD) (Unit: Kelvin) values estimated from overlapping AMSR and MWRI records for well correlated (R > 0.95) pixels, and the percentage of the global land domain with R > 0.95.
Table 3. Summary of sensor Descending orbit Tb mean biases (Unit: Kelvin) and Standard Deviation (STD) (Unit: Kelvin) values estimated from overlapping AMSR and MWRI records for well correlated (R > 0.95) pixels, and the percentage of the global land domain with R > 0.95.
Pairwise ComparisonsH_PolV_Pol

BiasSTDPercentageBiasSTDPercentage
10 GHz

MWRI-AMSR-E−1.392.5336.10−2.231.5252.46
MWRI-AMSR2−4.042.7333.95−4.341.6748.92

18 GHz

MWRI-AMSR-E0.842.2541.761.071.4354.68
MWRI-AMSR2−0.932.2141.12−1.071.5253.26

23 GHz

MWRI-AMSR-E−1.162.2960.21−1.521.5661.53
MWRI-AMSR2−3.001.9760.30−3.191.4060.62

36 GHz

MWRI-AMSR-E−1.412.2847.22−2.841.3756.53
MWRI-AMSR2−3.001.8150.56−3.941.3057.71

89 GHz

MWRI-AMSR-E−0.201.6234.20−0.391.1434.42
MWRI-AMSR2−0.841.8534.57−1.271.4436.96
Table 4. Summary of sensor ascending (PM) orbit Tb mean biases (Unit: Kelvin), Root Mean Square Error (RMSE) differences and Correlation Coefficients across all frequencies and polarizations estimated from overlapping AMSR and MWRI records before and after Tb calibration for five regions representing major Köppen-Geiger global climate classification zones.
Table 4. Summary of sensor ascending (PM) orbit Tb mean biases (Unit: Kelvin), Root Mean Square Error (RMSE) differences and Correlation Coefficients across all frequencies and polarizations estimated from overlapping AMSR and MWRI records before and after Tb calibration for five regions representing major Köppen-Geiger global climate classification zones.
Climate TypeBefore CalibrationAfter Calibration

BiasRMSECorrelationBiasRMSECorrelation
Polar−4.014.650.960.191.110.99
Continental−3.884.080.97−0.180.970.98
Tropical−2.592.880.900.040.460.98
Dry−2.342.740.97−0.010.580.99
Temperate−3.003.180.970.050.580.99
Table 5. Summary of sensor descending (AM) orbit Tb mean biases (Unit: Kelvin), Root Mean Square Error (RMSE) differences and Correlation Coefficients across all frequencies and polarizations estimated from overlapping AMSR and MWRI records before and after Tb calibration for five regions representing major Köppen-Geiger global climate classification zones.
Table 5. Summary of sensor descending (AM) orbit Tb mean biases (Unit: Kelvin), Root Mean Square Error (RMSE) differences and Correlation Coefficients across all frequencies and polarizations estimated from overlapping AMSR and MWRI records before and after Tb calibration for five regions representing major Köppen-Geiger global climate classification zones.
Climate TypeBefore CalibrationAfter Calibration

BiasRMSECorrelationBiasRMSECorrelation
Polar−3.053.920.940.031.120.98
Continental−2.512.830.97−0.040.630.99
Tropical−0.801.520.79−0.310.510.92
Dry−1.141.640.97−0.070.540.99
Temperate−2.012.330.96−0.060.470.99
Table 6. Summary of sensor ascending (PM) orbit Tb mean biases (Unit: Kelvin), Root Mean Square Error (RMSE) differences and Correlation Coefficients across all frequencies and polarizations estimated from overlapping AMSR and MWRI records before and after Tb calibration for major IGBP land cover type areas within the global domain.
Table 6. Summary of sensor ascending (PM) orbit Tb mean biases (Unit: Kelvin), Root Mean Square Error (RMSE) differences and Correlation Coefficients across all frequencies and polarizations estimated from overlapping AMSR and MWRI records before and after Tb calibration for major IGBP land cover type areas within the global domain.
Land Cover TypeBefore CalibrationAfter Calibration

BiasRMSECorrelationBiasRMSECorrelation
ENF *−3.935.600.85−0.493.510.87
EBF *−3.135.450.70−0.023.860.78
DNF *−3.214.300.94−0.042.180.95
DBF *−3.234.590.89−0.082.550.90
Mixed Forest−3.315.360.850.093.400.87
Shrublands *−3.244.810.89−0.152.840.91
Woody Savannas−3.074.530.84−0.182.600.88
Savannas−2.464.080.850.112.390.89
Grasslands−2.824.290.92−0.012.500.93
Croplands−3.215.330.870.033.450.88
Cropland/Natural Vegetation Mosaic−3.024.820.860.083.040.88
Barren or Sparsely Vegetated−2.474.080.83−0.262.480.88

*ENF for Evergreen Needleleaf Forest, EBF for Evergreen Broadleaf Forest, DNF for Deciduous Needleleaf Forest; DBF for Deciduous Broadleaf Forest; Shrublands include both Closed Shrubland and Open Shrubland.

Table 7. Summary of sensor descending (AM) orbit Tb mean biases (Unit: Kelvin), Root Mean Square Error (RMSE) differences and Correlation Coefficients estimated from overlapping AMSR and MWRI records before and after Tb calibration for major IGBP land cover type areas within the global domain.
Table 7. Summary of sensor descending (AM) orbit Tb mean biases (Unit: Kelvin), Root Mean Square Error (RMSE) differences and Correlation Coefficients estimated from overlapping AMSR and MWRI records before and after Tb calibration for major IGBP land cover type areas within the global domain.
Land Cover TypeBefore CalibrationAfter Calibration

BiasRMSECorrelationBiasRMSECorrelation
ENF *−2.924.740.84−0.303.120.88
EBF *−2.134.510.68−0.693.730.74
DNF *−2.693.870.92−0.142.000.94
DBF *−1.833.400.90−0.372.100.92
Mixed Forest−3.154.790.84−0.663.140.88
Shrublands *−2.193.940.88−0.292.540.91
Woody Savannas−1.573.480.79−0.342.460.84
Savannas−1.613.350.84−0.812.560.87
Grasslands−1.823.410.90−0.142.120.92
Croplands−2.484.520.86−0.563.280.88
Cropland/Natural Vegetation Mosaic−2.384.040.86−0.682.820.89
Barren or Sparsely Vegetated−1.383.590.83−0.102.430.88

*ENF for Evergreen Needleleaf Forest, EBF for Evergreen Broadleaf Forest, DNF for Deciduous Needleleaf Forest; DBF for Deciduous Broadleaf Forest; Shrublands include both Closed Shrubland and Open Shrubland.

Table 8. Summary of the accuracy of the satellite based daily minimum and maximum air temperature (Tmn, Tmx) estimates in relation to independent temperature measurements from 142 global WMO weather stations.
Table 8. Summary of the accuracy of the satellite based daily minimum and maximum air temperature (Tmn, Tmx) estimates in relation to independent temperature measurements from 142 global WMO weather stations.
SensorR2 TmxRMSE Tmx (Kelvin)R2 TmnRMSE Tmn(Kelvin)
AMSRE0.863.370.783.46
AMSR20.833.640.753.54
MWRI0.804.120.694.11

Share and Cite

MDPI and ACS Style

Du, J.; Kimball, J.S.; Shi, J.; Jones, L.A.; Wu, S.; Sun, R.; Yang, H. Inter-Calibration of Satellite Passive Microwave Land Observations from AMSR-E and AMSR2 Using Overlapping FY3B-MWRI Sensor Measurements. Remote Sens. 2014, 6, 8594-8616. https://doi.org/10.3390/rs6098594

AMA Style

Du J, Kimball JS, Shi J, Jones LA, Wu S, Sun R, Yang H. Inter-Calibration of Satellite Passive Microwave Land Observations from AMSR-E and AMSR2 Using Overlapping FY3B-MWRI Sensor Measurements. Remote Sensing. 2014; 6(9):8594-8616. https://doi.org/10.3390/rs6098594

Chicago/Turabian Style

Du, Jinyang, John S. Kimball, Jiancheng Shi, Lucas A. Jones, Shengli Wu, Ruijing Sun, and Hu Yang. 2014. "Inter-Calibration of Satellite Passive Microwave Land Observations from AMSR-E and AMSR2 Using Overlapping FY3B-MWRI Sensor Measurements" Remote Sensing 6, no. 9: 8594-8616. https://doi.org/10.3390/rs6098594

Article Metrics

Back to TopTop