Integrated Satellite System for Fire Detection and Prioritization

Several studies have shown the relevance of satellite systems in detecting, monitoring, and characterizing fire events as support to fire management activities. On the other hand, up to now, only a few satellite-based platforms provide immediately and easily usable information about events in progress, in terms of both hotspots, which identify and localize active fires, and the danger conditions of the affected area. However, this kind of information is usually provided through separated layers, without any synthetic indicator which, indeed, could be helpful, if timely provided, for planning the priority of the intervention of firefighting resources in case of concurrent fires. In this study, we try to fill these gaps by presenting an Integrated Satellite System (ISS) for fire detection and prioritization, mainly based on the Robust Satellite Techniques (RST), and the Fire Danger Dynamic Index (FDDI), an original re-structuration of the Índice Combinado de Risco de Incêndio Florestal (ICRIF), for the first time presented here. The system, using Moderate Resolution Imaging Spectroradiometer (MODIS), Advanced Very High Resolution Radiometer (AVHRR), and Spinning Enhanced Visible and InfraRed Imager (SEVIRI) data, provides near real-time integrated information about both the fire presence and danger over the affected area. These satellite-based products are generated in common formats, ready to be ingested in Geographic Information System (GIS) technologies. Results shown and discussed here, on the occasion of concurrent winter and summer fires in Italy, in agreement with information from independent sources, demonstrate that the ISS system, operating at a regional/national scale, may provide an important contribution to fire prioritization. This may result in the mitigation of fire impact in populated areas, infrastructures, and the environment.

Clouds and the Earth’s Radiant Energy System (CERES) FluxByCldTyp Edition 4 Data Product

The Clouds and the Earth’s Radiant Energy System (CERES) project has provided the climate community 20 years of globally observed top of the atmosphere (TOA) fluxes critical for climate and cloud feedback studies. The CERES Flux By Cloud Type (FBCT) product contains radiative fluxes by cloud-type, which can provide more stringent constraints when validating models and also reveal more insight into the interactions between clouds and climate. The FBCT product provides 1° regional daily and monthly shortwave (SW) and longwave (LW) cloud-type fluxes and cloud properties sorted by 7 pressure layers and 6 optical depth bins. Historically, cloud-type fluxes have been computed using radiative transfer models based on observed cloud properties. Instead of relying on radiative transfer models, the FBCT product utilizes Moderate Resolution Imaging Spectroradiometer (MODIS) radiances partitioned by cloud-type within a CERES footprint to estimate the cloud-type broadband fluxes. The MODIS multi-channel derived broadband fluxes were compared with the CERES observed footprint fluxes and were found to be within 1% and 2.5% for LW and SW, respectively, as well as being mostly free of cloud property dependencies. These biases are mitigated by constraining the cloud-type fluxes within each footprint with the CERES Single Scanner Footprint (SSF) observed flux. The FBCT all-sky and clear-sky monthly averaged fluxes were found to be consistent with the CERES SSF1deg product. Several examples of FBCT data are presented to highlight its utility for scientific applications.

Estimating Hourly Full-Coverage PM2.5 Concentrations Based On MODIS Data Over The Northeast of Thailand.

Particulate matter (PM2.5) pollutants are a significant health issue with impacts on human health; however, monitoring of PM2.5 is very limited in developing countries. Satellite remote sensing can expand spatial coverage, potentially enhancing our ability in a specific area for estimating PM2.5; however, some have reported poor predictive performance. An innovative combination of MODIS AOD was developed to fulfill all missing aerosol optical depth (AOD) data obtained from the Moderate Resolution Imaging Spectroradiometer (MODIS). Therefore, hourly PM2.5 concentrations were obtained in Northeastern Thailand. A Linear mixed-effects (LME) model was used to predict location-specific hourly PM2.5 levels. Hourly PM2.5 concentrations measured at 20 PM2.5 monitoring sites and 10- fold cross-validation were addressed for model validation. The observed and predicted concentrations suggested that LME obtained from MODIS AOD data and other factors are a potentially useful predictor of hourly PM2.5 concentrations (R2 >0.70), providing more detailed spatial information for local scales studies. Interestingly, PM2.5 along the Mekong River area was observed higher than in the plain area. The finding can infer that the monsoon wind brings polluted air into the province from sources outside the region. The results will be helpful to analyze air pollution-related health studies.

Analysis of the MODIS above-cloud aerosol retrieval algorithm using MCARS

The Multi-sensor Cloud and Aerosol Retrieval Simulator (MCARS) presently produces synthetic radiance data from Goddard Earth Observing System version 5 (GEOS-5) model output as if the Moderate Resolution Imaging Spectroradiometer (MODIS) were viewing a combination of atmospheric column inclusive of clouds, aerosols, and a variety of gases and land–ocean surface at a specific location. In this paper we use MCARS to study the MODIS Above-Cloud AEROsol retrieval algorithm (MOD06ACAERO). MOD06ACAERO is presently a regional research algorithm able to retrieve aerosol optical thickness over clouds, in particular absorbing biomass-burning aerosols overlying marine boundary layer clouds in the southeastern Atlantic Ocean. The algorithm's ability to provide aerosol information in cloudy conditions makes it a valuable source of information for modeling and climate studies in an area where current clear-sky-only operational MODIS aerosol retrievals effectively have a data gap between the months of June and October. We use MCARS for a verification and closure study of the MOD06ACAERO algorithm. The purpose of this study is to develop a set of constraints a model developer might use during assimilation of MOD06ACAERO data. Our simulations indicate that the MOD06ACAERO algorithm performs well for marine boundary layer clouds in the SE Atlantic provided some specific screening rules are observed. For the present study, a combination of five simulated MODIS data granules were used for a dataset of 13.5 million samples with known input conditions. When pixel retrieval uncertainty was less than 30 %, optical thickness of the underlying cloud layer was greater than 4, and scattering angle range within the cloud bow was excluded, MOD06ACAERO retrievals agreed with the underlying ground truth (GEOS-5 cloud and aerosol profiles used to generate the synthetic radiances) with a slope of 0.913, offset of 0.06, and RMSE=0.107. When only near-nadir pixels were considered (view zenith angle within ±20∘) the agreement with source data further improved (0.977, 0.051, and 0.096 respectively). Algorithm closure was examined using a single case out of the five used for verification. For closure, the MOD06ACAERO code was modified to use GEOS-5 temperature and moisture profiles as an ancillary. Agreement of MOD06ACAERO retrievals with source data for the closure study had a slope of 0.996 with an offset of −0.007 and RMSE of 0.097 at a pixel uncertainty level of less than 40 %, illustrating the benefits of high-quality ancillary atmospheric data for such retrievals.

A Multi-Year Study of GOES-13 Droplet Effective Radius Retrievals for Warm Clouds over South America and Southeast Pacific

Geostationary satellites can retrieve the cloud droplet effective radius (re) but suffer biases from cloud inhomogeneities, internal retrieval nonlinearities, and 3-D scattering/shadowing from neighboring clouds, among others. A 1-D retrieval method was applied to Geostationary Operational Environmental Satellite 13 (GOES-13) imagery, over large areas in South America (+5∘ to −30∘ N5∘ N–30∘ S; −20∘ to −70∘E20∘–70∘ W), the Southeast Pacific (+5∘ to −30∘ N5∘ N–30∘ S; −70∘ to −120∘E70∘–120∘ W), and the Amazon (+2∘ to −7∘ N2∘ N–7∘ S; −54∘ to −73∘E54∘–73∘ W), for four months in each year from 2014–2017. Results were regressedcompared against in situ aircraft measurements and the Moderate Resolution Imaging Spectroradiometer cloud product for Terra and Aqua satellites. Monthly regression parameters approximately followed a seasonal pattern. With up to 108,009 of matchups, slope, intercept, and correlation for Terra (Aqua) ranged from about 0.71 to 1.17, −2.8 to 2.5 μm, and 0.61 to 0.91 (0.54 to 0.78, −1.5 to 1.8 μm, 0.63 to 0.89), respectively. We identified evidence for re overestimation (underestimation) correlated with shadowing (enhanced reflectance) in the forward (backscattering) hemisphere, and limitations to illumination/ and viewing configurations accessible by GOES-13, depending on the time of day and season. A proposition is hypothesized to ameliorate 3-D biases by studying relative illumination and cloud spatial inhomogeneity.

Grounding Event of Iceberg D28 and Its Interactions with Seabed Topography

Iceberg D28, a giant tabular iceberg that calved from Amery Ice Shelf in September 2019, grounded off Kemp Coast, East Antarctica, from August to September of 2020. The motion of the iceberg is characterized herein by time-series images captured by synthetic aperture radar (SAR) on Sentinel-1 and the moderate resolution imaging spectroradiometer (MODIS) boarded on Terra from 6 August to 15 September 2020. The thickness of iceberg D28 was estimated by utilizing data from altimeters on Cryosat-2, Sentinel-3, and ICESat-2. By using the iceberg draft and grounding point locations inferred from its motion, the maximum water depths at grounding points were determined, varying from 221.72 ± 21.77 m to 269.42 ± 25.66 m. The largest disagreements in seabed elevation inferred from the grounded iceberg and terrain models from the Bedmap2 and BedMachine datasets were over 570 m and 350 m, respectively, indicating a more complicated submarine topography in the study area than that presented by the existing seabed terrain models. Wind and sea water velocities from reanalysis products imply that the driving force from sea water is a more dominant factor than the wind in propelling iceberg D28 during its grounding, which is consistent with previous findings on iceberg dynamics.

High altitudinal vegetation dynamics including treeline ecotone in Langtang National Park, Nepal

Identification of high altitudinal vegetation dynamics using remote sensing is important because of the complex topography and environment in the Himalayas. Langtang National Park is the first Himalayan park in Nepal representing the best area to study vegetation change in the central Himalaya region because of the high altitudinal gradient and relatively less disturbed region. This study aimed at mapping vegetation in Langtang National Park and its treeline ecotone using Moderate Resolution Imaging Spectroradiometer (MODIS), Normalized Difference Vegetation Index (NDVI). Two treeline sites with an altitude of 3927 and 3802 meters above sea level (masl) were selected, and species density was measured during the field survey. The linear slope for each pixel and the Mann Kendall test to measure significant trends were used. The results showed that NDVI has significantly increased at the rate of 0.002yr-1 in Langtang National Park and 0.003yr-1 in treeline ecotone during 2000-2017. The average 68.73% equivalents to 1463 km2 of Langtang National Park are covered by vegetation. At the same time, 16.45% equivalents to 350.43 km2 are greening, and 0.25%, i.e., 5.43 km2 are found browning. In treeline ecotone, the vegetation is mostly occupied by grasses, shrublands and small trees where the NDVI was found from 0.1 to 0.5. The relative changes of NDVI in barren lands are negative and vegetative lands above 0.5 NDVI are positive between 2000 and 2017. The dominant treeline vegetation were Abies spectabilis, Rhododendron campanulatum, Betula utilis and Sorbus microphyla, with the vegetation density of 839.28 and 775 individuals per hectare in sites A and B, respectively. The higher average NDVI values, significantly increased NDVI, and higher density of vegetation in both A and B sites indicate that the vegetation in treeline ecotone is obtaining a good environment in the Himalayas of Nepal.

Aerosol Characteristics and Their Impact on the Himalayan Energy Budget

The extensive work on the increasing burden of aerosols and resultant climate implications shows a matter of great concern. In this study, we investigate the aerosol optical depth (AOD) variations in the Indian Himalayan Region (IHR) between its plains and alpine regions and the corresponding consequences on the energy balance on the Himalayan glaciers. For this purpose, AOD data from Moderate Resolution Imaging Spectroradiometer (MODIS, MOD-L3), Aerosol Robotic Network (AERONET), India, and Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observation (CALIPSO) were analyzed. Aerosol radiative forcing (ARF) was assessed using the atmospheric radiation transfer model (RTM) integrated into AERONET inversion code based on the Discrete Ordinate Radiative Transfer (DISORT) module. Further, air mass trajectory over the entire IHR was analyzed using a hybrid Single-Particle Lagrangian Integrated Trajectory (HYSPLIT) model. We estimated that between 2001 and 2015, the monthly average ARF at the surface (ARFSFC), top of the atmosphere (ARFTOA), and atmosphere (ARFATM) were −89.6 ± 18.6 Wm−2, −25.2 ± 6.8 Wm−2, and +64.4 ± 16.5 Wm−2, respectively. We observed that during dust aerosol transport days, the ARFSFC and TOA changed by −112.2 and −40.7 Wm−2, respectively, compared with low aerosol loading days, thereby accounting for the decrease in the solar radiation by 207% reaching the surface. This substantial decrease in the solar radiation reaching the Earth’s surface increases the heating rate in the atmosphere by 3.1-fold, thereby acting as an additional forcing factor for accelerated melting of the snow and glacier resources of the IHR.

Effect of the Partitioning of Diffuse and Direct APAR on GPP Estimation

Accurate estimation of gross primary productivity (GPP) is necessary to better understand the interaction of global terrestrial ecosystems with climate change and human activities. Light use efficiency (LUE)-based GPP models are widely used for retrieving several GPP products with various temporal and spatial resolutions. However, most LUE-based models assume a clear-sky condition, and the influence of diffuse radiation on GPP estimations has not been well considered. In this paper, a diffuse and direct (DDA) absorbed photosynthetically active radiation (APAR)-based method is proposed for better estimation of half-hourly GPP, which partitions APAR under diffuse and direct radiation conditions. Firstly, energy balance residual (EBR) FAPAR, moderate resolution imaging spectroradiometer (MODIS) leaf area index (LAI) (MCD15A2H) and clumping index (CI) products, as well as solar radiation records supplied by FLUXNET2015 were used to calculate diffuse and direct APAR at a half-hourly scale. Then, an eddy covariance-LUE (EC-LUE) model and meteorological observations from FLUXNET2015 data sets were used for obtaining corresponding LUE values. A co-variation relationship between LUE and diffuse fraction was observed, and the LUE was higher under more diffuse radiation conditions. Finally, the DDA-based method was tested using the half-hourly FLUXNET GPP and compared with half-hourly GPP calculated using total APAR (GPP_TA). The results indicated that the half-hourly GPP estimated using the DDA-based method (GPP_DDA) was more accurate, giving higher R2 values, lower RMSE and RMSE* values (R2 varied from 0.565 to 0.682, RMSE ranged from 3.219 to 12.405 and RMSE* were within the range of 2.785 to 8.395) than the GPP_TA (R2 varied from 0.558 to 0.653, RMSE ranged from 3.407 to 13.081 and RMSE* were within the range of 3.321 to 9.625) across FLUXNET sites within different vegetation types. This study explored the effects of partitioning the diffuse and direct APAR on half-hourly GPP estimations, which demonstrates a higher agreement with FLUXNET GPP than total APAR-based GPP.

Estimation of the All-Wave All-Sky Land Surface Daily Net Radiation at Mid-Low Latitudes from MODIS Data Based on ERA5 Constraints

The surface all-wave net radiation (Rn) plays an important role in the energy and water cycles, and most studies of Rn estimations have been conducted using satellite data. As one of the most commonly used satellite data sets, Moderate Resolution Imaging Spectroradiometer (MODIS) data have not been widely used for radiation calculations at mid-low latitudes because of its very low revisit frequency. To improve the daily Rn estimation at mid-low latitudes with MODIS data, four models, including three models built with random forest (RF) and different temporal expansion models and one model built with the look-up-table (LUT) method, are used based on comprehensive in situ radiation measurements collected from 340 globally distributed sites, MODIS top-of-atmosphere (TOA) data, and the fifth generation of European Centre for Medium-Range Weather Forecasts Reanalysis 5 (ERA5) data from 2000 to 2017. After validation against the in situ measurements, it was found that the RF model based on the constraint of the daily Rn from ERA5 (an RF-based model with ERA5) performed the best among the four proposed models, with an overall validated root-mean-square error (RMSE) of 21.83 Wm−2, R2 of 0.89, and a bias of 0.2 Wm−2. It also had the best accuracy compared to four existing products (Global LAnd Surface Satellite Data (GLASS), Clouds and the Earth’s Radiant Energy System Edition 4A (CERES4A), ERA5, and FLUXCOM_RS) across various land cover types and different elevation zones. Further analyses illustrated the effectiveness of the model by introducing the daily Rn from ERA5 into a “black box” RF-based model for Rn estimation at the daily scale, which is used as a physical constraint when the available satellite observations are too limited to provide sufficient information (i.e., when the overpass time is less than twice per day) or the sky is overcast. Overall, the newly-proposed RF-based model with ERA5 in this study shows satisfactory performance and has strong potential to be used for long-term accurate daily Rn global mapping at finer spatial resolutions (e.g., 1 km) at mid-low latitudes.