Next Article in Journal
Assessment and Visualization of OSM Consistency for European Cities
Next Article in Special Issue
A GloVe-Based POI Type Embedding Model for Extracting and Identifying Urban Functional Regions
Previous Article in Journal
Examining the Effect of Squatter Settlements in the Evolution of Spatial Fragmentation in the Housing Market of the City of Buenos Aires by Using Geographical Weighted Regression
Previous Article in Special Issue
Simulation of Land-Use Changes Using the Partitioned ANN-CA Model and Considering the Influence of Land-Use Change Frequency
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Roman Aqueduct Flow Estimation Using Geomatic Measurement

1
Chemical, Materials and Environmental Engineering Department (DICMA), Sapienza Università di Roma, I-00184 Rome, Italy
2
Department of Civil, Constructional and Environmental Engineeering (DICEA), Sapienza Università di Roma, I-00184 Rome, Italy
3
MTA-ELTE Geological, Geophysical and Space Science Research Group, 1117 Budapest, Hungary
4
Alba Regia Technical Faculty, Institute of Geoinformatics, Óbuda University, 8000 Székesfehérvár, Hungary
*
Author to whom correspondence should be addressed.
ISPRS Int. J. Geo-Inf. 2021, 10(6), 360; https://0-doi-org.brum.beds.ac.uk/10.3390/ijgi10060360
Submission received: 2 April 2021 / Revised: 14 May 2021 / Accepted: 20 May 2021 / Published: 25 May 2021

Abstract

:
The aqueducts built by the ancient Romans are among the most impressive evidence of their engineering skills. The water inside the aqueducts was transported for kilometers, exploiting only the slight but constant differences in altitude throughout the route. To keep the differences in height constant, the aqueducts could proceed underground or aboveground on well-known arched structures that supported lead, ceramic or stone pipes. In order to reconstruct the characteristics of these structures, it is necessary to carry out an accurate survey of the orthometric heights, and therefore the most suitable technology is geometric levelling. In this case, however, it is not applicable, and therefore here we propose an alternative methodology. The final goal of this work was to estimate the flow of some sectors of these aqueducts preserved in the area south of the city of Rome. This has two main purposes: The first is to reconstruct the flow rate of these aqueducts for historical studies; the second is to check how much the orthometric heights have changed over the centuries, in order to reconstruct the movements from a geophysical and geodynamic point of view. The latter analysis will be developed in a following phase of this research. For this purpose, a high-precision geomatic survey was carried out in the area under study, partly retracing a survey already carried out in 1917 whose purpose and methodologies are not known. The area has been affected by a gradual subsidence over centuries, including since 1917. The observed sections of the aqueducts showed average inclinations, slightly lower than the 2 per thousand that is reported in the literature for similar aqueducts. The measurements carried out allowed the flow rate of the two specific aqueducts to be estimated more accurately, both as they were originally and in the presence of deposits that have accumulated during the years of use of the aqueducts. The reconstruction of the initial geometry will later be used as a reference to estimate how much the geodynamic deformations of the area have deformed the aqueducts themselves.

1. Introduction

The aqueducts built by the ancient Romans are generally considered their most original and impressive engineering achievements. Rome was founded at its location due to the presence of a river island (Tiber island), which made it possible to easily cross the Tiber river, which the Romans had a special emotional bond with. The ancient Romans managed to satisfy their water needs for centuries thanks to the river, and the springs and wells around it. In 312 BC, however, with the rapid development of the city, more water became necessary. This was due both to the increase in the number of inhabitants, but also to the increasing spread of the thermal baths that characterized the entire history of the Roman Empire [1]. For this reason, they started building aqueducts to bring to Rome water from all the neighboring regions (Figure 1), which was then used not only for the thermal baths but also for fountains and private houses. The water in the aqueducts was moved only by a slight but constant difference in altitude, which created pressure differences which were preserved through the use of ceramic, stone, and lead pipes [2].
The use of the latter for drinking water was widespread, even though the Romans were already aware of the health risk associated with the use of lead pipes for drinking water, and the possible occurrence of related diseases [3]. This is why, in reality, most of the aqueducts had stone and ceramic pipes, despite the greater technical difficulty of keeping the pipes watertight, and the pipes we examined were also all made of stone. To maintain a slight but constant decrease in the hilly morphology around Rome, the aqueducts proceeded underground (Figure 2) and along elevated stretches (Figure 3) on well-known arched structures. In some cases, the aqueducts were supplemented by sedimentation and distribution tanks, as well as sluice gates, to correctly regulate the flow according to specific local requirements.
In the present work, the path of some existing stretches of Roman aqueducts south of the city was studied. For this purpose, high precision levelling was carried out in the south-east area of the present city, repeating, in part, a survey carried out for the same purpose almost exactly one century ago (1917) by the Faculty of Engineering of the University of Rome [5].
The study of the characteristics of these engineering works requires the use of the techniques and methods of the historical metric survey, which has been described in detail in the literature [6,7].
The study for the survey and conservative documentation of monuments, buildings, and historical-archaeological remains is the subject of research and investigation by researchers all over the world. Generally, the aim is to document the remains in detail in order to be able to correctly plan their restoration or make it possible to visit them virtually.
Rarer is the survey for reverse engineering, i.e., to reconstruct the functioning of an instrument [8] or an infrastructure as in the case of this study.
Surveys for documentation purposes have the main objective of detecting in an efficient manner a large number of three-dimensional points that characterize the shape of the object. For this type of application, the most widespread methodology for many years has been terrestrial laser scanning, which allows reconstructing accurate point clouds in a relatively short time with high efficiency [9]. The only disadvantage of this technique is that it is usually not possible to acquire the whole object from a single position and the alignment of the various acquisitions can be time consuming and complex, which has stimulated the search for smart techniques to automate these processes [10]. Another technique often used for this type of survey is terrestrial photogrammetry [11], which has recently benefited from the increasing application of algorithms already known in other application fields such as structure from motion (SFM) [12]. SFMs, which often exploit acquisition redundancy to perform camera self-calibration [13], have made it possible to develop new techniques for indoor surveys [14] and, above all, to make the most of drone surveys [15,16,17].
On the other hand, in the case in which the survey of archaeological sites needs reverse engineering, it is necessary to reconstruct the three-dimensional coordinates of a number of points generally much lower than in the previous case, but with even greater three-dimensional accuracy [8].
In the specific case of the survey of aqueducts, the coordinate of greatest interest is obviously the height or its differences along the route. In particular, the height of some notable points that can characterize the height variations of the lowest part of the pipeline is of interest. This survey, in the design phase, can currently be carried out with great accuracy by means of geometric levelling. Similar methods were probably also used in ancient times, but unfortunately, they require access to the lower part of the conduit and this, once the aqueduct was completed, is generally no longer possible for obvious reasons of accessibility. Therefore, it is possible to use the already illustrated techniques of laser scanning or drone photogrammetry, but the accuracy is not entirely adequate for the specific investigation. In the case under study, advantage was taken of the numerous interruptions in the aqueduct itself, which made the lowest part of the conduits observable (although not stationable), allowing them to be surveyed using a total station.
The collapses that caused the interruptions to the aqueduct could be due to crustal movements affecting the area [18].
This paper aimed to evaluate the flow rate and to compare it with historical values, as well as to some results in the literature, including interpretations and experimental conclusions that can be drawn mainly with two main objectives:
  • To reconstruct the differences in historical elevation (comparing the recent survey with the original construction and the survey of 100 years ago), considering that the area is affected by crustal movements;
  • to investigate the geometry and the flow rates (also to understand if the aqueducts worked under pressure) of the original aqueducts.
This reconstruction work is part of a wider scientific project which is still in development, and whose main aim is to verify if and how much the area has been deformed both in the 2000 years since the construction of the aqueducts, and since the survey was carried out a century ago. To be able to hypothesize what the original geometry of the aqueducts was, and therefore estimate the deformations, the flow rate data is obviously of fundamental importance. The most accurate possible knowledge of the flow rates is also of historical importance, as a number of scholars are debating whether the currently known aqueducts were sufficient for the considerable water needs of Imperial Rome or how they intersected and combined with each other [4,19,20].
Where possible, the same points as used in the 1917 campaign were remeasured, considering the difference in height altimetric datum [21,22,23,24].
The two aqueducts we studied, built in the same period, followed their path one above the other taking the supplying of their waters in the distant Aniene valley so their complete original length is about 70 km, while their cross-sectional area is approximately 2 square meters. However, their geometry will be described in more detail in the paragraph on hydraulic modelling.
Here, we follow a path that, from the most extreme part of the current city (Roma Capannelle), arrives in the most urbanized area called ‘Porta Maggiore’. Specifically, the measurable points of the ‘Aqua Claudia’ end in the ‘Mura Latine’ area, while for the ‘Anio Novus’ aqueduct they end in the so-called ‘Porta Furba’ area. To correctly estimate the gradients, the aqueduct was divided into sections that were representative of the original trend, since the subsidence that has already been preliminarily studied in a first phase of this research project [2] reversed the gradient in other parts of the aqueduct, making them unreliable for the reconstruction of the flows.

2. Data Collection and Analysis

The two investigated aqueducts, built in the same historic period, followed the same path with one above the other (Figure 4), taking the supply of their water from the upper part of the Aniene river valley. As already mentioned, the water was only moved by gravity, for this reason, the water had to be taken from sources at a higher altitude than the city, in particular from the area east of the city. From there, the entire route of the aqueduct had to be designed with extreme care to make it compatible with the morphology of the terrain and, where possible, exploit the natural slopes. For these estimations, we used the ellipsoidal height differences measured with GNSS differential receivers and total stations, the area of transit sections of flow rates, reconstructed in previous studies in the literature, and the equations of flow in conditions of uniform motion in channels. Obviously, the ellipsoidal elevations were converted into orthometric elevations according to the procedure described in more detail in the section on surveying operations.
The reconstruction of the gradient profiles, through the measurement of the heights of several points on the bottom of the ‘specus’ along the two aqueducts of interest, allowed us to better study the aqueducts and the vertical movement of the surrounding area. We hypothesized, through the knowledge of the flow rates of the historical runoff, the placement of sections of the ancient channels, making possible a more accurate reconstruction of the ancient structures. It is well known that in order to reconstruct height differences for hydraulic calculations, it is necessary to use a technique capable of measuring orthometric height differences. The most accurate technique for this type of measurement is geometric levelling, with accuracies of up to a tenth of a millimeter on a single levelling section.
Unfortunately, this technique requires accessibility to the points to be measured, which in a historical aqueduct that has many collapsed sections, is not realistically feasible. This consideration required us to look into an alternative technique and raised the question as to which technique was used in 1917 by Reina et al. [5]. From the documentation available to us, it would appear that long geometric levelling was carried out at the ground level, from which a series of points were then ‘launched’ onto the internal structure of the pipelines where the various collapses made them visible. However, this second part described here is not extensively documented in the reports we have, and therefore the methods and relative accuracy are not known for sure.
To compare these surveys, a specific study was also necessary to reconstruct the original height datum, which is different from that presently used and that can be estimated by measuring unchanged benchmarks as they are described in contemporary maps [24], as we will explain in more detail below.

2.1. Geomatic Surveys

All the 1917 recognizable surveyed points were measured again, referring to modern height datum realized by the GNSS network managed by the Lazio regional administration. As it is well known, the heights and elevations calculated using a network based on GNSS permanent stations are ellipsoidal and not orthometric with the well-known physical and numerical differences. The levelling performed in the first part of this research allowed the reconstruction of the gradients of the aqueducts ‘Aqua Claudia’ and ‘Anio Novus’. For both aqueducts, the surveys made it possible to estimate, with sufficient approximation, the flow within the channels. These flows were subsequently compared with those used by the Romans.
To estimate the gradients, the elevations of the lower part of the channels (specus) and the distances along the route obtained from the plan-altimetric survey performed were taken into account.
To compare the elevations measured in 1917 and those of today on the same points, one must consider both the probable variation of the altimetric reference system (datum), and any altimetric modifications due to geodynamic deformations. In this regard, it must be remembered that the altimetric datum used in the course of a century can be profoundly varied. In particular today, in peninsular Italy, an altimetric datum is used which conventionally refers to 1 January 1942 [24,25], which is certainly not the one used in 1917. Unfortunately, from the official documentation of the time [5,26,27], it is not possible to reconstruct with certainty the datum used at that time, nor were useful parameters provided for the conversion between different datums. It was deduced that the only way to compare the two surveys carried out with different reference systems was to try to reconstruct the heights of some benchmarks already in existence at that time, and still preserved today, in order to estimate the relative variation of height. Unfortunately, only two points levelled in 1917 are still recognizable, but their undisturbed maintenance is not certain. In fact, the pillar at ‘Piazza San Giovanni’ is visibly tilted (Figure 5), and it is not certain whether it was also tilted in 1917.
Additionally, the last step at the church of Santa Croce in Gerusalemme is on a staircase that may have been subject to wear over the years (Figure 6).
The differences between the heights in 1917 and those of today, however, gave decidedly discordant values of −38 cm for the first benchmark and −59 cm for the second.
Estimating the datum shift on these two points alone is therefore not statistically significant, and an independent way to verify the estimated datum differences was needed. In order to do this, it was first assumed that since the survey was carried out by the Italian Military Geographical Institute (IGMI), they had used the same altimetric datum that they were using at the time for detailed maps of some urban centers.
A map of the IGMI, at a scale of 1:5000 from 1924, was available from previous research [27,28,29], on which the benchmarks used for its realization were reported. The cartography of 1924 contained both the heights of many benchmarks surveyed with geometric levelling, and therefore with millimetric/centimetric precision, and also a significantly greater number of points surveyed with lower precision. Therefore, it was decided to carry out a comparison between the elevations of the points on the 1924 map and a DEM obtained from the spot points of a modern (2005) map at the same scale, of which the altimetric datum was therefore certain (Genova, 1942).
The estimation was carried out on 139 points extracted from historical cartography which were compared with the DEM obtained from modern vector cartography. A 4-metre filter was imposed to eliminate any outliers due to errors in the interpolation of the DEM from vector mapping and 26 points were thus eliminated. The remaining 113 points (Figure 7) had a minimum difference value of −2.93 and a maximum value of 2.17 m. We were aware that the accuracy of the heights shown on the two maps was certainly not of the order of that obtainable with a geometric levelling, but we assumed that the large number of points compared could provide a statistically more reliable and therefore significant sample.
The difference value estimated in this way was −42 cm. On the one hand, this partly confirms the measurements made on the two benchmarks (although it does not coincide with either of them, but only with their average). On the other hand, it indicates that the altimetric variation between the two altimetric datums (1917 and 1941) is of the decimetric order, and cannot be ignored.
Considering that the uncertainty of the 1917 elevation datum could create the exact same shift at all points measured in 1917, it is clear that comparisons between differences in elevation measured in 1917 and those measured today are more meaningful than the differences between absolute elevations measured in 1917 and those measured today. This has been taken into account in this work, and therefore the trend differences between the two sets of measurements were considered more significant.
In order to acquire the heights of the remaining points, it was not possible (as already mentioned) to use geometric levelling, and therefore it was decided to employ GNSS levelling using dual-frequency geodetic receivers, and to ‘launch’ from the same stationary points the total station measurements of the visible points, effectively measuring trigonometric levelling.
The points were measured using a total station ‘Leica TCR703’ (Figure 8a), and the total station’s positions were surveyed using a Topcon Legacy-E double frequency GNSS receiver (Figure 8b).
The GNSS data were processed in a post-processing mode using RINEX data provided by the network of permanent stations managed by the Lazio region [10] (Figure 9).
The heights obtained from trigonometric levelling calculated from the GNSS heights of the same station points must obviously be considered as calculated with reference to the WGS84 ellipsoid in the same realization of the reference network (RDN2008) [21,23].
An average accuracy of less than 2 cm can be attributed to these heights, mainly due to the difficulty of collimating with certainty the main points of the specus from the ground.
These elevations were converted to orthometric elevations in the national elevation datum using data from the ITALGEO2005 geoid model, which has an estimated vertical accuracy of about 2.5 cm [30] when using the specific interpolation algorithms as implemented in the GeoTrasformer software [31]. We could identify and measure six levelled points from the 1917 measurement campaign: Points belonging to the ‘Marcio’, ‘Tepula’, and ‘Julia’ aqueducts at Porta Maggiore (Figure 10a), a point belonging to the ‘Aqua Claudia’ aqueduct at ‘Mura Latine Street’, and points belonging to the ‘Aqua Claudia’ and ‘Anio Novus’ aqueducts (Figure 10b) located in the Aqueduct Park (Figure 10c). These points are the only ones that were identified and measured with certainty.
In the work by Reina et al. [5], levelling campaigns that were performed on the aqueducts ‘Anio Novus’, ‘Anio Vetus’, ‘Marcia’, and ‘Claudia’ were described. They surveyed three different lines of levelling: From Rome to the ‘Capannelle’ area, from Anagnina street to the area ‘Pallavicina’, and from the ‘Pallavicina’ area to Cineto Romano.
The section of our interest, however, was only the first part of Rome-Capannelle, along which the aqueducts ‘Anio Novus’ and ‘Claudia’ are present (Figure 11 and Figure 12).

2.2. Hydraulics of the Roman Aqueducts

The Romans developed very high-level hydraulic engineering to support the development of the water supply to their main cities, especially Rome. Several aqueducts were realized, and their design was undertaken under consideration of the absence of pumping systems. The flow through them is entirely driven by gravity, flowing slowly from higher elevations. Thus, the aqueduct path was designed to respect the continuous decrease in the topographic level, considering a given gradient to avoid a high flowrate and its possible consequences. The Romans aimed for a slope of 0.02%, as reported in the Natural History by Plinius, who indicated ‘channel must be as solid as possible, and their slope should not be less than twenty centimetres per kilometre of length’. This value of 0.02% ensured the proper flow velocity without requiring a pumping system. The requirement of a solid channel was to guarantee the solidity of the structure, and its reliability during the water flow. The larger part of the Roman structures work at atmospheric pressure due to the difficulty, at that time, of creating piping systems that were solid enough. Roman aqueducts were usually built as buried masonry channels, with a rectangular internal profile and a semi-circular vault. The floor and the sidewalls were covered with a characteristic red waterproof cement, named opus signinum, with a special quarter round structure to strengthen the joints of the walls and floor. The channel was usually only half [33] to two-thirds full of water (Figure 13).

2.2.1. Flow Rate Calculation Method

The water flow into the aqueduct was assumed to be an open flow with a free surface. This assumption seems to be confirmed by Blackman, who cites ‘the language used by Frontinus and others implies a free surface’.
Following these assumptions, the flow rate was calculated with equations of flow in a uniform steady-state open channel flow.
The most common are Chézy’s equation and Manning’s equations [34,35]. To calculate the aqueduct discharge, the Manning equation and the cross-sectional area are used. Thus, the flow rate is calculated by the following formula:
Q = A 1 n R h 2 3 S
where Rh is the hydraulic radius (the ratio between the area and the wetted perimeter), A is the flow area, S is the hydraulic gradient or the bottom slope of the channel, and n is the Manning coefficient. To evaluate the capacity of the Aqua Claudia and the Aqua Anio Novus, the cross-sectional area characteristics, the slope of the channel, and the roughness of the walls are needed.

2.2.2. Reconstruction of the Outflow Sections and Flow Characteristics

The reconstruction of the duct section is not an easy task for two reasons:
-
The collapse and deformation of the masonry parts;
-
the formation of deposits of sediment inside the ducts themselves. The first indication of the cross section characteristics were reported by Ashby [36]. He undertook some surveys following the ones done by Reina et al. [5]. In the investigated stretch, four sections were identified on the Aqua Claudia and only one on the Aqua Anio Novus (Table 1).
Blackman [37] analyzed only the width of the specus, giving a distribution of all measurements. The mean values were 1.1 m for Aqua Claudia and 1.15 m for the Aqua Anio Novus. The approach followed by Blackman was to search for the depth of water having assumed a rectangular cross section.
More recent studies give some evidence of the real cross-section of the aqueducts. Keenan-Jones [38] and Motta et al. [39] highlighted the presence of travertine concretions on the bottom and walls of the Anio Novus aqueduct. They studied sketches near Tivoli and in Rome, near the “Casale Roma Vecchia” area, and focused on the travertine concretions to undertake a hydraulic evaluation of the discharge through the Anio Novus aqueduct.
For the present study, the reconstruction of the Anio Novus section near point I.4 from Reina et al. [5] was considered. The reconstruction of the section based on the caractheristics reported by Keenan-Jones et al. [38] is reported in Figure 14.
For the Aqua Claudia, the cross-section was derived from the available documentation, and its shape was determined to be rectangular due to the superposition of the Anio Novus (the roof of Claudia is the floor of the Anio Novus). In Figure 15 the geometry is shown. No information has been found regarding a travertine presence in the specus of Aqua Claudia.
The flow characteristic to be defined is the Manning coefficient. It depends on the roughness of the surfaces in contact with water. In case the water flowed in contact with the mortar lining, the estimate of the roughness coefficient was 0.014 s/m1/3, as suggested by Fahlbusch [40]. In other studies, lower values have been proposed.
The presence of travertine at the bottom of the specus was also well documented in the past. It was possible to survey the surface in contact with the flowing water, and its roughness was determined. Figure 16 shows the thickness of the travertine deposits and the surface created on the specus, such as a riverbed with cobbles of nearly 15 cm in diameter. The appearance of the surface suggests a Manning roughness coefficient around 0.035 s/m1/3. The actual surface did not respond to the status during ancient Rome, and was probably smoother due to the presence of a muddy layer, a roughness coefficient in agreement with Fahlbusch’s estimation (1982) of 0.019 s/m1/3 was used.

3. Results and Discussion

The first result of the study is the slope analysis of the aqueducts. In fact, to find the hydraulic gradient under the hypothesis of uniform open-channel flow, the slope of the aqueduct floor is required. To reconstruct the bottom profile of the two aqueducts, rather than using the well-known levels from Reina et al. [5], the data collected during the geomatics survey were used. The main objective was to verify the acquired level in 1917, and use them to evaluate the slope of the bottom of the channels of the ancient aqueducts.
The two aqueducts, as presented in Figure 17, follow the same path with one above the other, maintaining a distance between the channels that is almost constant with an average of 3 m and a standard deviation of 0.26 m. Ashby [36] reports that the distance between the bottom of the Anio Novus and the intrados of the Claudia is about 1.00 m.
Considering the height of the Claudia of 1.75 m and the thickness of the slab of 1 m, the distance between the two specus is 2.75 m. In our data, the values are always greater. In the three points taken as reference from the Reina et al. (1915) levelling survey (II.3); the distance between the specus is 2.85 m for the points I.1 and I.2 and 1.7 m for point I.3. It is clear that there is a mistake in the reported data. Comparing the level measurements of the bottom of the Claudia with the present survey, the difference of 1 m confirms the potential mistake.
In Figure 17, the reconstructed data, indicated with x, are reported as well as the connecting lines with dashed ones, having assumed a distance between the specus of 3 m. The point highlighted with a circle corresponds to an outlier of this distance. In the stretch, a high variability of the level and also inverse slopes appear. Those behaviors are probably due to ground or structural movements.
In fact, in some parts, due to the subsidence of the structures, the trend of the profile does not respect the hydraulic requirements for the outflow of waters. To estimate the gradients, it was decided to divide the route into sections, for which a positive slope of the aqueducts is guaranteed. In Table 2, the calculated slopes are given, along with the average between the positive values according to our calculations, and the average gradient between the extreme points.
Following the suggestion from Keenan-Jones on the incorrect positioning of benchmark 37, corresponding to point I.4, a verification of the more probable position was done. From the representation of Figure 18 and Figure 19 in Reina et al., is clear that point I.4 is about 350 m downstream from the indicated place. In this case, the slope of the two aqueducts over the entire stretch is 1.68‰ for Aqua Claudia and 1.64‰ for Aqua Anio Novus.
From historical documents and technical studies, is assumed that the ancient Romans made an effort to achieve an almost constant slope of 2‰. This slope could vary slightly concerning the routes that the aqueducts were to follow, the obstacles to be overcome, and the flow velocity desired.
The calculated average slope for the Anio Novus aqueduct in the analyzed section is 1.31‰, and for the Aqua Claudia aqueduct is 1.63‰. Those values agree with values previously cited in historical documents.

Evaluation of the Flow Capacity of Roman Aqueducts

To evaluate the volume capacity of the Claudia and Anio Novus aqueducts, two different calculations have been adopted. First, using the historical flow rate and the characterization of the aqueduct, the depth of water in the specific section is evaluated. The reference section lies in the stretch of the ‘Parco degli Acquedotti’ as previously reported.
The second calculation is the direct evaluation of the discharge flow rate based on the measured slopes and the sections in different flowing conditions.
The first stage is similar to the one followed by Blackman [37]. The basic assumption was the uniform channel flow, and the section was previously presented. Two different conditions were considered: A fully available section, and a restricted section due to travertine deposits.
From the calculation reported in Table 3, the water depth in Anio Novus is always greater than the maximum height of the specus (1.79 m) in the case of presence of travertine deposits. Without travertine deposits the conditions are in an open channel flow. On the contrary, Aqua Claudia always flows in an open channel condition. Of note, the water depth for the full section was in line with two-thirds of the height of the specus, as reported in the literature.
The second approach is based on direct calculation of the flow rate under the assumption of uniform steady-state open channel flow. For the Manning coefficient, different values have been assumed to have also a sensitivity analysis of discharge flow. The Mannings coefficients used were: 0.14 for the full section while for the case with travertine, 0.19 for low and 0.035 for high. The extreme value of 0.035 has been adopted only as reference for high roughness surfaces considering the limitation in the flow rate. The calculated discharge flow rates are reported in Table 4 for the three different conditions just reported. The assumption in the calculation of flow rates are: Maximum flow rate corresponding to the full section and the standard flow rate with the maximum water depth of 2/3 of the height section.
The presence of the travertine deposits can highly reduce the capacity of discharge of the aqueducts in the range of 30% to 68%. This was well known by Romans that were doing maintenance work by removing the scaling from the wall and from the basement.
The discharge capacity of standard flow rate highlights the requirements of higher water depth in the analyzed sections to accomplish the historical flow rate. This is not mandatory, but it gives an idea of the possible working range of the aqueducts.
The Anio Novus was designed without clearance for historical flow rates. The section restricted by travertine deposits has a reduced capacity of 1.27 m3/s (i.e., 109,454 m3/d). Thus, it will flow in pressure and to avoid pressure inside the channel, the flow should be less.

4. Conclusions

The techniques used here made it possible to reconstruct an orthometric height difference with an accuracy prudentially limited to more or less 4 cm, allowing the geometric characteristics of two Roman aqueducts to be reproduced for the observable sections. Since it was not possible to reconstruct the difference in the altimetric datum between the 1917 elevations and the current reference system with absolute certainty, comparisons of height differences measured in the two surveys must be considered more significant than the absolute differences in elevation. The comparison between the current survey and that of 1917 showed a difference of the orthometric heights datum of about 42 cm due to the different altimetric datum used in 1917 and presently (Genova, 1942).
In particular, the ‘Anio Novus’ was surveyed and reconstructed from the area of ‘Mandrione’ street to the area of the ‘Aqueduct Park’, where the aqueduct remains almost intact. The ‘Aqua Claudia’ aqueduct was taken over and rebuilt from the area of Via Latina up to the same area of the ‘Aqueduct Park’, where the survey of the other aqueduct was completed. In both cases, the average gradient was close to 2%, which is the value reported in the historical reference texts for the design of Roman hydraulic works. Specifically, ‘Anio Novus’ currently has a slope of 1.31‰, while ‘Aqua Claudia’ has a slope of 1.63‰. These slope values made it possible to estimate the flow of water, and consequently to verify the conditions of the channels in which the water flowed. To reconstruct the cross-section of both aqueducts, the study by Keenan-Jones et al. [39] and the study by Ashby [36] were used. Based on the comparison between the calculated flow rate and water heights and the data reported in the literature, two different conclusions can be reached for the two aqueducts. The Aqua Claudia has a calculated discharge flow rate higher than expected, and therefore the channel worked with a 70–95% filling of the section, as reported in the literature. The Anio Novus aqueduct, on the other hand, currently has a discharge flow rate lower than historically, suggesting that in the study area it was working under pressure at the time.
This work has highlighted the possibility to investigate the hydraulic of Roman aqueducts starting from a geodetic survey. The importance of the survey is clear when comparing the historical documents and data. The more difficult aspect is to understand after centuries which is the original path and relative slopes measured on ruins that can be exposed to ground movements.
A better knowledge of the hydraulic characteristics of the aqueduct will help in the subsequent phases of this research project where the current strain of the structure will allow confirming or refuting some theories on the deformations taking place in this area.

Author Contributions

The authors contributed equally to this work. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

The data presented in this study are available on request from the corresponding author. The data are not publicly available for privacy reasons.

Acknowledgments

The authors would like to thank all the people who helped carry out the surveys by allowing access to numerous private areas, and the facilities that enabled the research work to be carried out.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Alimonti, C.; Baiocchi, V.; Bonanotte, G.; Molnar, G. Measuring and leveling roman aqueducts to estimate their flows. In Proceedings of the 2018 Metrology for Archaeology and Cultural Heritage (MetroArchaeo), Cassino, Italy, 22–24 October 2018; pp. 527–531. [Google Scholar] [CrossRef]
  2. Bonanotte, G. Measuring the Recent Vertical Crustal deformation with Leveling of Old Roman Aqueducts. Master’s Thesis, Sapienza University of Rome, Rome, Italy, 2015. [Google Scholar]
  3. Vitruvius, M.P. De Architectura; Einaudi: Milano, Italy, 1997; VIII.6.10-11. [Google Scholar]
  4. Blanco, D.; Del Grande, P. Via Lucio Mariani (Gregna Sant’Andrea). Resti di acquedotti e tracciati viari d’epoca romana (Municipio X). Bull. Della Comm. Archeol. Comunale Roma 2016, 117, 380–385. [Google Scholar] [CrossRef]
  5. Reina, V.; Corbellini, G.; Ducci, G. Livellazione Degli Antichi Acquedotti Romani, Tip; Della Reale Accademia dei Lincei: Roma, Italy, 1917. [Google Scholar]
  6. Bitelli, G.; Balletti, C.; Brumana, R.; Barazzetti, L.; D’Urso, M.G.; Rinaudo, F.; Tucci, G. Metric documentation of cultural heritage: Research directions from the italian gamher project. ISPRS Int. Arch. Photogramm. Remote. Sens. Spat. Inf. Sci. 2017, XLII-2/W5, 83–90. [Google Scholar] [CrossRef] [Green Version]
  7. Balletti, C.; Guerra, F. The survey of cultural heritage: A long story. Rend. Lince 2015, 26, 115–125. [Google Scholar] [CrossRef]
  8. Baiocchi, V.; Brigante, R.; Del Pizzo, S.; Giannone, F.; Onori, M.; Radicioni, F.; Stoppini, A.; Tosi, G.; Troisi, S.; Baumgartner, M. Integrated Geomatic Techniques for Georeferencing and Reconstructing the Position of Underground Archaeological Sites: The Case Study of the Augustus Sundial (Rome). Remote. Sens. 2020, 12, 4064. [Google Scholar] [CrossRef]
  9. Costantino, D.; Angelini, M.G.; Caprino, G. Laser Scanner Survey of an Archaeological Site: Scala Di Furno (Lecce, Italy). Int. Arch. Photogramm. Remote Sens. Spat. Inf. Sci. 2010, 38, 178–183. [Google Scholar]
  10. Guarnieri, A.; Vettore, A.; Camarda, M.; Costantino, D. Automatic registration of large range datasets with spin-images. J. Cult. Herit. 2011, 12, 476–484. [Google Scholar] [CrossRef]
  11. Fazio, L.; Brutto, M.L.; Dardanelli, G. Survey and virtual reconstruction of ancient roman floors in an archaeological context. ISPRS Int. Arch. Photogramm. Remote. Sens. Spat. Inf. Sci. 2019, XLII-2/W11, 511–518. [Google Scholar] [CrossRef] [Green Version]
  12. Ullman, S. The interpretation of structure from motion. Proc. R. Soc. Lond. B 1979, 203, 405–426. [Google Scholar]
  13. Fraser, C.S. Digital camera self-calibration. ISPRS J. Photogramm. Remote Sens. 1997, 52, 149–159. [Google Scholar] [CrossRef]
  14. Masiero, A.; Fissore, F.; Guarnieri, A.; Pirotti, F.; Visintini, D.; Vettore, A. Performance evaluation of two indoor mapping systems: Low-cost UWB-aided photogrammetry and backpack laser scanning. Appl. Sci. 2018, 8, 416. [Google Scholar] [CrossRef] [Green Version]
  15. Ebolese, D.; Lo Brutto, M.; Dardanelli, G. Uav survey for the archaeological map of Lilybaeum (Marsala, Italy). ISPRS Int. Arch. Photogramm. Remote. Sens. Spat. Inf. Sci. 2019, XLII-2/W11, 495–502. [Google Scholar] [CrossRef] [Green Version]
  16. Alessandri, L.; Baiocchi, V.; Del Pizzo, S.; Di Ciaccio, F.; Onori, M.; Rolfo, M.F.; Troisi, S. A flexible and swift approach for 3D image–based survey in a cave. Appl. Geomat. 2020, 1–15. [Google Scholar] [CrossRef]
  17. Pepe, M.; Costantino, D. Techniques, tools, platforms and algorithms in close range photogrammetry in building 3D Model and 2D representation of objects and complex architectures. Comput. Des. Appl. 2020, 18, 42–65. [Google Scholar] [CrossRef]
  18. Riguzzi, F.; Atzori, S.; Obrizzo, F.; Pietrantonio, G.; Pingue, F. Recent terrestrial and satellite observations to model deformations at Colli Albani volcano (central Italy). Int. J. Earth Sci. 2012, 101, 1661–1671. [Google Scholar] [CrossRef]
  19. Deming, D. The Aqueducts and Water Supply of Ancient Rome. Ground Water 2020, 58, 152–161. [Google Scholar] [CrossRef] [Green Version]
  20. Hodge, A.T. Aqueducts and Water Supply. In A Companion to the Archaeology of the Roman Republic; Wiley & Blackwell: Hoboken, NJ, USA, 2013; pp. 285–295. [Google Scholar]
  21. Regione Lazio GNSS Web Page. Available online: http://62.149.194.137/SpiderWeb/frmIndex.aspx (accessed on 8 May 2021).
  22. EPSG Web Site. Available online: https://epsg.io/6708 (accessed on 8 May 2021).
  23. Barbarella, M. Digital technology and geodetic infrastructures in Italian cartography. Citta Stor. 2014, 9, 91–110. [Google Scholar]
  24. Surace, L. I sistemi di riferimento geotopocatografici in Italia. Boll. Geod. Sci. Affin. 1996, 57, 181–234. (In Italian) [Google Scholar]
  25. Pierozzi, M. Il Sistema Altimetrico Italiano, la liveLlazione, lo Zero Idrografico ed i Riflessi in Ambito Portuale. 2016. Available online: http://www.assoporti.it/sites/www.assoporti.it/files/eventiEsterni/Pierozzi.pdf (accessed on 8 May 2021).
  26. Mori, A. La cartografia ufficiale in Italia e l’Istituto Geografico Militare Nel cinquantenario Dell’istituto Geografico Militare (1872–1922); Istituto Geografico Militare, Stabilimento Poligrafico per l’Amministrazione della Guerra: Roma, Italy, 1922; 425p. [Google Scholar]
  27. Baiocchi, V.; Lelo, K. Assessing the accuracy of historical maps of cities: Methods and problems. Citta Stor. 2014, 9, 61–89. [Google Scholar]
  28. Loperfido, A. Survey for the Topographical Plan of Rome; Alfani e Venturi: Firenze, Italy, 1908. (In Italian) [Google Scholar]
  29. Baiocchi, V.; Lelo, K. Accuracy of 1908 high to medium scale cartography of Rome and its surroundings and related georeferencing problems. Acta Geod. Geophys. Hung. 2010, 45, 97–104. [Google Scholar] [CrossRef]
  30. Barzaghi, R.; Borghi, A.; Carrion, D.; Sona, G. Refining the estimate of the Italian quasigeoid. Boll. Geod. Sci. Affin. 2007, 3, 145–160. [Google Scholar]
  31. Baiocchi, V.; Camuccio, P.; Zagari, M.; Ceglia, A.; Del Gobbo, S.; Purri, F.; Cipollini, L.; Vatore, F. Development of a geographic database of a district area in open source environment Geoingegneria Ambientale e. Mineraria 2017, 151, 97–101. [Google Scholar]
  32. Baiocchi, V.; Alimonti, C.; Bonanotte, G.; Molnar, G. Geomatic measurement of “New Aniene” and “Claudia” roman aqueducts for flows estimation. IOP Conf. Ser. Mater. Sci. Eng. 2020, 949. [Google Scholar] [CrossRef]
  33. De Feo, G.; Angelakis, A.N.; Antoniou, G.P.; El-Gohary, F.; Haut, B.; Passchier, C.W.; Zheng, X.Y. Historical and Technical Notes on Aqueducts from Prehistoric to Medieval Times. Water 2013, 5, 1996–2025. [Google Scholar] [CrossRef] [Green Version]
  34. Chow, V.T. Open–Channel Hydraulics; McGraw–Hill: New York, NY, USA, 1959. [Google Scholar]
  35. White, F.M. Fluid Mechanics, 2nd ed.; McGraw–Hill: New York, NY, USA, 1986. [Google Scholar]
  36. Ashby, T. The Aqueducts of Ancient Rome; Clarendon Press: Rome, Italy, 1935. [Google Scholar]
  37. Blackman, D.R. The Volume of Water delivered by the Four Great Aqueducts of Rome. Pap. Br. Sch. Rome 1978, 46, 52–72. [Google Scholar] [CrossRef]
  38. Keenan-Jones, D.; Motta, D.; Garcia, M.H.; Fouke, B.W. Travertine-Based Estimates of the Amount of Water Supplied by Ancient Rome’s Anio Novus Aqueduct Journal of Archaeological Science: Reports; Elsevier: Amsterdam, The Netherlands, 2015. [Google Scholar]
  39. Motta, D.; Keenan-Jones, D.; Garcia, M.H.; Fouke, B.W. Hydraulic Evaluation of the Design and Operation of Ancient Rome’s Anio Novus Aqueduct. Archaeometry 2017, 59, 1150–1174. [Google Scholar] [CrossRef]
  40. Crouch, D.P.; Fahlbusch, H. Vergleich antiker griechischer und römischer Wasserversorgungsanlagen. Am. J. Archaeol. 1986, 90, 238. [Google Scholar] [CrossRef]
  41. Di Fenizio, C. Sulla Portata Degli Antichi Acquedotti Romani e Determinazione Della Quinaria Giornale del Genio Civile; Anno LIV; Luglio: Roma, Italy, 1916. [Google Scholar]
Figure 1. Main Roman aqueducts, and the two rivers (light blue), in the bold black circle are the stretches studied here.
Figure 1. Main Roman aqueducts, and the two rivers (light blue), in the bold black circle are the stretches studied here.
Ijgi 10 00360 g001
Figure 2. General view (a) and interior (b) of the Anio Novus specus (water channel) Reprinted with permission from ref. [4]. Copyright 2016, D. Blanco.
Figure 2. General view (a) and interior (b) of the Anio Novus specus (water channel) Reprinted with permission from ref. [4]. Copyright 2016, D. Blanco.
Ijgi 10 00360 g002
Figure 3. A section of an arched Roman aqueduct during the survey.
Figure 3. A section of an arched Roman aqueduct during the survey.
Ijgi 10 00360 g003
Figure 4. Anio Novus on top of the Aqua Claudia.
Figure 4. Anio Novus on top of the Aqua Claudia.
Ijgi 10 00360 g004
Figure 5. San Giovanni in Laterano square 1917 benchmark.
Figure 5. San Giovanni in Laterano square 1917 benchmark.
Ijgi 10 00360 g005
Figure 6. Santa Croce in Gerusalemme church 1917 (a) benchmark detail (b).
Figure 6. Santa Croce in Gerusalemme church 1917 (a) benchmark detail (b).
Ijgi 10 00360 g006
Figure 7. The 1924 map benchmarks on a 2005 map DEM.
Figure 7. The 1924 map benchmarks on a 2005 map DEM.
Ijgi 10 00360 g007
Figure 8. Survey phases using forced centerings with the total station (a) and GNSS antenna (b).
Figure 8. Survey phases using forced centerings with the total station (a) and GNSS antenna (b).
Ijgi 10 00360 g008
Figure 9. CORS network managed by Lazio regional administration.
Figure 9. CORS network managed by Lazio regional administration.
Ijgi 10 00360 g009
Figure 10. The Marcio, Tepula, and Julia aqueducts at Porta Maggiore square (a) and at Mura Latine street (b). The Aqua Claudia and Anio Novus aqueducts in the Aqueduct Park (c) Reprinted with permission from ref. [32]. Copyright 2020, V. Baiocchi.
Figure 10. The Marcio, Tepula, and Julia aqueducts at Porta Maggiore square (a) and at Mura Latine street (b). The Aqua Claudia and Anio Novus aqueducts in the Aqueduct Park (c) Reprinted with permission from ref. [32]. Copyright 2020, V. Baiocchi.
Ijgi 10 00360 g010
Figure 11. Points used to calculate the height difference along the entire route.
Figure 11. Points used to calculate the height difference along the entire route.
Ijgi 10 00360 g011
Figure 12. Points used to calculate the height difference in the area known as the “Parco degli Acquedotti”.
Figure 12. Points used to calculate the height difference in the area known as the “Parco degli Acquedotti”.
Ijgi 10 00360 g012
Figure 13. Typical cross-section of a Roman aqueduct channel. The channel is built of masonry and buried in the ground. The inside is covered with a red waterproof cement (opus signinum—in red), with characteristic quarter rounds to seal the edges. Carbonate deposits (gray) are common in some channels.
Figure 13. Typical cross-section of a Roman aqueduct channel. The channel is built of masonry and buried in the ground. The inside is covered with a red waterproof cement (opus signinum—in red), with characteristic quarter rounds to seal the edges. Carbonate deposits (gray) are common in some channels.
Ijgi 10 00360 g013
Figure 14. Cross-section geometry for the Anio Novus aqueducts. (a) Full section (b) with travertine deposits.
Figure 14. Cross-section geometry for the Anio Novus aqueducts. (a) Full section (b) with travertine deposits.
Ijgi 10 00360 g014
Figure 15. Cross-section geometry for the Aqua Claudia: Point II.3. (a) Full section (b) with travertine deposits (From: Ashby, 1935).
Figure 15. Cross-section geometry for the Aqua Claudia: Point II.3. (a) Full section (b) with travertine deposits (From: Ashby, 1935).
Ijgi 10 00360 g015
Figure 16. (a) A view of the Anio Novus aqueduct’s travertine deposits on the specus. (b) View from the top of the travertine deposits and their roughness.
Figure 16. (a) A view of the Anio Novus aqueduct’s travertine deposits on the specus. (b) View from the top of the travertine deposits and their roughness.
Ijgi 10 00360 g016
Figure 17. Specus level for the two aqueducts: On the x-axis the progressive distance along the aqueduct; on the y-axis the orthometric heights from levelling.
Figure 17. Specus level for the two aqueducts: On the x-axis the progressive distance along the aqueduct; on the y-axis the orthometric heights from levelling.
Ijgi 10 00360 g017
Figure 18. Slope profile of Aqua Anio Novus.
Figure 18. Slope profile of Aqua Anio Novus.
Ijgi 10 00360 g018
Figure 19. Slope profile of Aqua Claudia.
Figure 19. Slope profile of Aqua Claudia.
Ijgi 10 00360 g019
Table 1. Aqua Claudia (II) and Aqua Anio Novus (I) dimensions and levels (Ashby, 1935).
Table 1. Aqua Claudia (II) and Aqua Anio Novus (I) dimensions and levels (Ashby, 1935).
SiteTypeWidth
(m)
Height
(m)
Specus Level
(Orthometric Height, m)
I.4Vault1.141.2577.66
II.3Flat1.141.7565.76
II.4Flat1.141.7566.23
II.5Flat1.101.7969.78
II.6Flat1.101.7975.00
Table 2. Aqua Claudia and Aqua Anio Novus slopes.
Table 2. Aqua Claudia and Aqua Anio Novus slopes.
SlopesAqua Claudia
(‰)
Aqua Anio Novus
(‰)
Average value from data1.662.60
Over the entire stretch1.631.31
From Reina data1.571.54
Table 3. Aqua Claudia and Aqua Anio Novus water depths for historical flow rates.
Table 3. Aqua Claudia and Aqua Anio Novus water depths for historical flow rates.
ReferenceFlow Rate
(m3/d)
Full Section
(m)
With Travertine
(m)
Aqua ClaudiaDi Fenizio [41] (1916)191,1901.2672%1.7098%
De Feo [33] et al. (2013)184,2801.2270%1.6595%
Aqua Anio NovusDi Fenizio (1916)196,6271.5285%2.04115%
De Feo et al. (2013)189,5201.4685%2.05114%
Table 4. Aqua Claudia and Aqua Anio Novus estimated discharge capacity and historical flow rates.
Table 4. Aqua Claudia and Aqua Anio Novus estimated discharge capacity and historical flow rates.
Maximum Flow RateAqua ClaudiaAqua Anio Novus
Full section (m3/d)244,940198,156
With travertine–low (m3/d)166,714109,454
With travertine–high (m3/d)90,50259,418
Standard flow rate
Full section (m3/d)172,864157,709
With travertine–low (m3/d)126,448101,782
With travertine–high (m3/d)68,64355,253
Di Fenizio (1916) (m3/d)191,190196,627
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Alimonti, C.; Baiocchi, V.; Bonanotte, G.; Molnár, G. Roman Aqueduct Flow Estimation Using Geomatic Measurement. ISPRS Int. J. Geo-Inf. 2021, 10, 360. https://0-doi-org.brum.beds.ac.uk/10.3390/ijgi10060360

AMA Style

Alimonti C, Baiocchi V, Bonanotte G, Molnár G. Roman Aqueduct Flow Estimation Using Geomatic Measurement. ISPRS International Journal of Geo-Information. 2021; 10(6):360. https://0-doi-org.brum.beds.ac.uk/10.3390/ijgi10060360

Chicago/Turabian Style

Alimonti, Claudio, Valerio Baiocchi, Giorgia Bonanotte, and Gábor Molnár. 2021. "Roman Aqueduct Flow Estimation Using Geomatic Measurement" ISPRS International Journal of Geo-Information 10, no. 6: 360. https://0-doi-org.brum.beds.ac.uk/10.3390/ijgi10060360

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop