VieVS Analysis of a single session

Transcription

1 VieVS Analysis of a single session Hana Krásná, J. Böhm, M. Madzak, L. Plank 1, K. Teke 2, A. Hellerschmied, A. Hofmeister 1 University of Tasmania, Australia, 2 Hacettepe University, Turkey

2 VieVS Structure 2

3 3

4 Vie_setup GUI for Vienna VLBI Software Easy to use (windows style) Comes with plotting tool 4

5 Without previous knowledge: Visualize estimated parameters Visualization for analysts (clock breaks, exclude stations ) Plot world map, baseline length repeatability, covariance matrix Output EOP, baseline lengths new 5

6 Matlab interface.m file Code.fig file Objects Starts by running.m file 6

7 Minimum for VieVS Select session(s) Click Save + Run 7

8 More options A priori models (vie_mod) Estimation options (vie_lsm) Scheduling (vie_sched) Simulation (vie_sim) Global solution (vie_glob) 8

9 Plotting In menu Plotting Residuals Parameters Session information Output EOP, baselines new 9

10 Plot residuals First (only clock+zwd) Clock breaks Main solution Baseline wise Station wise Source wise 10

11 Plot residuals Outliers are shown in red Show station/source numbers new + Change reference clock + Clock breaks adding + Outlier selection 11

12 Plot parameters One/all session(s) in folder Parameters per station Comparison (up to three) possible Save values as text Print to any format 12

13 Plot session information Session network* Baseline length repeatabilities* Correlation matrix * Up to four sessions 13

14 Plot session information Session network* Baseline length repeatabilities* Correlation matrix * Up to four sessions 14

15 Plot session information Session network* Baseline length repeatabilities* Correlation matrix * Up to four sessions 15

16 Save and load GUI state Get GUI options from a parameter file Same parameterization Defaults Current before very latest action Save process_list 16

17 Program structure >> vievs opens current GUI version (2.2) >> vievs( xx ) opens GUI of version x.x >> vievs( batch ) runs batch version of vievs GUI Batch = Prepare processing = Do processing 17

18 vie_setup Save runp ( Last user action in GUI ) Creates input protocol Saves runp.mat Saves process_list.mat Saves parameter.mat file(s) Saves sched, sim and glob parameters 18

19 Input protocol input_protocol.txt in /WORK/ Saves GUI options (models, ) in textfile For user information (not used in processing!) 19

20 Runp.mat In /WORK/ Is overwritten (exists once) Gives: paths (LEVELx subfolders) 1 0 for all modules 20

21 Process_list.mat In /WORK/ Is overwritten (exists once) Gives VLBI experiments (one line = one session) 21

22 Vie_batch (1) Does (batch) processing E.g. vie_batch2_2.m Independent fromgui Does processing (vie_init, vie_mod, vie_lsm, ) 22

23 Vie_batch (2) Input Processing Output 23

24 24

25 What does VIE_INIT do? Reads observations from the NGS file Reads station coordinates and velocities from the superstation file Read source coordinates from supersource file Possible to: Remove outliers (specified in an outlier file) Exclude stations, sources, baselines (specified in OPT file) Introduce an elevation cut off angle 25

26 Superstation file Can be created/updated using the GUI Contain station coordinates, velocities, and additional antenna info (mount, axis offset, eccentricity, etc.) ITRF2005, ITRF2008, VTRF2008, VieTRF10a, VievsTrf If a station is not found in the chosen TRF, the coordinates from the backup VievsTrf are used.txt files in TRF directory: User defined TRF in ASCII format 26

27 Supersource file Can be created/updated using the GUI ICRF_Ext2, ICRF2, VieCRF10a, and VievsCRF Contain source coordinates If a source is not found in the chosen catalogue, the VievsCrf coordinates are used. 27

28 Outlier file Contains list of outliers for the session Created in VIE_LSM Outliers are removed in VIE_INIT. To detect and remove outliers you need to run VieVS twice: In the first run outliers are detected in VIE_LSM and saved it in an outlier file In the second run this file run is used in VIE_INIT for removing the outliers 28

29 OPT file Can be used to: Specify reference clock Introduce clock breaks Exclude stations, sources, and baselines Not use the cable cal of a station ASCII files in sub folders of DATA/OPT/, e.g. DATA/OPT/DEFAULT/2009/09JUN23XH.OPT 29

30 OPT file CLOCK REFERENCE: ONSALA60 CLOCK BREAKS: 3 WETTZELL WETTZELL KOKEE STATIONS TO BE EXCLUDED: 1 ONSALA85 SOURCES TO BE EXCLUDED: BASELINES TO BE EXCLUDED: 2 HOBART26 HOBART12 MATERA MEDICINA NO CABLE CAL: 1 FORTLEZA # This is a comment ONSALA60 is the reference clock Three clock breaks Exclude all observations including ONSALA85 Exclude all observation to sources , , and Exclude all observations of baselines HABART26-HOBART12 and MATERA-MEDICINA Do not use cable cal corrections for FORTLEZA 30

31 OPT file new STATIONS TO BE EXCLUDED: 1 HOBART a certain downtime for stations rather than excluding them for the whole session can be defined observations including HOBART12 between 01UTC 07UTC on Feb 8, 2013 will be excluded useful also for simulations 31

32 32

33 Command Window output These stations have no TRF coordinates 33

34 Command Window output Number of stations, sources, and scans 34

35 Command Window output Names of the stations 35

36 36

37 Theory computed delay τ τ = - r r b k c Vie_mod computed delay observed delay from NGS card file, corrected for ionosphere o c adjustment (lsm) 37

38 Theory models + EOP (oceanic & gravitational high frequency terms, nutation corrections) + solid Earth tides + troposphere delay (VMF / GMF / NMF) + gradients + ocean loading + thermal antenna deformation + atmosphere loading + hydrology loading partial derivatives τ dvar Vie_mod computed delay τ observed delay from NGS card file, corrected for ionosphere o c adjustment (lsm) 38

39 Vie_MOD usage TRF: stored at:../trf/ station coordinates (TRS) reference epoch station velocities Other: self created TRF../TRF/*.txt specified format 8 char. stat. name, column mode 39

40 Vie_MOD usage CRF: stored at:../crf right ascension ra declination de Other: self created CRF not finished yet 40

41 Vie_MOD usage Ephemerides: stored at:../ephem/jpl_421.mat Earth, Sun, Moon, planets remembered for 2nd run (jpl_421_05sep12xa_n004) not critical for normal VLBI default: gravitational bending (GRT) off all planets is taken into account. 41

42 Vie_MOD usage Pressure and temperature: Met data from NGS file GPT2 as backup External.trp files 42

43 Vie_MOD usage Mapping function: for atmospheric corrections on the delay VM1:../VM1/yyyy_VM1.mat regular update! GMF: applied when no VM1 available; no files necessary GPT2 43

44 Vie_MOD usage A priori troposphere gradients: 44

45 Vie_MOD usage Ionosphere correction: subtracted from the observed delay External ionosphere correction: must be created in advance stored at../ion/maps 45

46 Vie_MOD usage Station corrections: tidal & non tidal Thermdef: [A. Nothnagel, J. of Geodesy, 2008] 46

47 Vie_MOD usage Earth Orientation Parameters: c04 08 (daily) series from the IERS webpage;../eop/c04_08_1962_now.txt; updates might be necessary; predefined EOP:../EOP/*.txt; watch format finals: necessary for recent sessions dx, dy: include/exclude from c04 table; mainly contains FCN 47

48 Vie_MOD usage Earth Orientation Parameters: high frequency EOP ocean tides: interpf; table 8.2 (xp,yp) & 8.3 (ut) diurnal & semidiurnal libration: table 5.1a (xp,yp) & 5.1b (ut) IERS Conv 2010 Precession/Nutation Model 48

49 Vie_MOD usage Interpolation of EOP: lagrange: window of 4 data points linear: interpolation between the midnight before and after observation time. For a session from 18:00 to 18:00 this means, that there are 2 a priori lines and a break at midnight. tidal UT variations: removed before interpolation, re added 49

50 Vie_MOD usage Quality code limit: Accept observations < or = limit Good quality observations have quality code 0 Cut off Elevation angle: all observations below this local elevation [, degrees] are skipped 50

51 Vie_MOD usage For more information on the formalism, check the code (vie_mod.m) or DOC/vie_mod.pdf Partial derivatives are also set up Results are stored in DATA/LEVEL1/sesname_antenna _parameter _scan _sources DOC/structures.pdf 51

52 52

53 VIE_LSM_V22 VIE_LSM is a module of VieVS, which estimates geodetic parameters with least squares adjustment from VLBI observations. All the parameters are estimated as continuous piece wise linear offsets (CPWLO) in sub daily and daily temporal resolution. 53

54 VIE_LSM_V22 continuous piece wise linear offsets (CPWLO) t t x = x + ( x x ) i t2 t1 54

55 VIE_LSM_V22 Estimated parameters: Clocks (offset (cm), rate (cm/day), quadratic term (cm/day 2 ), CPWLO (cm)), Zenith wet delays (cm) as CPWLO, Troposphere gradients (cm) as CPWLO, EOP (mas and ms) as CPWLO, Antenna coordinates in TRF (cm) as one offset per session or as CPWLO, Source coordinates in CRF (declinations in mas and right ascensions in ms) as CPWLO. 55

56 Parametrization for Least Squares 56

57 Parametrization for relative clock errors Relative constraints on the clock offsets. Observation equations are added to the design matrix telling that the difference between two piecewise linear clock offsets is zero ± a certain standard deviation σ. (Mainly important to bridge gaps without observations to avoid singularity of the normal equation system.) 57

58 Troposphere delay H trop [ ] Δ τ = N s + N s ds 6 trop 10 h( ) w( ) 0 [ ] Δ τ ( αε, ) = ZHD m ( ε) + ZWD m ( ε) + m ( ε)cot( ε) G cos( α) + G sin( α) trop h w w n e reduced from observations a priori to the adjustment estimated estimated estimated 58

59 Troposphere delay zwd gradients 59

60 Parametrization for EOP To estimate one constant value per session strong relative constraints of 1e 4 m(a)s/day have to be set Example: The session is from 18 UT to 18 UT. Then, three piecewise linear offsets are set up for each EOP. (midnight before the session, during the session, and the midnight after the session) The strong constraints take care that all three estimates per session are the same. 60

61 VIE_LSM scan wise update 61

62 Scan wise update of normal equation system 62

63 VIE_LSM VIE_LSM estimates parameters as CPWLO from single session adjustment VIE_LSM corrects clock breaks and detects outlier observations VIE_LSM provides SINEX input and datum free normal equations for global solutions (VIE_GLOB) 63

64 Thank you for your attention!