The conversion of heights of the benchmarks of the detailed vertical reference network into the PL-EVRF2007-NH frame

The problem of calculating height of the benchmarks of the vertical reference network should be examined in the context of features of a given height frame and the type of the vertical reference network.

Pursuant to the Ordinance of the National Spatial Reference System, heights of the reference network points should be expressed in the Kronstad Vertical Reference Frame 1986 for Poland (PL-KRON86-NH) and the European Vertical Reference Frame 2007 for Poland (PL-EVRF2007-NH) (Regulation, 2012b). However, as of 1 January 2020, the PL-EVRF2007-NH frame with the reference level of the North Sea in Amsterdam (Somla, 2018) will be binding. Therefore, it is necessary to develop a new vertical reference frame. Due to the division of the surveying resource existing in Poland, including geodetic control networks into main (basic networks) and district (detailed networks) control networks, the task of calculating the height is performed in two stages, according to the classification of the given network (Regulation, 2012a).

During the first stage, heights of the basic vertical reference network benchmarks in the PL-EVRF2007-NH (Graszka et al., 2016; Kadaj, 2018) frame were calculated. The Head Office of Geodesy and Cartography provides data in both above-mentioned height frames, as well as a model of height conversions between these frames (PL-KRON86-NH and PL-EVRF2007-NH), developed for the area of the country (Head Office of Geodesy and Cartography, 2015). This article focuses on the second stage, that is, the conversion of height of the detailed reference network benchmarks. Although these works are already carried out in many districts, due to the diversified state of the detailed reference network, the use of different height frames and different financial means of the districts, the process of introducing the PL-EVRF2007-NH frame at the district level seems to be more complex. The article attempts to analyse selected aspects of height conversion in the detailed reference network. The proposed method based on height transformation in conditions of good saturation of the area with adjustment points may give results similar to those of the method of rigorous adjustment of levelling observations. The analysis is illustrated by an example of such conversion carried out for one of the districts.

The regulations concerning the classification of reference networks are contained in the Regulation on geodetic, gravimetric and magnetic control networks (Regulation, 2012a). The vertical reference network can be classified into basic fundamental, basic and detailed network. Such division corresponds to the accuracy of the network realised by classes 1, 2 and 3 accordingly. Surveyor General of Poland is responsible for the basic reference network, whereas the starost (head) of a district is responsible for the detailed reference network. The basic classification criterion is ‘the accuracy of determining the size (. . .) defined by the value of the average error of the determined size’ (Regulation, 2012a). For the basic vertical reference network, this value should not exceed 1.5 mm/km in spirit levelling. The results of the levelling measurement should be prepared in accordance with the comparative and thermal correction, the tidal correction and the levelling correction (normal) and permanent tide correction. Before aligning the observations, a detailed evaluation of the network should be carried out in order to calculate the values of average measurement errors (levelling line, random, systematic error and so on) (Regulation, 2012a). The calculation of benchmarks heights is, therefore, performed only by precision spirit levelling. Other methods of levelling are not allowed.

For the detailed reference network, the accuracy requirements are lower. First of all, apart from spirited levelling, satellite levelling is also allowed, which may be a convenient solution to the problem of levelling (e.g. in the case of terrain obstacles) (Banasik and Bujakowski, 2017). The average measurement error in this network should not exceed 4 mm. The above-mentioned corrections, used in precision levelling, are not introduced to observations in the detailed reference network.

The above synthesis of regulations allows us to claim that the basic (national) reference network is established by geodetic methods (e.g. taking into account the influence of geodynamic phenomena), whereas the detailed (district) reference network is established by land surveying methods.

The introduction of the new state elevation frames was the result of successive measurement campaigns of the basic vertical reference network in Poland. Over the years, in addition to the development of measurement technology, the way in which the tidal phenomenon was included was also subject to change. The mean-tide system was replaced by the tide-free system, and then by the EVRS2007 system, which by definition is zero-tide (Mäkinen, 2008). In the 1960s the error of the Potsdam system was observed (Barlik and Pachuta, 2007). The models of vertical movements of the Earth’s crust were also more detailed and the constant value of the whole country’s settlements of about -2 mm/year (Kowalczyk, 2006) was abandoned. The monitoring of the state precise satellite positioning network stations (called ASG-EUPOS) carried out since 2008 has provided information on the simultaneous existence of settlement and uplift areas (Kontny and Bogusz, 2012; Maciuk, 2016; Maciuk and Szombara, 2018). The result of the above changes is the lack of a constant conversion between the PL-KRON86-NH and PL-EVRF2007-NH elevation frames in Poland. According to the model made available by the Head Office of Geodesy and Cartography (HOGC, pl. GUGiK), the difference is within the range from +12.8 to +20.6 cm, with an average value of +16.6 cm (Borowski et al., 2017).

The synthesis of the heights frames is presented in Table 1. The previous state and local frames are not included in it. These frames were simultaneously used until recently in some cities (e.g. in Warsaw or Ełk): the former Amsterdam (1926–1937), Kronstadt’60 and Kronstadt’86.

One of the first works on the conversion of height between the Kronstadt’60 and Kronstadt’86 frames concerned the detailed reference network of Kraków district (Banasik et al., 2012). For this purpose, the authors used a polynomial transformation function, developed with the use of adjustment benchmarks coming from the basic reference network and with heights in both frames. The flat coordinates of the benchmarks necessary for transformation were obtained by digitalisation of the relevant cartographic materials.

Another study in this field concerned the detailed reference network in the Wrocław area (Osada et al., 2017; Osada and Gralak, 2017). The authors of this article applied a rigorous adjustment of archival observations from the levelling network, in relation to the benchmarks with known heights in the new frame. The whole procedure consisted of four stages, with each being adjusted in a different reference configuration. Most of the work preceding the calculation was related to the completion and preparation of observation data for adjustment.

In 2018, the first information on the conversion of the height of detailed reference networks to the PL-EVRF2007-NH frame (Graszka, 2018) appeared. Then, on the website of the GUGiK, the proposed methodology for height conversion (Head Office of Geodesy and Cartography, 2018) was published. The trainings in this field among district surveyors were also conducted (Królikowski, 2019).

The method recommended by GUGiK to calculate the height of the detailed reference network benchmarks is a rigorous adjustment of observations in the levelling network (Graszka, 2018; Head Office of Geodesy and Cartography, 2018; Somla, 2018).

Transformation is only permitted for the measurement reference network and situational details.

There may be doubts in case 2 and 3 (Table 2). At least some of the 380 districts may have difficulties in financing this project. The reason is the higher priority of the horizontal reference network on which cadastral data are based, as well as the establishment and maintenance of district databases. The introduction of a new heights frame enforces investments in the vertical reference network, and the above recommendations of the GUGiK increase its value. The approximate costs of case no. 2 given in Table 2 seem to be underestimated. For an ‘average’ district with 1000 points of vertical reference network, the financial outlays are as follows:

Collection of materials, digitalisation and preparation of data for adjustment – minimum 20,000 PLN

Thus, the real cost of elaboration may be even twice as high as estimated (case 2, Table 2). To confirm this, one can give an example of the reference network in Szczecinek district (of more than 900 points), where the cost of elaboration exceeded even the amount of 100 PLN per point (Fabisiak, 2019). Another example is the detailed reference network in Białystok, consisting of about 5400 points (Bialystok City Council, 2019). For this study, carried out without field-works and with a control measurement only for 1.5% of points, the cost exceeded the GUGiK’s estimates and amounted to over 60 PLN per point. Despite these significant costs, it is not certain whether the benchmarks of the reference network have not changed their location since the time of their measurement or have been destroyed (e.g. as a result of thermal modernisation of buildings). In order to eliminate the above factors, a full modernisation of the vertical reference network is necessary, which is a much more expensive operation (case 3, Table 2).

There is some misunderstanding about the statement concerning the preparation for the rigorous adjustment of the levelling network, contained in the recommendations (Chapter II, point 2, subitem 5 (Head Office of Geodesy and Cartography, 2018)):

Before the adjustment of the detailed vertical reference network, it is necessary to reduce archival observations to one epoch and one reference frame, which requires knowledge about methods and accuracy of the observations and methods of their reduction⁴⁾

Therefore, a question arises whether, in connection with the above, archival levelling observations in a detailed vertical reference network should be provided with appropriate corrections (e.g. lunosolar (tidal) and normal corrections), if they are to be used for re-alignment in the PL-EVRF2007-NH frame. This question can be answered after analysing the provisions of the ordinance on the vertical reference network, in the part concerning the detailed vertical reference network (Regulation, 2012a). However, it is worth using this opportunity to provide the justification for these provisions.

Observations in any levelling network (i.e. raw levelling measurement) results include, among others, the temporary effect of tidal phenomenon during the measurement. Thus, the levelling network, measured at various time intervals, in which the lunosolar correction was not taken into account, will be similar to the mean-tide model used in the Kronstadt’60 frame. The introduction of the lunosolar correction to the measured elevations ‘transfers’ those elevations to the ‘tide-free’ model that was used in the Kronstadt’86 frame, or to the ‘zero-tide’ model (EVRF2007 frame). The biggest differences between the above-mentioned tidal models occur between the mean- and zero-tide models, or, in the case of our country, between the Kronstadt’60 and PL-EVRF2007-NH frames (Mäkinen, 2008). Changing the tidal model used in these systems results in a difference in height of a few centimetres (Fig. 1). By converting this difference into a distance we get a value of a fraction of a millimetre per kilometre. Levelling (normal) correction will result in similarly small values. Therefore, we do not have to consider using tidal or normal correction in the detailed vertical reference network. In the use of archival measurement data, first of all, their completeness is important, as it enables the alignment of observations within the levelling network.

Figure 1

There is no doubt that the best solution to the problem of calculating the height of the benchmarks in detailed vertical reference network is its modernisation. Only this way a reference network will be fully functional, and its benchmarks will have current heights in the new elevation frame. For districts that fall under case 2 or 3 (Table 2), until the modernisation of the vertical reference network, a bridge solution would be useful, ensuring sufficient data quality and at the same time low cost of development. The proposal for such solution consists of the following stages:

A preliminary check whether it is possible to perform the transformation – an analysis of data (e.g. geoportal.gov.pl data (Head Office of Geodesy and Cartography, 2019)). A necessary condition is a good saturation of the area with adjustment points that is, having heights in both frames.

The above proposal assumes resignation from the rigorous adjustment and the need to perform field-work. In this manner, the most expensive elements of the development are minimised. The financial outlays for the ‘average’ district should amount to:

Collection of materials, digitisation and preparation of data – minimum 5000 PLN.

Thus, the cost of conversion of the heights should amount to approximately 25 PLN per point of the vertical reference network for office works only. Introduction of small measurement works doubles this number. The total amount of such conversion should not exceed 30,000 EUR (~130,000 PLN). The procedure for selection of the contractor may therefore be carried out on the basis of a request for quotation. In comparison with the tender procedure, this is a procedure with a shorter formal timeframe before the actual works are commenced. For districts that ‘postpone’ the conversion, the use of the tender procedure may result in exceeding of the deadline imposed by the regulation for the introduction of the new frame.

4.1

The description of the procedure – an example of Kętrzyn district

A detailed description of the proposed procedure is illustrated with an example of Kętrzyn district

License for scientific and educational purposes for the data of the reference network no. GKN-O.6642.1.1992.2017_2808_N issued by the Starost of Kętrzyn.

, where the Kronstadt’60 frame was used so far. The height conversion was done before the publication of the GUGiK recommendations.

In the method proposed above, rigorous adjustment was omitted. This significantly reduces the cost of data acquisition. Determining the coordinates of points through digitalisation is cheaper than scanning and laborious processing of the lists of elevations, especially in the case when there are several operators of the vertical reference network in a given district. The second factor that helps to reduce the cost is resignation from the field-works for all points. The example given in Section 3 indicates that this part of the work represents almost 50% of the order value. It seems that the field review of all points is an unnecessary expenditure, with little added value.

Before starting further works, it is recommended to supplement the geodetic network database with missing elements, such as non-digitalised topographic descriptions, corrections of entries in the database

There may be errors in geodetic database records, for example in the field of ‘frame height 1’ there are given heights from the Kronstadt’60 and PL-KRON86-NH frames.

, or conversion of coordinates (x, y) to the single, national coordinate system.

The initial control element of the reference network is the creation of profiles of differences of levelling lines (Banasik et al., 2012). Profiles with straight lines are particularly important. Four profiles of such characteristics can be identified in the Kętrzyn district area. The analysis of one of them is shown in Fig. 2. It is clearly visible that points 140 and 150 stand out (Fig. 2a). Excluding these points from the data set smoothens the profile (Fig. 2b). It would be advisable to use a function fit so that the adjustment of the profile to the function can be observed when further points are excluded. In this example, the second-degree polynomial (Eq. (2)) and the Pearson correlation coefficient (R²) are used. The value of 0.7 was assumed to be sufficient correlation to stop further analysis. Basing on this analysis, places to perform control measurements can be indicated. For example, points 140 and 150 have significantly changed their height in comparison to the others and a control measurement could explain the reason for this. For this reason, the authors recommend allocating part of the district funds to this type of measurements. However, this was not done in Kętrzyn district. The above analysis allowed the identification of outlying points for which the investigation should be carried out. The reason for the height discrepancy may be:

change in altitude between measurement campaigns (vertical displacement),

destruction of the benchmark and its reconstruction at a different altitude,

assigning an incorrect height in data-base.

Figure 2

The first thing to do is to check the available documentation on these points (topographic descriptions, orthophotomaps and so on). In this way, you can detect the points displacement due to change in topography or reproduction of damaged markers. Only those points for which no relevant information was found during the office works should be subjected to the field review. At the same time, the field-works should concern only the points of the detailed reference network. The GUGiK should be informed about possible height errors for the basic reference network points.

In Kętrzyn district the update measurements of detailed vertical network were completely omitted. Out of 108 points selected for the development of the conversion model, 104 belonged to the basic reference network. The preliminary analysis showed discrepancies on 14 points of the basic reference network. At these points, the differences between the frames significantly diverged from the others, in three cases with the opposite sense. The remaining points (94) were sufficient to develop the transformation function, with a relatively good distribution in the area of the district (Fig. 3). The only area without a basic reference network is located in the north-eastern part of the district. It is a densely forested area with a low degree of urbanisation.

Figure 3

The next stage of work was development of a transformation model. For this purpose, the second-degree polynomial with six parameters (Eq. (2)) was used.: (2) fΔH(x,y)=a00+a10x+a01y+a11xy+a20x2+a02y2 {f_{\Delta H}}\left( {x,y} \right) = {a_{00}} + {a_{10}}x + {a_{01}}y + {a_{11}}xy + {a_{20}}{x^2} + {a_{02}}{y^2} where: f_ΔH(x, y) – transformation function (heights conversion), a₀₀, a₁₀, . . . – function coefficients, x, y – flat coordinates of the points. This is a widely used function, among the others, to model a quasi-geoid in small areas (Tusat and Mikailsoy, 2018; Walo, 2000). Polynomial optimisation was performed with the use of one of the robust estimation methods – the Danish method. It gives a better fit than the least square adjustment method (Banaś, 2012; Borowski and Banaś, 2018). As a result, the distribution of height shifts improved (Table 3) (Borowski, 2017). The extreme values of displacement between the frames were approximately 2.7 cm, which gives a value below 1 mm/km.

Table 3

The assessment of the accuracy of the model ΔH_{Kronstadt′60–KRON86}: raw data (the points of the reference network) and after adjustments (the Danish method)

Parameter	Raw data	The Danish method
Average [cm]	–8.2	–8.2
Deviation [cm]	0.9	0.6
Maximum [cm]	–6.4	–6.9
Minimum [cm]	–10.5	–9.6
Curtosis	–0.11	–0.04
Skewness	–0.24	–0.01

The height difference ΔH_{Kronstadt′60–KRON86} in the grid nodes 0.1° × 0.1° was calculated on the basis of the model. The grid of this resolution is also used in the national elevation difference model (ΔH_{KRON86–EVRF2007}) (Head Office of Geodesy and Cartography, 2015). The final model ΔH_{Kronstadt′60–EVRF2007} was the sum of the values in the nodes of both grids (ΔH_{Kronstadt‘ 60–KRON86} + ΔH_{KRON86–EVRF2007}). The height at any point of the district resource was calculated by interpolation of the obtained grid, according to the formulas from instruction G-1.10 (Kadaj, 2001). Alternatively, it is also possible to fit the polynomial (2) into the state model. The calculation of the value at any point is then carried out on the basis of both models.

The control of matching the points of basic reference network with the GUGiK model was performed with the use of the author’s software, using interpolations according to the G-1.10 instruction. The average difference between KRON86 and EVRF2007 was less than 0.8 mm. This insignificant value made it possible to use the above-mentioned nodes grid without any additional processing.

The article presents an alternative method of heights conversion from the previous height frames to the PL-EVRF2007-NH frame, which is in use for the detailed reference network. This method uses height transformation, through a function developed for selected adjustment benchmarks. The adopted solution is based on the method used in Cracow district (Banasik et al., 2012), or local modelling of a quasi-geoid (Walo, 2000).

The main advantage of the proposed method is the lower cost of its implementation, with the accuracy acceptable for a detailed vertical reference network. The cost of such study should range from 25 PLN per point (without field works) to 40 PLN per point (limited field-works), for an ‘average’ district of about 1000 points. This is lower than the amount estimated by the GUGiK (case 2, Table 2). It should be assumed that the given amount of 50 PLN per point is clearly underestimated. Due to lower development costs the proposed method may be attractive especially for districts with a limited budget.