8.1. FENOLOGÍA DE FLORACIÓN 15 DE NOVIEMBRE DE 2012 Impartida por: José Miguel Muñoz Díaz
8.1. FENOLOGÍA DE FLORACIÓN
15 DE NOVIEMBRE DE 2012
Impartida por: José Miguel Muñoz Díaz
2009 (Part.)
2010
2011
2012
8.1. Fenología de floración
Captura de pantalla la base de datos y de Linaria Migrados los datos
de 2009 y 2010
Incluidos los datos de 2010 y 2011. Parte de 2012
Corrección de errores
Consultas. Propias de linaria o globales (.csv)
8.1. Fenología de floración
Análisis exploratorio de los resultados.
1.- Análisis previos. Regresión lineal
2.- Regresión no lineal. “Modelos de floración”
3. - Comentarios de los modelos
8.1. Fenología de floración
1.- Análisis previos. Regresión lineal.
8.1. Fenología de floración
Evolución de la floración a partir del valor medio del nivel de floración de la media
0 1 2 3 4 5
05
1015
20
Eje
mpl
ar
DJ 119
0 1 2 3 4 5
05
10
15
20
Eje
mpla
r
DJ 126
0 1 2 3 4 5
05
1015
20
Eje
mpl
ar
DJ 133
0 1 2 3 4 5
05
10
15
20
Eje
mpla
r
DJ 140
0 1 2 3 4 5
05
1015
20
Eje
mpl
ar
DJ 147
0 1 2 3 4 5
05
10
15
20
Eje
mpla
r
DJ 153
0 1 2 3 4 5
05
10
15
20
Eje
mpla
r
DJ 161
0 1 2 3 4 5
05
10
15
20
Eje
mpla
r
DJ 168
0 1 2 3 4 5
05
1015
20
Eje
mpl
ar
DJ 175
8.1. Fenología de floración
Una revisión mas detallada de los datos induce a pensar que el proceso de floración sigue un modelo no lineal logístico.
Abla Superior 1
0 50 100 150 200
0
1
2
3
4
5
6
Anthyllis cytisoides
Dia
Fase d
e f
lora
ció
n
2010
Abla Inferior 1
0 50 100 150 200
0
1
2
3
4
5
6
Anthyllis cytisoides
Dia
Fase d
e f
lora
ció
n
2010
Abla Superior 2
0 50 100 150 200
0
1
2
3
4
5
6
Anthyllis cytisoides
Dia
Fase d
e f
lora
ció
n
2010
Abla Inferior 2
0 50 100 150 200
0
1
2
3
4
5
6
Anthyllis cytisoides
Dia
Fase d
e f
lora
ció
n
2010
Abla Superior 1
0 50 100 150 200
0
1
2
3
4
5
6
Anthyllis cytisoides
Dia
Fas
e de
flo
raci
ón 2010
Abla Inferior 1
0 50 100 150 200
0
1
2
3
4
5
6
Anthyllis cytisoides
Dia
Fas
e de
flo
raci
ón 2010
Abla Superior 2
0 50 100 150 200
0
1
2
3
4
5
6
Anthyllis cytisoides
Dia
Fas
e de
flo
raci
ón 2010
Abla Inferior 2
0 50 100 150 200
0
1
2
3
4
5
6
Anthyllis cytisoides
Dia
Fas
e de
flo
raci
ón 2010
8.1. Fenología de floración
media ~ theta1/(1 + exp(-(theta2 + theta3 * dia.juliano1)))
Theta1. Asíntota superior (teóricamente 5)
Theta 2. Pendiente o tasa de crecimiento
Theta 3. Asíntota inferior (teóricamente 0)
8.1. Fenología de floración
Algunos ejemplos…
8.1. Fenología de floración
San Juan Inferior 2
0 100 200 300
0
1
2
3
4
5
6
Adenocarpus decorticans
Dia
Fas
e de
flo
raci
ón 2009
San Juan Inferior 2
0 100 200 300
0
1
2
3
4
5
6
Adenocarpus decorticans
Dia
Fas
e de
flo
raci
ón 2010
San Juan Inferior 2
0 100 200 300
0
1
2
3
4
5
6
Adenocarpus decorticans
Dia
Fas
e de
flo
raci
ón 2011
San Juan Inferior 2
0 100 200 300
0
1
2
3
4
5
6
Adenocarpus decorticans
Dia
Fas
e de
flo
raci
ón 2012
8.1. Fenología de floración
Abla Superior 2
0 100 200 300
0
1
2
3
4
5
6
Cytisus galianoi
DiaF
ase
de f
lora
ción 2009
Abla Superior 2
0 100 200 300
0
1
2
3
4
5
6
Cytisus galianoi
Dia
Fas
e de
flo
raci
ón 2010
Abla Superior 2
0 100 200 300
0
1
2
3
4
5
6
Cytisus galianoi
Dia
Fas
e de
flo
raci
ón 2011
0 100 200 300
0
1
2
3
4
5
6
Dia
Fas
e de
flo
raci
ón
San Juan Inferior 2
0 50 100 150 200
0
1
2
3
4
5
6
Adenocarpus decorticans
Dia
Fase d
e f
lora
ció
n
2009
San Juan Inferior 2
0 50 100 150 200
0
1
2
3
4
5
6
Adenocarpus decorticans
Dia
Fase d
e f
lora
ció
n
2010
San Juan Inferior 2
0 50 100 150 200
0
1
2
3
4
5
6
Adenocarpus decorticans
Dia
Fase d
e f
lora
ció
n
2011
0 50 100 150 200
0
1
2
3
4
5
6
DiaF
ase d
e f
lora
ció
n
8.1. Fenología de floración
San Juan Superior 2
0 50 100 150 200
0
1
2
3
4
5
6
Crataegus monogyna
Dia
Fase d
e f
lora
ció
n
2009
San Juan Superior 2
0 50 100 150 200
0
1
2
3
4
5
6
Crataegus monogyna
Dia
Fase d
e f
lora
ció
n
2010
San Juan Superior 2
0 50 100 150 200
0
1
2
3
4
5
6
Crataegus monogyna
Dia
Fase d
e f
lora
ció
n
2011
0 50 100 150 200
0
1
2
3
4
5
6
DiaF
ase d
e f
lora
ció
n
Con cambios de ubicación….
8.1. Fenología de floración
Resumen de los modelos de floración obtenidos. Por ajuste de los parámetros.
8.1. Fenología de floración
Anthyllis cytisoides
ABLARép. Año Theta1 Theta2 Theta3
CSR1
2009 No se hizo2010 4.48075*** -15.94181*** 0.11696***2011 4.84892*** -7.77491*** 0.05776***2012 Sin floración
CSR2
2009 No se hizo2010 4.61640*** -14.01098*** 0.09779***2011 4.73355** -7.05971* 0.05147*2012 Sin floración
CSI1
2009 No se hizo2010 4.14223*** -20.43161*** 0.15136***2011 4.88752*** -6.08300** 0.04857** 2012 Sin floración
CSI2
2009 No se hizo2010 4.33586*** -18.99564*** 0.14618***2011 4.61865*** -7.70632*** 0.06250***2012 Sin floración
8.1. Fenología de floración
ABLA CAÑARRép. Año Theta1 Theta2 Theta3 Theta1 Theta2 Theta3
CSR1
2009 5.7735 n.s. -21.6422 ns 0.1432 ns Sin datos2010 4.61022*** -25.23121*** 0.16760*** 4.64131*** -14.40223*** 0.10204***2011 No ajuste 4.94900** -8.39735** 0.06656*2012
CSR2
2009 4.8256 -32.7858 0.2187 Sin datos2010 4.6250*** -28.8621*** 0.1912*** 4.782616*** -11.799478*** 0.083896***2011 7.58388 -12.99996* 0.08222* 5.24347* -9.48723** 0.06931*2012
CIR1
2009 5.68389 -10.83944 0.08010 No se hizo2010 6.53710*** -12.67827*** 0.08617*** 4.93209*** -14.39371*** 0.10048*** 2011 4.98032** -11.07177 . 0.08793 . 7.41326** -8.43689** 0.05980**2012
CIR2
2009 5.60125 -13.26167 0.09686 No se hizo
2010 5.642972*** -17.928060*** 0.123668*** 5.23564*** -6.43997*** 0.05167**
2011 6.57377*** -11.89379*** 0.08254*** 5.07796*** -9.79800** 0.08450**2012
Adenocarpus decorticans (I)
8.1. Fenología de floración
Adenocarpus decorticans (II)
SAN JUANRép. Año Theta1 Theta2 Theta3
CSR1
2009 No se hizo No se hizo No se hizo2010 5.136861*** -12.957544*** 0.100022***2011 3.7020*** -34.1285 0.33252012
CSR2
2009 No se hizo No se hizo No se hizo2010 4.88539*** -13.61169*** 0.11067***2011 3.87433*** -22.04447** 0.22000**2012
CIR1
2009 No se hizo2010 4.96033*** -11.24885*** 0.09826***2011 4.38163 -13.59359 0.141052012
CIR2
2009 No se hizo No se hizo No se hizo2010 5.14877*** -9.83189*** 0.08836***2011 4.3643*** -26.7786* 0.2831*2012
8.1. Fenología de floración
SAN JUAN
Rép. Año Theta1 Theta2 Theta3
CSR1
2009
2010 4.85661*** -38.41434** 0.26355 **
2011 4.84598** -29.32714* 0.21983*
2012
CSR2
2009
2010 5.223269*** -25.484059*** 0.178465***
2011 4.81232*** -28.99210*** 0.22811***
2012
CSR1
2009
2010 5.32406*** -23.56673** 0.17049**
2011 5.04225*** -24.06907*** 0.19394***
2012
CIR2
2009
2010 5.33825*** -20.17816** 0.14829**
2011 5.18689*** -20.18140*** 0.16520***
2012
SAN JERÓNIMO
CSR1
2009 5.05739*** -28.97249* 0.19305*
2010
2011
2012
CSR2
2009 5.0524*** -34.4306* 0.2297*
2010
2011
2012
CIR1
2009 5.08958*** -30.05999** 0.20783**
2010 5.46509** -25.58646* 0.17210*
2011 4.1717*** -41.5238 . 0.3135 .
2012
CIR2
2009 4.98740*** -47.33341* 0.31936*
2010 5.47983*** -19.00110*** 0.13207**
2011 5.02086*** -30.56436*** 0.23471***
2012
Crataegus monogyna
8.1. Fenología de floración
Cytisus galianoi (I)ABLA
Rép. Año Theta1 Theta2 Theta3
CSR1
2009 4.97772*** -25.63315*** 0.18826***
2010 5.24619*** -20.29844** 0.14294**
2011 4.51984*** -28.49832*** 0.21241***
2012
CSR2
2009 4.95893*** -22.84882** 0.16940**
2010 4.98376*** -20.88339*** 0.14979***
2011 4.33635*** -28.11848*** 0.21168**
2012
CIR1
2009 5.02693*** -23.48046*** 0.17665***
2010 5.26247*** -13.13865** 0.09827*
2011 5.00250*** -20.61016* 0.15713*
2012
CIR2
2009 5.07069*** -19.08239** 0.14438**
2010 6.08384* -12.08893* 0.08795*
2011 5.28256*** -15.12228*** 0.11607***
2012
CAÑAR
Rép. Año Theta1 Theta2 Theta3
CSR1
2009 5.004517*** -49.837822* 0.353238*
2010 4.95924*** -21.14482** 0.15201*
2011 3.95380*** -39.58103* 0.30606*
2012
CSR2
2009 5.003139*** -57.389085* 0.402984*
2010 4.62179*** -27.97591** 0.19704**
2011 4.43013*** -26.74930** 0.20006**
2012
CIR1
2009 5.001755*** -69.670144** 0.483854**
2010 4.93567*** -19.09863*** 0.13413***
2011 5.55779 n.s -12.12015 n.s 0.08604 n.s.
2012
CIR2
2009 Mal ajuste
2010 4.98452*** -13.54839*** 0.10047***
2011 4.96342*** -12.09469*** 0.09228***
2012
8.1. Fenología de floración
Cytisus galianoi (II)
CAMARATE
Rép. Año Theta1 Theta2 Theta3
CSR1
2009 5.0432*** -30.0211** 0.2047**
2010 4.59822*** -29.68105** 0.20113**
2011
2012
CSR2
2009 5.09682*** -40.83671** 0.25877**
2010 4.53177*** -38.61141* 0.25663*
2011
2012
CIR1
2009 4.97898*** -33.21519*** 0.23454***
2010
2011 4.73789*** -26.03322** 0.19338**
2012
CIR2
2009 4.81521*** -24.44921* 0.17599*
2010 5.14324* -18.77941* 0.13475*
2011 4.60579*** -23.66798* 0.17802*
2012
8.1. Fenología de floración
Genista umbellataABLA
Rép. Año Theta1 Theta2 Theta3
CSR1
2009
2010 4.39157*** -34.50815** 0.21459**
2011 14.92788 -16.58578* 0.09005*
2012
CSR2
2009
2010 4.49673*** -32.24309** 0.20600**
2011 5.19324** -19.34780* 0.12216*
2012
CIR1
2009
2010 5.25181*** -16.42529*** 0.11344***
2011 5.1171*** -15.4793** 0.1055**
2012
CIR2
2009
2010 4.63313*** -32.28507*** 0.21882***
2011 5.01547*** -19.94244** 0.13636**
2012
CAÑAR
Rép. Año Theta1 Theta2 Theta3
CSR1
2009
2010 4.54613*** -29.35322*** 0.17956***
2011 5.97106* -18.11249* 0.10716*
2012
CSR2
2009
2010 4.62318*** -23.74208*** 0.14265***
2011 15.58081 -11.36791** 0.05812
2012
CIR1
2009
2010 4.72375*** -21.62465*** 0.13417***
2011 6.11205* -11.55821* 0.07203*
2012
CIR2
2009
2010 4.47403*** -23.97408*** 0.14522***
2011 6.24313 . -13.00364* 0.07687*
2012
8.1. Fenología de floración
Genista versicolor (I)
ABLA
Rép. Año Theta1 Theta2 Theta3
CSR1
2009 5.45468** -28.21357* 0.17719*
2010 5.26179*** -22.90763*** 0.14107***
2011 5.31505** -24.60752** 0.15336**
2012
CSR2
2009 5.43803** -28.62317* 0.17968*
2010 5.20710*** -22.67529*** 0.14057***
2011 5.29181** -23.95327* 0.14868*
2012
CIR1
2009 5.65458*** -20.93700** 0.13276**
2010 4.88599*** -25.71624** 0.16694**
2011 6.30073** -19.45700** 0.12006**
2012
CIR2
2009 5.08127*** -29.77492** 0.19896*
2010 4.88329*** -27.21640** 0.17863**
2011 5.1387** -24.9884* 0.1640*
2012
CAÑAR
Rép. Año Theta1 Theta2 Theta3
CSR1
2009 3.5972n.s. -32.1257n.s. 0.2031n.s.
2010 5.52116** -16.46200* 0.09518*
2011 4.12205. -22.48301. 0.13663.
2012
CSR2
2009 2.04407n.s. -33.56170n.s. 0.22082n.s.
2010 5.19141*** -17.59715** 0.10417**
2011 5.66392 -18.20944 0.10982
2012
CIR1
2009 4.63715n.s. -16.78564n.s. 0.10833n.s.
2010 5.05054*** -14.32705** 0.08770**
2011 8.79179 -11.36297* 0.06533
2012
CRI2
2009 No modelo logístico
2010 4.68985** -19.40928* 0.12274*
2011 4.29695*** -22.20777** 0.14563**
2012
8.1. Fenología de floración
CAMARATE
Rép. Año Theta1 Theta2 Theta3
CSR1
2009 5.26167** -60.33965* 0.35173*
2010 4.83724** -19.87372* 0.11989*
2011 10.6694 n.s. -25.7973 n.s. 0.1426 n.s.
2012
CSR2
2009 11.66803n.s. -27.56731 n.s. 0.15132 n.s.
2010 4.45921** -26.17383* 0.15581*
2011 Sin ajuste logístico
2012
CIR1
2009 5.28614** -33.33914* 0.20824*
2010 4.56810** -25.83016* 0.15887*
2011 5.53217*** -19.91565*** 0.12279***
2012
CIR2
2009 10.73147 n.s. -23.82990* 0.13670.
2010 4.37521*** -28.46899* 0.17820*
2011 5.80174 -17.15295 0.10454
2012
Genista versicolor (II)
8.1. Fenología de floración
Prunus avium
CAMARATE
Rép. Año Theta1 Theta2 Theta3
CSR1
2009
2010 4.068988** -56.780122** 0.482853**
2011 8.54614 -7.29547 0.06043
2012
CSR2
2009
2010 3.7617 -35.1122 0.3145
2011
2012
CIR1
2009
2010 5.4611 -13.7171 0.1244
2011 4.8416* -46.4035 0.4614
2012
CIR2
2009
2010 5.34208 . -13.36919 0.12412
2011 4.92714* -36.86035 . 0.36988 .
2012
8.1. Fenología de floración
SAN JUAN
Rép. Año Theta1 Theta2 Theta3
CSR1
2009
2010
2011 8.38517* -3.42119*** 0.04276*
2012
CSR2
2009
2010
2011 6.768436*** -5.155217*** 0.069196***
2012
CIR1
2009
2010
2011 5.44121*** -3.59079*** 0.06838***
2012
CIR2
2009
2010
2011 6.583183*** -3.192470*** 0.051907***
2012
Prunus dulcis
8.1. Fenología de floración
[JMMD1]Se queda la gráfica en 2
Retama sphaerocarpa
ABLA
Rép. Año Theta1 Theta2 Theta3
CSR1
2009
2010 4.73664*** -28.52959** 0.16270**
2011 2.23986 ** -82.66606* 0.50106*
2012
CSR2
2009
2010 5.3108*** -21.9996*** 0.1250**
2011 2.87993*** -57.64082** 0.34606**
2012
CIR1
2009
2010 3.6653*** -66.3598n.s. 0.4201n.s.
2011 5.84144n.s. -18.15712 n.s. 0.10855n.s.
2012
CSI2
2009
2010 4.67632*** -21.16414* 0.13065*
2011 4.99598n.s. -18.14388 n.s. 0.11220 n.s.
2012
[JMMD1]Se queda la gráfica en 2
8.1. Fenología de floración
Conclusiones
Los parámetros de los modelos tienen un alto nivel de significación para el año 2010
Los parámetros de los modelos tienen un bajo nivel de significación para el año 2011.
Esas alteraciones del modelo se producen en fases intermedias del proceso. Posiblemente se deba a cambios ambientales.
Es fácil estimar el inicio y final del proceso para establecer las series temporales. Además se pueden relacionar estos indicadores con otros factores como altitud, coordenadas geográficas, orientación, cobertura, etc….
8.1. Fenología de floración
Se han detectado cambios en la progresión del proceso para la misma réplica en años diferentes (cambios en theta 2).
El comportamiento de las especies ha sido muy dispar. Unas como Crataegus monogyna tienen un comportamiento muy regular. Por el contrario, otras como Anthyllis cytisoides ,son muy sensibles a factores ambientales. En esta última se puede dar el caso de que no florezca (como ha ocurrido 2012).
Los resultados son homologables a otras redes de seguimiento fenológico internacionales
8.1. Fenología de floración
8.1. Fenología de floración
Colaboran: