14-test3 real rainfall_meteo_bron_2008+solute

Upload: marcela-carolina-camacho-martinez

Post on 01-Nov-2015

229 views

Category:

Documents


0 download

DESCRIPTION

Documento de referencia entregado durante el curso de modelacion ambiental

TRANSCRIPT

  • MODELACIONDELTRANSPORTEDESOLUTOANUALPARAUN

    HIETOGRAMACONOCIDO

    EjemplosHYDRUS1DDatosmeteorolgicosdeBron,Francia

    Casodeunafuentefijadeunmetalpocosoluble

  • HyetogramBron,France(2008)Mat. 1

    Test3 =measurementofhourlyrainfall& estimation of potential&estimationofpotentialevapotranspiration

    Soil profile 1 m Soilprofile=1mMat.1=loamMat. 2 = sandyloamMat.2 sandy loam 4observationpoints l f l f h d Initialprofile,uniformpressurehead

    (hini =1m)Bottomb.c.=freedrainageMat. 2

    ottom b.c. free drainage

  • 5Hydraulic Properties: Head vs. Theta

    0

    Hydraulic Properties: log K vs. Theta

    3

    4

    -4

    -2

    0

    1

    2

    10

    -8

    -6

    00.0 0.1 0.2 0.3 0.4 0.5

    Theta [-]

    -100.0 0.1 0.2 0.3 0.4 0.5

    Theta [-]Hydraulic Properties: log K vs. log h

    0

    2 van Genuchten Mualem model (L = 0.5)

    r s (cm-1) n Ks (cm/h)

    6

    -4

    -2 1 0.078 0.43 0.036 1.56 1.042 0.065 0.41 0.075 1.89 4.42

    1 l (bl k)-8

    -6

    -4 -3 -2 -1 0 1 2log(|h|[m])

    1 = loam (black)2 = sandy-loam (bleu)

    log(|h|[m])

  • P & PET [cm/h]P&PET[cm/h]TOY = Time Of the Year, hora calendario (ver archivo Excel) TOY = 0 el 01/01/2008 a 0 horas

    0.0400 1000 2000 3000 4000 5000 6000 7000 8000 9000 10000

    TOY[h]

    0.03

    0.0350.5

    1

    0 015

    0.02

    0.0251.5

    2

    2 5

    E

    T

    P

    [

    c

    m

    /

    h

    ]

    Q

    i

    n

    p

    u

    t

    [

    c

    m

    /

    h

    ] 2008Qinput[cm/h]

    ETP[cm/h]

    0.005

    0.01

    0.0152.5

    3

    3.5

    0

    0 005

    0 1000 2000 3000 4000 5000 6000 7000 8000 9000 10000

    3 5

    4

    TOY[h]Jan Feb Mar Apr May June July Aug Sept Oct Nov DecJan Feb Mar Apr May June July Aug Sept Oct Nov Dec

  • P & PET cummulated100

    120P&PETcummulated

    80

    100

    60

    Q

    i

    n

    p

    u

    t

    [

    c

    m

    /

    h

    ]

    2008CumQinput[cm/h]

    C ETP [ /h]

    20

    40CumETP[cm/h]

    00 1000 2000 3000 4000 5000 6000 7000 8000 9000 10000

    TOY [h]TOY[h]

    Jan Feb Mar Apr May June July Aug Sept Oct Nov Dec

  • 0Profile Information: Water Content

    0

    Profile Information: Pressure Head

    -40

    -20

    -40

    -20

    -80

    -60

    -80

    -60

    -1000.10 0.15 0.20 0.25 0.30 0.35

    Theta [-]

    -100-300 -200 -100 0

    h [cm][ ]

    0.45

    Observation Nodes: Water Content

    0

    Observation Nodes: Pressure Heads

    0.30

    0.35

    0.40

    100

    -50

    0.15

    0.20

    0.25

    -150

    -100

    0.100 2000 4000 6000 8000 10000

    Time [hours]

    -2000 2000 4000 6000 8000 10000

    Time [hours]

  • All Cumulative FluxesTEST 3

    20Cum Potential surface flux

    10Cum. Potential surface flux

    -10

    0

    -20

    0

    C A t l

    Cum. Bottomflux

    -30Cum. Actual surface flux

    -400 2000 4000 6000 8000 10000

    Time [hours]Time [hours]

  • 2.0

    2.5

    Surface Run-Off

    3.0

    3.5

    4.0

    Cum. Surface Run-Off

    0 5

    1.0

    1.5

    1.0

    1.5

    2.0

    2.5

    0.0

    0.5

    0 2000 4000 6000 8000 10000Time [hours]

    0.0

    0.5

    0 2000 4000 6000 8000 10000Time [hours]

    Cum.Evaporation

    80

    100

    Cum. Infiltration

    50

    60

    70

    Cum. Evaporation

    20

    40

    60

    10

    20

    30

    40

    Soil Water Storage

    00 2000 4000 6000 8000 10000

    Time [hours]

    00 2000 4000 6000 8000 10000

    Time [hours]

    22

    24

    26

    16

    18

    20

    140 2000 4000 6000 8000 10000

    Time [hours]

  • Soluto condicindeNOequilibrio(MIM)ECI3-Solute

    Un soluto reactivo en la superficie del suelo constituye unaUnsolutoreactivoenlasuperficiedelsueloconstituyeunafuentedecontaminacinpermanentemobilisable acadalluvia.Laduracindecadapulse=mismotiempoquelaprecipitacincorrespondiente.

    c0=10mg/L=102 mg/cm3 sobreladuracindelasimulacind=1500mg/cm3Disp=1cmf=0.5im =0.1cm3/cm3Kd =0.005mg/cm3 =1s =0.12h1Condicindefrontera:Concentration BCNB.Elordendemagnituddec0 correspondealaconcentraciones

    deelementosmayores(Br,NO3,Cl,Ca,COT,).Para/metalespesadoslasconcentracionessondelordendelg/L.

  • 0.00030

    Observation Nodes: Concentration

    0

    Profile Information: Concentration

    0 00015

    0.00020

    0.00025

    -40

    -20

    0.00005

    0.00010

    0.00015

    -80

    -60

    0.000000 5000 10000

    Time [hours]

    -1000.000 0.002 0.004 0.006 0.008 0.010

    Conc [mg/cm3]

    P fil I f I bil C

    20

    0

    Profile Inform.: Immobile Conc.

    60

    -40

    -20

    100

    -80

    -60

    -1000.000 0.002 0.004 0.006 0.008 0.010

    Conc [mg/cm3]

  • 0.008

    0.010

    Surface Concentration

    1.5

    2.0

    Surface Solute Flux

    0.4

    0.5

    Cumulative Surface Solute Flux

    0 002

    0.004

    0.006

    0 0

    0.5

    1.0

    0.2

    0.3

    0.000

    0.002

    0 5000 10000Time [hours]

    -0.5

    0.0

    0 2000 4000 6000 8000 10000Time [hours]

    0.0

    0.1

    0 2000 4000 6000 8000 10000Time [hours]

    0 20

    All Solute Fluxes

    0 5

    All Cumulative Solute Fluxes

    0.15

    0.20

    0.4

    0.5

    0.05

    0.10

    0.2

    0.3

    -0.05

    0.00

    0.0

    0.1

    0 2000 4000 6000 8000 10000Time [hours]

    0 2000 4000 6000 8000 10000Time [hours]

  • ECI3-Solute 2C 3 So uteKd =0.001mg/cm3

    EnestecasosedisminuyeelvalordeKd (materialdel/suelomenosreactivoKd_inicial/5).Losotros

    parmetrosycondicionesdeflujosoninvariantes.

  • 0.008

    Observation Nodes: Concentration

    0

    Profile Information: Concentration

    0 004

    0.005

    0.006

    0.007

    -40

    -20

    0

    0.001

    0.002

    0.003

    0.004

    -80

    -60

    40

    0.0000 2500 5000 7500 10000

    Time [hours]-100

    -80

    0.000 0.002 0.004 0.006 0.008 0.010Conc [mg/cm3]Conc [mg/cm3]

    0

    Profile Inform.: Immobile Conc.

    -40

    -20

    -80

    -60

    -1000.000 0.002 0.004 0.006 0.008 0.010

    Conc [mg/cm3]

  • 0.30

    0.35

    Surface Solute Flux

    0.35

    0.40

    Cumulative Surface Solute Flux

    0.20

    Cumul. Non-Equil. Mass Transfer

    0.10

    0.15

    0.20

    0.25

    0.15

    0.20

    0.25

    0.30

    0.10

    0.15

    -0.05

    0.00

    0.05

    0 2000 4000 6000 8000 10000Time [hours]

    0.00

    0.05

    0.10

    0 2000 4000 6000 8000 10000Time [hours]

    0.00

    0.05

    0 2000 4000 6000 8000 10000Time [hours]Time [hours] Time [hours]

    0 40

    All Cumulative Solute Fluxes

    0.30

    0.35

    0.40

    0 15

    0.20

    0.25

    0.05

    0.10

    0.15

    0.000 2000 4000 6000 8000 10000

    Time [hours]

  • ECI3-Solute 3C 3 So ute 3Kd =0.0001mg/cm3

    EnestecasesedisminuyeelvalordeKd (materialdel/suelomuypocoreactivoKd_inicial/50).Losotros

    parmetrosycondicionesdeflujosoninvariantes.

  • 0.010

    Observation Nodes: Concentration

    0

    Profile Information: Concentration

    0.006

    0.008

    -40

    -20

    0.002

    0.004

    80

    -60

    0.0000 5000 10000

    Time [hours]

    -100

    -80

    0.000 0.002 0.004 0.006 0.008 0.010Time [hours] Conc [mg/cm3]

    0

    Profile Inform.: Immobile Conc.

    -40

    -20

    -80

    -60

    -1000.000 0.002 0.004 0.006 0.008 0.010

    Conc [mg/cm3]

  • Surface Solute Flux Cumulative Surface Solute Flux BottomConcentration

    0 06

    0.08

    0.10

    Surface Solute Flux

    0.25

    0.30

    0.35

    0.004

    0.005

    0.006

    Bottom Concentration

    0.00

    0.02

    0.04

    0.06

    0 05

    0.10

    0.15

    0.20

    0.001

    0.002

    0.003

    0.004

    -0.02

    0.00

    0 2000 4000 6000 8000 10000Time [hours]

    0.00

    0.05

    0 2000 4000 6000 8000 10000Time [hours]

    0.0000 2000 4000 6000 8000 10000

    Time [hours]

    0.16

    Cumul. Non-Equil. Mass Transfer

    0.000

    Cumulative Bottom Solute Flux

    0.000000

    Bottom Solute Flux

    0.08

    0.10

    0.12

    0.14

    -0.015

    -0.010

    -0.005

    -0.000020

    -0.000015

    -0.000010

    -0.000005

    0.00

    0.02

    0.04

    0.06

    0 2000 4000 6000 8000 10000-0.030

    -0.025

    -0.020

    0 2000 4000 6000 8000 10000-0.000035

    -0.000030

    -0.000025

    0 2000 4000 6000 8000 10000 0 2000 4000 6000 8000 10000Time [hours]

    0 2000 4000 6000 8000 10000Time [hours]

    0 2000 4000 6000 8000 10000Time [hours]

  • All Cumulative Solute Fluxes

    0.35

    0 25

    0.30

    0.20

    0.25

    0.10

    0.15

    0 00

    0.05

    -0.05

    0.00

    0 2000 4000 6000 8000 100000 2000 4000 6000 8000 10000Time [hours]