LAKE2K (o L2K) es un modelo que está diseñado para calcular las tendencias estacionales de la

calidad del agua en lagos estratificados. L2K se implementa dentro del entorno de Microsoft
Windows. Está programado en el lenguaje de macro de Windows: Visual Basic para Aplicaciones
(VBA). Excel se utiliza como la interfaz gráfica de usuario.

Esta documentación y el manual del usuario se han escrito para la versión beta del modelo. El
modelo ya ha sido probado alfa; es decir, ha sido depurado y sometido a pruebas numéricas para
confirmar que arroja resultados numéricamente precisos. Una versión beta representa un
software que está destinado a ser ejecutado por usuarios seleccionados antes de la distribución
general.

EMPEZANDO
 
Se requiere instalación para muchos modelos de calidad del agua. Este no es el caso
de LAKE2K porque el modelo está empaquetado como un libro de trabajo de Excel. El
programa está escrito en el lenguaje de macros de Excel: Visual Basic para
Aplicaciones o VBA. Las hojas de cálculo y los cuadros de trabajo del Libro de Excel se
usan para ingresar datos y mostrar resultados. En consecuencia, solo tiene que abrir el
Libro de trabajo para comenzar a modelar. Las siguientes son algunas instrucciones
paso a paso recomendadas sobre cómo obtener su primera ejecución del modelo.
 
Paso 1: cree una carpeta llamada LAKE2K para contener el libro de trabajo y sus
archivos de datos. Por ejemplo, en el siguiente ejemplo, se crea una carpeta llamada
LAKE2K en la unidad C: \. Crea un subdirectorio fuera de LAKE2K llamado DataFiles.
 
Paso 2: Copie el archivo L2KMaster (L2KMaster0102.xls) de su CD a C: \ LAKE2K.
 
Paso 3: abra Excel y asegúrese de que su nivel de seguridad de macros esté
configurado en medio. Esto se puede hacer utilizando los comandos del menú:
Herramientas  Macro  Seguridad. Asegúrese de que el botón de radio Medio esté
seleccionado como se muestra en la Figura 1.
Figura 1 El cuadro de diálogo de nivel de seguridad de macros de Excel. Para ejecutar
L2K, se debe seleccionar el nivel medio de seguridad.
 
Paso 4: abra L2KMaster0102.xls. Cuando haga esto, se mostrará el cuadro de diálogo
Macro Security Dialogue (Figura 2).

Figura 2 El cuadro de diálogo de seguridad Macro de Excel. Para ejecutar L2K, se debe
seleccionar el botón Habilitar macros.
 
Haga clic en el botón Habilitar macros.
 
Paso 5: guarde inmediatamente el archivo como L2K0102.xls. Este será el Libro de
Excel que usará de forma rutinaria. Si por alguna razón modificas L2K0102.xls de
manera que no se pueda usar, siempre puedes volver a L2KMaster0102.xls como copia
de seguridad.
 
Paso 6: en la hoja de cálculo de LAKE2K, vaya a la celda B10 e ingrese la ruta al
directorio de Archivos de datos: C: \ LAKE2K \ DataFiles como se muestra en la Figura
3.
Figura 3 La hoja de trabajo LAKE2K que muestra la entrada de la ruta del archivo en la
celda B10.

Paso 7: haz clic en el botón Ejecutar. LAKE2K comenzará a ejecutarse. Puede seguir el
progreso de la ejecución en la barra de estado que se muestra en la esquina inferior
izquierda de la hoja de cálculo (Figura 4).

Figura 4 La barra de estado de LAKE2K se encuentra en la esquina inferior izquierda de

la hoja de trabajo. Le permite seguir el progreso de una ejecución de modelo.

Si el programa se ejecuta correctamente, se mostrará el gráfico de temperatura. Si no

funciona correctamente, existen dos posibilidades:

Primero, puede estar usando una versión anterior de Microsoft Office. Aunque Excel es
compatible con versiones anteriores para algunas versiones anteriores, L2K no
funcionará con versiones muy antiguas.

En segundo lugar, es posible que haya cometido un error al implementar los pasos
anteriores. Un error común es haber escrito mal la ruta del archivo que ingresó en la
celda B10. Si este es el caso, recibirá un mensaje de error (Figura 5).

Figura 5 Un mensaje de error que aparecerá si escribe la ruta de archivo incorrecta en

la celda B10 en la hoja de trabajo LAKE2K.

Si esto ocurre, haga clic en el botón de finalización. Esto terminará la ejecución y lo

llevará nuevamente al Libro de Excel. Luego debe volver a la hoja de trabajo LAKE2K y
corregir la entrada de la ruta del archivo. Tenga en cuenta que se producirá el mismo
error si no configuró sus directorios con los nombres correctos como se especificó
anteriormente.
Paso 8: en la hoja de trabajo LAKE2K, haga clic en el botón Abrir archivo antiguo.
Busque el directorio: C: \ LAKE2K \ DataFiles. Debería ver que se ha creado un nuevo
archivo con el nombre especificado en la celda B9 y la extensión .l2k. Haga clic en el
botón Cancelar para regresar a L2K.

Tenga en cuenta que cada vez que se ejecuta L2K, se creará un archivo de datos con el
nombre de archivo especificado en la celda B9 en la hoja de trabajo LAKE2K (Figura 3).
El programa fija automáticamente la extensión
.l2k al nombre del archivo y almacena el archivo en el directorio especificado en la
celda B10. Como esto sobrescribirá el archivo, asegúrese de cambiar el nombre del
archivo cuando realice una nueva aplicación.

3 MANUAL DE USUARIO

3.1 VISIÓN GENERAL

El código de computadora utilizado para implementar los cálculos para LAKE2K está
escrito en Visual Basic para Aplicaciones (VBA). Excel sirve como la interfaz de usuario.
Se usa una serie de hojas de trabajo para ingresar los parámetros del modelo y los
datos del lago. Los resultados y datos del modelo se muestran en hojas de trabajo y
como una serie de gráficos de Excel.

El color se usa para indicar si la información debe ser ingresada por el usuario o
producida por el programa:

• Pale Blue designa los valores de variables y parámetros que debe ingresar el usuario.
• Pale Yellow designa los datos que ingresa el usuario. Los datos se muestran en
gráficos generados por L2K. Por lo tanto, tales datos son opcionales.
• Pale Green designa los valores de salida generados por L2K.
• Los colores sólidos oscuros se usan para las etiquetas y no se deben cambiar.

Todas las hojas de trabajo incluyen dos botones (Figura 6):

• Abrir archivo antiguo. Cuando se hace clic en este botón, el explorador de archivos se
abrirá automáticamente para permitirle acceder a un archivo de datos que
previamente generó L2K. Todos los archivos de datos LAKE2K tienen la extensión, *
.l2k.
• Correr. Este botón hace que L2K se ejecute y cree un archivo de datos que
contenga los valores de entrada. Luego, se puede acceder al archivo de datos usando

el botón Abrir archivo antiguo.

Figura 6. Los botones utilizados en L2K.

Las siguientes secciones describen las hojas de trabajo que componen la interfaz L2K.
  3.2 HOJAS DE TRABAJO DE PARÁMETROS MODELO Se usa una serie de hojas de trabajo
para ingresar la información requerida para ejecutar el modelo. En las versiones actuales de
Microsoft Office, estas hojas de trabajo se distinguen por pestañas turquesas.

3.2.1 Hoja de trabajo LAKE2K La hoja de trabajo LAKE2K (Figura 7) se usa para ingresar
información general sobre una aplicación de modelo en particular

Figura 7. La hoja de trabajo LAKE2K.

Nombre del lago Nombre del lago que se modela. Después de ejecutar el programa,
este nombre junto con la fecha se muestra en todas las hojas y cuadros.
Nombre del archivo. Este es el nombre del archivo de datos generado cuando se
ejecuta L2K.
Directorio donde se guardó el archivo. Esto especifica la ruta completa al directorio
donde se guarda el archivo de datos.
Latitud. El usuario debe ingresar la latitud del lago en grados decimales.
Elevación. El usuario debe ingresar la elevación de la superficie del lago sobre el nivel
del mar en metros. Paso de cálculo Este es el paso de tiempo utilizado para el cálculo.
Debe ser menor o igual a un día. Si el usuario ingresa un valor superior a 1 día, el
programa establece automáticamente el paso de cálculo en 1 día.
Paso de impresión Este es el paso del tiempo para el cual se muestran los resultados.
Debe ser mayor o igual que el paso de cálculo. Si el usuario ingresa un valor inferior al
paso de cálculo, el programa establece automáticamente el paso de impresión en 1
día.
Hora inicial. Este es el tiempo de inicio del cálculo ingresado como datos usando el
formato MM / DD / YR.
Tiempo final. Este es el tiempo de finalización del cálculo ingresado como datos usando
el formato MM / DD / YR.
Hora del último cálculo (salida). La computadora muestra automáticamente el tiempo
de la computadora requerido para la simulación.
 
3.2.2 Hoja de trabajo de área de elevación

Esta hoja de trabajo (Figura 8) se usa para ingresar la información hipsográfica del lago
en forma de una curva de área de elevación. Esta información se ingresa desde la
parte más profunda del lago hasta la elevación de superficie más alta posible. Cuando
se ejecuta el programa, el volumen debajo de cada elevación se calcula y se muestra
en la columna C. También se muestra un gráfico del área y el volumen frente a la
elevación.

Figura 8. La hoja de trabajo del área de elevación que se utiliza para ingresar los datos
hipsográficos del lago.

3.2.3 Hoja de trabajo de datos de entrada

Esta hoja de trabajo se usa para ingresar series de tiempo de los datos de entrada
como series de fechas y valores. El ingreso del lago se ingresa en las columnas A y B,
mientras que las series de tiempo de la temperatura de entrada y las concentraciones
se ingresan en las columnas restantes (Figura 9). El programa usa interpolación lineal
para determinar valores intermedios. Tenga en cuenta que cada serie de tiempo debe
abarcar la serie temporal de simulación como se especifica en las celdas B16 y B17 de

la hoja de trabajo LAKE2K

Figura 9. Hoja de trabajo de datos de entrada.

3.2.4 Hoja de trabajo de salida

Esta hoja de trabajo se usa para ingresar la información requerida para calcular la
salida del lago (Figura 10).

Figura 10. Hoja de trabajo de salida.

Se pueden usar siete enfoques para calcular el flujo de salida. El método deseado se
elige usando el menú desplegable en la celda B6.

1. Flujo (epi) = Flowin. El flujo de salida se establece igual a la entrada. Para este caso,
la salida sale del lago desde el epilimnion. Tenga en cuenta que si el Modo de salida se
deja en blanco, esta es la opción predeterminada.

2. Flujo (meta) = Flowin. El flujo de salida se establece igual a la entrada. Para este
caso, la salida sale del lago desde el metalimnion.

3. Flujo (hypo) = Flowin. El flujo de salida se establece igual a la entrada. Para este
caso, la salida sale del lago desde el hipolimnion.

4. Serie temporal (epi). El usuario puede ingresar una serie temporal de salidas en las
columnas A y B. El programa usa la interpolación lineal para determinar los valores
entre las fechas ingresadas. Para este caso, la salida sale del lago desde el epilimnion.

5. Serie temporal (meta). El usuario puede ingresar una serie temporal de salidas en
las columnas A y B. El programa usa la interpolación lineal para determinar los valores
entre las fechas ingresadas. Para este caso, la salida sale del lago desde el
metalimnion.

6. Serie de tiempo (hipo). El usuario puede ingresar una serie temporal de salidas en
las columnas A y B. El programa usa la interpolación lineal para determinar los valores
entre las fechas ingresadas. Para este caso, la salida sale del lago desde el hipolimnion.
 
7. Etapa-descarga. El usuario puede ingresar una curva de descarga de etapa en forma
tabular en las columnas D y E. El programa usa interpolación lineal para determinar el
intermedio.
Tenga en cuenta que si se elige cualquiera de las opciones de series temporales (4 a
6), la serie debe abarcar la serie temporal de simulación tal como se especifica en las
celdas B16 y B17 de la hoja de trabajo LAKE2K.

3.2.5 Hoja de trabajo de condiciones iniciales

Esta hoja de trabajo se usa para ingresar las condiciones iniciales para la elevación del
lago y las variables del modelo. Como en la Figura 11, los valores se ingresan para
epilimnion, metalimnion e hypolimnion. La elevación corresponde a la profundidad de
la parte superior de cada capa medida hacia arriba desde el fondo del lago (por

definición cero).

Figura 11. Hoja de trabajo de condiciones iniciales.

3.2.6 Hoja de trabajo de mezcla vertical

Esta hoja de trabajo se usa para ingresar información relacionada con la mezcla
vertical entre las capas del lago (Figura 12). Se especifican dos opciones (modos de
mezcla) usando el menú desplegable en la celda b6.
Estos son:

1. Munk-Anderson. Este enfoque emplea un modelo para calcular la mezcla vertical

entre las capas en función de la velocidad del viento y la densidad del agua. Tenga en
cuenta que esta es la opción predeterminada.
 
2. Series de tiempo (abajo). Los coeficientes de difusión turbulenta verticales entre las
capas se pueden ingresar como series de tiempo. El programa usa interpolación lineal
para determinar los coeficientes de mezcla entre las fechas ingresadas.
Figura 12. La hoja de trabajo de Mezcla Vertical.

Tenga en cuenta que si se elige la opción de series de tiempo, la serie debe abarcar la
serie temporal de simulación como se especifica en las celdas B16 y B17 de la hoja de
trabajo LAKE2K.

3.2.1 Hoja de trabajo de meteorología

Esta hoja de trabajo se usa para ingresar series temporales de datos meteorológicos
para el lago (Figura 13). El programa usa interpolación lineal para determinar valores
entre las fechas ingresadas. Tenga en cuenta que cada serie de tiempo debe abarcar la
serie temporal de simulación tal como se especifica en las celdas B16 y B17 de la hoja
de trabajo LAKE2K. Además, tenga en cuenta que el programa calcula internamente y
muestra valores para la radiación solar diaria promedio y el fotoperíodo en las

columnas H e I.
Figura 13. Hoja de trabajo de datos meteorológicos.  

3.2.8 Hoja de trabajo de hielo Esta hoja de trabajo ingresa valores que se usan para calcular el
espesor del hielo y para parametrizar su impacto en la transmisión de la luz al lago (Figura 14).

Figura 14. Hoja de trabajo de los parámetros de hielo.

3.2.9 Hoja de trabajo de luz y calor

Esta hoja de trabajo se utiliza para ingresar información relacionada con el régimen de
luz del lago y el balance de calor.

Figura 15. Hoja de trabajo de luz y calor utilizada para ingresar parámetros relacionados con la luz.
  3.2.9.1 Parámetros relacionados con la luz Radiación fotosintéticamente disponible.

Esta es la fracción de radiación solar entrante que está disponible para la fotosíntesis. Se
recomienda que este valor se establezca en 0.47. Fracción de dispersión hacia adelante. Esta es la
fracción de luz que se dispersa hacia adelante. Se recomienda que este valor se establezca en
0,94 / m. Coeficiente de absorción de agua Este parámetro representa la cantidad de luz absorbida
por el agua. Se recomienda que este valor se establezca en 0.012 / m. Coeficiente de absorción de
color Este parámetro representa la cantidad de luz absorbida por el color. Este valor es específico
del sistema. Constante de absorción de clorofila Este parámetro explica la cantidad de luz
absorbida por la clorofila de fitoplancton. Se recomienda que este valor se establezca en 0.03 m2 /
mgA. Constante de absorción POC. Este parámetro representa la cantidad de luz absorbida por
POC expresada como peso seco. Se recomienda que este valor se establezca en 0.16 m2 / gC.
Constante de absorción de sólidos suspendidos inorgánicos. Este parámetro representa la cantidad
de luz absorbida por los sólidos suspendidos inorgánicos. Se recomienda que este valor se
establezca en 0.031 m2 / gD. Coeficiente de dispersión de agua. Este parámetro representa la
cantidad de luz dispersada por el agua. Se recomienda que este valor se establezca en 0.0027 / m.
Constante de dispersión de clorofila. Este parámetro representa la cantidad de luz dispersada por
la clorofila de fitoplancton. Se recomienda que este valor se establezca en 0.01 m2 / mgA.
Constante de dispersión POC. Este parámetro representa la cantidad de luz dispersada por POC
expresada como peso seco. Se recomienda que este valor se establezca en 0.24 m2 / gD. Sólidos
en suspensión inorgánicos dispersos constantemente. Este parámetro representa la cantidad de
luz dispersada por los sólidos suspendidos inorgánicos. Se recomienda que este valor se
establezca en 0,37 m2 / gD. 3.2.9.2 Parámetros relacionados con el calor Conducción / constante
de evaporación. Coeficiente de viento de conducción / evaporación. Conducción / Evaporación
eólica. 3.2.10 Hoja de cálculo de tarif
Figura 15. Hoja de trabajo de luz y calor utilizada para ingresar parámetros relacionados con la luz.

3.2.9.1 Parámetros relacionados con la luz

Radiación fotosintéticamente disponible. Esta es la fracción de radiación solar entrante que está
disponible para la fotosíntesis. Se recomienda que este valor se establezca en 0.47.
Fracción de dispersión hacia adelante. Esta es la fracción de luz que se dispersa hacia adelante.
Se recomienda que este valor se establezca en 0,94 / m.
Coeficiente de absorción de agua Este parámetro representa la cantidad de luz absorbida por el
agua. Se recomienda que este valor se establezca en 0.012 / m.
Coeficiente de absorción de color Este parámetro representa la cantidad de luz absorbida por el
color. Este valor es específico del sistema.
Constante de absorción de clorofila Este parámetro explica la cantidad de luz absorbida por la
clorofila de fitoplancton. Se recomienda que este valor se establezca en 0.03 m2 / mgA.
Constante de absorción POC. Este parámetro representa la cantidad de luz absorbida por POC
expresada como peso seco. Se recomienda que este valor se establezca en 0.16 m2 / gC.
Constante de absorción de sólidos suspendidos inorgánicos. Este parámetro representa la cantidad
de luz absorbida por los sólidos suspendidos inorgánicos. Se recomienda que este valor se
establezca en 0.031 m2 / gD.
Coeficiente de dispersión de agua. Este parámetro representa la cantidad de luz dispersada por el
agua. Se recomienda que este valor se establezca en 0.0027 / m.
Constante de dispersión de clorofila. Este parámetro representa la cantidad de luz dispersada por
la clorofila de fitoplancton. Se recomienda que este valor se establezca en 0.01 m2 / mgA.
Constante de dispersión POC. Este parámetro representa la cantidad de luz dispersada por POC
expresada como peso seco. Se recomienda que este valor se establezca en 0.24 m2 / gD.
Sólidos en suspensión inorgánicos dispersos constantemente. Este parámetro representa la
cantidad de luz dispersada por los sólidos suspendidos inorgánicos. Se recomienda que este valor
se establezca en 0,37 m2 / gD.

3.2.9.2 Parámetros relacionados con el calor

Conducción / constante de evaporación. Coeficiente de viento de conducción / evaporación.
Conducción / Evaporación eólica.

3.2.10 Hoja de cálculo de tarifas

Esta hoja de trabajo se usa para ingresar los parámetros de velocidad del modelo. Debido a que
consta de muchas entradas, se divide en varias partes (Figura 16 a Figura 19).as Esta hoja de
trabajo se usa para ingresar los parámetros de velocidad del modelo. Debido a que consta de
muchas entradas, se divide en varias partes (Figura 16 a Figura 19).

Figura 16. La parte de la hoja de cálculo utilizada para ingresar la estequiometría y los
parámetros de velocidad para sólidos suspendidos inorgánicos, POC y DOC.

Estequiometría:

El modelo supone una estequiometría fija de plantas y materia detrítica. Los valores
recomendados para estos parámetros se enumeran en la Tabla 1. La clorofila es la más
variable de estos valores con un rango de aproximadamente 0,5 a 2 gA. Se sugiere que
el peso en seco se establezca en 100 gD y los otros valores se escalen en
consecuencia.

Tabla 1 Valores recomendados para estequiometría.

Peso seco 100 gD

Carbono 40 gC
Nitrógeno 7.2 gN
Fósforo 1 gP
Clorofila 1 gA
.
Sólidos suspendidos inorgánicos:

Velocidad de asentamiento

Carbono orgánico particulado:

Tasa de hidrólisis
Parámetro de temperatura para velocidad de hidrólisis Velocidad de sedimentación
 
Carbono orgánico disuelto:

Tasa de oxidación
Parámetro de temperatura para la tasa de oxidación del DOC
Modelo de inhibición de oxígeno para la oxidación de DOC. Se utiliza un menú
desplegable para seleccionar entre dos modelos de inhibición: media saturación y
exponencial.
Parámetro de inhibición de oxígeno

Figura 17. La parte de la hoja de cálculo usada para ingresar parámetros para
nitrógeno orgánico, nitrógeno amoniacal, nitrógeno nitrado, fósforo orgánico, sílice
orgánica y oxígeno disuelto.

Nitrógeno orgánico:
Tasa de hidrólisis
Parámetro de temperatura para la tasa de hidrólisis de nitrógeno orgánico Velocidad de
sedimentación

Nitrógeno amoniacal

Tasa de oxidación
Parámetro de temperatura para la tasa de nitrificación
Modelo de inhibición de oxígeno para la nitrificación. Se utiliza un menú desplegable
para seleccionar entre dos modelos de inhibición: media saturación y exponencial.
 
Parámetro de inhibición de oxígeno Nitrato de nitrógeno:
Tasa de desnitrificación
Parámetro de temperatura para la desnitrificación
Modelo de gatillo de oxígeno para desnitrificación. Se utiliza un menú desplegable para
seleccionar entre dos modelos de inhibición: media saturación y exponencial.
Parámetro de disparo de oxígeno
Constante de media saturación de carbono para desnitrificación Fósforo orgánico:
Tasa de hidrólisis
Parámetro de temperatura para la tasa de hidrólisis de fósforo orgánico Velocidad de
sedimentación

Sílice orgánico:

Tasa de disolución
Parámetro de temperatura para la tasa de hidrólisis de sílice orgánica Velocidad de
sedimentación

Oxígeno disuelto:

Parámetro de temperatura para la reacción Oxígeno por C oxidado

Oxígeno por N nitrificado
Fórmula de Reaeración Se utiliza un menú desplegable para seleccionar entre tres
modelos de transferencia de gas: Banks-Herrera, Wanninkhof y O'Connor. Si se deja en
blanco, el modelo se predetermina a la opción O'Connor.

Figura 18. La parte de la hoja de cálculo utilizada para ingresar los parámetros de la tasa de
fitoplancton.   Fitoplancton: Tasa de crecimiento máxima Modelo de temperatura para el
crecimiento. Se utiliza un menú desplegable para elegir entre el modelo theta o un modelo óptimo.
Theta o Topt. Dependiendo de la elección del modelo de temperatura (fila 55), esto mantendrá el 

para el modelo theta o la temperatura óptima para el modelo óptimo. Kappa1. Este parámetro
controla la frecuencia de disminución del modelo óptimo para temperaturas inferiores a la
temperatura óptima. Kappa2. Este parámetro controla la frecuencia de disminución del modelo
óptimo para temperaturas superiores a la temperatura óptima. Ritmo respiratorio Parámetro de
temperatura para la respiración Tasa de mortalidad Parámetro de temperatura para la muerte Mitad
de saturación de nitrógeno Mitad de saturación de fósforo Media saturación de sílice Modelo ligero
Se utiliza un menú desplegable para seleccionar entre tres modelos ligeros: media saturación,
Steele y Smith. Parámetro de luz Velocidad de asentamiento Parámetro de preferencia de
amoníaco Grazabilidad. Este parámetro especifica si el fitoplancton es comestible por el
zooplancton herbívoro. Puede ir desde 0 (totalmente incomestible) a 1 (totalmente comestible).
Contenido de sílice Esta es la proporción estequiométrica de sílice para el tipo particular de
fitoplancton

Figura 19. La parte de la hoja de cálculo utilizada para ingresar parámetros de tasas
para el zooplancton herbívoro y carnívoro.
 
Zooplancton herbívoro:

Tasa máxima de pastoreo

Parámetro de temperatura para el pastoreo de herbívoros Tasa de respiración
Parámetro de temperatura para la respiración de herbívoros Tasa de mortalidad
Parámetro de temperatura para la muerte de herbívoros Eficiencia de pastoreo
Las algas se sentaron medio Zooplancton carnívoro constante:
Tasa máxima de pastoreo
Parámetro de temperatura para el pastoreo de carnívoros Tasa de respiración
Parámetro de temperatura para la respiración de carnívoros Tasa de mortalidad
Parámetro de temperatura para la muerte de carnívoros Eficiencia de pastoreo
Hierba medio sentada constante

3.3 HOJAS DE TRABAJO DE DATOS

Estas hojas de trabajo se utilizan para ingresar los datos medidos que se muestran en
los gráficos. Como tal, son opcionales. En las versiones actuales de Microsoft Office,
estas hojas de trabajo se distinguen por pestañas de color amarillo pálido.

Tenga en cuenta que estos datos se pueden organizar en cualquier orden. Sin
embargo, no se deben incluir filas en blanco.
Si el programa encuentra una fila en blanco debajo de la primera fila, todos los datos
en las filas inferiores no se guardarán cuando se ejecute el programa.

3.3.1 Hoja de trabajo de datos físicos

Esta hoja de trabajo se usa para ingresar una serie temporal de profundidades del lago
y radiación solar medidas (Figura 20)

Figura 20. Hoja de trabajo de datos físicos. 3.3.2 Hojas de trabajo de datos de calidad del agua
Hay tres hojas de trabajo que se utilizan para ingresar datos de modelo para epilimnion,
metalimnion e hypolimnion. La Figura 21 muestra la hoja de datos de epilimnion.

Figura 21. Hoja de trabajo de datos de Epilimnion.

3.3.3 Hoja de trabajo de datos misceláneos

Esta hoja de trabajo se usa para ingresar datos para las profundidades de Secchi y los
coeficientes de extinción de la luz (Figura 22).
Figura 22. Hoja de trabajo de datos misceláneos.

3.4 Hojas de trabajo de salida

Estas son una serie de hojas de trabajo que presentan tablas de resultados numéricos
de series de tiempo generadas por L2K. La mayoría de esta información se muestra en
forma gráfica en gráficos posteriores. En las versiones actuales de Microsoft Office,
estas hojas de trabajo se distinguen por pestañas de color verde pálido.

Hoja de trabajo de salida hidráulica. Esta hoja de trabajo resume las series temporales
de flujo, volumen, profundidad y espesor del hielo generadas por el modelo (Figura 23).

Figura 23. Hoja de trabajo de salida de hidráulica. Hoja de trabajo de salida de variables de calidad
del agua. Esta hoja de trabajo presenta las series temporales de todas las variables del modelo
para cada capa del lago (Figura 24).
Figura 25. La hoja de trabajo de variables varias.

Hoja de trabajo de salida de variables de agua de sedimento. Esta hoja de trabajo

presenta las series de tiempo de las variables sedimento-agua que corresponden a los
sedimentos debajo de cada capa. Estos incluyen SOD y flujos de nutrientes. Tenga en
cuenta que un SOD positivo significa una pérdida de oxígeno de la capa a los
sedimentos. Por el contrario, un flujo de nutrientes positivo significa una transferencia
de nutrientes desde los sedimentos al agua.

Figura 26. Hoja de trabajo de salida de flujo de sedimento.

 
Documentación

Esta sección proporciona una descripción de la teoría del modelo.

 
4 FÍSICA
El modelo actualmente simula el lago como un sistema unidimensional que consta de
tres capas verticales (Figura 27). Los volúmenes de las dos capas más profundas
(metalimnion e hypolimnion) se mantienen fijos, mientras que el epilimnion se permite
cambiar en función del equilibrio entre las entradas y salidas.

Qp Qe
A0
Qout
epi (1) H1
Qin Temp
 Qout H2 H
meta (2)

A2

hypo (3)
Qout H3
Figura 27 El balance hídrico LAKE2K y el esquema de segmentación vertical. El flujo de entrada se
enruta hacia una de las capas en función de la densidad del flujo de entrada y la capa de acuerdo
con el algoritmo: • Si la densidad de entrada es menor que el promedio de las densidades de
epilimnion y metalimnion, la afluencia ingresa al epilimnion. • Si la densidad de entrada es mayor
que el promedio de las densidades de metalimnion e hipolimnion, el flujo entra al hipolimnion. • De
lo contrario, la entrada entra al metalimnion. 4.1 Equilibrio del agua Se calcula un balance hídrico
dinámico para el lago como en

dV in
 Q dt

donde V = volumen del lago [m3], t = tiempo (d), Qin = flujo de entrada [m3 / d], Qp =
flujo de precipitación [m3 / d], Qe
= flujo de evaporación [m3 / d], y Qout = flujo de salida [m3 / d]. Esta ecuación está
integrada para simular cómo cambia el volumen del lago en función del tiempo. Cada
término en el saldo se describe a continuación.

4.1.1 Entrada

El flujo de entrada representa la suma de todos los flujos que ingresan al lago desde los
afluentes, y las fuentes directas puntuales y no puntuales. El usuario ingresa estos
valores como una serie única de fechas y flujos.

4.1.2 Precipitación

 
El flujo de precipitación se ingresa como una serie temporal de fechas y tasas de
precipitación. El programa luego calcula el flujo de precipitación como

Qp  PA0
 

(2)
 

donde P = la tasa de precipitación [m / d], y A0 = la superficie del lago [m / d].

4.1.3 Evaporación

El modelo calcula internamente la pérdida de agua debido a la evaporación. Este flujo

se calcula como

 
Q  J e A0 
Le
 
m 100 cm
 

(3)
 
 
donde Je = el flujo de calor debido a la evaporación [cal / cm2 / d],  = la densidad del
agua [= 1 g / cm3], y
Le = el calor latente de vaporización [cal / g]. El calor latente está relacionado con la
temperatura

Le  597.3  0.57T1
donde T1 = la temperatura del epilimnion [oC]. El flujo de calor debido a la evaporación
se calcula como una función de la velocidad del viento, la temperatura del punto de
rocío y la temperatura del epilimnio como se describe a continuación en la Sec. 5.1.5.

4.1.4 Salida

Siete opciones están disponibles para calcular el flujo de salida. El método deseado se
elige usando el menú desplegable en la celda B6 en la hoja de trabajo de salida (Figura
10).

• Flowout (epi) = Flowin. El flujo de salida se establece igual a la entrada. Para este
caso, la salida sale del lago desde el epilimnion. Tenga en cuenta que si el Modo de
salida se deja en blanco, esta es la opción predeterminada.

• Flujo (meta) = Flowin. El flujo de salida se establece igual a la entrada. Para este
caso, la salida sale del lago desde el metalimnion.

Flujo (hypo) = Flowin. El flujo de salida se establece igual a la entrada. Para este caso,
la salida sale del lago desde el hipolimnion.

• Serie de tiempo (epi). El usuario puede ingresar una serie temporal de salidas en las
columnas A y B. El programa usa la interpolación lineal para determinar los valores
entre las fechas ingresadas. Para este caso, la salida sale del lago desde el epilimnion.

• Serie de tiempo (meta). El usuario puede ingresar una serie temporal de salidas en
las columnas A y B. El programa usa la interpolación lineal para determinar los valores
entre las fechas ingresadas. Para este caso, la salida sale del lago desde el
metalimnion.

• Serie de tiempo (hipo). El usuario puede ingresar una serie temporal de salidas en las
columnas A y B. El programa usa la interpolación lineal para determinar los valores
entre las fechas ingresadas. Para este caso, la salida sale del lago desde el hipolimnion.

• Etapa-descarga. El usuario puede ingresar una curva de descarga de etapa en forma

tabular en las columnas D y E. El programa usa interpolación lineal para determinar
valores intermedios.

Tenga en cuenta que si se elige cualquiera de las opciones de series temporales, la

serie debe abarcar la serie temporal de simulación tal como se especifica en las celdas
B16 y B17 de la hoja de trabajo LAKE2K

4.2 Segmentación vertical y mezcla El lago se divide verticalmente en tres capas, como se muestra
en la Figura 27. El metalimnion y el hypolimnion se establecen en volúmenes fijos, mientras que el
volumen de epilimnion puede variar. El ingreso se ingresa en una de las capas en función de la
densidad del flujo de entrada y de la densidad de las capas. Si es necesario (es decir, si el flujo
entra en el hypolimnion o metalimnion), el flujo se enruta a través de las otras capas y sale a través
del epilimnion o el hypolimnion dependiendo de la elección del modo de flujo de salida del usuario.
La densidad está relacionada con la temperatura del agua y los sólidos disueltos mediante una
ecuación de estado como se describe en el Apéndice A. El modelo usa difusión turbulenta para
mezclar agua entre las capas. Por ejemplo, para una sustancia conservadora sin flujo, se puede
escribir un balance de masa para el metalimnion como   ' dc2  E         'E'          1 c1 c2 dt V2  
2 c3 c2 V2   (5)   donde ci = la concentración en la capa i [g / m3], i = 1, 2 y 3 para el epi-, meta- e
hypolimnion, respectivamente, Ei '= el coeficiente de difusión turbulenta a granel a través del límite
inferior de la capa i [ m3 / d], y Vi = el volumen de la capa i [m3]. El coeficiente de dispersión
masiva está relacionado con cantidades más fundamentales por   E '   (Hola   Ei Ai  Hola   ) / 2  
(6)   donde Ei = el coeficiente de difusión turbulento a través del límite inferior de la capa i [m2 / d],
Ai = el área superficial del límite inferior de la capa i [m2], y Hi = el espesor de la capa i [m] (
4.2 Segmentación vertical y mezcla

El lago se divide verticalmente en tres capas, como se muestra en la Figura 27. El metalimnion y el
hypolimnion se establecen en volúmenes fijos, mientras que el volumen de epilimnion puede variar.
El ingreso se ingresa en una de las capas en función de la densidad del flujo de entrada y de la
densidad de las capas. Si es necesario (es decir, si el flujo entra en el hypolimnion o metalimnion),
el flujo se enruta a través de las otras capas y sale a través del epilimnion o el hypolimnion
dependiendo de la elección del modo de flujo de salida del usuario. La densidad está relacionada
con la temperatura del agua y los sólidos disueltos mediante una ecuación de estado como se
describe en el Apéndice A.

El modelo usa difusión turbulenta para mezclar agua entre las capas. Por ejemplo, para una
sustancia conservadora sin flujo, se puede escribir un balance de masa para el metalimnion como

'
dc2  E  

 'E'   

1 c1 c2
dt V2

2 c3 c2
V2

(5)

donde ci = la concentración en la capa i [g / m3], i = 1, 2 y 3 para el epi-, meta- e hypolimnion,

respectivamente, Ei '= el coeficiente de difusión turbulenta a granel a través del límite inferior de la
capa i [ m3 / d], y Vi = el volumen de la capa i [m3]. El coeficiente de dispersión masiva está
relacionado con cantidades más fundamentales por

E '

(Hola

Ei Ai
 Hola
)/2

(6)

donde Ei = el coeficiente de difusión turbulento a través del límite inferior de la capa i [m2 / d], Ai =
el área superficial del límite inferior de la capa i [m2], y Hi = el espesor de la capa i [m] (vea la
figura 27).

Hay dos opciones (modos de mezcla) disponibles para determinar los coeficientes de difusión
turbulenta tal como se especifica utilizando el menú desplegable en la celda b6 de la Hoja de
trabajo de mezcla vertical (Figura 12).
Estos son:

Munk-Anderson. Este enfoque emplea un modelo para calcular la mezcla vertical entre capas
basada en la velocidad del viento y la densidad del agua desarrollada originalmente por Munk y
Anderson (1948). El coeficiente de difusión turbulenta vertical se calcula comovea la figura 27).  
Hay dos opciones (modos de mezcla) disponibles para determinar los coeficientes de difusión
turbulenta tal como se especifica utilizando el menú desplegable en la celda b6 de la Hoja de
trabajo de mezcla vertical (Figura 12). Estos son: Munk-Anderson. Este enfoque emplea un modelo
para calcular la mezcla vertical entre capas basada en la velocidad del viento y la densidad del
agua desarrollada originalmente por Munk y Anderson (1948). El coeficiente de difusión turbulenta
vertical se calcula como

E0
E

(1  aRi )0.5
donde E0 = el coeficiente de difusión máximo [m2 / d], a = un coeficiente de sintonía,
y Ri = el número de Richardson.

 
El coeficiente de difusión máximo es proporcional a la velocidad de corte (Sundaram y
Rehm 1973),

E0  cU * (8)
donde c = un coeficiente de ajuste y U * = la velocidad de corte [m / s]. La velocidad
de corte se calcula como

u=

donde air = la densidad del aire (se supone igual a 1.2 kg / m3), Cd = un coeficiente
de resistencia adimensional (se supone igual a 1.3103), Uw = la velocidad del viento
[m / s], y  = la densidad del agua [kg / m3]. En nuestro caso,  es la densidad del agua
en el epilimnio del lago.
El número de Richardson es un número adimensional que representa la relación entre
la flotabilidad y las fuerzas de corte. Ford y Johnson (1986) han propuesto la siguiente
fórmula para calcular el número de Richardson,

 
R   (g / ) ( / z)
  
 

(10)
 
i U * / (z  z) 2

donde g = la aceleración debida a la gravedad (= 9.81 m / s2),  = la densidad del

agua en la profundidad a la que se calcula el número de Richardson [kg / m3],  / z =
el gradiente de densidad con profundidad, zs
= la elevación de la superficie del agua yz = la elevación a la que se calcula el número
de Richardson [m].

 
El gradiente de densidad en la interfaz entre las capas se calcula con una diferencia
dividida finita

  i  i1
  
 

(11)
 
Zz zi  zi1
 
donde i = la densidad para la capa i [kg / m3] y zi = la elevación en el punto medio de
la capa i [m].
El modelo se ajusta ajustando los coeficientes a y c hasta que haya una coincidencia
aceptable entre las temperaturas simuladas y las medidas para cada capa del modelo.
El valor predeterminado para a es 30 con un rango de 10 a 60. Un valor
predeterminado para c se estima como

 
c  H max
34
 

(12)
 

donde Hmax = la profundidad máxima del lago [m].

Series de tiempo (abajo). Los coeficientes de difusión turbulenta verticales entre las
capas se pueden ingresar como series de tiempo. El programa usa interpolación lineal
para determinar los coeficientes de mezcla entre las fechas ingresadas.

5 MODELO DE TEMPERATURA Se escribe un balance de calor para cada capa vertical. Para el
epilimnion para el caso donde el flujo de entrada y salida pasa a través del epilimnion y no hay
cubierta de hielo, el saldo se escribe como
(13)
donde Heat1 = calor del epilimnion [cal], t = tiempo [d],  = la densidad del agua [g / cm3], Cp = el
calor específico del agua [cal / (g oC)], Tin = temperatura de la afluencia [oC], Ti = temperatura en
la capa i [oC], Tair = temperatura del aire [oC], E'i = el coeficiente de dispersión masiva a través del
límite inferior de la capa i [m3 / d], Ai = el área en la parte inferior de la capa i [m2], y Jh = el flujo de
calor aire-agua [cal / (cm2 d)]. Se escriben balances similares para las otras capas con la notable
excepción de que se omite el flujo de calor superficial. 5.1 Flujo de calor superficial Como se
muestra en la Figura 28, el intercambio de calor superficial se modela como una combinación de
cinco procesos:   J h    I (0)  J an  J br  Jc  Je   (14)   donde I (0) = radiación de onda corta
solar neta en la superficie del agua, Jan = radiación de onda larga atmosférica neta, Jbr = radiación
de retorno de onda larga del agua, Jc = conducción y Je = evaporación. Todos los flujos se
expresan como cal / cm2 / d.
Figura 28 Componentes del intercambio de calor superficial

1 radiación solar

El modelo calcula la cantidad de radiación solar que ingresa al agua en una latitud
particular (Lat) en la superficie de la tierra. Esta cantidad es una función de la radiación
en la parte superior de la atmósfera terrestre que se atenúa por la transmisión
atmosférica, la cobertura de nubes y la reflexión.

 
I (0)  I 0
 
en ac
 
(1  Rs)
 

radiación extraterrestre
 

atenuación atmosférica
 

atenuación de la nube
 

reflexión
 
(15)
 
donde I (0) = radiación solar en la superficie del agua [cal / cm2 / d], I0 = radiación
extraterrestre (es decir, en la parte superior de la atmósfera terrestre) [cal / cm2 / d],
at = atenuación atmosférica, CL = fracción del cielo cubierto con nubes, y Rs = albedo
(fracción reflejada).

 
Radiación extraterrestre La radiación extraterrestre se calcula como (TVA 1972)

Yo  W0 pecado 
0r2
donde W0 = la constante solar [1367 W / m2 o 2823 cal / cm2 / d], r = radio
normalizado de la órbita de la tierra (es decir, la relación de la distancia tierra-sol real a
la distancia media tierra-sol), y  = la altitud del sol [radianes], que se puede calcular
como

pecado   pecado  pecado Lat  cos cos Lat

donde  = declinación solar [radianes], y Lat = latitud local [radianes].
El radio normalizado se puede estimar como
 

(dieciséis)
Ecuaciones 17.18

donde Dy = día juliano (1 de enero = 1, 2 de enero = 2, etc.). La declinación solar puede estimarse
como
ecuación 19

El fotoperíodo f [horas] se calcula como