High resolution limb observations of clouds by the CRISTA-NF experiment during the SCOUT-O3 tropical aircraft campaign

The Cryogenic Infrared Spectrometers and Telescopes for the Atmosphere – New Frontiers (CRISTA-NF) experiment on board the Russian research aircraft Geophysica measures limb emission spectra with an unprecedented vertical and horizontal resolution in the 4–15 μm wavelength region. The IR spectra measured during the SCOUT-O3 Tropical Aircraft Campaign have been analysed with respect of cloud occurrence, cloud vertical and horizontal extent, cloud spatial structures and their utilisation for trace gas retrievals. In addition indicators for ice water content and optical thickness of the clouds have been adopted. These new kinds of measurements in the upper troposphere/lower stratosphere region are especially valuable for the design and development of future space borne high resolution limb sounders.     

A feasibility study for the detection of the diurnal variation of tropospheric NO2 over Tokyo from a geostationary orbit

We have conducted a feasibility study for the geostationary monitoring of the diurnal variation of tropospheric NO₂ over Tokyo. Using NO₂ fields from a chemical transport model, synthetic spectra were created by a radiative transfer model, SCIATRAN, for summer and winter cases. We then performed a Differential Optical Absorption Spectroscopy (DOAS) analysis to retrieve NO₂ slant column densities (SCDs), and after converting SCDs into vertical column densities (VCDs), we estimated the precision of the retrieved VCDs. The simulation showed that signal-to-noise ratio (SNR) ? 500 is needed to detect the diurnal variation and that SNR ? 1000 is needed to observe the local minimum occurring in the early afternoon (LT13–14) in summer. In winter, the detection of the diurnal variation during LT08–15 needs SNR ? 500, and SNR ? 1000 is needed if early morning (LT07) and early evening (LT16) are included. The currently discussed sensor specification for the Japanese geostationary satellite project, GMAP-Asia, which has a horizontal resolution of 10 km and a temporal resolution of 1hr, has demonstrated the performance of a precision of several percent, which is approximately corresponding to SNR = 1000–2000 during daytime and SNR ? 500 in the morning and evening. We also discuss possible biases caused by the temperature dependence of the absorption cross section utilized in the DOAS retrieval, and the effect of uncertainties of surface albedo and clouds on the estimation of precisions.     

Semi-stochastic modification of second-order radial derivative of Abel–Poisson’s formula for validating satellite gravity gradiometry data

The geoid can be used to validate the satellite gravity gradiometry data. Validation of such data is important prior to their downward continuation because of amplification of the data errors through this process. In this paper, the second-order radial derivative of Abel–Poisson’s formula is modified stochastically to reduce the effect of the far-zone geoid and generate the second-order radial derivative of geopotential at 250 km level. The numerical studies over Fennoscandia show that this method yields the gradients with an error of 10 mE and when the long wavelength of geoid is removed from the estimator and restored after the computations (remove–compute–restore) the error will be in 1 mE level. We name this method semi-stochastic modification. The best case scenario is found when the degree of modification of the integral formula is 200 and the long wavelength geoid to degree 100 is removed and restored. In this case the geoid should have a resolution of 15′ × 15′ and the integration should be performed over a cap size of 3°.     

High vertical resolution water vapour profiles in the upper troposphere and lower stratosphere retrieved from MAESTRO solar occultation spectra

An algorithm has been developed that retrieves water vapour profiles in the upper troposphere and lower stratosphere from optical depth spectra obtained by the Measurements of Aerosol Extinction in the Stratosphere and Troposphere Retrieved by Occultation (MAESTRO) instrument onboard the SCISAT satellite as part of the Atmospheric Chemistry Experiment (ACE) mission. The retrieval relies on ro-vibrational absorption of solar radiation by water vapour in the 926–970 nm range. During the iterative inversion process, the optical depth spectra are simulated at the spectral resolution and sampling frequency of MAESTRO using the correlated-k approximation. The Chahine inversion updates the water vapour volume mixing ratio (VMR), adjusting all retrieval layers simultaneously, to match the observed differential optical depth due to absorption by water vapour and ozone at each tangent height. This approach accounts for significant line saturation effects. Profiles are typically obtained from ∼22 km down to the cloud tops or to 5 km, with relative precision as small as 3% in the troposphere. In the lower stratosphere, the precision on water vapour VMR is ∼1.3 μmol/mol in an individual retrieval layer (∼1 km thick). The spectral capability of MAESTRO allows for the clear separation of extinction due to water vapour and aerosol, and for the fitting quality to be quantified and used to determine an altitude-dependent convergence criterion for the retrieval. In the middle troposphere, interhemispheric differences in water vapour VMR are driven by oceanic evaporation whereas in the upper troposphere, deep convection dominates and a strong seasonal cycle is observed at high latitudes.     

Radiometric estimation of water vapor content over Brazil

A multi-channel microwave radiometre (make: Radiometrics Corporation) is installed at Instituto Nacional de Pesquisas Espaciais–INPE, Brazil (22°S). The radiometric output of two channels of the radiometer in the form of brightness temperature at 23.834 GHz and 30 GHz, initially, were used to find out the ambient water vapor content and the non-precipitable cloud liquid water content. The necessary algorithm was developed for the purpose. The best results were obtained using the hinge frequency 23.834 GHz and 30 GHz pair having an r.m.s. error of only 2.64. The same methodology was then adopted exploiting 23.034 GHz and 30 GHz pair. In that case the r.m.s. error was 3.42. These results were then compared with those obtained over Kolkata (22°N), India, by using 22.234 GHz and 31.4 GHz radiometric data. This work conclusively suggests the use of a frequency should not be at the water vapor resonance line. Instead, while measuring the vapor content for separation of vapor and cloud liquid, one of them should be a few GHz left or right from the resonance line i.e., at 23.834 GHz and the other one should be around 30 GHz.     

A comparison of different regression models for downscaling Landsat and MODIS land surface temperature images over heterogeneous landscape

Remotely sensed high spatial resolution thermal images are required for various applications in natural resource management. At present, availability of high spatial resolution (<200 m) thermal images are limited. The temporal resolution of such images is also low. Whereas, coarser spatial resolution (∼1000 m) thermal images with high revisiting capability (∼1 day) are freely available. To bridge this gap, present study attempts to downscale coarser spatial resolution thermal image to finer spatial resolution using relationships between land surface temperature (LST) and vegetation indices over a heterogeneous landscape of India. Five regression based models namely (i) Disaggregation of Radiometric Temperature (DisTrad), (ii) Temperature Sharpening (TsHARP), (iii) TsHARP with local variant, (iv) Least median square regression downscaling (LMS_DS) and (v) Pace regression downscaling (PR_DS) are applied to downscale LST of Landsat Thematic Mapper (TM) and Terra MODIS (Moderate Resolution Imaging Spectroradiometer) images. All the five models are first evaluated on Landsat image aggregated to 960 m resolution and downscaled to 480 m and 240 m resolution. The downscale accuracy is achieved using LMS_DS and PR_DS models at 240 m resolution at 0.61 °C and 0.75 °C respectively. MODIS data downscaled from 1000 m to 250 m spatial resolution results root mean square error (RMSE) of 1.43 °C and 1.62 °C for LMS_DS and PR_DS models, respectively. The LMS_DS model is less sensitive to outliers in heterogeneous landscape and provides higher accuracy when compared to other models. Downscaling model is found to be suitable for agricultural and vegetated landscapes up to a spatial resolution of 250 m but not applicable to water bodies, dry river bed sand sandy open areas.     

Synergetic retrieval of terrestrial AOD from MODIS images of twin satellites Terra and Aqua

Aerosol optical depth (AOD) is one of the most important indicators of atmospheric pollution. It can be retrieved from satellite imagery using several established methods, such as the dark dense vegetation method and the deep blue algorithm. All of these methods require estimation of surface reflectance prior to retrieval, and are applicable to a certain pre-designated type of surface cover. Such limitations can be overcome by using a synergetic method of retrieval proposed in this study. This innovative method is based on the fact that the ratio K of surface reflectance at different angles/geometries is independent of wavelength as reported by Flowerdew and Haigh (1995). An atmospheric radiative transfer model was then established and resolved with the assistance of the ratio K obtained from two Moderate Resolution Imaging Spectroradiometer (MODIS) spectral bands acquired from the twin satellites of Terra and Aqua whose overpass is separated by three hours. This synergetic method of retrieval was tested with 20 pairs of MODIS images. The retrieved AOD was validated against the ground observed AOD at the Taihu station of the AErosol RObotic NETwork (AERONET). It is found that they are correlated with the observations at a coefficient of 0.828 at 0.47 μm and 0.921 at 0.66 μm wavelengths. The retrieved AOD has a mean relative error of 25.47% at 0.47 μm and 24.3% at 0.66 μm. Of the 20 samples, 15 and 17 fall within two standard error of the line based observed AOD data on the ground at the 0.47 μm and 0.66 μm, respectively. These results indicate that this synergetic method can be used to reliably retrieve AOD from the twin satellites MODIS images, namely Terra and Aqua. It is not necessary to determine surface reflectance first.     

Sensitivity analysis for air temperature profile estimation methods around the tropopause using simulated Aqua/AIRS data

A sensitivity analysis is performed to investigate potential improvements to the accuracy of air temperature profile retrievals near the tropopause. A simple inversion method is employed to identify and remove redundant spectral channels from the retrievals using simulated data for the high-spectral resolution sounder AIRS (Atmospheric Infrared Sounder) on the Aqua satellite. Bayesian optimal theory and inverse technique are applied for the atmospheric temperature profile retrievals, and the 15 μm CO₂ absorption bands (620–750 cm⁻¹) are chosen for this study. Sequentially elimination of redundant channels is directly integrated into the inverse scheme for the temperature profile, in order to accurately retain the valuable channels and remove all the redundant channels, for accurate retrieval of the temperature profile. Also, the tropopause and troposphere are treated differently in the inverse scheme to improve the retrieval accuracy in the tropopause. Results of a sensitivity analysis based on this method, for the Tropical and Middle-Latitude Summer models simulated by MODTRAN4.0, show that the estimated accuracies are improved by 2 K around the tropopause, and are only changed by less than 0.2 K in the troposphere.     

Kinetic temperature and carbon dioxide from broadband infrared limb emission measurements taken from the TIMED/SABER instrument

The Sounding of the Atmosphere using Broadband Emission Radiometry (SABER) experiment is one of four instruments on NASA’s Thermosphere–Ionosphere–Energetics and Dynamics (TIMED) satellite. SABER measures broadband infrared limb emission and derives vertical profiles of kinetic temperature (Tk) from the lower stratosphere to approximately 120 km, and vertical profiles of carbon dioxide (CO₂) volume mixing ratio (vmr) from approximately 70 km to 120 km. In this paper we report on SABER Tk/CO₂ data in the mesosphere and lower thermosphere (MLT) region from the version 1.06 dataset. The continuous SABER measurements provide an excellent dataset to understand the evolution and mechanisms responsible for the global two-level structure of the mesopause altitude. SABER MLT Tk comparisons with ground-based sodium lidar and rocket falling sphere Tk measurements are generally in good agreement. However, SABER CO₂ data differs significantly from TIME-GCM model simulations. Indirect CO₂ validation through SABER-lidar MLT Tk comparisons and SABER-radiation transfer comparisons of nighttime 4.3 μm limb emission suggest the SABER-derived CO₂ data is a better representation of the true atmospheric MLT CO₂ abundance compared to model simulations of CO₂ vmr.     

Comparison between African equatorial station ground-based inferred vertical E × B drift,Jicamarca direct measured drift,and IRI model

The Incoherent Scatter Radar measurement over Jicamarca, together with the IRI model-2007 measurements were compared with ground-based digisonde inferred E × B drift over Ilorin in the African region during year of solar minima (F10.7 = 81). Seasonally, Ilorin pre-reversal enhancement (PRE) had peak drift velocities of 7.2, 3.7 and 7.9 m/s for March equinox, September equinox and December solstice respectively, while Jicamarca drifts indicated 13.0, 10.5 and 5.2 m/s; as well as the IRI model with 14.3, 8.4 and 0.7 m/s in similar order. PRE value was insignificant during June solstice. The PRE magnitude of the IRI-model during the equinoxes is twice the value obtained at Ilorin. The daytime E × B drift peaked over Ilorin 1–2 h earlier than both the modeled and Jicamarca observations. This could be due to the difference in sunset time at the conjugate points corresponding to the altitude of the observation. During the evening time PRE, the respective correlation coefficients (R) for Vz–F10.7 relation over Jicamarca, Ilorin and the modeled observations are −0.5559, 0.4796 and −0.4979. Similarly, the Vz–Ap relation exhibit excellent anti-correlation coefficient (R = −0.8637) for the IRI-model, −0.4827 over Jicamarca and 0.3479 for Ilorin. Annual mean drift velocities over Jicamarca, Ilorin and IRI model measurements respectively are 10, 5.6 and 10 m/s for the peak PRE observation; 15, 16 and 21 m/s for the daytime pre-sunrise peak values; and −21, −9 and −16 m/s for the nighttime downward reversals. The root-mean square (RMS) deviation between IRI-model and the Ilorin drift between 2000 and 0500 h is 4.37, 2.03, 3.71 and 2.42 m/s for March equinox, June solstice, September equinox and December solstice respectively. For Jicamarca–Ilorin drift relation, RMS deviation is 5.48, 2.30, 3.47 and 1.27 m/s in the same order respectively. Annual hmF2 inferred drift over Ilorin during daytime is higher by a factor of ≈2 and 3 at Jicamarca and IRI model measurements respectively; and by a factor of ≈5 for both during the night-time period. The limitations in using hmF2 to infer drifts are discussed.     

Analysis of precipitable water vapor from GPS measurements in Chengdu region: Distribution and evolution characteristics in autumn

The rainfall process of Chengdu region in autumn has obvious regional features. Especially, the night-time rain rate of this region in this season is very high in China. Studying the spatial distribution and temporal variation of regional atmospheric precipitable water vapor (PWV) is important for our understanding of water vapor related processes, such as rainfall, evaporation, convective activity, among others in this area. Since GPS detection technology has the unique characteristics, such as all-weather, high accuracy, high spatial and temporal resolution as well as low cost, tracking and monitoring techniques on water vapor has achieved rapid developments in recent years. With GPS–PWV data at 30-min interval gathered from six GPS observational stations in Chengdu region in two autumns (September 2007–December 2007 and September 2008–December 2008), it is revealed that negative correlations exist between seasonally averaged value of GPS–PWV as well as its variation amplitude and local terrain altitude. The variation of PWV in the upper atmosphere of this region results from the water vapor variation from surface to 850 hPa. With the help of Fast Fourier Transform (FFT), it is found that the autumn PWV in Chengdu region has a multi-scale feature, which includes a seasonal cycle, 22.5 days period (quasi-tri-weekly oscillation). The variation of the GPS–PWV is related to periodical change in the transmitting of the water vapor caused by zonal and meridional wind strengths’ change and to the East Asian monsoon system. According to seasonal variation characteristics, we concluded that the middle October is the critical turning point in PWV content. On a shorter time scale, the relationship between autumn PWV and ground meteorological elements was obtained using the composite analysis approach.     

Study of daytime vertical E × B drift velocities inferred from ground-based magnetometer observations of ΔH,at low latitudes under geomagnetically disturbed conditions

In this study, 30 storm sudden commencement (SSC) events during the period 2001–2007 for which daytime vertical E × B drift velocities from JULIA radar, Jicamarca (geographic latitude 11.91°S, geographic longitude 283.11°E, 0.81°N dip latitude), Peru and ΔH component of geomagnetic field measured as the difference between the magnitudes of the horizontal (H) components between two magnetometers deployed at two different locations Jicamarca (geographic latitude 11.91°S, geographic longitude 283.11°E, 0.81°N dip latitude) and Piura (geographic latitude 5.21°S, geographic longitude 279.41°E, 6.81°N dip latitude), in Peru, were considered. It is observed that a positive correlation exists between peak value of daytime vertical E × B drift velocity and peak value of ΔH for the three consecutive days of SSC. A qualitative analysis made after selecting the peak values of daytime vertical E × B drift velocity and ΔH showed that 57% of the events have daytime vertical E × B drift velocity peak in the magnitude range 20–30 m/s and 63% of the events have ΔH peak in the range 80–100 nT. The maximum probable (45%) range of time of occurrence of peak value for both vertical E × B drift velocity and ΔH during the daytime hours were found to be the same, i.e., 10:00–12:00 LT. A strong positive correlation was also found to exist between the daytime vertical E × B drift velocity and ΔH for all the three consecutive days of SSC, for all the events considered. To establish a quantitative relationship between day time vertical E × B drift velocity and ΔH, linear and polynomial (order 2 and 3) regression analysis (Least Square Method (LSM)) were carried out, considering the fully disturbed day after the commencement of the storm as ‘disturbed period’ for the SSC events selected for analysis. The formulae indicating the relationship between daytime vertical E × B drift velocity and ΔH, for the ‘disturbed periods’, obtained through the regression analysis were verified using the JULIA radar observed E × B drift velocity for 3 selected events. Root Mean Square (RMS) error analysis carried out for each case suggest that polynomial regression (order 3) analysis provides a better agreement with the observations from among the linear, polynomial (order 2 and 3) analysis.     

A new window on the cosmos: The Stratospheric Observatory for Infrared Astronomy (SOFIA)

The Stratospheric Observatory for Infrared Astronomy (SOFIA) is a joint US/German Project to develop and operate a gyrostabilized 2.5-m telescope in a Boeing 747-SP. This observatory will allow astronomical observations from 0.3 μm to sub-millimeter wavelengths at stratospheric altitudes as high as 45,000 ft where the atmosphere is not only cloud-free, but largely transparent at infrared wavelengths. The dynamics and chemistry of interstellar matter, and the details of embedded star formation will be key science goals. In addition, SOFIA’s unique portability will enable large-telescope observations at sites required to observe transient phenomena and location specific events. SOFIA will offer the convenient accessibility of a ground-based telescope for servicing, maintenance, and regular technology upgrades, yet will also have many of the performance advantages of a space-based telescope. Initially, SOFIA will fly with nine first-generation focal plane instruments that include broad-band imagers, moderate resolution spectrographs that will resolve broad features from dust and large molecules, and high resolution spectrometers capable of studying the chemistry and detailed kinematics of molecular and atomic gas. First science flights will begin in 2010, leading to a full operations schedule of about 120 8–10 h flights per year by 2014. The next call for instrument development that can respond to scientifically exciting new technologies will be issued in 2010. We describe the SOFIA facility and outline the opportunities for observations by the general scientific community with cutting edge focal plane technology. We summarize the operational characteristics of the first-generation instruments and give specific examples of the types of fundamental scientific studies these instruments are expected to make.     

Total electron content monitoring using triple frequency GNSS: Results with Giove-A/-B data

Triple frequency GNSS will be fully operational within the next decade, opening opportunities for new applications. Dual frequency GNSS already allow to study the ionosphere through the estimation of Total Electron Content (TEC). However, the precision is limited by the ambiguity resolution process. This paper studies a triple frequency TEC monitoring technique in which the use of Geometry-Free and Iono-Free linear combinations improves the ambiguity resolution process and therefore the precision of TEC. We have tested it on a set of triple frequency Giove-A/-B data from January and December 2008. The conclusions achieved are (1) TEC values are affected by an error of about 2–2.5 TECU produced through the ambiguity resolution process; (2) the error caused by the Geometric Free phase combination delays (hardware, multipath, noise, antenna phase center) on TEC is about 0.2 TECU; (3) the total error on TEC approximately reach 2–3 TECU.     

Quantum dosimetry and online visualization of X-ray and charged particle radiation in commercial aircraft at operational flight altitudes with the pixel detector Timepix

We investigate the application of the hybrid semiconductor pixel detector Timepix for precise characterization, quantum sensitivity dosimetry and visualization of the charged particle radiation and X-ray field inside commercial aircraft at operational flight altitudes. The quantum counting capability and granularity of Timepix provides the composition and spectral-characteristics of the X-ray and charged-particle field with high sensitivity, wide dynamic range, high spatial resolution and particle type resolving power. For energetic charged particles the direction of trajectory and linear energy transfer can be measured. The detector is operated by the integrated readout interface FITPix for power, control and data acquisition together with the software package Pixelman for online visualization and real-time data processing. The compact and portable radiation camera can be deployed remotely being controlled simply by a laptop computer. The device performs continuous monitoring and accurate time-dependent measurements in wide dynamic range of particle fluxes, deposited energy, absorbed dose and equivalent dose rates. Results are presented for in-flight measurements at altitudes up to 12 km in various flights selected in the period 2006–2013.     

Mesospheric doppler wind measurements from Aura Microwave Limb Sounder (MLS)

This paper describes a microwave limb technique for measuring Doppler wind in the Earth’s mesosphere. The research algorithm has been applied to Aura Microwave Limb Sounder (MLS) 118.75 GHz measurements where the O₂ Zeeman lines are resolved by a digital autocorrelation spectrometer. A precision of ∼17 m/s for the line-of-sight (LOS) wind is achieved at 80–92 km, which corresponds to radiometric noise during 1/6 s integration time. The LOS winds from Aura MLS are mostly in the meridional direction at low- and mid-latitudes with vertical resolution of ∼8 km. This microwave Doppler technique has potential to obtain useful winds down to ∼40 km of the Earth’s atmosphere if measurements from other MLS frequencies (near H₂O, O₃, and CO lines) are used. Initial analyses show that the MLS winds from the 118.75 GHz measurements agree well with the TIDI (Thermosphere Ionosphere Mesosphere Energetics and Dynamics Doppler Interferometer) winds for the perturbations induced by a strong quasi 2-day wave (QTDW) in January 2005. Time series of MLS winds reveal many interesting climatological and planetary wave features, including the diurnal, semidiurnal tides, and the QTDW. Interactions between the tides and the QTDW are clearly evident, indicating possible large tidal structural changes after the QTDW events dissipate.     

Equatorial vertical E × B drift velocities inferred from ionosonde measurements over Ouagadougou and the IRI-2007 vertical ion drift model

F-region vertical plasma drift velocities were deduced from the hourly hmF2 values acquired from ionogram data over a near dip equatorial station Ouagadougou (12.4°N, 358.5°E, dip angle 5.9°N) in Africa. Our results are compared against the global empirical model of Scherliess and Fejer (1999) incorporated in the IRI model (IRI-2007) for 1600 to 0800 LT from 1 year of data during sunspot maximum year of 1989 (yearly average solar flux intensity, F10.7 = 192) corresponding to the peak phase of solar cycle 22, under magnetically quiet conditions. The drifts are entirely downward between 2000 and 0500 LT bin for both techniques and the root mean square error (RMSE) between the modeled and the ionosonde vertical plasma drifts during these periods is 3.80, 4.37, and 4.74 m/s for June solstice, December solstice and equinox, respectively. Ouagadougou average vertical drifts show evening prereversal enhancement (PRE) velocity peaks (V_ZP) of about 16, 14, and 17 m/s in June solstice, December solstice, and equinox, respectively, at 1900–2000 LT; whereas global empirical model average drifts indicate V_ZP of approximately 33 m/s (June solstice), 29 m/s (December solstice), and 50 m/s (equinox) at 1800 LT. We find very weak and positive correlation (+0.10376) between modeled V_ZP versus F10.7, while ionosonde V_ZP against F10.7 gives worst and opposite correlation (−0.05799). The results also show that modeled V_ZP–Ap indicates good and positive correlation (+0.64289), but ionosonde V_ZP–Ap exhibits poor and negative correlation (−0.22477).     

Carbon monoxide spatial gradients over source regions as observed by SCIAMACHY: A case study for the United Kingdom

Carbon monoxide (CO) is an important air pollutant whose emissions and atmospheric concentrations need to be monitored. The measurements of the SCIAMACHY instrument on ENVISAT are sensitive to CO concentration changes at all atmospheric altitude levels including the boundary layer. The SCIAMACHY CO measurements therefore contain information on CO emissions. Until now no studies have been published where the SCIAMACHY CO measurements have been used to quantify CO emissions by applying, for example, inverse modelling approaches. Here we report about a step in this direction. We have analysed three years of CO columns to investigate if spatial gradients resulting from United Kingdom (UK) CO emissions can be observed from space. The UK is an interesting target area because the UK is a relatively well isolated CO source region. On the other hand, the UK is not the easiest target as its emissions are only moderate and because the surrounding water has low reflectivity in the 2.3 μm spectral region used for CO retrieval. We determined horizontal CO gradients from seasonally and yearly averaged CO during 2003–2005 over the UK taking into account daily wind fields. We show that the measured CO longitudinal (downwind) gradients have the expected order of magnitude. The estimated 2σ error of the gradients depends on time period and applied filtering criteria (e.g., land only, cloud free) and is typically 10–20% of the total column. The gradients are barely statistically significant within the 2σ error margin. This is mainly because of the relatively high noise of the SCIAMACHY CO measurements in combination with a quite low number of measurements (∼100) mainly due to cloud cover.     

Evaluation of MODIS GPP over a complex ecosystem in East Asia: A case study at Gwangneung flux tower in Korea

Moderate Resolution Imaging Radiometer (MODIS) gross primary productivity (GPP) has been used widely to study the global carbon cycle associated with terrestrial ecosystems. The retrieval of the current MODIS productivity with a 1 × 1 km² resolution has limitations when presenting subgrid scale processes in terrestrial ecosystems, specifically when forests are located in mountainous areas, and shows heterogeneity in vegetation type due to intensive land use. Here, we evaluate MODIS GPP (MOD17) at Gwangneung deciduous forest KoFlux tower (deciduous forest; GDK) for 2006–2010 in Korea, where the forests comprise heterogeneous vegetation cover over complex terrain. The monthly MODIS GPP data overestimated the GDK measurements in a range of +15% to +34% and was fairly well correlated (R = 0.88) with the monthly variability at GDK during the growing season. In addition, the MODIS data partly represented the sharp GPP reduction during the Asian summer monsoon (June–September) when intensive precipitation considerably reduces solar radiation and disturbs the forest ecosystem. To examine the influence of subgrid scale heterogeneity on GPP estimates over the MODIS scale, the individual vegetation type and its area within a corresponding MODIS pixel were identified using a national forest type map (∼71-m spatial resolution), and the annual GPP in the same area as the MODIS pixel was estimated. This resulted in a slight reduction in the positive MODIS bias by ∼10%, with a high degree of uncertainty in the estimation. The MODIS discrepancy for GDK suggests further investigation is necessary to determine the MODIS errors associated with the site-specific aerodynamic and hydrological characteristics that are closely related to the mountainous topography. The accuracy of meteorological variables and the impact of the very cloudy conditions in East Asia also need to be assessed.     

Radiance calibration of CRISTA-NF

The CRISTA-NF instrument is the airborne version of the CRISTA satellite infrared limb sounder. It has been successfully flown on the Geophysica research airplane during a test campaign in July 2005, during the SCOUT-O3 Tropical Aircraft Campaign in November/December 2005 and during the AMMA campaign in August 2006. Radiance calibrations of the airborne instrument are more complex compared to the satellite instrument because the vacuum shell of CRISTA-NF is confined by a ZnSe (zinc–selenide) window and the detectors can thermally drift during measurement flights. By comprehensive radiance calibrations with a blackbody source the window’s emissivity and transmissivity are determined and the dependence of the instrument sensitivity on the detector temperature is characterized. Taking these effects into account, the remaining radiance error of the calibration is smaller than 3%.