Metodo para el calculo de Datos de Actividad (DA)

La gerencia de Monitoreo Reporte y Verificación (MRV) de la Comisión Nacional Forestal empleó la metodología propuesta por el Panel Intergubernamental de Cambio Climatico (IPCC) para generar sus inventarios de gases de efecto invernadero. Se pone en practica el volumen 4 de las directices del IPCC. Como formula general:

\[ Emisiones=DA\times FE \]

Para comenzar se generarón Datos de Actividad, para ello se estimó una serie temporal de transiciones de uso de suelo del año 2001 al año 2022 en el estado de Jalisco, México. Se analizó el insumo principal, la malla nacional densificada del estado de Jalisco, un data frame o tabla llamada bdCambios con una serie temporal de cambios en el uso de suelo del año 2000 al 2022.

Así mismo se debe leer otro imsumo importante, una tabla con las 36 posibles transiciones de uso de suelo llamado transición. La funcionalidad de esta tabla es contabilizar los valores únicos para cada transición y eventualmente generar los conteos de cambios en el uso de suelo en el estado de Jalisco.

data.frame(bdCambios)

data.frame(transicion)

Una vez que se analizó la malla densificada, en este caso del estado de Jalisco, se extrajerón los valores unicos de la columna ClaseMalla del data frame bdCambios y se almacenarón en el vector idClaseMalla. Esto fue útil para identificar los diferentes tipos de clase malla presentes en los datos. Se imprimierón los valores en la consola, para verificar los datos antes de continuar.

idClaseMalla <- unique(bdCambios$ClaseMalla)
print(idClaseMalla)

## [1] "Nacional&Densificada"

Se extrajerón los valores únicos de la columna Condición del mismo data frame, bdCambios, se convirtierón en caracteres (si es que no lo estan) y se almacenarón en el vector condiciones. De igual manera se imprimierón los valores unicos de este vector.

condiciones <- c(unique(as.character(bdCambios$Condicion)))
print(condiciones)

##  [1] "Permanencia TF"                                                                                                  
##  [2] "Permanencia P"                                                                                                   
##  [3] "Permanencia TA"                                                                                                  
##  [4] "Deforestacion"                                                                                                   
##  [5] "Recuperacion Pradera"                                                                                            
##  [6] "Perdida de Pradera"                                                                                              
##  [7] "Permanencia AH"                                                                                                  
##  [8] "Otros Cambios"                                                                                                   
##  [9] "Permanencia H"                                                                                                   
## [10] "Recuperacion Forestal"                                                                                           
## [11] "Permanencia P (no alcanzo umbral Bosque)"                                                                        
## [12] "Permanencia OT"                                                                                                  
## [13] "Permanencia (DFA)"                                                                                               
## [14] "Permanencia P (Recuperacion Forestal: es posible que alcance los umbrales de la definicion Nacional de Bosques)" 
## [15] "Recuperacion de Pradera"                                                                                         
## [16] "Recuperacion Forestal: alcanzo los umbrales de la definicion Nacional de Bosques"                                
## [17] "Permanencia OT (Recuperacion Forestal: es posible que alcance los umbrales de la definicion Nacional de Bosques)"

Misma dinamica, pero ahora se extrajerón los valores unicos de la columna Equidist_Malla_Nacional del mismo data frame, bdCambios y se almacenarón en el vector malla Se imprimierón los valores de este vector en la consola.

malla <- c(unique(as.character(bdCambios$Equidist_Malla_Nacional)))#Selecciona las mallas de superficies
print(malla)

## [1] "5 x 5"     "2.5 x 2.5" "10 x 10"

Se extrajerón los valores unicos de la columna DESECON del mismo data frame bdCambios; se convirtierón en caracteres (si no lo estan) y fuerón almacenados en el vector ecorregiones. Se debierón imprmir los valores unicos de este vector en la consola.

ecorregiones <- c(unique(as.character(bdCambios$DESECON1)))#selecciona las ecorregiones
print(ecorregiones)

## [1] "Selvas Calido-Secas"                 "Sierras Templadas"                  
## [3] "Elevaciones Semiaridas Meridionales" "Selvas Calido-Humedas"

Se procedió a generar los conteos, primero fue necesarió crear la dirección de trabajo donde serán almacenados los conteos de transición de uso de suelo.

direccion <- "C:/Users/ABAX/Desktop/CONAFOR/Tasa_Jliasco_2024/Conteos_Jalisco_2000-2022"

Se ejecutó la función for para crear un bucle que haga el recuento de las transiciones del uso de suelo.

for(ClaseMalla in idClaseMalla ){#ejecuta los filtros y guarda el resultado e
  setwd(direccion)
  filtroClaseMalla <- bdCambios[bdCambios[,"ClaseMalla"]==ClaseMalla,]
  ClaseMalla <- unique(as.character(filtroClaseMalla$ClaseMalla))
  
  dir.create(ClaseMalla)
  setwd(ClaseMalla)
  
  ecorregiones <- c(unique(as.character(filtroClaseMalla$DESECON1)))#Cambiar aqui
  for(ecorregion in ecorregiones){
    filtroEcorregion <- filtroClaseMalla[filtroClaseMalla[,"DESECON1"]==ecorregion,]#cambiar aqui
    dir.create(ecorregion)
    setwd(ecorregion)
    malla <- c(unique(as.character(filtroEcorregion$Equidist_Malla_Nacional)))#cambiar aqui
    for(area in malla){
      filtroMalla= filtroEcorregion[filtroEcorregion[,"Equidist_Malla_Nacional"]==area,]#cambiar aqui
      tabla <- matrix(nrow=0,ncol = 26)
      for(j in 1:nrow(transicion)){
        numeroColumna = 5#cambiar si cambia el numero de columnas
        columna1 <-as.character(transicion[j,1])
        columna2 <- as.character(transicion[j,2])
        concat <- paste(columna1,"-",columna2)
        datosFiltrados <- c(ClaseMalla,ecorregion,area,concat)
        for(i in 1:26){
          numeroColumna = numeroColumna + 1
          numeroFilas <- nrow((filtroMalla[filtroMalla[,numeroColumna] == columna1 & filtroMalla[,numeroColumna +1] == columna2,]))
          datosFiltrados <- c(datosFiltrados, numeroFilas)
        }
        tabla<-rbind(tabla,as.character(datosFiltrados))
        datosFiltrados <- c()
      }
      dfTabla <-as.data.frame(tabla)
      names(dfTabla)<-c('ClaseMalla','Ecor','malla','tansicion','00-01','01-02','02-03','03-04',
                        '04-05','05-06','06-07','07-08','08-09','09-10','10-11','11-12',
                        '12-13','13-14','14-15','15-16','16-17','17-18','18-19','19-20','20-21','21-22')
      print(dfTabla)
      nameCsv = paste("Equidistan",area,".csv")
      write.csv(dfTabla, file=nameCsv)
    }
    direccion2<-paste(direccion,"/",ClaseMalla,sep = "")
    #print(direccion2)
    setwd(direccion2)
  }
}

Conteo de muestras

Ahora R almacenó los conteos de muestras en la carpeta que fue establecida anteriormente. Existe un conteo especifico para cada ecorregión y sus respectivas equidistancias. El estado de Jalisco tiene 4 ecorregiones nivel 1:

Elevaciones Semiardias Meridionales

2.5 x 2.5
5 x 5
10 x 10

Selvas Calido Secas

2.5 x 2.5
5 x 5
10 x 10

Selvas Calido Humedas

2.5 x 2.5
5 x 5
10 x 10

Sierras Templadas

2.5 x 2.5
5 x 5
10 x 10

Lo siguiente fue unir en una sola matriz, en otra hoja de cálculo todas las otras hojas de cálculo generadas por el Script de R. Si son 4 ecorregiones con 3 equidistancias tendremos 12 hojas de cálculo que deben ir integradas en una única hoja de cálculo. La cual se llamó (Conteos_Jalisco_2000_2022). Lo que se hizó fue copiar los datos de cada hoja de cálculo en orden y pegarlo en la hoja indicada, hasta tener todos los conteos por orden en una sola tabla.

data.frame(conteos_Jalisco_2000_2022)

Obtenidos los conteos de transiciones, se contarón las proyecciones de las mismas, para determinar eventualmente mediante las muestras:

Matriz de proporción
Matriz de Áreas estimadas
Matriz de desviación estandar de áreas
Matriz de incertidumbres

La estimación de las áreas deforestadas se realizó siguiendo el enfoque de superficies por medio de proporciones indicado en el capitulo 3 del Volumen 4 de las directrices del IPCC (2006) para lo cual, primero se generó un matriz en la que se identifica la evolución de las clases de cada una de las parcelas (bdCambios) a lo largo de los 21 periodos anuales de anslisis, dede el año 2000 hasta el año 2022.

Identificación de parcelas y años con cambios de clase

La base de datos con los registros de las clses INEGI y fechas de cmabios de clases de las parcelas (bdCambios) incluyerón varias vaiables. Sin embargo, para la evaluación de las tasas de deforestación solo se requirieron las variables sobre la calse INEGI de la parcela al año base 2001 al año 2022 y las correspondientes a aquellos casos en lo que se presentarón multiples cambios, ademas de las fechas de cambio (CONAFOR, 2020)

Las clases reportadas en el año 2022, 2000 y los diferentes años de cambio se etiquetarón a nivel de las clases de INEGI y fueron traducidas a las clases IPCC. una vez obtenidas las clases IPCC del año 2000, 2022 y de los diferentes eventos de cambio, fue posible identificar las parcelas que permanerierón con el mismo tipo de clase, las que presentarón cambio, así como su respectiva dirección de cambio. Se consideró como deforestación las siguientes transiciónes:

Tierra Forestal a Pradera (FL-GL)
Tierra Forestal a Agricultura (FL-CL)
Tierra Forestal a Asentamiento Humano (FL-SL)
Tierra Forestal a Cuerpo de Agua (FL-WL)
Tierra Forestal a Otras Tierras (FL-OL)

Con base en los criterios descritos se desarrolló una matriz anidada con otras matrices en la que para cada parcela se registró la evolución de su clase IPCC de manera anual desde el año 2001 hasta el año 2022. De manera que, si la parcela registró un cambio de clase en un año en particular, entonces se registró un evento de cambio.

Matriz de conteo de muestras

La primer matriz, que contabilizó el numero de parcelas totales, es el vaciado del data frame: (Conteos_Jalisco_2000_2022). Por lo que indicó las parcelas y su historial de cambios o permanencias en la serie temporal. Es indispensable contabilizar al final de los vaaciados que el numero total de muestras coincida con las muestras del data frame: bdCambios para cercionarse de que se estan llevando los conteos correctamente y evitar errores más adelante.

Una vez vaciados los conteos, las matrices de excell ya tienen las formulas predeterminadas para realizar las operaciones correspondientes para la estimación de los siguientes parametros. Sin embargo, se desglozó cada matriz para un óptimo entendimiento. Como se comentaba, la primera es la matriz: Conteo de muestras:

Para este caso de estudio se utilizó la ecorregión elevaciones semiaridas meridionales y la intensidad de muestro 2.5 x 2.5. Tengase en cuenta que para generar las tasas de deforestación que evebtualmente determinan las emisiones de co2, es forzoso determinar estos pasos para cada ecorregión y las tres intensidades de muestreo de la entidad federativa que se busca estimar sus emisiones.

Matriz conteo de muestras
Eleveaciones Semiaridas-Meridionales
Intensidad de muestreo: 2.5 x 2.5
Superficie: 169,495 hectareas
total de parcelas muestreadas (n): 263
Trancicion_IPCC	00-01_(n1,i)	01-02_(n2,i)	02-03_(n3,i)	03-04_(n4,i)	04-05_(n5,1)	05-06_(n6,i)	06-07_(n7,i)	07-08_(n8,i)	08-09_(n9,i)	09-10_(n10,i)	10-11_(n11,i)	11-12_(n12,i)	12-13_(n13,i)	13-14_(n14,i)	14-15_(n15,i)	15-16_(n16,i)	16-17_(n17,i)	17-18_(n18,i)	18-19_(n19,i)	19-20_(n20,i)	20-21_(n21,i)	21-22_(n22,i)
FL - FL	168	168	168	167	167	167	167	167	167	167	167	167	167	167	166	166	166	166	166	166	166	166
FL - GL	0	0	0	0	0	0	0	0	0	0	0	0	0	0	1	0	0	0	0	0	0	0
FL - CL	0	0	0	1	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
FL - WL	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
FL - SL	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
FL - OL	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
GL - FL	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
GL - GL	80	80	80	80	80	80	80	80	80	79	79	79	79	79	78	79	78	78	78	79	79	79
GL - CL	0	0	0	0	0	0	0	0	0	1	0	0	0	0	1	0	1	0	0	0	0	0
GL - WL	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
GL - SL	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
GL - OL	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
CL - FL	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
CL - GL	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	1	0	0	0
CL - CL	11	11	11	11	12	12	12	12	12	12	13	13	13	13	13	14	14	15	14	14	14	14
CL - WL	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
CL - SL	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
CL - OL	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
WL - FL	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
WL - GL	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
WL - CL	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
WL - WL	2	2	2	2	2	2	2	2	2	2	2	2	2	2	2	2	2	2	2	2	2	2
WL - SL	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
WL - OL	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
SL - FL	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
SL - GL	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
SL - CL	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
SL - WL	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
SL - SL	2	2	2	2	2	2	2	2	2	2	2	2	2	2	2	2	2	2	2	2	2	2
SL - OL	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
OL - FL	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
OL - GL	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
OL - CL	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
OL - WL	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
OL - SL	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
OL - OL	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0

Matriz de proporciones

Para estimar el área de cambio fue necesario obtener la proporción \((p_i = \dfrac{ni}{n})\) de las parcelas de cambio (\(n_i\)) con respecto al numero de parcelas muestreadas \((n)\) en un estrato especifico, (en nuestro caso, Elevaciones Semiaridas-Meridionales, equidistancia o intensidad de muestreo 2.5 x 2.5). Por ejemplo, para el el caso FL-FL 00-01, n1: \[ (p = \dfrac{168}{263}) \quad \therefore \quad p = 0.6388 \] Por lo tanto, como se ha ido mencionando, para cada estrato de “Ecorregión-Equidistancia” se generó una matriz de proporciones obtendida de la razón de las parcelas de cambio respecto al total de parcelas muestreadas. Posteriormente, esa matriz de proporciones fue ponderada por el área total de cada “Ecorregión-Equidistancia” y con ello se obtuvó el área total de cambio para cada uno de los periodos anuales y transiciones del IPCC.

Matriz conteo de Proporciones
Eleveaciones Semiaridas-Meridionales
Intensidad de muestreo: 2.5 x 2.5
Superficie: 169,495 hectareas
total de parcelas muestreadas (n): 263
Trancicion_IPCC	00-01_(n1,i)	01-02_(n2,i)	02-03_(n3,i)	03-04_(n4,i)	04-05_(n5,1)	05-06_(n6,i)	06-07_(n7,i)	07-08_(n8,i)	08-09_(n9,i)	09-10_(n10,i)	10-11_(n11,i)	11-12_(n12,i)	12-13_(n13,i)	13-14_(n14,i)	14-15_(n15,i)	15-16_(n16,i)	16-17_(n17,i)	17-18_(n18,i)	18-19_(n19,i)	19-20_(n20,i)	20-21_(n21,i)	21-22_(n22,i)
FL - FL	0.6388	0.6388	0.6388	0.6350	0.6350	0.6350	0.6350	0.6350	0.6350	0.6350	0.6350	0.6350	0.6350	0.6350	0.6312	0.6312	0.6312	0.6312	0.6312	0.6312	0.6312	0.6312
FL - GL	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0038	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000
FL - CL	0.0000	0.0000	0.0000	0.0038	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000
FL - WL	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000
FL - SL	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000
FL - OL	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000
GL - FL	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000
GL - GL	0.3042	0.3042	0.3042	0.3042	0.3042	0.3042	0.3042	0.3042	0.3042	0.3004	0.3004	0.3004	0.3004	0.3004	0.2966	0.3004	0.2966	0.2966	0.2966	0.3004	0.3004	0.3004
GL - CL	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0038	0.0000	0.0000	0.0000	0.0000	0.0038	0.0000	0.0038	0.0000	0.0000	0.0000	0.0000	0.0000
GL - WL	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000
GL - SL	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000
GL - OL	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000
CL - FL	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000
CL - GL	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0038	0.0000	0.0000	0.0000
CL - CL	0.0418	0.0418	0.0418	0.0418	0.0456	0.0456	0.0456	0.0456	0.0456	0.0456	0.0494	0.0494	0.0494	0.0494	0.0494	0.0532	0.0532	0.0570	0.0532	0.0532	0.0532	0.0532
CL - WL	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000
CL - SL	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000
CL - OL	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000
WL - FL	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000
WL - GL	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000
WL - CL	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000
WL - WL	0.0076	0.0076	0.0076	0.0076	0.0076	0.0076	0.0076	0.0076	0.0076	0.0076	0.0076	0.0076	0.0076	0.0076	0.0076	0.0076	0.0076	0.0076	0.0076	0.0076	0.0076	0.0076
WL - SL	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000
WL - OL	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000
SL - FL	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000
SL - GL	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000
SL - CL	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000
SL - WL	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000
SL - SL	0.0076	0.0076	0.0076	0.0076	0.0076	0.0076	0.0076	0.0076	0.0076	0.0076	0.0076	0.0076	0.0076	0.0076	0.0076	0.0076	0.0076	0.0076	0.0076	0.0076	0.0076	0.0076
SL - OL	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000
OL - FL	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000
OL - GL	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000
OL - CL	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000
OL - WL	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000
OL - SL	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000
OL - OL	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000	0.0000

Matriz de áreas estimadas

Para estimar el área de cambio fue necesarió obtener la proporción de las parcelas de cambio con respecto al numero de las parcelas muestreadas en un estrato especifico, ponderado por el área total del estrato especifico. El área deforestada \((AD_i)\) de un estrato específico se obtuvó de acuerdo con la siguiente formula \[ \widehat{AD}_i = \dfrac{ni}{n} \times A_T \] Donde:

\(\widehat{AD}_i\) representa el área todoal estimada de deforestación
\(ni\) es el i-ésimo valor de la proporción
\(AT\) es el área total del segmento evaluado, su equidistancia y ecorregión especifica.

Para calular la matriz de áreas estimadas se debe añadir la superficie en hectáreas total de cada ecorregión según su equidistancia. Tomamos la tabla SuperficiesEcorregionEstado y se debe utilizar el área totoal de la ecorregión y equidistancia corresponiente para poder realizar el caluclo de proporción del área.

data.frame(Superficie)

Matriz areas estimadas (ha)
Eleveaciones Semiaridas-Meridionales
Intensidad de muestreo: 2.5 x 2.5
Superficie: 169,495 hectareas
total de parcelas muestreadas (n): 263
Trancicion_IPCC	00-01_(n1,i)	01-02_(n2,i)	02-03_(n3,i)	03-04_(n4,i)	04-05_(n5,1)	05-06_(n6,i)	06-07_(n7,i)	07-08_(n8,i)	08-09_(n9,i)	09-10_(n10,i)	10-11_(n11,i)	11-12_(n12,i)	12-13_(n13,i)	13-14_(n14,i)	14-15_(n15,i)	15-16_(n16,i)	16-17_(n17,i)	17-18_(n18,i)	18-19_(n19,i)	19-20_(n20,i)	20-21_(n21,i)	21-22_(n22,i)
FL - FL	108270.5	108270.5	108270.5	107626.0	107626.0	107626.0	107626.0	107626.0	107626.0	107626.0	107626.0	107626.0	107626.0	107626.0	106981.5	106981.5	106981.5	106981.5	106981.5	106981.5	106981.5	106981.5
FL - GL	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	644.5	0.0	0.0	0.0	0.0	0.0	0.0	0.0
FL - CL	0.0	0.0	0.0	644.5	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0
FL - WL	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0
FL - SL	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0
FL - OL	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0
GL - FL	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0
GL - GL	51557.4	51557.4	51557.4	51557.4	51557.4	51557.4	51557.4	51557.4	51557.4	50912.9	50912.9	50912.9	50912.9	50912.9	50268.4	50912.9	50268.4	50268.4	50268.4	50912.9	50912.9	50912.9
GL - CL	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	644.5	0.0	0.0	0.0	0.0	644.5	0.0	644.5	0.0	0.0	0.0	0.0	0.0
GL - WL	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0
GL - SL	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0
GL - OL	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0
CL - FL	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0
CL - GL	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	644.5	0.0	0.0	0.0
CL - CL	7089.1	7089.1	7089.1	7089.1	7733.6	7733.6	7733.6	7733.6	7733.6	7733.6	8378.1	8378.1	8378.1	8378.1	8378.1	9022.5	9022.5	9667.0	9022.5	9022.5	9022.5	9022.5
CL - WL	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0
CL - SL	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0
CL - OL	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0
WL - FL	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0
WL - GL	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0
WL - CL	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0
WL - WL	1288.9	1288.9	1288.9	1288.9	1288.9	1288.9	1288.9	1288.9	1288.9	1288.9	1288.9	1288.9	1288.9	1288.9	1288.9	1288.9	1288.9	1288.9	1288.9	1288.9	1288.9	1288.9
WL - SL	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0
WL - OL	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0
SL - FL	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0
SL - GL	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0
SL - CL	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0
SL - WL	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0
SL - SL	1288.9	1288.9	1288.9	1288.9	1288.9	1288.9	1288.9	1288.9	1288.9	1288.9	1288.9	1288.9	1288.9	1288.9	1288.9	1288.9	1288.9	1288.9	1288.9	1288.9	1288.9	1288.9
SL - OL	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0
OL - FL	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0
OL - GL	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0
OL - CL	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0
OL - WL	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0
OL - SL	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0
OL - OL	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0	0.0

Matriz de desviación estandar de áreas

Varianza: \[ AT \times \dfrac{p1(1-p1)}{n-1} \] Desviación estandar: \[ AT \times \sqrt{\dfrac{p1(1-p1)}{n-1}} \quad \therefore \quad Area total \times \sqrt{\dfrac{valor proporcion(1-valor proporción)}{tamaño de la muestra-1}} \]

AT: Representa el área total de cada equidistancia
P1: Es la proporción muestral, la fracción de la muestra que posee cierta caracteristica o satisface un criterio (cambio de suelo o permanencia)
N: Es el tamaño de la muestra , en este caso para su equidistancia y ecorregión especifica
La expresión \((1-p1)\) representa la proporción de la muestra que no posee la caracteristica en cuestión
El termino \(p1 \times (1-p1)\) calcula la varianza de ña proporción muestral, es decir, mide cuanto varia la proporción muestral de la verdadera proporción de la población.
La raíz cuadrada de esta varianza denota la desviación estandar de la proporción muestral Esta matriz es importante ya que será un insumo clave para el calculo de incertidumbres en su primera etapa. Ya que las incertidumbres son estimadas en dos etapas.

Matriz desviacion estandar
Eleveaciones Semiaridas-Meridionales
Intensidad de muestreo: 2.5 x 2.5
Superficie: 169,495 hectareas
total de parcelas muestreadas (n): 263
Trancicion_IPCC	00-01_(n1,i)	01-02_(n2,i)	02-03_(n3,i)	03-04_(n4,i)	04-05_(n5,1)	05-06_(n6,i)	06-07_(n7,i)	07-08_(n8,i)	08-09_(n9,i)	09-10_(n10,i)	10-11_(n11,i)	11-12_(n12,i)	12-13_(n13,i)	13-14_(n14,i)	14-15_(n15,i)	15-16_(n16,i)	16-17_(n17,i)	17-18_(n18,i)	18-19_(n19,i)	19-20_(n20,i)	20-21_(n21,i)	21-22_(n22,i)
FL - FL	5030	5030	5030	5041	5041	5041	5041	5041	5041	5041	5041	5041	5041	5041	5052	5052	5052	5052	5052	5052	5052	5052
FL - GL	0	0	0	0	0	0	0	0	0	0	0	0	0	0	644	0	0	0	0	0	0	0
FL - CL	0	0	0	644	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
FL - WL	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
FL - SL	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
FL - OL	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
GL - FL	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
GL - GL	4817	4817	4817	4817	4817	4817	4817	4817	4817	4800	4800	4800	4800	4800	4783	4800	4783	4783	4783	4800	4800	4800
GL - CL	0	0	0	0	0	0	0	0	0	644	0	0	0	0	644	0	644	0	0	0	0	0
GL - WL	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
GL - SL	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
GL - OL	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
CL - FL	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
CL - GL	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	644	0	0	0
CL - CL	2096	2096	2096	2096	2185	2185	2185	2185	2185	2185	2270	2270	2270	2270	2270	2351	2351	2428	2351	2351	2351	2351
CL - WL	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
CL - SL	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
CL - OL	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
WL - FL	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
WL - GL	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
WL - CL	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
WL - WL	910	910	910	910	910	910	910	910	910	910	910	910	910	910	910	910	910	910	910	910	910	910
WL - SL	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
WL - OL	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
SL - FL	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
SL - GL	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
SL - CL	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
SL - WL	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
SL - SL	910	910	910	910	910	910	910	910	910	910	910	910	910	910	910	910	910	910	910	910	910	910
SL - OL	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
OL - FL	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
OL - GL	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
OL - CL	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
OL - WL	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
OL - SL	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
OL - OL	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0

Matriz de primera etapa de propagación de incertidumbres

El capitulo 3 de incertidumbres del volumen 1 de las directrices del IPCC (2006) indica que es una buena practica cuantificar las fuentes de incertidumbres. Para este caso las incertidumbres fuerón estimadas en dos etapas. La primera consistió en la estimación de las incertidumbres a nivel de las áreas de cambio a nivel de “Ecorregión-Equidistancia” para cada transición del IPCC. Como se muestra en la siguiente matriz, sin embargo, la ejemplificaición se reduce a una sola ecorregión que es Elevaciones Semiaridas-Meridionales. La sgunda etapa consistió en la propagación de las incertidumbres a travéz del metodo 1 de las directrices del IPCC.

Las incertidumbres de las áreas de cambio (para cada Ecorregión Qquidistancia-Transición IPCC) obtenidas mediante método de proporciones se estimarón de acuerdo con lo siguiente:

\[ U(AD_i) \dfrac{z_{\alpha / 2} * s(AD_i)}{AD_i} \times 100 \] Donde \(z_{\alpha / 2}\) es el percentil 95% del modelo de distribución empírico que ajusta los datos, y

\[ s(AD_i) = A_T * \sqrt{\dfrac{pi(1-pi)}{n-1}} \] \(S(A_i)\) es la desviación estandar del área deforestada estimada en la subsección anterior con los insumos ya definidos. De esta manera se estimarón las incertidumbres en la primera etapa. La segunda etapa se estimará en la sección de conteos anuales, ya que el área total deforestada es la suma de las áreas de cambio por Ecorregión-Equidistancia para las transiciones:

FL-GL
FL-CL
FL-WL
FL-SL
FL-OL

La segunda etapa de propagación de incertidumbres fue porpagada a través del metodo 1 del IPCC por adición:

\[ U(AD_T) = \sqrt{\frac{(U(\widehat{AD_1}) \times \widehat{AD_1})^2 + (U(\widehat{AD_2}) \times \widehat{AD_2})^2 + \dots + (U(\widehat{AD_n}) \times \widehat{AD_n})^2}{\mid \widehat{AD_1}) + \widehat{AD_2} +...+\widehat{AD}_n \mid}} \]

Matriz de incertidumbres
Eleveaciones Semiaridas-Meridionales
Intensidad de muestreo: 2.5 x 2.5
Superficie: 169,495 hectareas
total de parcelas muestreadas (n): 263
Trancicion_IPCC	00-01_(n1,i)	01-02_(n2,i)	02-03_(n3,i)	03-04_(n4,i)	04-05_(n5,1)	05-06_(n6,i)	06-07_(n7,i)	07-08_(n8,i)	08-09_(n9,i)	09-10_(n10,i)	10-11_(n11,i)	11-12_(n12,i)	12-13_(n13,i)	13-14_(n14,i)	14-15_(n15,i)	15-16_(n16,i)	16-17_(n17,i)	17-18_(n18,i)	18-19_(n19,i)	19-20_(n20,i)	20-21_(n21,i)	21-22_(n22,i)
FL - FL	9.11	9.11	9.11	9.18	9.18	9.18	9.18	9.18	9.18	9.18	9.18	9.18	9.18	9.18	9.26	9.26	9.26	9.26	9.26	9.26	9.26	9.26
FL - GL	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	196.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00
FL - CL	0.00	0.00	0.00	196.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00
FL - WL	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00
FL - SL	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00
FL - OL	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00
GL - FL	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00
GL - GL	18.31	18.31	18.31	18.31	18.31	18.31	18.31	18.31	18.31	18.48	18.48	18.48	18.48	18.48	18.65	18.48	18.65	18.65	18.65	18.48	18.48	18.48
GL - CL	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	196.00	0.00	0.00	0.00	0.00	196.00	0.00	196.00	0.00	0.00	0.00	0.00	0.00
GL - WL	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00
GL - SL	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00
GL - OL	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00
CL - FL	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00
CL - GL	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	196.00	0.00	0.00	0.00
CL - CL	57.96	57.96	57.96	57.96	55.38	55.38	55.38	55.38	55.38	55.38	53.10	53.10	53.10	53.10	53.10	51.07	51.07	49.24	51.07	51.07	51.07	51.07
CL - WL	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00
CL - SL	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00
CL - OL	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00
WL - FL	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00
WL - GL	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00
WL - CL	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00
WL - WL	138.33	138.33	138.33	138.33	138.33	138.33	138.33	138.33	138.33	138.33	138.33	138.33	138.33	138.33	138.33	138.33	138.33	138.33	138.33	138.33	138.33	138.33
WL - SL	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00
WL - OL	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00
SL - FL	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00
SL - GL	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00
SL - CL	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00
SL - WL	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00
SL - SL	138.33	138.33	138.33	138.33	138.33	138.33	138.33	138.33	138.33	138.33	138.33	138.33	138.33	138.33	138.33	138.33	138.33	138.33	138.33	138.33	138.33	138.33
SL - OL	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00
OL - FL	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00
OL - GL	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00
OL - CL	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00
OL - WL	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00
OL - SL	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00
OL - OL	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00

Matriz de incertidumbre ajustada

\[ (A \times U)^2 \] \[area total \times incertidumbres^2\] Multiplicar el área total de cada ecorregión o transición por la incertidumbre al cuadrado tiene un fundamento en el caluclo de errores en variables dependientes del tamaño o área. La razón de usar el cuadrado de la incertidumbre es que el error total o incertidumbre total se porpaga en función de la variabilidad de los datos en cada muestra. Al elevar al cuadrado, las diferencias mas significativas (mayor incertidumbre) pesan más, refñejando los errores de manera proporcional a su magnitud. Adicioalmente, la matriz de incertidumbre, que es una matriz auxiliar, ayuda a evitar la subestimación o sobreestimación de la incertidumbre total.

Matriz de incertidumbre ajustada
Eleveaciones Semiaridas-Meridionales
Intensidad de muestreo: 2.5 x 2.5
Superficie: 169,495 hectareas
total de parcelas muestreadas (n): 263
Trancicion_IPCC	00-01_(n1,i)	01-02_(n2,i)	02-03_(n3,i)	03-04_(n4,i)	04-05_(n5,1)	05-06_(n6,i)	06-07_(n7,i)	07-08_(n8,i)	08-09_(n9,i)	09-10_(n10,i)	10-11_(n11,i)	11-12_(n12,i)	12-13_(n13,i)	13-14_(n14,i)	14-15_(n15,i)	15-16_(n16,i)	16-17_(n17,i)	17-18_(n18,i)	18-19_(n19,i)	19-20_(n20,i)	20-21_(n21,i)	21-22_(n22,i)
FL - FL	971953000000	971953000000	971953000000	976338000000	976338000000	976338000000	976338000000	976338000000	976338000000	976338000000	976338000000	976338000000	976338000000	976338000000	980600000000	980600000000	980600000000	980600000000	980600000000	980600000000	980600000000	980600000000
FL - GL	0	0	0	0	0	0	0	0	0	0	0	0	0	0	15955615375	0	0	0	0	0	0	0
FL - CL	0	0	0	15955615375	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
FL - WL	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
FL - SL	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
FL - OL	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
GL - FL	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
GL - GL	891566000000	891566000000	891566000000	891566000000	891566000000	891566000000	891566000000	891566000000	891566000000	885232000000	885232000000	885232000000	885232000000	885232000000	878777000000	885232000000	878777000000	878777000000	878777000000	885232000000	885232000000	885232000000
GL - CL	0	0	0	0	0	0	0	0	0	15955615375	0	0	0	0	15955615375	0	15955615375	0	0	0	0	0
GL - WL	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
GL - SL	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
GL - OL	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
CL - FL	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
CL - GL	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	15955615375	0	0	0
CL - CL	168813000000	168813000000	168813000000	168813000000	183429000000	183429000000	183429000000	183429000000	183429000000	183429000000	197923000000	197923000000	197923000000	197923000000	197923000000	212295000000	212295000000	226545000000	212295000000	212295000000	212295000000	212295000000
CL - WL	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
CL - SL	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
CL - OL	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
WL - FL	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
WL - GL	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
WL - CL	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
WL - WL	31789432160	31789432160	31789432160	31789432160	31789432160	31789432160	31789432160	31789432160	31789432160	31789432160	31789432160	31789432160	31789432160	31789432160	31789432160	31789432160	31789432160	31789432160	31789432160	31789432160	31789432160	31789432160
WL - SL	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
WL - OL	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
SL - FL	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
SL - GL	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
SL - CL	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
SL - WL	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
SL - SL	31789432160	31789432160	31789432160	31789432160	31789432160	31789432160	31789432160	31789432160	31789432160	31789432160	31789432160	31789432160	31789432160	31789432160	31789432160	31789432160	31789432160	31789432160	31789432160	31789432160	31789432160	31789432160
SL - OL	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
OL - FL	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
OL - GL	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
OL - CL	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
OL - WL	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
OL - SL	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
OL - OL	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0

Esimación de áreas deforestadas anuales por transición y ecorregión

Los datos anteriormente calulcados se deben ejecutar para cada ecorregión y sus tres respectivas equidistancias por estado, sin embargo se realizan de manera automatica en excel. Lo unico que se debe de hacer es tomar la tabla Conteos_Jalisco_2000_2022 y vaciar los datos en la matriz conteo de muestras, como se ha mencionado anteriormente.

Podemos calcular la superficie anual de deforestación bruta por ecorregión y por tipo de transición en hectareas. Ya que la transición de tierra más común es la de tierra forestal a praderas (FL-GL) se midió esta observación de la matriz para la ecorregión elevaciones semiaridas-mediterraneas. Los valores los obtuvimos anteriormente, estan en la matriz de áreas estimadas.

Elevaciones semiaridas-meridionales

\[ FL-GL = \sum_{i=2001}^{2022} (A_i 2.5\times2.5+A_i5\times5+A_i10\times 10) \] Donde:

\(A_i 2.5 \times 2.5\) es el área de deforestación en el año \(i\) para la equidistancia \(2.5 \times 2.5\)
\(A_i 5\times 5\) es el área de deforestación en el año \(i\) para la equidistancia \(5 \times 5\)
\(A_i 10\times 10\) es el área deforestada en el año \(i\) para la equidistancia \(10 \times 10\)

Deforestacion por transicion FL-GL para Elevaciones semiaridas-meridionales
	2001	2002	2003	2004	2005	2006	2007	2008	2009	2010	2011	2012	2013	2014	2015	2016	2017	2018	2019	2020	2021	2022
FLGL	0	2586.42	2586.42	2586.42	5172.85	0	10219.03	2586.42	0	2586.42	12805.45	10219.03	0	2586.42	13449.92	2586.42	0	0	2586.42	10219.03	2586.42	0

La sumatoria de las áreas deforestadas para cada intensidad de muestreo de la ecorregión Elevaciones Semiaridas-Meridionales es solo una de las cuatro observaciones para cada transción de uso de suelo, ya que el estado de Jalisco cuenta con 4 ecorregiones nivel 1. Se deben obtener los mismos valores para cada ecorregión faltante:

Selvas Calido-Humedas
Selvas Calido-Secas
Sierras Templadas

FL-GL

Deforestacion por transicion FL-GL para ESM
	2001	2002	2003	2004	2005	2006	2007	2008	2009	2010	2011	2012	2013	2014	2015	2016	2017	2018	2019	2020	2021	2022
Selvas Calido-Humedas	0.00	0.00	1875.78	0.00	625.26	625.26	0.00	0.00	0.00	0.00	625.26	0.00	0.00	0.00	0.00	0.00	0.00	625.26	0.00	0.00	0.00	0.00
Selvas Calido-Secas	0.00	0.00	0.00	629.13	4455.49	0.00	5765.31	3197.22	7704.26	3197.22	3145.67	2516.54	629.13	0.00	16479.92	18367.33	5765.31	5765.31	24184.18	1258.27	3197.22	629.13
Sierras Templadas	635.13	5716.17	5081.04	2540.52	7621.55	9516.46	3165.17	8235.73	9495.51	5705.69	8860.38	8235.73	2540.52	6340.82	5705.69	3175.65	3810.78	0.00	8940.87	6975.95	6351.29	1270.26

Se debe repetir el calculo para obtener la estimación de áreas deforestadas para cada transición. Por ejemplo, se estimó la transición FL-GS, sin embargo faltan las permanencias y las demas transiciones, el total de transiciones para estimar deforestación son:

Matriz de permanencia FL-FL
Transición	Ecorregiones	2001	2002	2003	2004	2005	2006	2007	2008	2009	2010	2011	2012	2013	2014	2015	2016	2017	2018	2019	2020	2021	2022
FL-FL	Elevaciones Semiaridas-Meridionales	601570.47	598984.05	596397.63	590580.32	585407.47	585407.47	575188.44	572602.02	570015.59	567429.17	554623.7	544404.7	544404.7	541818.3	528368.3	525781.9	525781.9	515562.89	512976.47	502757.44	500171.01	500171.01
	Selvas Cálido_Húmedas	93794.16	93794.16	91918.38	91918.38	91293.12	90667.86	90667.86	90667.86	90667.86	90667.86	90042.6	90042.6	90042.6	90042.6	90042.6	90042.6	90042.6	89417.34	89417.34	89417.34	89417.34	89417.34
	Selvas Cálido-Secas	1096173.73	1096173.73	1096173.73	1095544.59	1085952.93	1085952.93	1080187.62	1076990.40	1069286.13	1066088.91	1062943.2	1060426.7	1057229.5	1057229.5	1040749.6	1022382.2	1014048.8	1008283.52	994184.82	992926.55	989100.19	988471.06
	Sierras Templadas	255742.98	2550026.82	2542415.74	2539875.23	2532253.67	2522102.08	2517666.65	2507525.54	2498030.03	2490418.95	2480288.3	2472052.6	2468241.8	2461265.9	2455560.2	2451749.4	2445398.1	2444762.98	2473092.37	2429481.30	2422494.87	2418059.45

Matriz de deforestación para la transición FL-GL
Transición	Ecorregiones	2001	2002	2003	2004	2005	2006	2007	2008	2009	2010	2011	2012	2013	2014	2015	2016	2017	2018	2019	2020	2021	2022
FL-GL	Elevaciones Semiaridas-Meridionales	0.00	2586.42	2586.42	2586.42	5172.85	0.00	10219.03	2586.42	0.00	2586.42	12805.45	10219.03	0.00	2586.42	13449.92	2586.42	0.00	0.00	2586.42	10219.03	2586.42	0.00
	selvas Cálido-Húmedas	0.00	0.00	1875.78	0.00	625.26	625.26	0.00	0.00	0.00	0.00	625.26	0.00	0.00	0.00	0.00	0.00	0.00	625.26	0.00	0.00	0.00	0.00
	Selvas Cálido-Secas	0.00	0.00	0.00	629.13	4455.49	0.00	5765.31	3197.22	7704.26	3197.22	3145.67	2516.54	629.13	0.00	16479.92	18367.33	5765.31	5765.31	24184.18	1258.27	3197.22	629.13
	Sierras Templadas	635.13	5716.17	5081.04	2540.52	7621.55	9516.46	3165.17	8235.73	9495.51	5705.69	8860.38	8235.73	2540.52	6340.82	5705.69	3175.65	3810.78	0.00	8940.87	6975.95	6351.29	1270.26

Matriz de deforestación para la transición FL-CL
Transición	Ecorregiones	2003	2004	2005	2006	2007	2008	2009	2010	2011	2013	2014	2016	2017	2018	2019	2020	2021	2022
FL-CL	Elevaciones Semiaridas-Meridionales	0.00	3230.89	0.00	0.00	0.00	0.00	2586.42	0.00	0.00	0.00	0.00	0.00	0.00	10219.03	0.00	0.00	0.00	0.00
	Selvas Cálido-Húmedas	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00
	Selvas Cálido-Secas	0.00	0.00	5136.17	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	0.00	2568.09	0.00	0.00	0.00	629.13	0.00
	Sierras Templadas	2530.04	0.00	0.00	635.13	1270.26	1905.39	0.00	1905.39	1270.26	1270.26	635.13	1270.26	2540.52	635.13	1905.39	635.13	0.00	3165.17

Matriz de deforestación para la transición FL-WL
Transición	Ecorregiones	2013
FL-WL	Elevaciones Semiaridas-Meridionales	0.00
	Selvas Cálido-Húmedas	0.00
	Selvas Cálido-Secas	2568.09
	Sierras Templadas	0.00

Matriz de deforestación para la transición FL-SL
Transición	Ecorregiones	2001	2002	2003	2004	2005	2006	2007	2008	2009	2010	2011	2012	2013	2014	2015	2016	2017	2018	2019	2020	2021	2022
FL-SL	Elevaciones Semiaridas-Meridionales	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
	Selvas Cálido-Húmedas	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
	Selvas Cálido-Secas	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
	Sierras Templadas	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0

Matriz de deforestación para la transición FL-OL
Transición	Ecorregiones	2021
FL-OL	Elevaciones Semiaridas-Meridionales	0.00
	Selvas Cálido-Húmedas	0.00
	Selvas Cálido-Secas	0.00
	Sierras Templadas	635.13

Superficie anual de deforestación bruta por transición en hectáreas 2000-2022

Anteriormente se estimó la deforestación para cada ecorregión y para cada transición de uso de suelo en la serie temporal con la sumatoria de las áreas estimadas según su respectiva intensidades de muestreo (equidistancia). Una vez obtenidos las sumatorias de los valores por equidistancia se obtuvierón los valores totales por ecorregión. Anteriormente estimamos unicamente la transición FL-GL para Elevaciones Semiaridas-Meridionales, y despues para las demas ecorregiones. Con esos datos obtenidos, se procedió a estimar la misma transición, pero total, sumando las áreas estimadas de las cuatro ecorregiones. Primero, se analizó unicamente la transición FL-GL para una mejor visualización y entendimiento. Después, se estimo la deforestación total por transición y eventualmente por ecorregión. \[ FL-GL = \sum_{i=2001}^{2022} (A_iES+A_iSCH+A_iSCS+A_iST) \] Donde:

\(A_iES\) es el área deforestada en el año \(i\) para la ecorregión Elevaciones Semiaridas-Meridionales
\(A_iSCH\) es el área deforestada en el año \(i\) para la ecorregión Selvas Cálido-Húmedas
\(A_iSCS\) es el área deforestada en el año \(i\) para la ecorregión Selvas Cálido-Secas
\(A_iST\) es el área deforestada en el año \(i\) para la ecorregión Sierras Templadas

Deforestacion total para la transición FL-GL
	2001	2002	2003	2004	2005	2006	2007	2008	2009	2010	2011	2012	2013	2014	2015	2016	2017	2018	2019	2020	2021	2022
FLGL	635.13	8302.59	9543.24	5756.08	17875.15	10141.72	19149.51	14019.37	17199.77	11489.33	25436.77	20971.3	3169.65	8927.24	35635.53	24129.4	9576.09	6390.57	35711.47	18453.25	12134.94	1899.39

\[ \text{Deforestación total por transición}_i = \sum_{t=2001}^{2022} \left( \text{ESM}_{i,t} + \text{SCH}_{i,t} + \text{SCS}_{i,t} + \text{ST}_{i,t} \right) \]

Donde:

Deforestación Total por transición: Es la suma total de la deforestación en la transición \(i\) (por ejemplo FL-FL, FL-GL, FL-CL, FL-WL, FL-SL, FL-OL) a lo largo de los años \(t\) (de 2001 a 2022).
Ecorregiones: \(ESM_{it}\), \(SCH_{it}\), \(SCS_{it}\), \(ST_{it}\) son las áreas deforestadas en cada ecorregión para la transición \(i\) y el año \(t\)
Indices
- \(i\): Transición (FL-FL,FL-GS etc)
- \(t\): Año en la serie temporal

Superficie anual de la deforestación bruta por transición en hectáreas 2000-2022
Conversión	2001	2002	2003	2004	2005	2006	2007	2008	2009	2010	2011	2012	2013	2014	2015	2016	2017	2018	2019	2020	2021	2022
FL-FL	4347281	4338979	4326905	4317919	4294907	4284130	4263711	4247786	4228000	4214605	4187898	4166927	4159919	4150356	4114721	4089956	4075271	4058027	4033671	4014583	4001183	3996119
FL-GL	635	8303	9543	5756	17875	10142	19150	14019	17200	11489	25473	20971	3170	8927	35636	24129	9576	6391	35711	18453	12135	1899
FL-CL	0	0	2530	3231	5136	635	1270	1905	2586	1905	1270	0	1270	635	0	1270	5109	10854	1905	635	629	3165
FL-WL	0	0	0	0	0	0	0	0	0	0	0	0	2568	0	0	0	0	0	0	0	0	0
FL-SL	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
FL-OL	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	635	0

Deforestación por Ecorregión:

\[ \text{Deforestación Total}_{\text{Eco}} = \sum_{t=2001}^{2022} \sum_{i=1}^{n} A_{i, \text{Eco}, t} \] Donde

\(A_{i, \text{Eco}, t}\) es el área deforestada en la transición \(i\) dentro de la ecorregioón en el año \(t\)
\(i\) representa cada una de las transiciones posibles (FL-GL, FL-CL, FL-WL, etc).
\(t\) cubre los años de la serie temporal de 2001 a 2002

Superficie anual de la deforestación bruta or ecorregión en hecráreas 2000-2022
Conversión	2001	2002	2003	2004	2005	2006	2007	2008	2009	2010	2011	2012	2013	2014	2015	2016	2017	2018	2019	2020	2021	2022
Elevaciones Semiaridas-Meridionales	0	2586	2586	5817	5173	0	10219	2585	2586	2586	12805	10219	0	2586	13450	2586	0	10219	2586	10219	2585	0
Selvas Cálido Húmedas	0	0	1876	0	625	625	0	0	0	0	625	0	0	0	0	0	0	625	0	0	0	0
Selvas Cálido-Secas	0	0	0	629	9592	0	5765	3197	7704	3197	3146	2517	3197	0	16480	18367	8333	5765	24184	1258	3826	629
Sierras Templadas	635	5716	7611	2541	7622	10152	4435	10141	9496	7611	10131	8236	3811	6976	5706	4446	6351	635	10846	7611	6986	4435

(Espacio para agregar los datos siguientes): es imprecinsibles que agruegue los indicadores, como la tasa de deforestación y los indicadores de la inforgrafía de CONAFOR.

Propagación de incertidumbres

\[ \text{Incertidumbre}_{\text{i,j}} = \sum_{r=1}^{m} \sum_{k=1}^{n} (\frac{U_{\text{i,j,k,r}}}{\triangle A_{\text{i,j,r}}} ) \] Donde

\(i\) es la transición (FL-FL, FL-GL, FL-CL, FL-WL, FL-SL, FL-OL)
\(J\) es el año de la serie temporal \(2001 \le j \le 2002\)
\(k\) representa cada intensidad de muestro usada para la matriz de incertidumbres ajustada \((A \times U)^2\)
\(r\) es la ecorregión en cuestión (ESM, SCH, SCS, ST)
\(U_{\text{i,j,k,r}}\) es el valor de incertidumbre ajustada de la trnaisción \(i\) para el año \(j\), la intensidad de muestreo \(k\) y la ecorregión \(r\)
\(\mid \triangle A_{\text{i,j,r}} \mid\) es el valor absoulto de la deforestación total calculada para la transición \(i\) en el año \(j\) y ecorregión \(r\)

Esta matriz es un insumo indispensable para la segunda etapa de propagación de incertidumbres a travéz del metodo 1 del IPCC. Como se mencionó en la formula anterior, para obtener el valor de cada observación se sumó cada valor nivel equidistancia-ecorregión de la matriz obtnida anteriormente: Matriz de incertidumbre ajustada \((A \times U)^2\) El método de este calculo es igual al que se utilizó anteriormente para obtener las estimaciones de áreas deforestadas por transición. Se contabilizo una sumatoria de tres valores, por ejemplo:

la transición FL-GL \(2.5 \times 2.5\) + FL-GL \(5 \times 5\) + FL-GL \(10 \times 10\) para Elevaciones Semiaridas Meridionales para cada año
la transición FL-GL \(2.5 \times 2.5\) + FL-GL \(5 \times 5\) + FL-GL \(10 \times 10\) para Selvas Calido Humedas para cada año

Se utilizó el mismo metodo para las demas ecorregiones, y despues para las demas transiciones hasta que se obtuvó la siguiente matriz.

Matriz de incertidumbres en porcentajes FL-FL
Transición	Ecorregiones	2001	2002	2003	2004	2005	2006	2007	2008	2009	2010	2011	2012	2013	2014	2015	2016	2017	2018	2019	2020	2021	2022
FL-FL	Elevaciones Semiaridas-Meridionales	16.58	16.66	16.73	16.91	17.06	17.06	17.11	17.19	17.27	17.35	17.49	17.53	17.53	17.62	17.77	17.86	17.86	17.89	17.98	18.01	18.11	18.11
	Selvas Cálido Húmedas	4.67	4.67	5.19	5.19	5.36	5.52	5.52	5.52	5.52	5.52	5.68	5.68	5.68	5.68	5.68	5.68	5.68	5.84	5.84	5.84	5.84	5.84
	Selvas Cálido-Secas	6.65	6.65	6.65	6.66	6.75	6.75	6.81	6.84	6.92	6.95	6.97	6.99	7.03	7.03	7.05	7.07	7.16	7.23	7.38	7.39	7.43	7.43
	Sierras Tmepladas	3.60	3.61	3.63	3.63	3.65	3.67	3.68	3.70	3.72	3.73	3.76	3.77	3.78	3.79	3.81	3.81	3.83	3.83	3.81	3.83	3.84	3.85

Matriz de incertidumbres en porcentajes FL-GL
Transición	Ecorregiones	2001	2002	2003	2004	2005	2006	2007	2008	2009	2010	2011	2012	2013	2014	2015	2016	2017	2018	2019	2020	2021	2022
FL-GL	Elevaciones Semiaridas-Meridionales	0	196.00	196.00	196.00	138.12	0.00	196.00	196.00	0.00	196.00	161.34	196.00	0.00	196.00	153.90	196.00	0.00	0.00	196.00	196.00	196.00	0.00
	Selvas Cálido Húmedas	0	0.00	112.39	0.00	196.00	196.00	0.00	0.00	0.00	0.00	196.00	0.00	0.00	0.00	0.00	0.00	0.00	196.00	0.00	0.00	0.00	0.00
	Selvas Cálido-Secas	0	0.00	0.00	196.00	122.68	0.00	125.12	162.09	112.82	162.09	87.27	97.68	196.00	0.00	121.14	109.31	125.12	125.12	89.95	138.44	162.09	196.00
	Sierras Templadas	196	65.26	69.23	97.96	56.49	67.77	161.53	91.36	92.72	99.65	98.36	91.36	97.96	91.79	99.65	87.60	79.96	0.00	182.61	85.31	61.90	138.57

Matriz de incertidumbres en porcentajes FL-CL
Transición	Ecorregiones	2003	2004	2005	2006	2007	2008	2009	2010	2011	2013	2014	2016	2017	2018	2019	2020	2021	2022
FL-CL	Elevaciones Semiaridas-Meridionales	0	161.7	0.00	0	0.00	0.00	196	0.00	0.00	0.00	0	0.00	0.00	196	0.00	0	0	0.00
	Selvas Cálido-Húmedas	0	0.0	0.00	0	0.00	0.00	0	0.00	0.00	0.00	0	0.00	0.00	0	0.00	0	0	0.00
	Selvas Cálido Secas	0	0.0	138.38	0	0.00	0.00	0	0.00	0.00	0.00	0	0.00	196.00	0	0.00	0	196	0.00
	Sierras Templadas	196	0.0	0.00	196	138.57	113.13	0	113.13	138.57	138.57	196	138.57	97.96	196	113.13	196	0	161.53

Matriz de incertidumbres en porcentajes FL-WL
Transición	Ecorregiones	2013
FL-WL	Elevaciones Semiaridas-Meridionales	0
	Selvas Cálido-Húmedas	0
	Selvas Cálido-Secas	196
	Sierras Templadas	0

Matriz de incertidumbres en porcentajes FL-SL
Transición	Ecorregiones	2001	2002	2003	2004	2005	2006	2007	2008	2009	2010	2011	2012	2013	2014	2015	2016	2017	2018	2019	2020	2021	2022
FL-SL	Elevaciones Semiaridas-Meridionales	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
	Selvas Cálido-Húmedas	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
	Selvas Cálido-Secas	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
	Sierras Templadas	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0

Matriz de incertidumbres en porcentajes FL-OL
Transición	Ecorregiones	2021
FL-OL	Elevaciones Semiaridas-Meridionales	0
	Selvas Cálido-Húmedas	0
	Selvas Cálido-Secas	0
	Sierras Templadas	196

La matriz obtuvida es el insumo para poder calcular las incertidumbres ajustadas ya que es la \(U\) que multiplicara el área \(A\) y se elevará al cuadrado.

Las incertidumbres por transición y ecorregión anuales fueron propagadas según la metodología de la primera etapa, esto quiere decir que la tabla anterior se obtuvo mediante la sumatoria de los valores de incertidumbre ajustada por Ecorregión-Equidistancia para cada transición, como se mencionó anteriormente. Estas incertidumbres permitierón la propagación de la matriz de incertidumbre ajustada anualizada ya que se puede estimar el producto de el área total deforestada por la incertidumbre al cuadrado:

Matriz de incertidumbres ajustadas FL-FL
Transición	Ecorregiones	2001	2002	2003	2004	2005	2006	2007	2008	2009	2010	2011	2012	2013	2014	2015	2016	2017	2018	2019	2020	2021	2022
FL-FL	Elevaciones Semiaridas-Meridionales	9.951507e+13	9.956436e+13	9.961012e+13	9.969548e+13	9.976236e+13	9.976236e+13	9.689684e+13	9.692501e+13	9.694965e+13	9.697077e+13	9.407700e+13	9.111979e+13	9.111979e+13	9.113387e+13	8.814563e+13	8.815267e+13	8.815267e+13	8.510375e+13	8.510727e+13	8.201251e+13	8.201251e+13	8.201251e+13
	Selvas Cálido-Húmedas	1.916681e+11	1.916681e+11	2.275292e+11	2.275292e+11	2.390743e+11	2.504150e+11	2.504150e+11	2.504150e+11	2.504150e+11	2.504150e+11	2.615514e+11	2.615514e+11	2.615514e+11	2.615514e+11	2.615514e+11	2.615514e+11	2.615514e+11	2.724834e+11	2.724834e+11	2.724834e+11	2.724834e+11	2.724834e+11
	Selvas Cálido-Secas	5.317430e+13	5.317430e+13	5.317430e+13	5.318303e+13	5.373252e+13	5.373252e+13	5.408239e+13	5.425918e+13	5.475508e+13	5.492562e+13	5.496627e+13	5.499761e+13	5.516595e+13	5.516595e+13	5.378762e+13	5.230944e+13	5.277064e+13	5.307264e+13	5.379696e+13	5.381038e+13	5.396103e+13	5.396751e+13
	Sierras Templadas	8.485299e+13	8.495671e+13	8.516063e+13	8.520622e+13	8.534212e+13	8.558604e+13	8.572576e+13	8.601799e+13	8.633349e+13	8.650958e+13	8.679930e+13	8.700262e+13	8.706781e+13	8.721328e+13	8.733384e+13	8.739807e+13	8.750441e+13	8.751499e+13	8.642315e+13	8.657134e+13	8.668651e+13	8.677845e+13

Matriz de incertidumbres ajustadas FL-GL
Transición	Ecorregiones	2001	2002	2003	2004	2005	2006	2007	2008	2009	2010	2011	2012	2013	2014	2015	2016	2017	2018	2019	2020	2021	2022
FL-GL	Elevaciones Semiaridas-Meridionales	0	256987036525	256987036525	256987036525	510453702686	0	4.011729e+12	256987036525	0	256987036525	4.268716e+12	4.011729e+12	0	256987036525	4.284671e+12	2.569870e+11	0	0	2.569870e+11	4.011729e+12	256987036525	0
	Selvas Cálido-Húmedas	0	0	44443297043	0	15018769346	15018769346	0.000000e+00	0	0	0	1.501877e+10	0.000000e+00	0	0	0.000000e+00	0.000000e+00	0	15018769346	0.000000e+00	0.000000e+00	0	0
	Selvas Cálido Secas	0	0	0	15205463544	298774296967	0	5.203731e+11	268561702012	755434150065	268561702012	7.536621e+10	6.042519e+10	15205463544	0	3.985838e+12	4.030661e+12	520373085367	520373085367	4.732003e+12	3.034482e+10	268561702012	15205463544
	Sierras Templadas	15496609543	139150333742	123724646898	61933246147	185374202252	415879749111	2.614012e+11	566167809998	775069566516	323299032032	7.595907e+11	5.661678e+11	61933246147	338751314888	3.232990e+11	7.739439e+10	92846677197	0	2.665634e+12	3.541947e+11	154567155250	30984353748

Matriz de incertidumbres ajustadas FL-CL
Transición	Ecorregiones	2003	2004	2005	2006	2007	2008	2009	2010	2011	2013	2014	2016	2017	2018	2019	2020	2021	2022
FL-CL	Elevaciones Semiaridas-Meridionales	0	272942651900	0	0	0	0	2.56987e+11	0	0	0	0	0	0	4.011729e+12	0	0	0	0
	Selvas Cálido-Húmedas	0	0	0	0	0	0	0.00000e+00	0	0	0	0	0	0	0.000000e+00	0	0	0	0
	Selvas Cálido-Secas	0	0	505167621823	0	0	0	0.00000e+00	0	0	0	0	0	253356238468	0.000000e+00	0	0	15205463544	0
	Sierras Templadas	245904637691	0	0	15496609543	30984353748	46463232616	0.00000e+00	46463232616	30984353748	30984353748	15496609543	30984353748	61933246147	1.549661e+10	46463232616	15496609543	0	261401247234

Matriz de incertidumbres ajustadas FL-wL
Transición	Ecorregiones	2013
FL-WL	Elevaciones Semiaridas-Meridionales	0
	Selvas Cálido-Húmedas	0
	Selvas Cálido-Secas	253356238468
	Sierras Templadas	0

Matriz de incertidumbres ajustadas FL-SL
Transición	Ecorregiones	2001	2002	2003	2004	2005	2006	2007	2008	2009	2010	2011	2012	2013	2014	2015	2016	2017	2018	2019	2020	2021	2022
FL-SL	Elevaciones Semiaridas-Meridionales	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
	Selvas Cálido-Húmedas	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
	selvas Cálido-Secas	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
	Sierras Templadas	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0

Matriz de incertidumbres ajustadas FL-oL
Transición	Ecorregiones	2021
FL-OL	Elevacione Semiaridas-Meridionales	0
	Selvas Cálido-Húmedas	0
	Selvas Cálido-Secas	0
	Sierras Templadas	15496609543

Segunda etapa de poropagación de incertidumbres

Debido a que el área total deforestada es la suma de las áreas de cambio por Ecorregión-Equidistancia para las transiciones: “Tierra Forestal a Pradera”, “Tierra Forestal a Agricultura”, “Tierra Forestal a Asentamiento Humano”, “TierraForestal a Cuerpo de Agua” y “Tierra Forestal a Otras Tierras”, sus incertidumbres se propagaron a través del método 1 del IPCC por adición de aacuerdo a lo mostrado a continuación:

\(U(AD_T)\) es la incertidumbre del área total deforestada y expresada como porcentaje.
\(U(\widehat{AD}_i)_i \widehat{AD}_i\) es el porcentaje de incertidumbre asociado con el área deforestada estimafa de Ecorregión-Equidistancia y el área deforestada respectivamente

Incertidumbres totales anualizadas por transición

Incertidumbres en porcentajes totales por transición de uso de suelo 2000-2022
Conversión	2001	2002	2003	2004	2005	2006	2007	2008	2009	2010	2011	2012	2013	2014	2015	2016	2017	2018	2019	2020	2021	2022
FL-FL	4	4	4	4	4	4	4	4	4	4	4	4	4	4	4	4	4	4	4	4	4	0
FL-GL	196	76	68	100	56	65	114	75	72	80	89	103	88	86	82	87	82	114	77	114	68	113
FL-CL	0	0	196	162	138	196	139	113	196	113	139	0	139	196	0	139	110	185	113	196	196	0
FL-WL	0	0	0	0	0	0	0	0	0	0	0	0	196	0	0	0	0	0	0	0	0	0
FL-SL	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0
FL-OL	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	0	196	0

Incertidumbres totales anualzadas por ecorregión

Incertidumbres en porcentajes totales por Ecorregión de uso de suelo 2000-2022
Ecorregiones	2001	2002	2003	2004	2005	2006	2007	2008	2009	2010	2011	2012	2013	2014	2015	2016	2017	2018	2019	2020	2021	2022
Elevaciones Semiaridas-Meridionales	0	196	196	125	138	0	196	196	196	196	161	196	0	196	154	196	0	196	196	196	196	0
Selvas Cálido-Húmedas	0	0	112	0	196	196	0	0	0	0	196	0	0	0	0	0	0	196	0	0	0	0
Selvas Cálido-Secas	0	0	0	196	93	0	125	162	113	162	87	98	162	0	121	109	106	125	90	138	139	196
Sierras Templadas	196	65	80	98	56	65	122	77	93	80	88	91	80	85	100	74	62	196	152	80	59	122

Estimación de Datos de Actividad

Jesús Arturo Fragozo Quesada

2024-09-24