Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “MODIS”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,783 records · Page 99

Improvement of aerosol optical depth data for localized solar resource assessment

Solar irradiance, especially for direct normal irradiance (DNI), is sensitive to the atmospheric aerosols, which can strongly extinguish DNI over areas with elevated aerosol loadings. The National Solar Radiation Database (NSRDB), developed by the National Renewable Energy Laboratory (NREL), uses AOD from MERRA-2 reanalysis and is further downscaled from 0.5° to 2-km based on elevation. However, the elevation-based downscaling may not accurately represent AOD, especially over the areas with large AOD gradients. This study examined whether the 1-km MODIS MAIAC satellite-retrieved AOD product can better represent the AOD distribution and be used to improve solar irradiance assessment. We focused on areas with relatively high AOD over North America (California, and New York City, US, and Mexico City in Mexico, in particular). The evaluation of MAIAC AOD and MERRA-2 AOD against ground-truth AERONET AOD shows that MAIAC AOD exhibits smaller RMSE by about 0.05 and higher correlation coefficient by 0.1-0.6 than MERRA-2 AOD. In addition, MERRA-2 AOD shows larger negative bias. The simulated DNI using MAIAC AOD (DNIMAIAC) and using MERRA-2 AOD (DNINSRDB) were evaluated with DNI observations. The results show the performance of DNIMAIAC is better than that of DNINSRDB with smaller RMSE by 0.5-1.5%, smaller positive mean bias by 0.8-3.1% and comparable correlation for the 3 sites analyzed. Overall, 1-km MAIAC AOD shows higher accuracy than 2-km elevation-based MERRA-2 AOD, leading to better performance of the simulated DNI using MAIAC AOD. Therefore, 1-km MAIAC AOD can be used to improve the accuracy of solar resource assessment.

14 SOLAR ENERGY↗

Radio-echo sounding of Crary Ice Rise reveals abundant marine ice in former ice shelf rifts and basal crevasses

Crary Ice Rise formed after the Ross Ice Shelf re-grounded ~1 kyr BP. We present new ice-penetrating radar data from two systems operating at center frequencies of 7 and 750 MHz that confirm the ice rise is composed of a former ice shelf buried by subsequent accumulation. Stacks of englacial diffraction hyperbolas are present almost everywhere across the central ice rise and extend up to ~350 m above the bed. In many cases, bed reflections beneath the diffraction hyperbolas are obscured for distances up to 1 km. Waveform modeling indicates that the diffraction hyperbolas are likely caused by marine ice deposits in former basal crevasses and rifts. The in-filling of rifts and basal crevasses may have strengthened the connection between the ice rise and the surrounding ice shelf, which could have influenced local and regional ice dynamics. Three internal reflection horizons mark the upper limit of disturbed ice and diffraction hyperbolas in different sections of the ice rise, indicating at least three stages of flow stabilization across the ice rise. A surface lineation visible in MODIS imagery corresponds spatially to deepening and strong deformation of these layers, consistent with the characteristics of former grounding lines observed elsewhere in Antarctica.

54 ENVIRONMENTAL SCIENCES↗

Seasonal Variability in Local Carbon Dioxide Biomass Burning Sources Over Central and Eastern US Using Airborne In Situ Enhancement Ratios

We present observations of local enhancements in carbon dioxide (CO 2 ) from local emissions sources over three eastern US regions during four deployments of the Atmospheric Carbon Transport-America (ACT-America) campaign between summer 2016 and spring 2018. Local CO 2 emissions were characterized by carbon monoxide (CO) to CO 2 enhancement ratios (i.e., ΔCO/ΔCO 2 ) in air mass mixing observed during aircraft transects within the planetary boundary layer. By analyzing regional-scale variability of CO 2 enhancements as a function of ΔCO/ΔCO 2 enhancement ratios, observed relative contributions to CO 2 emissions were separated into fossil fuel and biomass burning (BB) regimes across regions and seasons. CO 2 emission contributions attributed to biomass burning (ΔCO/ΔCO 2 > 4%) were negligible during summer and fall in all regions but climbed to ∼9%–11% of observed combustion contributions in the South during winter and spring. Relative CO 2 fire emission trends matched observed winter and spring BB contributions, but conflictingly predicted similar levels of BB during the fall. Satellite fire data from MODIS and VIIRS suggested the use of higher spatial resolution fire data that might improve modeled BB emissions but were not able to explain the bulk of the discrepancy.

54 ENVIRONMENTAL SCIENCES↗

Liquid Containing Clouds at the North Slope of Alaska Demonstrate Sensitivity to Local Industrial Aerosol Emissions

Abstract Cloud condensation nucleus control alter cloud solar albedo through cloud droplet size. Here, we leverage anthropogenic emissions at the North Slope of Alaska as a natural laboratory to study relationships between aerosols and Arctic liquid‐containing clouds. Averaging 14 years of MODIS observations, we found a reduction in temporally averaged cloud effective radius ( ) of up to 1.0 μm related to localized pollution. Pronounced regional gradients in cloud frequency of occurrence and liquid water path prohibit the detection of potential changes of these variables. Observed changes of alter radiative fluxes and increase cloud‐reflected shortwave radiation by up to 0.8 W m − 2 in the Prudhoe Bay area for the period covered by observations (April–September). Due to the frequent occurrence of liquid‐containing clouds, this implies that enhanced local emissions in Arctic regions can impact climate processes.

54 ENVIRONMENTAL SCIENCES↗

Cloud Morphology Evolution in Arctic Cold‐Air Outbreak: Two Cases During COMBLE Period

Abstract Cloud feedbacks play an important role in Arctic warming. Cloud morphology, for example, cloud size and spatial distributions, is among key factors that directly impact their radiative effects. In this work, we use two cases observed during the Cold‐air Outbreak (CAO) in the Marine Boundary Layer Experiment (COMBLE) to study the evolution of cloud size distributions as an air mass is advected from the Arctic over a comparatively warm ocean and cloud mesoscale organization changes from rolls to cells. Cloud objects are identified from Moderate Resolution Imaging Spectroradiometer (MODIS) reflectance images through an object segmentation procedure and roll breakup is identified by homogeneities in cloud water path (CWP). Roll breakup is found to be accompanied by a local minimum in wind shear and local maxima in cloud size and marine cold air outbreak index. The mean cloud horizontal aspect ratio has weak fetch dependency and is around 2 in roll, transition, and cell regimes. Regardless of distance from the ice edge, smaller clouds (<10 km 2 ) dominate the population number but not cloud cover. Cloud size distributions show bimodality in transition and cell regimes. For clouds with comparable sizes, mean nearest neighbor distances normalized by equivalent cloud radius converge to a single value for all regimes and for all but the smallest clouds, suggesting that clouds of comparable sizes in CAOs are separated by distance proportional to their sizes. The presented statistical results pave the way to evaluating model simulated cloud organizations during CAO events.

54 ENVIRONMENTAL SCIENCES↗

The Global Atmosphere‐aerosol Model ICON‐A‐HAM2.3–Initial Model Evaluation and Effects of Radiation Balance Tuning on Aerosol Optical Thickness

Abstract The Hamburg Aerosol Module version 2.3 (HAM2.3) from the ECHAM6.3‐HAM2.3 global atmosphere‐aerosol model is coupled to the recently developed icosahedral nonhydrostatic ICON‐A (icon‐aes‐1.3.00) global atmosphere model to yield the new ICON‐A‐HAM2.3 atmosphere‐aerosol model. The ICON‐A and ECHAM6.3 host models use different dynamical cores, parameterizations of vertical mixing due to sub‐grid scale turbulence, and parameter settings for radiation balance tuning. Here, we study the role of the different host models for simulated aerosol optical thickness (AOT) and evaluate impacts of using HAM2.3 and the ECHAM6‐HAM2.3 two‐moment cloud microphysics scheme on several meteorological variables. Sensitivity runs show that a positive AOT bias over the subtropical oceans is remedied in ICON‐A‐HAM2.3 because of a different default setting of a parameter in the moist convection parameterization of the host models. The global mean AOT is biased low compared to MODIS satellite instrument retrievals in ICON‐A‐HAM2.3 and ECHAM6.3‐HAM2.3, but the bias is larger in ICON‐A‐HAM2.3 because negative AOT biases over the Amazon, the African rain forest, and the northern Indian Ocean are no longer compensated by high biases over the sub‐tropical oceans. ICON‐A‐HAM2.3 shows a moderate improvement with respect to AOT observations at AERONET sites. A multivariable bias score combining biases of several meteorological variables into a single number is larger in ICON‐A‐HAM2.3 compared to standard ICON‐A and standard ECHAM6.3. In the tropics, this multivariable bias is of similar magnitude in ICON‐A‐HAM2.3 and in ECHAM6.3‐HAM2.3. In the extra‐tropics, a smaller multivariable bias is found for ICON‐A‐HAM2.3 than for ECHAM6.3‐HAM2.3.

54 ENVIRONMENTAL SCIENCES↗

Performance of Cloud 3D Solvers in Ice Cloud Shortwave Radiation Closure Over the Equatorial Western Pacific Ocean

Abstract For retrieving cloud optical properties from satellite images or computing these properties from climate model output, computationally efficient treatments of cloud horizontal inhomogeneity include the Monte Carlo Independent Column Approximation (McICA) and the Tripleclouds method. Computationally efficient treatment of cloud horizontal radiation exchanges includes the SPeedy Algorithm for Radiative TrAnsfer through CloUd Sides (SPARTACUS). As a test to derive properties from satellite images, we collocate Moderate Resolution Imaging Spectroradiometer (MODIS) cloud retrievals with near‐nadir Cloud and the Earth's Radiant Energy System (CERES) footprints in July 2008 over an equatorial western Pacific Ocean region to compare the performance of the McICA, Tripleclouds, and SPARTACUS solvers to the conventional plane‐parallel homogeneous (PPH) treatment. PPH overestimates cloud albedo, and the three solvers effectively reduce overestimation with root mean square error of shortwave upwelling irradiance decreasing between 15.72 and 18.53 W m −2 , or about 22%–25%. Although cloud top variability does not get fed into the simulations, all three solvers also reduce the effect of cloud top variability on cloud albedo. Entrapment (energy reflected downward from clouds) and horizontal radiation transfer have opposite effects on the SPARTACUS cloud albedo simulation. The net effect depends on the cloud vertical extent, the unawareness of which limits the performance of the SPARTACUS solver.

54 ENVIRONMENTAL SCIENCES↗

The Association Between Cloud Droplet Number over the Summer Southern Ocean and Air Mass History

The cloud properties and governing processes in Southern Ocean marine boundary layer clouds have emerged as a central issue in understanding the Earth's climate sensitivity. While our understanding of Southern Ocean cloud feedbacks have evolved in the most recent climate model intercomparison, the background properties of simulated summertime clouds in the Southern Ocean are not consistent with measurements due to known biases in simulating cloud condensation nuclei concentrations. This paper presents several case studies collected during the Capricorn 2 and Marcus campaigns held aboard Australian research vessels in the Austral Summer of 2018. Combining the surface–observed cases with MODIS data along forward and backward air mass trajectories, we demonstrate the evolution of cloud properties with time. These cases are consistent with multi–year statistics showing that long trajectories of air masses over the Antarctic ice sheet are critical to creating high droplet number clouds in the high latitude summer Southern Ocean. We speculate that secondary aerosol production via the oxidation of biogenically derived aerosol precursor gasses over the high actinic flux region of the high latitude ice sheets is fundamental to maintaining relatively high droplet numbers in Southern Ocean clouds during Summer.

54 ENVIRONMENTAL SCIENCES↗

Surface Energy Dynamics and Canopy Structural Properties in Intact and Disturbed Forests in the Southern Amazon

The Brazilian Amazon has been a focus of land development with large swaths of forests converted to agriculture. Forest degradation by selective logging and fires has accompanied the agricultural frontier and has resulted in significant impacts on Amazonian ecosystems. Changes in forest structure resulting from forest disturbances have large impacts on the surface energy balance, including on land surface temperature (LST) and evapotranspiration (ET). This study's objective is to assess the effects of forest disturbances on water fluxes and forest structure in a transitional forest site in the Southern Amazon. We used ET and LST products from MODIS and Landsat 8 and GEDI-derived forest structure data to address our research questions. We found that disturbances induced seasonal water stress, more pronounced in croplands/pastures than in forests (differences up to 20% in the dry season), and more pronounced in second-growth and recently burned areas than in logged and intact forests (differences up to 12% in the dry season). Moreover, ET and LST were negatively related, with more consistent relationships across disturbance classes in the dry season (R 2 : 0.41–0.87) than in the wet season (R 2 : 0.18–0.49). Forest and cropland or pasture classes showed contrasting relationships in the dry season. Finally, we found that forest structure exhibited stronger relationships with ET and LST in the most disturbed forests (R 2 : 0.01–0.43) than in the least disturbed forests (R 2 < 0.05). Our findings help to elucidate degraded forests functioning under a changing climate and to improve estimates of water and energy fluxes in Amazonian degraded forests.

54 ENVIRONMENTAL SCIENCES↗

Harnessing Satellite Data Alone for Mapping Global Thermal Anisotropy

Mapping thermal anisotropy across global lands is critical for advancing a wide range of Earth science studies. However, a comprehensive understanding of global thermal anisotropy intensity (TAI) and its governing factors remains missing. We introduce a novel data-driven methodology to quantify global TAI exclusively using multi-angle MODIS land surface temperature time series observations. Our analysis reveals distinct seasonal and diurnal TAI patterns, with global mean summertime TAI exceeding 2.9°C. Furthermore, we identify strong associations between TAI and key surface and atmospheric parameters, such as leaf area index and downward shortwave radiation. Our findings advocate for a paradigm shift from model-based to data-driven approaches in correcting thermal anisotropy, thereby addressing a critical bottleneck in Earth observation.

54 ENVIRONMENTAL SCIENCES↗

Integrating State Data Assimilation and Innovative Model Parameterization Reduces Simulated Carbon Uptake in the Arctic and Boreal Region

Model representation of carbon uptake and storage is essential for accurate projection of the response of the arctic-boreal zone to a rapidly changing climate. Land model estimates of LAI and aboveground biomass that can have a marked influence on model projections of carbon uptake and storage vary substantially in the arctic and boreal zone, making it challenging to correctly evaluate model estimates of Gross Primary Productivity (GPP). To understand and correct bias of LAI and aboveground biomass in the Community Land Model (CLM), we assimilated the 8-day Moderate Resolution Imaging Spectroradiometer (MODIS) LAI observation and a machine learning product of annual aboveground biomass into CLM using an Ensemble Adjustment Kalman Filter (EAKF) in an experimental region including Alaska and Western Canada. Assimilating LAI and aboveground biomass reduced these model estimates by 58% and 72%, respectively. The change of aboveground biomass was consistent with independent estimates of canopy top height at both regional and site levels. The International Land Model Benchmarking system assessment showed that data assimilation significantly improved CLM's performance in simulating the carbon and hydrological cycles, as well as in representing the functional relationships between LAI and other variables. Here, to further reduce the remaining bias in GPP after LAI bias correction, we re-parameterized CLM to account for low temperature suppression of photosynthesis. The LAI bias corrected model that included the new parameterization showed the best agreement with model benchmarks. Combining data assimilation with model parameterization provides a useful framework to assess photosynthetic processes in LSMs.

58 GEOSCIENCES↗

Remotely Sensed High‐Resolution Soil Moisture and Evapotranspiration: Bridging the Gap Between Science and Society

This paper reviews the current state of high‐resolution remotely sensed soil moisture (SM) and evapotranspiration (ET) products and modeling, and the coupling relationship between SM and ET. SM downscaling approaches for satellite passive microwave products leverage advances in artificial intelligence and high‐resolution remote sensing using visible, near‐infrared, thermal‐infrared, and synthetic aperture radar sensors. Remotely sensed ET continues to advance in spatiotemporal resolutions from MODIS to ECOSTRESS to Hydrosat and beyond. These advances enable a new understanding of bio‐geo‐physical controls and coupled feedback mechanisms between SM and ET reflecting the land cover and land use at field scale (3–30 m, daily). Still, the state‐of‐the‐science products have their challenges and limitations, which we detail across data, retrieval algorithms, and applications. We describe the roles of these data in advancing 10 application areas: drought assessment, food security, precision agriculture, soil salinization, wildfire modeling, dust monitoring, flood forecasting, urban water, energy, and ecosystem management, ecohydrology, and biodiversity conservation. We discuss that future scientific advancement should focus on developing open‐access, high‐resolution (3–30 m), sub‐daily SM and ET products, enabling the evaluation of hydrological processes at finer scales and revolutionizing the societal applications in data‐limited regions of the world, especially the Global South for socio‐economic development.

54 ENVIRONMENTAL SCIENCES↗

Representing Fine‐Scale Topographic Effects on Surface Radiation Balance in Hyper‐Resolution Land Surface Models

Land surface models are increasingly used to simulate land surface processes at hyper-spatial resolutions (e.g., ∼1 km). As model resolution increases, grid-scale topographic effects on surface radiation fluxes and their interactions between adjacent grids become more pronounced. However, current land surface models routinely neglect the fine-scale topographic effects on surface radiation balance. This study developed physically-based and computationally-efficient parameterizations (fineTOP) that explicitly resolve fine-scale topographic effects on downward shortwave and longwave radiation as well as land surface radiative properties. The newly developed parameterizations were implemented and tested in the Energy Exascale Earth System Model (E3SM) Land Model (ELM). Multi-decadal km-resolution ELM simulations over the California Sierra Nevada show that fine-scale topography significantly impacts the surface energy balance and snow processes across seasons. Slope determines the magnitude of topographic effects, while aspect controls their sign. For slopes larger than 30°, topography-induced change in annual surface temperature can be as large as 3.3 K. Regionally, the mean value and standard deviation of topography-induced changes in annual surface temperature are −0.22 ± 0.38 K and +0.25 ± 0.37 K over north-facing and south-facing slopes, respectively. Topography-induced changes in surface radiative properties account for 3.5% ± 13.8% of total topographic effects on annual net radiation. With fineTOP, ELM captures the aspect-dependence of snow cover fraction, snow water equivalent, and land surface temperature found in MODIS satellite observations and a snow reanalysis data set, while the default ELM fails to capture this phenomenon. The enhanced capability to represent fine-scale topographic effects on surface radiation balance can be used to advance understanding of the role of fine-scale topography in land surface processes and land-atmosphere interactions over mountainous regions.

Hao, Dalei [Pacific Northwest National Laboratory ↗

Seasonal stabilization effects slowed the greening of the Northern Hemisphere over the last two decades

Rising atmospheric CO₂ and warming spring temperatures increase vegetation growth and the terrestrial carbon sink. However, drought, heat stress, phenology, and resource limitations may stabilize or limit theses projected increases. We investigate the balance between these amplifying and stabilizing ecological factors by asking whether enhanced early-season growth leads to continued late-season growth. Using the Moderate Resolution Imaging Spectroradiometer (MODIS) leaf area index (LAI) dataset, we identify three seasonal growth patterns based on early- and peak-season positive LAI anomalies: (1) amplification, where late-season LAI anomalies exceed earlier ones; (2) weak stabilization, where late-season anomalies remain similar or slightly lower; and (3) strong stabilization, where late-season anomalies become negative. Weak and strong stabilization events dominate across 67% and 26% of Northern Hemisphere ecosystems above 30°N, respectively. The absence of any trend in amplifying or stabilizing events suggests stabilizing factors seasonally offset CO₂ and temperature-induced spring greening. Terrestrial biosphere models underestimate strong stabilization and overestimate amplification events. This inconsistency arises from the models’ underestimation late-season LAI sensitivity to precipitation in water-limited regions; overlook negative legacy effects of early enhanced LAI on late-season soil moisture via evapotranspiration losses in energy-limited regions. Our findings suggest water/heat stress and resource limitations limit greening and the land carbon sink.

Zhang, Wen↗

A biomass map of the Brazilian Amazon from multisource remote sensing

The Amazon Forest, the largest contiguous tropical forest in the world, stores a significant fraction of the carbon on land. Changes in climate and land use affect total carbon stocks, making it critical to continuously update and revise the best estimates for the region, particularly considering changes in forest dynamics. Forest inventory data cover only a tiny fraction of the Amazon region, and the coverage is not sufficient to ensure reliable data interpolation and validation. This paper presents a new forest above-ground biomass map for the Brazilian Amazon and the associated uncertainty both with a resolution of 250 meters and baseline for the satellite dataset the year of 2016 (i.e., the year of the satellite observation). A significant increase in data availability from forest inventories and remote sensing has enabled progress towards high-resolution biomass estimates. This work uses the largest airborne LiDAR database ever collected in the Amazon, mapping 360,000 km 2 through transects distributed in all vegetation categories in the region. The map uses airborne laser scanning (ALS) data calibrated by field forest inventories that are extrapolated to the region using a machine learning approach with inputs from Synthetic Aperture Radar (PALSAR), vegetation indices obtained from the Moderate-Resolution Imaging Spectroradiometer (MODIS) satellite, and precipitation information from the Tropical Rainfall Measuring Mission (TRMM). A total of 174 field inventories geolocated using a Differential Global Positioning System (DGPS) were used to validate the biomass estimations. The experimental design allowed for a comprehensive representation of several vegetation types, producing an above-ground biomass map varying from a maximum value of 518 Mg ha -1 , a mean of 174 Mg ha -1 , and a standard deviation of 102 Mg ha -1 . This unique dataset enabled a better representation of the regional distribution of the forest biomass and structure, providing further studies and critical information for decision-making concerning forest conservation, planning, carbon emissions estimate, and mechanisms for supporting carbon emissions reductions.

54 ENVIRONMENTAL SCIENCES↗

CubeSat constellations provide enhanced crop phenology and digital agricultural insights using daily leaf area index retrievals

Satellite remote sensing has great potential to deliver on the promise of a data-driven agricultural revolution, with emerging space-based platforms providing spatiotemporal insights into precision-level attributes such as crop water use, vegetation health and condition and crop response to management practices. Using a harmonized collection of high-resolution Planet CubeSat, Sentinel-2, Landsat-8 and additional coarser resolution imagery from MODIS and VIIRS, we exploit a multi-satellite data fusion and machine learning approach to deliver a radiometrically calibrated and gap-filled time-series of daily leaf area index (LAI) at an unprecedented spatial resolution of 3 m. The insights available from such high-resolution CubeSat-based LAI data are demonstrated through tracking the growth cycle of a maize crop and identifying observable within-field spatial and temporal variations across key phenological stages. Daily LAI retrievals peaked at the tasseling stage, demonstrating their value for fertilizer and irrigation scheduling. An evaluation of satellite-based retrievals against field-measured LAI data collected from both rain-fed and irrigated fields shows high correlation and captures the spatiotemporal development of intra- and inter-field variations. Novel agricultural insights related to individual vegetative and reproductive growth stages were obtained, showcasing the capacity for new high-resolution CubeSat platforms to deliver actionable intelligence for precision agricultural and related applications.

54 ENVIRONMENTAL SCIENCES↗

A Multi-Analysis Approach for Estimating Regional Health Impacts from the 2017 Northern California Wildfires

In the evening of October 8 and early hours of October 9, 2017, high winds in Northern California downed trees and power lines, igniting some of the most devastating wildfires the state had seen, and within hours unhealthy air quality impacted millions of people. We simulated these air quality conditions using fire detection information from the MODIS, VIIRS, and GOES-16 ABI satellite instruments, and applying a set of three WRF–CMAQ simulations, one data fusion, and three machine learning methods. We investigated using the 5-min available GOES-16 fire detection data to simulate timing of fire activity to allocate emissions hourly for the WRF-CMAQ air quality modeling system. Interestingly, this approach did not necessarily improve results compared to the baseline case, which used a default time profile. However, this approach was key to simulating the initial 12-hr explosive fire activity and smoke impacts. The WRF-CMAQ simulations compared well with observational data for the October 8-15 time period and tended to overestimate concentrations October 16-20. To improve these results, we applied three machine learning algorithms. We also had a unique opportunity to evaluate results with temporary monitors deployed specifically for wildfires, and performance was markedly different. For example, at the permanent monitoring locations, the WRF-CMAQ simulations had a Pearson correlation of 0.65, and the data fusion approach improved this (Pearson correlation = 0.95), while at the temporary monitor locations across the WRF-CMAQ, data fusion, and machine learning datasets, the best Pearson correlation was 0.5. The data fusion and machine learning results were biased low and WRF-CMAQ results were biased high. Finally, we applied the optimized PM2.5 exposure estimate in a short-term exposure-response function. Total estimated mortality attributable to PM2.5 exposure during the smoke episode was 83 (95% confidence interval: 0, 196) with 47% of these deaths attributable to wildland fire smoke.

O'Neill, Susan↗