Rendimiento hídrico estacional

Resumen

Existe una gran demanda de herramientas que estimen el efecto de la gestión del paisaje en el servicio de suministro de agua para usos como el riego, el consumo doméstico y la producción de energía hidroeléctrica. Aunque el modelo de rendimiento hídrico anual de InVEST proporciona una estimación del rendimiento hídrico total de una cuenca, muchas aplicaciones requieren conocer los caudales estacionales, especialmente durante la estación seca. Esto requiere la comprensión de los procesos hidrológicos en una cuenca, en particular la partición entre el flujo rápido (que se produce durante o poco después de los eventos de lluvia) y el flujo de base (que se produce durante el tiempo seco). En climas muy estacionales, es probable que el flujo de base tenga más valor que el flujo rápido, a menos que se disponga de un almacenamiento significativo (por ejemplo, un gran embalse). El modelo de rendimiento hídrico estacional de InVEST pretende ofrecer una orientación sobre la contribución de las parcelas a la generación de flujos de base y flujos rápidos. El modelo calcula índices espaciales que cuantifican la contribución relativa de una parcela de tierra a la generación tanto del caudal base como del caudal rápido. En la actualidad, no existen estimaciones cuantitativas del flujo de base, solo las contribuciones relativas de los píxeles; se está desarrollando una herramienta independiente para abordar esta cuestión.

Introducción

Comprender el efecto de la gestión del paisaje en el flujo estacional es de vital importancia para la gestión de las cuencas hidrográficas. La contribución de una parcela determinada al flujo de la corriente depende de una serie de factores ambientales, como el clima, el suelo, la vegetación, la pendiente y la posición a lo largo de la trayectoria del flujo (lo que determina si el píxel puede recibir agua de la ladera superior o si el agua recargada puede ser evapotranspirada posteriormente).

El agua que fluye por el paisaje se evapora, se transpira, se extrae de un pozo o sale de la cuenca en forma de flujo de aguas subterráneas profundas o de corrientes. Si consideramos un píxel individual, y su valor con respecto al rendimiento hídrico, podemos considerar dos enfoques:

  • El primero da crédito a la cantidad neta de agua generada en un píxel como igual a la precipitación entrante menos las pérdidas por evapotranspiración en ese píxel. En este esquema, es posible que la evapotranspiración real sea mayor que la precipitación si el agua se suministra al sitio desde la ladera. Así, la generación neta podría ser negativa. Este enfoque no tiene en cuenta la disposición final del agua generada en ese píxel; es decir, no considera si el agua aparece realmente como flujo de corriente o se evapora o se retira en algún punto de su recorrido.

  • El segundo enfoque da crédito al agua de una parcela que realmente aparece como flujo de corriente. Así, si una parcela genera agua que luego se evapora, la contribución se considera nula.

El primer enfoque hace más hincapié en el uso y la cobertura del suelo de un lugar, ya que se centra en la generación neta de ese píxel. Tiene en cuenta la subvención del agua de los píxeles pendiente arriba, pero no considera los efectos pendiente abajo. Representa el potencial de generar flujo de corriente (no una generación real de flujo).

El segundo enfoque pone más énfasis en la posición topográfica de un píxel, ya que eso determinará el potencial de que el agua generada en ese píxel se consuma antes de convertirse en caudal. Representa el caudal real generado por un píxel. Como el caudal real no puede ser inferior a cero, este enfoque, a diferencia del primero, dará lugar a índices mayores o iguales a cero.

Utilizamos ambos conceptos para desarrollar un conjunto de tres índices, uno para el flujo rápido, otro para la recarga (que representa el «flujo de base potencial») y otro para el flujo de base real. Aquí, el flujo de base se define como la generación de flujo de corriente con tiempos de residencia de la cuenca de meses a años, mientras que el flujo rápido representa la generación de flujo de corriente con tiempos de residencia de la cuenca de horas a días.

El modelo

Flujo Rápido

El flujo rápido (FR) se calcula con un enfoque basado en el número de curva (NC). Los eventos de lluvia mensuales hacen que la precipitación caiga sobre el paisaje. Las propiedades del suelo y de la cubierta del terreno determinan la cantidad de lluvia que se escurre rápidamente por la superficie del terreno (produciendo un flujo rápido) frente a la que se infiltra en el suelo (produciendo una recarga local). El número de curva es una forma sencilla de captar estas propiedades del suelo y de la cubierta del terreno: los valores más altos de NC tienen un mayor potencial de escorrentía (por ejemplo, los suelos arcillosos y la cubierta vegetal escasa), y los valores más bajos tienen más probabilidades de infiltrarse (por ejemplo, los suelos arenosos y la cubierta vegetal densa).

Para calcular el flujo rápido, utilizamos la profundidad media de los eventos, \(frac{P_{i,m}}{n_{i,m}}\) y asumimos una distribución exponencial de las profundidades de las precipitaciones diarias en los días con lluvia,

\[f\left( p \right) = \frac{1}{a_{i,m}}exp\left( - \frac{p}{a_{i,m}} \right)\]

Donde \(a_{i,m} = \frac{P_{i,m}}{n_{m}}/25.4\) y

  • \(a_{i,m}\) es la profundidad media de la lluvia en un día lluvioso en el píxel

    i on month m [in],

  • \(n_{i,m}\) es el número de eventos en el píxel i en el mes m

    [-],

  • \(P_{i,m}\) es la precipitación mensual para el píxel i en el mes

    m [mm].

El flujo rápido para los píxeles situados en las corrientes se ajusta a la precipitación en ese píxel, lo que supone que no hay infiltración, solo escorrentía.

\[\text{QF}_{stream,m} = \ P_{stream,m}\]

si no, se puede demostrar a partir de la distribución exponencial que la escorrentía mensual \(\text{QF}_{i,m}\) es

(85)\[\text{QF}_{i,m} = n_{m} \times \left( \left( a_{i,m} - S_{i} \right)\exp\left( - \frac{0.2S_{i}}{a_{i,m}} \right) + \frac{S_{i}^{2}}{a_{i,m}}\exp\left( \frac{0.8S_{i}}{a_{i,m}} \right)E_{1}\left( \frac{S_{i}}{a_{i,m}} \right) \right) \times \left( 25.4\ \left\lbrack \frac{\text{mm}}{\text{in}} \right\rbrack \right)\]

donde

  • \(S_{i} ` el el máximo de retención potencial, :math:\)frac{1000}{text{CN}_{i}} - 10` [in]

  • \(\text{CN}_{i}\) is the curve number for pixel i

    [in-1], tabulados en función del LULC local, y del tipo de suelo (el Apéndice I tiene una plantilla de esta tabla),

  • \(E_{1}\) es la función integral exponencial,

    \(E_{1}(x) = \int_{1}^{\infty}{\frac{e^{-xt}}{t}\text{dt}}\).

  • y \(25.4\) es un factor de conversión de pulgadas (utilizado por la ecuación) a milímetros (utilizado por el modelo)

(véase Guswa et al. 2018).

Algunos casos extremos se tratan de forma especial:

  • Cuando \(S_{i} = 0\), el término \(E_{1}\) llega a infinito. \(\text{QF}_{i,m}\) se pone a cero en este caso.

  • Para evitar problemas con la inestabilidad numérica cuando el resultado de exp se hace muy grande, cuando \(\frac{S_{i}}{a_{i,m}} > 100\), redondeamos \(\text{QF}_{i,m}\) a cero.

  • Con ciertas combinaciones de inputs, es posible que la ecuación \(\text{QF}_{i,m}\) anterior se evalúe a un pequeño número negativo. En estos casos \(\text{QF}_{i,m}\) se redondea a cero.

Así, el flujo rápido anual \(\text{QF}_{i}\), puede calcularse a partir de la suma de los valores mensuales de \(\text{QF}_{i,m}\),

(86)\[\text{QF}_{i} = \sum_{m = 1}^{12}{QF_{i,m}}\]

Recarga local

La recarga local, o la contribución potencial al flujo de base, de un píxel se calcula a partir del balance hídrico local. La precipitación que no se escurre como flujo rápido, y que no es evapotranspirada por la vegetación de un píxel, puede infiltrarse en el suelo para convertirse en recarga local. La recarga local puede ser negativa si un píxel no recibe suficiente agua propia para satisfacer sus necesidades de vegetación (determinadas por su factor de cultivo Kc), por lo que utiliza también el agua generada en la ladera del píxel (lo que se denomina «subsidio pendiente arriba»). El índice de recarga local se calcula en una escala temporal anual, pero utiliza valores derivados de los presupuestos hídricos mensuales.

Para un píxel i, la recarga local derivada del balance hídrico anual es (Figura 1):

(87)\[L_{i} = P_{i} - \text{QF}_{i} - \text{AET}_{i}\]

Donde la evapotranspiración real anual ETR es la suma de la ETR mensual:

(88)\[\text{AET}_{i} = \sum_{\text{months}}^{}\text{AET}_{i,m}\]

Para cada mes, \(\text{AET}_{i,m}\) está limitado por la demanda (evapotranspiración potencial - ETP) o por el agua disponible (de Allen et al. 1998):

(89)\[\text{AET}_{i,m} = min(\text{PET}_{i,m}\ ;\ P_{i,m} - \text{QF}_{i,m} + \alpha_{m}\beta_{i}L_{sum.avail,i})\]

donde \(\text{PET}_{i,m}\) es la evapotranspiración potencial mensual,

(90)\[\text{PET}_{i,m} = K_{c,i,m} \times ET_{0,i,m}\]

\(L_{sum.avail,i}\) se define recursivamente por (Figura 2),

(91)\[L_{sum.avail,i} = \sum_{j \in \{ neighbor\ pixels\ draining\ to\ pixel\ i\}}^{}{p_{\text{ij}} \cdot \left( L_{avail,j} + L_{sum.avail,j} \right)}\]

donde \(p_{text{ij}\) en lbrack 0,1rbrack` es la proporción de flujo de la celda i a j, y \(L_{avail,i}\) es la recarga disponible para un píxel, que es : math:L_{i} cuando \(L_{i}\) es negativo, y una proporción \(\gamma\) de \(L_{i}\) cuando es positivo (véase más abajo la definición de \(\gamma\)):

(92)\[L_{avail,i}\ = min(\gamma L_{i},L_{i})\]

En lo anterior:

  • \(P_{i}\) y \(P_{i,m}\) son la precipitación

    anual y mensual, respectivamente [mm].

  • \(text{QF}_{i}\) y \(text{QF}_{i,m}\) son los índices de flujo rápido

    definidos anteriormente [mm]

  • \(ET_{0,i,m}\) es la evapotranspiración de referencia para el mes m

    [mm]

  • \(K_{c,i,m}\) es el factor de cultivo mensual para el LULC del píxel

  • \(\alpha_{m}\) es la fracción de la recarga anual pendiente arriba

    disponible en el mes m (por defecto es 1/12)

  • \(\beta_{i}\) es la fracción del subsidio pendiente arriba que está

    disponible para la evapotranspiración pendienre abajo (por defecto es 1; véase el Apéndice II con más información)

  • γ es la fracción de la recarga del píxel que está disponible para los píxeles pendiente abajo

    (por defecto es 1)

Atribución de la recarga

El flujo de base total, \(Q_b\) (en mm), es la media de las recargas locales contribuyentes (negativas o positivas) en la cuenca,

(93)\[Q_{b} = \frac{\sum_{k \in \left\{ \text{pixels in catchment} \right\}}^{}L_{k}}{n_{\text{pixels in catchment}}}\]

El valor de atribución a un píxel es la contribución relativa de la recarga local \(L\) en ese píxel al flujo de base \(Q_b\):

(94)\[V_{R,i} = \frac{L_{i}}{{Q_{b} \times n}_{\text{pixels in catchment}}}\]

_images/fig1.png

Figura 1. Balance hídrico a escala de píxel para calcular la recarga local (Ecuación 3), donde Bsum es el flujo que realmente llega a la corriente.


_images/fig2.png

Figura 2. Recorrido a escala de ladera para calcular la evapotranspiración real (basada en las variables climáticas de cada píxel y en la contribución penduente arriba, véase la Ecuación 5) y el flujo de base (basado en Bsum, el flujo que realmente llega a la corriente, véanse las Ecuaciones 11-14).


Flujo base

El índice de flujo base representa la contribución de un píxel al flujo base (es decir, el agua que llega a la corriente durante la estación seca). Si la recarga local es negativa, entonces el píxel no contribuyó al flujo de base, por lo que \(B\) se establece en cero. Si el píxel contribuyó a la recarga de agua subterránea, entonces \(B\) es una función de la cantidad de flujo que sale del píxel y de la contribución relativa a la recarga de este píxel.

Para un píxel que no es adyacente al canal de la corriente, el flujo de base acumulado, \(B_{sum,i}\), es proporcional al flujo de base acumulado que sale de los píxeles adyacentes pendiente abajo menos el flujo de base acumulado que se generó en ese mismo píxel pendiente abajo (Figura 2):

(95)\[\begin{split}B_{sum,i} = L_{sum,i}\sum_{j \in \{\text{cells to which cell i pours}\}}^{}\begin{Bmatrix} p_{\text{ij}}\left( 1 - \frac{L_{avail,j}}{L_{sum,j}} \right)\frac{B_{sum,j}}{L_{sum,j} - L_{j}}\ \text{ if }j\text{ is a nonstream pixel} \\ p_{\text{ij}}\ \text{ if }j\text{ is a stream pixel} \\ \end{Bmatrix}\end{split}\]

En la salida de la cuenca (o en cualquier píxel adyacente a la corriente), la suma de la generación de flujo de base \(B_{sum,i}\) sobre todos los píxeles pendiente arriba es igual a la suma de la generación local sobre los mismos píxeles (porque no hay más oportunidad para que el flujo lento se consuma antes de llegar a la corriente):

(96)\[B_{sum,outlet} = L_{sum,outlet}\]

donde \(L_{sum,i}\) es la recarga acumulada aguas arriba definida por

(97)\[L_{sum,i} = L_{i} + \sum_{j,\ all\ pixels\ draining\ to\ pixel\ i}^{}{L_{sum,j} \cdot p_{\text{ji}}}\]

y el flujo de base, \(B_{i}\) puede derivarse directamente de la proporción del flujo de base acumulado que sale de la celda i, con respecto a la recarga disponible a la recarga acumulada aguas arriba:

(98)\[B_{i} = max\left(B_{sum,i} \cdot \frac{L_{i}}{L_{sum,i}}, 0\right)\]

Calibración

Se recomienda siempre validar con datos observados si es posible. Sin embargo, si bien el resultado del flujo rápido del modelo puede utilizarse como una medida cuantitativa, el flujo de base está pensado para ser utilizado como un índice, no como un valor absoluto. Por lo tanto, es difícil combinar el flujo rápido y el flujo de base y esperar obtener resultados realistas del modelo para validar contra el flujo observado. Una posibilidad es validar los valores relativos (es decir, la distribución de los valores en el paisaje). Esto requiere varios (al menos >3, más realista >5) medidores de corrientes, que pueden ser comparados con los resultados de flujo rápido y flujo base del modelo, agregados a los mismos puntos de medidores de corrientes. Alternativamente, los resultados pueden ser comparados con un modelo espacialmente explícito diferente, si está disponible.

Si intenta validar cuantitativamente el flujo rápido o una combinación de flujo rápido y flujo de base (de nuevo, no se recomienda, pero la gente lo intenta), tenga en cuenta que, dado que los resultados están en milímetros, si simplemente los sumamos en toda el área, es probable que el resultado sea de órdenes de magnitud demasiado grandes y no represente adecuadamente el volumen total de agua. En lugar de ello, utilice el valor medio de B o Qf en toda la cuenca, convierta los milímetros en metros, y luego multiplique por el área de la cuenca para obtener un valor en metros cúbicos, que puede ser comparado con los datos de flujo observados. Alternativamente, se puede calcular el volumen por píxel y sumarlos.

El documento de Hamel et al. (2020) posee un ejemplo de calibración del modelo de Rendimiento Hídrico Estacional contra los datos observados y otros modelos hidrológicos. Para una orientación más general sobre la evaluación de la incertidumbre en los análisis de los servicios ecosistémicoss, véase Hamel y Bryant (2017).

Limitaciones y simplificaciones

Como todos los modelos InVEST, aquí también se utiliza un enfoque simplificado para estimar el caudal rápido y el caudal base, y no incluye muchas de las complejidades que se producen cuando el agua se mueve a través de un paisaje. El caudal rápido se basa principalmente en el número de curva, que no tiene en cuenta la topografía. En cuanto al flujo de base, aunque el modelo utiliza un enfoque basado en la física, las ecuaciones están extremadamente simplificadas tanto a escala espacial como temporal, lo que aumenta significativamente la incertidumbre sobre las cifras absolutas producidas. Por ello, no sugerimos que se utilicen los valores absolutos, sino los valores relativos a través de los paisajes (donde suponemos que las simplificaciones importan menos, porque se aplican a todo el paisaje).

Necesidades de datos

Nota

Todos los inputs espaciales deben tener exactamente el mismo sistema de coordenadas proyectadas (con unidades lineales de metros), no un sistema de coordenadas geográficas (con unidades de grados).

Nota

Los inputs raterizados pueden tener diferentes tamaños de celda, y serán remuestreadas para que coincidan con el tamaño de celda del MDE. Por lo tanto, todos los resultados del modelo tendrán el mismo tamaño de celda que el MDE.

  • Workspace (directorio, requerido): The folder where all the model’s output files will be written. If this folder does not exist, it will be created. If data already exists in the folder, it will be overwritten.

  • Sufijo Del Archivo (texto, opcional): Sufijo que se añadirá a todos los nombres de archivos resultantes. Es útil para diferenciar entre las ejecuciones del modelo.

  • Directorio De Precipitación (directorio, condicionalmente requerido): Directory containing maps of monthly precipitation for each month. Only .tif files should be in this folder (no .tfw, .xml, etc files). Required if User-Defined Local Recharge is not selected.

    Se recomienda encarecidamente utilizar las mismas capas de precipitación que se utilizaron para crear los rásteres de input de evapotranspiración. Si se basan en diferentes fuentes de datos de precipitación, se introduce otra fuente de incertidumbre en los datos, y el desajuste podría afectar a los componentes del balance hídrico calculados por el modelo.

    Contents:

    • Precipitación (ráster, unidades: mm/month, requerido): Twelve files, one for each month. File names must end with the month number (1-12). For example, the filenames “precip_1.tif” and “precip1.tif” are both valid names for the month of January.

  • Directorio ET0 (directorio, condicionalmente requerido): Directory containing maps of reference evapotranspiration for each month. Only .tif files should be in this folder (no .tfw, .xml, etc files). Required if User-Defined Local Recharge is not selected.

    Se recomienda encarecidamente que los rásteres de input de evapotranspiración se basen en los misnos datos de precipitación como input para el modelo. Si se basan en diferentes fuentes de datos de precipitación, se introduce otra fuente de incertidumbre en los datos, y el desajuste podría afectar a los componentes del balance hídrico calculados por el modelo.

    Contenido:

    • Reference Evapotranspiration (ráster, unidades: mm/month, requerido): Twelve files, one for each month. File names must end with the month number (1-12). For example, the filenames “et0_1.tif” “evapotranspiration1.tif” are both valid for the month of January.

  • Modelo Digital De Elevación (ráster, unidades: m, requerido): Mapa de elevación sobre el nivel del mar.

  • Uso Del Suelo/Cobertura Del Suelo (ráster, requerido): Map of land use/land cover codes. Each land use/land cover type must be assigned a unique integer code. All values in this raster MUST have corresponding entries in the Biophysical Table.

  • Grupo Hidrológico Del Suelo (ráster, condicionalmente requerido): Mapa de los grupos hidrológicos del suelo. Los píxeles pueden tener los valores 1, 2, 3 o 4, correspondientes a los grupos hidrológicos de suelo A, B, C o D, respectivamente.

  • Área De Interés (vector, polygon/multipolygon, requerido): Un mapa de zonas sobre el que agregar y resumir los resultados finales.

  • Tabla Biofísica (CSV, requerido): Una tabla que asigna cada código LULC a las propiedades biofísicas de la clase LULC correspondiente. Todos los valores de la trama LULC deben tener los inputs correspondientes en esta tabla.

    Una tabla .csv (valores separados por coimas que contenga la información del modelo correspondiente a cada una de las clases de uso del suelo en el ráster LULC. Todas las clases de LULC en el ráster LULC DEBEN tener los valores correspondientes en esta tabla. Cada fila es una clase de uso del suelo/cobertura del suelo y las columnas deben nombrarse y definirse como sigue:

    Columnas:

    • lucode (entero, requerido): Código LULC del ráster LULC. Cada código debe ser un número entero único.

    • cn_[SOIL_GROUP] (número, unidades: sin unidades, requerido): Valores del número de la curva para cada combinación de grupo de suelos y clase LULC. Sustituya [GRUPO_SUELO] por el código de cada grupo de suelo A, B, C, D para que haya una columna para cada grupo de suelo. Los valores del número de la curva deben ser mayores que 0.

    • kc_[MONTH] (número, unidades: sin unidades, requerido): Valores del coeficiente de cultivo/vegetación (Kc) para esta clase LULC en cada mes. Sustituya [MES] por los números del 1 al 12 para que haya una columna para cada mes.

  • Tabla De Eventos De Lluvia (CSV, condicionalmente requerido): A table containing the number of rain events for each month. Required if neither User-Defined Local Recharge nor User-Defined Climate Zones is selected.

    Un evento de lluvia se define como una precipitación >0,1mm.

    Columnas:

    • month (número, unidades: sin unidades, requerido): Los valores son los números del 1 al 12 correspondientes a cada mes, de enero (1) a diciembre (12).

    • events (número, unidades: sin unidades, requerido): El número de eventos de lluvia en ese mes.

  • Umbral De Acumulación De Flujo (número, unidades: número de píxeles, requerido): El número de píxeles pendiente arriba que deben fluir en un píxel antes de que se clasifique como una corriente.

  • Parámetro Alpha_m (texto, condicionalmente requerido): La proporción de la recarga local anual disponible pendiente arriba que está disponible en cada mes. Requerido si no se selecciona la opción de usar la tabla mensual Alpha.

    Default value: 1/12.

  • Parámetro Beta_i (proporción, requerido): La proporción del subsidio ascendente que está disponible para la evapotranspiración descendente.

    Default value: 1.

  • Parámetro Gamma (proporción, requerido): La proporción de la recarga local del píxel que está disponible para los píxeles pendiente abajo.

    Default value: 1.

Opciones avanzadas del modelo

La tabla de eventos de lluvia mensuales es una forma sencilla de proporcionar datos de eventos de lluvia. Esto supone que hay un número de este tipo para toda la cuenca, lo que puede no ser cierto para grandes áreas o áreas con precipitaciones muy variables espacialmente.

Para representar la variabilidad en el número de eventos de lluvia, es posible introducir un mapa de zonas climáticas, y el número asociado de eventos de lluvia para cada zona.

Inputs

  • Zonas Climáticas (avanzado) (verdadero/falso): Utilizar datos de zonas climáticas definidas por quien usa el modelo en lugar de una tabla global de eventos de lluvia.

  • Tabla De Zonas Climáticas (CSV, condicionalmente requerido): Table of monthly precipitation events for each climate zone. Required if User-Defined Climate Zones is selected.

    Columnas:

    • cz_id (entero, requerido): Números de identificación de la zona climática, correspondientes a los valores del mapa de las zonas climáticas.

    • [MONTH] (número, unidades: sin unidades, requerido): El número de eventos de lluvia que se producen en cada mes en esta zona climática. Sustituya [MES] por las abreviaturas de los meses: ene, feb, mar, abr, may, jun, jul, ago, sep, oct, nov y dic, de modo que haya una columna para cada mes.

  • Mapa De Zonas Climáticas (ráster, condicionalmente requerido): Mapa de zonas climáticas. Todos los valores de este raster deben tener los inputs correspondientes en la Tabla de Zonas Climáticas.


El modelo calcula secuencialmente la capa de recarga local, y luego la capa de flujo base a partir de la recarga local. En lugar de que InVEST calcule la recarga local, esta capa podría obtenerse de un modelo diferente (por ejemplo, RHESSys.) Para calcular la contribución del flujo de base basada en su propia capa de recarga, es posible omitir la primera parte del modelo e introducir directamente un mapa de recarga local.

Inputs

  • Capa De Recarga Definida Por Quien Usa El Modelo (avanzado) (verdadero/falso): Utilizar los datos de recarga local definidos por quien usa el modelo en lugar de calcular la recarga local a partir de los otros datos proporcionados.

  • Recarga Local (ráster, unidades: mm, condicionalmente requerido): Mapa de datos de recarga local. Se necesita si se selecciona la opción de recarga local definida por quien usa el modelo.


El parámetro alpha representa la variabilidad temporal de la contribución del agua disponible en la ladera a la evapotranspiración en un píxel. En la parametrización por defecto, su valor se establece en 1/12, asumiendo que el suelo amortigua la liberación de agua y que la contribución mensual es exactamente 112th de la contribución anual.

Para permitir que la subvención de la pendiente sea temporalmente variable, usted puede proporcionar una tabla de valores mensuales de alpha.

Inputs

  • Utilizar La Tabla Alfa Mensual (avanzado) (verdadero/falso): Utilizar valores alfa mensuales en lugar de un único valor para todo el año.

  • Tabla Alfa Mensual (CSV, condicionalmente requerido): Table of alpha values for each month. Required if Use Monthly Alpha Table is selected.

Interpretación de los resultados

La resolución de los rásteres de resultados será la misma que la del MDE que se proporciona como input.

  • Carpeta [Espacio de trabajo]:

  • Registro de parámetros: Cada vez que se ejecute el modelo, se creará un archivo de texto (.txt) en el Espacio de Trabajo. El archivo enumerará los valores de los parámetros y los mensajes de salida para esa ejecución y se nombrará según el servicio, la fecha y la hora. Cuando se ponga en contacto con NatCap por errores en una ejecución del modelo, incluya el registro de parámetros.

  • B_[Sufijo].tif (tipo: ráster; unidades: mm, pero deben interpretarse como valores relativos, no absolutos): Mapa de valores de flujo base \(B\), la contribución de un píxel al flujo de liberación lenta (que no se evapotranspira antes de llegar a la corriente)

  • B_sum_[Sufijo].tif (tipo: ráster; unidades: mm, pero deben interpretarse como valores relativos, no absolutos): Mapa de valores \(B_{text{sum}\), el flujo que atraviesa un píxel, aportado por todos los píxeles de la ladera, que no es evapotranspirado antes de llegar a la corriente

  • CN_[Sufijo].tif (tipo: ráster): Mapa de los valores del número de la curva

  • L_avail_[Sufijo].tif (tipo: ráster; unidades: mm, pero deben interpretarse como valores relativos, no absolutos): Mapa de la recarga local disponible \(L_{\text{avail}}\)

  • L_[Sufijo].tif (tipo: ráster; unidades: mm, pero deben interpretarse como valores relativos, no absolutos): Mapa de valores de recarga local \(L\).

  • L_sum_avail_[Sufijo].tif (tipo: ráster; unidades: mm, pero deben interpretarse como valores relativos, no absolutos): Mapa de los valores de \(L_{text{sum.avail}\), el agua disponible para un píxel, aportada por todos los píxeles de ladera, que está disponible para la evapotranspiración de este píxel

  • L_sum_[Sufijo].tif (tipo: ráster; unidades: mm, pero deben interpretarse como valores relativos, no absolutos): Mapa de los valores de \(L_{text{sum}\), el flujo a través de un píxel, aportado por todos los píxeles pendiente arriba, que está disponible para la evapotranspiración a los píxeles pendiente abajo

  • QF_[Sufijo].tif (tipo: ráster; unidades: mm): Mapa de valores de flujo rápido (FR) anual

  • P_[Sufijo].tif (tipo: ráster; unidades: mm/año): La precipitación total en todos los meses de este píxel

  • Vri_[Sufijo].tif (tipo: ráster; unidades: mm): Mapa de los valores de recarga (contribución, positiva o negativa), a la recarga total

  • aggregated_results_swy_[Sufijo].shp: Tabla con los valores biofísicos de cada cuenca hidrográfica, con los campos siguientes:

    • qb (unidades: mm, pero deben interpretarse como valores relativos, no absolutos): Valor medio de la recarga local en la cuenca

    • vri_sum (units: mm): contribución total de recarga, (positiva o negativa) dentro de la cuenca. La suma de los valores de los píxeles de Vri_[Sufijo].tif dentro de la cuenca.

  • [Workspace]\intermediate_outputs folder:

  • aet_[Sufijo].tif (type: raster; units: mm): Mapa de evapotranspiración real (ETR)

  • qf_1_[Sufijo].tif…qf_12_[Suffix].tif (tipo: ráster; unidades: mm): Mapas de flujo rápido mensual (1 = enero… 12 = diciembre)

  • stream_[Sufijo.tif (tipo: ráster): Red de corrientes generada a partir del MDE de input y del Umbral de Acumulación de Flujo. Los valores de 1 representan corrientes, los valores de 0 son píxeles que no son corrientes.

Apéndice 1: Fuentes de datos y orientación para la selección de parámetros

Precipitation

Evapotranspiration

Digital Elevation Model

Land Use/Land Cover

Soil Groups

Watersheds

Curve Number

Kc

Rain Events

Threshold Flow Accumulation

Zonas climáticas

Los datos de las zonas climáticas están disponibles en el sitio de clasificación climática de Köppen-Geiger.

alpha_m

Por defecto: 1/12. Ver Apéndice 2

beta_i

Por defecto: 1. Ver Apéndice 2

gamma

Por defecto: 1. Ver Apéndice 2



Apéndice 2: Parámetros Alfa, beta y gamma - definición y valores alternativos

\(\alpha\) (alfa) y \(\beta_{i}\) (beta) representan la fracción anual de los píxeles de la ladera que está disponible para un píxel pendiente abajo para la evapotranspiración en un mes determinado. Se espera que el producto \(\alpha \times \beta_{i}\) sea <1, ya que parte del agua pendiente arriba puede no estar disponible, ya sea cuando sigue trayectorias de flujo profundas o cuando el calendario de la oferta y la demanda (de evapotranspiración) no está sincronizado.

\(\alpha\) es una función de la estacionalidad de las precipitaciones: la recarga de un mes determinado puede ser utilizada por las zonas de bajada durante los meses posteriores, dependiendo de los tiempos de recorrido del subsuelo. En la parametrización por defecto, su valor se establece en 1/12, asumiendo que el suelo amortigua la liberación de agua y que la contribución mensual es exactamente un 12th de la contribución anual. Un supuesto alternativo es fijar los valores a los valores de precipitación mensual antecedente, en relación con la precipitación total: Pm-1/Panual

\(beta_{i}\) es una función de la topografía local y de los suelos: para una cantidad dada de recarga de ladera, la cantidad de agua utilizada por un píxel es una función de la capacidad de almacenamiento. También depende de las características de la zona pendiente arriba: el uso del subsidio pendiente arriba está condicionado por la forma y la superficie de la zona de contribución (es decir, la recarga del píxel justo por encima del píxel de interés tiene menos probabilidades de perderse que los píxeles mucho más alejados).

En la parametrización por defecto, \(\beta\) se establece en 1 para todos los píxeles. Una alternativa es establecer \(\beta_{i}\) como TI, el índice de humedad topográfica para un píxel, definido como \(ln(\frac{A}{text{tan}\beta}\)) (u otra formulación que incluya el tipo de suelo y la profundidad).

γ (gamma) representa la fracción de la recarga del píxel que está disponible para los píxeles pendiente abajo. Es una función de las propiedades del suelo y posiblemente de la topografía. En la parametrización por defecto, γ es constante en todo el paisaje y juega un papel similar al de \(\alpha\).

En la práctica

Las opciones anteriores se ofrecen principalmente con fines de investigación. En la práctica, sugerimos que para los climas altamente estacionales, alpha debería ser establecida a los valores de precipitación mensual precedentes, en relación con la precipitación total: Pm-1/Panual .

A continuación, ofrecemos dos opciones para abordar la incertidumbre en torno a los valores de los parámetros:

  1. Verificación de la evapotranspiración real con las observaciones

El modelo da como resultado la evapotranspiración real en la escala temporal anual: ustedes pueden ajustar los parámetros para satisfacer la evapotranspiración real observada (por ejemplo, a partir de MODIS, https://www.ntsg.umt.edu/project/modis/mod16.php). En lo que sigue, «_mod» significa ETR (AET) modelizada, «_obs» significa ETR (AET) observada.

  • Si AET_mod > AET_obs, el modelo predice en exceso la evapotranspiración, lo que puede corregirse reduciendo los valores Kc, o reduciendo los valores gamma, y/o los valores beta (para que haya menos agua disponible para cada píxel).

  • Si AET_mod < AET_obs, el modelo no predice la evapotranspiración, lo que puede corregirse aumentando los valores de Kc (y aumentando los valores de gamma o beta si no están en su máximo de 1).

Si se dispone de valores mensuales de ETR, puede realizarse una calibración más fina cambiando el parámetro estacional alfa.

  1. Modelización de conjuntos

El modelo puede ejecutarse bajo diferentes supuestos y los resultados pueden compararse para estimar el efecto del error de los parámetros. Los rangos de los parámetros pueden determinarse a partir de las hipótesis sobre la proporción de subvención pendiente arriba disponible para un píxel determinado; pueden establecerse en los límites máximos (0 y 1) para obtener resultados preliminares.

Referencias

Allen, R.G., Pereira, L.S., Raes, D., Smith, M., 1998. Crop evapotranspiration - Guidelines for computing crop water requirements, FAO Irrigation and drainage paper 56. Roma, Italia.

Hamel, P. y Bryant, B. (2017). Uncertainty assessment in ecosystem services analyses: Seven challenges and practical responses. Ecosystem Services, Volume 24. https://doi.org/10.1016/j.ecoser.2016.12.008.

Hamel, P., Valencia, J., Schmitt, R., Shrestha, M., Piman, T., Sharp, R.P., Francesconi, W., Guswa, A.J., 2020. Modeling seasonal water yield for landscape management: Applications in Peru and Myanmar. Journal of Environmental Management 270, 110792.

NRCS-USDA, 2007. National Engineering Handbook. United States Department of Agriculture, https://www.nrcs.usda.gov/wps/portal/nrcs/detailfull/national/water/?cid=stelprdb1043063.