Spatio-temporal modeling of fine particulate matter

Size: px
Start display at page:

Download "Spatio-temporal modeling of fine particulate matter"

Transcription

1 Spatio-temporal modeling of fine particulate matter Sujit Sahu a University of Southampton Lisbon: March 2004 a Joint work with Alan Gelfand and David Holland

2 Fine particulate matter Particles with diameter less than 2.5 m, PM. PM has been linked to respiratory and cardiovascular problems. Contribute to other air-pollutants, stress to vegetation and ecosystems. In 1997, EPA gave new regulations and setup a network of monitoring sites.

3 PM formation PM can be emitted directly or formed in the atmosphere: road-dust, wood burning. It contains particles from gaseous emissions such as sulfates and nitrates. Power plant, automobile. PM concentration levels are affected by emissions, topography, land cover, and other factors, e.g. air movement.

4 Our data set Mid-western United states: Illinois, Indiana and Ohio. Data on 114 sites in the year Data are 24-hr averages. Sampling frequencies are varied: daily, every three days, weekly etc. We work with weekly averages. Have 5928 (114 times 52) data points. Only 2% missing observations.

5 114 sites in Illinois, Indiana, and Ohio cements

6 Objectives in modeling Investigate rural/urban differences Measure effect of population/density Produce weekly prediction surfaces Obtain annual predictions

7 Outline of the talk Exploratory analysis Methods Modeling details Prediction details Analysis Model choice Parameter estimates Prediction Discussion

8 Scale of Data Variance increases as mean increases. Model square root of data, but report predictions on original scale. Mean Variance Mean Variance Figure 1: Top: original; bottom: square root scale.

9 Rural/Urban separation cements PM Rural Suburban and urban 10 0 Site Figure 2: First 10 for rural; last 104 for suburban and urban.

10 Rural/Urban mean variance Rural Suburban Urban Will combine suburban and urban into one category.! " # $ if # rural, =1 otherwise. Problem in prediction: don t have these for a fine grid. Will develop a spatial model for! %. '&

11 Population density Have density for the year Can pop-density separate the means and vars? No! see below. mean variance low-pop high-pop

12 * 60 () PM cements Week cements PM Jan Feb Mar Apr May Jun Jul Aug Sep Oct Nov Dec Month

13 Temporal effects Higher in summer due to higher sulfate particle concentrations Higher in some winter months due to increased levels of nitrate concentrations, possible re-intrainment of road salt We shall use monthly indicators to model. '

14 Geodetic Distance Earth as a sphere of radius 6371 kilometers. Let,.- and, denote the latitudes and /0- and /1 denote the longitudes (in radian). d # :<;=;?>A@B!C (km) C I, - J@EDKFH I,4 LNM ;?>O@P I, - J;?>O@P I,4 QJ;?>O@P R/ - S /1 L.

15 Variograms Obtain variation in the residuals after fitting a linear model. There is spatial variation. Probably no sill. Points to non-stationary spatial variation. variogram All sites Urban Sites distance +

16 Hierarchical model T UWVYX[Z\] # PM ' ^ at location?v, week Z\. data = systematic process + white noise _!`VaX[Zb\] # c!`vdx[zb\] M ef!`vyx[zb\] systematic = mean + rural/background process + urban indicator * urban process c!wvyx[z\] #!`VdX[Z\P+M g!`vhx[z\]+m!`vijik!wvyx[z\]ql Mean process = pop-density + urban + urban * pop-density + month '

17 t v Spatio-temporal process Assume that g U?VYX[Zb\] and ik UWVaX[Z\P are zero mean spatio-temporal processes. Adopt separable covariance structure. Cov mnik!`vdx[zdopx ik Urq%X[Zb\]Bs # t u v w.x u!`v Syqzy{ x u v wa u }Zao S Z\~zy{ u X Covm g! V X[Z o X g! q X[Z \ s # w x!wv Syq%zy{ x v wa Zao S Zb\ zy{ l Take w s to be exponential: w! ƒzy{ # J S { 8 l Decay parameters: { x Xy{ Xy{ x u Xy{ u X are chosen by model choice.

18 Posterior Assume vague prior for t mean parameters. Have Š! ŒXEt and Š ŒXEt ˆ XEt u XEt x Ž, u x u Ž u. and the The joint posterior is a long expression but provides conjugate full conditional distributions for Gibbs sampling. '

19 Predicting urban indicators Want to predict at 3210 new locations. Haven t got the urban indicator h % for those. Use a spatial probit model. Set h % # 7 if U $ and 0 otherwise.! % is a spatial process with mean 6 M -! and variance 1, correlation function w! S = pzy{ š4. Select {ƒš using model choice. Run a Gibbs sampler for this first and save the predictions for!. '

20 Weekly predictions Model is: _!` X[Zd K Š!` X[ZY GHM g!` œx[zd G M UW žik!` X[Zd KX t ˆ Ÿl Bayesian predictive distributions: T < I T d rª # T y G I «9: 9: T ä yª 9: l Need to simulate from: g! X[Z and ik U X[Z l Distributions are normal but parameters are expensive to calculate. Q&

21 _ t t Annual predictions On the sqrt scale. By averaging over 52 weeks: U Š U HM g! M! ik U X ˆ Ō l To obtain the predictive, we need the distributions: g!` K and ik UW K. By sufficiency g!? g!` K <± g U - X=l=l=lfX g!w²6y³. Can show that: g!` K Š µ x U SW x - ' º¹A»½¼¾ d a is equivalent to X ± 7 S x! S` K - x x! S` Ky³PÀ l P

22 _ # È - Ä- Ä - _ M l Original scale orig UW # _ a d a \<Á - m d _ \<Á - UW px[z\] UW px[z\] S!` Â s _!` K is the square of the predictions on the square-root scale. UW Âl For the first, use: - d d \<Á - m q[á - m!` X[Z\P S!` X[Zb\hÅ S Take Æ # Ç weeks (equi-spaced) in four Ã!W Â s _ Ä UW s quarters. Do these predictive calculations for each of É# 24$O$4$6 MCMC iterations. The results for Æ # 2 look similar.

23 3 Ã We obtain annual mean surface, standard deviation and an estimate of the probability that annual PM ^ greater than 15. Thus we store, read and manipulate 7f$ v v 24$O$4$ 96 million numbers. Were we to average over all 52 weeks we had to work with 1000 million or a billion numbers! Q

24 Ê l Analysis: Model Choice Define the predictive model choice criterion Ë Ë (PMCC): # Var T < G I 'M T ä yª S T < I PMCC = P(enalty) + G(oodness of fit) Find the optimal decay parameters. Assume { x # { x u # { x and { # { u # {. Search for { x among (0.015, 0.03, 0.06, 0.12 and 0.30). Search for { among (3, 1.5, 1, 0.75 and 0.6). Optimal { x # $lì$43 (range = 100 km = about 1/10th of the maximal distance). Optimal { # 74l (range = 2 weeks).

25 Analysis: Model comparison Compare the proposed model with the basic model which is not spatio-temporal. PMCC values: Model P G P + G Base model Our model There are other intermediate models. We are not interested in those.

26 Analysis: Parameter estimates mean sd 95% interval Grand mean (4.150, 4.401) Pop-density (0.010, 1.023) Urban (0.244, 0.391) interaction ( 1.124, 0.123) Feb ( 0.788, 0.477) Mar ( 0.864, 0.560) Apr ( 1.365, 1.033) May ( 1.023, 0.706) Jun ( 0.288, 0.046) Jul ( 0.875, 0.526) Aug ( 0.800, 0.496) Sep ( 1.318, 0.977) Oct ( 1.357, 1.018) Nov ( 1.093, 0.784) Í Î Dec ( 1.244, 0.905) Ï (0.057, 0.070) Rural Var (0.363, 0.405) Urban Var (0.001, 0.003) Q

27 Discussion Current work Proposed and fitted non-stationary spatio-temporal model. Positive shift in mean level for urban areas. Also larger variability in urban areas. Obtained both weekly and annual predictions. Future work Produce annual map for the entire US. Simultaneous modeling with other pollutants. See if rural/urban differences are there for other pollutants.

28 Ö!Ø'Ù%ÚºÛaÜRÚÞÝ`ßàHáâÝ äãåø Ø[æ fç è Ø é% Rê1ër Iì + 'í + Ð ÑÂÒpÑ Ó}ÑÕÔIÑ Q

29 Ö% Ø'Ù"ÚîÛaÜ ÚïÝ`ß ànáðýq HãåØ[Ø æ Q& çòñ ë[ì í" P Ð ÑÂÒpÑ Ó}ÑÕÔIÑ Q

30 Ö% Ø'Ù"ÚîÛaÜ ÚïÝ`ß ànáðýq HãåØ[Ø æ ç«ó ê"ôwê àü 'Eí] Q& Ð ÑÂÒpÑ Ó}ÑÕÔIÑ &

31 Ö%!Ø+Ù"ÚºÛžÜ ÚïÝ`ß ànáðýq HãåØ[Ø[æ ç õ Ý+ö`Ø+ éoø Yí4 Ð ÑÂÒpÑ Ó}ÑÕÔIÑ ]

32 à!ù Ýá¾Ö% Ø'Ù"ÚîÛaÜ ÚïÝ`ß àháðýq NãåØ[Ø æ øðgñ Ðùø L

33 à Ù Ýá¾Ö%!Ø+Ù"ÚºÛžÜ ÚïÝ`ß àháðýq NãåØ[Ø[æ Q& øðgñ Ðùø

34 à Ù Ýá¾Ö%!Ø+Ù"ÚºÛžÜ ÚïÝ`ß àháðýq NãåØ[Ø[æ øðgñ Ðùø

35 à Ù ÝLá¾Ö%!Ø+Ù"ÚºÛžÜ ÚÞÝWß ànáâý Hã Ø Ø[æ ç õ Ý+ö`Ø+ é4ø[ Yí øðgñ Ðùø L

36 ëqß ßBê1ëüû0Ö%!Ø+Ù"ÚºÛžÜ ÚïÝ`ß ànáðýq Q&&] Ð ÓaÐ ÔbÐùødÐGú

37 à Ù ÝLáòëQß ßBê1ëüû6Ö%!Ø+Ù"ÚºÛžÜ ÚÞÝWß à Ð ýþô1ò!ýÿñ1ò!ýÿú

38 !Ýé Ü ëyü ëqß ßBê1ëüû ñ + Ñ ý Ò8Ñ!ý ø8ñ ý

Fusing point and areal level space-time data. data with application to wet deposition

Fusing point and areal level space-time data. data with application to wet deposition Fusing point and areal level space-time data with application to wet deposition Alan Gelfand Duke University Joint work with Sujit Sahu and David Holland Chemical Deposition Combustion of fossil fuel produces

More information

GAMINGRE 8/1/ of 7

GAMINGRE 8/1/ of 7 FYE 09/30/92 JULY 92 0.00 254,550.00 0.00 0 0 0 0 0 0 0 0 0 254,550.00 0.00 0.00 0.00 0.00 254,550.00 AUG 10,616,710.31 5,299.95 845,656.83 84,565.68 61,084.86 23,480.82 339,734.73 135,893.89 67,946.95

More information

NRCSE. Misalignment and use of deterministic models

NRCSE. Misalignment and use of deterministic models NRCSE Misalignment and use of deterministic models Work with Veronica Berrocal Peter Craigmile Wendy Meiring Paul Sampson Gregory Nikulin The choice of spatial scale some questions 1. Which spatial scale

More information

SYSTEM BRIEF DAILY SUMMARY

SYSTEM BRIEF DAILY SUMMARY SYSTEM BRIEF DAILY SUMMARY * ANNUAL MaxTemp NEL (MWH) Hr Ending Hr Ending LOAD (PEAK HOURS 7:00 AM TO 10:00 PM MON-SAT) ENERGY (MWH) INCREMENTAL COST DAY DATE Civic TOTAL MAXIMUM @Max MINIMUM @Min FACTOR

More information

Time Series Analysis

Time Series Analysis Time Series Analysis A time series is a sequence of observations made: 1) over a continuous time interval, 2) of successive measurements across that interval, 3) using equal spacing between consecutive

More information

SYSTEM BRIEF DAILY SUMMARY

SYSTEM BRIEF DAILY SUMMARY SYSTEM BRIEF DAILY SUMMARY * ANNUAL MaxTemp NEL (MWH) Hr Ending Hr Ending LOAD (PEAK HOURS 7:00 AM TO 10:00 PM MON-SAT) ENERGY (MWH) INCREMENTAL COST DAY DATE Civic TOTAL MAXIMUM @Max MINIMUM @Min FACTOR

More information

Technical note on seasonal adjustment for M0

Technical note on seasonal adjustment for M0 Technical note on seasonal adjustment for M0 July 1, 2013 Contents 1 M0 2 2 Steps in the seasonal adjustment procedure 3 2.1 Pre-adjustment analysis............................... 3 2.2 Seasonal adjustment.................................

More information

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

The Importance of Ammonia in Modeling Atmospheric Transport and Deposition of Air Pollution. Organization of Talk: The Importance of Ammonia in Modeling Atmospheric Transport and Deposition of Air Pollution Organization of Talk: What is modeled Importance of NH 3 emissions to deposition Status of NH 3 emissions (model-based)

More information

Framework for functional tree simulation applied to 'golden delicious' apple trees

Framework for functional tree simulation applied to 'golden delicious' apple trees Purdue University Purdue e-pubs Open Access Theses Theses and Dissertations Spring 2015 Framework for functional tree simulation applied to 'golden delicious' apple trees Marek Fiser Purdue University

More information

3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 REDCLIFF MUNICIPAL PLANNING COMMISSION FOR COMMENT/DISCUSSION DATE: TOPIC: April 27 th, 2018 Bylaw 1860/2018, proposed amendments to the Land Use Bylaw regarding cannabis

More information

ENGINE SERIAL NUMBERS

ENGINE SERIAL NUMBERS ENGINE SERIAL NUMBERS The engine number was also the serial number of the car. Engines were numbered when they were completed, and for the most part went into a chassis within a day or so. However, some

More information

Spatiotemporal variations of alpine climate, snow cover and phenology

Spatiotemporal variations of alpine climate, snow cover and phenology Spatiotemporal variations of alpine climate, snow cover and phenology S. Asam, M. Callegari, M. Matiu, G. Fiore, L. De Gregorio, A. Jacob, A. Menzel, C. Notarnicola, M. Zebisch Asam et al., Spatiotemporal

More information

Bayesian dynamic modeling for large space-time weather datasets using Gaussian predictive processes

Bayesian dynamic modeling for large space-time weather datasets using Gaussian predictive processes Bayesian dynamic modeling for large space-time weather datasets using Gaussian predictive processes Sudipto Banerjee 1 and Andrew O. Finley 2 1 Biostatistics, School of Public Health, University of Minnesota,

More information

Eddy and Chlorophyll-a Structure in the Kuroshio Extension Detected from Altimeter and SeaWiFS

Eddy and Chlorophyll-a Structure in the Kuroshio Extension Detected from Altimeter and SeaWiFS 14th Symposium on Integrated Observing and Assimilation Systems for the Atmosphere, Oceans, and Land Surface (IOAS-AOLS), AMS Atlanta, January 17-21, 21 Eddy and Chlorophyll-a Structure in the Kuroshio

More information

Technical note on seasonal adjustment for Capital goods imports

Technical note on seasonal adjustment for Capital goods imports Technical note on seasonal adjustment for Capital goods imports July 1, 2013 Contents 1 Capital goods imports 2 1.1 Additive versus multiplicative seasonality..................... 2 2 Steps in the seasonal

More information

Who is polluting the Columbia River Gorge?

Who is polluting the Columbia River Gorge? Who is polluting the Columbia River Gorge? Final report to the Yakima Nation Prepared by: Dan Jaffe, Ph.D Northwest Air Quality, Inc. 7746 Ravenna Avenue NE Seattle WA 98115 NW_airquality@hotmail.com December

More information

ETIKA V PROFESII PSYCHOLÓGA

ETIKA V PROFESII PSYCHOLÓGA P r a ž s k á v y s o k á š k o l a p s y c h o s o c i á l n í c h s t u d i í ETIKA V PROFESII PSYCHOLÓGA N a t á l i a S l o b o d n í k o v á v e d ú c i p r á c e : P h D r. M a r t i n S t r o u

More information

Bayesian dynamic modeling for large space-time weather datasets using Gaussian predictive processes

Bayesian dynamic modeling for large space-time weather datasets using Gaussian predictive processes Bayesian dynamic modeling for large space-time weather datasets using Gaussian predictive processes Andrew O. Finley 1 and Sudipto Banerjee 2 1 Department of Forestry & Department of Geography, Michigan

More information

Annual Average NYMEX Strip Comparison 7/03/2017

Annual Average NYMEX Strip Comparison 7/03/2017 Annual Average NYMEX Strip Comparison 7/03/2017 To Year to Year Oil Price Deck ($/bbl) change Year change 7/3/2017 6/1/2017 5/1/2017 4/3/2017 3/1/2017 2/1/2017-2.7% 2017 Average -10.4% 47.52 48.84 49.58

More information

Bayesian dynamic modeling for large space-time weather datasets using Gaussian predictive processes

Bayesian dynamic modeling for large space-time weather datasets using Gaussian predictive processes Bayesian dynamic modeling for large space-time weather datasets using Gaussian predictive processes Andrew O. Finley Department of Forestry & Department of Geography, Michigan State University, Lansing

More information

Time series and Forecasting

Time series and Forecasting Chapter 2 Time series and Forecasting 2.1 Introduction Data are frequently recorded at regular time intervals, for instance, daily stock market indices, the monthly rate of inflation or annual profit figures.

More information

WHEN IS IT EVER GOING TO RAIN? Table of Average Annual Rainfall and Rainfall For Selected Arizona Cities

WHEN IS IT EVER GOING TO RAIN? Table of Average Annual Rainfall and Rainfall For Selected Arizona Cities WHEN IS IT EVER GOING TO RAIN? Table of Average Annual Rainfall and 2001-2002 Rainfall For Selected Arizona Cities Phoenix Tucson Flagstaff Avg. 2001-2002 Avg. 2001-2002 Avg. 2001-2002 October 0.7 0.0

More information

Climate also has a large influence on how local ecosystems have evolved and how we interact with them.

Climate also has a large influence on how local ecosystems have evolved and how we interact with them. The Mississippi River in a Changing Climate By Paul Lehman, P.Eng., General Manager Mississippi Valley Conservation (This article originally appeared in the Mississippi Lakes Association s 212 Mississippi

More information

= observed volume on day l for bin j = base volume in jth bin, and = residual error, assumed independent with mean zero.

= observed volume on day l for bin j = base volume in jth bin, and = residual error, assumed independent with mean zero. QB research September 4, 06 Page -Minute Bin Volume Forecast Model Overview In response to strong client demand, Quantitative Brokers (QB) has developed a new algorithm called Closer that specifically

More information

Interannual variation of MODIS NDVI in Lake Taihu and its relation to climate in submerged macrophyte region

Interannual variation of MODIS NDVI in Lake Taihu and its relation to climate in submerged macrophyte region Yale-NUIST Center on Atmospheric Environment Interannual variation of MODIS NDVI in Lake Taihu and its relation to climate in submerged macrophyte region ZhangZhen 2015.07.10 1 Outline Introduction Data

More information

Optimal Control of PDEs

Optimal Control of PDEs Optimal Control of PDEs Suzanne Lenhart University of Tennessee, Knoville Department of Mathematics Lecture1 p.1/36 Outline 1. Idea of diffusion PDE 2. Motivating Eample 3. Big picture of optimal control

More information

Disentangling Impacts of Climate & Land Use Changes on the Quantity & Quality of River Flows in Southern Ontario

Disentangling Impacts of Climate & Land Use Changes on the Quantity & Quality of River Flows in Southern Ontario Disentangling Impacts of Climate & Land Use Changes on the Quantity & Quality of River Flows in Southern Ontario by Trevor Dickinson & Ramesh Rudra, Water Resources Engineering University of Guelph Acknowledgements

More information

All work must be shown in this course for full credit. Unsupported answers may receive NO credit.

All work must be shown in this course for full credit. Unsupported answers may receive NO credit. AP Calculus. Worksheet All work must be shown in this course for full credit. Unsupported answers may receive NO credit.. Write the equation of the line that goes through the points ( 3, 7) and (4, 5)

More information

Record date Payment date PID element Non-PID element. 08 Sep Oct p p. 01 Dec Jan p 9.85p

Record date Payment date PID element Non-PID element. 08 Sep Oct p p. 01 Dec Jan p 9.85p 2017/18 Record date Payment date PID element Non-PID element 08 Sep 17 06 Oct 17 9.85p - 9.85p 01 Dec 17 05 Jan 18-9.85p 9.85p 09 Mar 18 06 Apr 18 9.85p - 9.85p Final 22 Jun 18 27 Jul 18 14.65p - 14.65p

More information

Variability and trends in daily minimum and maximum temperatures and in diurnal temperature range in Lithuania, Latvia and Estonia

Variability and trends in daily minimum and maximum temperatures and in diurnal temperature range in Lithuania, Latvia and Estonia Variability and trends in daily minimum and maximum temperatures and in diurnal temperature range in Lithuania, Latvia and Estonia Jaak Jaagus Dept. of Geography, University of Tartu Agrita Briede Dept.

More information

Predicting Long-term Exposures for Health Effect Studies

Predicting Long-term Exposures for Health Effect Studies Predicting Long-term Exposures for Health Effect Studies Lianne Sheppard Adam A. Szpiro, Johan Lindström, Paul D. Sampson and the MESA Air team University of Washington CMAS Special Session, October 13,

More information

Outline. The Path of the Sun. Emissivity and Absorptivity. Average Radiation Properties II. Average Radiation Properties

Outline. The Path of the Sun. Emissivity and Absorptivity. Average Radiation Properties II. Average Radiation Properties The Path of the Sun Larry Caretto Mechanical Engineering 83 Alternative Energy Engineering II March, 2 Outline Review radiation properties for solar collectors Orientation of earth and sun Earth-based

More information

What is the difference between Weather and Climate?

What is the difference between Weather and Climate? What is the difference between Weather and Climate? Objective Many people are confused about the difference between weather and climate. This makes understanding the difference between weather forecasts

More information

ESTIMATING THE MEAN LEVEL OF FINE PARTICULATE MATTER: AN APPLICATION OF SPATIAL STATISTICS

ESTIMATING THE MEAN LEVEL OF FINE PARTICULATE MATTER: AN APPLICATION OF SPATIAL STATISTICS ESTIMATING THE MEAN LEVEL OF FINE PARTICULATE MATTER: AN APPLICATION OF SPATIAL STATISTICS Richard L. Smith Department of Statistics and Operations Research University of North Carolina Chapel Hill, N.C.,

More information

Jackson County 2018 Weather Data 67 Years of Weather Data Recorded at the UF/IFAS Marianna North Florida Research and Education Center

Jackson County 2018 Weather Data 67 Years of Weather Data Recorded at the UF/IFAS Marianna North Florida Research and Education Center Jackson County 2018 Weather Data 67 Years of Weather Data Recorded at the UF/IFAS Marianna North Florida Research and Education Center Doug Mayo Jackson County Extension Director 1952-2008 Rainfall Data

More information

Evapo-transpiration Losses Produced by Irrigation in the Snake River Basin, Idaho

Evapo-transpiration Losses Produced by Irrigation in the Snake River Basin, Idaho Nov 7, 2007 DRAFT Evapo-transpiration Losses Produced by Irrigation in the Snake River Basin, Idaho Wendell Tangborn and Birbal Rana HyMet Inc. Vashon Island, WA Abstract An estimated 8 MAF (million acre-feet)

More information

Multivariate spatial-temporal modeling and prediction of speciated fine particles 1

Multivariate spatial-temporal modeling and prediction of speciated fine particles 1 Multivariate spatial-temporal modeling and prediction of speciated fine particles 1 Jungsoon Choi, Brian J. Reich, Montserrat Fuentes, and Jerry M. Davis Abstract Fine particulate matter (PM 2.5 ) is an

More information

Over the course of this unit, you have learned about different

Over the course of this unit, you have learned about different 70 People and Weather TA L K I N G I T O V E R Over the course of this unit, you have learned about different aspects of earth s weather and atmosphere. Atmospheric scientists, climatologists, hydrologists,

More information

Principal Secretary to Government Haryana, Town & Country Planning Department, Haryana, Chandigarh.

Principal Secretary to Government Haryana, Town & Country Planning Department, Haryana, Chandigarh. 1 From To Principal Secretary to Government Haryana, Town & Country Planning Department, Haryana, Chandigarh. The Director General, Town & Country Planning Department, Haryana, Chandigarh. Memo No. Misc-2339

More information

Standardized Anomaly Model Output Statistics Over Complex Terrain.

Standardized Anomaly Model Output Statistics Over Complex Terrain. Standardized Anomaly Model Output Statistics Over Complex Terrain Reto.Stauffer@uibk.ac.at Outline statistical ensemble postprocessing introduction to SAMOS new snow amount forecasts in Tyrol sub-seasonal

More information

Local Ctimatotogical Data Summary White Hall, Illinois

Local Ctimatotogical Data Summary White Hall, Illinois SWS Miscellaneous Publication 98-5 STATE OF ILLINOIS DEPARTMENT OF ENERGY AND NATURAL RESOURCES Local Ctimatotogical Data Summary White Hall, Illinois 1901-1990 by Audrey A. Bryan and Wayne Armstrong Illinois

More information

TILT, DAYLIGHT AND SEASONS WORKSHEET

TILT, DAYLIGHT AND SEASONS WORKSHEET TILT, DAYLIGHT AND SEASONS WORKSHEET Activity Description: Students will use a data table to make a graph for the length of day and average high temperature in Utah. They will then answer questions based

More information

Jackson County 2013 Weather Data

Jackson County 2013 Weather Data Jackson County 2013 Weather Data 61 Years of Weather Data Recorded at the UF/IFAS Marianna North Florida Research and Education Center Doug Mayo Jackson County Extension Director 1952-2008 Rainfall Data

More information

F O R SOCI AL WORK RESE ARCH

F O R SOCI AL WORK RESE ARCH 7 TH EUROPE AN CONFERENCE F O R SOCI AL WORK RESE ARCH C h a l l e n g e s i n s o c i a l w o r k r e s e a r c h c o n f l i c t s, b a r r i e r s a n d p o s s i b i l i t i e s i n r e l a t i o n

More information

Temporal Trends in Forest Fire Season Length

Temporal Trends in Forest Fire Season Length Temporal Trends in Forest Fire Season Length Alisha Albert-Green aalbertg@sfu.ca Department of Statistics and Actuarial Science Simon Fraser University Stochastic Modelling of Forest Dynamics Webinar March

More information

Changes in spatial distribution of chub mackerel under climate change: the case study using Japanese purse seine fisheries data in the East China Sea

Changes in spatial distribution of chub mackerel under climate change: the case study using Japanese purse seine fisheries data in the East China Sea Changes in spatial distribution of chub mackerel under climate change: the case study using Japanese purse seine fisheries data in the East China Sea Tohya Yasuda, Ryuji Yukami, Seiji Ohshimo Seikai National

More information

AREP GAW. AQ Forecasting

AREP GAW. AQ Forecasting AQ Forecasting What Are We Forecasting Averaging Time (3 of 3) PM10 Daily Maximum Values, 2001 Santiago, Chile (MACAM stations) 300 Level 2 Pre-Emergency Level 1 Alert 200 Air Quality Standard 150 100

More information

An Example file... log.txt

An Example file... log.txt # ' ' Start of fie & %$ " 1 - : 5? ;., B - ( * * B - ( * * F I / 0. )- +, * ( ) 8 8 7 /. 6 )- +, 5 5 3 2( 7 7 +, 6 6 9( 3 5( ) 7-0 +, => - +< ( ) )- +, 7 / +, 5 9 (. 6 )- 0 * D>. C )- +, (A :, C 0 )- +,

More information

Funding provided by NOAA Sectoral Applications Research Project CLIMATE. Basic Climatology Colorado Climate Center

Funding provided by NOAA Sectoral Applications Research Project CLIMATE. Basic Climatology Colorado Climate Center Funding provided by NOAA Sectoral Applications Research Project CLIMATE Basic Climatology Colorado Climate Center Remember These? Factor 1: Our Energy Source Factor 2: Revolution & Tilt Factor 3: Rotation!

More information

Chapter 1 Handout: Descriptive Statistics

Chapter 1 Handout: Descriptive Statistics Preview Chapter 1 Handout: Descriptive Statistics Describing a Single Data Variable o Introduction to Distributions o Measure of the Distribution Center: Mean (Average) o Measures of the Distribution Spread:

More information

Dust storm variability over EGYPT By Fathy M ELashmawy Egyptian Meteorological Authority

Dust storm variability over EGYPT By Fathy M ELashmawy Egyptian Meteorological Authority WMO WORKSHOP ON CLIMATE MONITORING INCLUDING THE IMPLEMENTATION OF CLIMATE WATCH SYSTEMS FOR ARAB COUNTRIES IN WEST ASIA, AMMAN, JORDAN, 27-29 MAY 2013 Dust storm variability over EGYPT By Fathy M ELashmawy

More information

Faisal S. Syed, Shahbaz M.,Nadia R.,Siraj I. K., M. Adnan Abid, M. Ashfaq, F. Giorgi, J. Pal, X. Bi

Faisal S. Syed, Shahbaz M.,Nadia R.,Siraj I. K., M. Adnan Abid, M. Ashfaq, F. Giorgi, J. Pal, X. Bi ICTP NCP International Conference on Global Change 15-19 November, 2006, Islamabad Climate Change Studies over South Asia Region Using Regional Climate Model RegCM3 (Preliminary Results) Faisal S. Syed,

More information

CHARACTERISTICS OF TRAFFIC ACCIDENTS IN COLD, SNOWY HOKKAIDO, JAPAN

CHARACTERISTICS OF TRAFFIC ACCIDENTS IN COLD, SNOWY HOKKAIDO, JAPAN CHARACTERISTICS OF TRAFFIC ACCIDENTS IN COLD, SNOWY HOKKAIDO, JAPAN Motoki ASANO Director Traffic Engineering Division Civil Engineering Research Institute of 1-3 Hiragishi, Toyohira-ku, Sapporo, 062-8602,

More information

BAYESIAN PROCESSOR OF ENSEMBLE (BPE): PRIOR DISTRIBUTION FUNCTION

BAYESIAN PROCESSOR OF ENSEMBLE (BPE): PRIOR DISTRIBUTION FUNCTION BAYESIAN PROCESSOR OF ENSEMBLE (BPE): PRIOR DISTRIBUTION FUNCTION Parametric Models and Estimation Procedures Tested on Temperature Data By Roman Krzysztofowicz and Nah Youn Lee University of Virginia

More information

Global Climates. Name Date

Global Climates. Name Date Global Climates Name Date No investigation of the atmosphere is complete without examining the global distribution of the major atmospheric elements and the impact that humans have on weather and climate.

More information

! " # $! % & '! , ) ( + - (. ) ( ) * + / 0 1 2 3 0 / 4 5 / 6 0 ; 8 7 < = 7 > 8 7 8 9 : Œ Š ž P P h ˆ Š ˆ Œ ˆ Š ˆ Ž Ž Ý Ü Ý Ü Ý Ž Ý ê ç è ± ¹ ¼ ¹ ä ± ¹ w ç ¹ è ¼ è Œ ¹ ± ¹ è ¹ è ä ç w ¹ ã ¼ ¹ ä ¹ ¼ ¹ ±

More information

UWM Field Station meteorological data

UWM Field Station meteorological data University of Wisconsin Milwaukee UWM Digital Commons Field Station Bulletins UWM Field Station Spring 992 UWM Field Station meteorological data James W. Popp University of Wisconsin - Milwaukee Follow

More information

THE RESUSPENSION MODEL. Jan Macoun Czech Hydrometeorological Institute (CHMI), Prague, Czech Republic

THE RESUSPENSION MODEL. Jan Macoun Czech Hydrometeorological Institute (CHMI), Prague, Czech Republic THE RESUSPENSION MODEL Jan Macoun Czech Hydrometeorological Institute (CHMI), Prague, Czech Republic INTRODUCTION The particles contained in the ambient air present one of the biggest problems. Their air

More information

United States Climate

United States Climate 25-1 Alabama Annual Average Temperature Data Source: http://www.wrcc.dri.edu/cgi-bin/divplot1_form.pl?0106 25-2 How does climate vary as we traverse the U.S. along 40 N? 25-3 Average Temperature Along

More information

Rob J Hyndman. Forecasting using. 3. Autocorrelation and seasonality OTexts.com/fpp/2/ OTexts.com/fpp/6/1. Forecasting using R 1

Rob J Hyndman. Forecasting using. 3. Autocorrelation and seasonality OTexts.com/fpp/2/ OTexts.com/fpp/6/1. Forecasting using R 1 Rob J Hyndman Forecasting using 3. Autocorrelation and seasonality OTexts.com/fpp/2/ OTexts.com/fpp/6/1 Forecasting using R 1 Outline 1 Time series graphics 2 Seasonal or cyclic? 3 Autocorrelation Forecasting

More information

Comparison of Particulate Monitoring Methods at Fort Air Partnership Monitoring Stations

Comparison of Particulate Monitoring Methods at Fort Air Partnership Monitoring Stations Comparison of Particulate Monitoring Methods at Fort Air Partnership Monitoring Stations Melanie Larsen Harry Benders RS Environmental (Tom Dann) March 13, 2014 Executive Summary Historically FAP has acquired

More information

Changing Hydrology under a Changing Climate for a Coastal Plain Watershed

Changing Hydrology under a Changing Climate for a Coastal Plain Watershed Changing Hydrology under a Changing Climate for a Coastal Plain Watershed David Bosch USDA-ARS, Tifton, GA Jeff Arnold ARS Temple, TX and Peter Allen Baylor University, TX SEWRU Objectives 1. Project changes

More information

Stream Discharge and the Water Budget

Stream Discharge and the Water Budget Regents Earth Science Unit 6: Water Cycle & Climate Name: Lab # Stream Discharge and the Water Budget Introduction: The United States Geological Survey (USGS) measures and publishes values for the daily

More information

Scarborough Tide Gauge

Scarborough Tide Gauge Tide Gauge Location OS: 504898E 488622N WGS84: Latitude: 54 16' 56.990"N Longitude: 00 23' 25.0279"W Instrument Valeport 740 (Druck Pressure Transducer) Benchmarks Benchmark Description TGBM = 4.18m above

More information

QUESTIONS ON QUARKONIUM PRODUCTION IN NUCLEAR COLLISIONS

QUESTIONS ON QUARKONIUM PRODUCTION IN NUCLEAR COLLISIONS International Workshop Quarkonium Working Group QUESTIONS ON QUARKONIUM PRODUCTION IN NUCLEAR COLLISIONS ALBERTO POLLERI TU München and ECT* Trento CERN - November 2002 Outline What do we know for sure?

More information

2016 Meteorology Summary

2016 Meteorology Summary 2016 Meteorology Summary New Jersey Department of Environmental Protection AIR POLLUTION AND METEOROLOGY Meteorology plays an important role in the distribution of pollution throughout the troposphere,

More information

PRELIMINARY DRAFT FOR DISCUSSION PURPOSES

PRELIMINARY DRAFT FOR DISCUSSION PURPOSES Memorandum To: David Thompson From: John Haapala CC: Dan McDonald Bob Montgomery Date: February 24, 2003 File #: 1003551 Re: Lake Wenatchee Historic Water Levels, Operation Model, and Flood Operation This

More information

Models for spatial data (cont d) Types of spatial data. Types of spatial data (cont d) Hierarchical models for spatial data

Models for spatial data (cont d) Types of spatial data. Types of spatial data (cont d) Hierarchical models for spatial data Hierarchical models for spatial data Based on the book by Banerjee, Carlin and Gelfand Hierarchical Modeling and Analysis for Spatial Data, 2004. We focus on Chapters 1, 2 and 5. Geo-referenced data arise

More information

Pre-Calc Chapter 1 Sample Test. D) slope: 3 4

Pre-Calc Chapter 1 Sample Test. D) slope: 3 4 Pre-Calc Chapter 1 Sample Test 1. Use the graphs of f and g to evaluate the function. f( x) gx ( ) (f o g)(-0.5) 1 1 0 4. Plot the points and find the slope of the line passing through the pair of points.

More information

Winter Season Resource Adequacy Analysis Status Report

Winter Season Resource Adequacy Analysis Status Report Winter Season Resource Adequacy Analysis Status Report Tom Falin Director Resource Adequacy Planning Markets & Reliability Committee October 26, 2017 Winter Risk Winter Season Resource Adequacy and Capacity

More information

INDIAN INSTITUTE OF SCIENCE STOCHASTIC HYDROLOGY. Lecture -12 Course Instructor : Prof. P. P. MUJUMDAR Department of Civil Engg., IISc.

INDIAN INSTITUTE OF SCIENCE STOCHASTIC HYDROLOGY. Lecture -12 Course Instructor : Prof. P. P. MUJUMDAR Department of Civil Engg., IISc. INDIAN INSTITUTE OF SCIENCE STOCHASTIC HYDROLOGY Lecture -12 Course Instructor : Prof. P. P. MUJUMDAR Department of Civil Engg., IISc. Summary of the previous lecture Data Extension & Forecasting Moving

More information

Climatography of the United States No

Climatography of the United States No Month (1) Min (2) Month(1) Extremes Lowest (2) Temperature ( F) Lowest Month(1) Degree s (1) Base Temp 65 Heating Cooling 100 Number of s (3) Jan 63.9 39.3 51.6 86 1976 16 56.6 1986 20 1976 2 47.5 1973

More information

Climatography of the United States No

Climatography of the United States No Temperature ( F) Month (1) Min (2) Month(1) Extremes Lowest (2) Lowest Month(1) Degree s (1) Base Temp 65 Heating Cooling 100 Number of s (3) Jan 32.8 21.7 27.3 62 1918 1 35.8 1983-24 1950 29 10.5 1979

More information

Climate Change Impact Assessment on Indian Water Resources. Ashvin Gosain, Sandhya Rao, Debajit Basu Ray

Climate Change Impact Assessment on Indian Water Resources. Ashvin Gosain, Sandhya Rao, Debajit Basu Ray Climate Change Impact Assessment on Indian Water Resources Ashvin Gosain, Sandhya Rao, Debajit Basu Ray Objectives of the Study To quantify the impact of the climate change on the water resources of the

More information

Estimation of speciated and total mercury dry deposition at monitoring

Estimation of speciated and total mercury dry deposition at monitoring SUPPORTING INFORMATION Estimation of speciated and total mercury dry deposition at monitoring locations in eastern and central North America L. Zhang, P. Blanchard, D.A. Gay, E.M. Prestbo, M.R. Risch,

More information

5 - Seasons. Figure 1 shows two pictures of the Sun taken six months apart with the same camera, at the same time of the day, from the same location.

5 - Seasons. Figure 1 shows two pictures of the Sun taken six months apart with the same camera, at the same time of the day, from the same location. ASTR 110L 5 - Seasons Purpose: To plot the distance of the Earth from the Sun over one year and to use the celestial sphere to understand the cause of the seasons. What do you think? Write answers to questions

More information

Introduction to Forecasting

Introduction to Forecasting Introduction to Forecasting Introduction to Forecasting Predicting the future Not an exact science but instead consists of a set of statistical tools and techniques that are supported by human judgment

More information

Exercise 6. Solar Panel Orientation EXERCISE OBJECTIVE DISCUSSION OUTLINE. Introduction to the importance of solar panel orientation DISCUSSION

Exercise 6. Solar Panel Orientation EXERCISE OBJECTIVE DISCUSSION OUTLINE. Introduction to the importance of solar panel orientation DISCUSSION Exercise 6 Solar Panel Orientation EXERCISE OBJECTIVE When you have completed this exercise, you will understand how the solar illumination at any location on Earth varies over the course of a year. You

More information

Long-term Water Quality Monitoring in Estero Bay

Long-term Water Quality Monitoring in Estero Bay Long-term Water Quality Monitoring in Estero Bay Keith Kibbey Laboratory Director Lee County Environmental Laboratory Division of Natural Resource Management Estero Bay Monitoring Programs Three significant

More information

Simulation of Air Quality Using RegCM Model

Simulation of Air Quality Using RegCM Model Simulation of Air Quality Using RegCM Model The Regional Climate Model (RegCM) The Regional Climate Model (RegCM) is one of the RCMs that was originally developed at the National Center for Atmospheric

More information

Comparison of NASA AIRS and MODIS Land Surface Temperature and Infrared Emissivity Measurements from the EOS AQUA platform

Comparison of NASA AIRS and MODIS Land Surface Temperature and Infrared Emissivity Measurements from the EOS AQUA platform Comparison of NASA AIRS and MODIS Land Surface Temperature and Infrared Emissivity Measurements from the EOS AQUA platform Robert Knuteson, Steve Ackerman, Hank Revercomb, Dave Tobin University of Wisconsin-Madison

More information

Climatography of the United States No

Climatography of the United States No Climate Division: AK 5 NWS Call Sign: ANC Month (1) Min (2) Month(1) Extremes Lowest (2) Temperature ( F) Lowest Month(1) Degree s (1) Base Temp 65 Heating Cooling 90 Number of s (3) Jan 22.2 9.3 15.8

More information

TIME SERIES ANALYSIS AND FORECASTING USING THE STATISTICAL MODEL ARIMA

TIME SERIES ANALYSIS AND FORECASTING USING THE STATISTICAL MODEL ARIMA CHAPTER 6 TIME SERIES ANALYSIS AND FORECASTING USING THE STATISTICAL MODEL ARIMA 6.1. Introduction A time series is a sequence of observations ordered in time. A basic assumption in the time series analysis

More information

Climatic and Ecological Conditions in the Klamath Basin of Southern Oregon and Northern California: Projections for the Future

Climatic and Ecological Conditions in the Klamath Basin of Southern Oregon and Northern California: Projections for the Future Climatic and Ecological Conditions in the Klamath Basin of Southern Oregon and Northern California: Projections for the Future A Collaborative Effort by: CLIMATE LEADERSHIP INITIATIVE INSTITUTE FOR A SUSTAINABLE

More information

Bayesian data analysis in practice: Three simple examples

Bayesian data analysis in practice: Three simple examples Bayesian data analysis in practice: Three simple examples Martin P. Tingley Introduction These notes cover three examples I presented at Climatea on 5 October 0. Matlab code is available by request to

More information

Climatography of the United States No

Climatography of the United States No Climate Division: CA 4 NWS Call Sign: Month (1) Min (2) Month(1) Extremes Lowest (2) Temperature ( F) Lowest Month(1) Degree s (1) Base Temp 65 Heating Cooling 100 Number of s (3) Jan 55.6 39.3 47.5 77

More information

Climatography of the United States No

Climatography of the United States No Climate Division: CA 1 NWS Call Sign: Month (1) Min (2) Month(1) Extremes Lowest (2) Temperature ( F) Lowest Month(1) Degree s (1) Base Temp 65 Heating Cooling 100 Number of s (3) Jan 57.9 38.9 48.4 85

More information

Climatography of the United States No

Climatography of the United States No Climate Division: CA 5 NWS Call Sign: Month (1) Min (2) Month(1) Extremes Lowest (2) Temperature ( F) Lowest Month(1) Degree s (1) Base Temp 65 Heating Cooling 100 Number of s (3) Jan 56.6 36.5 46.6 81

More information

Climatography of the United States No

Climatography of the United States No Climate Division: CA 5 NWS Call Sign: Month (1) Min (2) Month(1) Extremes Lowest (2) Temperature ( F) Lowest Month(1) Degree s (1) Base Temp 65 Heating Cooling 100 Number of s (3) Jan 44.8 25.4 35.1 72

More information

Climatography of the United States No

Climatography of the United States No Climate Division: CA 4 NWS Call Sign: Month (1) Min (2) Month(1) Extremes Lowest (2) Temperature ( F) Lowest Month(1) Degree s (1) Base Temp 65 Heating Cooling 100 Number of s (3) Jan 49.4 37.5 43.5 73

More information

Climatography of the United States No

Climatography of the United States No Climate Division: CA 6 NWS Call Sign: Month (1) Min (2) Month(1) Extremes Lowest (2) Temperature ( F) Lowest Month(1) Degree s (1) Base Temp 65 Heating Cooling 100 Number of s (3) Jan 69.4 46.6 58.0 92

More information

Climatography of the United States No

Climatography of the United States No Climate Division: CA 4 NWS Call Sign: Month (1) Min (2) Month(1) Extremes Lowest (2) Temperature ( F) Lowest Month(1) Degree s (1) Base Temp 65 Heating Cooling 1 Number of s (3) Jan 58.5 38.8 48.7 79 1962

More information

Climatography of the United States No

Climatography of the United States No Climate Division: CA 6 NWS Call Sign: Month (1) Min (2) Month(1) Extremes Lowest (2) Temperature ( F) Lowest Month(1) Degree s (1) Base Temp 65 Heating Cooling 1 Number of s (3) Jan 67.5 42. 54.8 92 1971

More information

Climatography of the United States No

Climatography of the United States No Climate Division: CA 1 NWS Call Sign: Month (1) Min (2) Month(1) Extremes Lowest (2) Temperature ( F) Lowest Month(1) Degree s (1) Base Temp 65 Heating Cooling 1 Number of s (3) Jan 57.8 39.5 48.7 85 1962

More information

Jackson County 2014 Weather Data

Jackson County 2014 Weather Data Jackson County 2014 Weather Data 62 Years of Weather Data Recorded at the UF/IFAS Marianna North Florida Research and Education Center Doug Mayo Jackson County Extension Director 1952-2008 Rainfall Data

More information

SO2 observation at the summit of Mt. Fuji

SO2 observation at the summit of Mt. Fuji SO observation at the summit of Mt. Fuji Yasuhito Igarashi 1, Yosuke Sawa 1, Katsuhiro Yoshioka, Hiroshi Takahashi 3, Hidekazu Matsueda 1, Yukiko Dokiya 1 Geochemical Research Department, Meteorological

More information

2018 Annual Review of Availability Assessment Hours

2018 Annual Review of Availability Assessment Hours 2018 Annual Review of Availability Assessment Hours Amber Motley Manager, Short Term Forecasting Clyde Loutan Principal, Renewable Energy Integration Karl Meeusen Senior Advisor, Infrastructure & Regulatory

More information

Climatography of the United States No

Climatography of the United States No Climate Division: ND 8 NWS Call Sign: BIS Month (1) Min (2) Month(1) Extremes Lowest (2) Temperature ( F) Lowest Month(1) Degree s (1) Base Temp 65 Heating Cooling 100 Number of s (3) Jan 21.1 -.6 10.2

More information

Supplementary appendix

Supplementary appendix Supplementary appendix This appendix formed part of the original submission and has been peer reviewed. We post it as supplied by the authors. Supplement to: Lowe R, Stewart-Ibarra AM, Petrova D, et al.

More information

Determine the trend for time series data

Determine the trend for time series data Extra Online Questions Determine the trend for time series data Covers AS 90641 (Statistics and Modelling 3.1) Scholarship Statistics and Modelling Chapter 1 Essent ial exam notes Time series 1. The value

More information