Demo of GeoPEARL

Download Report

Transcript Demo of GeoPEARL

A spatially-distributed model of pesticide leaching Aaldrik Tiktak

Outline of presentation

• Why spatially distributed modelling?

• What is GeoPEARL?

• How is a GeoPEARL assessment carried out?

• Which are the main model outputs?

• Selected Results and Examples • Questions?

• Overview of the User Interface http://www.pearl.alterra.nl

Why spatially distributed modelling in registration?

Explicit consideration of the variability of the leaching process: • where are the vulnerable regions?

• Are specific areas (e.g. drinking water abstraction areas) sufficiently protected?

Leaching in drinking water abstraction areas http://www.pearl.alterra.nl

Why spatially distributed modelling in registration?

Explicit calculation of ‘realistic worst-case conditions’: • calculation of percentiles of the leaching concentration in space • more realistic choice of standard scenario’s, no need for safety factors in registration Standard Scenario http://www.pearl.alterra.nl

What is GeoPEARL?

• GeoPEARL is a spatially-distributed model of pesticide leaching to the groundwater and surface waters in the Netherlands. The core of the GeoPEARL consists of PEARL itself.

http://www.pearl.alterra.nl

The PEARL model - a 1D-column SWAP hydrology

preci pitat ion irrigati on t rans pirati on c rop c alendar s oil ev aporat ion ev aporat ion of int erc ept ed water t hroughf all ponding s urf ac e run-off

PEAR L pesticide fate

dis s ipation at the c rop c anopy v olati lis ation pes tic ide upt ake heat f low water uptak e by plant root s s oil wat er f luxes depos iti on applic at ion was h of f t rans form ation s olid-liquid gas partit ioning c onvec t ion dis persion dif fus ion injec t ion f luct uati ng groundwater lev el lat eral dis c harge to dit ches and f ield-drai ns s eepage leac hing http://www.pearl.alterra.nl

Position of column in landscape PEARL 1D column field ditches water courses < 3m water courses > 3 m Groundwater table Drainage into local surface waters Leaching into the regional groundwater

http://www.pearl.alterra.nl

Why a complex model??????

http://www.pearl.alterra.nl

The plot appraoch

• PEARL is run in a GIS context for a large number of unique combinations (plots) of soil, climate, land-use and hydrogeological units.

http://www.pearl.alterra.nl

Derivation of the unique combinations:

http://www.pearl.alterra.nl

Obtaining the lower boundary of SWAP:

• Lower boundary condition is a prescribed flux condition (so-called ‘Neumann condition’); • Lower boundary of SWAP is upper boundary of National Groundwater Model; • Start with NAGROM and a simple unsaturated zone model; • Adjust the boundary flux by running NAGROM and SWAP iteratively; • Wait until convergence is reached.

http://www.pearl.alterra.nl

Plots are represented by a single vertical column: A plot

http://www.pearl.alterra.nl

Basic soil properties in GeoPEARL Slightly acidic sandy soils (‘higher part’)

pH-KCl!!!!!

Clay soils with near neutral pH-values (‘lower part’) Acidic peat soils

http://www.pearl.alterra.nl

How are GeoPEARL assessments carried out?

• GeoPEARL is run for 250-6405 plots, depending on spatial detail wanted; • The model is run for the area of potential usage only. This area is assumed equivalent to the area of the crop(s) for which a registration is requested; • For each plot and crop-type, the model is run for 26, 46 or 66 years (FOCUS procedure); • For each plot the median leaching concentration in time is extracted and put into a map.

http://www.pearl.alterra.nl

To be specified by the user:

• Substance properties (start with standard dossier) • Application schedule • Number of plots (i.e. spatial detail wanted) • One or more crops for which a registration is requested http://www.pearl.alterra.nl

More attention for relationship between soil properties and substance properties:

• DT

50

in top-soil specified as a function of clay content, pH and organic matter content

DT DT

50 , 50 ,

plot plot

 

DT

50 , min if

DT

50 , min 

DT

50 ,

plot

DT

50 , max

DT

50 , max

DT

50 ,

plot

DT

50 , min

DT

50 ,

plot DT

50 ,

r

 

f om

(

m om

m om

.

r

) 

f l

(

m l

m l

,

r

) 

f p

H (

p

H 

p

H

r

)

DT

50 ,

plot

DT

50 , max • K

F

specified as a function of organic matter and pH:

K F,eq

m om K om

.

eq, ac 1  

K om

,

eq

,

ba M M ba

ac 10

p

H -

p

Ka 

p

H

M M ba

ac 10 pH -

p

Ka 

p

H or as a function of other properties:

K F

,

eq

K F

,

eq

, 0 

m om K om

,

eq

m sand K sand

,

eq

m silt K silt

,

eq

m clay K clay

,

eq

c sox K sox

,

eq

http://www.pearl.alterra.nl

More attention for relationship between soil properties and substance properties:

http://www.pearl.alterra.nl

The target variable produced by the model:

• The 90 th percentile in space and the 50 th percentile in time of the leaching concentration at 1 m depth at the area of

(potential) usage

http://www.pearl.alterra.nl

Further model outputs:

• Maps and frequency diagrams of terms of the water balances • Maps and frequency diagrams of terms of the substance balances, including the leaching, drainage and volatilisation fluxes • Maps of percentiles and frequency diagrams of the leaching concentration in time • Tables of balance terms http://www.pearl.alterra.nl

Example 1: Comparison of leaching of two substances

Ordinary substance pH-dependent sorption High pH Vulnerable sandy soils

Opposite spatial patterns!

http://www.pearl.alterra.nl

Example 2: Comparison of drainage and leaching

Drainage fluxes are highest in clay and peat soils...

Leaching fluxes are highest in sandy soils...

http://www.pearl.alterra.nl

Some conclusions from previous examples:

• Drainage fluxes are generally higher than leaching fluxes.

• For ordinary substances, drainage dominates in the peat and clay soils and leaching in sandy soils.

• The leaching pattern can be quite different for non-ordinary substances (e.g. substances with pH dependent sorption) http://www.pearl.alterra.nl

GeoPEARL considers the area of potential usage:

• Pesticides are usually authorised for particular crops; • The new decision tree therefore considers the area of potential usage, approximated by the area of crops for which a registration is requested • GeoPEARL is distributed with a database of areas for 24 major crops • All results are presented for the area of potential

usage

http://www.pearl.alterra.nl

GeoPEARL crops are not hydrologically different

• During the derivation of the UC’s, only four crops were included (grass, mais, potatoes and forest); • For each plot, the dominant crop is used. Hydrological parameters are assigned to these crops; • The 24 GeoPEARL crops are used for pre processing only; i.e. for calculating the area of use.

http://www.pearl.alterra.nl

Effect of area of use on pesticide authorisation (1)

Maize is grown on slightly acidic sandy soils Potatoes are grown on light sandy clay soils and in ‘the Veenkolonien’ http://www.pearl.alterra.nl

Effect of area of use on pesticide authorisation (2)

Slightly acidic sandy soils Light sandy clay soils and ‘the Veenkolonien’ http://www.pearl.alterra.nl

Substance with pH-dependent sorption behaviour

Effect of area of usage on authorisation (cntd.)

Unsafe in potatoes Safe in maize http://www.pearl.alterra.nl

Some conclusions from previous section:

• GeoPEARL explicitly considers the variation of the leaching process, so that it is now possible to evaluate the leaching for different crop types • Consideration of the area of usage can have a strong effect on the authorisation decision • It is extremely important to investigate the relationship between pesticide properties and soil properties http://www.pearl.alterra.nl

GeoPEARL assessments versus assessment items

• A notifier may wish to register for several crops; • For each crop, GeoPEARL must be run separately; • Combination of crop and application scheme is an assessment item; • The three items together form an assessment; • The 90 th percentile is calculated for the combination of the three crops, i.e. for the assessment.

http://www.pearl.alterra.nl

The GeoPEARL database

http://www.pearl.alterra.nl

Reduction of computation time by using zones:

• To reduce the computation time, plots are taken together; • Clusters of plots are called ‘zones’.

Plots Zones 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 SLICE 1 2 2 3 5 4 5 1 4 5 1 2 3 4 1 2 3 4 Plot without relevant crop Plot with relevant crop http://www.pearl.alterra.nl

Slicing is done on the basis of ‘vulnerability’

• The normal Arc/Info slice operation clusters neighbouring cells; • We cluster plots with adjacent vulnerability indexes; • This will give a more reliable estimation of the 90 th percentile, because we sample in each class of the cumulative frequency distribution.

http://www.pearl.alterra.nl

The vulnerability index for ordinary substances:

http://www.pearl.alterra.nl

How many plots must be included in an assessment?

these two are different!

almost no difference http://www.pearl.alterra.nl

How many plots must be included in an assessment? (cntnd.)

More plots = more detail, but 250 plots is generally enough...

http://www.pearl.alterra.nl

250 versus 6405 plots

10 1 0.1

0.01

spring 0.001

0.0001

0.00001

0.00001

0.0001

10 1 0.001

0.01

0.1

0.1

GeoPEARL_6405 (ug/L)

0.01

0.001

1 10 0.0001

autumn 0.00001

0.00001

0.0001

0.001

0.01

0.1

GeoPEARL_6405 (ug/L)

250 plots is generally enough...

http://www.pearl.alterra.nl

1 10

Some conclusions from previous section:

• The maximum number of zones is 6405 (which is equal to the number of plots); • A first impression of the leaching pattern can be obtained with less plots (e.g. the 30 as used in this course); • At least 250 plots are required for an accurate estimation of the 90 th percentile. This is the

minimum number for a final submission.

http://www.pearl.alterra.nl

Other uses of GeoPEARL

• GeoPEARL has been combined with a database on pesticide usage in order to calculate the actual pesticide loads into surface waters and groundwaters.

• Evaluation of policy plans... • A European version of GeoPEARL (called EuroPEARL) has been developed • EU registration in zones... http://www.pearl.alterra.nl

EuroPEARL example: 80 th percentile of the leaching concentration of FOCUS substance D

High leaching risk

Hotspots occur everywhere in Europe

http://www.pearl.alterra.nl

Possible scales…

Catchment http://www.pearl.alterra.nl

EU

Final remarks

• A model is now available that can be used to assess the leaching into the groundwater and the drainage into surface waters at a regional-scale • The model is operational for: • the Dutch registration procedure • evaluation of policy plans http://www.pearl.alterra.nl

The Graphical User Interface

GeoPEARL User http://www.pearl.alterra.nl

First install the program….

Download geopearl_1_1_1.exe

from the address http://www.pearl.alterra.nl.

Archive this application Logon with administrator rights Run geopearl_1_1_1.exe to install GeoPEARL Launch GeoPEARL to install Interbase If FOCUS PEARL 2.2.2 is installed: launch FOCUS PEARL as well Logon with your own user name Start using GeoPEARL; register at the PEARL website http://www.pearl.alterra.nl

To be specified by the user:

• Substance properties (start with standard dossier) • Application schedule • Number of plots (i.e. spatial detail wanted) • One or more crops for which a registration is requested http://www.pearl.alterra.nl

Getting started with the User Interface

Add o r e dit a substa nce .

Par ameters must be gi ven for pare nt substances an d for metabol ites Add o r e dit app licati on schemes.

Each assessmen t ite m may have i ts own a pplication sche me.

Use the wizard to ge nerate assessments.

An assessmen t consists o f a sin gle substa nce a nd one or asse ssment items.

Sel ect asse ssments for executi on and push the ca lcu late b utton to start the simulations http://www.pearl.alterra.nl

GUI step 1: Main form

http://www.pearl.alterra.nl

GUI step 2: Edit substance properties

http://www.pearl.alterra.nl

GUI step 3: Edit application schedules

http://www.pearl.alterra.nl

GUI step 4: Run wizard, select substance and resolution

http://www.pearl.alterra.nl

GUI step 5: Run wizard, select crops and applications (i.e. add items)

http://www.pearl.alterra.nl

GUI step 6: Run GeoPEARL

http://www.pearl.alterra.nl

The model is running….

http://www.pearl.alterra.nl

GUI step 7: Select output - Assessment as a whole

http://www.pearl.alterra.nl

GUI step 7: View output - Assessment as a whole

http://www.pearl.alterra.nl

GUI step 8: Select output - Assessment item

http://www.pearl.alterra.nl

GUI step 8: View output - Assessment item

http://www.pearl.alterra.nl

GUI step 9: View report (assessment as a whole)

http://www.pearl.alterra.nl

GUI step 10: Analyse the output

http://www.pearl.alterra.nl

GUI step 10: Analyse the output

http://www.pearl.alterra.nl