Source: https://www.youscribe.com/catalogue/documents/autres/seguimiento-del-estado-ecologico-de-las-aguas-continentales-2012758
Timestamp: 2020-04-05 15:47:55+00:00

Document:
SEGUIMIENTO DEL ESTADO ECOLÓGICO DE LAS AGUAS CONTINENTALES SUPERFICIALES ESPAÑOLAS MEDIANTE IMÁGENES MERIS (Monitoring the ecological state in Spanish inland water using MERIS images) - J. - Autres
SEGUIMIENTO DEL ESTADO ECOLÓGICO DE LAS AGUAS CONTINENTALES SUPERFICIALES ESPAÑOLAS MEDIANTE IMÁGENES MERIS (Monitoring the ecological state in Spanish inland water using MERIS images)
La teledetección tiene como uno de sus objetivos poder proporcionar información de utilidad en el menor tiempo posible para la toma de decisiones. En esta línea se han realizado multitud de estudios de teledetección en calidad de agua llegando a la conclusión de la necesidad de disponer de un sensor con bandas espectrales estrechas y centradas en las longitudes de ondas apropiadas, con una periodicidad alta y una adecuada resolución espacial. El sensor MERIS reúne casi todas estas características, siendo su resolución espacial (300 m) la más crítica. Aprovechando las facilidades que proporciona la ESA
para la realización de investigaciones con imágenes de los sensores de ENVISAT, durante los últimos años se han realizado las investigaciones necesarias y colaboraciones en todos las áreas (radiometría de campo, corrección atmosférica, desarrollo de algoritmos y software) para poder proporcionar en el menor tiempo (1 hora desde la recepción de la imagen) una cartografía fiable en calidad de agua. A continuación se presenta un resumen de todas las investigaciones así como de la operatividad del sistema.
Remote sensing has as an aim to give useful information in the shortest possible time in order to make the taking of decisions easier. This is why so much research on water quality using remote sensing has been done, the conclusion of which implies the necessity of disposing of narrow spectral bands in the appropriate wavelenghts, with a high periodicity and the right spatial resolution. MERIS has most ot these
characteristics, with a critical spatial resolution of 300 metres. Taking advantage of ESA convenience
to supply researchers with ENVISAT sensor imagery, in the last years there has been a lot of research in every area (field radiometry, atmospheric correction, algorhithm development and software) in order to produce in the shortest possible time (an hour since the reception of imagery) an accurate and trustworthy
water quality cartography. A summary of this research and the working of this system are presented below.
Revista de Teledetección. ISSN: 1988-8740. 2010. 34: 29-35
Seguimiento del estado ecológico de las aguas
continentales superficiales españolas
mediante imágenes MERIS
Monitoring the ecological state in Spanish
inland water using MERIS images
1 2 1 1 1 2J. A. Domínguez , C. Arancón , C. Alonso , A. Alonso , C. de Hoyos y A. Quesada
jose.a.dominguez@cedex.es jadomin@dfmf.uned.es
1 Centro de Estudios Hidrográficos (CEDEX). P.º Bajo Virgen del Puerto, 3. 28005 Madrid
2 Universidad Autónoma de Madrid. Departamento de Biología
Recibido el 03 de marzo de 2010, aceptado el 20 de mayo de 2010
La teledetección tiene como uno de sus objeti- Remote sensing has as an aim to give useful in-
vos poder proporcionar información de utilidad formation in the shortest possible time in order
en el menor tiempo posible para la toma de de- to make the taking of decisions easier. This is
cisiones. En esta línea se han realizado multitud why so much research on water quality using
de estudios de teledetección en calidad de agua remote sensing has been done, the conclusion
llegando a la conclusión de la necesidad de dis- of which implies the necessity of disposing of
poner de un sensor con bandas espectrales estre- narrow spectral bands in the appropriate wave-
chas y centradas en las longitudes de ondas apro- lenghts, with a high periodicity and the right
piadas, con una periodicidad alta y una adecuada spatial resolution. MERIS has most ot these
resolución espacial. El sensor MERIS reúne ca- characteristics, with a critical spatial resolution
si todas estas características, siendo su resolu- of 300 metres. Taking advantage of ESA con-
ción espacial (300 m) la más crítica. Aprove- venience to supply researchers with ENVISAT
chando las facilidades que proporciona la ESA sensor imagery, in the last years there has been
para la realización de investigaciones con imá- a lot of research in every area (field radiometry,
genes de los sensores de ENVISAT, durante los atmospheric correction, algorhithm develop-
últimos años se han realizado las investigacio- ment and software) in order to produce in the
nes necesarias y colaboraciones en todos las áre- shortest possible time (an hour since the recep-
as (radiometría de campo, corrección atmosfé- tion of imagery) an accurate and trustworthy
rica, desarrollo de algoritmos y software) para water quality cartography. A summary of this
poder proporcionar en el menor tiempo (1 hora research and the working of this system are pre-
desde la recepción de la imagen) una cartogra- sented below.
fía fiable en calidad de agua. A continuación se
presenta un resumen de todas las investigacio-
nes así como de la operatividad del sistema.
PALABRAS CLAVE: MERIS, cartografía de ca- KEY WORDS: MERIS, water quality carto-
lidad de aguas, aguas continentales, estado tró- graphy, inland water bodies, trophic state, eco-
fico, estado ecológico, Directiva Marco del Agua. logical state, Water Framework Directive.
la Planificación Hidrológica a los cambios in-INTRODUCCIÓN
troducidos en el texto refundido de la Ley de
La orden ministerial ARM/2656/2008 en la Aguas, transposición de la Directiva Marco del
cual se produce la adaptación del Reglamento de Agua 2000/60/CE que establece un marco co-
Revista de Teledetección. ISSN: 1988-8740. 2010. 34: 29-35 29J. A. Domínguez et al.
munitario en política de aguas, define la tipolo- carto su utilización para el seguimiento de pará-
gía de la masas de agua, así como los indicado- metros de calidad de agua, por lo que se decidió
res de calidad hidromorfológica, calidad físico- utilizar nivel L1b de MERIS (Peña et al., 2004).
química y presentación del estado ecológico de La metodología para implementar el sistema de
las aguas continentales superficiales. El Centro seguimiento del estado ecológico de las aguas
de Estudios Hidrográficos durante el siglo XXI continentales podemos resumirla en:
viene realizando diversos proyectos con el fin a) Generación de librería espectral de eco-
último de disponer de un sistema operacional de sistemas acuáticos continentales. Los datos a ad-
seguimiento del estado ecológico de las masas quirir en campo corresponderían a medidas de
de aguas continentales eficaz para la gestión de reflectividad, medidas de perfiles con sonda mul-
las aguas continentales españolas. Para ello era ti-paramétrica y muestras del primer espesor óp-
necesario disponer de un satélite con un resolu- tico para analizar pigmentos mediante HPLC y
ción espectral, temporal, radiométrica y espacial fitoplancton (Ruiz-Verdú et al., 2008, Domín-
adecuada. El único sensor que cumple la mayo- guez et al., 2008). La reflectividad medida des-
ría de los requerimientos necesarios para reali- de una embarcación a una determinada altura
zar este seguimiento es el sensor MERIS insta- respecto a la superficie del agua es:
lado en el satélite ENVISAT; la única deficiencia R = (L –ρ·L )/Ers sfc Sky s
de este sensor es su resolución espacial (300 m),
donde E es la irradiancia solar, L es la radian-s sfc aunque es suficiente para poder realizar el se-
cia procedente del agua y L es la radiancia pro-Sky guimiento del 80% de las masas de agua espa-
cedente del cielo.ñolas (Peña et al., 2004).
Estas variables se midieron con un espectro-
rradiómetro ASD-FR un espectralón y con un
panel de referencia lambertiano gris del 25%.METODOLOGÍA
A partir de los datos de campo (velocidad del
La idea inicial era utilizar imágenes L2, co- viento) y la información auxiliar (ángulos ceni-
rrespondientes a imágenes de reflectividad pro- tales solares) se determinan los valores de Rrs a
porcionadas por la ESA, a las cuales se les apli- partir del modelo de transferencia radiativa (Mo-
carían modelos desarrollados a partir de datos bley 1999).
de campo. La diferencia existente entre los va- La sonda multi-paramétrica consta de: CTD
lores de reflectividad de las imágenes L2 de ME- y fluorómetros de clorofila-a, materia orgánica
RIS con los valores de reflectividad de agua, en disuelta coloreada (CDOM), ficocianina y fico-
masas de aguas continentales, era tal que se des- eritrina, así como un turbidímetro y un sensor
0,09 –3[Cla] (mg m )
412 442 490 510 560 620 665 681 705 753 760 775 865 890 990
1 5 15 40 97 390
Figura 1. Espectros de campo para diferentes valores de concentración de clorofila-a.
30 Revista de Teledetección. ISSN: 1988-8740. 2010. 34: 29-35
ReflectanciaSeguimiento del estado ecológico de las aguas continentales superficiales españolas mediante imágenes MERIS
los correspondientes a la concentración de clo-
Pigmentos Grupos taxonómicos
rofila-a [Cla] y ficocianina [PC].
Peridinina Dianoflagelados La concentración de clorofila-a como indica-
Fucoxantina Diatomeas (algo en dinoflagelados) dor de biomasa y la concentración de ficociani-
Neo Euglenofíceas (algo en cloroficeas) na, pigmento característico de las cianofíceas,
Violaxantina Clorofíceas debido a que algunas de ellas son tóxicas (Peña
Alloxantina Criptofíceas
y Domínguez, 2006; Ruiz-Verdú et al., 2008).Luteina
Los algoritmos se han determinado a partir deZeaxantina Cianobacterias
datos de campo, en el mismo punto se han medi-Clorofila-b Clorofíceas (algo en uglonofíceas)ila-a Todos do espectros de agua, que posteriormente se han
Feofitina-a Todos transformado en espectros de agua correspon-
dientes a las bandas del sensor MERIS, concen-
Tabla 1. Los pigmentos analizados y grupos taxonómi- tración de clorofila-a en muestras posteriormen-
cos de algas a los que corresponden.
te determinadas en el laboratorio mediante HPLC
y medidas «in situ» mediante fluorómetro, pre-
PAR. El proceso de filtrado de las muestras pa- viamente calibrado, en el caso de la ficocianina.
ra el análisis de pigmentos fotosintéticos se re- El algoritmo de concentración de clorofila-a
aliza mediante una bomba de vacío de baja pre- en función de las bandas MERIS (Fig. 2) corres-
sión que absorbe el agua de una rampa en la cual ponde a una función por partes, necesaria debi-
se colocan los recipientes con filtros de fibra de do a que en los espectros de campo en todos se
vidrio Whatman GF/F de 25 mm de diámetro y observa el primer mínimo de absorción de la clo-
0,7 µm de tamaño medio de poro. El análisis me- rofila-a, pero no en todos se observa el segun-
diante cromatografía líquida de alta resolución do, sólo en aquellos en la concentración de clo-
(HPLC) consta de 2 fases, la fase de extracción rofila-a es elevada:
y el análisis de pigmentos mediante un equipo — Para valores altos de concentración de clo-
de cromatografía de líquidos de alta resolución rofila-a y en la cuales se observan con claridad
(HPLC), con muestreador automático y detector los valores para la cual se utiliza la banda de ab-
de Diodo Array, de la casa Hewlett-Packard (ac- sorción de la clorofila-a (banda 7) y la fluores-
tual Agilent), modelo 1050. cencia de la clorofila-a (banda 9).
b) Desarrollo de algoritmos. Los algoritmos — La segunda corresponde a aquellos valo-
que hasta el momento se han desarrollado son res de concentración de clorofila-a bajos y para
6,1257x 2y = 19,34e R = 0,72 x = {[(B5 – B2) / (B5 + B2)] – 0,662955} ≤ 0
5,2044x 2 R = 0,93 x = [(B9 – B7) / (B9 + B7)] – ≥ 0
–0,7 –0,6 –0,5 –0,4 –0,3 –0,2 –0,1 0 0,1 0,2 0,3 0,4 0,5 0,6 0,7
Figura 2. Analizando los espectros de campo se observa la necesidad de encontrar una función compuesta pero con-
tinua para la concentración de clorofila-a (p < 0,001, n = 78).
Revista de Teledetección. ISSN: 1988-8740. 2010. 34: 29-35 31
[Cla] (mg/m )J. A. Domínguez et al.
los cuales no hay respuesta espectral represen- ra esta banda se mejoró la corrección en el agua,
tativa para las longitudes de ondas de 680 nm y denominándola (SCAPE-M_B2), a partir de la
705 nm, En este caso se utilizaron la banda co- comparación de los espectros de campo y la co-
rrespondiente a la absorción de la clorofila-a, rrección SCAPE-M se determinó un procedi-
(banda 2) y el máximo de reflectividad del ver- miento de doble corrección atmosférica que con-
de (banda 5). siste en comparar la versión β de SCAPE-M y la
El algoritmo utilizado para la concentración definitiva para la banda 2 y restar este último va-
de la ficocianina [PC] en función de las bandas lor menos el primero.
MERIS (Peña y Dominínguez, 2006): [PC] = d) Validación de algoritmos en las imágenes
46678exp{5.1864·[(B9-B6)/(B9+B6)]} MERIS. La validación de los algoritmos de cloro-
2R = 0.9211 fila-a y demás pigmentos fotosintéticos se reali-
c) Corrección atmosférica de imágenes L1b. zó con la comparación de datos de campo corres-
Una buena corrección atmosférica supone una pondiente a las campañas de campo realizadas en
cartografía temática con precisión; supuesta- el 2006, para el seguimiento de la evolución tem-
mente esta correspondería al nivel L2 de ME- poral de la Perellonada (inundación de las arroza-
RIS, sin embargo la ESA no ha podido propor- les que circundan el lago de la Albufera) y calidad
cionar el procesamiento de precisión requerido de las aguas de La Albufera de Valencia (CEDEX,
(Peña et al., 2004). Esta fue la razón por la cual 2007) y las realizadas para la campaña 2007, pro-
se analizaron los diversos tipos de corrección at- yecto de un procesador de aguas continentales eu-
mosférica existentes, comparando los espectros tróficas a nivel europeo para MERIS (Koponen et
obtenidos en campo con los correspondientes a al., 2008). Los algoritmos obtenidos se aplican a
cada una de las correcciones atmosféricas utili- las imágenes L1b de MERIS con la corrección pre-
zadas en el mismo punto: el nivel L2 de MERIS, via atmosférica SCAPE-M-B2 y se comparan con
ATCOR, SCAPE-M. Se observó que la mejor los datos de campo, obteniendo unos errores cua-
corrección atmosférica correspondía a la reali- dráticos medios aceptables: RMSE[Cla] = 8,60
–3 –3zada por SCAPE-M. Sin embargo, la banda 2 co- mgm y RMSE[PC] = 12,92 mgm
rrespondiente a 442 nm, banda de absorción de e) Automatización de la cartografía temáti-
la clorofila-a que en los datos de campo corres- ca. La distribución de imágenes L1b de MERIS
ponde con un mínimo de reflectividad, intere- para los proyectos de investigación de la ESA se
sante para poder utilizarla en la determinación realiza mediante petición de la imagen a través
de valores bajos de concentración de clorofila- del programa EOLISA. Este software sólo per-
a, presentaba un error excesivo y no deseado. Pa- mite adquirir cuarto de imagen total y el tiempo
R(%)*10
413 442 490 510 560 620 665 681 708 753 762 778 865 885 900
ab060624a_Field ab060624a_ATCOR3 ab060624a_SCAPE-M_B2 ab060624a_MERIS L2
Figura 3. Comparación de los diferentes modelos de corrección atmosférica (L2, ATCOR, SCAPE-M_B2 y datos de
campo), en el punto ab060624a, correspondiente al muestreo realizado en el lago de la Albufera el 24 de junio de 2006.
32 Revista de Teledetección. ISSN: 1988-8740. 2010. 34: 29-35
R(%)*100Seguimiento del estado ecológico de las aguas continentales superficiales españolas mediante imágenes MERIS
entre la solicitud y la llegada de la imagen, que aplican los algoritmos de los pigmentos y se ob-
corresponde al cuarto de la imagen adquirida por tiene la cartografía temática.
MERIS, suele ser de unas 4 horas. Para tener la f) Generación de datos a multi-escala espa-
imagen completa hay que realizar un pedido de cial: local y estatal. La realización de los mapas
4 cuartos de imágenes. Los cuartos de imagen temáticos es el primer paso, sin embargo en la
se georreferencian a partir de los datos orbitales toma de decisiones se necesita disponer de un
y así se puede realizar el mosaico de la imagen. dato para toda la masa de agua, este valor dese-
El mosaico realizado a partir de estas imágenes ado, para poder tomar una rápida decisión, se ob-
presenta errores de continuidad. Este problema tendrá como media de los valores correspondien-
se resuelve corrigiendo los cuartos de imagen, tes a cada masa de agua. Además es interesante
con puntos de control respecto a una imagen de conocer los valores mínimo, máximo y desvia-
referencia generada a partir del mosaico de las ción típica de cada masas de agua para cada va-
imágenes TM del satélite Landsat para toda la riable. Para ello se ha diseñado y programado en
península. Una vez que las imágenes están co- IDL un programa (estadisticasporgrupo.pro),
rregidas atmosférica y geométricamente están que a partir de la imagen de cada pigmento y la
preparadas para trabajar con ellas. Los algorit- imagen de agua realiza la segmentación de la
mos se aplicaran en las masas de agua y para de- imagen de agua y determina las coordenadas ge-
terminar éstas se realiza el análisis del histogra- ográficas del centroide de cada grupo de píxe-
ma de una banda del infrarrojo cercano para les, y calcula el mínimo, máximo, media y des-
generar una máscara de agua. La banda del in- viación típica, para cada grupo de píxeles
frarrojo seleccionada es la banda 14. El análisis correspondiente a una masa de agua.
del histograma de la banda 14 se ha realizado pa-
ra más de 60 imágenes y se ha llegado a la con-
clusión de considerar que se realizará una bue- RESULTADOS
na máscara de agua en más del 99% de las
imágenes para valores inferiores al 9% de reflec- El trabajo desarrollado durante los últimos
tividad, evitando las sombras generadas por las años nos ha permitido, entre otras cosas, cono-
nubes. A las imágenes se les aplica la máscara cer las características espectrales de las masas
de agua para disponer de una imagen de agua pa- de agua continentales españolas, colaborar con
ra cada fecha. La imagen de agua tiene en cada otros investigadores para poder obtener una bue-
banda los valores correspondientes a reflectivi- na corrección atmosférica y desarrollar un sis-
dad sólo en las masas de agua, a las cuales se les tema de seguimiento del estado de la calidad de
Fecha de imagen
Superficie perellonada Superficie del lago de La Albufera
Figura 4. Seguimiento de la masa de agua del Parque de La Albufera desde
2/09/2007 hasta 22/12/2008.
Revista de Teledetección. ISSN: 1988-8740. 2010. 34: 29-35 33
Superficie inundada (Ha)
02/12/08J. A. Domínguez et al.
Fechas imágenes
3 3 3 3[Cla](mg/m )_Media [Cla](mg/m )_Max [Cla](mg/m )_Min [Cla](mg/m )_Desv
Figura 5. Seguimiento de los valores mínimo, máximo, media y desviación típica de la concentración de clorofila-a
3(mg/m ).
las aguas continentales, eficaz en un tiempo que lista para ser difundida mediante Internet. Du-
podemos considerar cuasi-real. El día 18 de mar- rante los primeros meses del 2008 se fueron re-
zo de 2009 por la tarde se pidieron a la ESA las solviendo los últimos problemas del sistema de
4 imágenes correspondientes a dicho día. A pri- seguimiento de la calidad de las aguas continen-
mera hora del día 19, la ESA nos envío el correo tales mediante teledetección y para ver su utili-
para descargar el fichero con las imágenes; 30 dad se empezó a solicitar imágenes a la ESA L1b
minutos después teníamos las imágenes y 2 ho- de MERIS desde septiembre de 2007. El análi-
ras más tarde la cartografía temática definitiva sis de las imágenes sin nubes de la península no
3[Cla](mg/m )_Media
–3Figura 6. Evolución temporal de concentración media de clorofila-a (mg · m ) en el lago de La Albufera entre (2/09/2007)
y (22/12/2008).
34 Revista de Teledetección. ISSN: 1988-8740. 2010. 34: 29-35
[Cla]_Media(mg/m )
[Cla] (mg/m )
14/01/2008 01/01/2008
24/12/2008Seguimiento del estado ecológico de las aguas continentales superficiales españolas mediante imágenes MERIS
fue tan satisfactorio como se deseaba y por lo Una vez desarrollado la metodología es nece-
tanto la periodicidad teórica de las imágenes sario su puesta en funcionamiento a un nivel de
MERIS de 3 días no siempre se cumple. La po- producción eficaz y útil para la gestión de las
tencia de la metodología desarrollada se obser- masas de aguas españolas. Para ello esta meto-
va en el seguimiento del humedal de La Albufe- dología se ha transferido al Ministerio de Medio
ra de Valencia entre septiembre de 2007 y Ambiente Rural y Marino con el objetivo de la
diciembre de 2008, una vez determinada la con- puesta en marcha del sistema de producción ade-
centración de clorofila-a se realizó el seguimien- cuado para su eficaz utilización.
to de la zona inundada y se determinó la super-
ficie del lago de La Albufera en 2.350 Ha,
comparado con la superficie determinada (2.027 REFERENCIAS
Ha) con sensores de mayor precisión (CEDEX,
2007). La ausencia de imágenes durante el mes CEDEX. 2007. Seguimiento de la evolución tempo-
de febrero del 2008 debido a la cubierta nubosa ral de la Perellonada y calidad de las aguas de La
Albufera de Valencia.no impide observar los valores bajos a primeros
DOMINGUEZ. J. A., ALONSO, C. & ALONSO, A.de marzo de 2008, correspondientes a la fase cla-
2008. Remote Sensing as a Basic Toolbox for Mo-ra que se produce en la Albufera.
nitoring Water Quality Parameters and as a System
of Surveillance of Cyanobacterial Harmful Algae
Blooms (SCYANOHABS). 2008 IEEE International
CONCLUSIONES Geoscience & Remote Sensing Symposium. Bos-
ton Massachusetts.USA.
KOPONEN, S., RUIZ-VERDÚ, A., HEEGE, T,. HE-La teledetección a través de imágenes MERIS
BLINSKI, J., SORENSEN, K., KALLIO, K.,y siguiendo la metodología expuesta en este tra-
PYHÄLAHTI, T., DOFFER, R., BROCKMANNbajo proporciona mayor información de la ma-
C. & PETERS, M. 2008. Develompent of MERISsa de agua y de la cuenca, incorporando la va-
Lake Water Algorithms, Ciudad, ESRIN Contract
riable temporal a la Directiva Marco del Agua. Nº 20436/06/I-LG,.http://www.space.tkk.fi/rese-
Los mapas obtenidos permitirían obtener car- arch/projects/water_quality_monitoring/MERIS_L
tografía de utilidad en la gestión de: akes.html. Frascati (Italia).
— Estado Trófico de las Aguas Continenta- MOBLEY, C.D. 1999. Estimation of remote sensing
reflectance from above-surface measurements. les, sencillamente delimitando los valores para
Applied Optics, 38: 7442-7455.cada nivel de eutrofia.
PEÑA, R., DOMÍNGUEZ, J. A., DE HOYOS, C. &— El estado ecológico (EQR) de las masas
RUIZ, A. 2004. Mapping of photosynthetic pig-
de aguas, ya que conocemos el nivel más bajo y
ments in Spanish reservoirs. Presented at MERIS
óptimo para cada tipología de masas de agua y user workshop http://www.esa.int/esapub/confe-
por lo tanto podemos determinar EQR al dividir rence/toc/tocSP549.pdf.
cada píxel por este valor. PEÑA, P., & DOMÍNGUEZ, J. A. 2006. Using re-
mote sensing for monitoring wetlands in Spain.— Un sistema de vigilancia de posible toxi-
Presented at GlobWetland: looking at wetlandscidad por cianobacterias con dos niveles: el ni-
from space. www.esa.int/esapub/conference/toc/ –3 –3vel-1 ([Cla] < 50 mg · m y [PC] > 30 mg · m ),
tocSP634.pdf
lugares donde hay posibilidad de encontrar cia-
RUIZ, A., SIMIS, S. G. H., DE HOYOS, C., GONS,–3nobacterias y el nivel-2 ([Cla] > 50 mg:m y H. J. & PEÑA, P. 2008. An evaluation of algo-
–3[PC]>30 mg·m ), a los que habría que ir a rithms for the remote sensing of cyanobacterial bio-
muestrear lo antes posible para analizar dichas mass. Remote Sensing of Environment, 112: 3996-
4008.muestras en laboratorio.
Revista de Teledetección. ISSN: 1988-8740. 2010. 34: 29-35 35

References: resolución 
 resolución 

resolución 
 resolución 
 resolución 
 resolución