3.1. Marco geodésico y gravimétrico

El marco de referencia geodésico de nuestra área de estudio está compuesto por las diferentes redes geodésicas establecidas en la Antártida, entre las que se encuentran: Red Geodésica de la Antártida Española (RGAE), Red Geodésica de Isla Decepción (REGID), Red de Nivelación de Isla Decepción (RENID) y la Red Gravimétrica de Isla Decepción (REGRID), esta última constituye el marco gravimétrico de la isla.

Todas las estaciones de estas redes están dotadas de coordenadas geocéntricas tridimensionales (X, Y, Z), transformadas a coordenadas geodésicas latitud (φ), longitud (λ) y altura elipsoidal (h) referidas al Marco Internacional de Referencia Terrestre ITRF2000.0 y ajustadas con un cierre de precisión milimétrico (^{Altamimi, Sillard y Boucher, 2002}), como veremos más adelante. Para los trabajos cartográficos se han utilizado las coordenadas s planas rectangulares UTM (Universal Transversa de Mercator) para la Zona UTM 20 sur, Meridiano Central 63º de longitud oeste.

3.1.1. Red Geodésica de la Antártida Española (RGAE)

La red RGAE fue la primera red diseñada en la Antártida por investigadores españoles, para definir un marco de referencia geodésico en la región de las Islas Shetland del Sur, el Mar de Bransfield y la Península Antártica (Figura 2).

Figura 2 Red RGAE y la ubicación de sus vértices geodésicos. 

Utilizando algunos vértices de esta red como puntos fundamentales, se determinaron las coordenadas de los vértices de la red REGID. Para este trabajo se utilizaron los siguientes: BEJC (BAE “Juan Carlos I”, isla Livingston), BEGC (BAE Gabriel de Castilla, isla Decepción) PEND (Caleta Péndulo, isla Decepción). Mayores detalles de esta red geodésica (^{Berrocoso et al., 2006}a, 2006b, 2006c, 2008).

3.1.2. Red Geodésica de Isla Decepción (REGID)

La Red REGID es la segunda red geodésica establecida en la Antártida, específicamente en isla Decepción, con la finalidad de constituir el marco de referencia geodésico para cualquier actividad científica o técnica y establecer los puntos de control para determinar los modelos de deformación horizontal debidos a la actividad volcano-tectónica. La red está constituida por 13 vértices geodésicos distribuidos alrededor de Puerto Foster, que se pueden ver y en la Figura 3 y coordenadas ajustadas de la red se muestran en la Tabla 1. El benchmark geodésico LN000 ha sido incluido en la red REGID, pues constituye el Punto Fundamental de Referencia Vertical de la red. Más detalles de la misma los podemos ver en Berrocoso et al. (2006a, 2006b, 2007, 2008).

Figura 3 Red de Nivelación de Islas Decepción (RENID). 

Tabla 1 Vértices red REGID. Coordenadas y precisiones (ITRF 2000.0) 

No	VERTICE	Latitud S (° ‘ “)			Longitud W (° ‘ “)			h (m)	σ φ (m)	σ λ (m)	σ h(m)
1	BEGC	62	58	43.6576	60	40	27.5304	82.070	0.001	0.001	0.002
2	BALL	62	58	38.5553	60	33	52.5089	25.968	0.003	0.001	0.006
3	FUMA	62	57	41.0170	60	42	59.3471	22.911	0.001	0.001	0.006
4	PEND	62	56	09.8456	60	35	34.3437	28.841	0.001	0.001	0.003
5	COLA	62	59	27.9180	60	37	31.8014	48.050	0.003	0.001	0.009
6	GLAN	62	57	58.3621	60	35	23.8528	27.541	0.003	0.003	0.012
7	GEOD	62	58	56.4127	60	39	11.7290	42.182	0.003	0.001	0.009
8	UCA1	62	56	28.4103	60	41	28.0891	28.667	0.003	0.001	0.009
9	CR70	62	55	23.6706	60	38	01.0013	23.600	0.003	0.003	0.015
10	TELE	62	55	37.9905	60	41	25.5485	23.791	0.003	0.003	0.012
11	BOMB	62	55	08.4201	60	39	33.8481	23.785	0.003	0.003	0.015
12	BARG	62	58	30.2679	60	41	53.3601	22.283	0.001	0.001	0.006

3.1.3. Red de Nivelación de Isla Decepción (RENID).

La red de nivelación RENID consta de 60 marcas de nivelación (benchmarks), divididas en seis líneas de nivelación independientes.

El punto fundamental de referencia de altitud es el LN00, situado junto a la Base Antártica Española "Gabriel de Castilla" y las líneas de nivelación que están vinculada al menos a uno de los vértices geodésicos de la red REGID, fueron establecidas en diferentes trabajos (^{Berrocoso et al. 2006}a, 2006b, 2007, 2008; ^{Vidal et al. 2012}). La Figura 4 muestra las diferentes líneas de nivelación que componen la red RENID y los vértices geodésicos de la red REGID a los que están vinculados. En la Tabla 2 se muestran los datos de nivelación obtenidos en todos los puntos pertenecientes a la red RENID, así como los errores asociados a cada uno de ellos.

Figura 4 a. Ubicación de las estaciones mareográficas, b. Isla Livingston, ubicación de LIVMAR y sus correspondientes TGBM, c. Isla Decepción, ubicación de DECMAR y sus correspondientes TGBM. 

Tabla 2 Puntos de la red RENID y con sus diferencias de nivel y sus errores asociados 

No	ESTACION	Δ n 2014 (m)	σ Δ n 2014 (m)		No	ESTACION	Δ n 2014 (m)	σ Δ n 2014 (m)
1	LN00	0.000	0.008		34	LN119	-1.94	0.020
2	BEGCC	57.025	0.003		35	LN120	-0.38	0.020
3	BALL	0.723	0.026		36	LN121	2.949	0.018
4	FUMA	-2.771	0.007		37	LN122	5.776	0.018
5	PEND	3.562	0.038		38	LN123	12.02	0.013
6	COLA	22.847	0.010		39	LN124	0.051	0.008
7	GLAN	2.213	0.033		40	LN125	-4.123	0.002
8	GEOD	17.018	0.014		41	LN201	17.845	0.009
9	UCA1	3.513	0.032		42	LN202	6.567	0.011
10	CR70	-1.531	0.042		43	LN203	12.550	0.022
11	TELE	-1.394	0.068		44	LN301	5.757	0.046
12	BOMB	-1.353	0.20		45	LN302	19.817	0.050
13	GRAV	23.456	0.000		46	LN303	-0.109	0.034
14	BARG	-2.883	0.008		47	LN401	4.645	0.003
15	TOCO	-3.640	0.008		48	LN402	12.409	-0.013
16	LN101	2.136	0.003		49	LN403	20.640	-0.020
17	LN102	8.606	0.004		50	LN404	32.454	-0.012
18	LN103	-0.756	0.008		51	LN501	-2.270	0.074
19	LN104	-3.845	0.014		52	LN502	19.396	0.056
20	LN105	4.960	0.009		53	LN503	26.238	0.061
21	LN106	1.167	0.006		54	BR-01	17.059	0.031
22	LN107	2.843	0.007		55	LN601	-4.347	0.054
23	LN108	0.266	0.000		56	LN602	1.785	0.050
24	LN109	3.758	0.004		57	LN603	8.167	0.043
25	LN110	20.362	0.003		58	LN604	-3.911	0.046
26	LN111	24.279	0.001		59	LN605	3.668	0.048
27	LN112	10.539	0.009		60	LN606	7.238	0.053
28	LN113	-1.488	0.004		61	LN607	12.875	0.074
29	LN114	-4.002	0.008		62	LN608	1.995	0.070
30	LN115	-4.655	0.011
31	LN116	-3.804	0.011		201	BEJC	0.000	0.000
32	LN117	-3.74	0.012		202	BEJ1	1.025	0.002
33	LN118	-2.082	0.012		203	TOJO	-9.482	0.000

3.2. Marco oceanográfico

El marco oceanográfico en el área de estudio está definido por las estaciones mareograficas, puntos de control de mareógrafos (TGBM, en inglés: Tide Gauge benchmark), puntos auxiliares de referencia de mareógrafos (TGAR, en inglés: Tide Gauge Auxiliary Reference) y los vértices de referencia. Se instalaron dos estaciones mareográficas en las islas Livingston y Decepción denominadas LIVMAR y DECMAR, que se muestran en la Tabla 5.

La Figura 4 muestra la ubicación de las estaciones mareográficas, con sus respectivos Puntos de Control de Mareógrafos (TGBM) y sus referencias auxiliares (TGAR).

3.3. Equipos utilizados para toma de datos

Los datos de marea se obtuvieron mediante dos sensores de presión y temperatura SAIV modelo TD304, CTD 204 SAIV SD y AQUAlogger PT520. Los sensores de presión fueron fondeados a unos 30 metros de la línea de playa y a profundidades aproximadas de 7.5 y 3.5 metros en LIVMAR y DECMAR respectivamente. Para el fondeo se utilizó una embarcación neumática Zodiac Pro 500.

Las estaciones geodésicas tienen coordenadas absolutas con precisión milimétrica en X, Y, Z, y en latitud (φ), longitud (λ) y altura elipsoidal (h) relativas al Sistema de Referencia ITRF2000.00 (^{Altamimi et al., 2002}). Para el posicionamento GPS se utilizaron receptores geodésicos Leica GX1230 y TRIMBLE 5700, en posicionamento estático (fase) con postprocesamiento y antena estándar, como se especifica en ^{Seeber (2003)}.

La nivelación geométrica se realizó utilizando un nivel óptico Leica modelo NA2 con una precisión de ± 0.7 mm para un kilómetro de doble recorrido. Para las nivelaciones trigonométricas y los enlaces geodésicos, utilizamos la Estación Total Leica modelo TPS 403. A estas medidas realizadas se le aplicaron las correcciones por esfericidad terrestre y efectos de refracción atmosférica.

Las observaciones gravimétricas se realizaron con un gravímetro relativo LaCoste & Romberg modelo D-203. Se aplicaron correcciones de mareas, altura y deriva a todo el conjunto de datos gravimétricos. Los detalles de los sensores y equipos utilizados los podemos ver en Jigena et al. (²⁰¹⁵, 2016).

4.1. Determinación de las constituyentes de marea

El período de registro de marea en ambas estaciones fue del 3 de febrero de 2011 al 11 de abril de 2013 (Figura 5). Todas las medidas de las series se normalizaron a 1 dato/hora. Las series temporales en cada estación fueron sometidas a un análisis armónico utilizando una aplicación desarrollada de acuerdo a ^{Foreman (1977)}. Durante todo el período del registro de datos de marea, se obtuvieron datos meteorológicos de presión atmosférica y temperatura del aire, tomados por la Agencia Española de Meteorología (AEMET), en las estaciones de “Gabriel de Castilla” (isla Decepción) y “Juan Carlos I” (isla Livingston).

Figura 5 Series de mareas obtenidas en LIVMAR y DECMAR. 

La presión atmosférica en la estación de isla Livingston se tomó como serie de referencia, al ser más estable y completa. La presión atmosférica enLivingston osciló entre un mínimo de 943,35 mb y un máximo de 1 025,59, con un promedio de 988,27 mb utilizado como referencia. Para convertir la presión hidrostática a una altura equivalente del nivel del mar, usamos (Ecuación 1):

donde, P es la presión registrada por el sensor del mareógrafo, Pa es la presión atmosférica de referencia, utilizando el valor constante de 988,27 mb, g es la aceleración debida a la gravedad cuyo valor medio en Isla Decepción es 9,822956 m/s2 y ρ la densidad del agua en el área de estudio, cuyo valor medio es de 1025 kg/m3, siguiendo las especificaciones de la Organización de las Naciones Unidas para la Educación, la Ciencia y la Cultura (UNESCO, ¹⁹⁸¹, ¹⁹⁹⁴) y Fofonot & Millard (¹⁹⁸³). Se realizó la corrección por efecto de barómetro invertido de acuerdo con Muñoz & Abarca (²⁰⁰⁹).

Las series de marea tienen un total de 19 143 datos horarios, que supone un registro de más de dos años de observación (798 días) en ambas estaciones. En el análisis armónico obtuvimos una condición de matriz de 0.82 y una bondad de ajuste de 0.993, teniendo como resultado 70 constituyentes de marea, siendo diecinueve (19) las más importantes, con amplitud mayor a un centímetro y un aporte de energía del 93% del total de la onda. En todos los análisis armónicos realizados se obtuvo que la mayor cantidad de energía de la onda de marea (85%) es aportada por ocho componentes, cuatro semidiurnas (M2, S2, K2, N2) y cuatro diurnas (K1, O1, P1, Q1), con un aporte similar para cada grupo de aproximadamente un 42,5%. En lo referente a las amplitudes de las series, los resultados muestran que las amplitudes máximas y mínima son de 4 986 m. y 2 224 m con un rango de 2 762 m para DECMAR, y de 8 334 m y 6 169 m con un rango de 2 165 m para LIVMAR. El factor de forma de marea o índice de Courtier (C) obtenido es 0,91 y 0.90 para LIVMAR y DECMAR respectivamente, determina que las mareas tienen un régimen mixto, con una componente predominantemente semidiurnas (0.25 < C < 1.50), definidas por ^{Defant (1961)}, ^{Dragani et al. (2004)} y ^{Vidal et al. (2012)}.

En general, los resultados en amplitud, desfase y régimen de marea, no difieren significativamente de los obtenidos por (^{Smithson, 1992}; SCAR, 1993; ^{López et al., 1993}; ^{García, 1994}; ^{Schöne et al., 1998}; ^{López et al., 1999}; ^{Speroni et al., 2000}; ^{Padman et al., 2002}; ^{D’onofrio et al., 2003}; ^{Dragani et al., 2004}; ^{Vidal et al., 2012}), sin embargo, se observan variaciones significativas en la constituyente K1 con respecto a los valores dados por (^{Dragani et al., 2004}, ^{Vidal et al., 2012}). El resultado del ajuste entre los datos reales y los datos predichos es de un 99.3%, que es un parámetro excelente y muy representativo para realizar predicciones en base a los datos obtenidos y utilizados en el análisis. Los resultados obtenidos en amplitud y desfase de las principales componentes de marea en ambas estaciones se muestran en la Tabla 6 y mayores detalles se pueden consultar en ^{Jigena et al. (2015)}.

4.2. Determinación del Nivel Medio del Mar y georreferenciación vertical

Para el control vertical, los datos de marea fueron referidos altimétricamente a un TGBM, que a su vez está enlazado con otras marcas auxiliares de referencia de nivel TGAR y con el vértice geodésico más cercano, STGPS (en inglés: Station GPS). Toda la referenciación altimétrica se realizó mediante nivelación geométrica y el esquema de la misma lo podemos ver en la Figura 6. Las mediciones de control, entre el PTG (en inglés: Principal Tide Gauge) y el TSG (en inglés: Tide Staff Gauge) se realizaron en días escogidos, sin viento y sin olas (mar en calma) tomando las lecturas simultáneas al menos 1 hora, con una frecuencia de 10 minutos.

Fuente: adaptada de ^{Jigena et al. (2015})

Figura 6 Vinculación altimétrica entre los TGBM (marca de referencia de mareógrafos), TGAR (marca auxiliar de referencia), PTG (sensor principal) y TSG (sensor móvil). 

Se realizó un ajuste lineal entre la lectura del sensor principal PTG, transformada a profundidad y corregida por presión atmosférica (D), contra los datos de altura (H) tomados por el sensor de control TSG, instalado en la regla de marea y corregidos por la referencia al TGBM mediante nivelación geométrica. Este ajuste nos proporciona un valor medio del nivel de marea para cada instante de la lectura referido al TGBM y los valores medios de profundidad de fondeo de los sensores y cuyos resultados se aprecian en la Figura 7. Los detalles y la metodología se encuentran explicados en ^{Vidal et al. (2012)} y Jigena et al. (²⁰¹⁴, 2015).

Figura 7 Ajuste lineal entre las medidas instantáneas del nivel medio del mar. 

En la Tabla 7 vemos la relación entre las estaciones mareográficas con respecto a los vértices geodésicos de referencia, con los resultados obtenidos en los vértices de BEJC, COLA y LN00, con los errores estimados asociados alos mismos. El valor obtenido del Nivel Medio del Mar (NMM) se corresponde con el valor de la altura ortométrica (H) para el TGBM de referencia con su error correspondiente en la determinación del NMM para cada estación, para DECMAR (CINF) y para LIVMAR (CIJO).

Se obtuvo una altura ortométrica muy precisa en todos los puntos de control de mareógrafos (TGBM) y en los vértices geodésicos de referencia (BEJC, COLA y LN00), corregida por efectos gravimétricos. La nivelación geométrica entre los TGBM, TGAR y los vértices geodésicos de referencia (STGPS) cumple con los estándares para nivelación geométrica de Primer Orden Geodésico. El benchmark LN00, fue definido como Punto Fundamental de Referencia Altimétrica de la red REGID y RENID, y su cota ortométrica queda determinada en 6 200 msnm. En la Tabla 7 se muestran los resultados finales del cálculo de los valores del NMM, donde el valor calculado del error de H (σH) incluye el error por medida y lectura del NMM. Además en la tabla h = altura alipsoidal; Δn = Diferencia de nivel (nivelación geométrica); N = Ondulación del geoide; H = Altura Ortométrica.

Las nuevas cotas ortométricas obtenidas tienen una diferencia + 0,06 m para BEJC y de + 0.20m para COLA, respecto a los valores obtenidos por (^{Vidal et al., 2012}). En general, los resultados obtenidos en este trabajo han mejorado los trabajos anteriores realizados en la zona (SHN, 1970; ^{Dragani et al., 2004}; Vidal et al., 2012).

4.3. Determinación del Geoide GEIODEC14

Otro de los resultados de este trabajo, fue la determinación de un geoide experimental de precisión para la isla Decepción, denominado GeoiDEC14, además de tres puntos dotados con ondulación del geoide en isla Livingston, en el área aledaña a la BAE “Juan Carlos I” y a Caleta Johnsons. Para la determinación del geoide de precisión se utilizó la metodología GPS/Nivelación/Gravedad, debido a que la misma computacionalmente es fácil de implementar, ya que simplemente se manejan vectores del tamaño de la cantidad de datos disponibles y el método se basa en la resolución de ecuaciones de segundo grado.

El método desarrollado está basado en la determinación de la ondulación del geoide (N) a partir de la diferencia entre la altura elipsoidal (h), obtenida directamente mediante observaciones GPS, y la altura ortométrica (H) obtenida por nivelación geodésica, geométrica o trigonométrica y corregida por efectos gravimétricos mediante medidas de la gravedad (^{Heiskanen & Moritz, 1985}; ^{Torge, 2001}; ^{Berrocoso et al., 1996}), para ello utilizamos la fórmula que aproximadamente es igual a (Ecuación 2):

donde h es la altura elipsoidal, H la altura ortométrica y N es la ondulación del geoide. Sabiendo que en un punto P de la superficie terrestre se conocen la cota ortométrica (H) y la altura elipsoidal (h), se puede determinar de forma muy aproximada la ondulación del geoide (N). La altura elipsoidal (h) la obtenemos directamente por observaciones GPS y partimos de un punto de altura ortométrica (H) conocida, el problema para la aplicación del método se reduce a obtener la altura ortométrica en los demás puntos. Para calcular la altura ortométrica en un punto Pj (HPj), obtenemos el siguiente desarrollo que nos permita calcular el de H partiendo de la diferencia de alturas ortométricas ΔHPi;j entre el punto de cálculo Pj y el punto de referencia Pi, y que nos da como resultado la siguiente Ecuación 3:

Es una ecuación de segundo grado, donde la incógnita es la altura ortométrica Hpj ) del punto Pj que puede ser calculada a partir de la altura ortométrica Hpi ) en el punto de referencia Pi, la diferencia de nivel ∆npi,j entre los puntos Pi y Pj y los valores de la gravedad absoluta en estos puntos, gpi y gpj. (Tabla 8)

Con estos datos se obtuvo el mapa de curvas de N aplicando el método interpolación para toda la isla (Figura 8 a). Igualmente, se obtuvo el mapa de errores en metros de la ondulación del geoide (Figura 8 b).

Figura 8 a. Valores de N obtenidos para GeoiDEC14, b. Errores estimados de N agrupados por zonas. 

Observamos que los mayores errores se obtuvieron en las zonas Cráteres del 70, Caleta Péndulo y Punta Collins, mientras que en la zona de Bahía Balleneros se observó un cambio de tendencia. En general, los errores estimados en la determinación de GeoiDEC14 no exceden los 20 cm y pueden agruparse en tres tipos y están relacionados geográficamente a un área concreta de la isla como se muestra en la Tabla 9.

En la isla Livingston, área de la BAE “Juan Carlos I” y caleta Johnsons, se han establecido tres puntos provistos de ondulación del geoide (N), altura ortométrica y coordenadas absolutas GNSS. En la Tabla 10 se muestran los resultados de estas observaciones en la isla Livingston y la ubicación de las estaciones se muestran en la Figura 4.

La importancia de este trabajo radica en que es el primer modelo de un geoide local de precisión obtenido en una zona volcánica activa, como es el volcán Decepción y en esta zona de la Antártida. Este modelo permitirá realizar trabajos técnicos y científicos y producir resultados con un alto nivel de precisión. Como es un modelo preciso puede técnicamente tener la capacidad de ser utilizado para:

Para obtener el modelo GeiDEC14, hemos utilizado la metodología GPS/nivelación/gravedad realizando el cálculo sobre un total de 108 puntos disponibles en toda la isla, distribuidos homogéneamente en la parte interior de la isla alrededor de Puerto Foster. Estos puntos tienen un valor promedio de ondulación del geoide (N) de 18,83 metros con un rango de valores entre 18,05 y 19,49 metros.

Se obtuvieron tres puntos en isla Livingston, todos se encuentran ubicados en la zona entre la BAE “Juan Carlos I” y caleta Johnson. Estos tres puntos tienen un valor promedio de 18,87 metros, muy similar y coherentes con los obtenidos en isla Decepción. Debido a que estos puntos son muy cercanos su obtención no es relevante ni suficiente para la determinación de un geoide local en esta área; sin embargo, sirven para contrastar los datos obtenidos de GEO-DEC14 en Decepción.

Se ha realizado una comparación entre las ondulaciones del geoide obtenidas con modelos globales de geoide como el EGM96, EGM2008, AUIB-Grace2003s y ITSG-Grace2014s, además del geoide experimental de isla Decepción 2007 determinado por ^{Berrocoso et al. (2007}, 2008), contra las obtenidas con el GEOIDEC14, que nos han servido para comparar y analizar sus diferencias con respecto a nuestro geoide experimental. Los resultados se observan en la Figura 9 c y Figura 9 f, los valores de GeoiDEC14 se han introducido en cuadricula del EGM08 y GRACE2014s para mostrar su influencia en la generación de los modelos combinados. Podemos observar que GeoiDEC14 tiene valores más pequeños de N, llegando a los 18,05 metros en algunas zonas, valor no alcanzado por ninguno de los geoides globales estudiados.

Figura 9 Valores de N obtenidos para a) EGM96, b) EGM08, d) AIUB_GRACE03s, e) ITSG-GRACE2014s, c) and f) muestran los resultados de la inserción de GeoiDEC14 dentro de EGM08 y GRACE2014. 

La comparación muestra una diferencia máxima del valor de N, en los puntos observados de alrededor de 2,5 m con respecto a los antiguos modelos globales EGM96, GRACE03 y EGM2008 (Figuras 9a, 9b y 9d). Respecto a los modelos DEC2007 y GRACE2014 (Figura 9e) la diferencia es de 1,9 m.

En la Tabla 11 se presenta un análisis estadístico realizado, denominado Tamaño de Efecto (ES, en inglés Effect Size), que nos sirve para comparar diferentes modelos tomando uno de ellos como referencia, en este caso GEOIDEC14 (^{Ledesma et al., 2008}). Los modelos donde ES es mínimo serían aquellos que se acerquen más a nuestro modelo patrón. Los resultados indican que los valores mínimos corresponden a los modelos GRACE2014 y DEC2007. Con este análisis nuevamente se marca una separación entre los modelos y se confirma que estos son los que tienen la mayor precisión y resolución entre los comparados. Este análisis validaría nuestra hipótesis de que el modelo de geoide local de precisión GeoiDEC14 es el modelo más preciso obtenido hasta la fecha en el área de la isla Decepción.

Los resultados comparativos con otros modelos de geoide indican que los valores mínimos corresponden a los modelos GRACE2014 y DEC2007. Este análisis marca una separación entre los modelos y confirma que estos son los que tienen la mayor precisión y resolución entre los comparados. Este análisis validaría nuestra hipótesis de que el modelo GeoiDEC14 es el más preciso obtenida hasta la fecha en el área de la isla Decepción. La razón de la existencia de una diferencia más pequeña respecto al modelo experimental DEC2007 y al modelo GRACE2014, podemos decir que se deben a que el primero fue obtenido con datos experimentales in situ, mientras que el modelo GRACE2014 es el modelo global más preciso obtenido hasta la fecha, resolución grado 200, que se traduce en una mayor precisión y por lo tanto sus datos se acercan más a modelos obtenidos con datos reales tomados in situ. Más detalles sobre esta comparativa de presentan en ^{Jigena et al. (2016)}.

El estudio de la distribución espacial de los valores de N en GeoiDEC14, muestra que algunas de sus anomalías coinciden con los puntos calientes de la isla, i.e., los valores mínimos se ubican en bahía Fumarolas y bahía Balleneros, que son zonas con presencia de fumarolas y puntos calientes. Por otro lado, los valores máximos los podemos encontrar en zonas donde existen restos de flujo de lavas, i.e., Colatinas, Glaciar Negro o Punta Murature. Estas anomalías pueden deberse a las propiedades de la corteza o a ladensidad del terreno en la zona (^{Crescentini y Amoruso, 2007}), que se pueden apreciar en la geología de la isla.