The management of a UNESCO Global Geopark (UGGp) requires a vast wealth of miscellaneous scientific knowledge that can be successfully organised using a Geographical Information System (GIS). This paper presents a pragmatic GIS database to assist in the suitable governance of the Courel Mountains UGGp (2017) in Northwest Spain. The database is structured in 66 coverages compiled from public sources and previous works or produced through traditional mapping (combining fieldwork and photointerpretation) and GIS tools. The acquired data was later homogenised and validated by a multidisciplinary team and archived in independent coverages. Forty thematic maps illustrate the broad range of cartographic information included in the GIS database. Among them, 25 basic maps provide an overview of the UGGp and 15 new maps focus on crosscutting and technical issues. All maps illustrate the huge potential of GIS to create new resources combining coverages and adapting the legend according to their purpose and audience. The database facilitates the suitable publishing of consistent outputs (e.g., brochures, books, panels, webpages, web serves), as well as the elaboration of technical data to assist the park management. The database furnishes information on the design of education actions, touristic routes, activities and Geopark facilities. The GIS database is also a supportive tool for scientific research and provides the necessary knowledge to conduct geoconservation actions based on land use, geological hazards and the occurrence of natural and cultural heritages. Altogether, the GIS database constitutes a powerful instrument for policy-making, facilitating the identification and evaluation of alternative strategy plans.

Las Loras UNESCO Global Geopark (UGGp) is geologically diverse, particularly in relation to water-derived features: springs, karst springs, travertine buildings, waterfalls, caves. In this work, the interactions between geology, geomorphology, structures and hydrogeology are analyzed. As a result of this study, a first conceptual model of the hydrogeological functioning at Las Loras UGGp is presented. The most plausible hypothesis is that the system is formed by two superimposed aquifer systems, separated by an aquitard formed by Lower Cretaceous material. The deep lower aquifer formed by the Jurassic limestones only outcrops on the northern and southern edges of the Geopark and in a small arched band to the south of Aguilar de Campoo. It forms a basement subject to intense deformation. The upper aquifer system, formed by outcropping materials from the Upper Cretaceous, is a free aquifer. It is formed by a multilayered aquifer system that is highly compartmentalized, constituting individual moorland and lora units acting as a separate recharge–discharge system. This model explains the base level of the permanent rivers and the abundant springs, important components of the water cycle and representing a contribution to the rich geological heritage of the location.

The eruption of Cumbre Vieja (also known as Tajogaite volcano, 19 September–13 December 2021, Spain) is an example of successful emergency management. The lessons learnt are yet to be fully disclosed as is whether the response can be further improved. The latter may include tools to predict lava flow inundation rheological characteristics, amongst other issues related to volcanic eruptions (i.e., ash fall and gas emission). The aim of this study was to explore if a scientific open-source, readily available, lava-flow-modelling code (VolcFlow) would suffice for lava emplacement forecasting, focusing on the first seven days of the eruption. We only the open data that were released during the crisis and previously available data sets. The rheology of the lava, as well as the emission rate, are of utmost relevance when modelling lava flow, and these data were not readily available. Satellite lava extent analysis allowed us to preliminarily estimate its velocity, the average flow emitted, and flow viscosity. These estimates were numerically adjusted by maximising the Jaccard morphometric index and comparing the area flooded by the lava for a simulated seven-day advance with the real advance of the lava in the same timescale. The manual search for the solution to this optimization problem achieved morphometric matches of 85% and 60%. We obtained an estimated discharge rate of about 140 m³/s of lava flow during the first 24 h of the eruption. We found the emission rate then asymptotically decreased to 60 m³/s. Viscosity varied from 8 × 10⁶ Pa s, or a yield strength of 42 × 10³ Pa, in the first hours, to 4 × 10⁷ Pa s and 35 × 10³ Pa, respectively, during the remainder of the seven days. The simulations of the lava emplacement up to 27 September showed an acceptable distribution of lava thickness compared with the observations and an excellent geometrical fit. The calculations of the calibrated model required less time than the simulated time span; hence, flow modelling can be used for emergency management. However, both speed and accuracy can be improved with some extra developments and guidance on the data to be collected. Moreover, the available time for management, once the model is ready, quasi-linearly increases as the forecasting time is extended. This suggests that a predictive response during an emergency with similar characteristics is achievable, provided that an adequate rheological description of the lava is available.

The application of simulation software has proven to be a crucial tool for tsunami hazard assessment studies. Understanding the potentially devastating effects of tsunamis leads to the development of safety and resilience measures, such as the design of evacuation plans or the planning of the economic investment necessary to quickly mitigate their consequences. This article introduces a pseudo-probabilistic seismic-triggered tsunami simulation approach to investigate the potential impact of tsunamis in the southwestern coast of Spain, in the provinces of Huelva and Cádiz. Selected faults, probabilistic distributions and sampling methods are presented as well as some results for the nearly 900 Atlantic-origin tsunamis computed along the 250 km-long coast.