Metodología para la cuantificación económica de daños tangibles directos provocados por inundaciones en estaciones depuradoras de agua residual. Aplicación en el Área Metropolitana de Barcelona

Methodology for the economic quantification of direct tangible flood damage in wastewater treatment plants: Application to the Metropolitan Area of Barcelona

Guillem Flor Tey ^a1*, Eduardo Martínez-Gomariz ^a2, Beniamino Russo ^{a3, b}, Joaquín Bosque Royo ^c

^aInstitut Flumen. Universitat Politècnica de Catalunya – BarcelonaTech, Campus Nord, C/ Gran Capità 6, 08034 Barcelona, España.

^bAQUATEC (AGBAR Group), Unidad de Resiliencia y Cambio Climático, Passeig de la Zona Franca, 08028 Barcelona, España.

^cAigües de Barcelona, Empresa Metropolitana de la Gestió del Cicle Integral de l’Aigua, C/General Batet 08028 Barcelona, España.

E-mail: ^a1 guillem.flor@upc.edu, ^a2 eduardo.martinez-gomariz@upc.edu, ^a3 beniamino.russo@upc.edu, ^c joaquin.bosque@aiguesdebarcelona.cat

*Autor para correspondencia

Recibido: 22/07/2025Aceptado: 19/09/2025Publicado: 31/10/2025

Citar como: Flor-Tey, G., Martínez-Gomariz, E., Russo, B., Bosque-Royo, J. 2025. Methodology for the economic quantification of direct tangible flood damage in wastewater treatment plants: Application to the Metropolitan Area of Barcelona. Ingeniería del agua, 29(4), 263-279. https://doi.org/10.4995/ia.2025.24361

RESUMEN

El presente trabajo desarrolla una metodología para la cuantificación económica de daños tangibles directos provocados por inundaciones en Estaciones Depuradoras de Agua Residual (EDAR). El objetivo principal pasa por estimar el impacto físico y económico de las inundaciones fluviales y pluviales sobre estas infraestructuras del ciclo del agua, utilizando una metodología que combina la modelización hidrodinámica, el desarrollo de curvas de daño específicas y un proceso metódico para la estimación del impacto. Se han modelizado diferentes caudales asociados a periodos de retorno de 10, 100 y 500 años para los dominios de seis plantas de tratamiento del Área Metropolitana de Barcelona gestionadas por Aigües de Barcelona. El impacto ha sido definido mediante la combinación de la exposición, la peligrosidad y la vulnerabilidad. Los resultados muestran daños inferiores al 6% del valor total de las instalaciones para todos los casos, siendo la obra civil la tipología de activo más afectada y la EDAR de Montcada y Reixac la planta más expuesta y vulnerable.

Palabras clave | Resiliencia climática urbana; evaluación económica de daños; estaciones depuradoras de aguas residuales; sistemas de saneamiento; curvas de daño; modelización hidrodinámica.

ABSTRACT

This study presents a methodology for the economic quantification of direct tangible damage caused by flooding in Wastewater Treatment Plants (WWTP). The main objective is to estimate the physical and economic impact of fluvial and pluvial floods on these key water cycle infrastructures, using a methodological framework that combines hydrodynamic modelling, the development of asset-specific damage curves, and a systematic approach to impact assessment. Flood scenarios corresponding to 10-, 100-, and 500-years return periods have been simulated for six WWTP in the Metropolitan Area of Barcelona, all operated by Aigües de Barcelona. The impact has been assessed based on the integration of exposure, hazard, and vulnerability. Results show damages below 6% of the total asset value in all cases, with civil works being the most affected asset category and the Montcada i Reixac WWTP identified as the most exposed and vulnerable facility.

Key words | Urban climatic resilience; economic damage assessment; wastewater treatment plant; sanitation systems; damage curves; hydrodynamic modelling.

INTRODUCCIÓN

Las Estaciones Depuradoras de Agua Residual (EDAR) son infraestructuras estratégicas en las áreas urbanas modernas, esenciales para el tratamiento del agua residual. Estas instalaciones suelen ubicarse en zonas bajas inundables y cercanas a cuerpos de agua como ríos y mares (Burian et al., 2013), lo que facilita el transporte del agua por gravedad y reduce los costes energéticos, además de permitir una descarga eficiente del agua tratada en el medio natural (Friedrich and Kretzinger, 2012). Sin embargo, esta ubicación incrementa altamente su exposición a eventos climáticos extremos, especialmente inundaciones fluviales, pluviales y costeras (Tavakol-Davani et al., 2013). Las inundaciones representan una amenaza significativa para el ciclo del agua y las infraestructuras del sistema de saneamiento. Debido a ello, cada vez son más recurrentes desbordamientos, obstrucciones y roturas de la red de drenaje urbano, los colectores y las estaciones de bombeo de agua residual (Hughes et al., 2021), provocando a su vez una aumento en la contaminación ambiental del entorno (Olds et al., 2018).

Específicamente en las EDAR, las inundaciones pueden generar diversos impactos directos, daños físicos relacionados con el contacto del agua con la instalación; e indirectos, principalmente de carácter medioambiental, asociados a la interrupción del servicio de la planta. Focalizando el análisis en los impactos tangibles y directos, estos están directamente asociados al daño físico en las instalaciones, incluyendo daño en los edificios de oficinas y el equipamiento de su interior o complicaciones vinculadas a la corrosión y humedades en la obra civil (Hughes et al., 2021). Específicamente en el interior de las instalaciones, se destacan los daños relacionados con la inundación de equipos mecánicos, motores y equipamiento eléctrico en general (Spirandelli et al., 2018). Además, no son despreciables los costes derivados de la eliminación y limpieza de fangos, escombros y maleza de los sistemas, interiores de edificios y urbanización exterior de la planta (Martínez-Gomariz et al., 2019a).

Entre los años 2000 y 2023, el Consorcio de Compensación de Seguros (CCS) ha indemnizado más de 15 millones de euros en equipamientos e instalaciones del sistema de saneamiento distribuidos en 97 eventos de inundación alrededor de España, según datos de la propia organización. De todos ellos, más del 60% fueron dedicados específicamente a la cobertura de daños en EDAR, siendo una cantidad significativa que justifica un estudio específico sobre el impacto de las inundaciones a este tipo de infraestructuras. Por lo general, la evaluación de los impactos en las instalaciones de tratamiento de agua residual ha sido tratada en la literatura de forma cualitativa y poco precisa. La mayoría de investigaciones se limitan a citar y clasificar los impactos climáticos potenciales en el sistema de saneamiento (Hughes et al., 2021; Sangsefidi et al., 2023) y proponer medidas de adaptación (Burian et al., 2013). En el marco de los Planes de Gestión del Riesgo de Inundación, promocionados por el Ministerio para la Transición Ecológica y el Reto Demográfico del Gobierno de España, se desarrolló un programa piloto para la estimación del impacto económico en dos EDAR de la Comunidad de Madrid (PGRI, 2015, 2016). Este estudio cuantificó de forma simplificada el impacto de las avenidas de 100 y 500 años de período de retorno a partir de la modelización hidrodinámica mediante simulaciones en HEC-RAS 2D y el uso de curvas de daño generalistas para el sector, estableciéndose como el estudio de partida para el desarrollo de la investigación que se presenta en este documento.

El presente estudio se enmarca dentro del proyecto ICARIA (Improving Climate Resilience of Critical Assets), un proyecto europeo que tiene como objetivo la promoción del desarrollo de modelos de evaluación de impactos climáticos de forma específica para diferentes infraestructuras estratégicas de alto valor para la sociedad (Russo et al., 2023). De esta forma, la definición del impacto en este estudio se caracteriza por la estimación económica de los daños físicos en las EDAR asociados con los eventos de inundación fluvial y pluvial. La conceptualización del daño en esta investigación se define en consonancia con las especificaciones del proyecto ICARIA (Leone et al., 2025), donde se define el impacto a partir de la combinación de la peligrosidad, la exposición y la vulnerabilidad de la instalación.

El presente trabajo expone el proceso metodológico desarrollado para la cuantificación económica de los daños de las inundaciones en las EDAR y aplica la metodología para las instalaciones depuradoras del área metropolitana de Barcelona gestionadas por Aigües de Barcelona. Para ello, se estudia la peligrosidad y la exposición de cada una de las instalaciones mediante modelizaciones hidrodinámicas de las contribuciones fluviales (ríos) y pluviales (rieras) que circulan por los dominios modelados de las plantas de tratamiento. Las modelizaciones han sido realizadas mediante IBER, herramienta de simulación hidráulica bidimensional que utiliza el método de volúmenes finitos para resolver las ecuaciones de aguas poco profundas (shallow water equations) y que tiene módulos incorporados con capacidad para resolver múltiples aplicaciones en la ingeniería del agua como el transporte de sedimentos, la calidad del agua, los procesos hidrológicos a nivel de cuenca rural y urbana o el transporte de troncos, entre muchos otros (Sanz-Ramos et al., 2025).

Paralelamente, se describe la vulnerabilidad de las estaciones a partir del uso de curvas de daño que relacionan el calado de inundación con el porcentaje de daño respecto al valor total del activo, desarrolladas específicamente para esta investigación y resultando totalmente novedosas para esta tipología de activo. Además, se describe una secuenciación de pasos que permiten la evaluación del impacto en términos económicos con suficiente detalle para replicar el uso de la metodología en otras infraestructuras fuera del caso de estudio presentado. Finalmente, se presentan los resultados de esta evaluación para las instalaciones depuradoras del área metropolitana de Barcelona y se detalla la validación de la metodología a partir de los datos proporcionados por el CCS para inundaciones alrededor de España.

CASO DE ESTUDIO

El Área Metropolitana de Barcelona (AMB) constituye una de las aglomeraciones urbanas más relevantes del sur de Europa. Situada en el noreste de la península ibérica, abarca una superficie aproximada de 636 km² e integra 36 municipios que, en conjunto, concentran más de 3.2 millones de habitantes (AMB, 2024), lo que representa cerca del 43% de la población catalana en apenas el 2% del territorio autonómico. Este entorno metropolitano está marcado por un elevado grado de urbanización, incluyendo también espacios agrícolas como el delta del río del Llobregat, amplias zonas forestales en los macizos prelitorales y los dos principales ríos de la región: el Llobregat y el Besòs. La interacción entre esta morfología compleja y la urbanización ha provocado una fuerte impermeabilización del suelo, lo cual, sumado a la intensidad y brevedad de las lluvias mediterráneas, favorece una dinámica torrencial en ríos y rieras, exponiendo las zonas bajas de Área Metropolitana a situaciones de riesgo de inundación (Martínez-Gomariz et al., 2019b).

El sistema de saneamiento del AMB está estructurado en cinco subsistemas independientes y es propiedad de la administración supramunicipal. La operativa del sistema está a cargo de Aigües de Barcelona, empresa responsable de la gestión del transporte y el tratamiento del agua residual. El sistema cuenta con más de 1600 km de redes de alcantarillado y grandes colectores, 39 estaciones de bombeo y una capacidad de transporte y tratamiento superior al millón de metros cúbicos diarios (Aigües de Barcelona, 2024). Las estaciones depuradoras del AMB presentan una gran heterogeneidad tanto en capacidad como en tipología de tratamiento, reflejando la diversidad de los entornos a los que dan servicio. Las plantas principales se sitúan estratégicamente junto a los ríos Llobregat y Besòs (EDAR de Sant Feliu, EDAR del Baix Llobregat, EDAR de Montcada i Reixac y EDAR Besòs), mientras que las instalaciones de menor capacidad de tratamiento, como la EDAR de Begues y la EDAR de Vallvidrera, se sitúan en zonas naturales de menor escala y se operan mayoritariamente en remoto. Por último, la EDAR de Gavà-Viladecans se localiza en el delta del Llobregat, rodeada de rieras y canales de riego. La ubicación de las estaciones y principales ríos y rieras de la zona metropolitana se presentan en la Figura 1.

La presencia de ríos y rieras efímeras en la región responde a su orografía, que combina sistemas montañosos cercanos a la costa con una urbanización densa. Aunque no existen registros históricos significativos de daños provocados por inundaciones en las EDAR del caso de estudio, su proximidad a cauces naturales implica una exposición potencial que justifica su análisis. Queda fuera del ámbito de estudio la EDAR del Besòs, puesto que no se encuentra expuesta a este peligro climático: se ubica aproximadamente a 1.5 km del río Besòs dentro del tramado urbano de Barcelona y Sant Adrià del Besòs, y no existen rieras cercanas a la instalación. Además, el mapa de zonas inundables de la Agencia Catalana del Agua (ACA, 2024) para avenidas de 500 años de periodo de retorno tampoco muestra exposición en la zona de la depuradora.

MATERIAL Y MÉTODOS

Evaluación de la peligrosidad y la exposición: modelización en IBER

La evaluación de la peligrosidad de la inundación y la exposición de las plantas se ha realizado utilizando modelización numérica hidrodinámica bidimensional, con el objetivo de obtener mapas de inundación con información sobre diferentes variables relacionadas con la peligrosidad, específicamente la altura de agua media de inundación y el porcentaje de área inundada. Para ello se han definido seis dominios alrededor de las plantas analizadas. Para cada uno de los dominios se han realizado tres simulaciones numéricas, añadiendo las contribuciones de ríos y rieras del territorio en forma de caudales, calculados a partir de los periodos de retorno de 10, 100 y 500 años de las curvas Intensidad-Duración-Frecuencia (IDF). Este proyecto ha utilizado los resultados de diferentes modelos hidrológicos proporcionados por la Agencia Catalana del Agua (ACA, 2024), introduciendo los caudales de ríos y rieras en forma de hidrograma como condiciones de contorno en los dominios de estudio, y ha modelado posteriormente el comportamiento hidrodinámico de los cauces de las diferentes contribuciones. Para evitar la sobreestimación de la inundación debido a la superposición de los caudales máximos de cada contribución, la modelización se ha desarrollado en flujo no permanente, introduciendo los diferentes flujos de agua utilizando un hidrograma triangular simplificado, caracterizado por alcanzar el caudal pico en el tiempo de concentración de cada una de las cuencas. Se ha asumido un régimen subcrítico como condición de contorno de salida, aplicando condición de vertedero con coeficiente de descarga entre 1.2 y 1.6. Se impone condición de calado 0 en todo el dominio como condición inicial.

El estudio hidrodinámico se ha realizado utilizando IBER (Bladé et al., 2014), software de modelización numérica desarrollado en colaboración entre el Grupo de Ingeniería del Agua y Medio Ambiente (GEAMA) y el Centro Internacional de Métodos Numéricos en Ingeniería (CIMNE). Las estaciones depuradoras son amplios complejos compuestos por diferentes equipamientos y zonas, usualmente urbanizados en su interior. Esta característica favorece que la inundación en su interior no sea homogénea y sea necesario sectorizar la instalación para conocer el detalle de la mancha de agua para cada una de las divisiones. Para esta investigación se ha utilizado el módulo hidrodinámico de IBER (Sanz-Ramos et al., 2025), desarrollando un mallado no estructurado de alta resolución (tamaño de elemento entre 1 y 5 metros, en función de la planta, en las áreas de cauce principal e interior de las plantas depuradoras) que permite conocer el comportamiento del flujo del agua en el interior de las instalaciones. La elevación del terreno ha sido incorporada utilizando un modelo de elevación del terreno (DEM) con una resolución de 2x2 m proporcionado por el Institut Cartogràfic i Geològic de Catalunya (ICGC, 2024). Además, se ha refinado el mallado utilizando información catastral de los edificios del interior del dominio estudiado, incorporando las diferentes canalizaciones identificadas, y añadiendo todo tipo de detalle sobre las instalaciones interiores de la EDAR a partir de planos y proyectos constructivos facilitados por Aigües de Barcelona. Por su lado, la rugosidad del terreno ha sido incorporada en el modelo mediante la asignación automática de usos del suelo de IBER, utilizando la información de la misma base de datos del ICGC y la correlación de coeficientes de rugosidad propuesta por el Hydraulic Reference Manual from the Hydrologic Engineering Centre (Brunner, 2016), y su adaptación para el territorio catalán (ACA, 2023a). Un ejemplo de la identificación de elementos de la instalación, junto con los resultados de las tres simulaciones realizadas en la planta de Montcada i Reixac, se muestra en la Figura 2. Información sobre los dominios de la modelización y las simulaciones se muestran en la Tabla 1.

La modelización ha permitido obtener imágenes ráster con información del calado (m) en cada elemento del mallado para las diferentes simulaciones de los seis dominios estudiados. La EDAR del Baix Llobregat no queda expuesta para ningún periodo de retorno del río Llobregat gracias a un recubrimiento de escollera en el lateral del cauce del mismo río y queda excluida del estudio de impacto económico (sin exposición, el impacto es nulo). Los mapas de inundación muestran que las EDAR de Vallvidrera y Montcada i Reixac emergen como las instalaciones con una mayor exposición, debido a la influencia de las rieras de su entorno. Además, la planta de Montcada i Reixac se encuentra expuesta al desbordamiento del río Besòs para caudales de 500 años de periodo de retorno. No existen registros históricos de inundaciones en las inmediateces de las instalaciones de depuración en el área metropolitana de Barcelona. En su ausencia, se han validado los modelos hidrodinámicos de forma simple a partir de la comparativa de la mancha de inundación de los grandes ríos con los mapas de inundabilidad para los mismos periodos de retorno simulados proporcionados por la Agencia Catalana del Agua (ACA, 2024).

Evaluación de la vulnerabilidad: curvas de daño

Las curvas de daño o curvas de vulnerabilidad son herramientas que permiten relacionar la peligrosidad y la exposición con los impactos económicos, ya que relacionan las variables climáticas con el daño tangible recibido, en forma de porcentaje o coste relativo a una unidad de referencia (€/m², €/m³, etc.) (McGrath et al., 2019). Esta metodología ha sido extensamente utilizada para evaluación de daños alrededor del mundo y en especial para estudios sobre inundaciones, como sintetiza Martínez-Gomariz (2017, 2020) en su estudio sobre curvas de daño para zonas urbanas en España.

No existen curvas de daño específicas para EDAR en la literatura existe. Algunos estudios previos han desarrollado curvas para tipologías de instalaciones como ‘Almacenes’ o ‘Industrias’, como tipología de tratamiento más similar a las plantas de depuración. Además, en el manual técnico de modelos de inundación HAZUS 6.1 (FEMA, 2024) existe una pequeña propuesta de curva de daño de forma generalista para ‘wastewater facilities’ que relaciona el % de daño físico para diferentes incrementos de calados de inundación. Así pues, esta investigación propone curvas de daño específicas y novedosas para diferentes tipologías de activo que se pueden encontrar en una EDAR, como son la obra civil, los equipos mecánicos, los equipos eléctricos y la urbanización de la instalación (entendiéndola como jardinería, señalización, alumbrado o pavimentación). Estas curvas, que se muestran en la Figura 3, han sido desarrolladas basadas en el conocimiento experto de un perito especializado, contacto preferente del Consorcio de Compensación de Seguros en materia de peritaje de inundaciones, y quien ya había participado en otros proyectos de evaluación de daños a infraestructuras (Martínez-Gomariz et al., 2020), en el marco del proyecto RESCCUE (Velasco et al., 2018). Las curvas desarrolladas relacionan la altura media de inundación en contacto con el activo (m) con el porcentaje (%) de daño en relación al valor total estimado para cada tipología de activo. Con el fin de facilitar el análisis, no se han tenido en cuenta otras variables definidoras del peligro, como por ejemplo la velocidad del agua en el momento de impacto con la infraestructura. Esta simplificación es habitual en la literatura existente para el desarrollo de curvas de vulnerabilidad para otras tipologías de infraestructuras e instalaciones (Martínez-Gomariz et al., 2021).

Las curvas presentadas han sido justificadas por el perito de inundaciones y validadas utilizando datos internos de Aigües de Barcelona. Son transferibles para cualquier instalación de depuración para más de 5000 habitantes equivalentes (he), según el análisis de validación aplicado que se presenta en el capítulo de resultados. La tipología de activo de obra civil hace referencia a daños estructurales y limpieza de escombros y superficies embarradas de interiores y equipos abiertos (p. ej. decantadores), de forma similar a las curvas definidas para almacenes industriales (Martínez-Gomariz et al., 2020). Por su parte, los equipos mecánicos refieren a los trabajos de mantenimiento y reparación de motores, generadores de emergencia, turbocompresores, bombas hidráulicas, compresores de aire o equipos relacionados con el tratamiento y secado de fangos. Los equipos eléctricos hacen referencia a conexiones eléctricas, subestaciones, cableado, cuadros eléctricos y paneles de control. Finalmente, los daños relacionados con la urbanización se caracterizan básicamente en la limpieza y la retirada de fangos sedimentados en el entorno urbanizado del interior de las instalaciones.

Metodología para la cuantificación económica del impacto de las inundaciones

La cuantificación económica de los daños tangibles, directos y físicos provocados por las inundaciones en las EDAR se define a partir de una propuesta metodológica (Figura 4) compuesta por seis etapas secuenciadas, y utiliza los resultados anteriores del análisis de la exposición, la peligrosidad y la vulnerabilidad. Las diferentes etapas de la metodología se detallan en la Figura 4.

1ª ETAPA. Zonificación

Las EDAR son instalaciones de gran extensión (p. ej., la planta del Baix Llobregat cuenta con una parcela de más de 350 000 m²), característica que hace que la exposición frente a las inundaciones no sea siempre homogénea en toda la finca. Para resolver esta situación, la primera etapa del proceso corresponde a una sectorización de los edificios y equipos de la planta. Para cada una de las zonas, se conocen los equipos mecánicos principales, los edificios y los elementos eléctricos asociados a dichas instalaciones a partir de su identificación en los planos y proyectos constructivos de Aigües de Barcelona. A modo de ejemplo, la Figura 5 y la Tabla 2 muestra esta distribución para la EDAR de Gavà-Viladecans.

Para cada una de las zonas analizadas, se ha identificado también el coste estimado de construcción a partir de proyectos constructivos de las EDAR realizados históricamente. Conociendo el valor total de la construcción de las instalaciones, después de un proceso de normalización de precios entre diferentes años con el fin de poder comparar los costes, se define el % sobre el valor total de la planta para cada una de las zonas y para cada uno de las tipologías de activo analizadas (obra civil, equipos mecánicos, equipos eléctricos y urbanización), cuya distribución es habitual en los presupuestos de los proyectos constructivos para las estaciones depuradoras.

Se destaca que la definición del valor o coste de la tipología de activo ‘Urbanización’ se ha definido mayoritariamente con un 95% del valor relativo a obra civil y un 5% a equipos eléctricos. En el estudio de vulnerabilidad, la zona ZU se analiza bajo la curva de daño específica para la urbanización.

2ª ETAPA. Valorización

Obtenido el % sobre el valor total de cada una de las zonas y, en cada una de ellas, el % específico referido a las diferentes tipologías de activo, el segundo paso pasa por la estimación económica de estos porcentajes. Para ello, se ha utilizado como referencia el Pliego de Prescripciones Técnicas que rige el Contrato de Servicio de Seguro para la Cobertura de Pérdidas y Daños Materiales de la Agencia Catalana del Agua (ACA, 2023b), un documento técnico que estima el valor total de la instalación a partir de su capacidad de tratamiento. Este procedimiento ha sido efectuado para las cinco instalaciones que han presentado algún tipo de exposición a la inundación: Begues, Vallvidrera, Montcada i Reixac, Sant Feliu de Llobregat y Gavà-Viladecans.

3ª ETAPA. Peligrosidad y Exposición

Se identifican las manchas de inundación para los diferentes dominios y avenidas (periodos de retorno de 10, 100 y 500 años), que corresponden a los mapas de máximos de inundación (calado máximo alcanzado en cada una de las celdas durante toda la simulación numérica). Se extraen capas de tipología ráster (.tiff) para un post proceso en entorno SIG. Se cruzan los mapas de inundación con la localización de las zonas para cada caso.

El objetivo de la tercera fase es conocer, para cada modelización, el porcentaje de zona (en área) inundada y la altura media del agua en el interior de la sección. Este último valor corresponde al valor medio del mapa de máximos de las celdas del interior de cada subdivisión. Para obtener el % de área inundada en cada zona, se extrae el área ocupada por edificios y equipos, conociendo en este caso el área potencialmente inundable. Finalmente, se compara el área del ráster de inundación con el área potencialmente inundable para obtener el % de área inundada. Este postproceso de la evaluación de la peligrosidad se ha realizado mediante herramientas SIG.

4ª ETAPA. Vulnerabilidad

Se aplican las curvas de daño, específicas por tipología de activo y simulación para cada una de las zonas identificadas y para todas las instalaciones. Este paso permite conocer el % de daño asociado a las alturas de agua identificadas en la etapa anterior. Para las curvas de daño para la obra civil y la urbanización se utiliza el calado de inundación obtenido en el apartado anterior.

Para todas aquellas zonas donde predominan los edificios convencionales (instalaciones con paredes y puertas), se aplica un coeficiente de permeabilidad para definir el calado de inundación interior. Se asume que en episodios de inundación las puertas de los edificios están cerradas. Para todos aquellos equipos mecánicos y eléctricos que se encuentran en el interior de un edificio, se evaluará la vulnerabilidad de la tipología de activo con la altura de agua corregida mediante dicho coeficiente. La variación del coeficiente de permeabilidad utilizada en este caso corresponde a la curva propuesta para instalaciones tipo ‘Industria’ por Martínez-Gomariz (2019c, 2021) en su investigación. Este coeficiente varía entre 0 y 1 progresivamente entre calados exteriores de entre 0 y 1.5 metros. A partir de este umbral se asume que el calado interior se iguala con el exterior.

Para aquellas zonas donde no predominan los edificios convencionales (p. ej., en los decantadores), el coeficiente de permeabilidad tiene valor de 0, cuando el nivel del agua exterior no alcanza la altura del muro exterior de la infraestructura, y valor de 1, cuando este límite es superado y por lo tanto la totalidad de los equipos mecánicos y eléctricos son expuestos.

5ª ETAPA. Daño Físico

Se trata de la primera etapa de obtención de resultados. A partir de la integración de los resultados del % de daño asociado a cada tipología de activo de las zonas, se obtiene el % total de daño de cada una de las divisiones. Esta información es identificada para todas las modelizaciones realizadas y se presenta mediante mapas de impacto por zona.

Se desarrolla también una integración de resultados para la totalidad de la planta, la cual permite conocer el % de afectación total para cada simulación (general y por tipología de activo). Estos valores son relevantes ya que permiten validar el modelo con datos reales de inundaciones históricas en otras EDAR de España registrados por el CCS.

6ª ETAPA. Impacto Económico

De forma paralela a la quinta etapa, se obtiene el coste económico de la inundación en cada caso, multiplicando, para cada tipología de activo i, el % de daño asociado a cada zona (%D_Z,i) por el valor económico estimado de cada combinación (V_Z,i) obtenido en la etapa 2. Así pues, el impacto económico (IE_Z) para cada zona se define en la Ecuación (1):

$$I{E}_Z(€)=\sum _{i=1}^n\mathrm{\%}{D}_{Z,i}\times V_{Z,i}(€)\times {\mathrm{\%}}_{FA} (1)$$

Donde, i denota la tipología de activo y n las tipologías totales, y %_FA es el % de área inundada para la zona estudiada. Finalmente, se integran los daños de cada zona para conocer la afectación por planta y modelización (Figura 6). Los resultados que se obtienen en este caso son los valores absolutos de costes tangibles por evento de inundación para cada una de las plantas y simulaciones. Los resultados de esta etapa permiten comparar los costes entre tipologías de activo y entre instalaciones.

RESULTADOS Y DISCUSIONES

Los resultados de la aplicación de la metodología han permitido conocer los porcentajes de daño específicos para cada una de las zonas estudiadas. Este apartado expone los principales resultados de la implementación de la metodología en las instalaciones de Begues, Vallvidrera, Montcada i Reixac, Sant Feliu de Llobregat y Gavà-Viladecans. Inicialmente se descartó del estudio la EDAR de Besòs. Además, queda excluida de la avaluación de daños la EDAR del Baix Llobregat, debido a que el estudio hidrodinámico no muestra exposición al peligro de inundaciones para ninguna de las simulaciones realizadas.

En esta sección se presentan los mapas de daño por zona, la distribución del daño por las diferentes divisiones y sus patrones y tendencias generales, y la distribución de los daños en función de la tipología de activo analizado. Aunque el proceso ha sido aplicado a las cinco instalaciones que han mostrado exposición en alguna de las simulaciones, los resultados se ejemplificaran en la planta de Montcada i Reixac, puesto que ha resultado ser una de las instalaciones con un mayor impacto para las diferentes modelizaciones. La EDAR presenta una elevada exposición y vulnerabilidad, debido a su proximidad a cauces pluviales y fluviales cercanos, así como a la inexistencia de muros perimetrales opacos en la zona norte, principal vía de acceso del agua a la instalación.

Finalmente, se muestra una comparativa de los resultados integrados de los porcentajes de afectación total en la planta y el impacto económico absoluto en cada caso. Además, se comparan estos últimos resultados con datos históricos de inundaciones a estaciones de depuración alrededor de España (datos proporcionados por el CCS) a modo de validación de la metodología.

Mapas de daño

Los mapas de daño permiten cruzar, en una misma imagen, el mapa de inundación con la afectación en % respecto el valor total de la zona. Se presentan los datos ya agregados por las diferentes tipologías de activo, ofreciendo el % de daño total para cada sección. Se han desarrollado mapas para las 5 instalaciones expuestas, y se presenta el mapa de la EDAR de Montcada i Reixac (Figura 7) puesto que es una de las plantas que ha registrado mayor afectación.

La EDAR de Montcada i Reixac resulta impactada por la Riera de Vallençana en la parte noreste, con daños que superan el 2.5% del valor total de la planta para las zonas correspondientes la decantación primaria y secundaria, y el reactor biológico. El daño es generalizado (entre un 5 y un 10%) para la mayoría de las zonas para la inundación de 500 años de periodo de retorno, la cual es provocada por el desbordamiento del río Besòs.

Por su lado, la EDAR de Gavà-Viladecans solo muestra afectaciones que rozan el 1% de daño en la zona este de la planta, donde se encuentran las divisiones de la línea de fangos y un sencillo tratamiento terciario. Esta poca afectación es debido a que la planta se encuentra en su totalidad construida en un pequeño terraplén, reduciendo la exposición a la inundación del delta. La planta de Begues resulta también poco afectada, como consecuencia de un pequeño desbordamiento de la riera que transcurre por la parte sur de la instalación. El daño se acerca al 1% solo en el edificio de control para la inundación de 500 años de periodo de retorno.

La EDAR de Sant Feliu de Llobregat no es afectada por el río, debido a la existencia de grandes vías de comunicación viaria que separan la estación del cauce principal. Sí que se detecta la posibilidad de desbordamiento de los canales de regadío y de la zona norte y el canal-riera que rodea la planta en el borde oeste para las inundaciones de 100 y 500 años de periodo de retorno. Es significativo el daño cercano al 2% de la sección donde se encuentra el tratamiento terciario para la inundación de 500 años de periodo de retorno, siendo en este caso una zona que concentra el 20% del valor total de la planta. Finalmente, la EDAR de Vallvidrera también es considerada altamente afectada para las inundaciones de 100 y 500 años de periodo de retorno. La instalación se encuentra dentro de la Riera de Vallvidrera, lo que la hace vulnerable cuando las obras de drenaje transversal no son capaces de canalizar el caudal de la misma. Se registran afectaciones superiores al 10% del valor para las zonas de la parte este (línea de fangos, desodorización y edificio de control).

Distribución del daño por zonas

La distribución del daño total en las diferentes zonas permite identificar secciones de la planta críticas donde se concentra el impacto y entender qué edificios o equipos son los que generan un mayor daño e impacto si son expuestos a la inundación. De nuevo, se muestra como caso de ejemplo la EDAR de Montcada i Reixac en la Figura 8.

De forma generalizada para todas las instalaciones analizadas, se puede observar como los daños a la urbanización son mayoritarios para inundaciones poco significativas, como muestran los resultados para las simulaciones de 10 años de periodo de retorno. A medida que las inundaciones son mayores, y consecuentemente las alturas de agua aumentan, crecen las afectaciones en las zonas donde se concentran los equipos mecánicos y eléctricos.

Se identifican en este análisis las zonas donde se encuentran el reactor biológico, los tratamientos terciarios y las líneas de fangos como aquellas más vulnerables al peligro de inundación. Son las zonas donde se concentran mayor parte de los equipos mecánicos y suelen ser más afectadas con inundaciones de mayor periodo de retorno. Por otro lado, se identifican las zonas relacionadas con la decantación o el pretratamiento como aquellas con un menor impacto directo si son alcanzadas por la inundación. Esta clasificación permite seleccionar con mayor precisión posibles medidas de adaptación de protección para las zonas más expuestas, y para aquellas más vulnerables y con mayor valor económico.

Distribución del daño por tipología de activo

La distribución del daño entra las diferentes tipologías de activo también permite entender mejor la afectación de la inundación. La obra civil ha resultado ser la tipología de activo que sufre la mayor afectación. En la Figura 9, se ejemplifica la distribución para la EDAR de Montcada i Reixac. En este caso, la urbanización se ha integrado dentro de la tipología de obra civil, entendiéndola como una zona más y definida con un 95% de su totalidad obra civil.

La tipología de obra civil concentra el 87% del daño total de media en todas las simulaciones (media entre todas las EDAR y todas las modelizaciones). La afectación para las otras tipologías crece con las inundaciones de mayor periodo de retorno. Se destaca en este apartado la EDAR de Vallvidrera, la instalación donde se registran los mayores calados de inundación, donde la suma del daño en los equipos mecánicos y eléctricos supera el 50% para la inundación de 500 años de periodo de retorno. Esta distribución se explica como consecuencia de la mayor protección de los equipos en el interior de los edificios y la aplicación del coeficiente de permeabilidad para la aplicación de las curvas de daño en equipos mecánicos y eléctricos.

Impactos económicos totales

Finalmente se muestran los resultados integrados para la totalidad de las EDAR. Esta sección permite identificar las instalaciones que sufren de un mayor impacto. Para ello se exponen en la Figura 10 el porcentaje de daño integrado respecto el valor total de la planta, y su correspondiente valor absoluto en € totales. El valor absoluto se ha calculado a partir de la suma del daño (expresado en €) para cada una de las tipologías de activo y zonas de cada instalación.

El análisis del % de daño total muestra como las EDAR de Vallvidrera y Montcada i Reixac tienen una afectación aproximada del 2% y el 6% respectivamente para inundaciones de 100 y 500 años de periodo de retorno. Por otro lado, solo la planta de Sant Feliu de Llobregat muestra una afectación superior al 0.75% del valor total para la inundación de 500 años de periodo de retorno, como consecuencia de la exposición en esa simulación de la zona del tratamiento terciario. Todas las simulaciones muestran un % de daño total entre el 0 y el 6%, intervalo y orden de magnitud que se compra posteriormente con otras depuradoras en España afectadas por inundaciones en el apartado de validación.

Por otro lado, el estudio del impacto económico en términos absolutos identifica únicamente la EDAR de Montcada i Reixac como instalación mayormente afectada. Se esperan daños tangibles directos de 1 millón y 3.5 millones de euros para las inundaciones de 100 y 500 años de periodo de retorno respectivamente. Se estima que el valor total de la planta es de 57 millones de euros, valor 30 veces mayor que el de la planta de Vallvidrera. Por este motivo, la EDAR de Vallvidrera no es significativa en términos absolutos.

Es destacable añadir que las dos instalaciones de mayor valor no han resultado expuestas al peligro de inundación para las condiciones climáticas actuales. El valor de la EDAR del Besòs se estima en 383 millones de €, mientras que el de la EDAR del Baix Llobregat se estima en 306 millones, ambas con un valor aproximado 6 veces mayor a la instalación de Montcada (ACA, 2023b). Es trabajo futuro de esta investigación evaluar las posibles modificaciones en los caudales analizados como consecuencia del cambio climático y analizar las posibles afectaciones en las instalaciones de mayor valor.

Un aspecto relevante del estudio es su capacidad para monetizar el impacto de eventos climáticos extremos con suficiente detalle como para priorizar medidas de adaptación y protección. Esta monetización permite comparar el riesgo entre distintas plantas o entre diferentes tipologías de activo dentro de una misma instalación, facilitando una toma de decisiones más racional y basada en evidencias. Los resultados destacan que la obra civil es la tipología de activo más afectada para todas las modelizaciones, aunque en eventos más severos aumentan significativamente los daños en equipos eléctricos y mecánicos, que además son críticos para el funcionamiento de la instalación.

Validación

La metodología desarrollada en esta investigación es novedosa y muy específica para las Estaciones Depuradoras de Agua Residual. Existen limitaciones en la metodología, principalmente asociadas al desarrollo de las curvas de daño, que no han podido contar con una validación experimental y se han basado únicamente en el criterio experto. Además, no han sido registradas inundaciones históricas en las instalaciones del caso de estudio, ni se han identificado los daños en caso que hayan ocurrido en el pasado. Así pues, ha sido necesario diseñar un proceso de validación alternativo.

Para ello, el Consorcio de Compensación de Seguros (CCS) ha facilitado a esta investigación más de 200 siniestros de inundaciones en EDAR en el territorio español entre el año 2020 y 2023. Para cada uno de ellos, se ha contado con información sobre el calado medio de inundación estimado para toda la planta y el coste total de la indemnización. Se ha estimado el valor total de cada planta utilizando la misma metodología que la utilizada para las instalaciones del caso de estudio (ACA, 2023b). Así pues, se ha podido estimar el % de daño total relacionando el valor indemnizado (actualizado a costes de 2024) con el valor total de la instalación.

Paralelamente, se ha identificado la altura de agua promedio de las instalaciones analizadas para los diferentes episodios de inundación y se ha superpuesto los resultados con los datos históricos del CCS. De las 138 instalaciones de capacidad mayor a 5000 he que han sido analizadas, solamente una de ellas ha registrado un daño total mayor al 6% (un 8,8%) sobre el valor total de la planta. Todas las demás, independientemente del calado de inundación, comparten el mismo orden de magnitud que las plantas del caso de estudio, entre 0 y 6% de afectación. Las instalaciones de capacidad menor a 5000 he, usualmente plantas situadas en entornos montañosos, quedan fuera del rango de aplicabilidad de la metodología, ya que los daños calculados no coinciden con el patrón anterior. De los 121 registros de inundaciones con calados medios < 0.70 m, solo el 1.5% supera el 2.5% de daño respecto el valor total de la planta, de forma coincidente con los resultados de la metodología. Por otro lado, de los 16 registros históricos > 0.70 m, el 50% ofrecen daños entre el 2.5% y el 6% del valor total, resultados coincidentes con los 3 eventos simulados en el caso de estudio con alturas de agua medias en la instalación de 0.80, 1.00 y 1.40 m, e impactos estimados de 5.94, 2.18 y 5.88% respectivamente.

La validación confirma la solidez de la metodología propuesta para instalaciones de más de 5.000 he, ya que tanto el orden de magnitud como los valores promedio del porcentaje de daño en relación con el valor estimado de la planta son coherentes con los datos registrados en otras EDAR de España afectadas por inundaciones en las últimas dos décadas. Esta metodología constituye una primera aproximación para la identificación de daños tangibles directos en estaciones depuradoras de aguas residuales, y se plantea como una herramienta de apoyo a la toma de decisiones en la priorización de medidas de adaptación y protección frente a eventos de inundación en el Área Metropolitana de Barcelona.

CONCLUSIONES

Las inundaciones representan una amenaza creciente para las infraestructuras del ciclo urbano del agua, especialmente para las Estaciones Depuradoras de Aguas Residuales (EDAR), cuya ubicación en zonas bajas cercanas a ríos o rieras incrementa significativamente su exposición. Este trabajo ha desarrollado y aplicado una metodología específica y replicable para la cuantificación económica de daños tangibles directos provocados por inundaciones en EDAR, centrando el análisis en las instalaciones gestionadas por Aigües de Barcelona en el Área Metropolitana de Barcelona (AMB).

La metodología ha sido estructurada en seis etapas, desde la zonificación de la instalación hasta la estimación del impacto económico final. La investigación combina la modelización hidrodinámica de eventos extremos (con caudales definidos para periodos de retorno de 10, 100 y 500 años para el caso de estudio) con el uso de curvas de daño novedosas particulares para diferentes tipologías de activo (obra civil, equipos mecánicos, equipos eléctricos y urbanización), desarrolladas específicamente para este proyecto a partir de conocimiento experto.

El análisis ha puesto de manifiesto que los daños físicos en las EDAR del AMB se sitúan en una horquilla de entre el 0% y el 6% del valor total de las plantas, en consonancia con los datos históricos recopilados por el Consorcio de Compensación de Seguros (CCS) para instalaciones de más de 5000 he en toda España entre los años 2000 y 2023. Este rango de afectación da solidez a la propuesta metodológica, que se perfila como una herramienta útil para su aplicación en otras regiones y contextos. La EDAR de Montcada i Reixac, al encontrarse expuesta al desbordamiento del río Besòs y a la Riera de Vallençana, ha sido la instalación con mayores daños estimados en términos tanto relativos como absolutos (entre 1 y 3,5 millones de euros para las modelizaciones de 100 y 500 años de periodo de retorno, respectivamente). En contraste, otras EDAR como Gavà-Viladecans, Begues o Sant Feliu de Llobregat muestran impactos significativamente menores gracias a su ubicación, protección natural o diseño constructivo.

En definitiva, este trabajo aporta una herramienta metodológica innovadora para la gestión del riesgo por inundación en las EDAR, integrando elementos técnicos, económicos y espaciales. Su aplicación permite no solo dimensionar los impactos físicos y económicos asociados a las inundaciones, sino también orientar la planificación estratégica y las inversiones necesarias entorno a la resiliencia climática. Como trabajo futuro, se propone la replicación de la misma metodología bajo distintos escenarios de cambio climático, considerando la posible intensificación de caudales y eventos extremos, así como incorporar mecanismos de validación adicionales a partir de la aplicación de la misma metodología en estaciones de tratamiento donde se dispongan suficientes datos para replicar el proceso y se cuente con registros históricos de inundaciones y daños asociados.

AGRADECIMIENTOS

Los autores agradecen a Aigües de Barcelona, Empresa Metropolitana de la Gestió de l’Aigua, el soporte logístico y la provisión de datos e información. Los autores agradecen el apoyo del proyecto ICARIA (Improving Climate Resilience of Critical Assets), proyecto financiado por la Comisión Europea a través del Horizon Europe Programme (Grant Number 101093806). Los autores agradecen a Salvador Castán, perito especialista en inundaciones, el soporte en el desarrollo de las curvas de daño. Los autores agradecen al Consorcio de Compensación de Seguros, y a Juan Manuel Morales y José Narciso Pigem, de Vitaqua SA, su colaboración facilitando datos de siniestros de inundación en España. El primer autor también agradece el apoyo del programa de Doctorats Industrials de la Generalitat de Catalunya.

REFERENCIAS

ACA (Agència Catalana de l’Aigua). 2023a. Recomanacions tècniques per als estudis d’inundabilitat d’àmbit local. Guia tècnica. Agència Catalana de l’Aigua, Barcelona.

ACA (Agència Catalana de l’Aigua). 2023b. Plec de prescripcions tècniques que regeix el contracte de serveis d’assegurances per a la cobertura de les pèrdues o danys materials de l’obra civil de l’Agència Catalana de l’Aigua (CTN2300486). Departament d’Acció Climàtica, Alimentació i Agenda Rural, Spain.

ACA (Agència Catalana de l’Aigua). 2024. Cabals normalitzats, innundabilitat, espais fluvials, capes d’informació geogràfica [Acceso 11 de febrero de 2024]. Visor Cartogràfic Agència Catalana de l’Aigua General, Catalunya. https://sig.gencat.cat/visors/VISOR_ACA.html

Aigües de Barcelona. 2024. Memòria d’explotació 2024. Aigües de Barcelona, Barcelona.

AMB (Àrea Metropolitana de Barcelona). 2024. Getting to know the metropolitan area. Location and land uses (2024) [Acceso 11 de febrero de 2024]. https://www.amb.cat/en/web/area-metropolitana/coneixer-l-area-metropolitana

Bladé, E., Cea, L., Corestein, G., Escolano, E., Puertas, J., Vázquez-Cendón, E., Dolz, J., Coll, A. 2014. Iber: herramienta de simulación numérica del flujo en ríos. Revista Internacional de Métodos Numéricos para Cálculo y Diseño en Ingeniería, 30(1), 1–10. https://doi.org/10.1016/j.rimni.2012.07.004

Brunner, G.W. 2016. HEC-RAS River Analysis System: Hydraulic Reference Manual. U.S. Army Corps of Engineers, Hydrologic Engineering Center (HEC), Davis, California.

Burian, S.J., Walsh, T., Kalyanapu, A.J., Larsen, S.G. 2013. Climate vulnerabilities and adaptation of urban water infrastructure systems. In Climate vulnerability: Understanding and addressing threats to essential resources (pp. 87–107). Elsevier Inc. https://doi.org/10.1016/B978-0-12-384703-4.00509-8

FEMA (Federal Emergency Management Agency). 2024. HAZUS 6.1 Flood Model Technical Manual. U.S. Federal Emergency Management Agency.

Friedrich, E., Kretzinger, D. 2012. Vulnerability of wastewater infrastructure of coastal cities to sea level rise: A South African case study. Water SA, 38(5), 755–764. https://doi.org/10.4314/wsa.v38i5.15

Hughes, J., Cowper-Heays, K., Olesson, E., Bell, R., Stroombergen, A. 2021. Impacts and implications of climate change on wastewater systems: A New Zealand perspective. Climate Risk Management, 31, 100262. https://doi.org/10.1016/j.crm.2020.100262

ICGC (Institut Cartogràfic i Geològic de Catalunya). 2024. Bessons digitals – Elevacions: Model d’elevacions del terreny de 2x2 m [Acceso 22 de febrero de 2024]. Generalitat de Catalunya. https://www.icgc.cat/ca/Geoinformacio-i-mapes/Dades-i-productes/Bessons-digitals-Elevacions

Leone, M.F., Zuccaro, G., De Gregorio, D., Turchi, A., Tedeschi, A., Bügelmayer-Blaschek, M., Sfetsos, A., Zarikos, I., de la Cruz Coronas, A., Russo, B. 2025. A holistic asset-level modelling framework for a comprehensive multi-hazard risk/impact assessment: Insights from the ICARIA project. International Journal of Disaster Risk Reduction, 119, 105319. https://doi.org/10.1016/j.ijdrr.2025.105319

Martínez-Gomariz, E., Gómez, M., Russo, B., Sánchez, P., Montes, J.A. 2017. Metodología para la evaluación de daños a vehículos expuestos a inundaciones en zonas urbanas. Ingeniería del Agua, 21, 247. https://doi.org/10.4995/ia.2017.8772

Martínez-Gomariz, E., Gómez, M., Russo, B., Sánchez, P., Montes, J.A. 2019a. Methodology for the damage assessment of vehicles exposed to flooding in urban areas. Journal of Flood Risk Management, 12. https://doi.org/10.1111/jfr3.12475

Martínez-Gomariz, E., Locatelli, L., Guerrero, M., Russo, B., Martínez, M. 2019b. Socio-economic potential impacts due to urban pluvial floods in Badalona (Spain) in a context of climate change. Water (Switzerland), 11. https://doi.org/10.3390/w11122658

Martínez-Gomariz, E., Guerrero-Hidalga, M., Russo, B., Yubero, D., Gómez, M., Castán, S. 2019c. Desarrollo y aplicación de curvas de daño y estanqueidad para la estimación del impacto económico de las inundaciones en zonas urbanas españolas. Ingeniería del Agua, 23, 229. https://doi.org/10.4995/ia.2019.12137

Martínez-Gomariz, E., Forero-Ortiz, E., Guerrero-Hidalga, M., Castán, S., Gómez, M. 2020. Flood depth-damage curves for Spanish urban areas. Sustainability (Switzerland), 12. https://doi.org/10.3390/su12072666

Martínez-Gomariz, E., Forero-Ortiz, E., Russo, B., Locatelli, L., Guerrero-Hidalga, M., Yubero, D., Castán, S. 2021. A novel expert opinion-based approach to compute estimations of flood damage to property in dense urban environments: Barcelona case study. Journal of Hydrology, 598, 126244. https://doi.org/10.1016/j.jhydrol.2021.126244

McGrath, H., Abo El Ezz, A., Nastev, M. 2019. Probabilistic depth–damage curves for assessment of flood-induced building losses. Natural Hazards, 97, 1–14. https://doi.org/10.1007/s11069-019-03622-3

Olds, H.T., Corsi, S.R., Dila, D.K., Halmo, K.M., Bootsma, M.J., McLellan, S.L. 2018. High levels of sewage contamination released from urban areas after storm events: A quantitative survey with sewage specific bacterial indicators. PLoS Medicine, 15, e1002614. https://doi.org/10.1371/journal.pmed.1002614

PGRI (Planes de Gestión del Riesgo de Inundación). 2015. Informe diagnóstico de la situación del riesgo de inundación de las EDAR de Butarque (Término municipal de Madrid). Ministerio para la Transición Ecológica y el Reto Demográfico, Gobierno de España, Madrid.

PGRI (Planes de Gestión del Riesgo de Inundación). 2016. Informe diagnóstico de la situación del riesgo de inundación de las EDAR de La China (Término municipal de Madrid). Ministerio para la Transición Ecológica y el Reto Demográfico, Gobierno de España, Madrid.

Russo, B., de la Cruz Coronas, À., Leone, M., Evans, B., Brito, R.S., Havlik, D., Bügelmayer-Blaschek, M., Pacheco, D., Sfetsos, A. 2023. Improving climate resilience of critical assets: The ICARIA Project. Sustainability (Switzerland), 15, 14090. https://doi.org/10.3390/su151914090

Sangsefidi, Y., Bagheri, K., Davani, H., Merrifield, M. 2023. Data analysis and integrated modeling of compound flooding impacts on coastal drainage infrastructure under a changing climate. Journal of Hydrology, 616, 128823. https://doi.org/10.1016/j.jhydrol.2022.128823

Sanz-Ramos, M., Sañudo, E., López-Gómez, D., García-Feal, O., Bladé, E., Cea, L. 2025. Evolución de la modelización numérica bidimensional del flujo en lámina libre a través del software Iber. Ingeniería del Agua, 29, 114–131. https://doi.org/10.4995/ia.2025.23259

Spirandelli, D., Babcock, R., Shen, S. 2018. Assessing the vulnerability of coastal wastewater infrastructure to climate change. University of Hawai‘i Sea Grant College Program.

Tavakol-Davani, H., Burian, S., Apul, D., Devkota, J. 2013. Developing urban water infrastructure modeling approach to control the flood damages and reduce life cycle impacts. In World Environmental and Water Resources Congress, Cincinnati, Ohio.

Velasco, M., Russo, B., Martínez, M., Malgrat, P., Monjo, R., Djordjevic, S., Fontanals, I., Vela, S., Cardoso, M.A., Buskute, A. 2018. Resilience to cope with climate change in urban areas: A multisectorial approach focusing on water – The RESCCUE project. Water (Switzerland), 10, 1356. https://doi.org/10.3390/w10101356