Remote Sensing and Environmental Monitoring Analysis of Pigment Migrations in Cave of Altamira’s Prehistoric Paintings

: The conservation of Cultural Heritage in cave environments, especially those hosting cave art, requires comprehensive conservation strategies to mitigate degradation risks derived from climatic influences and human activities. This study, focused on the Polychrome Hall of the Cave of Altamira, highlights the importance of integrating remote sensing methodologies to carry out effective conservation actions. By coupling a georeferenced Ground Penetrating Radar (GPR) with a 1.6 GHz central-frequency antenna along with photogrammetry, we conducted non-invasive and high-resolution 3D studies to map preferential moisture pathways from the surface of the ceiling to the first 50 cm internally of the limestone structure. In parallel, we monitored the dynamics of surface water on the Ceiling and its correlation with pigment and other substance migrations. By standardizing our methodology, we aim to increase knowledge about the dynamics of infiltration water, which will enhance our understanding of the deterioration processes affecting cave paintings re-lated to infiltration water. This will enable us to improve conservation strategies, suggesting possible indirect measures to reverse active deterioration processes. Integrating remote sensing techniques with geospatial analysis will aid in the validation and calibration of collected data, allowing for stronger interpretations of subsurface structures and conditions. All of this puts us in a position to contribute to the development of effective conservation methodologies, reduce alteration risks, and promote sustainable development practices, thus emphasizing the importance of remote sensing in safeguarding Cultural Heritage.


Introduction
The Cave of Altamira stands as a masterpiece of universal art, offering exceptional insight into a vanished cultural tradition, landscape, and technological heritage that illuminates a crucial period in human history.Recognized for its significance, the cave was designated as a UNESCO World Heritage Site in 1985.Situated within the Iberian Peninsula, specifically in the autonomous community of Cantabria, between the municipalities of Santillana del Mar and Reocín, The Cave of Altamira extends a length of 290 m (Figure 1).It holds a distinguished status as one of the most significant repositories of prehistoric art, with remarkable examples of Palaeolithic paintings dating back thousands of years.Particularly noteworthy are the depictions of signs and animals adorning the ceiling of the Polychrome Hall, offering invaluable insights into the beliefs and lifestyles of prehistoric hunter-gatherer societies.After the discovery of the Cave of Altamira in 1868, it underwent countless transformations, some aimed at preparing the cave for public visits and others to try to solve the serious stability problems affecting the cave ceilings.In parallel, to prevent water infiltration problems in the Polychrome Hall, a series of more or less fortunate interventions were carried out, among which was the injection of cement mortar over the central fracture that crosses the Ceiling from west to east, which caused significant changes in the water courses that had access through this large fracture, affecting the conservation of the paintings in areas near it [1].So, these and many other transformations affected and still affect the Cave of Altamira today.Its preservation, within an unstable conservation environment, is a major challenge for the cave's conservators and managers.
The main active deterioration factors in Altamira are directly related both to its internal microclimate and to its anthropic activities [2][3][4][5][6].The criteria used for the management and conservation of the Cave of Altamira have led to the application of a series of measures and protocols aimed at protecting the cave, its paintings, and its environment [1,7].Its preservation requires the implementation of global strategies that address the deterioration problems generated in its cave art due to the complex interaction between surrounding waters and the characteristics of the air in the cave (humidity, temperature, and CO2) [8].The variable action of these parameters produces alterations such as microcorrosion, dissolution, migration, and finally, detachment of the pigments constituting the cave art.In addition to this, the Cave of Altamira is a senile karst complex characterized by the presence of an important network of discontinuities (fissures, fractures, joints, detachments, stratification planes, and voids) from the upper part of the cave to its ceilings [9][10][11].These structural elements affect the cave's stability and condition significant changes in the infiltration water circulation, generating relevant alterations on the paintings such as flaking, washing/erosion, and carbonate concretions [12].
GPR is widely used in investigations of geological media, such as surveys of geological contact, groundwater, the detection of voids, fractures, and seepage, and even nondestructive testing of stone quality in quarries [13][14][15][16].Discontinuities in a rock medium (fractures, joints, voids, and bedding planes) represent electromagnetically active contacts/interfaces for GPR, where part of the radar pulse energy is reflected back to the surface, producing a radar record as a reflection profile (radargram) [17][18][19].GPR is used to efficiently characterize discontinuities in a karst system [20].Compared to other non-invasive techniques, GPR has been the most successful technique for high-resolution imaging of karst elements [21][22][23][24][25]. Additionally, GPR attribute analysis can provide further information from radar data [26][27][28][29][30].In this context, the integration of Ground Penetrating Radar (GPR) and photogrammetry has emerged as an indispensable tool to preserve rock art, offering non-invasive and high-resolution techniques to map and monitor cave environments.This integration presents a promising approach, allowing detailed 3D surveys and geospatial analysis crucial to understanding the factors influencing the degradation of Cultural Heritage, enabling the implementation of measures for its preservation [24,[31][32][33][34][35][36][37][38][39].In the GPR campaigns carried out in 2017 and 2018 in the Cave of Altamira, antennas with central frequencies of 100 MHz, 400 MHz, and 900 MHz were used, and through the integration of Global Navigation Satellite System (GNSS) data, a 3D Terrestrial Laser Scanner (3DTLS), and Unmanned Aerial Vehicle (UAV) photogrammetry, it was possible to map moisture zones and discontinuities in the overlying layer located internally right above the surface of the Polychrome Hall ceiling (Ceiling).This mapped set of moisture zones and discontinuities determined the main water infiltration pathways and related them to dripping points on the Polychrome Hall ceiling.In parallel, this study allowed for defining the main connection pathways for material and energy exchange between the outer surface and the inside of the Polychrome Hall.Additionally, the importance of the Dolomitic Layer and Polychrome Layer interface in the conservation status of the Ceiling was determined [35][36][37].
In 2023, a very high-resolution 1.6 GHz antenna was used to further study the internal structure of the Ceiling.This choice of a high-frequency antenna (>1 GHz) was intended to improve the resolution in detecting internal reflectors and changes/disruptions within the Polychrome Layer as close to the Dolomitic Layer as can be achieved, providing a reliable and non-invasive approach.
The proposed objective of the GPR technique is the high-resolution mapping of the water flow path in the Polychrome Layer from the Dolomitic Layer.The GPR data have allowed the estimation of discontinuities and moisture zones in the first 50 cm depth from the surface of the Ceiling.In this study, we refer to a control area ALT1 in the center of the Ceiling (Figure 2a).The ALT points depicted on the map represent sensitive areas where either pigment migration and carry-over processes or areas with bacterial growth have been detected.These locations undergo high-resolution microphotogrammetric monitoring to assess their evolution objectively.The control area under study is located on a large red clavicle placed near the legs of one of the large polychrome bison (Figure 2b).Claviforms are abstract signs dating from the Upper Palaeolithic period characterized by their elongated shape with a central bulge, although their precise meaning remains unclear; they are commonly in many caves, especially in the Cantabrian region, suggesting a cultural interaction between regional populations during that era.
The stratigraphy (Figure 2c) of the Polychrome Hall ceiling has been detailed in previous work [9][10][11].The orange layer, characterized by ochre and orange-colored limestone slabs with intercalated layers of loams and clays, is situated between the fissured layer and the upper level.The fissured layer comprises tableted limestone benches showing distinct morphologies resulting from dissolution processes.Additionally, the dolomitic layer, distinguished by a high degree of dolomitization and ochre coloration transitioning to brown upon alteration, is separated from the Polychrome layer by a fine clayey loam intercalation.These clarifications aim to enhance the reader's understanding of the stratigraphic composition.The geological layers directly influencing the processes affecting the ALT1 control area are the Polychrome Layer and the Dolomitic Layer [35][36][37].The Polychrome Layer, a 60-70 cm thick limestone layer, holds particular significance as Palaeolithic paintings were executed on its basal surface, using load-bearing structures and hydroplastic deformation morphology.The Dolomitic Layer, with a thickness ranging from 10-25 cm, exhibits located disintegration processes.The Dolomitic Layer overlies the Polychrome Layer and acts as a protective barrier, providing a degree of impermeability due to its thin intercalation of marls and clays in contact with the Polychrome Layer.Therefore, the Dolomitic Layer safeguards against water infiltration at the Ceiling surface.These two layers provide crucial insights concerning the site.
The results obtained in the Polychrome Layer using the 1.6 GHz antenna and its 3D data integration have allowed:

•
Obtaining a higher-resolution response in the first centimeters of depth from the surface of the Ceiling, where the mentioned red claviform sign is located.This has been complemented with earlier 2D records obtained with the central antenna of 900 MHz frequency.

•
Mapping the spatial configuration of moisture flows present in the study area directly involved in the Polychrome Layer.For this purpose, preferential moisture pathways have been mapped, as well as monitoring of the surface water dynamics, and the main dripping points associated with the migration and detachment processes of pigment.On the basis of these considerations, the main objective of this study is to provide valuable information that contributes to proposing reliable solutions to improve the conservation of the cave and its art.
Despite significant advances in remote sensing technologies, there are still problems when effectively applying these methodologies as basic tools for conservation in cave environments.Existing research on remote sensing has mainly focused on characterizing cave structures, with limited emphasis on correlating these findings with dynamics related to water infiltration and pigment migrations [40].The existence of divergent hypotheses about the effectiveness of remote sensing in adopting specific conservation measures highlights the need for comprehensive and standardized methodologies [1].Key publications in this field have highlighted the importance of remote sensing in obtaining useful basic information to implement measures aimed at safeguarding Cultural Heritage [41][42][43][44], emphasizing the need for interdisciplinary approaches that combine geophysical techniques with geospatial analysis.Controversies arise regarding the reliability and accuracy of remote sensing data, especially in complex cave environments where factors such as rock composition and moisture content vary significantly.As for moisture content, variations arise from different periods throughout the year, influenced by factors such as the access of infiltration water through fractures and fissures or changes in these water dynamics due to punctual variations in the order of moisture pockets or fillings, among others.
In this work, we show how remote sensing methodologies offer a valuable perspective, letting us understand the interaction between degradation agents (such as water/moisture) originating on the surface and, by infiltration, coming into contact with the paintings on the ceilings and walls of the Cave of Altamira.Taking advantage of the capabilities of remote sensing technologies, such as high-frequency GPR and photogrammetry, we aim to evaluate microclimatic environmental factors affecting the stability and integrity of the cave's ecosystem.
In our study of the Cave of Altamira, our aim was to understand the complex interplay between climatic variations and anthropogenic activities.While our primary focus centered on assessing damage related to moisture and pigment migration, we recognized the significant influence of human interventions on the environmental dynamics within the cave.Anthropogenic activities, such as applying cement mortar to seal fractures and constructing artificial walls, have changed the moisture circulation and microclimate, resulting in notable changes in the ecosystem of the cave.
Historical changes, including the installation of entrance doors and the construction of corridors with artificial lighting, have disrupted the natural state of the cave.These interventions over the past century have transformed the cave of Altamira from an oligotrophic environment to one significantly influenced by external factors related to anthropogenic action.Looking ahead, our research aims to dig deeper into these relationships to develop more comprehensive conservation strategies for safeguarding this invaluable cultural heritage site.
The main conclusions of this study emphasize the importance of integrating remote sensing techniques into conservation practices to preserve Cultural Heritage sites such as the Cave of Altamira.Thanks to high-resolution three-dimensional studies and geospatial analysis, we have improved our understanding of the complex interaction between environmental factors and the degradation of prehistoric paintings.Our findings highlight the essential role of remote sensing in helping with preventive conservation strategies involving decision-making and, in parallel, contributing to the long-term sustainability of these invaluable cultural assets.

Study Site
The study focused on the control zone designated as ALT1 located in the central part of the Polychrome Hall ceiling (Figure 3a).This zone covers an area of 1.2 m 2 and is arranged around a large red claviform sign located under the hind legs of one of the large polychrome bison.The detail of the sign corresponding to our control area ALT1 affected by various processes of migration and pigment washing associated with drip spots is depicted in the images.In the left image (Figure 3b), the sign seems scarcely affected by migration processes and paint washing.However, the right image (Figure 3c) shows the same sign eighteen years later, exhibiting numerous drip spots that have occasionally washed away the paint.These images collectively illustrate the progression of pigment loss and degradation.
In recent years, ALT1 has been affected by multiple infiltrating water drips that have caused dragging of pigment, leading to migration processes and, ultimately, the deterioration of the painting.This research poses significant challenges aimed at improving the conservation of this area and, by extension, other areas of the Ceiling affected by these same deterioration processes associated with water infiltration and consequent pigment dragging [45].

General Workflow Diagram
Integrating data into the model (Figure 4) begins with establishing a reference frame using the precision of the Global Navigation Satellite System (GNSS) [46,47].At first, in 2013, Topcon Hyper II receivers [48] were used to establish a geodetic reference system aligned with the European Terrestrial Reference System 1989 (ETRS89).Evaluating observed phase changes, aided by Topcon Tools 8.2 software [49], yielded a reference frame with an average accuracy of 1.7 cm for determining the coordinates of exterior cave vertices by using the method described in [50].
To ensure an accurate cartographic model of the cave, precise 3DTLS methods were employed [42].The process began from a closed traverse adjusted and compensated through the cave using a Topcon GPT-7503 (Topcon Corporation, Tokyo, Japan) total station [51].Sixteen traverse stations were established, covering a total length of 430 m, with an angular closure error of 0.0218 g and a linear closure error of X = −0.001m, Y = −0.005m, Z = 0.A total of 66 checkerboards, distributed throughout the cave, served as reference points for the 3DTLS, ensuring precise georeferencing of scans.The field campaign, conducted in December 2013, used a Faro Focus X-130 (Faro Technologies Inc., Lake Mary, USA) [52], requiring 300 scans.Calibrated spheres were employed as tie points, while traverse checkerboards served as references, achieving an adjustment accuracy of 2.7 mm for 95% of the points [53].In a previous photogrammetry campaign conducted in 2014, the Polychrome Hall Ceiling was documented using a Hasselblad H4 D-200 MS (Gothenburg, Sweden) camera, achieving a minimum resolution of 16 pixels/mm 2 .The photogrammetric process involved detecting homologous points extracted from the 3DTLS point clouds, strategically placing eighty evenly distributed points throughout the hall, with half serving as control points for model validation.These images were collectively adjusted to generate a high-density point cloud containing approximately 11 billion points, which was then filtered and generalized to produce another point cloud containing around 3.5 billion points.From this refined dataset, a high-resolution 3D digital model of the Ceiling was generated, comprising approximately 200 million polygons.This model facilitated the creation of a 6 Gigapixel orthoimage and a simplified version with 4 million polygons, enabling accurate integration.Utilizing the geometric models, hydrology calculations were conducted to comprehend the dynamics of water movement across the overlying layer and Polychrome Ceiling surface, analyzing surface water flow patterns and identifying potential entry points for water infiltration towards the Polychrome Hall [35][36][37].Finally, a dual-axis slider employing three 1.2 m sliders aimed to cover a 1 m 2 control area (ALT-1), with precise positioning achieved through scanning and 5 cm equidistance interval profiling.

Photogrammetry of the Polychrome Ceiling
The documentation of horizontal surfaces associated with gravitational issues poses challenges due to complexity and inaccessibility, prompting the use of advanced remote sensing techniques [53,54].In 2014, a photogrammetry campaign was conducted to document the Polychrome Hall ceiling [55], complementing laser scanning with a sub-millimeter resolution.The Hasselblad H4 D-200 MS camera, chosen for its R-G-B color capture and geometric accuracy, used HC 3.5/35 and HC 4/28 lenses, capturing frames at 8176 × 6132 pixels resolution.At least 16 pixels/mm² was maintained, with two 4800 Kelvin-degree LED screens providing 980 lux illumination, diffused with Rosco Cinegel 3000 Tough Rolux plastic.
Photogrammetric support involved homologous point detection from the 3DTLS point clouds.Eighty evenly distributed points, half serving as controls, were placed throughout the Hall.These images were collectively adjusted to generate a high-density point cloud of approximately 11 billion points.Subsequent filtering and generalization produced a refined point cloud with around 3.5 billion points.From this dataset, a highresolution 3D digital model was generated, boasting approximately 200 million polygons.This model helped with the creation of a 6 Gigapixel orthoimage and a simplified 4 million polygon version, ensuring seamless integration for the study of rock art and the associated processes of alteration.Once the geometric models were available, hydrological calculations were carried out, as understanding the dynamics of water flow through the overlying layer and the ceiling surface of the Polychrome Hall plays a fundamental role in understanding the active alteration processes associated with the ceiling surface paintings.
The data collected in 2014 maintain their relevance and reliability for our analysis, even considering potential changes in the cave's condition.The photogrammetric campaign conducted in 2014 provided valuable insights into the cave's micro-morphological particularities, particularly regarding pigment conditions and surface alterations.While minor changes may have occurred since then related to microcorrosion processes of the supporting rock, these are at the micron level and are very subtle changes.These are unlikely to significantly affect the integrity of the data or the accuracy of our analysis.Additionally, ongoing monitoring efforts, including annual high-resolution photogrammetry of key zones like ALT1, ensure that any changes are captured and accounted for in our assessments.So, we affirm that the data from 2014 remain a robust foundation for our hydrological calculations and the understanding of alteration processes within the cave.
The aim of integrating the geometric models with hydrological analyses was to gain insights into how water interacts with the internal structure of the Ceiling over time, revealing patterns of erosion, disintegration, concretion, and moisture infiltration, potentially revealing patterns.

GPR of ALT1 Control Area
The GPR technique consists of the propagation of short electromagnetic pulses (t ≤ 1 ns) in a medium.These pulses reflect when there are changes in the electromagnetic properties of the medium (magnetic permittivity, electrical conductivity, and dielectric permittivity).So, the characteristics of the GPR signal are strongly influenced by the electromagnetic properties of the medium and their contrasts with each other.The contributions of electrical conductivity and magnetic permeability to the recorded signal cannot be distinguished in typical GPR studies using frequencies between 8 MHz and 2.5 GHz [56,57].In low-loss media, the dielectric permittivity (εr) can be estimated by the equation [58,59]: where h is the depth of the layer, t is the two-way travel time (elapsed time) in the GPR profile, and c is the velocity of light (0.299 m/ns).The mean velocity (v) is derived from: The amplitude of a registered signal is a function of several factors: the amplitude of the transmitted signal, the geometrical dispersion (and thus the travel distance), the reflection coefficients between the interfaces (due to targets causing reflection and refraction of the signal), the signal scattering (due to small targets causing scattering), and the absorption of the medium.
The absorption of the medium is defined by the attenuation coefficient α.The attenuation coefficient is dependent on several factors (the angular frequency (ω), the electrical conductivity (σ), the dielectric permittivity (ε), and the magnetic permeability (μ)), as given in Equation (3) [60], with the average conductivity and the dielectric permittivity of the medium being key influences: Since most of the GPR studies are carried out in non-magnetic media (low-loss media), in this, case Equation (3) simplifies to Equation (4) [61]: Equations ( 3) and ( 4) show the dependence on the dielectric parameters of the media that depend on their physical properties, with the water content being one of the important factors.This attenuation coefficient α represents the amplitude behavior during the propagation of the electromagnetic wave through the medium as a result of energy absorption.
The reflection amplitudes at the two-material interface can be estimated using the reflection coefficient (R).GPR signal reflection amplitudes and intensities are proportional to the contrasts of the dielectric permittivities at the reflector boundaries and depend on the reflection coefficients.In non-magnetic media (low-loss media), expressions for the approximate electromagnetic impedance (Z) and reflection coefficient (R) for two media are given by [59,60,62]: where εr1 is the dielectric permittivity of the upper layer and εr2 is the dielectric permittivity of the lower layer; Z1 is the electromagnetic impedance of the upper layer and Z2 is the electromagnetic impedance of the lower layer.This reflection coefficient equation implies that the greater the contrast between the two media, the higher the percentage of incident energy that will be reflected at the boundary/discontinuity of the reflector.If the moisture/water content in a medium is significant, this moisture/water content becomes an influential factor for the values of dielectric permittivities of materials/media, as the dielectric permittivity of water ( 81) is higher than that of any other geological material   [59,63].
Thus, a higher water content in the interior of the medium results in two effects on the amplitude of the electromagnetic wave:

•
The effect of attenuation on the amplitude is a derivative of the attenuation coefficient.The amplitude of the wave propagating through the materials/media decreases drastically as a function of the degree of moisture.Consequently, the wave is attenuated and the wave incident on the targets within the medium generates a low-amplitude anomaly in the radar records, and in some cases, the anomalies may not even be visible [22,[64][65][66][67][68].

•
An additional effect on the amplitude is a function of the reflection coefficient.In the case of surface moisture in areas under dry material/medium, the contrast in dielectric permittivity between the wet and dry media is greater the higher the water content is, increasing the reflection coefficient.This means that the greater the contrast between the electromagnetic parameters of two materials/media that are in contact, the greater the percentage of incident energy that will be reflected at the boundary/discontinuity of the reflector.The strength (amplitude) of the reflected fields is proportional to the magnitude of the dielectric constant.Radar reflections of higher amplitude will occur mainly at interfaces within the same geological layer and at material/media interfaces that have significantly different electromagnetic properties.Thus, the anomaly due to the existence of moist areas has a higher amplitude in the reflection profile, as has been reported by different authors [20,37,60,[68][69][70][71][72][73][74].
As shown in Figure 2c, the rock medium studied with the 1.6 GHz antenna has been limited to the limestone Polychrome Layer due to the depth of penetration achieved.The surface tension water and active drip points observed on the ceiling are evidence of internal moisture zones in this limestone layer.In the dry zones of the limestone, the relative dielectric permittivity is lower than in the moist zones [59,69].The existence of wateraffected (moister) zones within the dry limestone medium produces interfaces with different relative dielectric permittivities.In the dry zones, the relative dielectric permittivity is lower than in the water-affected (moister) zones.As a consequence, when the wave reaches these interfaces, part of the energy is reflected.According to expression (5), highamplitude radar reflections are also generated at the interfaces of the water-affected (moister) and dry zones of the Polychrome Layer.It is relevant to notice that the existence of moisture in the interior of the limestone medium (Polychrome Layer) causes the limestone-water interface to have a higher dielectric contrast than the limestone-air interface and, consequently, a stronger reflection [69,70,75].
GPR signal penetration and detection resolution depend on the antenna frequency used.Higher-frequency antennas offer better detection resolution but also increase signal attenuation, which results in lower penetration into rocks [59].The penetration depth at which the GPR technique can record discontinuities within a rock mass is strongly influenced by its electromagnetic parameters (dielectric permittivity, electrical conductivity, and magnetic permeability).Variations in the electromagnetic parameter values can indicate alterations in the rock and conductive minerals, porosity, joints, fractures, water content, and fillings within the discontinuities.In geological media, the attenuation factor increases as the conductivity of the medium increases.In addition, the propagation of electromagnetic waves in a rock mass is also affected by attenuation due to the degree of water saturation, water mineralization, distance, porosity variation, and the number of discontinuities [60,[76][77][78][79].The high-frequency antenna choice (>1 GHz) provides a higher resolution in detecting the dimensions and geometries of reflectors/elements in a medium, which may be associated with internal changes/discontinuities [58,80,81].However, the penetration depth with these high-frequency antennas is lower than the 900 and 400 MHz central-frequency antennas used in earlier GPR campaigns.Due to the range in penetration depth of high-frequency antennas (>1 GHz), this survey was focused on the Polychrome Layer.
The comparison between antennas in GPR studies highlights a trade-off between resolution and penetration depth.The high-frequency antenna offers superior resolution, capable of detecting features as small as 2 cm, but its penetration depth is limited to 0.4 m.However, the low-frequency antenna provides a lower resolution but can penetrate up to 4.5 m, making it more suitable for deeper subsurface investigations [56,58,[82][83][84][85][86].The choice of antenna is closely related to the resolution of a specific conservation problem as accurately as possible: the washing of paint, which is directly linked to water accumulated in the first centimeters from the basal surface of the Ceiling.This justifies the use of a higher-frequency antenna (higher resolution and lower penetration), as shown in Table 1.Additionally, addressing deeper accumulations indirectly requires an antenna with greater penetration.Ultimately, the objective is to respond to a conservation need, balancing the requirement for resolution against the need for depth penetration to ensure optimal results in subsurface imaging and detection applications [56,58,[82][83][84][85][86].The digital processing of time-domain GPR data generally relies on the amplitude values of the signal for interpretation.Furthermore, the transformation of the GPR data into its instantaneous magnitude, phase, and frequency information requires additional processing of the GPR data.The Hilbert transform has been applied to complex rock masses to enhance the reflection and absorption contrasts of electromagnetic wave propagation.Thus, GPR attributes can reveal or highlight features or patterns that are not clear in standard amplitude data [86][87][88].
In this study, using Hilbert transform attributes in the interpretation of GPR data can help us to increase our knowledge of the internal structure of the Polychrome Layer by characterizing the electromagnetic response to the moisture zones within it.

GPR Data Acquisition
The high-resolution GPR survey was performed using an SIR 4000 system manufactured by Geophysical Survey Systems, Inc. (GSSI, Nashua, NH, USA).The GPR system was equipped with a 1.6 GHz central frequency antenna.The choice of a high-frequency antenna (>1 GHz) was intended to improve the resolution in detecting internal reflectors and changes/disruptions at a 50 cm depth in the Polychrome Hall ceiling.This antenna configuration provided depth penetration ranging from 0 cm to approximately 50-60 cm, depending on the medium and its electromagnetic properties, while maintaining a highresolution capability [58,80].
Due to the presence of Palaeolithic paintings on the Ceiling, the GPR data were collected with the air-coupled antenna on the Polychrome Hall Ceiling to avoid any contact between the antenna and the Ceiling that could cause damage to the cave paintings.It had to be assumed that the antenna height would affect the data recording process.To reduce this effect and consider the minimum allowable distance to the rock paintings on the Ceiling, data acquisition was performed at a height depending on the 1.6 GHz central-frequency antenna of 10 cm (Figure 5).This height ensured the maximum possible amount of energy transmitted to the rock mass and reduced the effects of coupling to the air to obtain accurate and high-quality GPR records [89,90].This antenna configuration provided a penetration depth of up to approximately 50 cm into the Polychrome Layer, while maintaining a high-resolution capability [58,80,81].A 1 m-long slider was chosen due to operative limitations inside the cave, on which the 1.6 GHz antenna was placed, focusing on the Ceiling (Figure 5).An odometer wheel (encoder) was used to calculate trace-interval distance.Further, the reflection profiles were georeferenced through 3D scanning of the sliders where the antenna was fit.
The survey aimed to spatially map the internal distribution of moisture in ALT1 and its immediate surroundings.To achieve this, a 5 cm × 5 cm grid pattern was designed specifically for ALT1, including 23 longitudinal and 23 transverse profiles.These profiles were strategically laid out and numbered for the study's purposes, as illustrated in Figure 6.Additionally, a square grid measuring 1 × 1 m was established, using equidistant orthogonal 2D lines spaced every 5 cm.This grid was centered on the ALT1 control point, as depicted in Figure 5, with the 2D lines oriented transversally and parallel to the central fracture of the Polychrome Hall ceiling.This approach helped with comprehensive data collection and analysis of moisture distribution within the designated area.

GPR Data Processing
The GPR wave mean velocity was determined by the hyperbolic fitting method on a set of hyperbolas recorded in different reflection profiles, resulting in a mean velocity value of 0.111 m/ns.The dielectric permittivity (ɛ) was calculated to be 7.2, according to Equation (2).To calculate the processed depths in the investigated polychromatic layer, this value was applied.This mean velocity was also used in the Kirchhoff migration for the 2D processing of the GPR data.
The raw field GPR dataset needs to have a post-acquisition processing procedure applied to it before generating amplitude maps and rendering 3D images with all reflection profiles on a grid.There are several techniques for post-acquisition processing of radar signals.In this study, a basic data-processing procedure and a Hilbert transform attribute (instantaneous magnitude) were applied to the raw dataset, using RADAN 7.6 software (Geophysical Survey Systems, Inc., GSSI, Nashua, NH, USA) [91].The first step was 1D processing consisting of zero-time correction (time-zero adjustment).After this, the 2D processing was applied for the reflection profiles according to the following steps: (i) background removal filter; (ii) bandpass filter; (iii) linear amplitude gains; (v) Kirchhoff migration filter using the calculated average velocity (0.111 m/ns); (vi) a time-to-depth conversion performed based on a calculated velocity mean velocity of 0.111 m/ns for the limestone of the Polychrome Layer; (vii) Hilbert transform attribute: instantaneous magnitude.Next, each processed reflector profile was aligned within the grid to generate horizontal amplitude 3D maps and isosurface (isoattribute-surface) images that allowed us to pinpoint the presence of significant moisture/water zones.Several amplitude slice maps were created at different levels for mapping the subsurface up to depths of 36 cm.Amplitude depth-slice maps from the 3D model were used to identify moisture/water zones at a constant depth.Once the depth-slice maps had been generated, the highest amplitudes in all profiles within the grid were rendered in an isosurface image.This isosurface image was colored to show certain amplitude values (the highest ones) while others were made transparent in order to best represent only moisture/water zones and simplify locating moisture/water zones and their interpretation.This transparent visualization of the 3D GPR dataset was performed to reveal the main moisture/water zones in this studied volume of the Polychrome Layer.

Polychrome Ceiling Hydrology
Modeling water movement within the cave's interior poses a significant challenge, particularly in the Polychrome Hall ceiling where water primarily moves under the influence of surface tension.This phenomenon occurs as water infiltrates through the porous and fractured limestone rock, guided by capillary forces and the cohesive properties of water molecules [92][93][94].GIS technologies play an important role in accurately modeling this process.Identifying the drip points involved that induce processes of erosion, flaking, disintegration, and washing by dissolution of some of the mineral parts, including the constituent pigments of the Altamira art (Figure 7), holds profound implications for cave conservation.Understanding the flow paths and points of ingress is essential for predicting the potential damage and degradation of rock paintings.The workflow involves using a 3D high-resolution photogrammetric model to conduct a hydrological analysis [95] in QGIS [96], aimed at elucidating the sources and destinations of water within the area.Various analytical functions, including the Strahler method [97], are applied to accurately represent water flow dynamics over the Ceiling.The Strahler method assigns numerical values to segments of a raster, reflecting branches in a linear network.It classifies streams based on their tributary counts, with an increase in stream order at intersections.The "Basin" function delineates drainage basins by identifying ridge lines that separate them, which is important for understanding spatial water distribution [98].Using input flow direction raster data, the analysis identifies connected cells within the same basin.The "Flow Distance" analysis calculates horizontal or vertical downslope distances along flow paths to converging stream cells, helping with the computation of various flow distance metrics.
The workflow includes the addition of drip-point calculation, which concerns specific locations within the cave Ceiling where accumulated water reaches a critical saturation point and drips down and constitutes an important feature to map and understand.They can significantly affect rock art deterioration, cave formations, water quality, and ecological dynamics (Figure 8).By simulating water movement and saturation levels within the cave Ceiling, GIS aids in predicting the precise areas of dripping.Identifying these drip points is important for geological and hydrological research and conservation efforts.It allows for tracking potential erosion, calcification growth, and the transport of water chemistry and nutrients to the cave floor, influencing unique cave ecosystems.All this information allows the managers of the Cave of Altamira to make decisions to anticipate damage to reduce the deterioration of the cave and the paintings it houses.

Results
The results obtained with the 1.6 GHz antenna shown in Section 3.1 indicate the distribution of moisture in the internal part and for the first 36 cm from the surface of the Polychrome Hall ceiling in the ALT1 control zone.Hydrographic calculations based on a three-dimensional model reveal the dynamics of water flow on the Ceiling, aiding in understanding the active fissures involved in pigment drag processes, as well as related drip points (Section 3.2.).Six relevant moisture zones with varying elevations and characteristics are observed.These zones are linked to surface microchannels and the large central fracture mapped in previous campaigns (2017 and 2018) as shown in Section 3.3...Despite encountering some issues with data acquisition using the antenna, such as noise interferences and depth limitations, the results allow for the interpretation of water and moisture access dynamics to the Ceiling surface.This information lays the groundwork for applying conservation strategies to minimize deteriorations related to infiltration water on the Polychrome Hall ceiling of the Cave of Altamira.

Moisture Mapping in ALT1
In the ALT1 control zone, after processing GPR data, six high-moisture zones were detected, each with different characteristics (Figure 9b,c):

•
Moisture Zone 1 is strongly associated with the geometry of the central fracture mapped in the 2017/2018 campaigns.This moisture accumulation is parallel to the central fracture and reaches a depth of 36 cm (Figure 9a).Moisture Zone 1 results from the moisture concentration in areas adjacent to the central fracture due to the dam effect generated by the injection of cement mortar applied over the entire development of this large fracture crossing the Ceiling from west to east (Figure 2a).The antenna information shows the existence of discontinuities in the mortar, indicating partial degradation, disintegration, or deterioration at certain points.Additionally, the antenna indicates that the depth of this injected filling from the early last century presents a variable thickness.

•
Moisture Zone 2 is closely associated with drip points ALT1_1, ALT1_2, and ALT1_4 as shown in Figure 7b, between 3 cm and 5 cm deep.Its highest moisture concentration is at a depth of 3 cm.These described moisture zones correlate with both the large fractures (mapped in the 2017 and 2018 campaigns) that develop vertically from the first meters of the overlying layer towards the Ceiling surface and the small fissures or fractures observed on the

Hydrology of the Polychrome Ceiling
Hydrographic calculations have been carried out using a detailed 3D model with a ground sampling distance (GSD) of less than 200 microns to reveal critical geological characteristics.The main objective is to identify the contribution basins, streams, and drip points, providing information on the dynamics of water flow in the studied environment.The obtained basins constitute an important metric that allows an understanding of the extent and activity of the fissures.In parallel, channels are accurately mapped to discern the trajectory of water flow, including possible pathways for pigment drag and other potential deposits such as clays or limestone particles, in some cases related to microcorrosion processes induced by another significant alteration agent present in caves such as CO2 ( [99][100][101][102].This high-resolution model helps to thoroughly examine the hydrographic network, contributing to a global understanding of the existing hydrogeological dynamics on the Polychrome Hall ceiling (Figure 11).To ensure the accuracy and validity of hydrographic calculations, a rigorous validation process has been executed.The results are compared with real field observations or through orthoimages, thus verifying the accuracy of the calculated drip points.This verification step is essential to affirm the fidelity of the 3D model predictions, aligning virtual representations with real geological features.In this verification process, information obtained through the so-called "witnesses" located on the floor of the Polychrome Hall is also important (Figure 12).These witnesses, covering an area of 1 × 1 m, are made of rigid material formed from two thin aluminum sheets joined to a solid polyethylene core and lined with a 100% synthetic fabric made of high-density polyethylene fibers, which are lightweight, durable, breathable, and resistant to water, abrasion, and bacterial penetration.The witnesses are georeferenced, allowing the documentation of any dripping from a specific area of the Ceiling surface.Their information, combined with the aforementioned hydrological study, determines the origin of clay, limestone, and pigment drag deposits.This valuable information is transferred to a series of tables that inform us of the dynamics and periodicity of these events [103].The importance of these hydrographic calculations goes beyond mere visualization as they provide valuable information about water distribution and pathways in the inner part of the Ceiling (see Figure 11).Drip points, indicating areas where water ultimately precipitates, are important indicators of the hydrological behavior of the environment.This detailed analysis not only improves the understanding of water access through active fissures and water channels but also sheds light on specific areas where water interacts with the cave floor, offering broader insight into the geological processes at play in the studied area.

The Central Fracture
The presence of a higher degree of superficial moisture, which is concentrated in the central part of the Polychrome Hall ceiling in the form of small laminar water flows on both sides of the large fracture sealed with cement, critically influences some of the pigment drag and washing processes involved in our ALT1 control zone (Figure 13).Geophysical prospecting studies conducted in the Cave of Altamira in recent years have revealed that the large central fracture has vertical development that, from a depth of 1.4 m relative to the exterior surface, penetrates, reaching the actual surface of the Polychrome Hall ceiling.On the other hand, the aforementioned injection with cement mortar acts partially as a 'plug', diverting water flows towards adjacent areas, precisely where our ALT1 control zone is located (Figure 11b).

Limitations and Challenges
The problems regarding data acquisition are derived, among other issues, from the difficulty in locating the areas of interest, as they are situated on the Ceiling, which obviously does not allow any contact of our antennas with its surface.On the other hand, the Ceiling surface is bathed in a thin but constant layer of water, which also conditions data collection.Additionally, there are significant conservation constraints related to both the length of stay and the use of materials inside the Polychrome Hall [1].
Noise interferences and signal attenuation have limited the depth penetration of data acquisition beyond the first 50 cm, making subsurface characterization difficult.The combined use of different antennas and the application of a novel methodology by using antennas from the exterior surface in combination with those applied from the interior of the Polychrome Hall has allowed us, after correlating all the data, to obtain more extensive and more precise information on internal fractures and the moisture involved.
The results reveal complex relationships between the distribution of moisture and its associated fractures and fissures and highlight how certain environmental factors decisively influence the conservation of this type of Cultural Heritage site.

Discussion
The results of our study shed light on the complex dynamics of moisture distribution in the Polychrome Hall of the Cave of Altamira.Thanks to the integration of high-resolution GPR data with previous campaigns and geospatial analysis, we have been able to better understand the underlying processes related to the dynamics of infiltrating water affecting the conservation of this Cultural Heritage.
The geophysical prospecting campaign with the 1.6 GHz antenna in the ALT1 control zone and up to the first 50 cm depth allowed us to conclude that the highest amount of internal moisture is found between 10 and 20 cm from the surface, with the least moisture load found at greater depths of above 30 cm.Additionally, the presence of this moisture is mainly related to two lateral contributions, one on the southern side, clearly associated with the large central fracture, and another on the northern side coming from a large fracture of similar development to the central one that even begins at a higher elevation and, discontinuously, reaches almost the same depth as the central one (Figures 10 and  11b).
The greatest water percolation is located towards the northern side, taking as a reference the large central fracture that crosses the Ceiling from west to east.This direction of moisture flow could be explained by: a) The dip of the strata of the Altamira karst itself.
b) The greater percolation of moisture in the western area of the studied quadrant is due to the sealing with cement mortar applied over this large fracture.

Interpretation of Results in Light of Previous Studies
Our observations are based on studies carried out in the Cave of Altamira, particularly those using GPR technology with variable frequencies.The use of a 1.6 GHz antenna allowed for improved resolution in detecting wet areas and structural characteristics in the upper layers of the cave Ceiling.A comparison with previous campaigns using lower frequencies helped to better understand subsurface conditions and temporal changes.

Implications for Conservation Strategies
The identification of different moisture zones, areas, and channels and the origin of deposits detached from the Ceiling, especially those of pigment and limestone and dolomite particles, provides important information as they speak of the activation of a series of deterioration processes associated with microclimatic variations.The correlation between moisture distribution and active drip points underscores the need to adopt proactive conservation measures to mitigate the potential damage caused to rock art by the access and circulation of this infiltrating water on the Ceiling surface.

Broad Context and Future Research Directions
Our results contribute to the broader discourse on the conservation of Cultural Heritage in cave environments, highlighting the importance of integrating remote sensing methodologies for effective conservation.Future research directions may include:

•
Further studying moisture dynamics and its correlation with seasonal variations to refine conservation strategies.

•
Continuous monitoring of structural stability to assess long-term changes and the effectiveness of possible conservation interventions.
Considering these lines of research would advance our understanding of the complex interactions that shape the preservation of such valuable Cultural Heritage, which in turn would help implement new and more effective preventive conservation strategies.Establishing sustainable practices based on knowledge of these active deterioration processes will enable the activation of new control and monitoring protocols.Ultimately, this will put us in a better position to understand the processes involved in the deterioration of rock art, allowing us to anticipate damage and foresee deterioration in top-tier Heritage Sites such as the Cave of Altamira.

Conclusions
Our study represents a significant advancement in understanding the deterioration factors involved in the conservation of art in the Cave of Altamira, especially in the Polychrome Hall.Leveraging advanced remote sensing techniques, particularly high-frequency Ground-Penetrating Radar (GPR), along with geospatial analysis, we have achieved a comprehensive characterization of moisture distribution, and the constraints affecting the structural stability of the Ceiling have been better understood.
The data obtained provide relevant insights into the factors influencing the degradation and conservation of the cave's prehistoric paintings, elucidating the complex interaction between climatic influences, anthropogenic activities, and geological characteristics.
By delineating moisture zones and identifying detachment areas, we have established a foundation for implementing new preventive conservation strategies aimed at mitigating risks and proposing measures to safeguard this invaluable Cultural Heritage.
This study underscores the importance of interdisciplinary collaboration among professionals from different fields to address the multiple challenges posed by the conservation of this type of Heritage.By integrating research for conservation purposes, we can develop holistic approaches focused on the long-term sustainability of these heritage sites.
The research has allowed us to determine the trajectory of water movement in the first 50 cm from the ceiling surface, enabling the identification of the main internal structural discontinuities involved in water access to areas where cave art is distributed.This study was repeated at different times of the year, allowing observations of variations in water behavior and movement within the ALT1 control zone located in the center of the Polychrome Hall ceiling.
Looking ahead, our work lays the groundwork for exciting avenues of research, such as continuous monitoring of the environmental dynamics involved in rock support and pigment detachment processes, the refinement of conservation strategies, and the exploration of innovative technologies for the non-invasive characterization of cave environments [104].Additionally, our findings underscore the need for ongoing commitment to local communities and stakeholders to ensure effective management and protection of Cultural Heritage resources.
The conservation of the Cave of Altamira faces significant challenges, especially regarding deterioration processes related to its geological stability and the infiltration of water.To mitigate these threats posed by active deterioration processes, it is essential to develop effective conservation strategies that improve our understanding of the current conservation status of the cave.
The integration of geomatics remote sensing technologies, such as 3DTLS, UAVs, and GPR, has helped to rapidly, accurately, reliably, and non-invasively document complex elements such as fracture systems involved in the main access routes for water infiltration into the cave.We have correlated these fractures with several hydrographic basins dividing the surface of the Polychrome Hall ceiling, water streams, and drip points, inducing erosion, scaling, disaggregation, and dissolution processes affecting mineral components, including the pigments constituting Altamira art.With this knowledge, Cave of Altamira managers will be able to make proactive conservation decisions to prevent damage and reduce the risk of cave and painting deterioration [105].

Data Availability Statement:
The research data supporting this publication are not publicly available.The data were collected by GIM Geomatics as part of research and conservation studies of the Cave.These data are kept in the Museo Nacional y Centro de Investigación de Altamira.

Figure 1 .
Figure 1.Location map of the Cave of Altamira: projection of the Cave of Altamira on the topographic surface, showing the Polychrome Hall and the orthoimage of its Ceiling.

Figure 2 .
Figure 2. (a) Location of control area ALT1 in Polychrome Hall ceiling.The dashed white line indicates the position of the central fracture.(b) Picture of the general view of ALT1 study area.(c) Schematic stratigraphy of the Polychrome Hall ceiling according to [9-11], showing the depth range reached with the 1.6 GHz antenna in Polychrome Layer.

Figure 3 .
Figure 3. (a) Control zone ALT1 and comparative pictures between (b) 2004 and (c) 2022 where migration and pigment dragging and pigment loss processes associated with drip points occur.

Figure 4 .
Figure 4. General workflow scheme followed in this study.

Figure 5 .
Figure 5. Photograph showing the device designed (a) in which the 1.6 GHz antenna was installed (b) with the odometer (c) maintaining a minimum distance in relation to the study area of the Polychrome Hall ceiling of 10 cm.

Figure 6 .
Figure 6.Detail of the control area ALT1, with layout of the longitudinal and transverse profiles made in ALT1 and image of the field configuration.

Figure 7 .
Figure 7. (a) Location of the ALT-1 active drip zone near the hind legs of the great bison figure.Drip points are indicated with blue polygons and streams with blue lines.The thicker the blue lines (streams), the higher the hierarchy in the ALT-1 active drip zone.(b) Active drip basins highlighted in green involved in the recorded pigment migrations in recent years.Black circles indicate areas of paint drop and loss and white circles indicate areas where pigment has been transported by the water film with no pigment fall recorded in recent years.

Figure 8 .
Figure 8. Image showing the limestone rock stripped of paint due to the process of dragging and dripping of water.The white circle highlights a new drip that has already caused the loss of pigment adhesion to the supporting rock.

•Figure 9 .
Figure 9. (a) Reflection profile L5 (greyscale amplitudes) showing high-amplitude and strong reflections due to water, interpreted as the presence of moisture within the Polychrome Layer (highlighted in red rectangles: zones 1 and 6).(b) Areas of water/moisture: Superimposition of the sections of the first 36 cm.In blue the areas with higher presence and in white the areas with lower presence.(c) The 6 areas of high moisture mapped in the 3D isosurface model in ALT1 (highlighted with red numbers).

Figure 10 .
Figure 10.(a) General view of the main surface fissures and fractures in the ALT-1 control area.(b) View of some of the small fissures or fractures that act as small contributions of moisture affecting, as highlighted in the box, the drip point ALT1_4.(c) Detail view of ALT1_4 control area.

Figure 11 .
Figure 11.(a) Representation of hydrographic basins and active drip points over a wide surface of the Polychrome Hall ceiling.(b) 3D view showing the moisture distribution within the large central fracture and the control zone ALT1.It showcases the surface-level features like basins, streams, and drip points on the Ceiling area, with varying shades showing the degree of moisture accumulation.Lighter colors mean lower moisture accumulation, while blue indicates the highest accumulation of moisture.The green box highlights the internal moisture close to the Ceiling surface that is directly involved in the processes of paint carry-over, migration, and loss in the ALT1 control area.

Figure 12 .
Figure 12.Location of the different witnesses located on the floor of the Polychrome Hall.Witness covering a series of drips located on the hump and head of one of the large polychrome bison.

Figure 13 .
Figure 13.(a) Large central fracture sealed with cement mortar (indicated by the white arrow) and location of the nearby claviform signal (white box).(b) Drip point ALT1_1 and its associated main water inputs (blue arrows indicate the location and orientation of the marked water flows).

Table 1 .
Comparison of parameters for the 400 MHz, 900 MHz, and 1.6 GHz antennas used on Polychrome Hall ceiling.