A Comprehensive Gridded Dataset Associated to the Climate Change Effect on the Water Resources in the Grand Est Region, France

Water resources and environment quality are nowadays under high pressure because of climate change, land use practices, as well as human actions. A comprehensive gridded dataset becomes a necessary instrument to assess the risk level at regional scale, and also for territorial planning, the defining strategies to address future natural and anthropological challenges. In order to obtain a complete database with the most important parameters at spatial scale, this study is constructed as a preparation of layers used for various environmental risks, but mostly with the climate change effect on the water resources from the Grand Est region, France. In addition, geological formations, terrain data, and land cover were harmonized as grid format for the study area. Thus, the temperature and precipitation parameters, related to the 1961–1990 (1990s), 2011–2040 (2020s), and 2041–2070 (2050s), become useful data for evapotranspiration, water availability, and effective precipitation calculations. The geology layer indicates the composition and types of aquifers and it contributes to the potential infiltration map (PIM). The morphology of the terrain contributes to the slope angle and PIM. Through the typology of land cover, the pollution load index (PLI) was estimated. The findings indicate intense aridization and the depletion of the effective precipitation (below 650 mm) during the present and future periods. With respect to these concerns, the surface waters and groundwater resources from the Grand Est region are experiencing the negative effects of climate change on runoff and aquifers recharge respectively. In addition, the high PLI in the industrial and agricultural areas contribute to the possibility of the increasing water resources vulnerability. The affected areas extend mainly in the western, north-central, and north-eastern parts of the region, mainly in the Rhine, Aube, and Marne Valleys. Considered as a precious resource in the region, the water management should follow best practices for vulnerability and risk assessment, and further to delineate the protection areas. As a comprehensive gridded dataset, the calculations and original maps presented in this paper represent a complex product with main environmental parameters processed at spatial scale of 1 km2 in ArcGIS. This product has the purpose to integrate the geospatial data for the Grand Est region of France.


Introduction
Climate change and intensive land use for agriculture and resources extraction have shown to have, in most cases, a negative influence on the water resources and environment, in various places of the globe [1][2][3]. Often, the global and regional changes are related to improper management of the territory and to limitations in the action plans for natural resources conservation [4][5][6]. Climate change has definitely a negative impact on glaciers and water systems in many regions of the globe. In addition, the human influence on the environment contributes to a decline, in both quantity and quality, of the resources available for immediate consumption [7]. The agricultural and industrial practices conduct to higher pressure on natural resources and to many other environmental stresses (i.e., deforestation, increased water demand, pollution etc.,) [8][9][10].
Among the consequences of climate change, there are also health risks. An example of direct effect of the climate on health is the increased incidence of skin cancer [11] because of prolonged exposure to sunlight. However, most of the time the effects on health are mediated; such an example may be the increased contamination of groundwater with various pathogens because of higher frequency and intensity of floods. This results in enteric disease [12], often caused by the spread of pathogens like Campylobacter, Cryptosporidium, norovirus, or vibrio [13].
Climate change may extend the habitat of various parasites, like the liver fluke, which is a proven causal agent of cholangiocarcinoma, a highly aggressive form of liver cancer [14]. Equally dangerous, although insidious and more difficult to identify is the progressive salinization of groundwater in some inhabited areas close to sea, because of gradual sea levels rise, sometimes combined with changes in land use, like farming seawater animal species for human consumption. The higher content of salt in the drinking water has been connected with an increasing risk of arterial hypertension in the populations exposed [15,16].
Climate change influences the quality of the environment in numerous ways. Among these, a highly significant one is the mean air temperature, which influences the evapotranspiration of the hydrogeological basins, along with the land cover pattern. Starting from the beginning of this century, steady increases in the mean annual temperature were observed [17], and several models predict a continued rise in temperature during the 21st century [17][18][19][20][21][22]. The melting of glaciers represents a highly convincing proof of global warming [23][24][25][26][27][28].
Additionally, it was argued that the climate changes have a negative effect on groundwater recharge and ecosystems [29][30][31]. It was also found that the changes in land use impact on watersheds, surface water, land evapotranspiration, and groundwater [29]. Further, seawater intrusion due to the rising sea level, negatively influences the quality of the coastal fresh waters [6,[30][31][32]. Thus, the climate change effect on the water resources is divided into quantity and quality impacts. Through the evapotranspiration and water availability, the quantity impact could be assessed as part of the water balance. In the specialty literature, there are numerous methods for the evapotranspiration calculation [33][34][35][36]. Incorporating the land cover, the crop evapotranspiration (ETc) is estimated using various methods such as Penman-Monteith method and 'New Implemented Spatial-Temporal On Regions-Climate Effect on Groundwater' (NISTOR-CEGW) [37,38].
The geological formations and the type of lithology play an important role in the aquifers typology, productivity of the aquifers, intrinsic vulnerability of the aquifers, and infiltration process. Thus, the effective precipitation becomes a significant parameter that is influenced by the climate and geology of the area. Because of these reasons, this work treats in detail both climate and territory aspects with respect to the water resources and risk assessment.
The development of various regions should take into account the current threats coming from changes in climate and in land use. In order to support the preservation of the environment and create a holistic image of the geographical region, the objectives of this work are focused to calculate the main parameters that influence the water resources in the Grand Est region and to integrate the geospatial data at 1 km 2 of these parameters for long-term periods of 1990s, 2020s, and 2050s. The product is important for the assessment of climate change effects on the water resources, groundwater vulnerability determination, and future plans of water management. Through this study, the drought areas and high land use impact from the Grand Est region could be easily identified and used for advance analyses. The complexity of the survey is much consistent with integration of the terrain morphology and geology. This dataset is prepared for a risk assessment framework and it provides both spatial and temporal scales, especially for past, current, and mid 21st century.

Study Area
The Grand Est region is located in the northeastern part of France ( Figure 1) and it presents a variety of geological and geomorphological features, which are reflected by the morphology of the territory and water resources. The Vosges Mountains are located in the southeastern part of the Grand Est region and the Ardennes Mountains in the south-central and in northwest. The main rivers that cross the region are Rhine, Aube, and Marne [39]. Along these rivers are extending the largest valleys and lowlands of the region.
Atmosphere 2020, 11, x FOR PEER REVIEW 3 of 20 identified and used for advance analyses. The complexity of the survey is much consistent with integration of the terrain morphology and geology. This dataset is prepared for a risk assessment framework and it provides both spatial and temporal scales, especially for past, current, and mid 21st century.

Study Area
The Grand Est region is located in the northeastern part of France ( Figure 1) and it presents a variety of geological and geomorphological features, which are reflected by the morphology of the territory and water resources. The Vosges Mountains are located in the southeastern part of the Grand Est region and the Ardennes Mountains in the south-central and in northwest. The main rivers that cross the region are Rhine, Aube, and Marne [39]. Along these rivers are extending the largest valleys and lowlands of the region. The geology of the Grand Est region presents a wide variety of formations. Thus, the chalkstones and limestones formations in the west and central parts, the dolomites are extended in the westcentral, southern, and eastern parts, the jointed volcanic rocks, sandstones, and shales are present in the east, while the marls, clays, and sands are sparsely over the region and in the main valleys.
The region has a temperate climate with four seasons per year. The mean annual temperature ranges from 4 to 10 °C, while the precipitation regime exceeds 2000 mm/year in the Vosges Mountains [40]. The main climatic influence is coming from the Atlantic Ocean. According to the Koppen-Geiger climate classification, the region has a fully humid warm temperate climate with warm summers (Cfb class) [41].
The vegetation pattern is adapted to the relief and climate conditions of the region. The land use indicates large areas with agriculture cultivations in the eastern part, forest and natural grasslands in the center and western parts, while the main urban areas are located along the valleys, mainly in the eastern, central, northern, and western parts. The diversity of the land cover typology influences with The geology of the Grand Est region presents a wide variety of formations. Thus, the chalkstones and limestones formations in the west and central parts, the dolomites are extended in the west-central, southern, and eastern parts, the jointed volcanic rocks, sandstones, and shales are present in the east, while the marls, clays, and sands are sparsely over the region and in the main valleys.
The region has a temperate climate with four seasons per year. The mean annual temperature ranges from 4 to 10 • C, while the precipitation regime exceeds 2000 mm/year in the Vosges Mountains [40]. The main climatic influence is coming from the Atlantic Ocean. According to the Koppen-Geiger climate classification, the region has a fully humid warm temperate climate with warm summers (Cfb class) [41].
The vegetation pattern is adapted to the relief and climate conditions of the region. The land use indicates large areas with agriculture cultivations in the eastern part, forest and natural grasslands in the center and western parts, while the main urban areas are located along the valleys, mainly in the eastern, central, northern, and western parts. The diversity of the land cover typology influences with different ranges the evapotranspiration and also the contamination of the phosphorus through pollution load of each land cover type.

Climate Data
The temporal frames used in this study for the climate data are 1961-1990 (for the 1990s), 2011-2040 (for the 2020s), and 2041-2070 (for the 2050s). The annual temperature, precipitation, and ET0 data were processed at high-spatial resolution (1 km 2 ). The main reason for the climate data analysis in such an investigation on water resources and risk is that the aridity-humidity status of the area is directly correlated with the water availability. In addition, the climate data and water availability variation in time over the Grand Est region is closely related to the assessment of risk of runoff leading to flooding, and the risk of infiltration, important for civil engineering projects and for predicting landslides and slope instability. Both aridity index-here De Martonne Aridity Index (Figure 1j-l)-and water availability (Figure 1m-o) were prepared at 1 km 2 resolution for the Grand Est region. Ref. [42] prepared the monthly climate models for Europe, from where we have extracted and processed the data for this study. He has used the Parameter Regression of Independent Slopes Model (PRISM) to compute the climate models of precipitation, while the ANUSplin interpolation method was used for the temperature models computation. A moderate climate change projection was applied, with emission of Representative Concentration Pathway (RCP) 4.5, which implies a global prediction of +1.4 • C (±0.5) [43,44]. The climate models of temperature and precipitation represent an ensemble average of 15 Atmosphere Ocean Global Climate Models (AOGCMs), which were processed according to the Coupled Model Intercomparison Project 5 (CMIP5) multi-model dataset, following the Intergovernmental Panel on Climate Change (IPCC) Assessment Report 5 [45]. The climate data models were performed into ClimateEU v4.63 environment, available on the website [42].
Previous works [38,[46][47][48] described the data processing methods and the details on how the models were validated. Instead of bias corrections, a change factor (CF) method was applied to the raw GCM results. These climate models served as input for several studies and investigations at continental and regional scale [49][50][51].

Annual ETc and Annual AETc
The potential evapotranspiration (ET0) phenomena incorporating various patters of vegetation and land use is called crop evapotranspiration (ETc). In order to calculate the ETc, the ET0 in the study area was determined and a specific crop coefficient (Kc) was assigned to each land cover type. Thus, the annual ETc is calculated by multiplying the annual ET0 with the Kc (Equation (1)). For the annual ET0, the Thornthwaite formula [36] was used, while for the Kc values, the standard values were used [37]. This approach was previously used in several studies from the temperate zone [49,50]. The Kc values for the artificial and urban areas published by Ref. [51] were used in this work. Tables 1 and 2 summaries the Kc values for with respect to the land cover of the Grand Est region. The annual actual crop evapotranspiration (AETc) is the closer estimation of the real ETc, that incorporates in the calculation the annual precipitation and aridity index. For the AETc determination, the Budyko approach [52] was used because of its applicability on long-term period and large area. The Budyko formula (Equation (2)) is a useful method for the hydrological studies and water balance calculation because of the consideration of heat energy and precipitation which are taken into account [53]. Recently, numerous studies based on the Budyko equation indicate good performance [2,49]. The mathematical operations between raster layers were completed using ArcGIS environment.

Geological Data and Aquifers
Jointed and karstified rocks such as limestones, chalkstones, sandstones, and conglomerates, are representative for the geology of the Grand Est region. In addition, volcanic and plutonic rocks, and also gneisses, marls, clays, shales, mica schists, quartzites, and gravels could be found in the study area. The lithological layers (Figure 2a), together with the relief morphology, control the water flow direction and infiltration rate into the soil. This amount of infiltration is important also for the aquifers' recharge and the soil water content [54]. Thus, the typology of rocks plays a significant role in determining aquifers types, permeability, saturation of the soil, and infiltration process [55,56]. The last two processes could take varied time and are dependent on lithology and rainfall intensity.
According to the geology, six types of productivity aquifers were defined in the International Hydrogeological Map of Europe, dating from 2013, at 1:1,500,000 scale [57]. The aquifer types include fractured and porous media, and also rocks which are not considered aquifers. The main aquifers types, identified in the International Hydrogeological Map of Europe according to the productivity type are: highly productive fissured aquifers, highly productive porous aquifers, low and moderately productive fissured aquifers, low and moderately productive porous aquifers, locally aquiferous rocks-porous or fissured, and practically non-aquiferous-porous or fissured (Figure 2b).

Terrain Data and Potential Infiltration Map
The digital elevation model (DEM) of the Grand Est region was used for the slope angle and potential infiltration map calculation. The Rhine Valley in the eastern part and the Marne and Aube Valleys in the western part represent the lowlands in the Grand Est region. The highest elevations are recorded in the south-east side of the region, where the Vosges Mountains are present. The altitudes in the study area vary from 51 m to 1298 m above sea level (Figure 2e), thus the relief morphology influences greatly the runoff flow. By combining relief morphology and lithology, the groundwater flow directions and accumulation could be estimated at spatial scale.
The DEM layer is important in water resources and environmental risks studies, not only because of the influence of variations of relief topography and slopes' angle gradient on the water flow direction and infiltration control, but also because of their role in flooding, flow accumulation, and slope stability. With that in mind, using the slope angle layer (Figure 2f) together with the PIC layer, which reflects the lithology permeability, a potential infiltration map (PIM) was drawn (Figure 2g). The calculation is based on the reasoning that where the PIC is higher and the slope angle is lower, the infiltration values will be higher. This approach was developed using the ratio between PIC and slope angle (normalized values). The procedure was also applied before in studies of groundwater [9,50]. ArcGIS environment served as the software to analyze and manage the data assignation of aquifers' features in the attribute table.

Land Cover Data and Pollution Load Index
The land cover has a two-fold significance for water resources and environmental risk mapping. First, the quantitative aspect is related to the evapotranspiration, which determines the final water availability amount. Second, the qualitative aspect includes the pollution load because of phosphorus transfers, corresponding to each type of land cover. For these two aspects, the spatial distribution of crop coefficient (Kc) and pollution load index (PLI) [58] in different periods over the Grand Est region, were determined by using the CORINE Land Cover database and Hercules projections. The CORINE Land Cover database data used were from 1990 to 2012, the source being Copernicus Land Monitoring Services [59]. The Hercules projections indicate four land cover models related to 2040. These future projections were developed within the framework of "Sustainable futures for Europe's HERitage in CULtural landscapES" (HERCULES) project [60]. The methodology applied in this project is based on the macro-economic and land use models, following fourteen trajectories in the land cover changes trend [60].

Effective Precipitation
The effective precipitation represents the amount of water that infiltrates and contributes to the groundwater recharge. For the calculation of effective precipitation, the water availability and potential infiltration map were considered at spatial scale of the Grand Est region (Figure 2a-f). The precipitation, lithology permeability, soil saturation, and terrain configuration are the main drivers for the infiltration process. For this reason, all these factors were used for effective precipitation by adopting the formula from equation (Equation (4)). Being a study for a long-term period (averages of 30 years), the storage yield was neglected. where: EP, Effective precipitation (mm) WA, Water availability (mm) PIM, Potential infiltration map (dimensionless) The calculations were completed at spatial scale using ArcGIS environment and the normalized raster of PIM, following the normalization procedure [50]. The procedures of the calculation with respect to the surface and subsurface hydrologic properties, soil infiltration capacity, lithology, slope, and land use implication for the effective precipitation are in line with the previous studies [49,50].

Normalization of Data
In order to obtain a homogeneous database, all datasets were normalized between 0 and 1. Thus, the effective precipitation (EF, for 1990s, 2020s, and four scenarios for 2050s), the pollution load index (PLI) (CLC1990, CLC2012, and HERCULES Projections 2040), the potential infiltration coefficient (PIC), the hydrogeological vulnerability factor (HVF), and the potential infiltration map (PIM) were normalized and processed in a unique scale. Data normalization was done in conformity with the standard formula (Equation (5)). In this way, the datasets are ready for risk mapping analysis and other various spatial interference.
where, X, the values in series Xnorm, the value after normalization Xmax, the maximum value sin the series Xmin, the minim value in the series The data presented in this work were validated through publications [49,50] regarding the climate change effect on groundwater resources in the Grand Est region, but also through spatial analysis in ArcGIS regarding the groundwater vulnerability [50].

Climate Models of Temperature and Precipitation
The climatological maps indicate a mean annual air temperature from 4.1 • C to 10 • C (Figure 1a) (1990s), but with predicted increases up to 12.4 • C (2020s) and 13.2 • C (2050s) (Figure 1b,c). During the 1990s, the annual precipitation values range was from 568 mm to 2414 mm (Figure 1d). The higher temperatures are in the western and eastern parts of the region, mainly in the lowlands and valleys. The low temperatures were identified in the Vosges Mountains and central parts of the region. The climate models indicate annual precipitation values between 582 mm and 2452 mm during the 2020s and between 560 mm and 2442 mm in the 2050s (Figure 1e,f). The higher amount of precipitation is depicted in the Vosges Mountains and in the South-central part of the Grand Est region. The potential evapotranspiration (ET0) varied in the study area from 392 mm to 548 mm in the past period (1990s) (Figure 1g), while the current period (2020s) indicates values from 444 mm to 630 mm (Figure 1h). During the 2050s, the ET0 predicted values ranging from 468 mm to 656 mm (Figure 1i). The higher values of ET0 were identified in the western and eastern parts, while the lower values of ET0 were found in the Vosges Mountains.

Variation of the De Martonne Aridity Index
The spatial variation of the De Martonne Aridity Index in the Grand Est region during 1990s indicates extremely humid climate (index higher than 55) in the Vosges Mountains. The very humid (values from 35 to 55) and humid (values above 28 and below 35) climates were found in major parts of the territory (Figure 1j). In 2020s and 2050s, the climate will be semi-humid in some locations from east and west of the Grand East region, with lower values (below 28) of the De Martonne Aridity Index (Figure 1k-l). This means that during the 2020s and 2050s, the De Martonne Aridity Index decreases and indicates gradual aridization. In the future periods, the very humid and extremely humid climates are depicted in the mountain areas and in the central parts of the Grand Est region.

Variation of Water Availability and Evapotranspiration
The spatial distribution of water availability during the 1990s period in the Grand Est region shows values between 122 mm and 1972 mm. The higher values (above 1000 mm) were depicted in the south-central, south-eastern, and north-western areas (Figure 1m). The lower values (below 300 mm) of water availability were identified in the western and eastern parts, especially in the Marne Valley, Aube Valley, and Rhine Valley. In 2020s and 2050s, the water availability ranges from 103 mm to 1989 mm and from 96 mm to 1810 mm (Figure 1n-o). In the future periods, the high water availability values are extended in the Vosges Mountains while the low values are located in the Rhine, Marne and Aube Valleys.
Combining the climate models and land cover data, the annual ETc and AETc were calculated at spatial scale of the Grand Est region. During the 1990s, the annual ETc varied from 92 mm to 777 mm ( Figure 3a). The higher values were depicted in the western, eastern, and central parts of the region, while the lower values were found in the mountains and urban areas. In the period of 2020s, the annual ETc indicates values between 102 mm and 892 mm (Figure 3b). The future scenarios of the 2050s period show values between 138 mm and 873 mm of the ETc (Figure 3c-f). The future scenario has almost the same pattern, with few modifications that follow the land cover projections.
The spatial distribution of the annual AETc indicates values from 92 mm to 650 mm during 1990s, with high values in the central and south-eastern parts of the region (Figure 3g). The lower values of the AETc were found in the west and east, and also in the urban areas. In the 2020s, the annual AETc shows values between 100 mm and 737 mm (Figure 3j). The projections for 2050s show values between 136 mm and 723 mm. It is interesting that the AETc increases the values for the 2050s in the mountains area. For the future periods, the lower values of AETc were depicted in the valleys, where the annual precipitation is low.

Variation of Water Availability and Evapotranspiration
The spatial distribution of water availability during the 1990s period in the Grand Est region shows values between 122 mm and 1972 mm. The higher values (above 1000 mm) were depicted in the south-central, south-eastern, and north-western areas (Figure 1m). The lower values (below 300 mm) of water availability were identified in the western and eastern parts, especially in the Marne Valley, Aube Valley, and Rhine Valley. In 2020s and 2050s, the water availability ranges from 103 mm to 1989 mm and from 96 mm to 1810 mm (Figure 1n-o). In the future periods, the high water availability values are extended in the Vosges Mountains while the low values are located in the Rhine, Marne and Aube Valleys.
Combining the climate models and land cover data, the annual ETc and AETc were calculated at spatial scale of the Grand Est region. During the 1990s, the annual ETc varied from 92 mm to 777 mm (Figure 3a). The higher values were depicted in the western, eastern, and central parts of the region, while the lower values were found in the mountains and urban areas. In the period of 2020s, the annual ETc indicates values between 102 mm and 892 mm (Figure 3b). The future scenarios of the 2050s period show values between 138 mm and 873 mm of the ETc (Figure 3c-f). The future scenario has almost the same pattern, with few modifications that follow the land cover projections.
The spatial distribution of the annual AETc indicates values from 92 mm to 650 mm during 1990s, with high values in the central and south-eastern parts of the region (Figure 3g). The lower values of the AETc were found in the west and east, and also in the urban areas. In the 2020s, the annual AETc shows values between 100 mm and 737 mm (Figure 3j). The projections for 2050s show values between 136 mm and 723 mm. It is interesting that the AETc increases the values for the 2050s in the mountains area. For the future periods, the lower values of AETc were depicted in the valleys, where the annual precipitation is low.

Variation of Kc and PLI
In the Grand Est region, the land cover pattern is complex, with a multitude of vegetation type and artificial areas. The broad-leaved, coniferous, and mixed forests are located predominantly in the mountains and hilly areas. The arable and croplands extend in the lowlands and the main valleys of Rhine, Aube, and Marne. The cities and villages are spread over the region, and their areas are

Variation of Kc and PLI
In the Grand Est region, the land cover pattern is complex, with a multitude of vegetation type and artificial areas. The broad-leaved, coniferous, and mixed forests are located predominantly in the mountains and hilly areas. The arable and croplands extend in the lowlands and the main valleys of Rhine, Aube, and Marne. The cities and villages are spread over the region, and their areas are expected to increase by the mid-21 century. Each land cover type has a particular capacity for evapotranspiration [54]. Specific annual Kc values (Tables 1 and 2) for each type of vegetation and land use were assigned to the land cover layer. The method follows the Ref. [Haidu and Nistor and Mindrescu] approach and this procedure contributes to the annual ETc and AETc. Figure 4a-f shows the spatial variation of the Kc over the Grand Est region. The higher values were depicted in the agriculture and forest areas, while the lower values overlap with the industrial units, road and rail networks, and construction sites.
Atmosphere 2020, 11, x FOR PEER REVIEW 11 of 20 expected to increase by the mid-21 century. Each land cover type has a particular capacity for evapotranspiration [54]. Specific annual Kc values (Tables 1 and 2) for each type of vegetation and land use were assigned to the land cover layer. The method follows the Ref. [Haidu and Nistor and Mindrescu] approach and this procedure contributes to the annual ETc and AETc. Figure 4a-f shows the spatial variation of the Kc over the Grand Est region. The higher values were depicted in the agriculture and forest areas, while the lower values overlap with the industrial units, road and rail networks, and construction sites. Regarding the quality of the environment relative to land cover, the literature [2,55] proposed a particular PLI for each land cover type. Thus, we assigned the specific PLI for each land cover type in the ArcGIS environment. Figure 4g-l indicates the spatial variation of pollution load index in the study area. These results are highly recommended for use in the spatial analysis of water resources and other environmental risks evaluation. In the Grand Est region, the PLI has higher values in the agricultural areas, mainly in the Rhine, Marne, and Aube Valleys. The low values of PLI are associated with the pastures, sparsely vegetated areas, bare rocks, peat bogs, salt marshes, and water courses.
The values of the PLI used in the present study are reported in the Tables 3 and 4.  Regarding the quality of the environment relative to land cover, the literature [2,55] proposed a particular PLI for each land cover type. Thus, we assigned the specific PLI for each land cover type in the ArcGIS environment. Figure 4g-l indicates the spatial variation of pollution load index in the study area. These results are highly recommended for use in the spatial analysis of water resources and other environmental risks evaluation. In the Grand Est region, the PLI has higher values in the agricultural areas, mainly in the Rhine, Marne, and Aube Valleys. The low values of PLI are associated with the pastures, sparsely vegetated areas, bare rocks, peat bogs, salt marshes, and water courses.
The values of the PLI used in the present study are reported in the Tables 3 and 4.

Variation of PIC, HVF, and PIM
The PIC and HVF values are related to the geological formations and lithology, while the PIM is related to the geology and terrain morphology. Thus, based on the geological data and type of aquifers, the most important outputs, consisting in the potential infiltration coefficient (PIC) [42,46] and in aquifers' vulnerability factor were assigned [2] to the aquifers layer. In the Grand East region, the higher value of PIC are located in the areas with gravels, sands, limestones, and dolomitic formations. The lower values of PIC are mainly identified in the areas with clays and marls. The determination of PIC ( Figure 2c) at spatial scale plays a crucial role in the effective infiltration calculation [47].
The HVF layer (Figure 2d) is significant for the groundwater vulnerability mapping [50]. In the last years, the specialists in hydrogeology and environmental studies have used successfully the PIC and vulnerability factors in European aquifers into CC-WARE European Project [2] and DRINK Adria Project [61]. Table 5 reports the aquifer vulnerability factors used in the Grand Est region. The higher values of HVF are found in the dolomitic, limestones, calcarenites, and sands. The lower HVF is related to the areas with marls, clays, gneiss, mica, schists, plutonic rocks, shales, and claystones.
The PIM of the Grand Est region was generated using the DEM, slope angle, and PIC layers. High values of the PIM were depicted in the lowlands, mainly in the valleys with gravels and sands deposits. These areas show higher infiltration values because of higher PIC and low slope angle. As we can observe in Figure 4g, the areas with high values of PIM overlap with the Rhine, Marne, and Aube Valleys. The major part of the Grand Est territory appears with low infiltration rate, because of the low PIC values and multitude of slopes. Vosges Mountains, where the plutonic rocks are, represent a territory with an impermeable lithology layer.

Variation of Effective Precipitation
One of the most complex layers from the calculation point of view and with major importance for the water resources and groundwater vulnerability assessment is the effective precipitation layer. In the Grand Est region, the effective precipitation varied from 0 mm to 636 mm (1990s and 2020s). The high values were found in the Rhine Valley and along main rivers banks (Figure 5a). The lower values of effective precipitation could be depicted in major parts of the territory. During the 2050s, the effective precipitation ranges from 0 mm to about 572 mm. The higher values of the effective precipitation were found in the large valleys, but with less amount of this parameter. The significant changes between past and future could be observed in the central and north-central parts (Figure 5c

Variation of Effective Precipitation
One of the most complex layers from the calculation point of view and with major importance for the water resources and groundwater vulnerability assessment is the effective precipitation layer. In the Grand Est region, the effective precipitation varied from 0 mm to 636 mm (1990s and 2020s).
The high values were found in the Rhine Valley and along main rivers banks (Figure 5a). The lower values of effective precipitation could be depicted in major parts of the territory. During the 2050s, the effective precipitation ranges from 0 mm to about 572 mm. The higher values of the effective precipitation were found in the large valleys, but with less amount of this parameter. The significant changes between past and future could be observed in the central and north-central parts (Figure 5c Figure 6 shows the normalized maps of effective precipitation, PLI, PIC, HVF, and PIM in the Grand Est region, France. These maps represent a comprehensive gridded dataset that serves for the climate and pollution risk evaluation in the study area with respect to the water resources. The data grids are homogenous layers ready for different kind of analyses, such as spatial analysis by weights,  Figure 6 shows the normalized maps of effective precipitation, PLI, PIC, HVF, and PIM in the Grand Est region, France. These maps represent a comprehensive gridded dataset that serves for the climate and pollution risk evaluation in the study area with respect to the water resources. The data grids are homogenous layers ready for different kind of analyses, such as spatial analysis by weights, inference matrix, fuzzy logic, or any other approach that is focused on the quality and quantity assessment of nowadays and future factors. In addition, our database includes the future scenarios for the 2020s and 2050s, that could help for the environmental planning and delineation of protection zones.

Normalized Layers
Atmosphere 2020, 11, x FOR PEER REVIEW 15 of 20 for the 2020s and 2050s, that could help for the environmental planning and delineation of protection zones.

Discussion on the Usage Notes
The comprehensive gridded dataset based on GIS technology and complex approaches was set-up for the Grand Est region, France. The gridded dataset is associated with the climate models, geology, terrain morphology, and land cover and it represents a useful cartographical support for water resources analysis and environmental risk evaluation in the Grand Est region of France. The uniform representation of climate data, geological and aquifers data, terrain and potential infiltration map data, effective precipitation, land cover data and pollution load index, as well as annual crop evapotranspiration and annual actual crop evapotranspiration were combined and described using a regular grid cell of 1 km 2 .
The results indicate a negative impact of climate change on the water resources for the future periods. Analyzing the De Martonne Aridity Index variation, the index decreased during the 2020s and 2050s. This fact clearly indicates slight aridization and it suggests the negative effect of climate change in the region. The annual ETc and AETc values increase for the future, while the water availability decreases for the 2050s.
The land cover pattern indicates a large variety in the Grand Est region. Thus, the values of Kc and PLI show a wide range starting from the artificial areas, to the agriculture, pasture, and forest. For this reason, the importance of the land cover for the water resources has two-fold impact: a quantitative one, in terms of evapotranspiration and water availability amount, and a qualitative one, related to the phosphorous and load pollution variation in different kinds of land cover.
The PIC, HVF, and PIM are essential layers with respect to the substrate data and infiltration process. Thus, these three layers are useful for the aquifers recharge and groundwater vulnerability calculation.
The effective precipitation layer illustrates directly the areas with high amounts of water that infiltrate in the fractured and porous media. Considering that the geological formations do not change in a temporal scale of 100 years, the decline of the effective precipitation in the Grand Est region is due to climate change.
The originality of this study consists in providing the first regional high-resolution datasets on various patterns of distribution of the variables studied, over the Grand Est region. The work and layers prepared here are essential for water resources management and future predictions, but could be also useful for different environmental hazards, such as flooding, landslides, or droughts. These layers are appropriate files for development of the groundwater vulnerability map using spatial analyst tools in ArcGIS, on one hand, but also for the flow accumulation and flood areas determination, on other hand. This comprehensive dataset of the Grand Est region can be combined with additional datasets (e.g., groundwater table, soil properties) and used as a driver for hazard maps in the region. In this sense, TRIGRS and Scoops 3D are highly used for slope stability because of the large GIS data reading and results at spatial scale. Thus, the work presented here completes the previous studies of the Grand Est region, offering a much more complex homogenous dataset. An example may be provided by the spatial analysis, using GIS technology for groundwater vulnerability mapping, over the Grand Est region, performed by Ref. [50].
To support the development of climate change mitigation planning for water and natural resource management in the Grand Est region, and also in the scientific community, we provide the gridded data layers through the zenodo platform [62].

Conclusions
In this work, the main layers of temperature, precipitation, geology, terrain configuration, and land cover of the Grand Est region were used to analyze the annual ETc, AETc, water availability, PIC, HVF, PIM, and effective precipitation. In order to serve as a comprehensive gridded dataset for the climate change and water resources risk assessment, all these layers were processed at 1 km 2 spatial resolution and normalized. The strength of comprehensive gridded dataset processed in this study lies in the future datasets, which are of high interest for both the scientific community and the administrative authorities of the region. This region represents a study case of the relationship between climate change and agriculture. Our comprehensive gridded dataset provides the "stand alone" maps and derived maps, as well as the normalization of all layers.
The findings presented in this paper illustrate the decreasing of the De Martonne Aridity Index, which suggests negative impact of climate change on the Grand Est territory through the aridization. The annual ETc and AETc increase and the water availability and effective precipitation decrease for the future periods in the study area. The high values of the water availability could be identified in the Vosges and Ardennes Mountains, where the annual precipitation amount is higher than the evapotranspiration. On the contrary, in the mountain areas, the lithological strata are much impermeable and the infiltration is very low. For this reason, the HVF is very low in the area of Vosges and Ardennes Mountains.
The pollution load index (PLI) indicates higher values in the agricultural areas, more precise in the Rhine, Marne, and Aube Valleys and it has lower values in the natural areas such pastures and vegetated areas, bare rocks, marshes, and water courses.
Thus, the environmental planning for the near future and mid-century time period could consider the original maps that we have developed in this work. Moreover, the comprehensive gridded dataset contributes to environmental risk evaluation and to the calibration of various models (e.g., groundwater flow models, slope stability analyses). The technical validation of the presented data and methodology is supported by the previous publications [42,43], while future work may develop a more extensive dataset including hydrological models based on runoff and groundwater flow, spring discharge monitoring, and field research data including soil sampling and geotechnical investigations.