Research Article | Open Access
Kleanthis Simyrdanis, Ian Moffat, Nikos Papadopoulos, Jarrad Kowlessar, Marian Bailey, "3D Mapping of the Submerged Crowie Barge Using Electrical Resistivity Tomography", International Journal of Geophysics, vol. 2018, Article ID 6480565, 11 pages, 2018. https://doi.org/10.1155/2018/6480565
3D Mapping of the Submerged Crowie Barge Using Electrical Resistivity Tomography
This study explores the applicability and effectiveness of electrical resistivity tomography (ERT) as a tool for the high-resolution mapping of submerged and buried shipwrecks in 3D. This approach was trialled through modelling and field studies of Crowie, a paddle steamer barge which sunk at anchor in the Murray River at Morgan, South Australia, in the late 1950s. The mainly metallic structure of the ship is easily recognisable in the ERT data and was mapped in 3D both subaqueously and beneath the sediment-water interface. The innovative and successful use of ERT in this case study demonstrates that 3D ERT can be used for the detailed mapping of submerged cultural material. It will be particularly useful where other geophysical and diver based mapping techniques may be inappropriate due to shallow water depths, poor visibility, or other constraints.
During the last few years there has been an increasing trend of employing the electrical resistivity method in subaqueous areas for geological and environmental studies [1, 2] and underwater geological mapping . A number of methodological studies on the use of ERT in the aquatic environment have been undertaken, including investigating the most appropriate electrode position and array type. Orlando  used numerical simulation modelling to investigate the resolution of underwater resistivity surveys, employing both floating and submerged electrodes. Her results showed that floating cables give low resolution data when the contrast between the resistivity of water and sediment layer is too small. Submarine electrical resistivity data can be collected utilising different electrode arrays such as Dipole-Dipole or Pole-Dipole which are the most widely used ones due to their simple geometry and high data density. In addition, Orlando discusses the Schlumberger array which gives comparable results  and Rucker et al. recommend a Pole-Pole array for the collection of marine resistivity data . Baumgartner and Christensen  describe a particular array where probes are vertically aligned to increase the investigation depth for mapping resistivity anomalies located at greater water depth. Chiang et al. [6, 7] studied the efficiency of a marine towed electrical resistivity method and its sensitivity using 2D numerical simulation modelling as a well-suited method for methane hydrate exploration at a shallow depth beneath the seafloor.
However the use of electrical resistivity method in submarine archaeology prospection studies is not common. Ranieri et al. (: 11) studied the development and application of an integrated methodology for the 3D rendering of geoelectrical data of buried and submerged archaeological features, in complex environments such as coastal areas. A comprehensive feasibility study was also undertaken, investigating the efficiency of ERT in reconstructing submerged archaeological relics (e.g., building foundations) in shallow seawater environment through numerical modelling and field experimentation . The study showed that floating and submerged survey modes can be used with equal success in cases of relatively shallow marine environments when the water depth does not exceed 1 m. In deeper marine environments, the submerged mode survey is recommended for outlining isolated targets. Valid a priori information, in terms of the seawater resistivity and thickness, is important and can greatly improve the inversion results for the data captured with the submerged ERT mode. On the other hand, erroneous information can cause severe distortions in the inversion ERT models and misleading interpretations. Passaro et al. used multibeam bathymetry, seismic, geoelectric, and magnetic methods to investigate archaeological targets protruding from very shallow water . Passaro (: 3) also used electrical resistivity survey to detect a shipwreck in Salerno, southern Italy.
Despite this interest, ERT so far appears not to have been used to map a shipwreck in detail. ERT is particularly useful in cases when the water is too turbid for visual surveys and the water depth is too shallow for acoustic techniques. Furthermore, since many ships are metal-based structured, they become easily detectable with ERT when the target is located in an electrically high resistivity environment such as sand.
2. Crowie Shipwreck
The present study undertook a high-resolution survey of the barge Crowie, which is submerged in the Murray River (South Australia). Crowie was a commercial cargo barge, built in 1911 by Captain J. G. Arnold. It is claimed to be one of the largest barges ever built for the Murray River (: 63; in : 137), and its enormous size is recorded in historical accounts as causing several accidents when traveling up the river. Crowie was eventually superseded by other transport methods and was abandoned near the Morgan wharf where it sank at its mooring, sometime between 1946 and 1950 (: 137) (Figure 1).
The dimensions of Crowie and its structural composition are well recorded in historical records as well as from previous multibeam and side-scan sonar surveys (summarised by ) (Figure 2). The barge is 45 m long, with a 9 m beam and depth from the base of the hold to the deck of 2.5 m (: 140). It was built using a bottom-based construction technique. This led to a flat bottom and straight sides amenable to carry the maximum possible amount of cargo (: 141). The vessel had one deck, a carvel built sharp stern, straight stem, and a composite construction consisting of iron frames and topsides overlain by timber. It had nine internal bulkheads (: 139).
The vessel sits upright in a known location immediately upstream of the Morgan wharf and downstream of the wreck of the PS Corowa. The vessel slightly protrudes from the water during low water periods and extends to a depth of several meters. The extremely turbid water of the Murray makes a visual assessment of the wreck impossible; however geophysical data and inspection of the wreck by touch suggest it is largely intact below the deck-level for its entire length (: 139). The frames of the vessel are intact and easily identifiable. Any iron topsides that were on the vessel at the time of wrecking have corroded away; however metal below a depth of 15 cm below the waterline is believed to be largely intact (: 139).
The large amount of metal used in the construction of the Crowie vessel is likely to give it a highly conductive electrical signature, which should contrast significantly the fresh water and sand which host the wreck. For that reason, electrical resistivity tomography can be considered an appropriate method to apply to this survey, in order not only to detect the specific “target” and define its dimensions (as the previous methods successfully accomplished) but also to estimate its burial depth and confirm the nature of the material that the barge is composed of.
Initially the viability of the proposed field survey was validated using data produced from a three-dimensional numerical simulation model. The ERT numerical experiments were performed using 3D forward and inversion algorithms within the “Res3DMod” and “Res3DInvx64” software . A finite-element modelling subroutine was used to calculate the apparent resistivity values, and a nonlinear smoothness-constrained least-squares optimisation technique was used to calculate the resistivity of the model blocks . The model was constructed with parameters approximating those expected for the field site.
The indicative model used for the numerical simulations included a water layer m thick (under the air half-space) and three homogenous medium layers m, m, and m as shown in Figure 3. The resistivity value of the water layer was set to ohm·m and for the homogeneous layers was set to ohm·m, ohm·m, and ohm·m, respectively. The above resistivity values were estimated based on the local geological map and descriptions of the lithology of units within the Murray Basin [16–18]. A conductive target, which simulated a metallic ship structure (shape was defined according to the multibeam image results in : 141), was embedded in the first layer with dimensions 5.5 m width and 12.5 m length. Internal divisions were also included denoting the different compartments of the upper part of the barge. The resistivity value for the target’s structure was set to ohm·m.
The grid (10 × 15 m) for the numerical simulation was constructed using 672 electrodes (21 lines of 32 electrodes) with spacing m in both directions . The electrodes were submerged in the water bottom 6 m below the water surface. A Pole-Dipole (including the combination of forward and reverse measurements where the electrode B was placed at a theoretically infinite distance) and maximum separation distance (“” unit electrode spacing) was selected and 29618 apparent resistivity measurements were simulated. Data were corrupted with 10% Gaussian noise so that the simulation data were approaching the real noise levels that would be expected in the field. A priori information for constraining the inversion procedure was included, taking into consideration the water depth and the resistivity of the water. The inversion algorithm was terminated after seven iterations unless some other criteria were met (e.g., slow convergence rate of less than 3%, RMS error smaller than the noise level).
The outer structure of the ship was reproduced quite accurately within the dimensions of the constitutive model’s dimensions (Figure 4). The reconstructed resistivity depth sections clearly depict the outer structure of the metal target, with resistivity values close to 0.05–0.2 ohm·m from 0 m to the depth of 1 m below water bottom (depth layers D1–D4). The inner divisions of the barge are not well resolved except for the main axis of the barge which is more obvious at depth slices D3 to D5. Although the target extends no more than 1 m in depth (D4) below the water bottom (according to the model), the inversion result depicts the target (with distorted shape) continuing to 5 m below the water bottom (D5 to D11). The upper background highly resistive layer (0 m to 3.5 m, D1–D10) is represented with a range of gradually increased resistivity values from 150 to 300 ohm·m. The deeper layers (D11-D12) are represented with values close to 1600 ohm·m. The numerical limitations of the modelling and inversion procedures cannot efficiently cope with the magnitude of the resistivity contrast between the conductive target and the resistive background. Hence, the RMS error was high with values reaching almost 30% after 6 iterations.
4. Field Case Study
Following the theoretical modelling of this feature, field survey was undertaken on Crowie. The barge Crowie was submerged near the bank of the Murray River in Morgan (South Australia) (Figure 5). The exact position of the target is described in Roberts et al. .
The study area is located in the Murray Basin, epicratonic sedimentary basin of Paleocene to Quaternary age . This basin contains a sedimentary sequence up to 600 m thick including freshwater, marine, coastal, and continental sediments [17, 18]. At the survey location the stratigraphy includes the Oligocene aged Ettrick Formation, the Early-Middle Miocene aged Mannum Limestone, and a number of Quaternary aged aeolian units (; see Figure 6). The Quaternary aeolian sediments are estimated to have a resistivity value of ohm·m, the Mannum Limestone has a value of ohm·m, and the Marly Ettrick Formation has a value of ohm·m.
4.1. Data Acquisition (Grid, Equipment)
The survey grid was set as a 60 by 15.5 m rectangle approximately parallel to the riverbank as shown in Figure 7. The four corners of the grid were established using heavy rocks as anchors, making sure that they would be unmovable and stable throughout the fieldwork, with taut ropes running vertically from the anchor to a float on the surface. At the extremities of the grid, measuring tapes were floated as guidelines between these anchors for determining the position of each survey line.
The data were acquired using 58 parallel lines (length of each line 15.5 m) oriented perpendicular to the bank and equally spaced ( m apart). The sensors (32 electrodes), equally spaced ( m apart) on each survey line, were submerged on the water bottom on top of the sunken ship perpendicular to its length.
A ZZ FlashRES-Universal resistivity meter was used for injecting and measuring the potential throughout the survey. The water conductivity value and temperature were measured daily with a YSI Pro-1030 hand-held conductivity meter. The average value of conductivity ( ohm·m) was used as a priori information for the inversion procedure. A Pole-Dipole array was used for data collection with a total of #1498 measurements (consisting of 749 forward and 749 reverse measurements). The distance “” between the dipole “MN” varied between 0.5 and 4 m. The distance factor “” between electrode “A” and dipole “MN” varied from 1 to 8.
4.2. Data Acquisition (Bathymetry)
A total station survey was conducted over the wreck site using a Leica TS16 to measure the river depth in the location of the ERT survey. The position of the total station was georectified based on two static GPS points collected with a CHC 90+ GPS postprocessed using the Auspos service. The position of all bathymetry points were recorded using a prism mounted on a 7.37-meter long pole which was placed on the river bottom by an operator in a small boat. The bathymetric survey recorded 116 points on the bottom of the river around the study area. The points were distributed around the survey area spaced out to ensure reasonable site coverage. The water level was recorded using a prism held at the water surface at 8 different locations. To calculate the bathymetry of the study area, a surface was interpolated using the Inverse Distance Weighting (IDW) surface interpolation function in ESRI ArcGIS software package. The IDW interpolation is a local neighborhood interpolation method, which estimates unknown values based on the distance and values of nearby known points. This interpolation method was chosen due to the limited number of points selected and the heterogeneity of the measured values (which is thought to be based on debris on the river bottom). Unlike kriging or other geostatistical surface interpolations, IDW should be less affected by the potential outliers in the data set . The surface was calculated using a variable search radius with 12 points and no maximum distance selected.
The surface produced can be seen in Figure 8(a). No bounding box was used for the surface interpolation, so the values produced are only accurate within the 60 × 15.5 m survey grid. The resulting surface provided elevations based on the Australian Height Datum for a continuous surface covering the study area. To calculate the water depth, the raster calculator tool was used to create a new raster surface, by subtracting the elevation of the river bottom from the average water surface elevation. This produced a continuous surface describing the depth of the water for the site.
The surface describing water depth for the site can be seen in Figure 8(b). Individual ERT lines were surveyed with a total station marking the start and end position of each line. The individual sensor locations were estimated by evenly spacing a number of points along the length of the measured line, using the “Points Along Line” tool in ESRI ArcGIS software package. An attribute was calculated for all sensors based on the water depth for use in later ERT processing. The elevation (above sea level) for the survey area was between −0.31 m and −4.13 m which corresponds to water depths between 3.47 m and 7.29 m.
4.3. Data Processing
In total, 112,984 raw data measurements were acquired during the field survey from survey lines L00 until L57. Data were exported after despiking outlier values (due to high geometric factor, small potential values, or values collected with insufficient current) and removing data with high % RMS error using a spreadsheet. The resistance values varied from 0.01 to 30 V/I, as depicted in Figure 9(a). Specifically, the first 12 lines (L00–L12) had much higher group of resistance values (from 0.5 to 30 V/I) and the lines that were over the wreck (L13–L57) mainly have values less than 0.5 V/I. Regarding the apparent resistivity values (Figure 9(b)), the first lines (L00–L12) have values ranging from 50 to 400 ohm·m, and the rest (L13–L57) have lower values from 1 to 20 ohm·m, reflecting the presence of the wreck in these lines.
After the preprocessing of the raw data, the a priori information was incorporated in each line. Before adding the water column depth in each survey line, the slope was smoothed in both directions : On direction the vertical distance should be less than the electrode spacing ( m) and on direction the vertical distance should be less than the line spacing ( m). The a priori information was completed after the resistivity value of the river water ( ohm·m) was added at the end of each survey line file.
High initial value of the damping factor (factor = 10) was selected (close to the surface) in order to cope with the high electrical contrast between the high conductive metallic target and the highly resistive background. The inversion for each line was completed using L2-norm (smooth) method for the data and L1-norm (robust) for the model since the target consists of sharp boundaries. Each 2D survey line was filtered individually (according to the %RMS error) and final 3D model images were produced after inverting all the 58 merged lines.
4.4. Inversion Results
In Figure 10, 12 different depth slices (D1–D12) represent the resistivity distribution at different depths below the river bottom where the first slice D1 is at ~0.10 m and the last one D12 at ~5.5 m, below the river bottom. The overall resistivity values (displayed on a logarithmic scale) varied from 0.05 to 1500 ohm·m.
At the NNE side of the grid from m to m, high resistivity values of 1000 ohm·m are present in all depth slices (D1–D12) to a maximum depth of 6 m. From m to m the background resistivity value varies from 30 to 50 ohm·m, from 0 to 1 m depth (D1–D4), and gradually increases to 1500 ohm·m into deeper layers till 6 m depth (D5–D12).
Until the first 0.5 m of depth (D1, D2), the exterior of the hull is difficult to discriminate. After this depth the boundary becomes more resolved to a depth of 1.5 m (D3–D6) with very low resistivity values that vary from 0.05 to 0.2 ohm·m. These boundaries are confined between m and m and between m and m presenting a structure of 44 m length and 10 m width, respectively. From the depth of 1.5 to 4 m (D7–D10) the conductive area becomes wider, presumably showing the ship’s deck. After the depth of 4 m (D11) the conductive area gradually disappears until it completely fades away after 6 m depth (D12). This is an indication either that the metallic part of the barge is no longer present or that the maximum depth of the barge is 6 m.
Apart from the external boundaries of the ship, internal sections are also resolved in the inversion depth images. At the depth of ~0.6 m (D3), a bulkhead (long metallic bar) can be seen running parallel to the centre line of the ship from bow to stern from m to m and at m. Furthermore, at m, 32 m, and 41 m at a depth of ~1.7 m (D7) internal sections, perpendicular to the main axis of the ship, are visible as well. A highly resistive anomaly is observed at a depth of 2.5 to 5 m (D9–D11) located at m to m and at m which is interpreted as a damaged part of the vessel, as validated by the multibeam data (Figure 2).
At the ESE part of the survey grid from m to m and m to m from a depth of 3 to 6 m (D10–D12) a low resistivity area parallel to the main ship axis can be seen with resistivity values close to 0.1 ohm·m. This has been interpreted as a detached part of the ship lying nearby although it may also be an artefact of the inversion process given its location near the boundary of the survey area. The RMS error is estimated as less than 30%, which may be attributed to the high contrast between the resistivity values of the anomaly related to Crowie and the surrounding substrate.
The inversion model values were exported into commercial 3D graphics software (Voxler). A three-dimensional representation of the barge can be seen in Figure 11, which also includes the bathymetry of the area. The resistivity values are mapped on a logarithmic scale. In Figure 11(a), the water layer (light blue) and the topography (light brown) are shown, using the volume rendered plot mode. In Figure 11(b) selected orthoimages (slices through the model) in all directions (, , and ) are depicted, indicating with dark blue colours the ships’ lower resistive areas. The overall volume with the target in it is shown in Figure 11(c) with highly resistive areas denoted by red colour and the target itself is shown in Figure 11(d) depicted with green colour using the isosurface mode of plotting. The isovalue used for ships’ metal structure representation is ohm·m. Apart from the outer boundaries of the ship, the internal bulkhead can be seen as well, across the long axis of the ship.
Overall, this study demonstrates that ERT is a suitable technique for investigating shipwrecks in aquatic environments. Despite this potential, some aspects of underwater ERT are more challenging than terrestrial surveys. For example, establishing the survey grid is a time consuming procedure as it is more difficult to set the four corners of the grid at the right position and keep them fixed throughout the survey period. Furthermore, placing the cable on the right position at each different line (if it is floating) is not an easy operation because it depends mainly on the current and wave conditions and any obstacles that might be protruding from the bottom. Any movement of the cable during the data acquisition can increase the noise levels, so extra stabilisation technique is required for some surveys. On the other hand, no metal probes are required since the cable itself is sufficient to propagate the current into the water and the river sediments. After the data acquisition, extra time for the data processing routine is needed since the a priori information including bathymetry and water resistivity should be included for each different survey line. The same arrays and equipment can be used for terrestrial and underwater ERT surveys; however aquatic surveys are more demanding on batteries due to the very conductive environment.
From an archaeological point of view, underwater ERT provides a new means of investigating material culture items that would otherwise be impossible due to turbid water (visual methods) or shallow water depths (acoustic methods).
The shipwreck Crowie was successfully mapped using ERT deployed underwater in the Murray River. Data from a model simulating the barge structure were obtained and analysed before undertaking the field survey, so that the optimum survey parameters (such as the appropriate array, probe spacing) could be selected according to the particular characteristics of this target and survey area. A priori information, including water column depth and water conductivity, were incorporated throughout the data processing for more accurate model representation.
The field survey involved 58 parallel lines being placed underwater on top of the wreck. Following processing and inversion, the resistivity distribution in the 3D inversion images clearly shows a conductive target shaped like the structure of a ship. The external boundaries of the ship, due to the low resistivity values (metallic parts), can be distinguished from the highly resistive background (sand and limestone sediments) and they are in good correlation with the historical record. The outer formation of the barge is accurately defined (in accordance with the modelling results) indicating the efficacy of the survey methodology. Apart from the successful imaging of the exterior of the ship, internal compartments of the barge are also visible where they present a suitable resistivity contrast.
Using the ERT method in this novel way to map a submerged barge provides important new information about this geophysical technique and about the wreck of Crowie. This result demonstrates that several important characteristics of a submerged ship can be defined, such as the burial depth of the wreck, the outer limits, overall dimensions, and crucially the shape of the ship. The excellent results obtained from Crowie suggest many future marine archaeological applications of the ERT method.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
This project was supported by an Endeavour Fellowship to Dr. Kleanthis Simyrdanis funded by the Australian Government. Lee Rippon, John Naumann, Celeste Jordan, Belinda Duke, Shannon Sullivan, and Anika Johnstone contributed to the field work for the project. Thanks are due to Lisa and Barry from the Morgan Waterfront Marina and ZZ Resistivity Imaging for their support of this research. Dr. Ian Moffat is the recipient of Australian Research Council Discovery Early Career Award funded by the Australian Government (DE160100703).
- C. Colombero, C. Comina, F. Gianotti, and L. Sambuelli, “Waterborne and on-land electrical surveys to suggest the geological evolution of a glacial lake in NW Italy,” Journal of Applied Geophysics, vol. 105, pp. 191–202, 2014.
- G. Tassis, P. Tsourlos, J. Rønning, and T. Dahlin, “Detection and Characterization of Fracture Zones in Bedrock: Possibilities and Limitations,” in Proceedings of the Near Surface Geoscience 2014 - First Applied Shallow Marine Geophysics Conference, Athens, Greece, September 2014.
- D. F. Rucker, G. E. Noonan, and W. J. Greenwood, “Electrical resistivity in support of geological mapping along the Panama Canal,” Engineering Geology, vol. 117, no. 1-2, pp. 121–133, 2011.
- L. Orlando, “Some considerations on electrical resistivity imaging for characterization of waterbed sediments,” Journal of Applied Geophysics, vol. 95, pp. 77–89, 2013.
- F. Baumgartner and N. B. Christensen, “Analysis and application of a non-conventional underwater geoelectrical method in Lake Geneva, Switzerland,” Geophysical Prospecting, vol. 46, no. 5, pp. 527–541, 1998.
- C.-W. Chiang, T.-N. Goto, C.-C. Chen, and S.-K. Hsu, “Efficiency of a marine towed electrical resistivity method,” Terrestrial, Atmospheric and Oceanic Sciences, vol. 22, no. 4, pp. 443–446, 2011.
- C.-W. Chiang, T.-N. Goto, H. Mikada, C.-C. Chen, and S.-K. Hsu, “Sensitivity of deep-towed marine electrical resistivity imaging using two-dimensional inversion: A case study on methane hydrate,” Terrestrial, Atmospheric and Oceanic Sciences, vol. 23, no. 6, pp. 725–732, 2012.
- G. Ranieri, F. Loddo, A. Godio et al., “Reconstruction of archaeological features in a mediterranean coastal environment using non-invasive techniques,” 2009.
- K. Simyrdanis, N. Papadopoulos, J.-H. Kim, P. Tsourlos, and I. Moffat, “Archaeological investigations in the shallow seawater environment with electrical resistivity tomography,” Near Surface Geophysics, vol. 13, no. 6, pp. 601–611, 2015.
- S. Passaro, F. Budillon, S. Ruggieri et al., “Integrated geophysical investigation applied to the definition of buried and outcropping targets of archaeological relevance in very shallow water,” Italian Journal of Quaternary Sciences, vol. 22, no. 1, pp. 33–38, 2009.
- S. Passaro, “Marine electrical resistivity tomography for shipwreck detection in very shallow water: A case study from Agropoli (Salerno, southern Italy),” Journal of Archaeological Science, vol. 37, no. 8, pp. 1989–1998, 2010.
- R. Parsons, Ships of the inland rivers: An Outline History And Details of All Known Paddle Ships, Barges And Other Vessels Trading on The Murray-Darling System, Gould Genealogy, Modbury, Australia, 3rd edition, 2005.
- A. Roberts, W. van Duivenvoorde, M. Morrison et al., “‘They call ’im Crowie’: an investigation of the Aboriginal significance attributed to a wrecked River Murray barge in South Australia,” International Journal of Nautical Archaeology, vol. 46, no. 1, pp. 132–148, 2017.
- M. H. Loke and J. W. Lane, “Inversion of data from electrical resistivity imaging surveys in water-covered areas,” Exploration Geophysics, vol. 35, no. 4, pp. 266–271, 2004.
- C. deGroot-Hedlin and S. Constable, “Occam's inversion to generate smooth, two-dimensional models from magnetotelluric data,” Geophysics, vol. 55, no. 12, pp. 1613–1624, 1990.
- J. B. Firman et al., Renmark 1:250000 Geological Map, Geological Survey of South Australia, 1971.
- P. A. Rogers, “Continental Sediments of the Murray Basin,” in The geology of South Australia, J. F. Drexel and W. V. Preiss, Eds., vol. 54, pp. 252–254, The Phanerozoic, South Australia, Australia, 1995.
- P. A. Rogers, J. M. Lindsay, N. F. Alley, S. R. Barnett, K. Lablack, and G. LKwitko, “Murray Basin,” in The geology of South Australia, J. F. Drexel and W. V. Preiss, Eds., vol. 2, pp. 157–163, The Phanerozoic, South Australia, Australia, 1995.
- V. Chaplot, F. Darboux, H. Bourennane, S. Leguédois, N. Silvera, and K. Phachomphon, “Accuracy of interpolation techniques for the derivation of digital elevation models in relation to landform types and data density,” Geomorphology, vol. 77, no. 1-2, pp. 126–141, 2006.
Copyright © 2018 Kleanthis Simyrdanis et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.