A Physical Analysis of Changes in Australian Pan Evaporation

Technical Report prepared as part of Land & Water Australia Project No. ANU49 titled Agro-ecological implications of change to the terrestrial water balance

Michael L. Roderick & Graham D. Farquhar CRC for Greenhouse Accounting Research School of Biological Sciences The Australian National University

Principal Investigator: Dr Michael L. Roderick CRC for Greenhouse Accounting Research School of Biological Sciences The Australian National University Canberra ACT 0200

File: PhysicEpan_2.doc Date: 10 March 2006

______Roderick & Farquhar CRC for Greenhouse Accounting The Australian National University 1

Summary

Roderick and Farquhar (2004) found that evaporative demand as measured by pan evaporation has, on average, declined throughout Australia since the 1970s. In this study we used Bureau of Meteorology measurements of temperature, humidity, wind speed and solar radiation, and an existing model called the PenPan model to determine why pan evaporation had decreased. We follow a Penman-style approach and calculated the radiative and aerodynamic components of pan evaporation separately. This allowed us to “attribute” the observed trend to one or more contributing factors.

Pan evaporation rate was calculated at 7 sites having long term solar radiation and aerodynamic measurements extending back 30 years. The linear trend in pan evaporation as determined by the model-based calculations agreed with the observed trends within acceptable limits. Using that as a basis we calculated the aerodynamic component at 41 sites having long-term measurements (wind, humidity, temperature) for 1975-2004. The radiative component was estimated as the difference between the observed and aerodynamic trend.

The results showed that in general, there was very little change over time in the overall humidity regime. There was some evidence for a decline in solar radiation in the north-west of Australia that contributed to the observed decline in pan evaporation in that region. Overall, the most important factor was a decline in the wind speed.

Whether the “stilling” is a local phenomenon (e.g. local obstructions causing changes in air flow over the pans) or a regional (or larger scale) change remains unknown.

______Roderick & Farquhar CRC for Greenhouse Accounting The Australian National University 2

Introduction

Roderick and Farquhar (2004) reported that evaporative demand, as measured by pan evaporation rate had, on average, declined in Australia from the early 1970s until

2002. This came as a surprise to many in the climate and natural resource management community because, until that research, the dominant paradigm was that evaporative demand would increase with warming. One possibility to explain the decline was that a mistake had been made in the measurements. However, a special forum hosted by the Australian Academy of Science, involving leading Australian and international scientists, did not identify any problems of that nature (Gifford, 2005).

From physical principles we know that evaporative demand should increase slightly with increases in near-surface air temperature provided that all else is constant.

Assuming that the pan evaporation measurements have been made correctly, then the above-noted observed decline in evaporative demand means that all else is not constant.

In fact, it would have been a surprise to find that Australian pan evaporation had on average increased, because the dominant worldwide trend has been for declines. For example, declines have been reported in the US and former Soviet Union (Peterson et al., 1995; Golubev et al., 2001), China (Liu et al., 2004; Chen et al., 2005), India

(Chattopadhyay and Hulme, 1997), Thailand (Tebakari et al., 2005), Canada (Hesch and Burn, 2005) and New Zealand (Roderick and Farquhar, 2005).

The decline in Australian pan evaporation rate, and hence evaporative demand, is of fundamental importance for understanding changes in water availability. The results

______Roderick & Farquhar CRC for Greenhouse Accounting The Australian National University 3 raise two important questions; (i) is the decline happening all year around or is it specific to a particular season/s?, and (ii) what are the physical reasons for the decline?

The first question has been addressed in a separate report (Roderick and Farquhar,

2006). The results showed that pan evaporation rate was, on average, declining in all seasons. This report addresses the second question – why?

The approach is to compile existing meteorological data and to calculate how pan evaporation rate has changed as a function of changes in the radiative (e.g. solar radiation) and aerodynamic (e.g. humidity, wind speed) regimes. To do that we used data from the Bureau of Meteorology (BoM) and an existing model describing evaporation from a class A pan called the PenPan model (Linacre, 1994). The model- based calculations are then compared to observations.

Theory

The PenPan Model

Linacre (1994) used a Penman-style approach to formulate a combination equation describing the mass and energy balance of a class A pan. He called it the “PenPan” model. The PenPan model is based on a steady-state heat balance. For a class A pan this means that it can be applied for periods of a week or more. Linacre’s formulation applies to a month. The derivation (along with a worked example) is described in

Appendix A – only the main results are listed here. The PenPan equation is (see Table

1 for definition of terms),

______Roderick & Farquhar CRC for Greenhouse Accounting The Australian National University 4

 s  λE pan =   (RN , pan + 6 u F (Ta − Td )) (1)  s + 2.4γ 

Note that Eqn (1) can be rewritten as the sum of radiative ( λE pan,R ) and aerodynamic

(λE pan,A ) terms,

λE pan = λE pan,R + λE pan,A (2) where the respective components are,

 s  λE pan,R =   (RN , pan ) (3a)  s + 2.4γ  and,

 s  λE pan,A =   (6 uF (Ta −Td )) (3b)  s + 2.4γ 

Symbol Units Comments λ J kg-1 Latent heat of vaporization of liquid water (~ 2.4 MJ kg-1) {Estimate in J kg-1 using λ = 106 × (2.5 − 0.0024T) with T in deg C} -2 -1 Epan kg m s Evaporation rate from class A pan -2 RN,pan W m Net radiation of the class A pan s Pa K-1 Derivative of saturation vapour pressure with respect to temperature {Estimate in Pa K-1 using s = 100*(0.5 + 0.01T + 0.0019T 2 ) with T in deg C} γ Pa K-1 Psychrometric constant (~ 67 Pa K-1) {Estimate in Pa K-1 using γ = 100*(0.67 − z 7.2×10−5 ) with z (elevation) in metres} Ta K Air temperature Td K Dew point temperature of air F - Correction for change in air density with altitude {Estimate using F =1.0 − z 8.7×10−5 with z (elevation) in metres} u m s-1 Monthly mean wind speed

Table 1 Definition of terms in the PenPan model ______Roderick & Farquhar CRC for Greenhouse Accounting The Australian National University 5

To calculate the aerodynamic component of pan evaporation requires estimates of

monthly mean air temperature (Ta, Eqn 3b, also used to estimate s, see Table 1), monthly

mean dew point temperature of the air (Td, Eqn 3b), monthly wind run expressed as a monthly mean wind speed (u, Eqn 3b) and elevation (z, used to calculate γ and F). To

calculate the radiative component, RN,pan is estimated using the global solar irradiance,

the site latitude and the albedo of the ground surrounding the pan. See Appendix A for

details.

Two key points about the PenPan model (see Appendix A for details) need to considered

when interpreting the results. First, the PenPan model assumes that the aerodynamic

component varies linearly with wind speed. Second, the PenPan model as currently

implemented assumes that the net long wave irradiance of a pan is constant over space

and time. This appears to be a plausible first approximation for many climatological

purposes. However, for use in studying the enhanced greenhouse effect over long

periods (e.g. 100 years) the assumption may not be realistic. For the 30 year period

(1975-2004) considered here, we assume, in the absence of contrary evidence, that it is a reasonable approximation.

Methods

Data Sources & Site Selection

As a basis for the project we used the 61 BoM sites identified previously as having near-continuous class A pan evaporation (Epan) records since 1975 (Roderick and

______Roderick & Farquhar CRC for Greenhouse Accounting The Australian National University 6

Farquhar, 2004). A further 9 potential sites were identified from the list compiled by the BoM (Jovanovic et al., 2005). This gave 70 sites that were potentially suitable for analysis of changes in pan evaporation and the database was constructed for those sites.

Database Construction

Daily data for pan evaporation, precipitation, radiation, air temperature, dew point temperature, relative humidity and wind speed were purchased in a digital format from the BoM. These data were used to construct daily databases for each separate element. These daily databases were then summarised into monthly databases for the

70 sites over the period 1967-2004. In doing that, we calculated monthly mean daily values of the relevant variables for all months having at least 25 daily observations with a ‘Y’ quality control flag (i.e. ‘quality controlled and considered acceptable’) in the daily BoM database. The final database record was left NULL (i.e. blank) for any month that did satisfy this criterion. The final database was constructed by combining all relevant monthly data (Database Name: LWA_1967_2004_70stns) and contains

31920 records (38 years X 12 months X 70 stns = 31920 records). The database is available (as a text file, as a MS Excel spreadsheet and as a MS Access database) on the CD-ROM that accompanies this report. Details of the database are listed in Table

2.

______Roderick & Farquhar CRC for Greenhouse Accounting The Australian National University 7

Database Data Source Description Field STN_NUM BOM Site Number Year Month nDays Actual number of days in the month Epan_raw IDCJDCO5.200506 (mm d-1) - Monthly mean daily Epan_raw (this are the data as they appear on the BOM CDs and are not adjusted for bird guards) Epan_nobs Actual number of obs used to calculate the mean Pre MADAM (mm mth-1) - Total monthly precipitation Solar_Num_Rs The number of the solar stn (BoM have different site numbers for climate and solar stns) – should be the same as Solar_Num_Rl Rs NCCSOL Vers. 2.209 (J m-2 d-1) - Monthly mean daily global solar irradiance Rs_nobs Actual number of obs used to calculate the mean Solar_Num_Rl The number of the solar stn (BoM have different site numbers for climate and solar stns) – should be the same as Solar_Num_Rs Rl NCCSOL Vers. 2.209 (W m-2) - Monthly mean daily incoming longwave irradiance Rl_nobs Actual number of obs used to calculate the mean Ta IDCJHCO2.200506 (degC) - Monthly mean air temperature Ta_nobs Actual number of obs used to calculate the mean Td IDCJHCO2.200506 (degC) - Monthly mean dew point temperature Td_nobs Actual number of obs used to calculate the mean Rh IDCJHCO2.200506 (%) - Monthly mean relative humidity Rh_nobs Actual number of obs used to calculate the mean U IDCJDCO5.200506 (m s-1) - Monthly mean wind speed U_nobs Actual number of obs used to calculate the mean Tw_max IDCJDCO5.200506 (degC) - Monthly mean daily max T of water in the Class A Pan Tw_max_nobs Actual number of obs used to calculate the mean Tw_min IDCJDCO5.200506 (degC) - Monthly mean daily min T of water in the Class A Pan Tw_min_nobs Actual number of obs used to calculate the mean

Table 2 Details of the database. The data source column lists the BoM product code (CD-ROM) used for each data element. All CD-ROM’s were dated June 2005, except for the radiation data. The radiation database (CD-ROM titled NCCSOL Vers. 2.209, dated 2001) is no longer being updated. The BoM provided an update CD (ASCII files, dated June 2005).

______Roderick & Farquhar CRC for Greenhouse Accounting The Australian National University 8

Description of the Meteorological Data

The pan evaporation, wind speed and pan water temperature data were all taken from the daily pan evaporation database (BoM product code: IDCJDCO5.200506, see

Table 2). Those data were already in a daily format.

The air temperature (Ta), dew point temperature (Td) and relative humidity data (Rh) data were taken from the so-called daily THP (Temperature, Humidity, Pressure) CD-

ROM (BoM product code: IDCJHCO2.200506, see Table 2). This CD-ROM contains temperature, humidity and pressure measurements taken at particular instants of time.

At some sites over some time periods, there may be as few as 2 measurements per day, while at other sites/times, there are up to 8 measurements (i.e. one every 3 hours) per day. In all cases, whenever there is an air temperature measurement there is also a dew point estimate – they come together. The database used here was constructed by averaging all of the available measurements each month. Thus some database elements (Ta, Td, Rh) are averages of 60 measurements, while others are averages of

240 measurements. Further, this is not specific to individual stations, i.e. in a few cases some stations have 60 measurements per month for a period, then 240, then back to 60, and so on, as the local circumstances changed. These situations can be distinguished using the relevant “nobs” fields in the database which are counts of the number of measurements used to calculate the average (Ta_nobs, Td_nobs, Rh_nobs; see Table 2). Similarly, standard BoM protocol is to take measurements based on local standard time. Hence during periods of daylight saving, the local solar time of the observations will show a step change of 1 hour. Because of these issues, any

______Roderick & Farquhar CRC for Greenhouse Accounting The Australian National University 9 individual trends derived from the Ta and Td fields must be considered uncertain until the data are made more uniform with respect to observing practice. The BoM are currently working on the humidity database and future improvements are expected

(Coughlan et al., 2005). However, in the evaporation equations the key factor is the difference between Ta and Td. The database estimates, while by no means perfect, should be a much better estimate of that difference, and are considered sufficient for this initial analysis.

The quality of the solar radiation data is uncertain prior to c. 1993, and especially prior to c. 1988 (Bruce Forgan, pers. comm.). Most measurements made post-1993 have used new and more precise instrumentation and the BoM expect an accuracy of better than c. 1 % for those data. Those data can be identified in the database by having a three digit identifier (in the Solar_Num_Rs field, see Table 2). Prior to 1993, the network used less accurate instrumentation and the BoM expect lower accuracy.

An additional complication is that prior to 1988, all solar radiation measurements were calibrated using a statistical (instead of a physical) approach. The calibration was based on the assumption that the clear sky global solar irradiance on a particular day of the year should remain constant over time (see description in Frick et al., 1987) and that any observed variations were due to sensor drift. The measurements were subsequently adjusted by this assumed sensor drift. It is not possible to re-process the original observations using physically based calibration procedures because the original observations are no longer available (Bruce Forgan, pers. comm.).

The solar calibration procedure is implicitly based on the assumption that there would be no trends in non-cloud atmospheric elements, of which the most important are

______Roderick & Farquhar CRC for Greenhouse Accounting The Australian National University 10 aerosols and water vapour. The aerosol assumption may be plausible in Australia but cannot be evaluated because the original measurements are no longer available (Bruce

Forgan, pers. comm.). The water vapour assumption is most likely to be unreasonable because an increase in atmospheric water vapour is expected due to the enhanced greenhouse effect and should lead to increased absorption of solar radiation (Arking,

1999). Hence any changes in global solar irradiance due to changes in water vapour would have been inadvertently removed, although they should be small over a 30 year period. Any cloud-based effect on global solar irradiance should still be present in the data because the calibration procedure was based on clear-sky irradiance.

The net effect is that the accuracy of solar radiation measurements made prior to 1993, and especially those made prior to 1988, are uncertain. Any downward (or upward) trend due to non-cloud effects (e.g. aerosols, water vapour) would have been removed.

While not entirely satisfactory, this is the best that we can do at the moment.

Bird Guard Corrections

All current BoM evaporation pans are fitted with a standardised bird guard. In the early years, the BoM retro-fitted bird guards on some pans. The installation dates for bird guards were available in the BoM database. All pan evaporation measurements made prior to the bird guard installation were reduced by 7% in line with empirical findings (Hoy and Stephens, 1979; van Dijk, 1985).

Calculated estimates of pan evaporation

______Roderick & Farquhar CRC for Greenhouse Accounting The Australian National University 11

For all months having estimates of global solar irradiance (Rs) and mean air

temperature (Ta), the radiative component of pan evaporation was calculated at each site using the modified PenPan model. There were 5843 such estimates (out of a

maximum possible 31920) from 27 sites. For all months having estimates of air

temperature (Ta), dew point temperature (Td) and wind speed (u), the aerodynamic component of pan evaporation was calculated. There were 21044 such estimates from

58 sites. Where the calculated radiative and aerodynamic components were both available in a given month at a given site, they were added to give a total calculated estimate of monthly pan evaporation. There were 5264 such monthly estimates that came from a total of 26 stations. Unfortunately, of those 26 sites, only 7 have a more or less continuous time series from the 1970s until 2004. The underlying reason for this is that many of the solar radiation measurement stations have been shut down.

The calculated values (see Table 3) are stored in the final database (available on the

CD-ROM).

Field Description STN_NUM BOM Site Number Year Month nDays Actual number of days in the month Epan_adj (mm d-1) - OBSERVATION - observed monthly mean daily Epan that have been adjusted where necessary for the Bird Guard (reduce by 7%) Epan_rad (mm d-1) - CALCULATION - Radiative component of the calculated monthly mean daily Epan. All calculations per a modified PenPan Model. Epan_aero (mm d-1) - CALCULATION - Aerodynamic component of the calculated monthly mean daily Epan. All calculations per a modifed PenPan Model. Epan_calc (mm d-1) - CALCULATION - Epan_rad + Epan_aero when both are available.

Table 3 Details of the pan evaporation calculations held in the database.

______Roderick & Farquhar CRC for Greenhouse Accounting The Australian National University 12

Trend Estimation

Linear regression (ordinary least squares) was used to describe the changes over time.

The slope of the linear regression is denoted the “trend” throughout this report. The regressions were calculated separately at each site for each month, and because a linear regression was used, the annual trend at the site could be calculated as the sum of the monthly trends. The principal advantage of using this procedure is that we avoided the need to estimate any missing data. Where appropriate we have listed the actual number of observations used to estimate the trend. This gives an indication of the amount of missing data.

Results

Evaluation of the modified PenPan model (26 sites)

The modified PenPan model was evaluated at the 26 sites having at least one model- based estimate of pan evaporation. In making the model-based estimate, the same parameters were used at all sites. In particular, the albedo of the pan surroundings was assumed to be 0.22 (typical value for short green grass) at all sites, although we recognise that this would in fact vary between sites, and over time. The subsequent comparisons are summarised in Table 4 and plots showing the comparison at each of the 26 sites are available in Fig. B1 (Appendix B). At some sites, there is evidence of consistent bias (i.e. slope not equal to unity, intercept not equal to zero, see Table 4).

That bias could be easily removed by “tuning” the model parameters - the two main parameters that could be “tuned” are the albedo of the ground surrounding the pan (α,

______Roderick & Farquhar CRC for Greenhouse Accounting The Australian National University 13 see Eqn A10 in Appendix A) and the heat transfer coefficient (the numerical value

“6” in Eqn 3b, also see Eqn A2 in Appendix A). While such “tuning” would have been physically realistic, no attempt was made to do that because improving the fit between calculations and observations was not the aim of the study. The RMS varied across the 26 sites from 12 to 49 mm mth-1 with a mean of about 26 mm mth-1. The

RMS is a measure of disagreement between calculations and observations. In general, in the absence of gross measurement error, the observations should be much more precise than the model-based calculations. In that context, the RMS is a measure of the disagreement due to errors in the meteorological measurements (temperature, humidity, wind, solar irradiance) and also any errors in the model parameters and/or model formulation. The RMS errors reported here are in the same range as those reported by Linacre (1994) and are typical of results in the environmental physics field.

For most of the 26 sites, the model-based calculations explained more than 90% of the variation in the observations but at three sites, the percentage was closer to 80% (see

Table 4); Site 3003 (Broome , R2 = 0.84), Site 14015 (Darwin Airport, R2 =

0.81), Site 32040 (Townsville Aero, R2 = 0.84). There were no “glaring” problems with the meteorological data at these sites and we concluded that this was a reasonable representation of the results. In contrast, there were several sites where the model explained more than 98% of the variation in the observations. These were mostly in relatively arid environments.

______Roderick & Farquhar CRC for Greenhouse Accounting The Australian National University 14

Stn Slope Intercept R2 N RMS mm mth-1 mm mth-1 2012 1.03 -24.8 0.91 235 26.9 3003 1.00 -22.4 0.84 174 28.3 4032 0.96 -7.1 0.90 245 28.2 5007 1.20 -45.8 0.97 190 26.2 6011 1.19 -41.7 0.97 92 19.1 7045 1.16 -46.6 0.97 101 34.5 8051 1.08 -25.9 0.96 363 23.4 9741 0.97 -10.3 0.97 216 16.3 9789 0.99 -19.1 0.99 95 21.2 12038 1.14 -31.9 0.99 142 20.4 14015 0.85 -8.2 0.81 299 42.9 15135 1.13 -61.6 0.90 96 33.4 15590 1.11 -43.4 0.93 374 37.7 16001 1.15 -35.3 0.96 180 33.0 26021 0.98 -8.0 0.98 340 13.2 31011 0.84 11.3 0.90 67 20.9 32040 0.82 -4.9 0.84 123 49.3 36031 1.18 -41.4 0.95 226 31.6 39083 0.94 -16.2 0.93 297 30.1 50052 1.09 -10.1 0.97 85 22.6 61078 0.97 -12.3 0.95 153 20.9 66037 1.11 -24.5 0.97 127 14.8 70014 0.91 -14.4 0.98 110 28.8 72150 1.04 -16.4 0.98 297 16.7 76031 1.12 -23.3 0.99 361 17.8 86071 1.04 -13.7 0.98 83 12.0

Table 4 Evaluation of the modified PenPan model at 26 BoM sites (see Fig. B1 for plots at each of the 26 sites). The slope and intercept are from a linear regression between Y (model-based calculation) and X (observed) pan evaporation rate. N denotes the number of months used in the comparison at each site. RMS is the root- mean-square of the difference between observed and calculated pan evaporation rate.

Detailed inspection of the plots (Fig. B1, Appendix B) shows that at many sites there

is a residual seasonal variation (e.g. Fig. B1, 6011 – Carnarvon Airport, 7045 –

Meekatharra Airport, etc.). This demonstrates that there is a “radiative” process

occurring that is not represented by the model. This is almost certainly due to the

geometric calculations (see Appendix A). In particular, in the original formulation,

Linacre calculated the irradiance of the pan walls relative to the water surface (the P parameter in Eqn A10) as a function of latitude. However, as Linacre noted in the discussion, this should also be a function of time (i.e. the apparent solar position ______Roderick & Farquhar CRC for Greenhouse Accounting The Australian National University 15

changes with latitude and time). Given that we are primarily interested in trends over

more than one year, this will not have an impact on the results. A more important

problem is that the residuals often show correlations through time (see lower right

panels in Fig. B1). Ideally, the residuals should be uncorrelated, i.e. “white noise”.

For example, at Darwin Airport (Site 14015, lower right panel in Fig. B1) there is a

“jump” in c. 1987, followed by a “trend” in the residuals. This indicates either; (a) a

problem with the measurements, and/or (b) that something in the surroundings of the

pan is changing in a systematic way, e.g. changes in albedo, that is not accounted for

in the model. For attribution purposes, it is important that the residuals are

uncorrelated over the period being considered. At Darwin Airport, attribution is still

possible over the period 1977-2004 because a linear regression of the residuals

(observed minus calculated) over that period would have a slope near zero. Further, at

Darwin, the magnitude of the residual trend (i.e. c. 60 mm a-1 from 1987-2000, i.e. c.

5 mm a-2) is much smaller than the trend in the observations (c. - 17 mm a-2). The

situation is similar at Alice Springs (Site 15590, Fig. B1, lower right panel) where

there is a “step” in c. 1980 separating two different trends in the residuals. As for

Darwin, the data at Alice Springs could be used for attribution over the period 1975-

2004 because the overall slope of the residuals over that period is near zero. However, in contrast to Darwin, this is not the case for other periods. For example, at Alice

Springs, there is a systematic trend in the residuals of about + 100 mm a-1 from 1980-

2000 (i.e. c. + 5 mm a-2) that is of the same order as the observed trend in that period.

Alternatively, at Mount Gambier (Site 26021, Fig. B1), the residuals are more or less

uncorrelated through time and the data at this site can be used for attribution purposes.

Future improvements in the model and the BoM data should remove some of these

______Roderick & Farquhar CRC for Greenhouse Accounting The Australian National University 16 residual trends. In the interim, it is prudent to check each station individually if attribution is being attempted.

Trend Calculations at Sites with Long Term Records (7 sites)

Site Name Period dEobs/dt dEcalc/dt dER/dt dEA/dt (mm a-2) (mm a-2) (mm a-2) (mm a-2) 8051 GERALDTON 1969-2004 -5.0 -2.0 (± 2.3) -0.4 -1.6 AIRPORT (n = 432) (n=427) (n=369) (n=429) 14015 DARWIN 1977-2004 -17.0 -16.0 (± 5.3) -6.0 -10.0 AIRPORT (n = 336) (n=331) (n=297) (n=334) 15590 ALICE SPRINGS 1969-2004 10.9 14.1 (± 3.6) 1.4 12.7 AIRPORT (n=432) (n=432) (n=370) (n=428) 26021 MOUNT GAMBIER 1973-2004 -5.4 -7.6 (± 1.4) -0.4 -7.2 AERO (n=384) (n=378) (n=344) (n=383) 39083 ROCKHAMPTON 1973-2004 10.0 8.9 (± 3.3) 3.6 5.3 AERO (n=384) (n=374) (n=307) (n=382) 72150 WAGGA WAGGA 1972-2004 -0.7 3.9 (± 1.8) 0.6 3.3 AMO (n=396) (n=394) (n=299) (n=394) 76031 MILDURA 1974-2004 -5.7 -9.1 (± 2.0) 1.4 -10.5 AIRPORT (n=372) (n=370) (n=363) (n=370)

Table 5 Comparison between observed (dEobs/dt) and calculated (dEcalc/dt) trends in pan evaporation rate at seven sites having long term meteorological records. Note that the calculated trend is the sum of the radiative (dER/dt) and aerodynamic (dEA/dt) components. The maximum number of possible monthly observations is shown in brackets under the period while the actual number of observations used in each calculation are listed under the trend estimates. The (1σ) error estimates (dEcalc/dt) are based on the RMS values (Table 4) at each site. The trends in the shaded columns should agree (within the error estimates). The calculations at Darwin Airport (Site 14015) were begun in 1977 to avoid a problem (discontinuity) in the wind speed measurements at that site.

Of the 26 sites having suitable radiative and aerodynamic observation, there were 7 sites having long term observations. At these 7 “elite” sites, the trends in observed and calculated pan evaporation rates were determined using the longest possible period at each site (Table 5). Using the best example - at Darwin (Site 14015), the calculations for 1977-2004 estimated the trends to be -6.0 mm a-2 and –10.0 mm a-2 for the radiative and aerodynamic component respectively, giving a total calculated trend of -

16.0 mm a-2. By comparison, the observed trend was -17 mm a-2 – excellent

______Roderick & Farquhar CRC for Greenhouse Accounting The Australian National University 17 agreement. The agreement at the seven “elite” sites is summarised in Table 5. In all cases the observed and calculated trends are within the error estimates. While one can say little about regional trends from these 7 sites, this was still a very important result because it gives an independent confirmation that the results based on the BoM measurements and the modified PenPan model are generally reliable for the period considered.

Estimates of the aerodynamic and radiative trends (41 sites)

The major limitation to the attribution study is the lack of long term solar radiation measurements. However, of the 66 sites, there are 41 that have suitable “near- continuous” measurements to calculate the aerodynamic component. At each of the 41 sites, the trend in the observations, and the model-based calculations of the aerodynamic component were estimated. The trend in the radiative component was then estimated as the difference between the observed trend and the aerodynamic trend at the 41 sites (Fig. 1).

One interesting feature (Fig. 1) in the trend maps is the two inland sites showing conspicuously large increasing trends. Those site locations are:

(133.89°, -23.80°, Site 15590) and Woomera Aerodrome (136.81°, -31.16°, Site

16001). As noted previously (Roderick and Farquhar, 2006), the large increase at

Alice Springs is due to the very low Epan values recorded during a very wet period in the mid 1970s – around the starting year for the calculations. If the calculations at

Alice Springs are begun at a different date (e.g. either 1970 or 1980) the sharp upward

______Roderick & Farquhar CRC for Greenhouse Accounting The Australian National University 18 trend is not nearly as strong. However, the trend at Woomera appears to be a realistic record of changes at that pan.

Site Map 1975-2004 Epan OBS

+ 10 mm a-2 - 10 mm a-2

1975-2004 Epan AERO 1975-2004 Epan RAD

+ 10 mm a-2 + 10 mm a-2 - 10 mm a-2 - 10 mm a-2

Figure 1 Results of the attribution analysis for 1975-2004. Locations of the 41 sites (top left) and the linear trends in: observed pan evaporation rate (top right), aerodynamic component of pan evaporation rate (lower left), radiative component of pan evaporation rate (lower right). Note that the radiative component is computed as the difference between the observed and aerodynamic trends. Black (grey) dots denote decreasing (increasing) trends with the dot diameter scaled per the legend in the bottom left.

______Roderick & Farquhar CRC for Greenhouse Accounting The Australian National University 19

The results show that at many sites, the observed trend is mostly due to changes in the

aerodynamic component, although declines in the radiative component are prominent

in the north-west region (Fig. 1).

Segmenting the Aerodynamic trend into Wind and VPD components

The decline in the aerodynamic component (Fig. 1) raises the question: is this due to

changes in the humidity regime and/or the wind regime? To address that question we

partitioned the aerodynamic trend into two components; one due to changes in wind

(u) and a second due to changes in VPD (here represented as Ta-Td). This was

estimated by differentiating Eqns 3a and 3b,

dE pan,A  s   du   s   d(Ta −Td )  λ ≈   6 F (Ta −Td )  +   6 F u  (4) dt  s + 2.4γ   dt   s + 2.4γ   dt 

where the first term is an estimate of the trend due to changes in wind speed while the

second term is that due to changes in the VPD. The linear trends were estimated from

the trends in the BoM observations ( d(Ta − Td ) dt , du dt ) (see Fig. 2) and the

calculations were completed separately for each month using appropriate monthly

mean values ( Ta − Td , u ) for the 1975-2004 period. Note that in partitioning the trends, we have assumed that s (s + 2.4γ ) (which is a function of temperature) has no

trend as well as neglecting higher order terms (Eqn 4). Because of those

approximations, the sum of the humidity and wind terms will not exactly equal the

calculated aerodynamic component. However, the differences are small and do not

alter the final interpretation (see Table C1, Appendix C). The results are summarised

in Fig. 3. A full listing for the 41 sites is also available in Table C1.

______Roderick & Farquhar CRC for Greenhouse Accounting The Australian National University 20

1975-2004 Ta-Td 1975-2004 U

+ 0.2 K a-1 + 0.05 m s-1 a-1 - 0.2 K a-1 - 0.05 m s-1 a-1

Figure 2 Observed trends in humidity (Ta-Td, left panel) and wind (u, right panel) regimes at 41 sites for 1975-2004.

1975-2004 Epan AERO

+ 10 mm a-2 - 10 mm a-2

1975-2004 AERO: Ta-Td 1975-2004 AERO: U

+ 10 mm a-2 + 10 mm a-2 - 10 mm a-2 - 10 mm a-2

Figure 3 Partition of the aerodynamic trend (top panel) into VPD (lower left) and wind (lower right) components at 41 sites for 1975-2004.

______Roderick & Farquhar CRC for Greenhouse Accounting The Australian National University 21

Summary and Synthesis of the Attribution Study

According to the modified PenPan model used here, the changes (mostly declines) in pan evaporation rate reported over Australia since the 1970s are mostly due to declines in the wind speed (Figs 1-3). At several sites, the trend in wind speed, according to the BoM observations can be as large as - 0.05 m s-1 a-1. Over a thirty

year period, this is equivalent to a drop in mean annual wind speed of c. 1.5 m s-1.

This is a relatively large decline given that a typical global mean wind speed is roughly 2 m s-1. The trend in the radiative component was generally smaller although

a decline in the radiative component was evident in the north-west of Australia (Fig.

1). Changes in pan evaporation rate due to changes in the humidity regime were

generally small (Fig. 3).

Discussion & Conclusion

The attribution reported here is dependent on the model-based calculations (i.e. model

and data). The results at the 7 “elite” sites provide support for the attribution with an

error margin of about 4 mm a-2 (see Table 5). With this context in mind, we calculated

the trend in the aerodynamic component of pan evaporation using monthly

observations of air temperature, humidity and wind speed that had been made by the

BoM at 41 sites. The radiative component was estimated as the difference between the

trend in the observed pan evaporation rate and the aerodynamic component. Note that

any error made in the aerodynamic calculations is compensated for in the radiative

component because we obtained that by difference. Ideally one should compute the

______Roderick & Farquhar CRC for Greenhouse Accounting The Australian National University 22 radiative trend from observations but there is a paucity of radiation measurements, both in Australia, and elsewhere. The encouraging results reported here suggest that evaporation pans might be useful “surrogate” radiometers. This is important because there are many more evaporation pans than radiometers, in Australia, and elsewhere.

The results of the above-noted approach showed that a widespread (but not universal) decline in wind speed was the main reason for the widespread (but not universal) decline in pan evaporation. Trends in the radiative component were generally smaller, except in the north-west of Australia, where there is evidence of a decline in solar radiation. This is consistent with BoM observations showing a decline in solar radiation at Darwin (Table C1). Presumably this is due to increased cloud cover as might be expected to accompany the well known increase in rainfall in the north-west of Australia over the last 30 years (e.g. see www.bom.gov.au & www.longpaddock.qld.gov.au). Trends in pan evaporation rate due to changes in the humidity were generally small. This was expected (Roderick and Farquhar, 2002).

A decline in wind speed at an individual site could be due to a change in local exposure (e.g. growing trees that are part of a wind break located upwind of the BoM class A pan). However, the fact that many sites show declines in wind speed (Fig. 2) suggests that it might be a more widespread climate-related phenomenon. If that is true, the reason for this “stilling” remains unknown.

______Roderick & Farquhar CRC for Greenhouse Accounting The Australian National University 23

Acknowledgements

We thank Alan Beswick, John Carter, Greg McKeon and David Rayner from the

Climate Impacts Group at the Queensland Department of Natural Resources and

Mines (QDNRM) for their generous and helpful advice. We also thank Randall

Donohue and Michael Hobbins for useful discussions and especially Alison Saunders for expert assistance in acquiring and processing the data.

______Roderick & Farquhar CRC for Greenhouse Accounting The Australian National University 24

References

Arking A. 1999. The influence of clouds and water vapour on atmospheric absorption. Geophysical Research Letters 26: 2729-2732.

Chattopadhyay N, Hulme M. 1997. Evaporation and potential evapotranspiration in India under conditions of recent and future climate change. Agricultural and Forest Meteorology 87: 55-73.

Chen D, Gao G, Xu C-Y, Guo J, Ren G. 2005. Comparison of the Thornthwaite method and pan data with the standard Penman-Monteith estimates of reference evapotranspiration in China. Climate Research 28: 123-132.

Coughlan M, Braganza K, Collins D, Jones D, Jovanovic B, Trewin B, 2005: Observed climate change in Australia. In: Greenhouse 2005, November 2005, Melbourne, Victoria.

Frick RA, Walsh PJ, Rice SP, Leadbeater M, 1987: Australian Solar radiation Data Handbook. Vol II. End of Grant Report No. 714, National Energy Research, Development and Demonstration Program, Canberra.

Gifford RM (ed), 2005: Pan evaporation: An example of the detection and attribution of trends in climate variables, Australian Academy of Science. Canberra. 86 pp.

Golubev VS, Lawrimore JH, Groisman PY, Speranskaya NA, Zhuravin SA, Menne MJ, Peterson TC, Malone RW. 2001. Evaporation changes over the contiguous United States and the former USSR: A reassessment. Geophysical Research Letters 28: 2665-2668.

Hesch NM, Burn DH, 2005: Analysis of trends in evaporation - Phase I, Technical Report, Agriculture and Agri-Food Canada: Prairie Farm Rehabilitation Administration, Waterloo, Ontario, 127 pp.

Hoy RD, Stephens SK, 1979: Field Study of Lake Evaporation - Analysis of Data from Phase 2 Storages and Summary of Phase 1 and 2, Technical Paper No. 41, Australian Water Resources Council, 177 pp.

Jovanovic B, Jones DA, Nicholls N, 2005: A historical monthly pan-evaporation dataset for Australia. In: Australian Institute of Physics 16th Biennial Congress, February 2005, Canberra, ACT.

Linacre ET. 1994. Estimating U.S. class A pan evaporation from few climate data. Water International 19: 5-14.

Liu B, Xu M, Henderson M, Gong W. 2004. A spatial analysis of pan evaporation trends in China, 1955-2000. Journal of Geophysical Research 109: D15102, doi:10.1029/2004JD004511.

Peterson TC, Golubev VS, Groisman PY. 1995. Evaporation losing its strength. Nature 377: 687-688. ______Roderick & Farquhar CRC for Greenhouse Accounting The Australian National University 25

Roderick ML. 1999. Estimating the diffuse component from daily and monthly measurements of global radiation. Agricultural and Forest Meteorology 95: 169-185.

Roderick ML, Farquhar GD. 2002. The cause of decreased pan evaporation over the past 50 years. Science 298: 1410-1411.

Roderick ML, Farquhar GD. 2004. Changes in Australian pan evaporation from 1970 to 2002. International Journal of Climatology 24: 1077-1090.

Roderick ML, Farquhar GD. 2005. Changes in New Zealand pan evaporation since the 1970s. International Journal of Climatology 25: 2031-2039.

Roderick ML, Farquhar GD, 2006: Seasonal changes in Australian pan evaporation. Land & Water Australia Technical Report, Canberra, February 2006.

Tebakari T, Yoshitani J, Suvanpimol C. 2005. Time-space trend analysis in pan evaporation over Kingdom of Thailand. Journal of Hydrologic Engineering 10: 205- 215.

van Dijk MH. 1985. Reduction in evaporation due to the bird screen used in the Australian class A pan evaporation network. Australian Meteorological Magazine 33: 181-183.

Wallace BB, 1994: Modelling evaporation from a U.S. Weather Bureau class A pan, B.E. (Hons) Thesis. Centre for Water Research, University of Western Australia.

______Roderick & Farquhar CRC for Greenhouse Accounting The Australian National University 26

Appendix A – Linacre’s “PenPan” Equation with minor modifications

Symbol Units Comments λ J kg-1 Latent heat of vaporization of liquid water (~ 2.4 MJ kg-1) {Estimate in J kg-1 using λ = 106 × (2.5 − 0.0024T ) with T in deg C} -2 -1 Epan kg m s Evaporation rate from class A pan -2 RN,pan W m Net radiation of the class A pan ρ kg m-3 Density of air (~ 1.2 kg m-3) C J kg-1 K-1 Specific heat of dry air (~ 1000 J kg-1 K-1) D Pa Vapour pressure deficit s Pa K-1 Derivative of saturation vapour pressure with respect to temperature`{Estimate in Pa K-1 using s = 100*(0.5 + 0.01T + 0.0019T 2 ) with T in deg C} -1 ra s m Resistance to sensible heat loss from water surface (~ 200/u with u {wind speed} in m s-1) γ Pa K-1 Psychrometric constant (~ 67 Pa K-1) {Estimate in Pa K-1 using γ = 100*(0.67 − z 7.2×10−5 ) with z (elevation) in metres} -1 rs s m Additional resistance to sensible heat loss due to walls of pan (~ 280/u with u {wind speed} in m s-1) Ta K Air temperature Td K Dew point temperature of air F - Density correction for change with altitude {Estimate using F =1.0 − z 8.7×10−5 with z (elevation) in metres} -2 RS,pan W m Solar irradiance of pan (i.e. water surface plus walls of the pan) -2 RS W m Global solar irradiance at surface (i.e. the quantity observed by the BoM) B - Fractional reduction in solar irradiance of the water surface due to the bird guard F - Ratio of direct to global solar irradiance (Eqn A11) P - A latitude dependent correction factor (Eqn A10) {Estimate using P = 1.32 + 4×10−4 φ + 8×10−5 φ 2 based on latitude (φ , decimal degrees). NB. φ is absolute value of the latitude} α - Albedo of surface surrounding the pan (0.22 used here)

Table A1 Definition of terms in the PenPan model

______Roderick & Farquhar CRC for Greenhouse Accounting The Australian National University 27

Theory

Linacre (1994) formulated a Penman-type equation for the evaporation from a class A pan (without a bird guard) that is known as the “PenPan” model. The model is based on a steady-state heat balance. In practice, this means that the equation is applicable over time periods longer than a week or so, and Linacre’s formulation was applied to monthly time scales. The basis of that formulation was a combination-type (i.e. Penman-type) equation as follows (see Table A1 for definitions); ρ c D (R + ) N , pan s r λE = a (A1) pan γ r + r 1+ ( a s ) s ra

-1 -1 Linacre used typical heat transfer coefficients, i.e., ra ~ 200/u s m , rs ~ 280/u s m where u is the wind speed. With these values, and keeping SI units throughout ( ρ ~ 1.2 kg m-3 ,c ~ 1000 J kg-1 K -1 ) we have, ρ c (1.2)(1000)u ~ = 6 u (A2) ra 200 Note that Linacre used 6uF (cf. Eqn A2) instead, where F is a (small) correction to account for the dependence of air density on elevation (z, metres), F =1.0 − z 8.7×10−5 (A3) Also, with the above-noted numerical values, we have, r + r 200/u + 280/u a s ~ = 480/ 200 = 2.4 (A4) ra 200/u Linacre also noted that D ~T − T (A5) s a d Combining the above equations, we have (R + 6uF(T − T )) λE = N , pan a d (A6) pan γ 1+ 2.4 s which can be rewritten as  s  λE pan =   (RN , pan + 6uF(Ta − Td )) (A7)  s + 2.4γ  Note that Eqn (A7) can be interpreted as the sum of radiative and aerodynamic terms,

λE pan = λE pan,R + λE pan,A (A8a) where the radiative component is  s  λE pan,R =   (RN , pan ) , (A8b)  s + 2.4γ  and the aerodynamic component is  s  λE pan,A =   (6uF(Ta −Td )) . (A8c)  s + 2.4γ 

______Roderick & Farquhar CRC for Greenhouse Accounting The Australian National University 28

In this project, these two terms were calculated for each month having the appropriate meteorological data. The aerodynamic calculation is straightforward and not described further. The radiative calculation is described in more detail below.

The Radiative Calculations

Linacre used an empirical equation for monthly net irradiance of the pan, -2 RN , pan = 0.71RS , pan − 40.0 (W m ) (A9) where RS , pan is the solar irradiance of the pan (i.e. solar irradiance of the water surface plus the pan walls). Note that Eqn A9 implies that the net long-wave irradiance is -40 W m-2 and is invariant over space and through time. (This could be modified with more research on the energy balance of a pan.) In turn, RS , pan is related to the global solar irradiance (Rs, i.e. the quantity measured by the BoM) by the following equation which Linacre derived based on the geometry of a Class A pan (4 ft diameter, 10 inches deep),

RS , pan = RS (1.42 − b + f (P −1.42) + 0.42 α) (A10) The b term in Eqn A10 is a modification added here to account for the presence of a birdguard (see below), α is the albedo of the ground surrounding the pan (set as 0.22 for all the calculations in this report), P is a latitude dependent factor (see Table A1) and f is the ratio of direct solar irradiance to global solar irradiance. The calculation of f was based on previous work (Roderick, 1999), i.e., R f = −0.11+1.31 S (A11) RO where RO is the solar irradiance at the top of the atmosphere (calculations for Ro are described by Roderick (1999)).

The Birdguard Correction

Linacre’s original formulation did not account for birdguards. Wallace (1994) noted that the shade cast by the standard BoM birdguard reduced the solar irradiance at the water surface by about 5.5% and that this was (more or less) seasonally invariant. It is very close to the 7% reduction due to the BoM birdguard that has been found empirically. The Wallace result was implemented here by setting b (Eqn A10) equal to 0.055. (Note that in the absence of a bird guard, b is 0).

Linacre’s augmented long-wave formulation for dry ground

Linacre’s original PenPan formulation included an “adjustment” that accounts for additional long-wave irradiance received by the pan walls when the ground around the pan is dry. The idea was that when the ground surrounding the pan is dry, evaporation from the ground is reduced, and the ground becomes warmer. However, the pan is never short of water, so when the ground is dry, the temperature of the ground is higher than the water in the pan, and presumably also higher than the metal walls of the pan. (NB. Presumably this only applies during daylight hours.) We originally implemented the formulation for an additional long-wave flux as given by Linacre.

We found that at high pan evaporation rates (i.e. usually when the ground surrounding the pan is dry) the calculated pan evaporation rate was systematically larger than

______Roderick & Farquhar CRC for Greenhouse Accounting The Australian National University 29 observations. The deviation was about the amount of the augmented long-wave irradiance. Hence, the formulation of the augmented long-wave irradiance was not implemented. Possible problems are that while the dry ground might be warmer than the walls of the pan during the day light hours, it can also be cooler at night and over a full diurnal cycle the effect might be diminished (at least partly). Secondly, Linacre’s formulation was based on relating the net long wave flux to the temperature difference between the surroundings and the pan walls and effectively assumed that both radiating surfaces were near black bodies (i.e. long-wave emissivity near unity). However, while the dry ground would be expected to have an emissivity near unity (e.g. 0.85 – 0.95 would be typical), the metal walls of the pan would have a much lower emissivity.

A Worked Example – Calculations using the modified PenPan

SITE DETAILS

Site Name: Site Number: 2012 (latitude -18.2292 decimal deg, elevation 422 m ASL) Birdguard installation: 10 November 1968

OBSERVATIONS

Date: May 1969 (hence bird guard was installed at time of observations) -1 Epan observed: Monthly mean daily average of 6.95 mm d (per database) (× 31 d mth-1 = 215.5 mm mth-1) -2 -1 Rs observed: Monthly mean daily average of 20537142.86 J m d (per database) (/ 86400 s d-1 = 237.7 W m-2) Ta: 24.94 °C (Monthly mean air temperature) (per database) Td: 5.49 °C (Monthly mean dew point temperature of the air) (per database) u: 1.23 m s-1 (Monthly mean wind speed) (per database)

PARAMETERS

α: 0.22 (assumed albedo of the surroundings, i.e. short grass surface) b: 0.055 (shade cast onto water surface by the birdguard, see main text)

CALCULATIONS

(a) Calculate s (T dependent); s = 100× (0.5 + 0.01T + 0.0019T 2 ) with T = 24.94 oC, s = 193.1 Pa K −1

(b) Calculate γ (elevation (z) dependent); γ = 100× (0.67 − z 7.2×10−5 ) with z = 422 metres, γ = 65 Pa K -1

______Roderick & Farquhar CRC for Greenhouse Accounting The Australian National University 30

(c) Calculate F (elevation (z) dependent); F =1.0 − z 8.7×10−5 with z = 422 metres, F = 0.963

(d) Calculate the aerodynamic component of pan evaporation:  s  λE pan,A =   ()6uF(Ta − Td )  s + 2.4γ  193.1 λE = (6×1.23× 0.963× (24.94 − 5.49)) pan,A 193.1+ 2.4× 65 -2 λE pan,A = 77 W m

(e) Now calculate radiative component. First, the top of atmosphere solar irradiance; -2 RO = 318.325 W m (calc per Roderick (1999); also see computer program solarpos.pro on the CD-ROM available with this report).

(f) Ratio of direct to global solar irradiance; R f = −0.11+1.31 S R O 237.7 f = −0.11+1.31 = 0.87 318.325

(g) The latitude factor, P; P = 1.32 + 4×10−4 φ + 8×10−5 φ 2 with φ = −18.2292

P = 1.354

(h) Solar irradiance of pan;

RS , pan = RS (1.42 − b + f (P −1.42) + 0.42α) with b = 0.055, f = 0.87,P =1.354,α = 0.22

RS , pan = 237.7()1.42 − 0.055 + 0.87(1.354 −1.42) + 0.42×0.22 -2 RS , pan = 332.8 W m

(i) Net irradiance of pan; -2 RN , pan = 0.71RS , pan − 40.0 = 0.71× 332.8 − 40.0 = 196.3 W m

(j) Calculate the radiative component of pan evaporation;  s  λE =   ()R pan,R  s + 2.4γ  N , pan   193.1 λE = ()196.3 = 109 W m -2 pan,R 193.1+ 2.4× 65

Calculations all finished.

______Roderick & Farquhar CRC for Greenhouse Accounting The Australian National University 31

(k) Unit conversions between energy-mass-depth; λ = 106 × (2.5 − 0.0024T ) with T = 24.94 oC, λ = 2.44×106 J kg -1 109 E = = 4.467 ×10-5 kg m -2 s -1 pan,R 2.44×106 The calculations are for a month so we need to express it over that period. −5 -2 -1 -1 E pan,R = 4.467 ×10 kg m s × ( 24× 60× 60× 31) s mth -2 -1 E pan,R = 119.6 kg m mth and with the density of liquid water set as 1000 kg m -3 , we convert to a depth equivalent, -1 -1 E pan,R = 0.1196 m mth = 119.6 mm mth

SUMMARY OF RESULTS

For this worked example (Halls Creek Airport, Site 2012, May 1969), the observed pan evaporation rate was 215.5 mm mth-1 (= 6.95 mm d-1 averaged over the month) while the calculated estimate was slightly lower at 204.1 mm mth-1 (= 6.6 mm d-1).

CALCULATION -2 -1 λEpan,A 77 W m Epan,A 84.5 mm mth -2 -1 λEpan,R 109 W m Epan,R 119.6 mm mth -2 -1 λEpan 186 W m Epan 204.1 mm mth OBSERVATION -2 -1 λEpan 196.2 W m Epan 215.5 mm mth

Table A2 Result summary for the worked example

______Roderick & Farquhar CRC for Greenhouse Accounting The Australian National University 32

Appendix B – Evaluation of the modifed PenPan model at 26 sites

Figure B1 (26 pages following) Comparison between observed and model-based calculations of pan evaporation rate at 26 sites. Each page shows: (top panel) observed (line) and calculated (×) pan evaporation rate; (lower left panel) comparison between observed and calculated pan evaporation with dotted line showing linear regression (see Table 4 for statistics) and 1:1 line also shown (full line); (lower right panel) variations in observed minus calculated pan evaporation rate as a function of time.

______Roderick & Farquhar CRC for Greenhouse Accounting The Australian National University 33

2012 - HALLS CREEK AIRPORT 500

400 ) -1 300

200 Epan (mm mth

100

0 1970 1980 1990 2000 Time (Yrs)

2012 - Epan 2012 - Residuals 600 100

500

50 ) 400 -1 ) -1

300 0

CALC (mm mth 200 OBS - CALC (mm mth -50

100

0 -100 0 100 200 300 400 500 600 1965 1970 1975 1980 1985 1990 1995 OBS (mm mth-1) Time (Yrs)

______Roderick & Farquhar CRC for Greenhouse Accounting The Australian National University 34

3003 - BROOME AIRPORT 400

300 ) -1

200 Epan (mm mth

100

0 1970 1980 1990 2000 Time (Yrs)

3003 - Epan 3003 - Residuals 500 150

400 100 ) -1 ) -1 300

50

200 CALC (mm mth OBS - CALC (mm mth 0 100

0 -50 0 100 200 300 400 500 1985 1990 1995 2000 2005 OBS (mm mth-1) Time (Yrs)

______Roderick & Farquhar CRC for Greenhouse Accounting The Australian National University 35

4032 - PORT HEDLAND AIRPORT 500

400 ) -1 300

200 Epan (mm mth

100

0 1970 1980 1990 2000 Time (Yrs)

4032 - Epan 4032 - Residuals 600 80

60 500

) 40 400 -1 ) -1 20 300 0

CALC (mm mth 200 OBS - CALC (mm mth -20

100 -40

0 -60 0 100 200 300 400 500 600 1965 1970 1975 1980 1985 1990 1995 OBS (mm mth-1) Time (Yrs)

______Roderick & Farquhar CRC for Greenhouse Accounting The Australian National University 36

5007 - 500

400 ) -1 300

200 Epan (mm mth

100

0 1970 1980 1990 2000 Time (Yrs)

5007 - Epan 5007 - Residuals 600 150

500 100 ) 400 -1 ) -1 50

300

0 CALC (mm mth 200 OBS - CALC (mm mth

-50 100

0 -100 0 100 200 300 400 500 600 1985 1990 1995 2000 2005 OBS (mm mth-1) Time (Yrs)

______Roderick & Farquhar CRC for Greenhouse Accounting The Australian National University 37

6011 - CARNARVON AIRPORT 400

300 ) -1

200 Epan (mm mth

100

0 1970 1980 1990 2000 Time (Yrs)

6011 - Epan 6011 - Residuals 500 40

400 20 ) -1 ) -1 300

0

200 CALC (mm mth OBS - CALC (mm mth -20 100

0 -40 0 100 200 300 400 500 1986 1988 1990 1992 1994 1996 OBS (mm mth-1) Time (Yrs)

______Roderick & Farquhar CRC for Greenhouse Accounting The Australian National University 38

7045 - MEEKATHARRA AIRPORT 800

600 ) -1

400 Epan (mm mth

200

0 1970 1980 1990 2000 Time (Yrs)

7045 - Epan 7045 - Residuals 800 100

600 50 ) -1 ) -1

400 0 CALC (mm mth OBS - CALC (mm mth 200 -50

0 -100 0 200 400 600 800 1984 1986 1988 1990 1992 1994 1996 OBS (mm mth-1) Time (Yrs)

______Roderick & Farquhar CRC for Greenhouse Accounting The Australian National University 39

8051 - 500

400 ) -1 300

200 Epan (mm mth

100

0 1970 1980 1990 2000 Time (Yrs)

8051 - Epan 8051 - Residuals 600 200

150 500

) 100 400 -1 ) -1 50 300 0

CALC (mm mth 200 OBS - CALC (mm mth -50

100 -100

0 -150 0 100 200 300 400 500 600 1960 1970 1980 1990 2000 2010 OBS (mm mth-1) Time (Yrs)

______Roderick & Farquhar CRC for Greenhouse Accounting The Australian National University 40

9741 - ALBANY AIRPORT 250

200 ) -1 150

100 Epan (mm mth

50

0 1970 1980 1990 2000 Time (Yrs)

9741 - Epan 9741 - Residuals 400 50

40

300 ) -1 30 ) -1

200 20

CALC (mm mth 10 OBS - CALC (mm mth 100

0

0 -10 0 100 200 300 400 1965 1970 1975 1980 1985 1990 OBS (mm mth-1) Time (Yrs)

______Roderick & Farquhar CRC for Greenhouse Accounting The Australian National University 41

9789 - ESPERANCE 300

250

200 ) -1

150

Epan (mm mth 100

50

0 1970 1980 1990 2000 Time (Yrs)

9789 - Epan 9789 - Residuals 400 50

40 300 ) -1 ) -1 30

200

20 CALC (mm mth OBS - CALC (mm mth 100 10

0 0 0 100 200 300 400 1986 1988 1990 1992 1994 1996 OBS (mm mth-1) Time (Yrs)

______Roderick & Farquhar CRC for Greenhouse Accounting The Australian National University 42

12038 - KALGOORLIE-BOULDER AIRPORT 500

400 ) -1 300

200 Epan (mm mth

100

0 1970 1980 1990 2000 Time (Yrs)

12038 - Epan 12038 - Residuals 600 60

500 40 ) 400 -1 20 ) -1

300 0

CALC (mm mth 200 -20 OBS - CALC (mm mth

100 -40

0 -60 0 100 200 300 400 500 600 1978 1980 1982 1984 1986 1988 1990 1992 OBS (mm mth-1) Time (Yrs)

______Roderick & Farquhar CRC for Greenhouse Accounting The Australian National University 43

14015 - DARWIN AIRPORT 300

200 ) -1

Epan (mm mth 100

0 1970 1980 1990 2000 Time (Yrs)

14015 - Epan 14015 - Residuals 400 120

100

300

) 80 -1 ) -1 60 200 40 CALC (mm mth

OBS - CALC (mm mth 20 100

0

0 -20 0 100 200 300 400 1960 1970 1980 1990 2000 2010 OBS (mm mth-1) Time (Yrs)

______Roderick & Farquhar CRC for Greenhouse Accounting The Australian National University 44

15135 - 600

500

400 ) -1

300

Epan (mm mth 200

100

0 1970 1980 1990 2000 Time (Yrs)

15135 - Epan 15135 - Residuals 800 100

600 50 ) -1 ) -1

400 0 CALC (mm mth OBS - CALC (mm mth 200 -50

0 -100 0 200 400 600 800 1996 1998 2000 2002 2004 2006 OBS (mm mth-1) Time (Yrs)

______Roderick & Farquhar CRC for Greenhouse Accounting The Australian National University 45

15590 - ALICE SPRINGS AIRPORT 600

500

400 ) -1

300

Epan (mm mth 200

100

0 1970 1980 1990 2000 Time (Yrs)

15590 - Epan 15590 - Residuals 800 100

50

600 ) -1 0 ) -1

400 -50

CALC (mm mth -100 OBS - CALC (mm mth 200

-150

0 -200 0 200 400 600 800 1960 1970 1980 1990 2000 2010 OBS (mm mth-1) Time (Yrs)

______Roderick & Farquhar CRC for Greenhouse Accounting The Australian National University 46

16001 - WOOMERA AERODROME 600

500

400 ) -1

300

Epan (mm mth 200

100

0 1970 1980 1990 2000 Time (Yrs)

16001 - Epan 16001 - Residuals 800 100

50 600 ) -1 ) -1 0

400

-50 CALC (mm mth OBS - CALC (mm mth 200 -100

0 -150 0 200 400 600 800 1965 1970 1975 1980 1985 1990 OBS (mm mth-1) Time (Yrs)

______Roderick & Farquhar CRC for Greenhouse Accounting The Australian National University 47

26021 - MOUNT GAMBIER AERO 300

250

200 ) -1

150

Epan (mm mth 100

50

0 1970 1980 1990 2000 Time (Yrs)

26021 - Epan 26021 - Residuals 400 40

300 20 ) -1 ) -1

200 0 CALC (mm mth OBS - CALC (mm mth 100 -20

0 -40 0 100 200 300 400 1970 1980 1990 2000 2010 OBS (mm mth-1) Time (Yrs)

______Roderick & Farquhar CRC for Greenhouse Accounting The Australian National University 48

31011 - CAIRNS AERO 300

200 ) -1

Epan (mm mth 100

0 1970 1980 1990 2000 Time (Yrs)

31011 - Epan 31011 - Residuals 400 50

40

300 ) -1 30 ) -1

200 20

CALC (mm mth 10 OBS - CALC (mm mth 100

0

0 -10 0 100 200 300 400 1996 1998 2000 2002 2004 2006 OBS (mm mth-1) Time (Yrs)

______Roderick & Farquhar CRC for Greenhouse Accounting The Australian National University 49

32040 - TOWNSVILLE AERO 400

300 ) -1

200 Epan (mm mth

100

0 1970 1980 1990 2000 Time (Yrs)

32040 - Epan 32040 - Residuals 500 200

400 150 ) -1 ) -1 300

100

200 CALC (mm mth OBS - CALC (mm mth 50 100

0 0 0 100 200 300 400 500 1980 1982 1984 1986 1988 1990 1992 1994 OBS (mm mth-1) Time (Yrs)

______Roderick & Farquhar CRC for Greenhouse Accounting The Australian National University 50

36031 - LONGREACH AERO 600

500

400 ) -1

300

Epan (mm mth 200

100

0 1970 1980 1990 2000 Time (Yrs)

36031 - Epan 36031 - Residuals 600 100

50 ) 400 -1 ) -1 0

-50 CALC (mm mth 200 OBS - CALC (mm mth

-100

0 -150 0 200 400 600 1965 1970 1975 1980 1985 1990 1995 OBS (mm mth-1) Time (Yrs)

______Roderick & Farquhar CRC for Greenhouse Accounting The Australian National University 51

39083 - ROCKHAMPTON AERO 400

300 ) -1

200 Epan (mm mth

100

0 1970 1980 1990 2000 Time (Yrs)

39083 - Epan 39083 - Residuals 500 80

60 400 ) -1 40 ) -1 300

20

200 CALC (mm mth 0 OBS - CALC (mm mth

100 -20

0 -40 0 100 200 300 400 500 1970 1980 1990 2000 2010 OBS (mm mth-1) Time (Yrs)

______Roderick & Farquhar CRC for Greenhouse Accounting The Australian National University 52

50052 - CONDOBOLIN AG RESEARCH STN 500

400 ) -1 300

200 Epan (mm mth

100

0 1970 1980 1990 2000 Time (Yrs)

50052 - Epan 50052 - Residuals 500 40

20 400 ) -1 0 ) -1 300

-20

200 CALC (mm mth -40 OBS - CALC (mm mth

100 -60

0 -80 0 100 200 300 400 500 1978 1980 1982 1984 1986 1988 OBS (mm mth-1) Time (Yrs)

______Roderick & Farquhar CRC for Greenhouse Accounting The Australian National University 53

61078 - WILLIAMTOWN RAAF 400

300 ) -1

200 Epan (mm mth

100

0 1970 1980 1990 2000 Time (Yrs)

61078 - Epan 61078 - Residuals 500 60

400 40 ) -1 ) -1 300 20

200 0 CALC (mm mth OBS - CALC (mm mth

100 -20

0 -40 0 100 200 300 400 500 1970 1975 1980 1985 1990 OBS (mm mth-1) Time (Yrs)

______Roderick & Farquhar CRC for Greenhouse Accounting The Australian National University 54

66037 - AMO 300

200 ) -1

Epan (mm mth 100

0 1970 1980 1990 2000 Time (Yrs)

66037 - Epan 66037 - Residuals 400 40

300 20 ) -1 ) -1

200 0 CALC (mm mth OBS - CALC (mm mth 100 -20

0 -40 0 100 200 300 400 1982 1984 1986 1988 1990 1992 1994 1996 OBS (mm mth-1) Time (Yrs)

______Roderick & Farquhar CRC for Greenhouse Accounting The Australian National University 55

70014 - 400

300 ) -1

200 Epan (mm mth

100

0 1970 1980 1990 2000 Time (Yrs)

70014 - Epan 70014 - Residuals 500 80

400 60 ) -1 ) -1 300 40

200 20 CALC (mm mth OBS - CALC (mm mth

100 0

0 -20 0 100 200 300 400 500 1982 1984 1986 1988 1990 1992 1994 1996 OBS (mm mth-1) Time (Yrs)

______Roderick & Farquhar CRC for Greenhouse Accounting The Australian National University 56

72150 - WAGGA WAGGA AMO 500

400 ) -1 300

200 Epan (mm mth

100

0 1970 1980 1990 2000 Time (Yrs)

72150 - Epan 72150 - Residuals 600 60

500 40 ) 400 -1 20 ) -1

300 0

CALC (mm mth 200 -20 OBS - CALC (mm mth

100 -40

0 -60 0 100 200 300 400 500 600 1970 1980 1990 2000 2010 OBS (mm mth-1) Time (Yrs)

______Roderick & Farquhar CRC for Greenhouse Accounting The Australian National University 57

76031 - 400

300 ) -1

200 Epan (mm mth

100

0 1970 1980 1990 2000 Time (Yrs)

76031 - Epan 76031 - Residuals 500 40

400 20 ) -1 ) -1 300 0

200 -20 CALC (mm mth OBS - CALC (mm mth

100 -40

0 -60 0 100 200 300 400 500 1970 1980 1990 2000 2010 OBS (mm mth-1) Time (Yrs)

______Roderick & Farquhar CRC for Greenhouse Accounting The Australian National University 58

86071 - MELBOURNE REGIONAL OFFICE 300

250

200 ) -1

150

Epan (mm mth 100

50

0 1970 1980 1990 2000 Time (Yrs)

86071 - Epan 86071 - Residuals 400 30

20 300 ) -1 ) -1 10

200

0 CALC (mm mth OBS - CALC (mm mth 100 -10

0 -20 0 100 200 300 400 1978 1980 1982 1984 1986 1988 OBS (mm mth-1) Time (Yrs)

______Roderick & Farquhar CRC for Greenhouse Accounting The Australian National University 59

Appendix C – Summary of Attribution Analysis Site Site Name OBS N RAD N RAD* AERO N U TaTd (mm a-2) (mm a-2) (mm a-2) (mm a-2) (mm a-2) (mm a-2) 2012 HALLS CREEK AIRPORT -12.5 355 -5.6 -6.9 356 -6.5 -0.7 2014 KIMBERLEY RES.STATION -22.0 315 -15.0 -7.0 298 -5.7 -1.4 3003 BROOME AIRPORT -8.9 355 -1.8 -7.1 355 -5.8 -1.0 4032 PORT HEDLAND AIRPORT 8.1 354 17.5 -9.4 352 -8.9 0.0 5007 LEARMONTH AIRPORT -5.9 354 -4.1 -1.7 353 -0.9 -0.6 6011 CARNARVON AIRPORT -0.5 359 0.7 -1.3 358 0.0 -0.9 7045 MEEKATHARRA AIRPORT -17.5 360 3.6 -21.1 360 -17.1 -2.8 8051 GERALDTON AIRPORT -4.1 355 0.0 298 -2.3 -1.8 359 -1.6 0.5 9592 PEMBERTON -7.4 321 -0.1 -7.2 327 -7.1 -0.1 9741 ALBANY AIRPORT 5.3 358 4.4 0.9 357 1.2 -0.2 9789 ESPERANCE 1.1 360 -2.1 3.2 356 1.3 1.7 12038 KALGOORLIE-BOULDER AIRPORT -4.3 359 -5.9 1.6 351 1.2 0.4 13017 GILES METEOROLOGICAL OFFICE 3.7 360 -2.7 6.5 347 11.4 -5.4 14015 DARWIN AIRPORT -16.4 337 -10.8 -5.6 357 -4.7 -0.4 15135 TENNANT CREEK AIRPORT -3.5 360 3.6 -7.0 360 -10.3 2.4 15590 ALICE SPRINGS AIRPORT 25.8 360 2.6 309 0.6 25.2 360 21.7 3.8 16001 WOOMERA AERODROME 22.8 357 -13.9 36.6 356 24.2 11.4 18012 CEDUNA AMO 2.4 360 -1.8 4.2 354 2.7 1.6 18139 POLDA (GUM VIEW) -11.2 350 2.6 -13.8 343 -10.9 -2.3 23343 ROSEDALE (TURRETFIELD RESEARCH -8.4 359 5.2 -13.7 345 -7.0 -6.3 26021 MOUNT GAMBIER AERO -6.1 360 0.0 320 1.5 -7.6 359 -6.7 -0.8 29127 MOUNT ISA AERO -3.9 353 8.4 -12.3 353 -13.7 2.8 31011 CAIRNS AERO 7.1 355 2.8 4.3 346 3.3 0.5 36031 LONGREACH AERO 3.6 360 -3.7 7.3 358 -0.1 5.9 39083 ROCKHAMPTON AERO 10.9 360 4.1 285 6.3 4.6 359 0.2 3.6 41359 OAKEY AERO 7.1 355 11.3 -4.2 348 -1.0 -3.0 44021 CHARLEVILLE AERO 8.7 360 11.0 -2.3 351 -1.5 -0.1 48027 COBAR MO -6.9 360 1.1 -8.0 348 -10.9 3.1 50052 CONDOBOLIN AG RESEARCH STN -3.3 339 12.1 -15.4 346 -2.6 -9.0 51049 TRANGIE RESEARCH STATION AWS 4.0 351 0.1 3.9 330 1.7 2.7 59040 COFFS HARBOUR MO -11.0 360 -0.8 -10.2 353 -9.7 -0.3 61078 WILLIAMTOWN RAAF -1.6 359 0.6 -2.3 360 -4.2 1.9 61089 SCONE SCS 0.7 356 3.8 -3.1 349 0.5 -2.8 61242 CESSNOCK (NULKABA) -7.4 316 4.4 -11.8 315 -10.5 -0.8 61288 LOSTOCK DAM -16.5 335 7.1 -23.6 333 -6.7 -11.3 63005 BATHURST AGRICULTURAL STATION 1.0 352 4.7 -3.7 358 -2.9 -0.6 66037 SYDNEY AIRPORT AMO 1.1 359 0.4 0.7 357 -1.7 1.6 70014 CANBERRA AIRPORT -2.3 356 -0.5 -1.8 353 -1.7 -0.2 72150 WAGGA WAGGA AMO -1.9 359 0.9 265 -3.6 1.7 359 -0.5 2.4 76031 MILDURA AIRPORT -8.8 358 1.0 351 4.4 -13.2 359 -14.4 2.1 82039 RUTHERGLEN RESEARCH -5.5 354 1.1 -6.6 327 -2.1 -2.1 82042 STRATHBOGIE 0.6 311 2.9 -2.3 313 3.7 -4.4 85072 EAST SALE AIRPORT -1.7 360 4.8 -6.5 355 -5.9 -0.4 88023 LAKE EILDON 0.1 359 3.0 -2.9 338 -1.6 -1.2 91104 COMPARISON -1.4 347 0.9 -2.3 356 -4.1 1.6 94069 GROVE (COMPARISON) -3.7 313 -7.4 3.7 316 2.7 0.7

Table 6 Overall summary of attribution analysis for 1975-2004 at 41 sites. OBS denotes observed trend in pan evaporation rate; N is the number of monthly observations used to calculate the trend; RAD denotes observed trend in the radiative component at 6 sites having (near) continuous records over the period; RAD* denotes the difference between OBS and AERO; AERO denotes trend in the aerodynamic component of pan evaporation rate. The AERO component is further partitioned into changes due to changes in wind (U) and and humidity (TaTd) regimes. U is the component of the aerodynamic trend due to changes in mean wind speed; TaTd is the component of the aerodynamic trend due to changes in the mean difference between air and dewpoint temperature (i.e. humidity) (see main text for details).

______Roderick & Farquhar CRC for Greenhouse Accounting The Australian National University