Skip to contents

Trophic State Index (Lamparelli)

Source: R/iet_lamparelli.R
iet_lamparelli.Rd

Computa o Indice do Estado Trofico de Lamparelli (2004), como adotado pela CETESB, a partir de fosforo total e clorofila-a. O indice agregado e a media de IET(PT) e IET(CL) disponiveis na linha.

Pode receber um data.frame como primeiro argumento (ver Detalhes).

Usage

iet_lamparelli(
  tp = NULL,
  chla = NULL,
  sd = NULL,
  ambiente = c("rio", "reservatorio"),
  .keep_ids = FALSE,
  add_status = TRUE,
  locale = c("pt", "en"),
  ...
)

Arguments

tp

Fosforo total (mg/L) ou um data.frame contendo colunas tp (ug/L) ou p_total (mg/L), chla ou clorofila (ug/L), e sd ou secchi (m). Se for data.frame, chla e sd devem ser NULL.

chla

Clorofila-a (ug/L).

sd

Profundidade do disco de Secchi (m).

ambiente

Tipo de ambiente: "rio" ou "reservatorio".

.keep_ids

Logico; quando data.frame, vincula colunas de ID (rio, ponto, data, lat, lon). Padrao FALSE.

add_status

Logico; se TRUE (padrao), adiciona a coluna TSI_status com a classificacao qualitativa (Lamparelli).

locale

Idioma de TSI_status: "pt" (padrao) ou "en".

...

Reservado para uso futuro (ignorado).

Value

Um data.frame com colunas (quando aplicavel):

  • IET_TP — componente de fosforo total (0-100).

  • IET_Chla — componente de clorofila-a (0-100).

  • IET_Secchi — componente auxiliar de Secchi (0-100), fora do agregado.

  • IET_Lamp — indice Lamparelli agregado: media de IET_TP e IET_Chla (0-100).

  • TSI_status — classe qualitativa (quando add_status=TRUE).

  • ambiente — tipo de ambiente informado.

Details

Equacoes (CETESB 2023, Apendice D, secao 4, p. 10; PT e CL em ug/L, ln = logaritmo natural):

  • Rios: IET(CL) = 10 (6 - ((-0.7 - 0.6 ln CL) / ln 2)) - 20; IET(PT) = 10 (6 - ((0.42 - 0.36 ln PT) / ln 2)) - 20.

  • Reservatorios: IET(CL) = 10 (6 - ((0.92 - 0.34 ln CL) / ln 2)); IET(PT) = 10 (6 - ((1.77 - 0.42 ln PT) / ln 2)).

Conferidas contra os cortes das Tabelas 5 e 6 da mesma fonte. A CETESB nao usa a transparencia no IET; quando sd e informado, IET_Secchi e devolvido como componente auxiliar (definicao de Carlson, 10 (6 - ln SD / ln 2)) e nao entra em IET_Lamp. Ate a v0.10.0 a funcao usava 10 + 10 log10(x) para PT e CL, igual para rio e reservatorio, o que classificava praticamente qualquer amostra como ultraoligotrofica.

Entradas com virgula decimal (ex.: "3,2") ou desigualdades (ex.: "<0,1") sao convertidas com seguranca por helpers internos. Se houver apenas p_total (mg/L), e convertida para tp (ug/L) via tp = p_total * 1000.

Os componentes e o indice agregado sao limitados ao intervalo [0, 100] para consistencia com as figuras e tabelas do pacote/artigo.

References

Lamparelli, M. C. (2004). Grau de trofia em corpos d'agua do Estado de Sao Paulo: avaliacao dos metodos de monitoramento. Tese (Doutorado), Instituto de Biociencias, Universidade de Sao Paulo.

CETESB (2023). Apendice D: Indices de Qualidade das Aguas e Criterios de Avaliacao da Qualidade dos Sedimentos. In: Qualidade das Aguas Interiores no Estado de Sao Paulo. Secao 4, pp. 9-11.

See also

iet_carlson(), iqa(), conama_check()