Replacement of Dry Deposition Velocities (Vd) at CASTNET Sites (cont d) Gary Lear George Bowker CAMD

Similar documents
Climate Dataset: Aitik Closure Project. November 28 th & 29 th, 2018

Northern New England Climate: Past, Present, and Future. Basic Concepts

Monthly probabilistic drought forecasting using the ECMWF Ensemble system

Temporal validation Radan HUTH

Gapfilling of EC fluxes

Reliability of Daily and Annual Stochastic Rainfall Data Generated from Different Data Lengths and Data Characteristics

Urban Forest Effects-Dry Deposition (UFORE D) Model Enhancements. Satoshi Hirabayashi

Water information system advances American River basin. Roger Bales, Martha Conklin, Steve Glaser, Bob Rice & collaborators UC: SNRI & CITRIS

AERMOD Sensitivity to AERSURFACE Moisture Conditions and Temporal Resolution. Paper No Prepared By:

Drought Monitoring in Mainland Portugal

Extremes Events in Climate Change Projections Jana Sillmann

A new sub-daily rainfall dataset for the UK

A sensitivity and uncertainty analysis. Ministry of the Walloon Region Agricultural Research Centre

APPENDIX 3.6-A Support Information for Newcastle, Wyoming Meteorological Monitoring Site

MMIF-processed WRF data for AERMOD Case Study: North ID mountain terrain

Ammonia Emissions and Nitrogen Deposition in the United States and China

DROUGHT ASSESSMENT USING SATELLITE DERIVED METEOROLOGICAL PARAMETERS AND NDVI IN POTOHAR REGION

Urban heat island in the metropolitan area of São Paulo and the influence of warm and dry air masses during summer

The Canadian ADAGIO Project for Mapping Total Atmospheric Deposition

Laboratory Exercise #7 - Introduction to Atmospheric Science: The Seasons and Daily Weather

A re-sampling based weather generator

The Importance of Ammonia in Modeling Atmospheric Transport and Deposition of Air Pollution. Organization of Talk:

Drought Criteria. Richard J. Heggen Department of Civil Engineering University of New Mexico, USA Abstract

AN INTERNATIONAL SOLAR IRRADIANCE DATA INGEST SYSTEM FOR FORECASTING SOLAR POWER AND AGRICULTURAL CROP YIELDS

Climate Change and Arizona s Rangelands: Management Challenges and Opportunities

The Latest Science of Seasonal Climate Forecasting

UPPLEMENT A COMPARISON OF THE EARLY TWENTY-FIRST CENTURY DROUGHT IN THE UNITED STATES TO THE 1930S AND 1950S DROUGHT EPISODES

FOREST FIRE HAZARD MODEL DEFINITION FOR LOCAL LAND USE (TUSCANY REGION)

Use of reanalysis data in atmospheric electricity studies

Operational MRCC Tools Useful and Usable by the National Weather Service

A re-sampling based weather generator

Indices and Indicators for Drought Early Warning

Graphical Techniques Stem and Leaf Box plot Histograms Cumulative Frequency Distributions

METADATA MANAGEMENT TOOL (LIO)

Variability of Reference Evapotranspiration Across Nebraska

HEC-HMS Lab 4 Using Frequency Storms in HEC-HMS

Atmospheric Mercury Deposition Modeling

Application and verification of ECMWF products 2017

Using Innovative Displays of Hydrologic Ensemble Traces

Ambient Ammonium Contribution to total Nitrogen Deposition

BETWIXT Built EnvironmenT: Weather scenarios for investigation of Impacts and extremes. BETWIXT Technical Briefing Note 7 Version 1, May 2006

YACT (Yet Another Climate Tool)? The SPI Explorer

Climate Variables for Energy: WP2

Table 13 Maximum dewpoint temperature (degrees Celsius), primary meteorological station, H. J. Andrews Experimental Forest, 5/10/72 through 12/31/84

Francina Dominguez*, Praveen Kumar Department of Civil and Environmental Engineering University of Illinois at Urbana-Champaign

WEATHER NORMALIZATION METHODS AND ISSUES. Stuart McMenamin Mark Quan David Simons

A STUDY ON EVAPORATION IN IOANNINA, NW GREECE

UPPLEMENT. PRECIPITATION VARIABILITY AND EXTREMES IN CENTRAL EUROPE New View from STAMMEX Results

Snow and glacier change modelling in the French Alps

Comparing the Relationships Between Heat Stress Indices and Mortality in North Carolina

Global Climate Change and the Implications for Oklahoma. Gary McManus Associate State Climatologist Oklahoma Climatological Survey

Dust Storm: An Extreme Climate Event in China

Daily Rainfall Disaggregation Using HYETOS Model for Peninsular Malaysia

DEVELOPMENT OF A LARGE-SCALE HYDROLOGIC PREDICTION SYSTEM

PRELIMINARY DRAFT FOR DISCUSSION PURPOSES

peak half-hourly Tasmania

5B.1 DEVELOPING A REFERENCE CROP EVAPOTRANSPIRATION CLIMATOLOGY FOR THE SOUTHEASTERN UNITED STATES USING THE FAO PENMAN-MONTEITH ESTIMATION TECHNIQUE

peak half-hourly New South Wales

THE METEOROLOGICAL DATA QUALITY MANAGEMENT OF THE ROMANIAN NATIONAL SURFACE OBSERVATION NETWORK

An index to indicate precipitation probability and to investigate effects of sub-grid-scale surface parameterizations on model performance

Approach Estimating Mercury Dry Deposition for AMNeT Leiming Zhang

3/3/2013. The hydro cycle water returns from the sea. All "toilet to tap." Introduction to Environmental Geology, 5e

Reduced Overdispersion in Stochastic Weather Generators for Statistical Downscaling of Seasonal Forecasts and Climate Change Scenarios

Memo. I. Executive Summary. II. ALERT Data Source. III. General System-Wide Reporting Summary. Date: January 26, 2009 To: From: Subject:

Heavier summer downpours with climate change revealed by weather forecast resolution model

Item 1: Edits to Select Appendix S Class I Air Dispersion Modeling Report text

Time Series Analysis -- An Introduction -- AMS 586

Validation of the Weather Generator CLIGEN with Precipitation Data from Uganda. W. J. Elliot C. D. Arnold 1

Exploration of Nitrogen Total Deposition Budget Uncertainty at the Regional Scale

Trevor Lee Director, Buildings. Grant Edwards PhD Department of Environment and Geography

Development of an ensemble-based volcanic ash dispersion model for operations at Darwin VAAC

Impact on Agriculture

Impact of aerosol on air temperature in Baghdad

PROPOSAL OF SEVEN-DAY DESIGN WEATHER DATA FOR HVAC PEAK LOAD CALCULATION

PREDICTING OVERHEATING RISK IN HOMES

Analysis of gross alpha, gross beta activities and beryllium-7 concentrations in surface air: their variation and statistical prediction model

Solar Radiation in Thailand:

SWIM and Horizon 2020 Support Mechanism

Improved Fields of Satellite-Derived Ocean Surface Turbulent Fluxes of Energy and Moisture

COURSE CLIMATE SCIENCE A SHORT COURSE AT THE ROYAL INSTITUTION

Who is polluting the Columbia River Gorge?

A quantile matching adjustment algorithm for Gaussian data series

Model Output Statistics (MOS)

Characteristics of long-duration precipitation events across the United States

Investigation of the January Thaw in the Blue Hill Temperature Record

Course Outline CLIMATE SCIENCE A SHORT COURSE AT THE ROYAL INSTITUTION. 1. Current climate. 2. Changing climate. 3. Future climate change

Appendix 1: UK climate projections

Chapter 2 Available Solar Radiation

Vermont Soil Climate Analysis Network (SCAN) sites at Lye Brook and Mount Mansfield

Ashraf S. Zakey The Egyptian Meteorological Autority

Weather Prediction Using Historical Data

Northern Rockies Adaptation Partnership: Climate Projections

Weather and climate outlooks for crop estimates

VARIABILITY OF SOLAR RADIATION OVER SHORT TIME INTERVALS

2/23/2015 GEOGRAPHY 204: STATISTICAL PROBLEM SOLVING IN GEOGRAPHY THE NORMAL DISTRIBUTION THE NORMAL DISTRIBUTION

Examining Changes to the Temporal Distribution of Seasonal-to-Annual Precipitation across the Mid-Atlantic Region,

Advancements and validation of the global CryoClim snow cover extent product

Table of Contents. Page

How to Run the Analysis: To run a principal components factor analysis, from the menus choose: Analyze Dimension Reduction Factor...

Laboratory Exercise #7 - Introduction to Atmospheric Science: The Seasons

Transcription:

Replacement of Dry Deposition Velocities (Vd) at CASTNET Sites (cont d) Gary Lear George Bowker CAMD

Introduction For CASTNET, hourly deposition velocity estimates are the output of the Multi-Layer Model (MLM) MLM utilizes numerous hourly meteorological parameters as inputs. Windspeed Sigma theta (turbulence) Temperature Precipitation Solar radiation Deposition Velocity (Vd) If any of the met parameters is missing, the hourly Vd estimate is also missing. Missing observations can skew aggregations and distort trends

Introduction CASTNET data completeness for total deposition is cumulative for all variables Missing data = $ EPA discontinued met measurements at our CASTNET sites so now it is essential to estimate deposition velocities without these measurements. Percent (%) 100 90 80 70 60 50 40 30 20 10 0 Concentrations Data Completeness Meteorological measurements Deposition velocity Total (wet + dry) deposition

Replacement of limited missing/invalid deposition velocity estimates. If longer-term (annual, seasonal) dry deposition estimates are most important for deposition assessments, it may be possible to impute (replace) the missing Vd s using historical values. We observed that replacing a missing/invalid Vd with a historical Vd could simplify the replacement process and could probably result in a moreaccurate long-term deposition estimate (since the historic variability of a longer-term average is just a few percent). This is based on the concept that we replace missing hours of Vd with a representative historical Vd from that same hour. Quality Assurance Decisions with Air Models: A Case Study of Imputation of Missing Input Data Using EPA?s Multi-layer Model. G.E. Bowker, D.B. Schwede, G.G. Lear, W.J. Warren-Hicks, P.L. Finkelstein. 2010. Water, air, and soil pollution. DOI: 10.1007/s11270-011-0808-7DOI: 10.1007/s11270-011-0808-7

Seasonal Variation in Vd Although meteorology dominates the magnitude of Vd at short timescales such as hours or days, longer time scales are dominated by site-specific seasonal trends Weekly deposition velocity for individual years Average Vd for all years (solid black)

Seasonal Variation in Vd

For each site: Seasonal variation is relatively consistent from year to year Variance of each week is similar from year to year Diurnal patterns are also very consistent with seasonal variations (not shown) Seasonal Variation in Vd

How do we quantify the error of these imputed data sets? Variable mixture of real measurements and imputed values Both real and imputed measurements have some but different! -- degree of autocorrelation Missing periods are not of truly random occurrence or duration and are difficult to model Classical statistics (e.g., mean and variance) are tricky in these situations!

Our solution: Quantify uncertainty of imputation scheme using known data sets replaced with imputed values and compare known with estimate From the entire CASTNET data record: Identify site-years with essentially complete data ( Master data). Create a complete Historical average time series (average hourly Vd) based on at least 5-years. Historical average for site Master data is a complete year at this site Vd Time (season)

Create Mask patterns Create Mask patterns from patterns of missing data from different sites The black line shows pattern of missing data. It is a Mask file. Vd Solid blue line is a representative year at a different site with missing data. It is used to make a Mask file. Time (season)

Vd Vd Master-Mask files Replace Master data with historical averages where mask pattern indicates a missing value to make Master-Mask Calculate the average Vd for the replaced file. Compare to the original Master file Time (season) Time (season)

Rinse and repeat 132 Master site-year files. 152 Mask files Result Many thousands of different combinations, each with a different percent missing.

Difference of Imputed Values Compared to Known Relative percent differences (RPD) of the annual aggregate of the Master-Mask imputed values compared to the known Master data The same figure, but looking at the absolute relative percent difference (ARPD) Annual Averages The red squares are the standard deviation of the decile

How robust is the scheme? We want to compare results across parameters, seasons and sites. So, we need to put them in the same units: Normalize values by dividing them by the standard deviation at 100% missing

How robust is the scheme? Normalize the figure, by dividing all values by the standard deviation at 100% missing. The result is the Standard Score, or Z-score, which is the number of standard deviations from the mean. These values are approximately three standard deviations from the mean. The red line has a value of 1 standard deviation at 100% missing.

Compare by Parameter Annual Averages

Compare by Season Summer

Compare by Season Spring

Compare by Season Fall

Compare by Season Winter

Uncertainty vs. Percent Replaced By parameter, how does uncertainty vary with amount of data replaced? Plot the Z-scores as a function of the percent of data replaced

Uncertainty vs. Percent Replaced When normalized, the patterns for all seasons appear to be the same. Same with parameters!

Uncertainty vs. Percent Replaced Conclusion: The patterns across parameters are very similar, and indicates the data replacement scheme is quite robust. y = -4.4E-05x 2 + 0.0142x

How do we estimate uncertainty for specific sites, seasons and parameters? Data needed: Mean and standard deviation for the site based on the historic CASTNET database for the season and parameter Percentage of missing data for your site-year. The uncertainty equation: y = -4.4E-05x 2 + 0.0142x (where x is the percent of missing data, and y is the estimated uncertainty in units of standard deviations for your site s season and parameter). Result: The average expected uncertainty expressed as a standard deviation (percent of the average value).

How do the imputed annual fluxes compare with previous estimates?

Explore how these changes affect the data http://epa.gov/castnet/javaweb/mcannual_replaced.html