The tutorial is prepared to use the J2000 hydrological model for hydrological system analysis of a river catchment. A test catchment and dataset of the Dudh Kosi river basin has been provided along with the tutorial. The Dudh Kosi river basin was used for the hydrological system analysis by using the J2000 hydrological model as a part of the PhD research (Nepal, 2012). The information provided here is largely based on this study. PhD Thesis. The motivation, objectives and methodological approach and the rational of using the J2000 model in the Dudh Kosi river basin are also presented. The users can use the test data to get familiar with the model application. At the same time, users can prepare their own dataset to understand the hydrological system dynamics of any river basin by following this tutorial.

Índice 1 Who can use the tutorial 2 Description of the test dataset 2.1 Motivation 2.2 Study area 2.3 Objectives and methods 3 O modelo J2000 4 Preparação do Conjunto de Dados 4.1 Arquivos de parâmetros de modelos 4.1.1 Arquivo de parâmetros do solo 4.1.2 Arquivo de parâmetros de cobertura de terra 4.1.3 Arquivo de parâmetros hidrogeológico 4.1.4 Arquivos de parâmetros de HRUs e trechos fluviais 4.2 Dados de entrada meteorológicos 4.3 Área de trabalho para os dados de entrada 4.3.1 Pastas e arquivos 4.3.2 Modelo xml 5 Exibir resultados do modelo e de saída 6 Módulos dentro do modelo hidrológico J2000 6.1 Cálculo da Evapotranspiração 6.2 Módulo de distribuição de Precipitação 6.3 Módulo de Intercepção 6.4 Módulo de Neve 6.5 Módulo de água do solo 6.6 Módulo de Águas Subterrâneas 7 Calibração do modelo 8 Criação de um novo modelo 9 Fórum

Who can use the tutorial

The tutorial is prepared in such a way that the J2000 hydrological model can be used independently without any techtical support from model developers. Therefore, it can be used by students, model developers and researchers for the hydrological system analysis of a catchment. The tutorial should be read in conjunction with other sub-tutorials which has been mentioned in different part of this tutorial. Additionally, the tutorial is supplied with test dataset of the Dudh Kosi river basin (Nepal, 2012) which users can use to get familiar with the different aspects of the J2000 model. Similarly, users can also create their own dataset of the catchment of interest to run the model.

Description of the test dataset

The tutorial is accompanied by the test dataset of the Dudh Kosi river basin. This hydro-meteorological data were provided by the Department of Hydrology and Meteorology (DHM), Government of Nepal. The DHM has provided permission to use the data along with the tutorial. The users are expected to use the tutorial along with the test data to understand different aspects of the J2000 modelling system and also aim to prepare their own dataset to run the model.

Motivation

This is the motivation of the study.

Study area

This is the Study area, Dudh Kosi river basin.

Objectives and methods

Objectives and methods

O modelo J2000

O modelo J2000 é um modelo hidrológico distribuído e orientado em processos para simulações hidrológicas de bacias de meso e macro-escala. Ele está incorporado no Jena Adaptable Modelling System (JAMS), que é uma plataforma de software baseada em componentes para o desenvolvimento e aplicação de modelos ambientais (Kralisch e Krause 2006, Kralisch et.al. 2007). A simulação de diferentes processos hidrológicos é realizada em módulos de processo encapsulados, os quais são em grande parte independentes uns dos outros. Isto permite alterar, substituir ou adicionar módulos individuais ou processos sem se ter que reestruturar novamente o modelo desde o início. Com esta flexibilidade, um módulo de geleira é integrado como parte do estudo realizado por Nepal (2012) na região do Himalaia.

O aplicativo de modelagem representa os processos hidrológicos importantes de uma bacia hidrográfica. O layout principal do modelo hidrológico J2000 é fornecido na figura abaixo. O layout também inclui o módulo de uma geleira, que foi aplicado na região do Himalaia. O sistema de modelagem diferencia-se entre os quatro componentes de escoamento diferentes de acordo com sua origem específica. O componente com maior dinâmica temporal é o escoamento rápido e direto (RD1) (fluxo terrestre). Ele consiste no escoamento de áreas fechadas e escoamento superficial vindos do escoamento de acesso saturado de acesso de infiltração. O componente de escoamento lento direto (RD2) (também conhecido como Interfluxo 1), que corresponde ao escoamento hipodérmico lateral na zona do solo, reage ligeiramente mais lentamente do que o RD1. Dois outros componentes de escoamento de fluxo básicos podem ser distinguidos. O escoamento de base relativamente rápido (RG1) (também conhecido como Interfluxo 2) simula o escoamento a partir da parte superior de um aquífero, que é mais permeável devido às intempéries, em comparação à zona inferior do aquífero. O componente de escoamento de base lento (RG2), o qual pode ser visto como o fluxo dentro de fracturas de rochas sólidas ou de matriz em aquíferos homogéneos não consolidados.

A descrição pormenorizada dos sistemas de modelagem pode ser encontrada em muitas publicações. Algumas das publicações importantes são: (Krause, 2001,; Krause 2002,; Krause, 2010,; Nepal, 2012; Kralisch e Krause 2006,; Kralisch et.al. 2007). Algumas das publicações também podem ser acessadas ​​a partir deste link: http://jams.uni-jena.de/index.php?id=5582&L=2

Preparação do Conjunto de Dados

Arquivos de parâmetros de modelos

Os requisitos de dados para executar o modelo hidrológico J2000 são discutidos aqui em detalhe. Dois tipos de dados são necessários i) arquivos de parâmetros de modelo e ii) os dados de entrada meteorológicos. Os primeiros são preparados e quantificados dentro do ambiente GIS e conhecidos como arquivos de parâmetros de modelo. Os arquivos de parâmetros e seus valores ficam estáticos na aplicação de modelagem.

Os usuários têm que preparar todos os dados de entrada (ou seja, solo, cobertura do solo, geologia, DEM) em formato raster com uma dada resolução. Ao delinear HRUs, todos os dados de entrada têm de ser fornecidos na mesma resolução. A resolução do conjunto de dados controla principalmente o número de HRUs a serem formadas sem se perder a heterogeneidade de uma bacia. Portanto, a resolução dos dados de entrada depende da bacia a ser modelada. Por exemplo, se a bacia é pequena (por exemplo, 1000 km ²), uma resolução entre 30-90 será adequada dependendo da resolução do conjunto de dados disponível. Da mesma forma, para bacias hidrográficas de meso-escala (por exemplo, 4000 km ²), uma resolução entre 250-500 m será adequado. Além disso, uma bacia com topografia plana (baixo gradiente, por exemplo) precisa de dados de alta resolução para caracterizar as traços de uma bacia hidrográfica.

A descrição detalhada para a derivaçao dos arquivos de parâmetros são fornecidas a seguir:

Arquivo de parâmetros do solo

A informação detalhada necessária para um arquivo de parâmetro do solo é apresentada na tabela abaixo.

Parâmetro	Descrição
SID	tipo de solo ID
depth	profundidade do solo
kf_min	coeficiente de permeabilidade mínima
depth_min	profundidade do horizonte acima do horizonte com o menor coeficiente de permeabilidade
kf_max	coeficiente de permeabilidade máxima
cap_rise	variável de Boole, que permite (1) ou restringe (0) a ascensão capilar
aircap	capacidade de ar (equivalente ao reservatório de poros grandes (LPS))
fc_sum	capacidade de campo utilisável (equivalente ao reservatório de poros médios (MPS))
fc_1 ...22	capacidade de campo utilisável por decímetro da profundidade do perfil

O arquivo de parâmetros do solo é um dos importantes, os quais necessitam de uma série de informações, como mostrado na tabela acima, para produzir uma caracterização completa sobre a capacidade de retenção de água de diferentes tipos de solo. Para tal, a informação sobre a textura de tipos de solo de diferentes horizontes de solo é necessária. Uma descrição detalhada de como produzir um arquivo de parâmetro do solo pode ser encontrada aqui:

[Soil_parameterization [| Como preparar arquivo de parâmetros do solo]]

Arquivo de parâmetros de cobertura de terra

O arquivo de parâmetro do uso da terra requer informações sobre o uso da terra e da cobertura da terra de uma bacia hidrográfica que controla os diferentes aspectos da hidrologia. Tal informação pode ser derivada a partir da literatura na qual a informação espacial acerca do uso do solo e da cobertura da terra é fornecida. Alternativamente, tal informação pode ser estimada por meio de imagens de satélites e de classificação subsequente. O modelo hidrológico J2000 requer uma classificação maior do uso do solo e da cobertura da terra, o que afeta a dinâmica hidrológica.

Como preparar a terra arquivo de parâmetros de cobertura

Arquivo de parâmetros hidrogeológico

As informações necessárias para o arquivo de parâmetros hidrogeológico estão disponíveis abaixo:

• Hgeo.par

Parâmetro	Descrição
GID	hidro-geologia ID
RG1_max	capacidade máxima de armazenamento do reservatório de água subterrânea superior
RG2_max	capacidade máxima de armazenamento do reservatório de água subterrânea inferior
RG1_k	coeficiente de armazenamento do reservatório de água subterrânea superior
RG2_k	Coeficiente de armazenamento da parte inferior do reservatório de água subterrânea

A capacidade de armazenamento do reservatório de água subterrânea superior (RG1) e inferior (RG2) pode ser estimada analisando-se a informação geológica da área. A capacidade de armazenamento é normalmente controlada pela formação geológica, os tipos de rochas, origem e natureza das rochas e permeabilidade. Estes valores são expressos como volume de armazenamento máximo em mm/dia de cada tipo de reservatório. Os valores de coeficiente de armazenamento (RG1_k e RG2_k) são usadas como um co-eficiente de recessão geral de dois reservatórios. Estes são expressos como tempo de retenção em dias no reservatório específico. A recessão é ainda controlada por um parâmetro de calibração flexível dentro do modelo.

A descrição detalhada do parâmetro hidrogeológico está disponível aqui: Como derivar hidro-geológica arquivo de parâmetro

Arquivos de parâmetros de HRUs e trechos fluviais

As Unidades de Resposta Hidrológicas (HRUs - Hydrological Response Units) são as unidades de modelagem para o modelo hidrológico J2000. As HRUs são "unidades de modelo espaciais que são unidades estruturadas distribuídas e heterogêneas com um clima, uso da terra, solos, geologia comuns , controlando seus dinâmica hidrológica" (Flügel, 1995). As áreas que compreendem propriedades semelhantes, como a topografia (declividade, outros aspectos), uso da terra, solo e geologia, e se comportam de forma semelhante em sua resposta hidrológica, são mescladas para se desenvolver uma HRU. A variação da dinâmica do processo hidrológico dentro da HRU deve ser relativamente pequena em comparação com a dinâmica de uma HRU diferente (Flügel 1995).

O processo de delineamento de HRUs é descrito no seguinte tutorial. GRAMA-HRU tutorial. Os usuários precisam preparar o seguinte arquivo para o delineamento de uma HRU.

• Modelo de Elevação Digital (DEM)
• Solo
• Uso da terra
• Hydro-geologia

Todos estes dados devem ser fornecidos em um arquivo formato de dados *.tiff com a mesma resolução. O delineamento de processos de HRUs acaba por fornecer arquivos de parâmetros de HRU e trechos fluviais.

• Arquivo de parâmetro HRU (* hru.par)

Parâmetro	Descrição
ID	HRU ID
x	Avanço para leste (easting) do ponto centróide
y	Avanço para o norte (northing) do ponto centróide
elevation	elevação média
area	área
type	tipo de drenagem: a HRU drena na HRU (2), a HRU drena na parte do canal (3)
to_poly	ID da HRU subjacente
to_reach	ID da parte do canal adjacente
slope	declive
aspect	aspecto
flowlength	Comprimento do fluxo
soilID	ID da classe do solo
landuseID	ID da classe do uso da terra
HgeoID	ID da classe hidrogeológica

Um exemplo de arquivo de parâmetro HRU é fornecido abaixo.

O arquivo de parâmetro HRU armazena os atributos espaciais da bacia hidrográfica onde a informação sobre a elevação, área, aspecto, as coordenadas (x, y), tipo de uso da terra (landuseID), hidrogeologia (hgeoID) e solo (soilID) são armazenados para cada HRU. As HRUs são topologicamente ligadas ao encaminhamento lateral para simular os processos de transporte de água a partir de uma HRU para outra HRU e foram ainda ligadas a um afluente próximo para o roteamento de trechos fluviais. A coluna (to_poly) define a HRU que passa água para a próxima HRU.

A conexão entre o arquivo de parâmetro HRU e outros arquivos de parâmetros é resolvida dentro da anterior. Por exemplo, no arquivo de parâmetro HRU, o id HRU 1 tem todas as informações necessárias, conforme exigido na tabela acima, incluindo o uso da terra, o solo e o tipo de geologia, as quais pertencem à HRU1:

• Arquivo de parâmetro de trechos (* reach.par)

Parâmetro	Descrição
ID	ID da parte do canal
length	comprimento
to-reach	ID da parte subjacente do canal
slope	declive
rough	Valor de rugosidade de acordo com MANNING
width	Largura

O arquivo de parâmetro de trechos armazena as informações sobre as características de fluxo, bem como a relação entre as redes de fluxo para realizar o roteamento de trechos. Ele ainda contém informação sobre a estrutura da topologia de fluxo observando o ID para cada trecho fluvial em que transfere. O parâmetro de trechos é produzido em conjunto com o processo de delimitação da HRU e compreende também a informação requerida na tabela acima.

Com relação à figura do arquivo de parâmetro HRU acima, o ID HRU1 contribui com água diretamente para o REACH ID1, enquanto o HRU ID 16 contribui com água para HRU ID 5 que então contribui para o REACH ID 2. As interações entre os arquivos de parâmetro foram resolvidos por uma relação entre solo, uso da terra e descritores hidrogeológicos no arquivo de parâmetro HRU e respectivos descritores nos arquivos de outros parâmetros.

Dados de entrada meteorológicos

O modelo J2000 hidrológico requer os seguintes dados de entrada para a inicialização do modelo:

Nome	Descrição	Unidade
Ahum.dat	| g / cm 3
Orun.dat	passagem do fluxo medida no segundo turno	m3/s
Rain.dat	quantidade medida de precipitação	milímetros
Rhum.dat	umidade relativa	%
Sunh.dat	duração do sol	h
Tmax.dat	temperatura máxima	°C
Tmean.dat	temperatura média do ar	°C
Tmin.dat	temperatura mínima	°C
Wind.dat	velocidade do vento	m/s

O sistema de modelagem J2000 usa ponderações de distâncias inversas (IDW - Inverse Distance Weightings) com o método de correção de altitude para a regionalização dos dados climáticos de entrada. A descrição detalhada da abordagem de regionalização pode ser acessada aqui:

[Regionalization approach of the J2000]

Nem todos os dados, como mostrados na figura acima, podem estar disponíveis em algumas bacias. Normalmente, os dados de temperatura e precipitação estarão disponíveis. Se há apenas poucas estações (menos de três) para alguns parâmetros, a IDW não funcionará corretamente. Neste caso, o mesmo valor de entrada é aplicado à bacia inteira. Para algumas variáveis ​​particulares, ex. a temperatura, esta abordagem poderia trazer grandes quantidades de erros / incertezas. Em tais casos, sugere-se a abordagem da regionalização com base em um gradiente térmico para a temperatura. Os detalhes desta abordagem são fornecidos em Nepal, 2012.

As horas com luz solar, a umidade relativa e o vento também não estão freqüentemente disponíveis em algumas bacias. Estes valores são utilizados para a estimativa da evapotranspiração, ao usar-se a abordagem de Penman-Monteith. As horas de sol e velocidade do vento podem ser consideradas como suficientes a partir de uma estação, caso nenhum dado de outra estação estiver disponível. Em tais casos, o mesmo valor de estação é aplicado a uma bacia inteira. Esse valor único de estação para a umidade relativa também traz alguns erros consigo, ao calcular a umidade relativa utilizando a umidade absoluta e a temperatura. No sistema de modelagem J2000, não é recomendada uma regionalização direta dos valores de humidade relativa. Os detalhes são fornecidos no cálculo da evapotranspiração sub-tutorial calculation of evapotranspiration.

No caso destes dados (rhum, sunh, wind) não estiverem disponíveis, a abordagem Pennmann-Monteith não pode ser utilizada. Em seu lugar uma abordagem mais empírica com base na temperatura, tais como a Hargreaves, poderá ser usada.

Uma amostra dos dados de entrada do arquivo de precipitação (rain.dat) pode ser encontrado abaixo:

Os dados de entrada devem ser salvos com a extensão dat (exemplo: rain.dat). Os dados em formato Excel podem ser salvos como 'Texto (guia delineado) (*.txt)' e alterando a extensão de *.txt para *.dat. No final de cada conjunto de dados, a coluna de dados deve ser terminada com # data.dat . Para mais detalhes, baixe o arquivo de amostra de dados:

Cada ficheiro de dados tem a seguinte estrutura (demonstrado aqui para os exemplos de "precipitação"):

Linha	Descrição
#rain.dat chuvas
@dataValueAttribs
rain 0 9999 mm	nome da série de dados, o menor valor possível, o maior valor possível, a unidade
@dataSetAttribs
missingDataVal -9999	Valor para marcar valores de dados em falta
DataStart 1979/01/01 07:30	data e hora do primeiro valor de dados
DataEnd 2000/12/31 07:30	data e hora do último valor de dados
tres D	resolução temporal dos dados (aqui: dias)
@statAttribVal
name Stat1 stat14	nomes das estações pluviométricas
ID 1574... 1309	ID numérica das estações pluviométricas (ID)
elevation 525... 321	estação de elevação 1 ... estação de elevação 14
X 4402310... 4406282	estação de avanço ao leste 1 ... estação de avanço ao leste 14
Y 5620906... 5644937	estação de avanço ao norte 1 ... estação de avanço ao leste 14
dataColumn 1... 14	número da coluna em particular na parte de dados
@dataVal	parte inicial de dados
1979/01/01 07:30 0,8... 0,3	data, hora, valor, estação 1 ... estação 14
...
2000/10/17 07:30 0,1..0.1	data, hora, valor, estação 1 ... estação 14
#end of rain.dat	fim da parte de dados

Os arquivos de exemplo dos dados de entrada podem ser baixados a partir do diretório de área trabalho aqui.

Área de trabalho para os dados de entrada

Os dados de entrada do modelo hidrológico J2000 tem de ser fornecidos em uma pasta específica considerada como um "diretório de área de trabalho". Este diretório contém todos os dados de entrada necessários para se executar o modelo, bem como os arquivos de saída do modelo.

Diretório da área de trabalho do modelo J2000 DudhKosi

Módulo de extensão de geleira para o modelo J2000

Nota Importante:

O J2000Himalaya.zip contém um arquivo J2000Himalaya.jar que é uma extensão do módulo geleira ao modelo hidrológico J2000 padrão. Portanto, este arquivo jar tem que ser copiado na pasta lib (\Arquivos de Programas\jams\lib) [se sua bacia tiver geleiras]. A pasta lib já contém o arquivo J2K.jar, quando usuários fazem o download do software JAMS, juntamente com o conjunto de dados de teste da bacia Gelberg. Eles também podem manter o arquivo J2000Himalaya.jar em locais diferentes, mas o caminho tem de ser definido quando o modelo for executado pela primeira vez usando-se os seguintes passos:

[JAMS Launcher (ou JAMS Builder) - >> Editar - >> Editar Preferências]. Uma nova janela "Preferências JAMS" aparecerá. Os usuários precisam localizar o local do arquivo jar * clicando no sinal +).

Pastas e arquivos

O diretório da área de trabalho tem três pastas principais: entrada, saída e parâmetros. A pasta de entrada tem todos os dados de entrada necessários para executar o modelo. A pasta “parâmetros” tem arquivos de parâmetros (hru.par, hgeo.par, landuse.par, soil.par e reach.par). A pasta de saída contém arquivos de saída de diferentes variáveis ​​depois que o modelo tiver sido executado com sucesso. Uma área de trabalho de amostra do conjunto de dados de teste é provida em anexo, o que dá uma idéia de como organizar o diretório de área de trabalho para o modelo hidrológico J2000.

Pasta: entrada (input)

A pasta de entrada tem 12 arquivos com extensão xml de todos os dados de entrada. Copie e cole esses arquivos já que eles são o conector para a entrada de dados reais que estão localizados dentro de uma pasta "local".

Subpasta: local

A pasta dentro da pasta de entrada contém dados de entrada para oito variáveis ​​(rain.dat, rhum.dat, sunh.dat, tmax.dat, tmean.dat, tmin.dat, wind.dat). ahum.dat é criada quando o modelo é executado pela primeira vez, utilizando os dados existentes.

Subpasta: gis

Algumas camadas GIS podem ser colocadas aqui para visualizar a distribuição espacial de algumas variáveis ​​de saída (por exemplo, a distribuição espacial da precipitação em uma bacia (2D e 3D). Para isso, os usuários precisam utilizar um arquivo DEM de uma bacia hidrográfica (formato de dados: *.asc). A resolução do DEM deve ser semelhante à entrada de DEM para o processo de delimitação de HRUs. Os usuários precisam copiar o arquivo styles.sld, que é necessário para se exibir um mapa. Além disso, as HRUs, córregos e arquivos de dados de estação (*.shp) podem ser colocados em uma pasta separada para se visualizar as variáveis ​​em um componente de mapa. Os nomes destes arquivos e pastas devem ser definidos em um arquivo de modelo xml [model xml file].

'subpasta: despejo (dump)'

Crie uma pasta com o nome "dump (despejo)", que será usado para armazenar os arquivos temporários.

Pasta: saída (output)

A pasta de saída tem dois arquivos xml (HRULoop.xml e TimeLoop.xml) e uma corrente de pastas. Estes aquivos *.xml definem as variáveis ​​para as quais os produtos de saída são criados. Da mesma forma, os dados de saída são colocados dentro da pasta atual (nomes de arquivo: TimeLoop.dat e HRULoop.dat). A relevância destes dados de saída é discutida no sub-tutorial "Produto do Modelo abaixo. [Subseção: visor numérico]

Pasta: parâmetros (parameter)

Essa pasta contém arquivos de parâmetros (hgeo.par, hru.par, landuse.par, reach.par, soils.par). Lembre-se que estes nomes devem ser os mesmos no arquivo de modelo xml.

Pastas: explorador e tmp (explorer and tmp)

A outra pasta dentro da área de trabalho é o explorador e o tmp que são utilizados para armazenar alguns arquivos temporários gerados durante a execução do modelo.

Modelo xml

O diretório da área de trabalho também contém um modelo de arquivo xml. Os arquivos de modelo podem ser lidos como *.xml ou *.jams. Estes arquivos de modelo são fornecidos em um conjunto de dados teste.

Um exemplo de alguns arquivos de modelo xml é fornecido em anexo. Descompacte o arquivo para usá-lo.

Arquivo de modelo xml da bacia do Gelberg: Este é um arquivo de modelo xml padrão do J2000 fornecido juntamente aos dados de teste do Gelberg. A bacia do Gelberg na Alemanha tem dados de entrada suficientes para utilizar o método IDW para a regionalização.

Arquivo de modelo xml da bacia do Gelberg: Este é um arquivo de modelo xml da bacia do rio Kosi Dudh na região do Himalaia. A região tem geleiras e apenas uma estação de temperatura. Portanto, módulo de taxa de lapso da temperatura é usado para a regionalização da mesma.

J2K com o módulo Hargreaves: Este é um arquivo de modelo xml para regiões de dados escassos (apenas temperatura e precipitação) e a evapotranspiração potencial é calculada usando o método Hargreaves Salami.

Através da aplicação de diferentes módulos, o requisito de dados para o modelo é alterado. Em tal situação, diferentes módulos são desativados ou removidos no arquivo xml. Por exemplo, para o módulo de Hargreaves Salami, o vento, as horas com luz solar, e a umidade relativa não são necessários. Portanto, o leitor de dados e regionalização destes parâmetros são desativados no J2k_Hargreaves xml. Se for necessário usar o módulo Hargreaves Salami, é aconselhável comparar os arquivos xml ou requisitos de dados e dados xml do Kosi Dudh ou do Gelberg.

O modelo xml contém a estrutura lógica da estrutura do modelo e os módulos utilizados no mesmo. Ele é organizado de forma sistemática, de modo que o produto de saída de um módulo é fornecido como entrada para a próxima (exemplo: Snowmelt - Fusão da neve – ((o output do módulo de neve serve de entrada para o módulo de água do solo). Este arquivo também contém informações sobre a exibição de diferentes variáveis ​​e saídas no âmbito JAMS. O modelo xml pode ser visto usando o JAMS Builder para entender o componente diferente de um modelo específico.

Pode-se seguir este tutorial para se familiarizar com o [Http://jams.uni-jena.de/ilmswiki/index.php/Tutorial_Advanced_Users # JAMS Builder]. Frasco Um exemplo de modelo do Kosi Dudh é fornecido no JAMS Builder na figura abaixo.

A janela da esquerda mostra a localização dos arquivos de modelo de código fonte, que são necessários para executar um modelo. Todos os códigos-fonte do modelo estão dentro do arquivo *.jar. Para o modelo J2000, é utilizado um arquivo J2000.jar. A janela do meio fornece informações sobre os módulos utilizados no arquivo de modelo xml do modelo Kosi Dudh. Estes são fornecidos em uma estrutura lógica em que os resultados de um módulo são fornecidos como dados de entrada para o próximo. Uma descrição detalhada destes módulos diferentes será fornecida na seção subsequente, ao passo que um exemplo do módulo de regionalização de temperatura máxima (TmaxLapseRate) é mostrada no JAMS Builder abaixo. O módulo utiliza dados de temperatura de uma única estação para regionalizar a temperatura máxima em uma bacia. Ao clicar-se no TmaxLapseRate na coluna do meio sob Regionalização, a informação detalhada do módulo é exibida nas janelas à direita, como mostrado na figura a seguir.

Todas as variáveis ​​utilizadas nos módulos são fornecidos na coluna "Nome", como mostrado na figura acima. A coluna 'Tipo' descreve as características das variáveis ​​em termos da informação (dados, qualidade), as quais essas variáveis armazenam. A coluna R / W determina a natureza dessas entradas e bem como dos resultados. O 'statElev' é a elevação de uma estação de temperatura, o que é indicado por R. Isto significa que a informação é a entrada para o módulo a partir de um módulo anterior e denotados por 'Read'. O 'W' denota escrever/gravar, o que significa o novo valor de saída de este módulo. Os parâmetros de calibração, se for o caso, são fornecidos na coluna "Valor". Esta informação pode ser alterada no JAMS Builder instantaneamente. Ao clicar-se na variável, a informação sobre a 'Configuração de Atributo' será preenchida. Esta informação pode ser alterada. Clique em 'set' para guardar as informações.

O modelo pode ser executado a partir do JAMS Builder, clicando-se no botão "Executar Modelo[1]" ou "Executar Modelo a partir do JAMS Launcher [2]", como indicado pela caixa vermelha na figura abaixo. Ao clicar-se na caixa "1", o modelo será inicializado diretamente, enquanto um clique na caixa "2" inicializará o JAMS Launcher. Este último irá fornecer opções para alterar os parâmetros do modelo (como valores de parâmetro, período de tempo) etc

Uma descrição mais detalhada do arquivo de modelo xml em relação aos diferentes módulos e as variáveis ​​utilizadas nos códigos de fonte do modelo podem ser encontrados em Krause (2011). Este documento pode ser gerado a partir do JAMS Builder instantaneamente. Clique em 'Modelo' no menu superior e em 'Gerar documentação do modelo". Esta será disponibilizada em formato pdf para download.

O modelo xml também pode ser visto e editado através de um editor de textos (como o PSPad) como mostrado abaixo para o 'módulo de regionalização de temperatura máxima'.

---

 <component class="org.unijena.j2k.regionalisation.TemperatureLapseRate1" name="TmaxLapseRate">

• Esta linha define o local do módulo TemperatureLapseRate1 no arquivo de biblioteca de modelos J2000.jar *.

 <var name="lapseRateWinter"value="0.6"/>  <var name="lapseRateSummer"value="0.55"/>

• O valor de lapseRateWinter e lapserateSummer é o parâmetro de calibragem, que é uma taxa de lapso de mudança de temperatura por 100 metros.

 <var attribute="elevation" context="HRULoop" name="entityElev"/>

• O atributo elevation define a elevação de uma HRU que é a variável de entrada para o módulo TemperatureLapseRate1. O modelo lê a elevação de cada HRU a partir do ficheiro de parâmetros HRU como explicado anteriormente.

 <var attribute="time" context="J2K" name="time"/>

• time define a resolução temporal do modelo (por exemplo, diariamente)

 <var attribute="tmax" context="HRULoop" name="outputValue"/>

• Tmax é a temperatura máxima como um valor de saída do módulo, que calcula uma temperatura máxima de cada HRU utilizando as variáveis ​​de entrada dentro do módulo.

 <var attribute="elevationTmax" context="J2K" name="statElev"/>

• ElevationTmax é a variável de entrada da estação de elevação de temperatura máxima. O modelo lê a partir da estação de elevação de temperatura em um arquivo de entrada (por exemplo, tmax.dat).

 <var attribute="dataArrayTmax" context="J2K" name="inputValue"/>

• DataArrayTmax é a variável de entrada da temperatura máxima em uma determinada data.

 <var attribute="tmaxOrder" context="HRULoop" name="statOrder"/>

• TmaxOrder define a ordem da estação em caso de mais de uma estação de temperatura estiver disponível com base na distância de uma HRU com esse tipo de estação. Nesse caso, o modelo reconhece a estação próxima da HRU em particular.

A ordem lógica das variáveis ​​utilizadas no módulo é muito importante no modelo de arquivo xml. Por exemplo, cada uma das variáveis ​​de entrada deve ser definida anteriormente, antes de ser usada. Por exemplo, a variável de entrada "dataArrayTmax" é definida anteriormente em um módulo chamado 'TmaxDataReader'. O módulo de leitura de todos os dados diários de temperatura máxima em um formato que o modelo pode reconhecer. Da mesma forma, a variável de saída "tmax" (temperatura máxima) é, então, usada numa seção posterior. Por exemplo: o valor tmax (temperatura máxima) de cada HRU é usado para calcular o seu degelo.

Os parâmetros de calibração podem ser exibidos na tela do JAMS antes de se executar o modelo. Para isso, pode-se utilizar o GUI Builder do JAMS Launcher (figura abaixo). Aqui foi mostrado o exemplo taxa de lapso de temperatura. Essa taxa pode ser mantida sob a regionalização de grupo.

1. Primeiramente, clique no GUI Builder e escolha "regionalização" de grupo e, em seguida, clique em Adicionar propriedades. Uma nova janela, Editor de parâmetros de Modelo aparecerá, dentro da qual a informação do componente de modelo pode ser alterada.

2. Por exemplo, na classe 'Componente', escolha TmaxLapseRate. Em Variável/atributo, escolha lapseRateSummer. Da mesma forma, o nome e descrição do parâmetro também têm de ser preenchidos, em conjunto com o seu limite inferior e superior. Os usuários podem escolher o parâmetro de calibração entre a faixa de fronteira. Eles também podem inserir informações extras sob Texto de ajuda para fornecer informações mais detalhadas sobre o parâmetro. Quando toda a informação estiver preenchida, um clique em "OK" encorporará o parâmetro na modelagem. Do mesmo modo, no passo seguinte, os utilizadores podem escolher o 'lapseRateWinter' para o 'TmaxLapseRate'. Isto se deve ao fato de que a temperatura máxima é regionalizada com duas taxas diferentes de lapso de verão e de inverno. Pode ser o caso de se necessitar preencher as outras informações, como mostrado na figura a seguir:

O mesmo processo pode ser repetido com o TmeanLapseRate, o TminLapseRate e por períodos de verão e de inverno.

A figura abaixo mostra os parâmetros de calibração para a taxa de lapso, que são exibidos na tela do JAMS. O usuário pode alterar o valor do JAMS Launcher antes de se executar o modelo.

Exibir resultados do modelo e de saída

O arquivo xml compreende componentes para exibir os resultados do modelo de certas variáveis. Os resultados de saída do modelo podem ser vistos na janela do JAMS após a execução do modelo. A distribuição espacial de certas variáveis ​​de saída (tais como evapotranspiração, precipitação) pode ser exibida no final da execução do modelo, incluindo o mapa 3D. Além disso, o modelo também fornece os resultados de eficiência de diferentes avaliação estatística no final da execução do modelo.

A descrição da saída dos resultados do modelo são fornecidos no seguinte sub-tutorial : [Output of model results]

Módulos dentro do modelo hidrológico J2000

Esta seção descreve os diferentes módulos dentro do modelo hidrológico J2000. Os parâmetros de calibração aplicáveis ​​a cada um dos módulos estão também descritos com foco na sua influência sobre a vazão. A seguir estão os módulos importantes no modelo hidrológico J2000.

• Módulo de distribuição de Precipitação
• Módulo de Intercepção
• Módulo de neve
• Módulo de geleira
• Módulo de Intercepção
• Módulo do solo
• Módulo de águas subterrâneas
• Roteamento lateral
• Alcance de roteamento

Cálculo da Evapotranspiração

A descrição detalhada de cálculo da evapotranspiração usando a abordagem Pennmann-Monteith (Allen, 1998) também é oferecida neste sub-tutorial.

[Calculation of evapotranspiration]

Módulo de distribuição de Precipitação

Parâmetros de calibração

Parâmetros	Descrição	Escala Global	Para o modelo do Dudh Kosi
Trans	temperatura limiar	0 to - 5	2
Trs	temperatura básica para chuva e neve	-5 to +5	0

Estes parâmetros são considerados como não-flexíveis e não necessariamente colocados na tela do JAMS como parâmetros ajustáveis.

No sistema de modelagem J2000, a precipitação é primeiramente distribuída entre chuva e neve, dependendo da temperatura do ar. Dois parâmetros de calibração (Trans, eTrs) são usados ​​onde Trs é a temperatura base e Trans é um intervalo de temperatura (limite superior e inferior) acima e abaixo da temperatura de base. A fim de se determinar a quantidade de neve e chuva, presume-se que uma precipitação abaixo de certas temperaturas limiares resulta em precipitação totalmente de neve, ao passo que acima de um outro limiar ela será interiamente composta de chuva. No gama (Trans) entre as referidas temperaturas limiares, a precipitação ocorre de forma mista.

Relevância na modelagem

Ao inserir-se os valores Trs abaixo de zero (por exemplo, -2), tem-se mais precipitação sob a forma de "chuva" do que de "neve".

A descrição detalhada deste módulo, juntamente com o algoritmo, tal como definido no código-fonte é fornecida [Precipitation distribution module]

Módulo de Intercepção

"Parâmetros de calibração '

Parâmetros (unidades)	Descrição	Escala Global	Para o modelo do Dudh Kosi
α_rain (mm)	capacidade de armazenamento (m²) de cobertura do solo em particular para chuva em mm	0 to +5	1.0
α_snow (mm)	capacidade de armazenamento (m²) de cobertura do solo em particular para neve em mm	0 to +5	1.28

Os parâmetros de calibração do módulo de intercepção no JAMS Builder é fornecido na figura abaixo:

A intercepção é um processo durante o qual a precipitação é armazenada em folhas, e outras superfícies abertas de vegetação. Este processo é identificado como componentes importantes de um ciclo hidrológico que podem afetar os componentes de balanço hídrico. Copa e intercepção de resíduos são consideradas perdas para o sistema, já que qualquer chuva interceptada por um desses componentes será posteriormente reduzida pelo processo de evapotranspiração. O módulo de intercepção utiliza uma abordagem simples de armazenamento de acordo com Dickinson (1984), que calcula a capacidade máxima de armazenamento de intercepção com base no Índice de área foliar (LAI), do tipo particular de cobertura do solo. O modelo recebe a informação de LAI de diferentes estações do arquivo de parâmetros da cobertura da terra. Quando o armazenamento máximo é atingido, o excedente é passado como interceptado ao módulo do solo.

Os parâmetros, como mostrado na figura acima tem um valor diferente, dependendo do tipo da precipitação interceptada (chuva ou neve). Isso é dessa forma devido a capacidade de interceptação máxima de neve ser ​​visivelmente maior do que a de precipitação líquida.

Relevância na modelagem

Ao manter-se o alto valor dos parâmetros (a_snow e a_rain), têm-se um armazenamento maior de água em superfícies de folhas, o que leva a uma evapotranspiração maior e menos água disponível para o próximo módulo (ou seja, módulo do solo), escoamento e vice-versa.

A descrição detalhada deste módulo, juntamente com o algoritmo, tal como definido no código-fonte pode ser encontrada aqui [Interception module]

Módulo de Neve

"Parâmetros de calibração '

Parâmetros	Descrição	Escala Global	Para o modelo do Dudh Kosi
SnowCritDens	densidade crítica da neve	0-1	0,381
SnowColdContent	conteúdo frio da neve	0-1	0,0012
BaseTemp	temperatura limite para degelo	-5 a 5	0
T_factor	fator de derretimento no caso de calor sensível	0-5	2,84
R_factor	fator de derretimento no caso de precipitação líquida	0-5	0,21
G_factor	fatorde derretimento pelo fluxo de calor no solo	0-5	3,73

Os parâmetros de calibração do módulo de neve no JAMS Builder é fornecido na figura abaixo:

O módulo de neve J2000 é realizado de acordo com o método descrito em Knauf (1976) incluindo algumas modificações e extensões. O módulo calcula variáveis ​​de estado diferentes do bloco da neve, tais como sua profundidade, densidade e equivalente de água, o que pode ser usado para o processo de validação orientada a múltiplas respostas. O módulo descreve principalmente as fases de fusão, acumulação e afundamento do bloco de neve.

Se a temperatura de fusão for superior ao valor limite de temperatura (baseTemp), o processo de degelo começa. O parâmetro snowColdContent leva em conta o teor de frio do bloco de neve. A temperatura desse bloco tem que ser cerca de cerca 0 para iniciar o processo de degelo. Temperaturas do ar negativas são acumuladas e diminúem apenas com temperaturas positivas e para realizar o processo de degelo.

A energia de fusão necessária para o cálculo de derretimento da neve potencial é fornecido sob a forma de calor sensível à temperatura do ar ( t_factor), a entrada de energia a partir da precipitação na forma de chuva ( r_factor) e de entrada devido a um fluxo de calor no solo (g_factor). Esta taxa de derretimento potencial é ainda modificada de acordo com a inclinação e aspecto da HRU. A cobertura de neve é ​​capaz de armazenar água líquida (tal como encontrada na forma de derretimento de neve potencial) nos seus poros, até um certo limite de densidade (CritDens). Esta capacidade de armazenamento é perdida quase que completamente e de forma irreversível ao atingir uma certa densidade crítica da neve acumulada (ou seja, a quantidade de água livre em relação ao equivalente da água da neve total, entre 40% e 45%) de acordo com Bertle (1966) e Herrmann (1976). O degelo da neve resultante é passado ao módulo do solo.

Relevância na modelagem

A baseTemp é um limite de temperatura para a fusão da neve. A fusão ocorre apenas se a temperatura do ar for superior a baseTemp. Ao se manter um valor alto, tem-se mais neve para armazenar e menos degelo ocorre e vice-versa. Este parâmetro é muito importante no caso da bacia ter capacidade de armazenamento e derretimento de neve sazonais. O melhor valor pode ser encontrado em torno de 0.

O SnowCritDens é um parâmetro importante de calibração, porque ele permite que a água do bloco de neve funda quando a densidade crítica for superior a este valor. Esta capacidade de armazenamento é perdida quase que completa e irreversivelmente quando uma certa quantidade de água líquida, proporcionalmente ao equivalente total de neve (cerca de 40%) é atingida. Mantendo-se esse valor de parâmetro baixo, ocasionará a liberação da água em estado líquido a partir do bloco de neve rapidamente (já que a densidade crítica da neve é ​​atingida mais rapidamente), com menos quantidade de neve e profundidade relacionada. O elevado valor deste parâmetro resulta em mais tempo para atingir a densidade crítica, a qual exigirá mais neve fresca para ocorrer.

O SnowColdContent contribui para alcançar o teor de frio de uma estreita camada de neve a zero para o início do processo de fusão. Um valor maior ajudará a levar o teor de frio de uma camada de neve a zero mais rapidamente do que um valor baixo.

O T_factor é um dos parâmetros mais importantes e sensíveis para a produção de escoamento em fusão a partir de uma camada de neve. Um valor maior produzirá maior degelo e vice-versa. Quanto maior o t_factor, juntamente com o valor r_factor e o g_factor irá fornecer energia de calor para derreter a neve e produzir o escoamento.

A descrição detalhada deste módulo, juntamente com o algoritmo, tal como definido no código-fonte pode ser acessada aqui: [Snow module]

=== Módulo de geleira ===

"Parâmetros de calibração '

Parâmetros	Descrição	Escala Global	Para o modelo do Dudh Kosi
MeltFactorIce	fator para o derretimento do gelo	0-5	2,5
Tbase	temperatura limite para derretimento	-5 a 5	-1
DebrisFactor	fator de detritos para derretimento do gelo	0-10	3
AlphaIce	fator de derretimento de radiação para o gelo	0-5	0,2
Kice	fator de derretimento por calor sensível	0-50	5
KSnow	fator de derretimento por precipitação líquida	0-50	5
Krain	fator de derretimento pelo fluxo de calor no solo	0-50	3

Os parâmetros de calibração do módulo de geleira no JAMS Builder é fornecido na figura abaixo. Os usuários podem escolher o método para calcular o derretimento do gelo da geleira entre um fator grau-dia simples [1] e um fator grau-dia aprimorado [2] na 'fórmula de derretimento'. O parâmetro de calibração na parte inferior 'ddfIce' é aplicável para o fator grau-dia simples [1].

A neve sazonal ou fresca nas geleiras é processada ​​usando-se o J2KProcessSnow como descrito no Módulo de neve na seção anterior. Quando o armazenamento de neve for ​​zero na HRU de geleira, o processo de derretimento do gelo da geleira começa a usar o fator grau-dia maior. A energia para a o derretimento do gelo de geleira é representada pelo fator grau-dia, que é a função geral de temperatura e radiação (meltFactorIce e alphaIce). A taxa de fusão é ainda adaptada para a inclinação e o aspecto da HRU específica. O modelo também segrega as geleiras de HRU cobertos por detritos ou não, utilizando para isso, a inclinação das HRUs. Se a inclinação for superior a 30 graus, a geleira da HRU é considerada como uma geleira não coberta por detritos. Na geleira de HRU coberta por detritos, o degelo é ainda controlado pelo parâmetro de calibração debrisFactor. O resultado da área de geleira é a neve derretida (de neve fresca), derretimento do gelo da geleira e escoamento da água da chuva (chuva nas geleiras). Todos estes três escoamentos de derretimento são encaminhados através das áreas glaciais usando-se três encaminhamentos diferentes para neve, gelo e chuva (kSnow, KIce e Krain). No final, o material fundido da geleira, que é o produto de três escoamentos diferentes (derretimento da neve, do gelo, e drenagem da chuva), é passado para o próximo trecho como escoamento superficial (RD1).

Relevância na modelagem:

MeltFactorIce é um dos parâmetros sensíveis do módulo glacial. Quanto maior o valor, mais elevado o derretimento de gelo glacial e vice-versa.

Tbase é um limite de temperatura para se derreter a neve. A fusão ocorre apenas se a temperatura do ar for superior a tbase. Mantendo-se o alto valor, mais neve para armazenar assim como menos degelo ocorrerá e vice-versa. O melhor valor pode ser encontrado em torno do valor 0.

DebrisFactor controla o derretimento do gelo. O valor mais elevado (por exemplo, 4) significa que o deretimento é reduzido em 40% nas geleiras cobertas por escombros .

AlphaIce é o fator de radiação para degelo. Quanto maior o valor, maior o derretimento do gelo da geleira.

Kice, kSnow e Krain são coeficientes de roteamento. Os valores menores representam o curto tempo residencial do escoamento de derretimento potencial. Assume-se que a chuva sobre a superfície glacial drena mais rápido do que derretimento de neve e gelo (snowmelt/icemelt).

A descrição pormenorizada do módulo de geleira, juntamente com o algoritmo, tal como definido no código de fonte do modelo, pode ser encontrada aqui: Módulo geleira # Glacier_module

Módulo de água do solo

Parâmetros de calibração

Parâmetros (unidades)	Descrição	Escala Global	Para o modelo do Dudh Kosi
soilMaxDPS	reservatório de depressão máxima	0-10	2
soilLinRed	co-eficiente de redução linear para AET	-5 a 5	-1
soilMaxInfSummer	infiltração máxima no verão	0-200	60
soilMaxInfWinter	infiltração máxima no inverno	0-200	75
soilMaxInfSnow	infiltração máxima em áreas de cobertura de neve	0-200	40
soilImpGT80	infiltração para áreas maiores que 80% de vedação	0-1	0,5
soilImpLT80	infiltração para áreas menores que 80% de vedação	0-1	0,5
soilDistMPSLPS	coeficiente de distribuição MPS-LPS	0-10	0,3
soilDiffMPSLPS	coeficiente de difusão MPS-LPS	0-10	0,5
soilOutLPS	coeficiente de vazão para LPS	0-10	0,3
soilLatVertLPS	coeficiente de distribuição lateral vertical	0-10	0,5
soilMaxPerc	Taxa de percolação máxima para água subterrânea	0-100	10
soilConcRD1Flood	coeficiente de recessão para a inundação	0-10	1,3
soilConcRD1Floodthreshold	valor limite para soilConcRD1Flood	0-500	300
soilConcRD1	coeficiente de recessão para escoamento superficial	0-10	2,8
soilConcRD2	coeficiente de recessão para interfluxo	0-10	3

Os parâmetros de calibração do módulo de água do solo no JAMS Builder são fornecidos na figura abaixo:

O módulo de solo é a parte mais complexa do modelo hidrológico J2000, o que reflete o papel central da zona de solo como um regulamento e sistema de distribuição. A entrada para o módulo de solo é o derretimento da neve e a precipitação na forma de chuva. O reservatório de poros médios (MPS) e o reservatório de poros grandes (LPS) representa a capacidade de retenção de água do solo. A água nos MPS representa a capacidade de campo, no qual a água é mantida contra a gravidade, mas pode ser subtraída de uma tensão ativa, por exemplo transpiração de plantas. As informações de entrada sobre o módulo do solo são primeiramente utilizadas para preencher os MPS e LPS, o que determina a umidade atual do solo. As condições de umidade do solo influenciam o processo de infiltração. Ao passo que a capacidade de infiltração de um solo em particular é calculada com base na umidade real do mesmo. Além disso, existem três diferentes parâmetros de infiltração (soilMaxInfSummer, soilMaxInfWinter e soilMaxInfsnow), que controlam a infiltração durante o inverno, o verão e a cobertura de neve. No caso das superfícies seladas, apenas determinada quantidade de água na superfície é capaz de se infiltrar, a qual é controlada por dois parâmetros (soilImpGT80 e soilImpLT80). A água que não é capaz de se infiltrar é armazenada na superfície da terra, como o reservatório de depressão, até um certo grau, tal como definido no parâmetro de calibração (soilMaxDPS) e qualquer excesso é tratado como escoamento superficial (RD1). A água infiltrada é então distribuída entre os MPS e LPS, o que é controlado por dois parâmetros de calibração (soilDistMPSLPS e soilDiffMPSLPS). A água disponível no MPS pode ser reduzida por evapotranspiração, para a qual a profundidade da raiz da respectiva tampa de terra do solo é importante. A água no LPS se distribui entre os componentes laterais e verticais (soilLatVertLPS). O fluxo lateral é responsável pela produção de interfluxo na zona insaturada (RD2), que pode ser controlado por soilOutLPS. O fluxo vertical (percolação) é fornecido para as águas subterrâneas da zona (zona saturada), na qual a taxa de percolação é controlada por soilMaxPerc. O excedente é passado aos LPS para serem liberadas como escoamento superficial (RD2). O escoamento superficial e Interfluxo 2 são controlados por coeficientes de retenção (soilConcRD1 e soilConcRD2).

No caso do escoamento superficial, o período de tempo de retenção pode ser diferente em períodos de fluxo elevados, devido a um comportamento não-linear de uma bacia. Para isso, são introduzidos novos parâmetros para poder-se representar o comportamento não-linear durante a estação das monções da região do Himalaia. Por conseguinte, um novo parâmetro (soilConcRD1Flood), como um coeficiente de recessão de escoamento superficial durante o período de pico de inundação, quando o volume do escoamento superficial ultrapassa o valor-limite no soilConcRD1Floodthreshold definido como um parâmetro de calibração.

Relevância na modelagem

SoilMaxDPS influencia a água armazenada nas áreas de depressão. Quanto maior o valor deste parâmetro, maior quantidade de água a ser armazenada como reservatório de depressão e menos água ficará disponível para escoamento superficial.

SoilLinRed reduz a quantidade de taxa de evapotranspiração. Quanto menor o valor maior a redução da evapotranspiração.

SoilMaxInfSummer,soilMaxInfWinter, e soilMaxInfSnow influenciam a taxa de infiltração máxima para zonas saturada (água do solo) e insaturada (água subterrânea).Um menor valor para estes parâmetros permitirá que apenas uma parte da chuva e neve derretidas entrem na zona insaturada (O valor 20 significa que apenas 20 mm de chuva e neve derretidas poderá passar através do módulo de água no solo por intervalo de tempo). Em tais casos, o escoamento superficial seria maior, já que a maior parte da entrada é drenada como escoamento superficial. Considera-se que o soilMaxInfSummer é ligeiramente mais baixo do que no inverno, porque o solo é mais saturado durante o período de verão chuvoso. O soilMaxInfSnow é considerado menor entre os três, por causa da condição do solo congelado no ambiente onde ocorre neve.

O parâmetro soilImpGT80 é ativado, se a impermeabilidade da cobertura da terra é alta, tal como definido na permeabilidade da cobertura da terra (sealedgrade) (tais como áreas urbanizadas) no arquivo de parâmetros de cobertura de terra. Os valores no sealedgrade atuam como um limite para a ativação dos parâmetros soilImpGT80 ou soilImpLT80. Um valor menor destes parâmetros indica que a quantidade mais baixa de entrada será capaz de se infiltrar, bem como que o resto fluxo escoerá superficialmente.

SoilDistMPSLPS e soilDiffMPSLPS são principalmente responsáveis ​​pela distribuição e difusão da água entre os MPS e LPS. São parâmetros menos sensíveis e têm um papel menor no interfluxo. Um menor valor destes parâmetros irá alocar menos entrada de fluxo para os reservatórios de poros médios MPS.

SoilOutLPS influencia a componente Interfluxo 2 (RD2). Um valor menor irá fazer com que mais água seja escoada a partir de um reservatório de poros grandes LPS e, assim, aumentando o componente RD2.

SoilLatVertLPS é um dos parâmetros mais sensíveis no módulo de água do solo. Ele distribui o fluxo de entrada (depois da infiltração) entre vertical (interfluxo 1) e percolação. Quanto maior o valor irá alocar maior quantidade de fluxo para o Interfluxo 2 (e menos fluxo a ser percolado para a água subterrânea). Um valor maior irá aumentar o interfluxo 2 componente e, ao mesmo tempo, a contribuição das águas subterrâneas (RG1 e RG2) será reduzida.

SoilMaxPerc controla a taxa de infiltração de água subterrânea por intervalo de tempo. Um valor mais elevado (por exemplo, 20) indica que no máximo 20 milímetros de água equivalente poderão passar para as águas subterrâneas. Um valor maior irá aumentar a contribuição das águas subterrâneas (RG1 e RG2) e ao mesmo tempo diminuir o componente RD2 em maior magnitude. Isto também diminuirá a RD1, mas em menor grau.

SoilConcRD1 e soilConcRD2 são coeficientes de recessão para o escoamento superficial (RD1) e o Interfluxo 1 (RD2) e são uns dos parâmetros sensíveis no módulo. O valor representa o tempo de retenção (por intervalo de tempo) para escoamento superficial. O valor mais elevado indica que o período de retenção é elevado e, portanto, menos água flui como escoamento superficial. Basicamente, o período de retenção do RD1 deve ser menor do que o do RD2.

O parâmetro soilConcRD1Flood é implementado no modelo hidrológico padrão J2000 para replicar o comportamento de escoamento, especialmente durante o período de alta inundação na região do Himalaia, dominada por monções. Este parâmetro é apenas ativado quando a entrada para o escoamento superficial atravessa o parâmetro soilConcRD1FloodThreshold, o qual é definido pelo usuário. Em princípio, o valor de soilConcRD1Flood deve ser inferior a soilConcRD1, já que o tempo de retenção do escoamento superficial é menor durante o tempo de forte inundação.

A descrição pormenorizada do módulo de água do solo, juntamente com o algoritmo, tal como definido no código fonte é fornecido no modelo Soil Water Module.

Módulo de Águas Subterrâneas

Parâmetros de calibração

Parâmetros (unidades)	Descrição	Escala Global	Para o modelo do Dudh Kosi
gwRG1RG2dist	coeficiente de distribuição RG1-RG2	0-10	2,1
gwRG1Fact	adaptação para fluxo de RG1	0-10	0,3
gwRG2Fact	adaptação para fluxo de RG2	0-10	0,5
gwCapRise	coeficiente de ascensão capilar	0-10	0,01

Os parâmetros de calibração do módulo de águas subterrâneas no JAMS Builder são fornecidos na figura abaixo:

O módulo de águas subterrâneas recebe entrada a partir da zona do solo não saturado (módulo de água no solo), num compartimento de dois armazenamentos de uma zona de águas subterrâneas ou seja, de uma zona de águas subterrâneas superior (RG1) e uma zona de águas subterrâneas inferior (RG2). A entrada está então entre as duas zonas, nas quais a sua distribuição é realizada pelo parâmetro de calibração gwRG1RG2dist. A descarga de água a partir das áreas de armazenamento superior e inferior (RG1 e RG2) é levada a cabo de acordo com a quantidade de armazenamento de corrente sob a forma de uma função linear, utilizando-se o coeficiente de retenção de armazenamento para dois reservatórios gwRG1Fact e gwRG2Fact. Existe também a possibilidade de que as águas subterrâneas sejam transferidas a partir da zona de água do solo através da ascensão capilar com o parâmetro gwCapRise.

"'Relevância na modelagem

GwRG1RG2dist distribui a água de entrada para RG1 e RG2. Um valor maior para esse parâmetro aumenta a proporção de água de entrada para a zona RG2.

GwRG1Fact e gwRG2Fact influenciam a saída das unidades de armazenamento RG1 e RG2. Os valores do parâmetro representam o tempo de retenção nessas unidades. Quanto maior o valor, maior será a condução a uma menor vazão e mais água permanece nelas.

GwCapRise influencia a distribuição de água entre a água do solo e o módulo de água do solo. Com um maior valor, uma maior quantidade de água fluirá a partir de águas subterrâneas para a zona de água do solo.

A descrição pormenorizada do módulo de água subterrânea juntamente com o algoritmo, tal como definido no código de fonte do modelo é fornecida em Groundwater Module.

=== Módulo de Roteamento ===

Parâmetros de calibração

Parâmetros	Descrição	Escala Global	Para o modelo do Dudh Kosi	-	FlowRouteTA	Parâmetro de calibração para adaptar a velocidade das ondas de fluxo	0-10	1,3	-

Os parâmetros de calibração do módulo de alcance de encaminhamento no JAMS Builder são fornecidos na figura abaixo:

O modelo hidrológico J2000 tem dois componentes de roteamento. O roteamento lateral (lateral routing) serve para simular os processos de fluxo lateral na área de captação de uma unidade de modelo (HRU) para o próximo até que a água finalmente chega a um trecho fluvial. O roteamento de trechos fluviais (reach routing) descreve processos de fluxo de um canal de fluxo, utilizando a abordagem comumente aplicada de onda cinemática e o cálculo de velocidade de acordo com Manning e Strickler (Krause, 2001).

Relevância na modelagem

FlowRouteTA influencia o tempo de execução de ondas de escoamento no canal de fluxo. Um valor mais alto aumenta a velocidade das ondas de escoamento e mais água flui do canal.

Uma descrição detalhada do módulo de encaminhamento, juntamente com o algoritmo, tal como definido no código de fonte do modelo pode ser encontrada em Routing module.

Calibração do modelo

A fim de aplicar modelos hidrológicos com sucesso, é necessário definir com precisão os parâmetros dos mesmos. A medição direta dos parâmetros não é possível na maioria das vezes, ou muito cara ou não há relação física clara. Por estas razões, os parâmetros são ajustados com um processo de tentativa e erro, na medida que os elementos (por exemplo, os escoamentos simulados) correspondem melhor aos valores medidos. Esta tarefa pode ser demorada e difícil se o modelo correspondente for complexo ou tiver um grande número de parâmetros.

O modelo J2000 fornece a plataforma para o processo de calibração offline e online. A calibração offline é realizada no âmbito do JAMS, enquanto na calibração online, os arquivos de modelo e parâmetros necessários são definidos através de ferramentas de calibração baseadas na web chamadas 'OPTAS'. Em seguida, a calibração é realizada no servidor da Universidade e os resultados podem ser baixados. Esse processo é eficiente e consome menos tempo, já que o cálculo é realizado em um servidor externo, dessa forma não mantando os computadores locais ocupados.

As informações detalhadas sobre a calibração de modelos são fornecidas em Model calibration

Nota Importante: Os parâmetros do modelo, os valores específicos incluindo os tipos de parâmetros do modelo Dudh Kosi foram explicados anteriormente em cada um dos módulos. Os valores destes parâmetros foram definidos por uma combinação de 'tentativa-e-erro' e utilizando métodos de otimização de parâmetros automática ou numérica. Além disso, as análises de sensibilidade e incerteza também foram realizadas na bacia do rio Kosi Dudh. A descrição destes métodos e processos estão descritos em Nepal (2012).

Criação de um novo modelo

Para configurar o modelo hidrológico J2000 para uma captação nova, requere-se dois passos importantes. Primeiro, prepare os arquivos do modelo de parâmetros e dados de entrada, como explicado nas seções anteriores Dataset_preparation # seções anteriores. Segundo, crie um novo modelo de arquivo XML que controla as variáveis ​​de entrada e de saída com base nos dados de entrada. O modelo de arquivo XML poderá ser diferente em aplicações de modelos diferentes, já que esse se baseará em dados de entrada e módulos diferentes usados ​​para calcular processos hidrológicos. (Por exemplo: a exigência de dados para estimar a evapotranspiração potencial pelo método Hargreave-Salami é diferente da abordagem Penmann-Monteith, e, portanto, a configurção do modelo xml será também diferente.)

Recomenda-se usar os arquivos de modelo XML existentes como base para a criação de um novo modelo, a qual foi explicada nas seções anteriores. Se a captação tem geleiras, pode-se usar o modelo xml da bacia hidrográfica do Kosi Dudh. Além de geleiras, você pode usar o modelo de arquivo xml da bacia de Gelberg que é baixada durante o download de software JAMS como um conjunto de dados de teste. Esses arquivos de modelo xml podem ainda ser alterados de acordo com suas necessidades de dados e preferências.

Além disso, é preciso levar em conta algumas características específicas da bacia nova no arquivo de modelo xml.

1. As coordenadas geográficas da área de estudo em UTM tem de ser definidas em componente 'CalcLatLong' para estimar o 'fator de correção de aspecto de inclinação', como mostrado na figura abaixo.

2. As coordenadas geográficas da área de estudo têm que ser definidas no componente 'Calcular radiação extraterrestre (ExtRad)', baseando-se na latitude e longitude. Isto pode ser feito através da edição de informações no JAMS Builder como mostrado na figura a seguir:

Alternativamente, a informação pode ser alterada diretamente no JAMS Launcher.

Por favor, não se esqueça de salvar as novas informações.

3. O sistema de modelagem J2000 é bastante flexível em termos de uso de diferentes componentes com base na disponibilidade de dados. No conjunto de dados de teste fornecido neste tutorial, para a regionalização da temperatura, os dados de uma estação são usados e regionalizados no verão e no inverno usando taxas de erro. Mas, se houver dados de mais de três estações, eles podem também usar ponderações de distância inversas (IDW). Para tal, é preciso substituir o módulo de taxa de lapso pelo módulo de regionalização IDW. Ambos os módulos são fornecidos nos dois XMLS diferentes acima.

Da mesma forma, para estimar a evapotranspiração usando Penmann-Monteith, são necessários muitos dados, tais como: umidade relativa, vento e horas de sol. Estes podem não estar disponíveis em muitas áreas, no qual caso, o módulo de alternativa para a evapotranspiração chamado 'Hargreaves' pode ser empregado no modelo.

É necessário ter cuidado para que, ao se alterar módulos nos arquivos XML padrão, os quais estão disponíveis aqui (Dudh Kosi e Gelberg), os requisitos dos dados também mudem. Por razão disso, alguns módulos podem não ser necessários no novo modelo criado. Por exemplo, se o método Hargreves-Salame for usado, o componente de leitura de dados e um componente de regionalização para a umidade relativa, vento e da hora do sol não serão necessários. Em tais casos, esses módulos precisam ser desativados.

4. A localização do diretório de área de trabalho e um arquivo de dados tem de ser definida após o carregamento do arquivo de modelo xml, como mostrado na figura abaixo.

O arquivo de modelo xml é aberto usando-se o JAMS Builder ou JAMS Launcher [Arquivo -> Carregar modelo].

Workspace directory: A localização do diretório de área de trabalho na máquina local.

Time Interval: O intervalo de tempo em que o modelo é executado.

Parameter file: A localização da HRU e arquivo de parâmetro de trecho fluvial.

Efficiency: Pode-se aqui fornecer um período de tempo diferente para a estimativa de eficiência.

Nota Importante: O modelo J2000 foi aplicado com sucesso na bacia do rio Kosi Dudh, como parte de uma pesquisa de doutorado. O modelo hidrológico J2000 calibrado e validado foi ainda usado para avaliar o impacto da mudança no uso da terra sobre o regime hidrológico. Dois cenários de mudança do uso da terra hipotéticos foram implementados e as informações do uso da terra do arquivo de parâmetro da HRU foram alteradas de acordo com a quantidade do impacto dessa mudança em diferentes processos hidrológicos. Além disso, o impacto da mudança climática sobre o regime hidrológico também foi analisado usando os dados climáticos regionais do modelo na bacia hidrográfica Kosi Dudh. A descrição destas análises são fornecidas em Nepal (2012).

Fórum

É provável que durante a utilização do modelo e tutorial, possa haver problemas e mensagens de erro. Em tais casos, sugere-se entrar em contato com o fórum de discussão ILMS, onde os usuários do modelo e comunidades de desenvolvedores podem ser contactados.

Sistema Integrado de Gestão de Terras (ILMS) Fórum de discussão:

http://ilms.uni-jena.de/ilms/board/index.php?sid=f545b1932d03a6781393eea0fed040e3

O fórum de discussão ILMS foi formulado para a discussão dos vários componentes do software ILMS. Para o modelo ILMS, o seguinte fórum foi criado:

http://ilms.uni-jena.de/ilms/board/viewforum.php?f=6

É necessário se registrar no fórum ILMS Registrar aqui para deixar uma nova pergunta e também para acompanhar as postagens.

=Bibliografia e Leitura para aprofundamento=