Systematic and random errors between collocated satellite ice water path observations

There remains large disagreement between ice-water path (IWP) in observational data sets, largely because the sensors observe different parts of the ice particle size distribution. A detailed comparison of retrieved IWP from satellite observations in the Tropics (!30 " latitude) in 2007 was made using collocated measurements. The radio detection and ranging(radar)/light detection and ranging (lidar) (DARDAR) IWP data set, based on combined radar/lidar measurements, is used as a reference because it provides arguably the best estimate of the total column IWP. For each data set, usable IWP dynamic ranges are inferred from this comparison. IWP retrievals based on solar reflectance measurements, in the moderate resolution imaging spectroradiometer (MODIS), advanced very high resolution radiometer–based Climate Monitoring Satellite Applications Facility (CMSAF), and Pathfinder Atmospheres-Extended (PATMOS-x) datasets, were found to be correlated with DARDAR over a large IWP range (~20–7000 g m -2 ). The random errors of the collocated data sets have a close to lognormal distribution, and the combined random error of MODIS and DARDAR is less than a factor of 2, which also sets the upper limit for MODIS alone. In the same way, the upper limit for the random error of all considered data sets is determined. Data sets based on passive microwave measurements, microwave surface and precipitation products system (MSPPS), microwave integrated retrieval system (MiRS), and collocated microwave only (CMO), are largely correlated with DARDAR for IWP values larger than approximately 700 g m -2 . The combined uncertainty between these data sets and DARDAR in this range is slightly less MODIS-DARDAR, but the systematic bias is nearly an order of magnitude.

The impact of a thermodynamic sea-ice module in the COSMO numerical weather prediction model on simulations for the Laptev Sea, Siberian Arctic

Previous versions of the Consortium for Small-scale Modelling (COSMO) numerical weather prediction model have used a constant sea-ice surface temperature, but observations show a high degree of variability on sub-daily timescales. To account for this, we have implemented a thermodynamic sea-ice module in COSMO and performed simulations at a resolution of 15 km and 5 km for the Laptev Sea area in April 2008. Temporal and spatial variability of surface and 2-m air temperature are verified by four automatic weather stations deployed along the edge of the western New Siberian polynya during the Transdrift XIII-2 expedition and by surface temperature charts derived from Moderate Resolution Imaging Spectroradiometer (MODIS) satellite data. A remarkable agreement between the new model results and these observations demonstrates that the implemented sea-ice module can be applied for short-range simulations. Prescribing the polynya areas daily, our COSMO simulations provide a high-resolution and high-quality atmospheric data set for the Laptev Sea for the period 14-30 April 2008. Based on this data set, we derive a mean total sea-ice production rate of 0.53 km3/day for all Laptev Sea polynyas under the assumption that the polynyas are ice-free and a rate of 0.30 km3/day if a 10-cm-thin ice layer is assumed. Our results indicate that ice production in Laptev Sea polynyas has been overestimated in previous studies.

Evaluation of a polynya flux model by means of thermal infrared satellite estimates

We test the ability of a two-dimensional flux model to simulate polynya events with narrow open-water zones by comparing model results to ice-thickness and ice-production estimates derived from thermal infrared Moderate Resolution Imaging Spectroradiometer (MODIS) observations in conjunction with an atmospheric dataset. Given a polynya boundary and an atmospheric dataset, the model correctly reproduces the shape of an 11 day long event, using only a few simple conservation laws. Ice production is slightly overestimated by the model, owing to an underestimated ice thickness. We achieved best model results with the consolidation thickness parameterization developed by Biggs and others (2000). Observed regional discrepancies between model and satellite estimates might be a consequence of the missing representation of the dynamic of the thin-ice thickening (e.g. rafting). We conclude that this simplified polynya model is a valuable tool for studying polynya dynamics and estimating associated fluxes of single polynya events.

Regional representativity of AERONET observation sites during the biomass burning season in South America determined by correlation studies with MODIS Aerosol Optical Depth

This paper presents an analysis of ground-based Aerosol Optical Depth (AOD) observations by the Aerosol Robotic Network (AERONET) in South America from 2001 to 2007 in comparison with the satellite AOD product of Moderate Resolution Imaging Spectroradiometer (MODIS), aboard TERRA and AQUA satellites. Data of 12 observation sites were used with primary interest in AERONET sites located in or downwind of areas with high biomass burning activity and with measurements available for the full time range. Fires cause the predominant carbonaceous aerosol emission signal during the dry season in South America and are therefore a special focus of this study. Interannual and seasonal behavior of the observed AOD at different sites were investigated, showing clear differences between purely fire and urban influenced sites. An intercomparison of AERONET and MODIS AOD annual correlations revealed that neither an interannual long-term trend may be observed nor that correlations differ significantly owing to different overpass times of TERRA and AQUA. Individual anisotropic representativity areas for each AERONET site were derived by correlating daily AOD of each site for all years with available individual MODIS AOD pixels gridded to 1 degrees x 1 degrees. Results showed that for many sites a good AOD correlation (R(2) > 0.5) persists for large, often strongly anisotropic, areas. The climatological areas of common regional aerosol regimes often extend over several hundreds of kilometers, sometimes far across national boundaries. As a practical application, these strongly inhomogeneous and anisotropic areas of influence are being implemented in the tropospheric aerosol data assimilation system of the Coupled Chemistry-Aerosol-Tracer Transport Model coupled to the Brazilian Regional Atmospheric Modeling System (CCATT-BRAMS) at the Brazilian National Institute for Space Research (INPE). This new information promises an improved exploitation of local site sampling and, thus, chemical weather forecast.

Geoprocessamento de dados meteo-oceonográficos (cor do oceano e temperatura da superfície do mar) aplicado ao monitoramento ambiental na costa setentrional do Rio Grande do Norte

Orbital remote sensing has been used as a beneficial tool in improving the knowledge on oceanographic and hydrodynamic aspects in northern portion of the continental shelf of Rio Grande do Norte, offshore Potiguar Basin. Aspects such as geography, temporal and spatial resolution combined with a consistent methodology and provide a substantial economic advantage compared to traditional methods of in situ data collecting. Images of the Moderate Resolution Imaging Spectroradiometer (MODIS) sensor onboard NASA's AQUA satellite were obtained to support systematic data collections related to the campaign of environmental monitoring and characterization of Potiguar Basin, held in May 2004. Images of Total Suspension Matter (TSM) and values of radiance standard were generated for the calculation of concentrations of total suspension matter (TSM), chlorophyll-a and sea surface temperature (SST). These data sets were used for statistical comparisons between measures in situ and satellite estimates looking validate algorithms or develop a comprehensive regional approach empirically. AQUA-MODIS images allowed the simultaneous comparison of two-dimensional water quality (total suspension matter), phytoplankton biomass (chlorophyll-a) variability and physical (temperature). For images of total suspension matter, the generated models showed a good correlation with the field data, allowing quantitative and qualitative analysis. The images of chlorophyll-a showed a consistent correlation with the in situ values of concentration. The algorithms adjusted for these images obtained a correlation coefficient fairly well with the data field in order that the sensor can be having an effect throughout the water column and not just the surface. This has led to a fit between the data of chlorophyll-the integration of the average sampling interval of the entire water column up to the level of the first optical depth, with the data generated from the images. This method resulted in higher values of chlorophyll concentration to greater depths, due to the fact that we are integrating more values of chlorophyll in the water column. Thus we can represent the biomass available in the water column. Images SST and SST measures in situ showed a mean difference DT (SST insitu - SST sat) around -0.14 ° C, considered low, making the results very good. The integration of total suspension matter, chlorophyll-a, the temperature of the sea surface (SST) and auxiliary data enabled the recognition of some of the main ways to fund the continental shelf. The main features highlighted were submerged canyons of rivers Apodi and Açu, some of the lines and beachrocks reefs, structural highs and the continental shelf break which occurs at depths around -60 m. The results confirmed the high potential for use of the AQUA-MODIS images to environmental monitoring of sea areas due to ease of detection of the field two-dimensional material in suspension on the sea surface, temperature and the concentration of chlorophyll-a

Large-Scale Water Productivity Assessments with MODIS Images in a Changing Semi-Arid Environment: A Brazilian Case Study

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Análise do potencial do sensor OLI (Operational Land Imager) a bordo do Landsat-8 em estimar a concentração de sólidos totais em suspensão no reservatório de Barra Bonita

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Recent shift from forest to savanna burning in the Amazon Basin observed by satellite

The numbers of fires detected on forest, savanna and transition lands during the 2002-10 biomass burning seasons in Amazonia are shown using fire count data and co-located land cover classifications from the Moderate Resolution Imaging Spectroradiometer (MODIS). The ratio of forest fires to savanna fires has varied substantially over the study period, with a maximum ratio of 0.65:1 in 2005 and a minimum ratio of 0.27:1 in 2009, with the four lowest years occurring in 2007-10. The burning during the droughts of 2007 and 2010 is attributed to a higher number of savanna fires relative to the drought of 2005. A decrease in the regional mean single scattering albedo of biomass burning aerosols, consistent with the shift from forest to savanna burning, is also shown. During the severe drought of 2010, forest fire detections were lower in many areas compared with 2005, even though the drought was more severe in 2010. This result suggests that improved fire management practices, including stricter burning regulations as well as lower deforestation burning, may have reduced forest fires in 2010 relative to 2005 in some areas of the Amazon Basin.

Calibration of a Species-Specific Spectral Vegetation Index for Leaf Area Index (LAI) Monitoring: Example with MODIS Reflectance Time-Series on Eucalyptus Plantations

The leaf area index (LAI) is a key characteristic of forest ecosystems. Estimations of LAI from satellite images generally rely on spectral vegetation indices (SVIs) or radiative transfer model (RTM) inversions. We have developed a new and precise method suitable for practical application, consisting of building a species-specific SVI that is best-suited to both sensor and vegetation characteristics. Such an SVI requires calibration on a large number of representative vegetation conditions. We developed a two-step approach: (1) estimation of LAI on a subset of satellite data through RTM inversion; and (2) the calibration of a vegetation index on these estimated LAI. We applied this methodology to Eucalyptus plantations which have highly variable LAI in time and space. Previous results showed that an RTM inversion of Moderate Resolution Imaging Spectroradiometer (MODIS) near-infrared and red reflectance allowed good retrieval performance (R-2 = 0.80, RMSE = 0.41), but was computationally difficult. Here, the RTM results were used to calibrate a dedicated vegetation index (called "EucVI") which gave similar LAI retrieval results but in a simpler way. The R-2 of the regression between measured and EucVI-simulated LAI values on a validation dataset was 0.68, and the RMSE was 0.49. The additional use of stand age and day of year in the SVI equation slightly increased the performance of the index (R-2 = 0.77 and RMSE = 0.41). This simple index opens the way to an easily applicable retrieval of Eucalyptus LAI from MODIS data, which could be used in an operational way.

PM(10) remote sensing from geostationary SEVIRI and polar-orbiting MODIS sensors over the complex terrain of the European Alpine region

The subject of this study is to investigate the capability of spaceborne remote sensing data to predict ground concentrations of PM10 over the European Alpine region using satellite derived Aerosol Optical Depth (AOD) from the geostationary Spinning Enhanced Visible and InfraRed Imager (SEVIRI) and the polar-orbiting MODerate resolution Imaging Spectroradiometer (MODIS). The spatial and temporal resolutions of these aerosol products (10 km and 2 measurements per day for MODIS, ∼ 25 km and observation intervals of 15 min for SEVIRI) permit an evaluation of PM estimation from space at different spatial and temporal scales. Different empirical linear relationships between coincident AOD and PM10 observations are evaluated at 13 ground-based PM measurement sites, with the assumption that aerosols are vertically homogeneously distributed below the planetary Boundary Layer Height (BLH). The BLH and Relative Humidity (RH) variability are assessed, as well as their impact on the parameterization. The BLH has a strong influence on the correlation of daily and hourly time series, whilst RH effects are less clear and smaller in magnitude. Despite its lower spatial resolution and AOD accuracy, SEVIRI shows higher correlations than MODIS (rSEV∼ 0.7, rMOD∼ 0.6) with regard to daily averaged PM10. Advantages from MODIS arise only at hourly time scales in mountainous locations but lower correlations were found for both sensors at this time scale (r∼ 0.45). Moreover, the fraction of days in 2008 with at least one satellite observation was 27% for SEVIRI and 17% for MODIS. These results suggest that the frequency of observations plays an important role in PM monitoring, while higher spatial resolution does not generally improve the PM estimation. Ground-based Sun Photometer (SP) measurements are used to validate the satellite-based AOD in the study region and to discuss the impact of aerosols' micro-physical properties in the empirical models. A lower error limit of 30 to 60% in the PM10 assessment from space is estimated in the study area as a result of AOD uncertainties, variability of aerosols properties and the heterogeneity of ground measurement sites. It is concluded that SEVIRI has a similar capacity to map PM as sensors on board polar-orbiting platforms, with the advantage of a higher number of observations. However, the accuracy represents a serious limitation to the applicability of satellites for ground PM mapping, especially in mountainous areas.

Validation of a modified AVHRR aerosol optical depth retrieval algorithm over Central Europe

The Advanced Very High Resolution Radiometer (AVHRR) carried on board the National Oceanic and Atmospheric Administration (NOAA) and the Meteorological Operational Satellite (MetOp) polar orbiting satellites is the only instrument offering more than 25 years of satellite data to analyse aerosols on a daily basis. The present study assessed a modified AVHRR aerosol optical depth τa retrieval over land for Europe. The algorithm might also be applied to other parts of the world with similar surface characteristics like Europe, only the aerosol properties would have to be adapted to a new region. The initial approach used a relationship between Sun photometer measurements from the Aerosol Robotic Network (AERONET) and the satellite data to post-process the retrieved τa. Herein a quasi-stand-alone procedure, which is more suitable for the pre-AERONET era, is presented. In addition, the estimation of surface reflectance, the aerosol model, and other processing steps have been adapted. The method's cross-platform applicability was tested by validating τa from NOAA-17 and NOAA-18 AVHRR at 15 AERONET sites in Central Europe (40.5° N–50° N, 0° E–17° E) from August 2005 to December 2007. Furthermore, the accuracy of the AVHRR retrieval was related to products from two newer instruments, the Medium Resolution Imaging Spectrometer (MERIS) on board the Environmental Satellite (ENVISAT) and the Moderate Resolution Imaging Spectroradiometer (MODIS) on board Aqua/Terra. Considering the linear correlation coefficient R, the AVHRR results were similar to those of MERIS with even lower root mean square error RMSE. Not surprisingly, MODIS, with its high spectral coverage, gave the highest R and lowest RMSE. Regarding monthly averaged τa, the results were ambiguous. Focusing on small-scale structures, R was reduced for all sensors, whereas the RMSE solely for MERIS substantially increased. Regarding larger areas like Central Europe, the error statistics were similar to the individual match-ups. This was mainly explained with sampling issues. With the successful validation of AVHRR we are now able to concentrate on our large data archive dating back to 1985. This is a unique opportunity for both climate and air pollution studies over land surfaces.

Lava discharge rate estimates from thermal infrared satellite data at Pacaya Volcano, Guatemala

Time-averaged discharge rates (TADR) were calculated for five lava flows at Pacaya Volcano (Guatemala), using an adapted version of a previously developed satellite-based model. Imagery acquired during periods of effusive activity between the years 2000 and 2010 were obtained from two sensors of differing temporal and spatial resolutions; the Moderate Resolution Imaging Spectroradiometer (MODIS), and the Geostationary Operational Environmental Satellites (GOES) Imager. A total of 2873 MODIS and 2642 GOES images were searched manually for volcanic “hot spots”. It was found that MODIS imagery, with superior spatial resolution, produced better results than GOES imagery, so only MODIS data were used for quantitative analyses. Spectral radiances were transformed into TADR via two methods; first, by best-fitting some of the parameters (i.e. density, vesicularity, crystal content, temperature change) of the TADR estimation model to match flow volumes previously estimated from ground surveys and aerial photographs, and second by measuring those parameters from lava samples to make independent estimates. A relatively stable relationship was defined using the second method, which suggests the possibility of estimating lava discharge rates in near-real-time during future volcanic crises at Pacaya.

Victorian Dominant Land Cover, 2009-2013

There is a demand for regularly updated, broad-scale, accurate land cover information in Victoria from multiple stakeholders. This paper documents the methods used to generate an annual dominant land cover (DLC) map for Victoria, Australia from 2009 to 2013. Vegetation phenology parameters derived from an annual time series of the Moderate Resolution Imaging Spectroradiometer Vegetation Indices 16-day 250 m (MOD13Q1) product were used to generate annual DLC maps, using a three-tiered hierarchical classification scheme. Classification accuracy at the broadest (primary) class level was over 91% for all years, while it ranged from 72 to 81% at the secondary class level. The most detailed class level (tertiary) had accuracy levels ranging from 61 to 68%. The approach used was able to accommodate variable climatic conditions, which had substantial impacts on vegetation growth patterns and agricultural production across the state between both regions and years. The production of an annual dataset with complete spatial coverage for Victoria provides a reliable base data set with an accuracy that is fit-for-purpose for many applications.

Chlorophyll a concentration, and proportion of time in the water column obtained by instrumented bowhead whales

Sixty hours of direct measurements of fluorescence were collected from six bowhead whales (Balaena mysticetus) instrumented with fluorometers in Greenland in April 2005 and 2006. The data were used to (1) characterize the three-dimensional spatial pattern of chlorophyll-a (Chl-a) in the water column, (2) to examine the relationships between whale foraging areas and productive zones, and (3) to examine the correlation between whale-derived in situ values of Chl-a and those from concurrent satellite images using the NASA MODIS (Moderate Resolution Imaging Spectroradiometer) EOS-AQUA satellite (MOD21, SeaWifs analogue OC3M and SST MOD37). Bowhead whales traversed 1600 km**2, providing information on diving, Chl-a structure and temperature profiles to depths below 200 m. Feeding dives frequently passed through surface waters ( >50 m) and targeted depths close to the bottom, and whales did not always target patches of high concentrations of Chl-a in the upper 50 m. Five satellite images were available within the periods whales carried fluorometers. Whales traversed 91 pixels collecting on average 761 s (SD 826) of Chl-a samples per pixel (0-136 m). The depth of the Chl-a maximum ranged widely, from 1 to 66 m. Estimates of Chl-a made from the water-leaving radiance measurements using the OC3M algorithm were highly skewed with most samples estimated as <1 mg/m**3 Chl-a, while data collected from whales had a broad distribution with Chl-a reaching >9 mg/m**3. The correlation between the satellite-derived and whale-derived Chl-a maxima was poor, a linear fit explained only 10% of the variance.

Daily MODIS composites of thin-ice thickness and ice-surface temperatures for the Southern Weddell Sea

Based upon high-resolution thermal-infrared Moderate-Resolution Imaging Spectroradiometer (MODIS) satellite imagery in combination with ERA-Interim atmospheric reanalysis data, we derived long-term polynya parameters such as polynya area, thin-ice thickness distribution and ice-production rates from daily cloud-cover corrected thin-ice thickness composites. Our study is based on a thirteen year investigation period (2002-2014) for the austral winter (1 April to 30 September) in the Antarctic Southern Weddell Sea. The focus lies on coastal polynyas which are important hot spots for new-ice formation, bottom-water formation and heat/moisture release into the atmosphere. MODIS has the capability to resolve even very narrow coastal polynyas. Its major disadvantage is the sensor limitation due to cloud cover. We make use of a newly developed and adapted spatial feature reconstruction scheme to account for cloud-covered areas. We find the sea-ice areas in front of Ronne and Brunt Ice Shelf to be the most active with an annual average polynya area of 3018 ± 1298 and 3516 ± 1420 km2 as well as an accumulated volume ice production of 31 ± 13 and 31 ± 12 km**3, respectively. For the remaining four regions, estimates amount to 421 ± 294 km**2 and 4 ± 3 km**3 (Antarctic Peninsula), 1148 ± 432 km**2 and 12 ± 5 km**3 (Iceberg A23A), 901 ± 703 km**2 and 10 ± 8 km**3 (Filchner Ice Shelf) as well as 499 ± 277 km**2 and 5 ± 2 km**3 (Coats Land). Our findings are discussed in comparison to recent studies based on coupled sea-ice/ocean models and passive-microwave satellite imagery, each investigating different parts of the Southern Weddell Sea.