Evaluation of Anthropogenic Substrate Variability Based on Non-Destructive Testing of Ground Anchors

The purpose of this paper is to describe the variability of soil rheological properties based on research carried out using load tests of ground anchors under complex geotechnical conditions. The heterogeneity of soil should always be considered when designing geotechnical constructions. In the present case, the earthwork created at the Warsaw Slope revealed an embankment of anthropogenic origin, located in a geologically and geomorphologically complex area of the Vistula valley slope. Excavation protection was anchored mainly in soils of anthropogenic origin. When the acceptance tests of the ground anchor were completed, the subsoil randomness was confirmed using nondirect, geostatistical methods. A standard solid rheological model with nonlinear fitting to the data was used. This model was established to describe the creeping activity of the ground anchor more accurately. The characteristics of man-made embankments were described using the parameters obtained with the rheological model of the subsoil.


Introduction
In this paper, the problem of recognizing substrate variability is investigated using a nondirect method. This task is often associated with difficult access to the adjacent subsoil, conservation or preservation of buildings. There are some relevant methods for these tests is the measurement of displacement or strain. One of the common applications of geodetic methods in relation to geotechnical structures is the measurement of the indirect determination of the material properties of a substrate. Nondestructive testing of geotechnical structures provides data on the structure-soil interaction. The scope of displacement and deformation of soil-structure systems [1] used for slow-variable phenomena. Deformation analysis based on multi-temporal terrestrial laser scanning (TLS) surveys has been applied for many years to investigate commercial and academic problems [2]. Close-range photogrammetry with image processing can be used to accurately measure ground and structure deformation [3]. Geodetic methods have found direct application in the design process (observational methods) [4], as well as in structure monitoring and maintenance control [5]. It is difficult to use conventional instruments to measure the shear stress (deformation) along geotechnical structures, such as pile foundations, soil nails and retaining walls, due to their inherent disadvantages, such as large size and electromagnetic interference. Difficult access and problems with securing wires and connections during construction also contribute to the complexity of this process. To overcome these obstacles, fiber Bragg grating (FBG) sensors are recommended for measuring the shear stress along a pile anchor in pull-out tests. The different types of fiber-optic sensing technology include FBG, low-coherence interferometry (LCI), optical time domain reflectometry (ODTR), Fabry-Perot interferometry (FPI) and Brillouin scattering. In [6], a series of pull-out tests were conducted on a model pile in soil with different initial water contents and different pile surface roughnesses, ranging from smooth to rough, to study their influences on pile skin friction. Changes in axial strain and skin friction were measured by FBG sensors using a laboratory test procedure. In this paper, the concept of an analytical diagnostic tool for soil anchors is presented. Anthropogenic ground, which is often difficult to identify in practice, with numerous overlaps and discontinuities, was treated together with anchorage material. Information about the variability of the characteristics of the homogenized medium was obtained. In addition, the authors noticed many possibilities of applying the diagnostic method to the inspection of historical objects. In the presented work, it is the object and its scale that drive the search for new methods to describe variability. Scientific information on both the spatial variability and distribution of anthropogenic soil properties is necessary to model soil-structure systems and to make decisions regarding the need for modifications, strengthening or replacement [7]. Geotechnical engineering has a considerable arsenal of methods for parameter modification. However, little is known about the spatial distribution and variability of soil properties on the scale of the study site, and historical data queries introduce even more uncertainty about quality. This information is confirmed by geological and geotechnical reconnaissance. The purpose of this study is to investigate the spatial relationships and variability of bedrock properties using geostatistical methods and by means of nondestructive measurements of rheological characteristics.
For decades, the area of the capital city of Warsaw has undergone intensive urbanization processes, including a constant influx of people, the development of processing industry, the extraction of mineral aggregate and damage from military conflicts, as well as subsequent reconstruction. These activities caused major changes in the natural environment, especially in the subsoil. Anthropogenic changes in the environment with respect to geological structure and terrain relief, climate, hydrographic conditions, soils and flora have considerable influence on civil structures. Areas with anthropogenic soils had already been marked in the geological documentation of Warsaw in 1937 [8]. Urbanization processes have been occurring since the establishment of the city at the turn of the 13th and 14th centuries, peaking in the 19th century and in the 1970s of the 20th century.
The contemporary environment is not free from anthropogenic transformations, and the earthwork in question is an example of this. Anthropogenic changes also take place on the surface of the terrain, where they have an impact on civil structures. Many geologists [9], construction engineers [10,11] and geographers [12] have examined anthropogenic soils in urbanized areas, as well as outside of them. The high interest in this matter is due to the highly diverse properties of soils that make up subsoils with altered characteristics owing to human activity. These soils are a lithological foundation, especially in urbanized areas, and a subsoil for engineering constructions. Their presence also indicates the transformation of the terrain relief. Morphodynamic processes, which are hazardous to the environment, occur in areas with anthropogenic soils and require the thorough design of geotechnical structures. The characteristics of anthropogenic embankments indicate their diverse composition in terms of both vertical and horizontal distribution. Some of these embankments are natural soils of various origins and compositions (obtained from different sites and deposited), and some are products of human activity. During the long process of city development, the material itself and the means of its accumulation underwent many changes. It was influenced, among other events, by fires, floods and damage that occurred over the centuries. The task of extracting anthropogenic soils from backgrounds of different origins is often complex. The relative proportions of anthropogenic material and the soils of natural origin reveal anisotropy in both the vertical and horizontal directions. The upper part of the anthropogenic layer usually has a higher level of heterogeneity of mechanical characteristics in the plan and contains more material with higher average grain size distribution. The classification of soils is specified in the PN-EN ISO 14688-1 standard [13]. A schematic division of the soils according to this standard is shown in Figure 1 below. higher average grain size distribution. The classification of soils is specified in the PN-EN ISO 14688-1 standard [13]. A schematic division of the soils according to this standard is shown in Figure 1 below. Figure 1. The division of anthropogenic soils according to standard [6].
Only the following types of soil can be considered for the construction subsoil: anthropogenic soils in tiers I and II, present in the form of post-mining and post-energy sites [14], the construction debris of the re-cultivated excavation sites in tier III, and macrolevel areas. It is recommended to increase the development of foundations with low load per unit area, with the type of foundation adjusted to the heterogeneous subsoil, namely, a strip foundation with enhanced reinforcement, grillage or slab foundation. Anthropogenic soils in tiers I and II have become increasingly popular as a construction material for engineering objects, such as railway and road embankments, dams with small damming, embankments of wet landfills or levees. These soils are used in technical recultivation, for example, to fill troughs above underground coal and sulfur excavation sites and to fill the excavations themselves. They are also used as an insulation layer or an ingredient of mineral insulation layers and as an ingredient of soil compositions to increase their endurance and physical characteristics. When it comes to soils in tier I, two types of embankments should be distinguished: a controlled and uncontrolled one, with the former being an embankment formed according to a particular technology, in which its geometry, density and the assumed forming techniques can be controlled.
Anthropogenic soils come in various forms, depending on the technology of their storage (dumping) and usage. At present, dumping (external and internal) is used to create landfills [15,16], heaps and earthwork structures [17][18][19] (e.g., road embankments, levees). Anthropogenic soils, as human-occupied layers, have also shaped the surface and subsoil in cities. Geological engineering research on these soils should be carried out in cooperation with archaeologists. Traces of human activity are often the criterion that shows the anthropogenic origin of geological layers. Due to the many possibilities of anthropogenic soil usage, the research aims and scopes are extremely diverse [20,21]. When it comes to engineering constructions, as with natural or native soils, it is crucial to investigate compaction parameters, as well as the choice of equipment and compaction technology [22], but it is also important to take into account the possible geochemical transitions of built-in soils and the degree of their aggressiveness towards concrete and steel constructions.
Soil transformation in processes I and II has a substantial influence on the engineering solutions used during the transitions and renovations of buildings. This requires a specific understanding of geodiversity, which synthesizes information on the abiotic layers and reveals the spatial interactions and mutual connections [23,24]. Only the following types of soil can be considered for the construction subsoil: anthropogenic soils in tiers I and II, present in the form of post-mining and post-energy sites [14], the construction debris of the re-cultivated excavation sites in tier III, and macrolevel areas. It is recommended to increase the development of foundations with low load per unit area, with the type of foundation adjusted to the heterogeneous subsoil, namely, a strip foundation with enhanced reinforcement, grillage or slab foundation. Anthropogenic soils in tiers I and II have become increasingly popular as a construction material for engineering objects, such as railway and road embankments, dams with small damming, embankments of wet landfills or levees. These soils are used in technical recultivation, for example, to fill troughs above underground coal and sulfur excavation sites and to fill the excavations themselves. They are also used as an insulation layer or an ingredient of mineral insulation layers and as an ingredient of soil compositions to increase their endurance and physical characteristics. When it comes to soils in tier I, two types of embankments should be distinguished: a controlled and uncontrolled one, with the former being an embankment formed according to a particular technology, in which its geometry, density and the assumed forming techniques can be controlled.
Anthropogenic soils come in various forms, depending on the technology of their storage (dumping) and usage. At present, dumping (external and internal) is used to create landfills [15,16], heaps and earthwork structures [17][18][19] (e.g., road embankments, levees). Anthropogenic soils, as human-occupied layers, have also shaped the surface and subsoil in cities. Geological engineering research on these soils should be carried out in cooperation with archaeologists. Traces of human activity are often the criterion that shows the anthropogenic origin of geological layers. Due to the many possibilities of anthropogenic soil usage, the research aims and scopes are extremely diverse [20,21]. When it comes to engineering constructions, as with natural or native soils, it is crucial to investigate compaction parameters, as well as the choice of equipment and compaction technology [22], but it is also important to take into account the possible geochemical transitions of built-in soils and the degree of their aggressiveness towards concrete and steel constructions.
Soil transformation in processes I and II has a substantial influence on the engineering solutions used during the transitions and renovations of buildings. This requires a specific understanding of geodiversity, which synthesizes information on the abiotic layers and reveals the spatial interactions and mutual connections [23,24].
The complex composition of the Masovian Lowland, which, as a whole, contains traces of glacier activity from the Riss glaciation, as well as signs of accumulative and erosive activity of river and glacial waters during the interglacial periods [25], makes geotechnical conditions more complicated. Previous geological documentation presents the consequences of processes that have been transpiring for millions of years, and only in the last few centuries have there been substantial changes to the structure, composition and morphology of this terrain. During the construction and the many years of use of the Warsaw Citadel, different substances were introduced to the environment, including ammunition, waste, food and construction materials [26]. With the development of fortification techniques and the changes taking place in the international arena, the fortress had been expanded by the erection of forts and external positions near the end of the nineteenth century, which required a huge scope of earthwork. It must be highlighted that anthropogenic interference is rather shallow, and on the profile scale, it concerns the surface layer only. Ground fortifications are still visible today and are a historic monument [27]. The study of the development of earthwork methods has led us to the conclusion that contemporary knowledge of the parameters and mechanics of soils is slowly changing the tradition of raising earth constructions using empirical experience and corrections during the time of use. When designing ground anchors and the soldier pile wall, both the geodiversity and the anthropogenic origin of soils should be taken into account [28]. The influence of anthropogenic soils on the design of geotechnical constructions was shown in [29], among other studies. Studies on deformations in retaining structures and retrospective analyses have shown that anthropogenic soils that had initially been classified as low bearing have presented average durability parameters.
The presence of cultural layers is highly visible when performing deep excavations that reveal the spatial heterogeneity of the ground and its interactions with contemporary constructions.

Research Object
The geotechnical material under investigation originates from the historical embankments of a citadel located in Warsaw, the capital city of Poland, located in the Mazovian Lowland. The trough is filled with sediments from the Cretaceous, Tertiary and Quaternary Periods and shaped by geomorphological processes that took place in the Quaternary. The main geomorphological units within the city of Warsaw are the Warsaw Plain, the Wołomin Plain and the Central Vistula Valley. The Warsaw and Wołomin Plains were created by glaciers in the Riss glaciation period and by the accumulative and erosive activity of river and glacier waters in the interglacial periods. The predominating materials in the area of the Warsaw and Wołomin Plains are the Pleistocene glacial tills (from sandy clays through clayey sands). The highest peak of the plain is about 115 m above sea level and is located in the city center. In the East, the plain is divided by the Vistula river, where the Warsaw Slope's cliff is located [30,31]. The view of the Warsaw Slope with the city panorama is characteristic of the city. The slope divides the city into north/south sides. Recreational areas, parks and palaces that survived WWII are located nearby. Currently, the Vistula valley is created by the riverbed and the following terraces: floodplain, fluvial, dune and lacustrine. Most of the right-bank part of Warsaw is situated on the terraces along the Vistula river. These terraces are made of alluvial sediments, such as gravels, sands, river sediments, fluvisols and organic sediments. Dune and eolian sands are present in the fluvial terraces and make vast areas of sorted sands and dunes. The citadel was raised in the left-bank part of Warsaw, directly by the river bed, in the area of the Warsaw Slope and the Warsaw Plain terraces. In 1830, the November Uprising against Russian domination broke out [32]. Warsaw became the most important and long-lasting point of resistance for the uprising. As a result, the construction of a stronghold that could be a resistance point for the Russian army guarding the capital city of the Kingdom of Poland became a necessity. It had an entirely non-metropolitan character and function; however, due to modifications to its military functions and the endurance of its technical substance, it has managed to survive to this day [33].
The Warsaw Citadel was raised in theŻoliborz district, which was located in the 16th-century Fawora manor farm and the Piarists Monastery from the 17th century [34].
Additionally, during archaeological work before the museum construction, a former Protestant cemetery was discovered [35].
The construction of the Citadel in place of the Royal Guard barracks was initiated by Tzar Nicholas I in 1832. The Citadel has the shape of a semicircle, whose chord is based on the river (Figure 2). The whole stronghold is surrounded by a scarp with a dry moat and counterscarp. There are three bastions and two half-bastions. The arms of each of them included a single-level cannoneer position. The solutions used in the Citadel were based on schemes popular in the 19th century [36,37]. There was neither a long-range defense system nor protection from the artillery, no ammunition bunker and no defensive barracks. The ceilings had no earthen ramparts, with only a thin roof covered with ceramic ones for larger areas. The Citadel was not meant to protect against the regular army, but only against insurgents with no artillery whom the Tzar expected to come from the conquered Warsaw. It was designed as a display of the power of the Tzar and to provide a vast and capacious space. On the first day of the January Uprising, the garrison had 16,000 soldiers. A crucial element of the Citadel was its artillery, whose aim was to keep the center of Warsaw under control: this threat worked well during the January Uprising, as the temporary government decided not to provoke a rebellion in the city [38]. due to modifications to its military functions and the endurance of its technical substance, it has managed to survive to this day [33].
The Warsaw Citadel was raised in the Żoliborz district, which was located in the 16th-century Fawora manor farm and the Piarists Monastery from the 17th century [34].
Additionally, during archaeological work before the museum construction, a former Protestant cemetery was discovered [35].
The construction of the Citadel in place of the Royal Guard barracks was initiated by Tzar Nicholas I in 1832. The Citadel has the shape of a semicircle, whose chord is based on the river ( Figure 2). The whole stronghold is surrounded by a scarp with a dry moat and counterscarp. There are three bastions and two half-bastions. The arms of each of them included a single-level cannoneer position. The solutions used in the Citadel were based on schemes popular in the 19th century [36,37]. There was neither a long-range defense system nor protection from the artillery, no ammunition bunker and no defensive barracks. The ceilings had no earthen ramparts, with only a thin roof covered with ceramic ones for larger areas. The Citadel was not meant to protect against the regular army, but only against insurgents with no artillery whom the Tzar expected to come from the conquered Warsaw. It was designed as a display of the power of the Tzar and to provide a vast and capacious space. On the first day of the January Uprising, the garrison had 16,000 soldiers. A crucial element of the Citadel was its artillery, whose aim was to keep the center of Warsaw under control: this threat worked well during the January Uprising, as the temporary government decided not to provoke a rebellion in the city [38]. Architects of the WCXA studio are transforming the Citadel into a tourist attraction and a memorial. The new museums that will be located there will change the defensive character of the land development. The walls and ramparts must be open for vehicle and pedestrian passage (Figures 2 and 3). Geotechnical works, as a stage of this project, are shown in Figure 4. Architects of the WCXA studio are transforming the Citadel into a tourist attraction and a memorial. The new museums that will be located there will change the defensive character of the land development. The walls and ramparts must be open for vehicle and pedestrian passage (Figures 2 and 3). Geotechnical works, as a stage of this project, are shown in Figure 4.
According to a design concept selected in an architectural contest in 2010, the following parts of the Warsaw Citadel will be used for museum purposes: The Caponier, which will be the head office of the Katyń Museum, with a permanent exhibition [39]; the Shoulder Battery; and the Nowomiejska Gate (the main entrance to the museum and the open-air area). In 2022, the Museum of the Polish Army will be fully relocated to the Warsaw Citadel. The new facilities will generate significant visitor traffic, especially car drivers, hence the concept of opening up the site and building a new tunnel through the earth embankment. ials 2021, 14, x FOR PEER REVIEW 6 of 21  According to a design concept selected in an architectural contest in 2010, the following parts of the Warsaw Citadel will be used for museum purposes:  According to a design concept selected in an architectural contest in 2010, the following parts of the Warsaw Citadel will be used for museum purposes:

Geotechnical Site Investigation
The geotechnical conditions in the substrate and the quality of the material become important in the construction and use of massive structures, including fortifications, in particular embankments and the protection of deep excavations. The variable layering visible during the earthworks ( Figure 5) is a confirmation of the former geotechnical investigation.
The Caponier, which will be the head office of the Katyń Museum, with a permanent exhibition [39]; the Shoulder Battery; and the Nowomiejska Gate (the main entrance to the museum and the open-air area). In 2022, the Museum of the Polish Army will be fully relocated to the Warsaw Citadel. The new facilities will generate significant visitor traffic, especially car drivers, hence the concept of opening up the site and building a new tunnel through the earth embankment.

Geotechnical Site Investigation
The geotechnical conditions in the substrate and the quality of the material become important in the construction and use of massive structures, including fortifications, in particular embankments and the protection of deep excavations. The variable layering visible during the earthworks ( Figure 5) is a confirmation of the former geotechnical investigation. The area of Warsaw Slope terraces includes regions that have active erosive processes and creeping movements, as well as places where there are springs, wetlands areas and depressions of suffosion origin. These factors adversely affect the foundations of buildings and engineering structures, in particular: The occurrence of an anthropogenic soil layer in the substrate; • The occurrence of a layer of low-bearing soils in the substrate; • An underground water table with a shallow level.
Geotechnical investigations of the substrate in the area of the designed cutout included the geotechnical drilling hole and cone penetration tests (CPT) [40]. As a result of anthropo-pressure, which especially intensified during the Citadel's erection, the primary terrain was leveled and additionally reinforced by building up a layer of anthropogenic embankments. The sediments lithologically comprise the following materials: clayey sands, sandy clays and silty clay with stiff-consistency gravels and cobbles. All of the above-mentioned sediments are of Quaternary origin. Older subsoil, i.e., from the Tertiary and Pliocene periods, was found at the bottom of the scarp.
For design purposes, the following basic geotechnical layers have been established ( Figure 6):  The area of Warsaw Slope terraces includes regions that have active erosive processes and creeping movements, as well as places where there are springs, wetlands areas and depressions of suffosion origin. These factors adversely affect the foundations of buildings and engineering structures, in particular: The occurrence of an anthropogenic soil layer in the substrate; • The occurrence of a layer of low-bearing soils in the substrate; • An underground water table with a shallow level.
Geotechnical investigations of the substrate in the area of the designed cutout included the geotechnical drilling hole and cone penetration tests (CPT) [40]. As a result of anthropopressure, which especially intensified during the Citadel's erection, the primary terrain was leveled and additionally reinforced by building up a layer of anthropogenic embankments. The sediments lithologically comprise the following materials: clayey sands, sandy clays and silty clay with stiff-consistency gravels and cobbles. All of the above-mentioned sediments are of Quaternary origin. Older subsoil, i.e., from the Tertiary and Pliocene periods, was found at the bottom of the scarp.
For design purposes, the following basic geotechnical layers have been established ( and sandy clays with gravels and cobbles with a stiff consistency and plasticity index IL = 0.1-0.2.

•
Layer IV: glaciolacustrine deposits of the Wartanian Glaciation formed as silts, sandy silts, silty clays and clays with a stiff consistency. • Layer V: consolidated cohesive glacial deposits of the Odranian Glaciation consisting of sandy clays and clayey sands with gravels and cobbles with a very stiff consistency. Figure 6. Cross-section along the road support (with courtesy of Geoteko data [40]).
For the design task under consideration, the layering between the two drilled test holes was adopted. Reconnaissance was carried out to confirm the adequate capacity of the soil anchors and to ensure their durability for a period of two years. As mentioned above, the embankments are made of local soils found close to the fortifications. The substrate classified as original can also be anthropogenic. The current layer layout is revealed only at the time of excavation. The carefully observed walls of the excavation make it possible to view the true arrangement of geological layers, often different from the one drawn on the basis of research drills. The attempt to adopt horizontal layer boundaries is visible in most geological engineering cross-sections included in the documentation. These basic assumptions work well in the case of a simple and complex geological structure, but in the case of a complicated structure, due to the glaciotectonic origin of ground displacements, they are not always correct. An example is the geological structure presented in the cross-sections generated for the needs of the renovation works under consideration.
The described adaptation requires a cross-cut in the embankment and protection of the excavation in the embankment. The first stage of implementation of the entrance to the designed museum object requires the demolition of approximately 54 m of the Carnot wall. The excavation under investigation was carried out in an uncontrolled anthropogenic substrate and in non-cohesive and cohesive subsoil. The task requires the identification of a substrate made of anthropogenic soils in which excavation and estimation of geotechnical parameters will be carried out. After the Carnot wall had been demolished and the soldier pile wall was completed, the work was in progress to create the excavation. The further stages of the structure are presented in Figures 5 and 7. Figure 6. Cross-section along the road support (with courtesy of Geoteko data [40]).
For the design task under consideration, the layering between the two drilled test holes was adopted. Reconnaissance was carried out to confirm the adequate capacity of the soil anchors and to ensure their durability for a period of two years. As mentioned above, the embankments are made of local soils found close to the fortifications. The substrate classified as original can also be anthropogenic. The current layer layout is revealed only at the time of excavation. The carefully observed walls of the excavation make it possible to view the true arrangement of geological layers, often different from the one drawn on the basis of research drills. The attempt to adopt horizontal layer boundaries is visible in most geological engineering cross-sections included in the documentation. These basic assumptions work well in the case of a simple and complex geological structure, but in the case of a complicated structure, due to the glaciotectonic origin of ground displacements, they are not always correct. An example is the geological structure presented in the cross-sections generated for the needs of the renovation works under consideration.
The described adaptation requires a cross-cut in the embankment and protection of the excavation in the embankment. The first stage of implementation of the entrance to the designed museum object requires the demolition of approximately 54 m of the Carnot wall. The excavation under investigation was carried out in an uncontrolled anthropogenic substrate and in non-cohesive and cohesive subsoil. The task requires the identification of a substrate made of anthropogenic soils in which excavation and estimation of geotechnical parameters will be carried out. After the Carnot wall had been demolished and the soldier pile wall was completed, the work was in progress to create the excavation. The further stages of the structure are presented in Figures 5 and 7.  Materials 2021, 14, x FOR PEER REVIEW 9 of 21  (Figure 6). The fixed anchors are positioned in the zone of occurrence of the Wartanian Glaciation's sediments, represented by fluvioglacial and glaciolacustrine sands (fine sands, silty sands and medium sands), glacial tills (clayey sands, sandy clays with gravels, cobbles and boulders) and locally lacustrine soils (silt, silty clay and clay). The total thickness of these sediments ranges from 5.0 to more than 10.0 m. The anchor capacity was predicted and designed as a deterministic value based on empirical formulation. The final verification of anchorages was subsequently performed using creep tests after pre-compression of anchors by a static pull-out load test.

The Principles of Design and Testing
As part of site investigations, excavations, ground anchor testing and other engineering works, boreholes are crucial for understanding Warsaw Citadel's geodiversity. The designers used anchors as temporary supports for soldier pile walls ( Figure 5). During construction, it was possible to look closer at the substrate and uncover some mechanical dependencies [41,42]. The standard [43] provides the general definition: an anchor is a structure capable of carrying tensile loads through a free length of a tendon onto a support layer. On the other hand, for an injection anchor, it is specified that an anchor with a bolt solidified by the injection of resin, cement grout or other material transfers the tensile force to the ground. The standard [44] applies to temporary or permanent soil anchors that are ground-bonded by injection, compressed and tested according to the strictest safety measure described. The field test of the anchor head is shown in Figure 8. The measuring set consists of a dial indicator sensor and an actuator (conventional hollow-piston cylinder).
Ground anchors mainly carry pulling forces, sometimes slightly transverse, and do not transfer compressive forces [45]. The main areas of their application include the stabilization of an excavation and other support structures. This stretch does not make the surrounding soil stronger. Typical ranges of anchor load capacity vary from 500 kN to 2.0 MN. The main document referenced for the testing procedure in this study is the standard [46]. It presents a complete range of design calculations and implementation procedures for both ultimate and serviceability limit states. It concerns only the design and testing of anchors with free length. However, at present, there is no effective method for monitoring the anchors. The  (Figure 6). The fixed anchors are positioned in the zone of occurrence of the Wartanian Glaciation's sediments, represented by fluvioglacial and glaciolacustrine sands (fine sands, silty sands and medium sands), glacial tills (clayey sands, sandy clays with gravels, cobbles and boulders) and locally lacustrine soils (silt, silty clay and clay). The total thickness of these sediments ranges from 5.0 to more than 10.0 m. The anchor capacity was predicted and designed as a deterministic value based on empirical formulation. The final verification of anchorages was subsequently performed using creep tests after pre-compression of anchors by a static pull-out load test.

Anchor Test Description The Principles of Design and Testing
As part of site investigations, excavations, ground anchor testing and other engineering works, boreholes are crucial for understanding Warsaw Citadel's geodiversity. The designers used anchors as temporary supports for soldier pile walls ( Figure 5). During construction, it was possible to look closer at the substrate and uncover some mechanical dependencies [41,42]. The standard [43] provides the general definition: an anchor is a structure capable of carrying tensile loads through a free length of a tendon onto a support layer. On the other hand, for an injection anchor, it is specified that an anchor with a bolt solidified by the injection of resin, cement grout or other material transfers the tensile force to the ground. The standard [44] applies to temporary or permanent soil anchors that are ground-bonded by injection, compressed and tested according to the strictest safety measure described. The field test of the anchor head is shown in Figure 8. The measuring set consists of a dial indicator sensor and an actuator (conventional hollow-piston cylinder).
ning and end of each load step.
After reaching the test tension, the creep rate α is determined by: where S2 and S1 are creep displacement at times t2 and t1, respectively. The anchors were installed using the Hütte HBR 609-3 set. This is a hydraulic drilling rig for micropiles and ground anchors, foundation construction, soil investigation and slope stabilization. The most important parameters are engine power 209 kW, extraction force 100 kN, crowd force 100 kN and working extraction speed 17.5 m/min. The diameter of the casing used was 133 mm. The device performs the entire technological process of creating the anchor quickly and efficiently. The anchors were drilled with a gimlet bit in pipe casing with water flush. Upon making a borehole, grout injection was carried out and then the anchor bearing element-a set of wires-was installed with double corrosion protection (DCP) to ensure a prolonged lifespan; the element was equipped with spacers, injection hoses and separator(s). Primary pressure grouting was applied through the casing every 2 m during its removal, and then the pressure applied was 1 MPa. After approximately 24 h, all anchors were grouted by tube-a-manchette with a pressure of 6-7 N/mm 2 . Ground anchors mainly carry pulling forces, sometimes slightly transverse, and do not transfer compressive forces [45]. The main areas of their application include the stabilization of an excavation and other support structures. This stretch does not make the surrounding soil stronger. Typical ranges of anchor load capacity vary from 500 kN to 2.0 MN. The main document referenced for the testing procedure in this study is the standard [46]. It presents a complete range of design calculations and implementation procedures for both ultimate and serviceability limit states. It concerns only the design and testing of anchors with free length. However, at present, there is no effective method for monitoring the anchors. The three coordinated standards [43,44,46] assessing the design, execution and testing of grouted anchors, respectively, in the EU, establish three types of tests: investigation tests, suitability tests and acceptance tests. The standard [46] includes 3 test methods: TM1, TM2 and TM3. All of these tests are designed to evaluate the global performance of the anchors by taking external measurements (e.g., applied force, relaxation and creep). In particular, the investigation test is a load test aimed at evaluating the ultimate geotechnical resistance of the anchor and its behavior at working loads.
This article discusses an acceptance test according to the TM3 model (Figure 9b). The procedure must be carried out on all anchors. The anchor is loaded in incremental steps from a datum to the maximum load (here, the design capacity of the anchor). The incremental steps were 10%, 30%, 50%, 80% and 100% of the design capacity. The displacement of the anchor points was measured under a maintained load at each loading step. The measurements summarize the anchor head displacement versus anchor load at the beginning and end of each load step. During primary and secondary pressure grouting, the grout used had w/c = 0.50 and was produced from CEMII 32.5R. Finally, 7 days later, the anchors were tested, compressed and stabilized with anchorage ( Figure 9a). System creep tests were performed after pre-compression of anchors.

Modeling Methodology for Time-Dependent Process
The genetic and geomorphological diversity of the subsoil material not only affects the physical characteristics of the substrate in which the fixed anchor tests are performed but also presents difficulties in describing its rheological characteristics [47]. There are many methods of creep testing for soils and soft rocks: these are methods for loading large-scale cylindrical loads [48]. They require a long time with a special apparatus and a constant climate (controlled pressure, temperature and humidity). Research procedures are available in [49]. As a different approach, forced testing after destruction on the direct surface of slip is possible in a modified rotary shear apparatus [50]. In the case of anchor tests lowered in the ground at the construction site, an attractive option is the in situ test presented in the article. The authors used the method of direct measurement of these properties based on creep tests [46]. The standard linear solid model of viscoelastic soil [51,52] was used. The Hook element (first spring 1 ) was connected in series with the Kevin-Voigt element (second spring 2 and first damper 2 ). The strains of the elementary model components After reaching the test tension, the creep rate α is determined by: where S 2 and S 1 are creep displacement at times t 2 and t 1 , respectively. The anchors were installed using the Hütte HBR 609-3 set. This is a hydraulic drilling rig for micropiles and ground anchors, foundation construction, soil investigation and slope stabilization. The most important parameters are engine power 209 kW, extraction force 100 kN, crowd force 100 kN and working extraction speed 17.5 m/min. The diameter of the casing used was 133 mm. The device performs the entire technological process of creating the anchor quickly and efficiently. The anchors were drilled with a gimlet bit in pipe casing with water flush. Upon making a borehole, grout injection was carried out and then the anchor bearing element-a set of wires-was installed with double corrosion protection (DCP) to ensure a prolonged lifespan; the element was equipped with spacers, injection hoses and separator(s). Primary pressure grouting was applied through the casing every 2 m during its removal, and then the pressure applied was 1 MPa. After approximately 24 h, all anchors were grouted by tube-a-manchette with a pressure of 6-7 N/mm 2 . During primary and secondary pressure grouting, the grout used had w/c = 0.50 and was produced from CEMII 32.5R. Finally, 7 days later, the anchors were tested, compressed and stabilized with anchorage ( Figure 9a).
System creep tests were performed after pre-compression of anchors.

Modeling Methodology for Time-Dependent Process
The genetic and geomorphological diversity of the subsoil material not only affects the physical characteristics of the substrate in which the fixed anchor tests are performed but also presents difficulties in describing its rheological characteristics [47]. There are many methods of creep testing for soils and soft rocks: these are methods for loading large-scale cylindrical loads [48]. They require a long time with a special apparatus and a constant climate (controlled pressure, temperature and humidity). Research procedures are available in [49]. As a different approach, forced testing after destruction on the direct surface of slip is possible in a modified rotary shear apparatus [50]. In the case of anchor tests lowered in the ground at the construction site, an attractive option is the in situ test presented in the article. The authors used the method of direct measurement of these properties based on creep tests [46]. The standard linear solid model of viscoelastic soil [51,52] was used. The Hook element (first spring γ 1 ) was connected in series with the Kevin-Voigt element (second spring γ 2 and first damper θ 2 ). The strains of the elementary model components (Hook and Kevin-Voigt) were summarized as a one-dimensional process, along the anchor tendon, to the total value. The calculation was conducted for a constant load value in the anchor head, where stress is σ = P 0 /A, with A as the cross-sectional area of the anchor tendon and P 0 as the axial force. According to the Code [46], the total time of the experiment was 30 minutes (Figure 9c), and the lack of long-term creep results was addressed by the assumption that σ is constant, where σ represents compressive stress. The boundary task has strain ε(t) equal to ε 0 at time t 0 = 0, and it allows us to determine the constant of elasticity of the first spring γ 1 (Hook's element) as γ 1 = σ/ε(0).

Results
Measured displacement u [mm] in this study was divided by the fixed anchor length L tot for calculating the anchor strain. This assumption was made because the anchortendon-soil system was analyzed. The Kevin-Voigt part has two sub-elements: liquid and solid body types. Their stress values are summarized with the output values. Simultaneously, the strains of these two submodels are equal. The first-order differential equation was obtained as: At the final time of the experiment t l , it was assumed that the second damper had disappeared . ε(t l ) = 0, which allows the second spring coefficient γ 2 to be deduced from (2): which allows the solution to be obtained with one unknown parameter θ 2 : To test the spatial variability of the model parameters, semivariograms were used and are presented graphically in Figure 10 and by direct statistical parameters in Table 1. The minimum distance between the testing points is compatible with the anchor spacing and is equal to 2.5 m. The semivariogram data were taken along the wall of the south and north wings of the excavation. The two sides constitute the cross-section of the embankment. The shape of the semivariograms does not allow for building a model for the empirical variability of the tested characteristics along the cross-section of the excavation. For both sides under consideration, the randomness of the model parameter value is quite high, and along the walls of the excavation, it can be regarded as white noise, with no functional correlation between the variance and the distance. This confirms the previous assumptions that the soils creating the embankment are highly heterogeneous, which also makes the distinction of homogeneous layers in the embankment difficult.

Both Wings Mean
Median SD  The correlation between model parameters was examined, and the results are p sented in Table 2. The model parameters 2 , 2 have a significant positive correlation. T statistically significant dependence of parameters describing the viscous part in the mo indicates a similar material characteristic for time-dependent strains. Basic statistical relationships between model parameters are presented graphically Figure 11.  The correlation between model parameters was examined, and the results are presented in Table 2. The model parameters γ 2 , θ 2 have a significant positive correlation. The statistically significant dependence of parameters describing the viscous part in the model indicates a similar material characteristic for time-dependent strains. Basic statistical relationships between model parameters are presented graphically in Figure 11. For the analysis of the variability of anchorage characteristics for the south and north wing sections, a geostatistical measure of variability was used to take into account the mutual position of the anchors. The semivariogram describes the degree of dependence of the feature (rheological and elastic parameters) as a function of distance in normalized Euclidean space ||h|| for isotropic phenomena, or as a function of distance and direction when anisotropy is assumed for phenomena in two or more dimensions. The semivariogram estimator has the form: In Equation (5), N(h) denotes the number of all pairs (z(si) − z(sj)) 2 that are separated by distance ||h||; z(s) denotes the value of the feature at location s; z denotes the set of all discrete values of that feature; and i, j are natural numbers in the interval [1,n], where n is the number of measured points. For practical reasons [53] and the symmetry of the phenomenon with respect to zero spacing, a semivariogram, defined as half of the variogram γ(h), was used. The semivariogram is a measure of dissimilarity between points observed For the analysis of the variability of anchorage characteristics for the south and north wing sections, a geostatistical measure of variability was used to take into account the mutual position of the anchors. The semivariogram describes the degree of dependence of the feature (rheological and elastic parameters) as a function of distance in normalized Euclidean space ||h|| for isotropic phenomena, or as a function of distance and direction when anisotropy is assumed for phenomena in two or more dimensions. The semivariogram estimator has the form: In Equation (5), N(h) denotes the number of all pairs (z(s i ) − z(s j )) 2 that are separated by distance ||h||; z(s) denotes the value of the feature at location s; z denotes the set of all discrete values of that feature; and i, j are natural numbers in the interval [1,n], where n is the number of measured points. For practical reasons [53] and the symmetry of the phenomenon with respect to zero spacing, a semivariogram, defined as half of the variogram γ(h), was used. The semivariogram is a measure of dissimilarity between points observed at a given location z(s i ) and z(s j ). Figure 12 shows points on the empirical semivariogram for the south wing, and Figure 13 shows points on the empirical semivariogram for the north wing. In both figures, the lines were fitted by the least square method (LSM). The high value of the nugget effect in the semivariogram analyzed in Figure 12 and Table 3 and its linear theoretical model indicates significant randomness and diversity of the parameter 1 values for the south wing, while the character of the variation in this parameter for the north wing is presented differently, as shown in Figure 13 and Table 4. Here, the semivariogram is described by a spherical model, and a smaller value of the nugget effect can be observed, which means greater homogeneity of the material. Differences in the origin of the material and methods used to erect the embankments surrounding the Warsaw Citadel are evident on both sides associated with adjacent masses of anthropogenic embankments subjected to the analysis of the variability of elastic features of the rheological model of the subsoil. The spatial variability of the feature 2 in Figures  12b and 13b and 2 in Figures 12c and 13c is supported by significant divergences between the analyzed semivariograms.  The high value of the nugget effect in the semivariogram analyzed in Figure 12 and Table 3 and its linear theoretical model indicates significant randomness and diversity of the parameter γ 1 values for the south wing, while the character of the variation in this parameter for the north wing is presented differently, as shown in Figure 13 and Table 4. Here, the semivariogram is described by a spherical model, and a smaller value of the nugget effect can be observed, which means greater homogeneity of the material. Differences in the origin of the material and methods used to erect the embankments surrounding the Warsaw Citadel are evident on both sides associated with adjacent masses of anthropogenic embankments subjected to the analysis of the variability of elastic features of the rheological model of the subsoil. The spatial variability of the feature γ 2 in Figures 12b and 13b and θ 2 in Figures 12c and 13c is supported by significant divergences between the analyzed semivariograms.  There was no significant correlation in the remaining combinations of parameter pairs. A lack of correlation with the element responsible for immediate deformations 1 allows for the presumption of different mechanical processes affecting capillaries or grains and influencing the behavior of the material. The human influence on soil formation is much deeper and more intense than we can perceive. For this reason, it is imperative that anthropogenic processes are clearly recognized as a factor in soil formation.  There was no significant correlation in the remaining combinations of parameter pairs. A lack of correlation with the element responsible for immediate deformations γ 1 allows for the presumption of different mechanical processes affecting capillaries or grains and influencing the behavior of the material. The human influence on soil formation is much deeper and more intense than we can perceive. For this reason, it is imperative that anthropogenic processes are clearly recognized as a factor in soil formation.
In the present case, the design of the anchors and their execution took into account the presence of anthropogenic soils with a smaller range than was shown by the acceptance test results. In the performed tests, the load was not increased until reaching the ultimate or serviceability limit state of the anchor; therefore, it is not possible to indicate the mechanism of failure of the grout or soil here. Numerous test programs have been carried out for anchors in concrete and have shown that three different failure modes are considered, including concrete cone failure at the unloaded end without interfacial debonding, interfacial debonding plus concrete cone failure at the same height as the interfacial shear crack tip, and interfacial debonding plus concrete cone failure at the unloaded end. In theoretical studies, linear elastic fracture mechanics is generally used to calculate the cone failure load. A similar study was carried out for soil nails [43], the function of which is identical to that of anchors. The shear strength between the soil nails and the soil is an important parameter for the design of the soil nails. This parameter is governed by a number of factors, such as the stress condition, normal stress acting on the nail-soil properties (strength, particle size, soil dilatancy, degree of saturation, etc.) and nail-soil boundary conditions (surface area and roughness). Multiple nail-soil tests were performed in a pull-out box. Full-scale tests are presented in [47]. All anchors were tested according to [45]. The investigation test using the cyclic load test procedure was conducted, and after reaching the maximum level of the prestressing force in every cycle, the force was kept constant, and the displacement of the tendon was measured and recorded at specified times. In the tests [47], all boreholes for anchors were drilled vertically. Rotary drilling without casing was used to drill the boreholes; however, this is not typical for anchors that are drilled diagonally or horizontally.
A similar process of behavior under the loading and failure mode of tested anchors was observed on all tested members. After reaching the ultimate resistance (maximum value of the prestressing force), a sudden decrease in the prestressing force occurred. This was accompanied by a visible displacement of the grout body out of the ground.
The ultimate capacity of ground anchors is mostly determined by using various empirical and semi-empirical methods [54,55]. These procedures are usually successful but highly simplified. A wide range of empirical factors is required, for which knowledge is limited only to certain types of soil. One of the major simplifications of these methods is the assumption that the bond stress at the soil-grout interface is mobilized uniformly along the whole fixed length. The result of these methods is usually only the ultimate capacity, with no additional information about the displacement required to mobilize this ultimate force. As outlined above, most of the scientific research on anchors involves a pull-out test. These are not common studies, but they are well documented in the literature. Acceptance tests are standardized and performed almost on a mass scale, but anchor displacement and creep velocity are only assessed in terms of compliance with standard values. Acceptance tests are of underestimated value, because analysis of displacements over time enables spatial evaluation of the rheological parameters of the soil.
In view of the numerous factors that affect the load capacity of anchors, the performance of acceptance tests for temporary structures is justified. The use of the full range of anchor acceptance test results is innovative; it allows the assessment of the variability of anthropogenic soil, which is not fully revealed by a typical discrete borehole investigation.

Conclusions
The authors focused on presenting an indirect research method to describe the heterogeneity of embankments. Based on the anchor creep model, the authors identified the heterogeneity of the material, the lack of spatial correlation between its properties and the large dispersion of the (rheological) parameters γ 1 , γ 2 and θ 2 . The variability of the substrate-anchor system was investigated indirectly. The directly obtained experimental results allow the anchors to be accepted as a temporary support element for a retaining wall, according to the standard.
Geostatistical methods were used to investigate the variability of the mechanical features of the subsoil in the cross-section of the south and north wings of the tunnel wall protection. The values of anthropogenic embankment features were determined to have high variability and heterogeneity. The high variability prevented us from describing the rheological features of homogenized subsoil material and anchors by one theoretical semivariogram model. Linear and spherical models were used, and in each case, the effect of nuggets describing the randomness of the phenomenon plays an important role. Even at small distances between anchors, a significant value of semivariance-differentiation of mechanical properties-can be observed. This confirms the thesis of the high variability of parameters and, indirectly, of substrate characteristics. An important conclusion from the presented analysis is the confirmation of the necessity of acceptance testing for each anchor in the studied area. The presented method allows us to reduce the number of necessary anchor tests in the case of more homogeneous soils. In place of in situ testing, numerical simulation (randomization of features based on a semivariogram) can be performed, and the probability of failure can be investigated computationally. Analysis of the rheological model does not support the claim that the anthropogenic and fluvioglacial layers distinguished in the geotechnical documentation are homogeneous. The results obtained indicate major material heterogeneity. The spatial variability of the rheological model parameters means that the thickness of heterogeneous anthropogenic soils is higher than initially assumed; soils of postglacial origin have been relocated to and deposited in their current location. The notion of an uncontrolled embankment should be broadened to encompass the body of the heap as well. The presented research procedure and analysis of results also confirm the suitability of the designer's (co-author of the article) decision to secure the temporary retaining wall with anchors. By adjusting the tension force and subjecting each anchor to an acceptance test, it was possible to achieve partial independence of the significant variability of soil characteristics. The spatial variability of the creep model parameters can be considered a new element in the process of determining the geodiversity of the subsoil. Given the variability of the mechanical features of the anthropogenic subsoil, the description of the geotechnical situation illustrates, first of all, the activities of humans and their omissions in the erection of embankments. It is unlikely that we will be able to discover the detailed reasons for this heterogeneity.