Advances in Ground-based Gravity for Hydrologic Studies

Transcription

1 Advances in Ground-based Gravity for Hydrologic Studies L. Longuevergne (1,2), C.R. Wilson 1, B.R. Scanlon 2 1 Department of Geological Sciences, UT Austin 2 Bureau of Economic Geology, UT Austin

2 What is Hydro-geodesy? Study of the continental part of the hydrological cycle Study of the Earth s deformations and the Earth s gravity field aet g Slope Water mass redistribution Gravity variations Load deformation

3 What is Hydro-geodesy? Conventional units 1 nm.s -2 (i.e. 0.1 Gal) gravity change is the attraction of an infinite 2.4 mm thickness water layer. Significant contribution of hydrology on gravity time series Among others: Lambert and Beaumont, JGR, 1977 Mäkinen and Tattari, IAG, 1990 Bower and Courtier, PEPI, 1998 Hinderer & Crossley, Surv. Geophys, 2000 Kroner et al., J. geodyn, 2004 Van Camp et al., JGR, 2006 Development of geodetic methods for hydrological studies Among others: Pool and Eychaner, Groundwater, 1995 Hasan et al., WRR, 2008 Longuevergne, PhD, 2008 Naujok, PhD, 2008 Krause et al., JoH, 2009 Creutzfeld et al., HESS, 2010

4 Gravity variations [ nm.s -2 ] Spraying test Spraying test in MOXA observatory (Germany) Spraying Natural rainfall 12 nm.s -2 Drainage Gravi Time From Kroner et al., J. Geodyn., 2004 Measurable effect 4cm fullwaterlayer 12 nms g

5 Another test Comparison to soil moisture measurements From Longuevergne et al., J. Geodyn., 2009

6 Outline Physical processes and spatial scales Instrument types Observation strategies and examples Application of gravity tools for hydrological studies

7 Physical processes 1. Newtonian Attraction Horizontal

8 Physical processes 1. Newtonian Attraction Information on water mass Horizontal Water mass

9 Physical processes 1. Newtonian Attraction Information on water mass 2. Elastic deformation Horizontal Surface

10 Physical processes 1. Newtonian Attraction Information on water mass 2. Elastic deformation Information on the pressure applied on the crust (water mass) (Farrell, RGSP, 1977) Water mass Horizontal Surface

11 Physical processes 1. Newtonian Attraction Information on water mass 2. Elastic deformation Information on the pressure applied on the crust (water mass) (Farrell, RGSP, 1977) Water mass Horizontal Surface Phenomenon Eq. water layer Vertical deformation Tide (Britany) 10 m 10 cm (Llubes et al., 2007) Amazon basin water 0.3 m 3 cm (VanDam et al., 2001)

12 Spatial scales Measurement methods Gravity variations Tilt variations Vertical deformation Main spatial scales of sensitivity Ground gravimeter Tiltmeter GPS Water mass Surface Llubes et al., J. geodyn, 2004 GRACE satellite 100 m 1 km 10km 100 km 1000 km Investigated Spatial scales

13 Spatial scales Measurement methods Gravity variations Tilt variations Vertical deformation Main spatial scales of sensitivity Ground gravimeter Tiltmeter GPS Water mass Surface Llubes et al., J. geodyn, 2004 GRACE satellite 100 m 1 km 10km 100 km 1000 km Investigated Spatial scales

14 Recorded contributions Local mass sensitivity g Surface

15 Sensitivity Recorded contributions Local mass sensitivity g space

16 Sensitivity Recorded contributions Local mass sensitivity g Surface h 66% at 6h 70% at 6h Relative height h between gravimeter & target 70% of sensitivity on a 6xh diameter disc space

17 Sensitivity Recorded contributions Local mass sensitivity g Surface h 66% at 6h 70% at 6h Relative height h between gravimeter & target 70% of sensitivity on a 6xh diameter disc space Sensitivity to global scales, amplitude ~ 30 nm.s -2, seasonal Especially linked to the elastic deformation of the Earth Typically corrected using LSM simulations Boy & Hinderer, J. geodyn, 2006 Wziontek et al., J. geodyn, 2009

18 Sensitivity Recorded contributions Local mass sensitivity g Surface h 66% at 6h 70% at 6h Relative height h between gravimeter & target 70% of sensitivity on a 6xh diameter disc space Sensitivity to global scales, amplitude ~ 30 nm.s -2, seasonal Especially linked to the elastic deformation of the Earth Typically corrected using LSM simulations Boy & Hinderer, J. geodyn, 2006 Wziontek et al., J. geodyn, 2009 Other contributions: tides (Wenzel, BIM 1997) atmospheric effect (Boy & Lyard, GJI, 2008) vertical displacement 1 cm ~ 30 nm.s -2

19 Outline Physical processes and spatial scales Instrument types Observation strategies and examples Application of gravity tools for hydrological studies

20 Hydrology Instrument improvement cm water Relative & absolute gravimeters Year Adapted from Torge, 1989

21 Hydrology Instrument improvement cm water Relative & absolute gravimeters 2000 Year 0.25 Superconducting gravimeter Adapted from Torge, 1989

22 Which gravimeter for hydrological studies? Micro g Lacoste Gphone ~ 100 k$ CG-5 ~ 100 k$ ZLS Burris ~ 50 k$ Instrument Resolution (eq. full water layer) Field gravimeter (relative) 25 mm Drifting instruments Rq Portable instrument, quick measure Ideal for network repetition

23 Which gravimeter for hydrological studies? Micro g Lacoste FG5 ~ 400 k$ A-10 ~ 400 k$ Instrument Resolution (eq. full water layer) Rq Field gravimeter (relative) 25 mm Drifting instruments Absolute gravimeter (free fall) 25 mm Absolute Transportable instrument Provide reference value

24 Which gravimeter for hydrological studies? GWR OSG ~ 450 k$ igrav ~ 250 k$ Instrument Resolution (eq. full water layer) Rq Field gravimeter (relative) 25 mm Drifting instruments Absolute gravimeter (free fall) 25 mm Absolute Superconducting (relative) < 2 mm Very low drift Accurate and stable instrument Ideal for continuous monitoring

25 Outline Physical processes and spatial scales Instrument types Observation strategies and examples Application of gravity tools for hydrological studies

26 Vertical integration Observation strategy Continuous observation : monitoring vertically integrated mass variations

27 Vertical integration Observation strategy Continuous observation : monitoring vertically integrated mass variations Unconfined aquifer g geom. Sy. h g h Gravity variation Well level variation Sy ~10 1 geom Sy Geometric factor Porosité de drainage

28 GRAVITY CHANGE, MICROGAL DEPTH TO WATER, METERS Vertical integration Observation strategy Example in Tucson, Arizona R = 0.94 S.Y = 0.27 Sep-97 Sep-98 Sep-99 Sep-00 Sep-01 Sep Good correlation unconfined aquifer From Pool, Geophysics, 2008

29 Vertical integration Observation strategy Continuous observation : monitoring vertically integrated mass variations Soil moisture Unconfined aquifer

30 GRAVITY CHANGE, MICROGAL DEPTH TO WATER, METERS Vertical integration Observation strategy Example in Tucson, Arizona Gravity Depth to Water Sep-97 Sep-98 Sep-99 Sep-00 Sep-01 Sep Significant vadose zone contribution From Pool, Geophysics, 2008

31 Vertical integration Observation strategy Continuous observation : monitoring vertically integrated mass variations Soil moisture Unconfined aquifer Confined aquifer

32 Vertical integration Observation strategy g geom. S. h g h Gravity variation Well level variation S ~10 3 geom S Geometric factor Storativity

33 GRAVITY CHANGE, MICROGAL DEPTH TO WATER, METERS Vertical integration Observation strategy g geom. S. h g h Gravity variation Well level variation S ~10 3 geom S Geometric factor Storativity Example in Tucson, Arizona Gravity Depth to Water Sep-97 Sep-98 Sep-99 Sep-00 Sep-01 Sep From Pool, Geophysics, 2008

34 Network repetition Observation strategy - No continental effect in spatial differences, focus on local effects - Investigate lateral variations, localize heterogeneities Soil Moisture Groundwater

35 Network repetition Observation strategy - No continental effect in spatial differences, focus on local effects - Investigate lateral variations, localize heterogeneities Soil Moisture Groundwater Example in MOXA Observatory (Germany) Geometry validation of a spatially-distributed model From Naujok, PhD, 2008, Naujok et al., J.Geod., 2008

36 Network repetition Observation strategy - No continental effect in spatial differences, focus on local effects - Investigate lateral variations, localize heterogeneities - Vertical gradient: g 2 times water storage in the layer in between Soil Moisture Groundwater Example in Durzon karst (France) Estimation of epikarst storage From Jacob et al., GJI., 2009

37 Multi-instrumental investigation Observation strategy 1. Precise local measurements 2. Inference at larger scale 3. Localization of water mass 4. Investigation of lateral & vertical heterogeneity Soil Moisture Groundwater

38 Multi-instrumental investigation Observation strategy 1. Precise local measurements 2. Inference at larger scale 3. Localization of water mass 4. Investigation of lateral & vertical heterogeneity Soil Moisture Groundwater Ongoing projects (among others) International project gathering several teams operating superconducting gravimeters GGP network (Global Geodynamic Project) Crossley et Hinderer, J. Geodyn, 2009

39 Multi-instrumental investigation Observation strategy 1. Precise local measurements 2. Inference at larger scale 3. Localization of water mass 4. Investigation of lateral & vertical heterogeneity Soil Moisture Groundwater Ongoing projects (among others) Ground and space gravity In AMMA / CATCH area GHYRAF project (Gravity and Hydrology in Africa) Hinderer et al. J. Geodyn, 2009 GGP network (Global Geodynamic Project) Crossley et Hinderer, J. Geodyn, 2009

40 Outline Physical processes and spatial scales Instrument types Observation strategies and examples Application of gravity tools for hydrological studies

41 Pumping test Hydrological applications Simulation of a homogeneous unconfined aquifer Q=5 m 3 /day over 7 days From Damiata et al., JoH; 2006

42 Pumping test Hydrological applications Simulation of a homogeneous unconfined aquifer Q=5 m 3 /day over 7 days T = m 2 /s Sy=0.25 T = m 2 /s Sy=0.15 T = m 2 /s Sy=0.25 T = m 2 /s Sy=0.35 T = m 2 /s Sy=0.25 From Damiata et al., JoH; 2006

43 Pumping test Hydrological applications T = m 2 /s Sy=0.15 T = m 2 /s Sy=0.25 T = m 2 /s Sy=0.35 Simulation of a homogeneous unconfined aquifer Q=5 m 3 /day over 7 days T = m 2 /s Sy=0.25 T = m 2 /s Sy=0.15 T = m 2 /s Sy=0.25 T = m 2 /s Sy= Complementary information - Heterogeneity investigation from surface measurements - Artificial recharge monitoring T = m 2 /s Sy=0.25 From Damiata et al., JoH; 2006

44 Model validation Hydrological applications Mass balance equation S Total water storage (catchment) d ( S) dt P aet Q s Q sub P aet Q Precipitations Actual evapotranspiration Surface runoff Q sub Subsurface runoff After time integration g ( t) S P Q s dt aet Q sub dt Measurable fluxes Fluxes to be modeled Any systematic error generate a drift-like signal Gravity brings a additional constrain on fluxes realism

45 Model validation Hydrological applications Example in Strasbourg observatory (France) SG CLSM model Residuals

46 Model validation Hydrological applications Example in Strasbourg observatory (France) SG CLSM model Residuals Main discrepancy linked to 2003 European heat wave Gravity data highlight a systematic error in ET CLSM does not stress vegetation enough in dry conditions

47 Adapting SGs for field Hydrological applications Idea: Keep SG characteristics for field studies Example in Austin, Texas From Wilson et al. IAG, 2009

48 Adapting SGs for field Hydrological applications Idea: Keep SG characteristics for field studies Example in Austin, Texas Austin, TX Worst drought in 50 years Results: original sensitivity for periods > 1 day Importance of pillar From Wilson et al. IAG, 2009

49 Adapting SGs for field Hydrological applications Idea: Keep SG characteristics for field studies Example in Austin, Texas Austin, TX Worst drought in 50 years Results: original sensitivity for periods > 1 day Importance of pillar From Wilson et al. IAG, 2009

50 Adapting SGs for field Hydrological applications Idea: Keep SG characteristics for field studies Next generation of SGs igrav, GWR ~ 50 kg Low energy consumption ( ~1.5 kw )

51 Conclusions Interest of gravity instruments Sensitive to water storage variations complementary information Integrative information (~ 100 m) Non-destructive & well calibrated information Multi-instrument study Precise local measurements (~ few mm) Inference at larger scale with network repetition Water mass localization in active hydrological systems Lateral & vertical heterogeneity investigation

52 Thank you for attention Thank you for attention