skip to main content


Title: Intraurban NO 2 hotspot detection across multiple air quality products
Abstract

High-resolution air quality data products have the potential to help quantify inequitable environmental exposures over space and across time by enabling the identification of hotspots, or areas that consistently experience elevated pollution levels relative to their surroundings. However, when different high-resolution data products identify different hotspots, the spatial sparsity of ‘gold-standard’ regulatory observations leaves researchers, regulators, and concerned citizens without a means to differentiate signal from noise. This study compares NO2hotspots detected within the city of Chicago, IL, USA using three distinct high-resolution (1.3 km) air quality products: (1) an interpolated product from Microsoft Research’s Project Eclipse—a dense network of over 100 low-cost sensors; (2) a two-way coupled WRF-CMAQ simulation; and (3) a down-sampled product using TropOMI satellite instrument observations. We use the Getis-OrdGi*statistic to identify hotspots of NO2and stratify results into high-, medium-, and low-agreement hotspots, including one consensus hotspot detected in all three datasets. Interrogating medium- and low-agreement hotspots offers insights into dataset discrepancies, such as sensor placement and model physics considerations, data retrieval caveats, and the potential for missing emission inventories. When treated as complements rather than substitutes, our work demonstrates that novel air quality products can enable researchers to address discrepancies in data products and can help regulators evaluate confidence in policy-relevant insights.

 
more » « less
NSF-PAR ID:
10464253
Author(s) / Creator(s):
; ; ; ; ; ;
Publisher / Repository:
IOP Publishing
Date Published:
Journal Name:
Environmental Research Letters
Volume:
18
Issue:
10
ISSN:
1748-9326
Page Range / eLocation ID:
Article No. 104010
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract. The role of clouds in the Arctic radiation budget is not well understood. Ground-based and airborne measurements provide valuable data to test and improve our understanding. However, the ground-based measurements are intrinsically sparse, and the airborne observations are snapshots in time and space. Passive remote sensing measurements from satellite sensors offer high spatial coverage and an evolving time series, having lengths potentially of decades. However, detecting clouds by passive satellite remote sensing sensors is challenging over the Arctic because of the brightness of snow and ice in the ultraviolet and visible spectral regions and because of the small brightness temperature contrast to the surface. Consequently, the quality of the resulting cloud data products needs to be assessed quantitatively. In this study, we validate the cloud data products retrieved from the Advanced Very High Resolution Radiometer (AVHRR) post meridiem (PM) data from the polar-orbiting NOAA-19 satellite and compare them with those derived from the ground-based instruments during the sunlit months. The AVHRR cloud data products by the European Space Agency (ESA) Cloud Climate Change Initiative (Cloud_CCI) project uses the observations in the visible and IR bands to determine cloud properties. The ground-based measurements from four high-latitude sites have been selected for this investigation: Hyytiälä (61.84∘ N, 24.29∘ E), North Slope of Alaska (NSA; 71.32∘ N, 156.61∘ W), Ny-Ålesund (Ny-Å; 78.92∘ N, 11.93∘ E), and Summit (72.59∘ N, 38.42∘ W). The liquid water path (LWP) ground-based data are retrieved from microwave radiometers, while the cloud top height (CTH) has been determined from the integrated lidar–radar measurements. The quality of the satellite products, cloud mask and cloud optical depth (COD), has been assessed using data from NSA, whereas LWP and CTH have been investigated over Hyytiälä, NSA, Ny-Å, and Summit. The Cloud_CCI COD results for liquid water clouds are in better agreement with the NSA radiometer data than those for ice clouds. For liquid water clouds, the Cloud_CCI COD is underestimated roughly by 3 optical depth (OD) units. When ice clouds are included, the underestimation increases to about 5 OD units. The Cloud_CCI LWP is overestimated over Hyytiälä by ≈7 g m−2, over NSA by ≈16 g m−2, and over Ny-Å by ≈24 g m−2. Over Summit, CCI LWP is overestimated for values ≤20 g m−2 and underestimated for values >20 g m−2. Overall the results of the CCI LWP retrievals are within the ground-based instrument uncertainties. To understand the effects of multi-layer clouds on the CTH retrievals, the statistics are compared between the single-layer clouds and all types (single-layer + multi-layer). For CTH retrievals, the Cloud_CCI product overestimates the CTH for single-layer clouds. When the multi-layer clouds are included (i.e., all types), the observed CTH overestimation becomes an underestimation of about 360–420 m. The CTH results over Summit station showed the highest biases compared to the other three sites. To understand the scale-dependent differences between the satellite and ground-based data, the Bland–Altman method is applied. This method does not identify any scale-dependent differences for all the selected cloud parameters except for the retrievals over the Summit station. In summary, the Cloud_CCI cloud data products investigated agree reasonably well with those retrieved from ground-based measurements made at the four high-latitude sites.

     
    more » « less
  2. Abstract

    Gridded monthly rainfall estimates can be used for a number of research applications, including hydrologic modeling and weather forecasting. Automated interpolation algorithms, such as the “autoKrige” function in R, can produce gridded rainfall estimates that validate well but produce unrealistic spatial patterns. In this work, an optimized geostatistical kriging approach is used to interpolate relative rainfall anomalies, which are then combined with long-term means to develop the gridded estimates. The optimization consists of the following: 1) determining the most appropriate offset (constant) to use when log-transforming data; 2) eliminating poor quality data prior to interpolation; 3) detecting erroneous maps using a machine learning algorithm; and 4) selecting the most appropriate parameterization scheme for fitting the model used in the interpolation. Results of this effort include a 30-yr (1990–2019), high-resolution (250-m) gridded monthly rainfall time series for the state of Hawai‘i. Leave-one-out cross validation (LOOCV) is performed using an extensive network of 622 observation stations. LOOCV results are in good agreement with observations (R2= 0.78; MAE = 55 mm month−1; 1.4%); however, predictions can underestimate high rainfall observations (bias = 34 mm month−1; −1%) due to a well-known smoothing effect that occurs with kriging. This research highlights the fact that validation statistics should not be the sole source of error assessment and that default parameterizations for automated interpolation may need to be modified to produce realistic gridded rainfall surfaces. Data products can be accessed through the Hawai‘i Data Climate Portal (HCDP;http://www.hawaii.edu/climate-data-portal).

    Significance Statement

    A new method is developed to map rainfall in Hawai‘i using an optimized geostatistical kriging approach. A machine learning technique is used to detect erroneous rainfall maps and several conditions are implemented to select the optimal parameterization scheme for fitting the model used in the kriging interpolation. A key finding is that optimization of the interpolation approach is necessary because maps may validate well but have unrealistic spatial patterns. This approach demonstrates how, with a moderate amount of data, a low-level machine learning algorithm can be trained to evaluate and classify an unrealistic map output.

     
    more » « less
  3. Abstract

    The13C/12C of dissolved inorganic carbon (δ13CDIC) carries valuable information on ocean biological C‐cycling, air‐sea CO2exchange, and circulation. Paleo‐reconstructions of oceanic13C from sediment cores provide key insights into past as changes in these three drivers. As a step toward full inclusion of13C in the next generation of Earth system models, we implemented13C‐cycling in a 1° lateral resolution ocean‐ice‐biogeochemistry Geophysical Fluid Dynamics Laboratory (GFDL) model driven by Common Ocean Reference Experiment perpetual year forcing. The model improved the mean of modernδ13CDICover coarser resolution GFDL‐model implementations, capturing the Southern Ocean decline in surfaceδ13CDICthat propagates to the deep sea via deep water formation. Controls onδ13CDICof the deep‐sea are quantified using both observations and model output. The biological control is estimated from the relationship between deep‐sea Pacificδ13CDICand phosphate (PO4). Theδ13CDIC:PO4slope from observations is revised to a value of 1.01 ± 0.02‰ (μmol kg−1)−1, consistent with a carbon to phosphate ratio of organic matter (C:Porg) of 124 ± 10. Model output yields a lowerδ13CDIC:PO4than observed due to too low C:Porg. The ocean circulation impacts deep modernδ13CDICin two ways, via the relative proportion of Southern Ocean and North Atlantic deep water masses, and via the preindustrialδ13CDICof these water mass endmembers. Theδ13CDICof the endmembers ventilating the deep sea are shown to be highly sensitive to the wind speed dependence of air‐sea CO2gas exchange. Reducing the coefficient for air‐sea gas exchange following OMIP‐CMIP6 protocols improves significantly surfaceδ13CDICrelative to previous gas exchange parameterizations.

     
    more » « less
  4. Abstract

    We present the analysis of ∼100 pc scale compact radio continuum sources detected in 63 local (ultra)luminous infrared galaxies (U/LIRGs;LIR≥ 1011L), using FWHM ≲ 0.″1–0.″2 resolution 15 and 33 GHz observations with the Karl G. Jansky Very Large Array. We identify a total of 133 compact radio sources with effective radii of 8–170 pc, which are classified into four main categories—“AGN” (active galactic nuclei), “AGN/SBnuc” (AGN-starburst composite nucleus), “SBnuc” (starburst nucleus), and “SF” (star-forming clumps)—based on ancillary data sets and the literature. We find that “AGN” and “AGN/SBnuc” more frequently occur in late-stage mergers and have up to 3 dex higher 33 GHz luminosities and surface densities compared with “SBnuc” and “SF,” which may be attributed to extreme nuclear starburst and/or AGN activity in the former. Star formation rates (SFRs) and surface densities (ΣSFR) are measured for “SF” and “SBnuc” using both the total 33 GHz continuum emission (SFR ∼ 0.14–13Myr−1, ΣSFR∼ 13–1600Myr−1kpc−2) and the thermal free–free emission from Hiiregions (median SFRth∼ 0.4Myr−1,ΣSFRth44Myr−1kpc−2). These values are 1–2 dex higher than those measured for similar-sized clumps in nearby normal (non-U/LIRGs). The latter also have a much flatter median 15–33 GHz spectral index (∼−0.08) compared with “SBnuc” and “SF” (∼−0.46), which may reflect higher nonthermal contribution from supernovae and/or interstellar medium densities in local U/LIRGs that directly result from and/or lead to their extreme star-forming activities on 100 pc scales.

     
    more » « less
  5. Abstract

    The partial pressure of carbon dioxide (pCO2) was surveyed across the Amazon River plume and the surrounding western tropical North Atlantic Ocean (15–0°N, 43–60°W) during three oceanic expeditions (May–June 2010, September–October 2011, and July 2012). The survey timing was chosen according to previously described temporal variability in plume behavior due to changing river discharge and winds.In situsea surfacepCO2and air‐sea CO2flux exhibited robust linear relationships with sea surface salinity (SSS; 15 < SSS < 35), although the relationships differed among the surveys. Regional distributions ofpCO2and CO2flux were estimated using SSS maps from high‐resolution ocean color satellite‐derived (MODIS‐Aqua) diffuse attenuation coefficient at 490 nm (Kd490) during the periods of study. Results confirmed that the plume is a net CO2sink with distinctive temporal variability: the strongest drawdown occurred during the spring flood (−2.39 ± 1.29 mmol m−2 d−1in June 2010), while moderate drawdown with relatively greater spatial variability was observed during the transitional stages of declining river discharge (−0.42 ± 0.76 mmol m−2 d−1in September–October 2011). The region turned into a weak source in July 2012 (0.26 ± 0.62 mmol m−2 d−1) when strong CO2uptake in the mid‐plume was overwhelmed by weak CO2outgassing over a larger area in the outer plume. Outgassing near the mouth of the river was observed in July 2012. Our observations draw attention to the importance of assessing the variable impacts of biological activity, export, and air‐sea gas exchange before estimating regional CO2fluxes from salinity distributions alone.

     
    more » « less