un modelo hidrodinÁmico-barotrÓpico para los … marin.pdf · son más recientes (salomons et...

Modelo Hidrodinámico - Barotrópico para fiordos australes 125

UN MODELO HIDRODINÁMICO-BAROTRÓPICO PARA LOSFIORDOS AUSTRALES DE CHILE ENTRE LOS 41° S Y LOS 46° S*

A HYDRODYNAMIC-BAROTROPIC MODEL FOR THEAUSTRAL FJORDS IN CHILE BETWEEN 41° S AND 46° S

VÍCTOR H. MARÍN1

FRANCISCO J. CAMPUZANO2

1Laboratorio de Modelación Ecológica,Departamento de Ciencias Ecológicas,

Facultad de Ciencias, Universidad de Chile;Casilla 653, Santiago, Chile.2Instituto Superior Técnico,Pavilhão de Mecânica I, 1º,

Av. Rovisco Pais, 1049-001 Lisboa,E-mail: [email protected]

Recepción: junio de 2006 – Versión corregida aceptada: julio de 2008

RESUMEN

La modelación de los ecosistemas estuarino-costeros característicos de la región de los fiordosaustrales de Chile debería realizarse teniendo como base modelos hidrodinámicos que den cuenta dela circulación y sus principales forzantes. En este trabajo presentamos un modelo barotrópico, for-zado por la marea (FIORDOS_NIVEL1), para la zona comprendida entre los 41° S y los 46° S,generado por medio del sistema numérico MOHID. FIORDOS_NIVEL1 puede ser usado comobase para la generación de modelos anidados de menor dimensión y mayor resolución en ámbitosdiversos como estudios de calidad de agua o simulaciones ecosistémicas.

Palabras claves: Modelo numérico, circulación, fiordos, MOHID.

ABSTRACT

Modelling coastal-estuarine ecosystems, characteristic of the southern fjords of Chile, shouldbe based on hydrodynamic models accounting for the water circulation and its forcing. In this workwe describe a barotropic, tidally forced, model (FIORDOS_NIVEL1) for the region between 41° Sand 46° S, based on the MOHID modelling system. FIORDOS_NIVEL1 can be used as the basis ofnested, smaller dimension and higher resolution, models for a variety of studies such as water qualityand ecosystem simulations.

Key words: Numerical model, circulation, fjords, MOHID.

* Proyecto CONA-C9F 03-06.

Cienc. Tecnol. Mar, 31 (2): 125-136, 2008

126 Revista Ciencia y Tecnología del Mar, Vol. 31 (1) - 2008

INTRODUCCIÓN

El desarrollo de modelos cuantitativos esimportante en la mayoría de las ciencias natu-rales. Estos pueden ir desde simples modelosempíricos (Marín et al., 2003) hasta modelosanalíticos (Caswell, 2001) y numéricos (He-aps, 1987). Los modelos pueden servir a unavariedad de objetivos tales como: síntesis delconocimiento sobre el sistema bajo estudio;hacer evidentes y explícitas las falencias endicho conocimiento; como plataformas con-ceptuales para poner a prueba hipótesis sobrela dinámica de los sistemas y como herramien-tas de predicción y construcción de escenariosde desarrollo futuro (Canham et al., 2003;Meadows et al., 2004). En el caso de la ocea-nografía en general, y la oceanografía costeraen lo particular, existe una larga tradición en lageneración y uso de modelos numéricos (Ni-houl, 1975; Bowden, 1983; Heaps, 1987;Dyke, 1996). Sin embargo, modelos acopladosfísico-ecológicos y físico-ecológico-socialesson más recientes (Salomons et al., 1999; Penget al., 2002; Arrequin-Sanchez et al., 2004;Fennel & Neumann, 2004). Estos últimos sonnecesarios para generar predicciones respectode los efectos potenciales de las sociedadeshumanas en los ecosistemas costeros, para po-ner a prueba hipótesis respecto de mecanismosde funcionamiento de estos ecosistemas y pararealizar simulaciones respecto de acciones demanejo. Lo anterior es especialmente válidobajo el esquema del manejo integrado de zo-nas costeras o en la elaboración de normas se-cundarias sobre la calidad de las aguas mari-nas (Salomons et al., 1999).

El proyecto ECOManage, financiado por laUnión Europea1 y desarrollado en Chile por laUniversidad de Chile2 y el Centro de EcologíaAplicada, tuvo como objetivo principal am-pliar la capacidad de los manejadores de re-cursos para el uso de un sistema integral deconocimiento basado en la implementación demodelos numéricos de ecosistemas coste-ros. El proyecto se desarrolló en tres áreasgeográficas de Latinoamérica: (1) el Estua-rio Santos en Brasil, (2) Bahía Blanca enArgentina y (3) el Fiordo Aysén en Chile.La zona de los fiordos australes de Chile espor estos días escenario de un acelerado pro-

ceso de inversiones y proyectos de desarrollo(e.g. salmonicultura, centrales hidroeléctricas)los cuales pueden generar modificaciones tan-to en la calidad del agua como en la ecologíade los fiordos. Sin embargo, aún cuando sehan desarrollado modelos numéricos sobre laCorriente de Humboldt (Batten et al., 1995;Marín, 1997; Escribano et al., 2004; Palma etal., 2006), no existe un modelo numérico de lacirculación para la zona de los fiordos. ECO-Manage ha generado una serie de modelos ani-dados de circulación, con el propósito de ana-lizar procesos ecológicos e influencias antro-pogénicas en el interior del Fiordo Aysén yBahía Chacabuco. La circulación es modeladapor medio de tres modelos anidados. El demayor dimensión (FIORDOS_NIVEL1) es unmodelo barotrópico forzado por la marea y queabarca la zona comprendida entre los 41°S ylos 46°S y entre la costa y los 76,5° W. El prin-cipal objetivo de este trabajo es describir estemodelo con el propósito de ponerlo a disposi-ción de la comunidad científica como platafor-ma para estudios y desarrollo de otros mode-los en zonas al interior de sus límites.

MÉTODOS

1. El sistema de modelación MOHID

MOHID3 es un sistema numérico para lamodelación tridimensional de ecosistemasacuáticos. Es desarrollado y mantenido por elCentro de Investigaciones Marinas y de Tec-nología Ambiental (MARETEC) del InstitutoSuperior Técnico perteneciente a la Universi-dad Técnica de Lisboa en Portugal. El sistemaestá compuesto de herramientas de pre-proce-samiento (e.g. MOHID GIS), de una interfasegráfica para la implementación de los modelos(MOHID GUI) y de herramientas para el post-procesamiento de los resultados (e.g. MOHIDGIS, MOHID POST, MOHID Time Series Edi-tor y MOHID Statistical Analyzer). Desde elpunto de vista de la modelación de ecosiste-mas marinos, MOHID está compuesto por másde 40 módulos los cuales interactúan y se aco-plan para el desarrollo de modelos en un vastocampo de aplicaciones. El sistema puede serinstalado en un Laptop PC con un procesadorde 1,5 GHz y una memoria RAM de 1 Gbyte o

1 http://www.ecomanage.info2 http://antar.uchile.cl3 http://www.mohid.com


más. Sin embargo, modelos anidados, espe-cialmente aquellos que requieren un modelo deNIVEL1 extenso (ver más abajo), deberíanejecutarse en estaciones de trabajo. En el casodel Laboratorio de Modelación Ecológica,MOHID corre en una estación de trabajoDELL Precision 670 con dos procesadores In-tel Xeon de 3,6 GHz y con 2 Gbyte de memo-ria RAM. MOHID, como sistema informático,es de libre acceso (freeware) y de código abier-to.

El módulo hidrodinámico de MOHID, co-rresponde a un modelo baroclínico 3D, imple-mentado en volúmenes finitos (Martins et al.,2001), que resuelve las ecuaciones primitivasincompresibles, suponiendo equilibrio hidros-tático y empleando la aproximación de Bous-sinesq. La viscosidad vertical es calculadapor medio del modelo General Ocean Tur-bulence Model (GOTM4). Una descripciónmás detallada del módulo hidrodinámico sepuede encontrar en Martins et al. (2001),Coelho et al. (2002) y en los manuales dis-ponibles en el portal Internet de MOHID(http://www.mohid.com).

MOHID permite la modelación de áreascosteras y estuarinas de compleja batime-tría y topografía, a través del anidamientode una vía (Braunsweig et al., 2004). Esteanidamiento implica que para cada nivel, lasolución de referencia para las condicionesde borde abierto corresponde al nivel jerár-quico superior. MOHID, desde la perspec-tiva de su desarrollo como software, no tie-ne limitaciones respecto de los niveles deanidamiento siendo solo limitado por lamemoria RAM y la velocidad del procesa-dor. Sin embargo, la experiencia muestraque para la mayoría de las aplicaciones tresniveles de anidamiento son suficientes.Para los propósitos de la modelación de losfiordos australes hemos implementado unsistema de tres niveles. El primer nivel,FIORDOS_NIVEL1, corresponde a una so-lución barotrópica 2D de la corriente demarea para la zona comprendida entre los41º S y los 46º S. Al interior de este nivelse han implementado en orden jerárquicodescendiente otros dos modelos: un modelo3D ba roc l ín i co pa ra e l F io rdo Aysén

(AYSÉN_NIVEL2) y un modelo más detalla-do (3D, baroclínico) para la Bahía de Chaca-buco (CHACABUCO_NIVEL3) al interior delfiordo. Este trabajo solo se concentra en el pri-mer nivel.

2. Implementación del modelo FIORDOS_NIVEL1

Batimetría

La batimetría, especialmente para zo-nas de geografía compleja como los fior-dos australes, es una variable importanteen modelos numéricos de circulación. Labatimetría del modelo FIORDOS_NIVEL1se generó a partir de tres fuentes princi-pales: (a) cartas náuticas provenientes delServicio Hidrográfico y Oceanográfico dela Armada de Chile, (b) la base de datosgeofísicos GEODAS (GEOphysical DataSystem; NOAA´s National GeophysicalData Center) y la base de datos GEBCO(Genera l Ba thymet r i c Char t o f theOceans; IOC, IHO and BODC, 2003). Enel primer caso, se digitalizó la informa-ción proveniente de las cartas náuticasdisponibles para zona de estudio por me-dio de ARCVIEW GIS 3.3 (Fig.1). Porotra parte, la base de datos GEODAS con-tiene información batimétrica, sísmica,magnética y de gravedad proveniente delos bitácoras de navegación de crucerosmarinos desde 1953 a la fecha (GEODAS,2006). La distribución espacial de los da-tos GEODAS usados en e l modelo semuestra en la Figura 2. Por último, la basede datos GEBCO consiste en una colec-ción de datos batimétricos que cubre todoel planeta con una resolución de un minu-to. La principal ventaja de esta base dedatos es su amplia cobertura, sin embargopresenta el inconveniente que en zonassomeras cercanas a la costa, estos datosno son fiables, por lo que las otras dosfuentes de datos complementan perfecta-mente a esta última. A partir de este con-junto de datos se generó, por medio del pro-grama MOHID GIS, la grilla numérica de ba-timetría. La resolución espacial para el mode-lo FIORDOS_NIVEL1 es de 0,02º x 0,02º (2,2km aprox.) lo que genera una grilla de 282 x212 elementos (Fig. 3).

4 http://www.gotm.net


Condiciones de borde para la frontera oceáni-ca

Los principales forzantes de la circulaciónsuperficial en la zona de fiordos son las co-rrientes de marea y las descargas de aguadulce desde las cuencas hidrográficas (Sa-linas & Hormazábal, 2004). El modeloFIORDOS_NIVEL1 es forzado en su fron-tera oceánica por la marea, incluyendo la con-dición de radiación descrita por Blumberg &Kantha (1985). Los componentes de la mareason generados a partir del modelo global

FES2004, que corresponde a una versión ac-tualizada del FES95.2 (Le Provost et al.,1998). Los componentes son impuestos en es-taciones ubicadas en la frontera oceánica a par-tir de las cuales MOHID genera las triangula-ciones para imponer la marea en cada una delas celdas (Fig. 3).

Datos usados para la validación del modelo

El modelo FIORDOS_NIVEL1 fue valida-do para los niveles de marea respecto de esta-ciones disponibles entre los 41° S y los 46° S.

Fig. 1: Cobertura de puntos de batimetría digitalizados a partir de las cartas náuticas del ServicioHidrográfico y Oceanográfico de la Armada de Chile.

Fig. 1: Bathymetric point coverage digitized from the nautical charts of the Hydrographic and OceanographicService of the Chilean Navy.


Se tuvo acceso a dos tipos de datos: (1) esta-ciones mareográficas permanentes, a cargo delServicio Hidrográfico y Oceanográfico de laArmada de Chile y (2) mediciones de cortoperíodo y alta resolución temporal (10 minu-tos) realizadas durante el crucero CIMAR 4Fiordos. La Tabla I muestra el nombre, la loca-lización y el período de observaciones para lasestaciones de marea. Para el caso de las esta-ciones permanentes del SHOA, se tuvo accesoa datos horarios para los años 2003 y 2004. LaFigura 3 muestra la ubicación de las principa-les estaciones mareográficas.

RESULTADOS

El análisis de los valores de la altura de lamarea para las distintas estaciones controlmostró que el modelo FIORDOS_NIVEL1captura de buena forma las fluctuaciones delnivel del mar en la zona de estudios. La Figura4 muestra un ejemplo de la relación entre losdatos de cuatro estaciones mareográficas (va-lores obtenidos a partir de las componentes demarea por análisis armónico realizado con elsoftware TASK-2000 POL/PSMSL TidalAnalysis Software Kit 2000; Bell et al, 1999)

Fig. 2: Batimetría GEODAS disponible para el área del modelo FIORDOS_NIVEL1 (GEODAS, 2006).Fig. 2: Available GEODAS bathymetry for the FIORDOS_NIVEL 1 model (GEODAS, 2006).


Fig. 3: Posición de las estaciones donde se genera el forzamiento de la marea en el borde oceánico (cuadradosrojos) sobre la batimetría final del modelo FIORDOS_NIVEL 1 (resolución espacial de 0,02º x 0,02º.)En cada estación se impone la solución del modelo de mareas FES2004. La señal de marea para cadagrilla de la frontera se genera por medio de la triangulación de las señales en las estaciones de marea.Los puntos geográficos identificados con nombres corresponden a las principales estacionesmareográficas usadas para validar los resultados de nivel del mar del modelo FIORDOS_NIVEL1.

Fig. 3: Position of the stations where the tidal forcing is generated in the oceanic boundary (red squares)superimposed to the final bathymetry of the FIORDOS_NIVEL 1 model (spatial resolution 0.02° x0.02°). The solution of the FES2004 tidal model is imposed in each station. The tidal signal for eachboundary grid is generated triangulating the signal from all stations. Geographic names identified bynames correspond to the main tide-gauge stations used to validate the sea level results from theFIORDOS_NIVEL 1 model.


y los predichos por el modelo. La bondad delajuste, expresada en términos del coeficientede determinación (r2) fluctuó entre un mínimode 0,92 para Ancud y un máximo de 0,99 paraPuerto Cuptana. De especial interés, para lospropósitos de desarrollo de modelos anida-dos de menor escala, resulta el buen ajuste(r2= 0,94) logrado entre la altura de marea delmodelo y los valores medidos en Bahía Cha-cabuco al interior del Fiordo Aysén (Fig. 3).La condición interna de esta Bahía unido alhecho que está localizada al Sur-Oeste de laConstricción de Meninea (la principal restric-ción topográfica del Canal Moraleda) sirvecomo una buena indicación para el uso delmodelo FIORDOS_NIVEL1 como base parafuturos modelos anidados.

La dinámica de la marea en la zona se pue-de apreciar gráficamente en la Figura 5. Estacorresponde a seis imágenes (cada 2 horas) delmodelo correspondiente a 12 horas para lacondición de sicigia. La onda de marea aúncuando ingresa a los fiordos a través de todoslos canales que se comunican con el océano,lo hace preferentemente por la Boca del Guafodesviándose principalmente hacia el Norte.Ello, unido al ingreso de la onda a través delCanal de Chacao, genera las altas mareas ob-servadas en la zona de Puerto Montt. Duranteel proceso, la onda viaja hacia el sur a travésdel Canal de Moraleda y a pesar de que granparte del agua abandona el sistema por la bocadel Guafo, una parte viaja hacia el sur unién-

dose al agua que ingresa a través de los cana-les que conectan el Canal Moraleda con elOcéano Pacífico. Como resultado se obtieneuna intensificación de la señal de marea en elextremo sur del canal (imágenes 5 y 6, Figura 5).Las mayores velocidades encontradas enFIORDOS_NIVEL1 se presentan en la zonadel Canal de Chacao. Concordante con las ob-servaciones realizadas por Salinas & Horma-zabal (2004), las corrientes meridionales en elCanal Moraleda son débiles, produciéndoseademás una clara división del Canal Moraledaal Norte y al Sur de la Constricción de Meni-nea. Finalmente, nuestra principal área de es-tudio (el Fiordo Aysén) recibe una marea mix-ta que es consecuencia del ingreso de la ondadirectamente desde el océano abierto a través,preferentemente, del Canal Darwin así comolos reflejos de la onda a ambos extremos delCanal Moraleda. En este último, las mayoresvelocidades se presentan en la Constricción deMeninea (imágenes 1, 3, y 6, Figura 5).

DISCUSIÓN

La modelación numérica es una herramien-ta importante para el análisis cuantitativo desistemas complejos. Para el caso de sistemasnaturales (eco-sociales) esta herramienta esuna de las pocas disponibles para la genera-ción de escenarios de desarrollo y para el dise-ño de experimentos a escala sistémica. El Mo-delo MOHID ha sido utilizado en numerosos

Estación

CF = Crucero CIMAR FIORDO

Latitud Longitud Tipo de

estación

Periodo

considerado

Ancud (AN) 41º 52' 06" S 73º 49' 53" W SHOA 2003-20042003-20042003-2004

P. Montt (PM) 41º 29' 05" S 72º 57' 39" W SHOAChacabuco (CH) 45º 28' 06" S 72º 49' 09" W SHOAP. Cuptana (PC) 44º 39' 50" S 73º 37' 33" W CF 4 Oct.’98-Feb. ’99

Oct.’98-Feb. ’99Oct.’98-Feb. ’99Oct.’98-Feb. ’99

-I. Isquiliac (II) 45º 24' 35" S 74º 17' 14" W CF 4I. Castillo (CI) 45º 18' 17" S 73º 42' 53" W CF 4C. Tronador (CT) 45º 31' 30" S 73º 34' 30" W CF 4

Tabla I. Nombre, localización geográfica y período considerado de las estaciones mareográficas usadas parala validación de los niveles de marea generados por el modelo FIORDOS_NIVEL1.

Table I. Name, geographic location and time period for the tide-gauge stations used in the validation of thesea levels generated by the FIORDOS_NIVEL 1 model.


estudios integrando una amplia variedad deprocesos y escalas. El modelo MOHID ha sidoaplicado satisfactoriamente en varios estuarios(Sado (Martins et al., 2001); Tajo (Leitão,2003; Braunschweig et al., 2003); Scheldt(Campuzano et al., 2004)) y áreas costeras(Ría de Pontevedra (Villarreal et al., 2002)) asícomo en estudios de océano abierto (MargenOceánico Europeo (Santos et al., 2002), CostaIberica (Coelho et al., 2002); Costa Brasileña(Leitão, 2004), además de las áreas de estudioque se incluyen en el presente proyecto.

En este trabajo hemos mostrado el desarro-llo de un modelo hidrodinámico, barotrópico,para la zona de los fiordos australes de Chile(41° S - 46° S). La comparación de las salidasde FIORDOS_NIVEL1 con la informacióndisponible respecto de la evolución del niveldel mar en la zona muestra que el modelo re-produce eficientemente las variaciones de estavariable en diversos lugares tanto en el inte-rior de los fiordos (e.g. Bahía Chacabuco)

como en zonas intermedias (Puerto Montt,Puerto Cuptana) y externas (Ancud). Los es-tudios sobre corrientes realizados en la zona(e.g. Salinas & Hormazabal, 2004) muestranque la marea es una señal dominante. Estemodelo puede, por tanto, ser usado como baseconceptual y numérica para el estudio y aná-lisis de sub-sistemas (por medio del anida-miento) al interior de sus límites. Si se consi-dera que FIORDOS_NIVEL1 está implemen-tado por medio del sistema numérico MO-HID, entonces se dispone de una amplia gamade herramientas y módulos capaces de serusados para la generación de modelos de ca-lidad de agua, lagrangianos, de ciclos de nu-trientes y para la incorporación de forzantesprovenientes de las cuencas hidrográficas.Sin embargo, aún cuando el modelo reprodu-ce de buena forma la señal de la marea, noreproduce algunas características de los flu-jos de superficie como por ejemplo la unidi-reccionalidad que caracterizan la circulacióndel Canal Moraleda (Fierro et al., 2000). Ello

Fig. 4: Relación entre los valores observados del nivel del mar con relación a los predichos por el modeloFIORDOS_NIVEL1 para las cuatro principales estaciones mareográficas.

Fig. 4: Relationships between observed sea level values and those predicted by the FIORDOS_NIVEL 1model for the four main tide-gauge stations.

Observed Values (m)

Observ

ed

Valu

es

(m)

Observed Values (m) Observed Values (m)

Predicted Values (m)

Pre

dic

ted

Valu

es

(m)

Pre

dic

ted

Valu

es

(m)

4.50 7.00

6.00

5.00

4.00

3.00

2.00

1.00

-1.00

-1.00

0.00

0.000.00

0.000.00

0.50

0.50

1.001.00

1.001.00

1.50

1.50

2.00

2.00

2.002.00

2.50

2.50

3.00

3.00

3.003.00

3.50

3.50

4.00

4.00

4.004.00

4.50

4.50

5.00

5.00

5.00

6.00

6.005.00

1.00 2.00 3.00 4.00 5.00 6.00 7.00

4.50

4.00

4.00

3.50

3.50

3.00

3.00

2.50

2.50

2.00

1.501.50 2.00

Pre

dic

ted

Valu

es

(m)

Ancud r = 0.94252

I. Castillo r = 0.94252

P. Chacabuco r = 0.93642

P. Montt r = 0.96322


Fig. 5: Secuencia de seis imágenes (ordenadas del 1 al 6) del nivel del mar y de la intensidad de las corrientespara un período de 12 horas en condiciones de sicigia.

Fig. 5: Six image sequence (ordered from 1 to 6) of the sea level and current intensity for a 12-hour periodduring spring tide condition.

-76.68

-41.15

-41.57

-42.00

-42.42

-42.84

-43.27

-43.69

-44.12

-44.54

-44.96

-45.39

-45.81

-45.24

-46.66

-41.15

-41.57

-42.00

-42.42

-42.84

-43.27

-43.69

-44.12

-44.54

-44.96

-45.39

-45.81

-45.24

-46.66-76.32 -75.90 -75.47 -74.62-75.05 -74.20 -73.78 -73.35 -72.93 -72.50

-76.26 -75.63 -75.41 -74.98 -74.56 -74.14 -73.71 -73.29 -72.86 -72.44

Wa

ter

Le

ve

l(m

)

Father Model HydrodynamicsForced with FES 2004

1998/10/6

21h 0min

6.50

4.75

3.00

1.25

-5.00

1)-76.68

-41.15

-41.57

-42.00

-42.42

-42.84

-43.27

-43.69

-44.12

-44.54

-44.96

-45.39

-45.81

-45.24

-46.66

-41.15

-41.57

-42.00

-42.42

-42.84

-43.27

-43.69

-44.12

-44.54

-44.96

-45.39

-45.81

-45.24

-46.66-76.32 -75.90 -75.47 -74.62-75.05 -74.20 -73.78 -73.35 -72.93 -72.50

-76.26 -75.63 -75.41 -74.98 -74.56 -74.14 -73.71 -73.29 -72.86 -72.44

Wa

ter

Le

ve

l(m

)


1998/10/6

23h 0min

6.50

4.75

3.00

1.25

-5.00

2)

-76.68

-41.15

-41.57

-42.00

-42.42

-42.84

-43.27

-43.69

-44.12

-44.54

-44.96

-45.39

-45.81

-45.24

-46.66

-41.15

-41.57

-42.00

-42.42

-42.84

-43.27

-43.69

-44.12

-44.54

-44.96

-45.39

-45.81

-45.24

-46.66-76.32 -75.90 -75.47 -74.62-75.05 -74.20 -73.78 -73.35 -72.93 -72.50

-76.26 -75.63 -75.41 -74.98 -74.56 -74.14 -73.71 -73.29 -72.86 -72.44

Wa

ter

Le

ve

l(m

)


1998/10/7

1h 0min

6.50

4.75

3.00

1.25

-5.00

3)-76.68

-41.15

-41.57

-42.00

-42.42

-42.84

-43.27

-43.69

-44.12

-44.54

-44.96

-45.39

-45.81

-45.24

-46.66

-41.15

-41.57

-42.00

-42.42

-42.84

-43.27

-43.69

-44.12

-44.54

-44.96

-45.39

-45.81

-45.24

-46.66-76.32 -75.90 -75.47 -74.62-75.05 -74.20 -73.78 -73.35 -72.93 -72.50

-76.26 -75.63 -75.41 -74.98 -74.56 -74.14 -73.71 -73.29 -72.86 -72.44

Wa

ter

Le

ve

l(m

)


1998/10/7

3h 0min

6.50

4.75

3.00

1.25

-5.00

4)

-76.68

-41.15

-41.57

-42.00

-42.42

-42.84

-43.27

-43.69

-44.12

-44.54

-44.96

-45.39

-45.81

-45.24

-46.66

-41.15

-41.57

-42.00

-42.42

-42.84

-43.27

-43.69

-44.12

-44.54

-44.96

-45.39

-45.81

-45.24

-46.66-76.32 -75.90 -75.47 -74.62-75.05 -74.20 -73.78 -73.35 -72.93 -72.50

-76.26 -75.63 -75.41 -74.98 -74.56 -74.14 -73.71 -73.29 -72.86 -72.44

Wa

ter

Le

ve

l(m

)


1998/10/7

5h 0min

6.50

4.75

3.00

1.25

-5.00

5)-76.68

-41.15

-41.57

-42.00

-42.42

-42.84

-43.27

-43.69

-44.12

-44.54

-44.96

-45.39

-45.81

-45.24

-46.66

-41.15

-41.57

-42.00

-42.42

-42.84

-43.27

-43.69

-44.12

-44.54

-44.96

-45.39

-45.81

-45.24

-46.66-76.32 -75.90 -75.47 -74.62-75.05 -74.20 -73.78 -73.35 -72.93 -72.50

-76.26 -75.63 -75.41 -74.98 -74.56 -74.14 -73.71 -73.29 -72.86 -72.44

Wa

ter

Le

ve

l(m

)


1998/10/7

7h 0min

6.50

4.75

3.00

1.25

-5.00

6)


es debido a que FIORDOS_NIVEL1 no in-corpora componentes baroclínicas, ni el in-greso de agua desde las cuencas hidrográfi-cas. Ellas deben ser incorporadas para re-producir la estructura vertical que caracte-riza la circulación en sistemas estuarinos.De hecho, la implementación de un mode-lo baroclínico para el fiordo Aysén (en de-sarrollo) muestra que es posible reproducirtal estructura vertical por medio de MO-HID.

Por otra parte, es importante enfatizar lanecesidad de contar con una buena base dedatos (de libre disposición) respecto de labatimetría en la zona. Experimentos numé-ricos realizados durante el desarrollo deFIORDOS_NIVEL1 mostraron que debidoa la condición mixta que caracteriza a la se-ñal de la marea en zonas tales como el fior-do Aysén (e.g. sub-sistema Boca del Gua-fo-Canal Moraleda; ingreso de otra señal através de los canales intermedios entre elocéano y los fiordos internos) la altura delnivel del mar en zonas internas es altamen-te sensible a la batimetría de los canales in-termedios. El mejoramiento de la capacidadpredictiva de la señal de mareas en diver-sos sectores al interior del modelo, depen-derá por tanto en gran medida de la incor-poración de nueva información batimétrica.

AGRADECIMIENTOS

Este trabajo forma parte de ECOManage,proyecto financiado por el Sexto ProgramaMarco de la Unión Europea bajo contrato N°INCO-CT-2004-003715, del cual ambos auto-res son co-investigadores y fue parcialmentefinanciado por fondos del Programa Bicente-nario de Ciencia y Tecnología de CONICYT/Banco Mundial y del Concurso de Proyectospara el Crucero CIMAR 9 Fiordos ambos otor-gados a V. Marín. Los autores desean agrade-cer al Servicio Hidrográfico y Oceanográficode la Armada de Chile por hacer disponibleslos datos de batimetría y de series de tiempode nivel del mar para la zona de estudios y alComité Ejecutivo del Programa CIMAR Fior-dos por el uso de la base de datos del cruceroCIMAR 4 Fiordos.

REFERENCIAS

ARREQUIN-SANCHEZ, F. ZETINA-RE-JON, M., MANICKCHAND-HEILEMAN,S. RAMIREZ-RODRIGUEZ, M. & L. VI-DAL. 2004. Simulated response to harves-ting strategies in an exploited ecosystem inthe southwestern Gulf of Mexico. Ecologi-cal Modelling 172:421-432.

BATTEEN, M., HU, C.-P., BACON, J. L., &C.S. NELSON. 1995. A numerical study ofthe effects of wind forcing on the Chile Cu-rrent system. Journal of Oceanography 51:585-614.

BELL, C., J.M. VASSIE & P.L. WOODWOR-TH. 1999. POL/PSMSL TIDAL ANALY-SIS SOFTWARE KIT 2000 (TASK-2000),Tech. Rep., Permanent Service for MeanSea Level, Proudman Oceanographic Lab.,Bidston Observatory, Birkenhead, Mersey-side.

BLUMBERG, A.F. & L.H. KANTHA. 1985.Open boundary condition for circulationmodels. J. of Hydraulic Engineering, ASCE,111: 237-255.

BOWDEN, K. F. 1983. Physical oceanographyof coastal waters. Ellis Horwood Limited,New York, USA.

BRAUNSCHWEIG, F., F. MARTINS, P. C.LEITÃO & R. NEVES. 2003. A methodo-logy to estimate renewal time scales in es-tuaries: the Tagus Estuary case, Ocean Dy-namics, Volume 53(3), 137 - 145.

BRAUNSCHWEIG, F., CHAMBEL, P., FER-NANDES, L., PINA, P. & R. NEVES2004. The object-oriented design of the in-tegrated modelling system MOHID, Com-putational Methods in Water Resources In-ternational Conference, Chapel Hill, NorthCarolina, USA.

CAMPUZANO, F.J., J.H. ALLEN & T. SCO-TT. 2004. The numerical modelling ofecosystem response to nutrients: Applica-tion to the Scheldt estuary and plume. IECS(Institute of Estuarine and Coastal Studies)


report, University of Hull, for the EuropeanCommission, contract nº EVK3-CT-2000-00040 "EUROTROPH" (report 2 of 4).

CANHAM, C., J.J. COLE & W. K. LAUEN-ROTH. 2003. Models in ecosystem scien-ce. Princeton University Press, Princeton,USA.

CASWELL, H.. 2001. Matrix population mo-dels. Construction, analysis, and interpreta-tion. Sinauer Assoc. Inc. Publishers, Sun-derland, Massachusetts, USA.

COELHO, H. S., NEVES, R. J.J., WHITE, M.,LEITÃO, P.C. & A. J. SANTOS. 2002. Amodel for ocean circulation on the Iberiancoast. Journal of Marine Systems 32:153-179.

DYKE, P. 1996. Modelling marine processes.Prentice Hall, New York, USA.

ESCRIBANO, R., ROSALES, S. A. & J.LBLANCO. 2004. Understanding upwellingcirculation off Antofagasta (northern Chi-le): A three-dimensional numerical-mode-ling approach. Continental Shelf Research24:37-53.

FENNEL, W. & T. NEUMANN. 2004. Intro-duction to the modelling of marine ecosys-tems. Elsevier Oceanography Series 72, El-sevier Scientific Pub. Co., New York, USA.

FIERRO, J., M. BRAVO & M. CASTILLO.2000. Caracterización del régimen de ma-reas y corrientes a lo largo del Canal Mora-leda (43° 54' S - 45° 15' S). Cienc. Tecnol.Mar, 23:3-14.

GEODAS, 2006. - http://www.ngdc.noaa.gov/mgg/geodas/geodas.html (Último acceso:Enero 2006).

HEAPS, S. 1987. Three-Dimensional coastalocean models. Coastal and Estuarine Scien-ce 4, American Geophysical Union, Was-hington D.C., USA.

IOC, IHO, AND BODC (2003), CentenaryEdition of the GEBCO Digital Atlas. Publis-

hed on CD-ROM on behalf of the Intergo-vernmental Oceanographic Commissionand the International Hydrographic Organi-zation as part of the General BathymetricChart of the Oceans; British OceanographicData Centre, Liverpool. http://w w w. n g d c . n o a a . g o v / m g g / g e b c o /gebco.html (Último acceso Mayo 2006).

LE PROVOST C., F. LYARD, J.M. MOLI-NES, M.L. GENCO & F. RABILLOUD.1998. A hydrodynamic ocean tide modelimproved by assimilating a satellite altime-ter derived data set, J. Geophys. Res., Vol.103: 5513-5529.

LE PROVOST, C., 2002. `FES2002. A NewVersion of the FES Tidal solution Series',Abstract Volume, Jason{1 Science WorkingTeam Meeting, Biarritz, France.

LEITÃO, P.C., R. NEVES, H. COELHO, F.BRAUNSCHWEIG & J.C. LEITÃO. 2003.3D Hydrodynamic Modelling of the TagusRegion of Fresh Water Influence. Procee-dings of the 4rd Symposium on the IberianAtlantic Margin. Vigo, 2003.

LEITÃO, P.C., J.C. LEITÃO, R. NEVES, G.BERZIN & A.J.R. SILVA. 2004. Hydrody-namics And Transport In The Coastal ZoneOf São Paulo - Brazil. Proceedings 29th In-ternational Conference of Coastal Enginee-ring, Vol. 3, pp 3316-3328, Lisbon 2004.

MARÍN, V. H. 1997. A simple-biology, stage-structured population model of the springdynamics of Calanus chilensis at Mejillonesdel Sur Bay, Chile. Ecological Modelling105:65-82.

MARÍN, V. H., L. E. DELGADO & R. ESCRI-BANO. 2003. Upwelling shadows at Meji-llones Bay (northern Chilean coast): a re-mote sensing in situ analysis. Invest. Mar.,Valparaíso, 31:47-55.

MARTINS, F., P. C. LEITÃO, A. SILVA & R.NEVES. 2001. 3D modelling of the SadoEstuary using a new generic vertical discre-tization approach, Oceanologica Acta, 24(1), 51-62.


MEADOWS, D., J. RANDERS & D. MEA-DOWS. 2004. Limits to growth. The 30-Year update. Chelsea Green Pub. Co., Whi-te River Junction Vermont, USA.

NIHOUL, J. C. J. 1975. Modelling of marinesystems. Elsevier Oceanography Series, 10.Elsevier Scientific Pub. Co., New York,USA.

PALMA, W., ESCRIBANO, R., & S.A. RO-SALES. 2006. Modeling study of seasonaland inter-annual variability of circulation inthe coastal upwelling site of the El Loa Ri-ver off northern Chile. Estuarine, Coastaland Shelf Science 67:93-107.

PENG, G., LESLIE, L. M. & Y. SHAO. 2002.Environmental modelling and prediction.Springer-Verlag, Berlín, Alemania.

SALINAS, S. & S. HORMAZÁBAL. 2004.Capacidad de transporte de la constricción

de Meninea para un flujo de dos capas yel efecto de la corriente de marea. Cienc.Tecnol. Mar., Chile 27:5-15.

SALOMONS, W., TURNER, R. K., DE LA-CERDA, L. D. & S. RAMACHAN-DRAN. 1999. Perspectives on integratedcoastal zone management. Springer-Ver-lag, Berlín.

SANTOS, A., H. MARTINS, H. COELHO,R.NEVES, P. C. LEITÃO. 2002. A circu-lation model for the European oceanmargin, Applied Mathematical Mode-lling, 26(5), 563-582.

VILLARREAL, M. R., P. MONTERO, J. J.TABUADA, R. PREGO, P. C. LEITÃO& V. PÉREZ-VILLAR. 2002. Hydrody-namic model study of the Ria de Ponte-vedra under estuarine conditions, Estua-rine, Coastal and Shelf Science, 54, 101-113.

un modelo hidrodinÁmico-barotrÓpico para los … marin.pdf · son más recientes (salomons et...

Documents