extremas en el Norte de Mozambique INTRODUCCIÓN La estimación de la magnitud

y

frecuencia de las láminas de precipitaciones máximas es un elemento clave en estudios

naturaleza , resultan raros , es decir , eventos con períodos de retorno de entre 10

y

500 años . Esto se debe fundamentalmente a la naturaleza destructiva de los mismos

y

al riesgo que suponen para la vida y las infraestructuras . Una estimación fiable

a la naturaleza destructiva de los mismos y al riesgo que suponen para la vida

y

las infraestructuras . Una estimación fiable de tales eventos requiere de series

la calidad requerida es uno de los principales retos en diferentes zonas del mundo

y

muy especialmente en regiones de África . La región de África del Sur es especialmente

frecuencia de los eventos extremos se basa exclusivamente en el uso de información local

y

presenta el inconveniente de la escasa longitud de las series de observaciones disponibles

carencia de datos en el tiempo con la abundancia en el espacio ” ( Hosking

y

Wallis , 1997 ) . El análisis regional es un método útil para caracterizar la

donde la información es escasa o simplemente no existe . A nivel práctico ,

y

a partir de los años 90 , la reformulación del análisis regional de frecuencia

extremos tomando como base los métodos del índice de avenida ( Dalrymple , 1960 )

y

los L - momentos ( Hosking , 1990 ; Hosking y Wallis , 1997 )

Dalrymple , 1960 ) y los L - momentos ( Hosking , 1990 ; Hosking

y

Wallis , 1997 ) se han empleado de forma habitual a nivel internacional con estudios

al . , 1997 ) , Italia ( Norbiato et al . , 2007 )

y

China ( Yang et al . , 2010 ) , entre otros . En la

2010 ) , entre otros . En la región de África del Sur , Parida

y

Moalafhi ( 2008 ) utilizaron el método de los L - momentos en el análisis

la General de Valores Extremos ( Jenkinson , 1955 ) . En Sudáfrica , Smithers

y

Schulze ( 2001 ) utilizaron los L - momentos para el análisis de tormentas de

aplicaron el análisis regional de frecuencia basado en el método del índice de avenida

y

los L - momentos a las series anuales de precipitaciones máximas de 23 estaciones

de 23 estaciones localizadas en la región meridional del país . En Mozambique ,

y

en particular en la región norte de país , no se han realizado estudios de

Algunas de las estimaciones realizadas hasta la fecha se centran en el análisis local

y

toman como base el ajuste de la función de distribución de Gumbel ( Gumbel ,

incalculables dimensiones que ha servido para agravar el proceso de recuperación económica

y

social en la que se encuentra inmerso el país . Se estima , por ejemplo ,

cuencas de los principales ríos del país fueron las peores de los últimos 150 años

y

costaron la vida a 800 personas y afectaron a otros 4 . 5 millones de

fueron las peores de los últimos 150 años y costaron la vida a 800 personas

y

afectaron a otros 4 . 5 millones de personas . Algunas estimaciones apuntan a que

, tomando como base los métodos del índice de avenida de Dalrymple ( 1960 )

y

los L - momentos ( Hosking , 1990 ; Hosking y Wallis , 1997 )

Dalrymple ( 1960 ) y los L - momentos ( Hosking , 1990 ; Hosking

y

Wallis , 1997 ) . La aplicación de ese método puede ser de utilidad práctica

método puede ser de utilidad práctica en estudios relacionados con las avenidas

y

la delimitación de zonas inundables en la región . MATERIAL Y MÉTODOS La región

con las avenidas y la delimitación de zonas inundables en la región . MATERIAL

Y

MÉTODOS La región de estudio La región objeto de este estudio está formada por

este estudio está formada por las cuencas internas de la provincia de Cabo Delgado

y

se localiza en la zona norte de Mozambique . Las principales cuencas son las de

de Mozambique . Las principales cuencas son las de los ríos Megaruma , Montepuez

y

Messalo , existiendo , además , otras cuencas de dimensiones menores que conforman

está marcada por toda la franja litoral de Cabo Delgado sobre el Océano Índico

y

discurre fundamentalmente en dirección norte - sur . La orientación de los ejes

principales de las cuencas que componen la región discurre en dirección oeste - este

y

vierten sus aguas en el Océano Índico . El clima de la región está

clima de la región está influenciado por el régimen de los monzones del Índico

y

por la corriente cálida del Canal de Mozambique . Presenta características de clima

Mozambique . Presenta características de clima templado en las regiones montañosas

y

tropical lluvioso en las restantes zonas . En toda la región se distinguen claramente

dos estaciones : la húmeda , de noviembre a marzo con vientos del Nordeste ,

y

la estación seca , más acentuada de abril a septiembre con temperaturas suaves

y

predominancia de vientos del Sudoeste . La precipitación anual alcanza los 1020

alcanza los 1020 mm . La temperatura media anual se sitúa en 25 . 2˚C

y

oscila entre 24 . 3˚C y los 26 . 1˚C . La humedad relativa es

temperatura media anual se sitúa en 25 . 2˚C y oscila entre 24 . 3˚C

y

los 26 . 1˚C . La humedad relativa es elevada y alcanza valores entre el

entre 24 . 3˚C y los 26 . 1˚C . La humedad relativa es elevada

y

alcanza valores entre el 70 y el 80 % . La población de las

26 . 1˚C . La humedad relativa es elevada y alcanza valores entre el 70

y

el 80 % . La población de las cuencas internas de Cabo Delgado en

expuestas a los desastres naturales relacionados con las precipitaciones intensas

y

las inundaciones . Disponibilidad de información Se han utilizado los registros

estaciones localizadas dentro de los límites de las cuencas internas de Cabo Delgado

y

sus cercanías . La disponibilidad de información pluviométrica en la región ,

región , cuya extensión supera los 5 . 6 · 104 km2 , resulta escasa

y

el acceso a la misma , difícil . La densidad de la red pluviométrica es

que se encuentran en el rango de una estación cada 250 km2 en zonas montañosas

y

cada 900 km2 para las zonas costeras . Los registros de las estaciones utilizadas

PRDA ) promovido por la Biblioteca Central de la Administración Nacional Oceánica

y

Atmosférica ( NOAA ) del Departamento de Comercio de los Estados Unidos . Dentro

anuarios de observaciones meteorológicas del Servicio Meteorológico de Mozambique

y

pueden ser descargados gratuitamente en formato digital ( http : / / docs . lib

series disponibles tienen registros generalmente cortos e incompletos de entre 20

y

34 años de observaciones , siendo la longitud de la serie promedio del orden de

28 años . La altitud de las estaciones oscila entre los 10 metros en Mecufi

y

847 metros sobre el nivel medio del mar en la meseta de Mueda .

partir de un conjunto de pruebas implementadas en el modelo estadístico TREND ( Chiew

y

Siriwardena , 2005 ; Grayson et al . , 1996 ) , que recoge las

, 2005 ; Grayson et al . , 1996 ) , que recoge las recomendaciones

y

métodos emanados del Taller de Expertos de la OMM / UNESCO ( Kundzewicz y Robson

recomendaciones y métodos emanados del Taller de Expertos de la OMM / UNESCO ( Kundzewicz

y

Robson , 2000 ) sobre pruebas estadísticas para la verificación de la calidad de

) las de Mann - Kendall ( Mann , 1945 ; Kendall , 1975 )

y

Spearman´s Rho ( Siegel y Castellan , 1988 ) para el análisis de tendencias ;

Kendall ( Mann , 1945 ; Kendall , 1975 ) y Spearman´s Rho ( Siegel

y

Castellan , 1988 ) para el análisis de tendencias ; ( II ) las de

II ) las de suma de rangos ( Hirsh et al . , 1982 )

y

la t de Student ( Student , 1908 ) para el análisis de la diferencia

1908 ) para el análisis de la diferencia de medianas / medias en dos períodos

y

, ( III ) las pruebas de cruce de la mediana ( Fisz , 1963

mediana ( Fisz , 1963 ) , diferencia de rangos ( Meacham , 1968 )

y

de auto - correlación ( Jenkins y Watts , 1968 ) para el análisis de

diferencia de rangos ( Meacham , 1968 ) y de auto - correlación ( Jenkins

y

Watts , 1968 ) para el análisis de la aleatoriedad . El nivel de significación

precipitaciones máximas se utilizó el algoritmo de los L - momentos de Hosking

y

Wallis ( 1997 ) y el método del índice de avenida ( Dalrymple , 1960

utilizó el algoritmo de los L - momentos de Hosking y Wallis ( 1997 )

y

el método del índice de avenida ( Dalrymple , 1960 ) . El uso de

, 1993 ) . Los L - momentos fueron introducidos por Hosking ( 1990 )

y

constituyen un sistema alternativo para describir las formas de las funciones de

alternativo para describir las formas de las funciones de distribución de probabilidades

y

se derivan de los momentos ponderados probabilísticamente ( MPP ) de Greenwood et

estimadores muestrales de los ratios de los L - momentos ( t , t3 ,

y

t4 ) se calculan reemplazando λ1 , λ2 , λ3 y λ4 por sus respectivos

t , t3 , y t4 ) se calculan reemplazando λ1 , λ2 , λ3

y

λ4 por sus respectivos estimadores muestrales l1 , l2 , l3 , y l4 en

, λ3 y λ4 por sus respectivos estimadores muestrales l1 , l2 , l3 ,

y

l4 en las Ecuaciones 3 - 6 . En Hosking y Wallis ( 1997 )

l2 , l3 , y l4 en las Ecuaciones 3 - 6 . En Hosking

y

Wallis ( 1997 ) se exponen detalladamente las ecuaciones para la estimación vía

muestras , la función de distribución subyacente ( Voguel y Fennessey , 1993 ; Hosking

y

Wallis , 1997 ) . Etapas del análisis regional de frecuencia El ARF

selección de la función de distribución regional de frecuencia de mejor ajuste

y

, d ) la estimación de los cuantiles de frecuencia . Estas etapas se describen

: la medida de discordancia D . La medida de discordancia D propuesta por Hosking

y

Wallis ( 1997 ) sirve para determinar la existencia de posibles estaciones “ discordantes

región . Los cambios de tendencias , así como los valores anómalos ( outliers )

y

demás errores muestrales de las series de observaciones se reflejan en sus respectivos

) tiene como objetivo evaluar el grado de heterogeneidad de la región propuesta

y

compara la variabilidad de los ratios de los L - momentos muestrales de las estaciones

los L - momentos muestrales de las estaciones de la región ( LCv , LCs

y

LCk ) con la variabilidad esperada en una región homogénea , obtenida a partir

una región homogénea , obtenida a partir de técnicas de simulación ( Hosking

y

Wallis , 1997 ) . La función de distribución Kappa - KAP ( Hosking ,

número elevado de M de simulaciones ( ~1000 ) de una región con N estaciones

y

número de registros por estación igual que la región original . Una región se

si Hn < 1 , “ posiblemente heterogénea ” si 1 < Hn < 2

y

“ definitivamente heterogénea ” en el caso que Hn ≥ 2 . Selección de

homogéneas corresponde seleccionar la función de distribución de mejor ajuste . Hosking

y

Wallis ( 1997 ) recomiendan la elección entre las siguientes cinco distribuciones

( GEV ) , Normal Generalizada ( GNO ) , Pareto Generalizada ( GPA )

y

la Pearson Tipo III ( PE3 ) . Para discernir entre las DIST candidatas se

aceptables , la de mejor ajuste será la que tenga el mínimo |ZDIST| ( Hosking

y

Wallis , 1997 ) y por tanto , la de mejor capacidad descriptiva . La obtención

ajuste será la que tenga el mínimo |ZDIST| ( Hosking y Wallis , 1997 )

y

por tanto , la de mejor capacidad descriptiva . La obtención de la medida ZDIST

descriptiva . La obtención de la medida ZDIST se basa en simulaciones de Monte Carlo

y

lleva el mismo proceso que el descrito en la etapa anterior para la medida de

estiman a partir de cada una de las distribuciones candidatas que resulten aceptables ,

y

de éstas , la DIST de mejor ajuste , siguiendo el método regional del índice

la ecuación anterior se define por el error medio cuadrático relativo ( RMSE )

y

se evalúa siguiendo las recomendaciones dadas por Hosking y Wallis ( 1997 ) , basadas

cuadrático relativo ( RMSE ) y se evalúa siguiendo las recomendaciones dadas por Hosking

y

Wallis ( 1997 ) , basadas en el empleo de simulaciones de Monte Carlo a

disponibles en la región de estudio ( Nº de estaciones , registros por estación

y

ratios medios regionales de los L - momentos muestrales ) y reproduciendo un número

registros por estación y ratios medios regionales de los L - momentos muestrales )

y

reproduciendo un número elevado de regiones artificiales ( M~1000 ) iguales que

cálculo seguidos en este estudio , se pueden encontrar con todo detalle en Hosking

y

Wallis ( 1997 ) . Desde el punto de vista computacional , la metodología del

metodología del ARF se encuentran implementada en los paquetes lmomRFA ( Hosking , 2014a )

y

lmom ( Hosking , 2014b ) que funcionan bajo R , entorno de programación para

2014b ) que funcionan bajo R , entorno de programación para análisis estadístico

y

gráfico ( http : / / www . r - project . org ) .

http : / / www . r - project . org ) . RESULTADOS

Y

DISCUSIÓN La calidad de las 12 series de precipitaciones máximas diarias empleadas

de regiones homogéneas . Dada la considerable extensión del área de la región

y

el limitado número de estaciones existente se consideró inicialmente una única

única región . Dentro de la región hay 7 estaciones que ocupan las partes media

y

alta de las cuencas internas de Cabo Delgado y están emplazadas entre los 200 y

estaciones que ocupan las partes media y alta de las cuencas internas de Cabo Delgado

y

están emplazadas entre los 200 y los 850 m de altitud , observándose , en

y alta de las cuencas internas de Cabo Delgado y están emplazadas entre los 200

y

los 850 m de altitud , observándose , en algunas de ellas , láminas precipitaciones

estaciones están localizadas en la franja litoral de la provincia de Cabo Delgado

y

tienen como característica común la altura a la que se encuentran emplazadas (

característica común la altura a la que se encuentran emplazadas ( ≤ 60 m )

y

láminas anuales de precipitación cercanas o inferiores en algunos casos a los 1000

. Para la región , fueron calculadas las respectivas medidas de discordancia Di

y

heterogeneidad Hn . Según se aprecia en la Tabla 3 , no existen a nivel

valores más elevados de la medida de discordancia recaen sobre las estaciones de Muite

y

Pemba , 2 . 14 y 2 . 51 respectivamente , pero quedan aún distantes

medida de discordancia recaen sobre las estaciones de Muite y Pemba , 2 . 14

y

2 . 51 respectivamente , pero quedan aún distantes del valor crítico de la región

función de distribución Kappa de una región con el mismo número de estaciones

y

mismas longitudes muestrales que la real . En la tabla 3 , se observa que

en la región de estudio , debido a la extensión del área de la misma

y

a la cercanía entre ciertas estaciones . Esta situación se refleja en el valor

correlación obtenido de la matriz de correlaciones de las 12 estaciones de la región

y

que es igual 0 . 171 . Por esta razón se recalcularon los correspondientes valores

cálculo muestra que su valor sube a 1 . 12 , mientras que para H2

y

H3 , los valores suben ligeramente a 0 . 39 y 0 . 24 respectivamente

mientras que para H2 y H3 , los valores suben ligeramente a 0 . 39

y

0 . 24 respectivamente . Puesto que H1 tiene un mayor peso en la evaluación

los cuantiles más precisos que aquellos obtenidos por el análisis local ( Hosking

y

Wallis , 1997 ) , especialmente en regiones con registros insuficientes , como es

5 distribuciones candidatas de tres parámetros ( GLO , GEV , GNO , GPA ,

y

PE3 ) , recomendadas por Hosking y Wallis ( 1997 ) . Los resultados de

GLO , GEV , GNO , GPA , y PE3 ) , recomendadas por Hosking

y

Wallis ( 1997 ) . Los resultados de la prueba se muestran en la siguiente

candidatas aceptadas según las simulaciones de Monte Carlo realizadas , fueron la GLO

y

la GEV respectivamente , dado que cumplían con el criterio ZDIST ≤ |1 . 64|

≤ |1 . 64| , mientras que las restantes tres distribuciones ( GNO , PE3

y

GPA ) fueron rechazadas por presentar un ZDIST >|1 . 64| . De las distribuciones

vecino con la región noroccidental de Mozambique , en las que las distribuciones GLO

y

GEV también mostraron el mejor ajuste a las series de máximos anuales de precipitaciones

diagrama de la siguiente figura se muestran las relaciones teóricas entre los LCs

y

LCk dadas por Hosking y Wallis ( 1997 ) para cada una de las funciones

figura se muestran las relaciones teóricas entre los LCs y LCk dadas por Hosking

y

Wallis ( 1997 ) para cada una de las funciones de distribución candidatas . En

máximas . Las distribuciones teóricas de tres parámetros se representan como líneas

y

las de 2 como puntos ( cuadrados negros ) . Clave de distribuciones : GLO

- Exponencial , G – Gumbel , L - Logística , N – Normal ,

y

U – Uniforme . El promedio regional de LCs y LCk se representa mediante un

, N – Normal , y U – Uniforme . El promedio regional de LCs

y

LCk se representa mediante un triángulo rojo . La función de distribución Logística

La función de distribución Logística Generalizada ( Hosking , 1986 ; Hosking

y

Wallis , 1997 ) se define por sus parámetros de posición ( u )

define por sus parámetros de posición ( u ) , escala ( α )

y

forma ( k ) estimados regionalmente a partir de los L - momentos y se

) y forma ( k ) estimados regionalmente a partir de los L - momentos

y

se expresa como : Los parámetros de la ecuación anterior se determinan en

ecuación anterior se determinan en función de los L - momentos muestrales ( Hosking

y

Wallis , 1997 ) En la siguiente tabla se da un resumen con la

de los cuantiles locales estimados para algunos períodos de retorno T de entre 10

y

500 años que pueden ser de interés práctico en estudios de hidrológicos de avenidas

años que pueden ser de interés práctico en estudios de hidrológicos de avenidas

y

delimitación de zonas inundables . Centrando el análisis de los resultados de

de los resultados de la tabla anterior , en los períodos de retorno de 100

y

500 años se observa que , para T igual a 100 años , la precipitación

máxima diaria en la región oscila entre los 237 mm en la estación de Quissanga

y

los 183 mm en la de Montepuéz , siendo el promedio en toda la región

años , en esas mismas estaciones , las láminas máximas estimadas alcanzan los 343

y

265 mm respectivamente , mientras que promedio regional sube a los 300 mm . Si

300 mm . Si se toma como referencia el período de retorno de 100 años

y

se comparan los resultados obtenidos en este estudio con los derivados por Ngongondo

. 4 ; Mangochi : 179 . 5 mm , Monkebay : 243 . 7

y

Ntaja : 219 . 4 mm ) , se observan resultados muy similares . El

de retorno de las estaciones de la región G2 se sitúa en los 208 mm

y

resulta coincidente con el obtenido en el presente estudio . Se destaca también

la citada región , los valores regionales de los coeficientes de L - asimetría

y

L - curtosis ( ~0 . 26 y 0 . 24 ) son cercanos a

de los coeficientes de L - asimetría y L - curtosis ( ~0 . 26

y

0 . 24 ) son cercanos a los obtenidos en la región de estudio (

) son cercanos a los obtenidos en la región de estudio ( 0 . 248

y

0 . 216 ) . CONCLUSIONES En este estudio se ha realizado el

de las precipitaciones máximas basado en el uso de datos locales exclusivamente

y

justifican la conveniencia del análisis regional como alternativa a la carencia

5 distribuciones candidatas de tres parámetros , de las que sólo dos , la GLO

y

la GEV , mostraron un ajuste aceptable . De ellas , se adoptó el modelo

los observados en el país vecino de Malawi , en las que las distribuciones GLO

y

GEV también mostraron el mejor ajuste . Para el período de retorno de 100

máxima diaria en la región oscila entre los 237 mm en la estación de Quissanga

y

los 183 mm en la de Montepuéz , siendo el promedio en toda la región

mismas estaciones las láminas máximas estimadas para 500 años alcanzan los 343

y

265 mm respectivamente , mientras que el promedio regional sube a los 300 mm .

Communication Technology ) . Estos sistemas están formados por centros de control

y

adquisición de datos ( SCADA , del inglés Supervisory Control and Data Acquisition

adquisición de datos ( SCADA , del inglés Supervisory Control and Data Acquisition )

y

sistemas de información geográfica ( SIG , del inglés Geographic Information System

recogida por el sistema no se traduce automáticamente en una mejora de la operación

y

el servicio al usuario . Es necesario transformar toda la información en conocimiento

de la red para que ésta sea útil para la toma de decisiones a corto

y

largo plazo , con la consiguiente mejora del funcionamiento del sistema . Con este

caracterizar el estado de grandes redes eléctricas ( Schweppe and Wildes , 1970 )

y

extendidas poco después a sistemas de abastecimiento de agua ( Coulbeck , 1977 )

permite deducir , en base a toda la información recopilada por el sistema SCADA

y

las ecuaciones que gobiernan el problema hidráulico , las condiciones de flujo en

gobiernan el problema hidráulico , las condiciones de flujo en cualquier instante

y

cualquier punto de la red . Concretamente , la EE se plantea de forma que

son el mínimo conjunto de variables que permite caracterizar el estado de la red ,

y

a partir de las cuales se determinarán el resto de variables del sistema gracias

lineal , x que resulta de la aplicación de las ecuaciones de conservación de masa

y

energía : donde ε representa los errores asociados con las medidas , típicamente

asociados con las medidas , típicamente asumidos como gaussianos con media cero

y

matriz varianza - covarianza Rz , la cual resulta ser una matriz diagonal puesto

en ese caso la matriz [ JT – 1J ( ν ) ] es invertible

y

el sistema resulta compatible determinado . Esta condición implica que no cualquier

, es un estudio previo que permite identificar qué variables ( flujos , demandas

y

presiones ) podrán conocerse al aplicar un posterior proceso de EE para una configuración

para llevar a cabo este análisis en redes de agua . Bargiela ( 1984 )

y

Carpentier and Cohen ( 1991 ) presentan métodos para el AO basados en la teoría

presentan una metodología que también evalúa la incertidumbre asociada a las medidas

y

los parámetros del sistema . Díaz et al . ( 2015 ) adaptan del sector

propia el sector eléctrico , un referente en la implementación de técnicas de EE

y

AO , se pueden encontrar una gran variedad de propuestas : métodos topológicos

red , el procedimiento para cuantificar la incertidumbre en el escenario inicial

y

dada una configuración de medidas determinada , y el consiguiente AO mediante el

incertidumbre en el escenario inicial y dada una configuración de medidas determinada ,

y

el consiguiente AO mediante el denominado método estocástico . Posteriormente ,

se presenta un Ejemplo ilustrativo que permite comprender el potencial del método

y

la información que aporta con respecto a los planteamientos tradicionales . A continuación

Conclusiones asociadas . El método presentado en este trabajo muestra un buen rendimiento

y

permite obtener información adicional de la red en base a los datos recopilados

cierto instante de tiempo , para luego proceder a analizar la incertidumbre inicial

y

condicionada del sistema , en base a las cuales se evaluará su observabilidad .

el valor esperado ( representado con un subíndice b que acompaña a la variable )

y

la desviación o incertidumbre de la variable ( representada con un apóstrofe '

n=nn+nd , donde nn se refiere al número de nodos que no son depósitos

y

nd a los depósitos existentes en el sistema . b ) Flujos en los

del sistema ( demandas ) , positivos cuando se produce una entrada ( aportaciones )

y

nulos cuando no existe ninguna salida ni entrada ( nodos de tránsito ) . Por su parte

en redes de distribución de agua son la ecuación de conservación de la energía

y

la ecuación de continuidad o conservación de masa . Estas expresiones pueden adaptarse

: a ) Ecuación de energía . Establece la relación entre los flujos q

y

los niveles piezométricos en los nodos h , teniendo en cuenta la pérdida de energía

pérdida de energía que se produce a lo largo del tubo debido a su geometría

y

rugosidad . Recogiendo en una matriz de conectividad CP la información relativa

matriz de conectividad CP la información relativa a los nodos de inicio ( +1 )

y

final ( – 1 ) de cada tubo del sistema , esta ecuación se puede

de forma general como : donde K es la matriz diagonal de resistencia hidráulica

y

c el exponente adoptado en función de la ecuación de pérdidas asumida ( c

la Ecuación ( 4 ) en términos de valor esperado e incertidumbre de las variables

y

linealizando la solución en torno al valor medio , se obtiene : Así

: b ) Ecuación de continuidad . Establece la relación entre la demanda b

y

el resto de caudales q que entran o salen de un nodo en cuestión .

en cuestión . Si la información relativa a los flujos que entran ( +1 )

y

salen ( – 1 ) de cada nodo del sistema se agrupa en una matriz

por un signo negativo para que los flujos sean positivos cuando entran al nodo

y

negativos cuando salen de él . Reescribiendo la Ecuación ( 7 ) en términos

de él . Reescribiendo la Ecuación ( 7 ) en términos del valor esperado

y

las desviaciones , y centrando el análisis únicamente en la incertidumbre de la

Reescribiendo la Ecuación ( 7 ) en términos del valor esperado y las desviaciones ,

y

centrando el análisis únicamente en la incertidumbre de la variable , la ecuación

variable , la ecuación de continuidad resulta : Reorganizando las expresiones ( 6 )

y

( 8 ) en una única matriz , las ecuaciones involucradas en la resolución de

en el modelo de incertidumbre de la red existen 2nn+nd+np variables o incógnitas

y

np+nn ecuaciones , por lo que es necesario definir nn+nd condiciones adicionales

apartado anterior . En concreto , se considera que existen pseudomedidas de demanda

y

de niveles en depósitos , ya que en la práctica es fácil que exista una

que en la práctica es fácil que exista una estimación aproximada de los consumos

y

del rango de niveles que se puede dar en los depósitos de la red .

refiere a las desviaciones del nivel piezométrico en los nodos que no son depósitos

y

hd ' a la incertidumbre de los niveles piezométricos en los propios depósitos .

sistema de ecuaciones , designada matriz M . Además , puesto que b '

y

hd ' son variables aleatorias asumidas , el sistema anterior se puede reagrupar en

mismo : que recoge en la parte izquierda de la ecuación los términos desconocidos

y

en la derecha los términos aleatorios conocidos . Nótese que las matrices MU y

y en la derecha los términos aleatorios conocidos . Nótese que las matrices MU

y

MA de esta expresión se refieren a las matrices que aparecen en el primer y

y MA de esta expresión se refieren a las matrices que aparecen en el primer

y

segundo miembro de la Ecuación ( 11 ) , respectivamente . De esta forma ,

es la matriz de varianzas - covarianzas de las distribuciones aleatorias asumidas

y

el resto pueden ser derivadas atendiendo a las siguientes expresiones : Nótese

piezométricos en el sistema ( ) conforme a las distribuciones aleatorias asumidas para b '

y

hd ' . Este valor inicial no proporciona información relevante por sí solo , pero

del que se conoce con certeza que la demanda es nula ( bb=0 ) siempre

y

cuando no existan fugas , no debe incluirse en las ecuaciones de la ( 15

observabilidad determinística . Por tanto , si se considera que existe un vector de medidas X

y

el resto de variables no medidas se denotan como Y , se pueden utilizar las

un vector de medidas X y el resto de variables no medidas se denotan como

Y

, se pueden utilizar las propiedades de una distribución multivariada normal condicionada

propiedades de una distribución multivariada normal condicionada : donde µY|X

y

RY|X representan , respectivamente , la media y la matriz varianza - covarianza de

normal condicionada : donde µY|X y RY|X representan , respectivamente , la media

y

la matriz varianza - covarianza de la distribución condicionada . De forma similar

condicionado implica asumir que la incertidumbre de la lectura es cero . Por tanto , tal

y

como sucede con el cálculo de la incertidumbre inicial , se deben utilizar versiones

se deben utilizar versiones reducidas de las matrices de varianza - covarianza RX

y

RY para calcular RY|X , siendo necesario prestar especial atención al cambio de

RY|X . Análisis de observabilidad Una vez obtenidas las incertidumbres inicial

y

condicionada de las variables del sistema , es posible definir índices que permitan

red para un conjunto de medidas disponible , las diagonales de las matrices RY|X

y

R [ A , U ] contienen la varianza asociada a cada una de las

otros métodos , como la propuesta algebraica de estos mismos autores . Además ,

y

como se verá en el ejemplo ilustrativo , el método estocástico complementa perfectamente

de la red no debería resultar complicado , puesto que se conoce la cota mínima

y

máxima que puede alcanzar el agua en ellos y se puede por tanto deducir la matriz

que se conoce la cota mínima y máxima que puede alcanzar el agua en ellos

y

se puede por tanto deducir la matriz diagonal asociada a la desviación típica de

matriz diagonal asociada a la desviación típica de cada nivel , RAhd . No obstante ,

y

con el objeto de establecer un marco común para la implementación de este método

consiguiente extracción de la matriz para la propagación de incertidumbre , tal

y

como se ha comentado en los apartados anteriores . EJEMPLO ILUSTRATIVO A continuación

como ejemplo está formada por 17 nodos ( 2 depósitos , 4 nodos de tránsito

y

11 nodos de demanda ) y 19 tubos , siendo su disposición en planta la

nodos ( 2 depósitos , 4 nodos de tránsito y 11 nodos de demanda )

y

19 tubos , siendo su disposición en planta la que aparece en la Figura 1 .

necesario partir de un estado hidráulico conocido para poder calcular las matrices MU

y

MA , definidas por las Ecuaciones ( 11 ) y ( 12 ) . Para

poder calcular las matrices MU y MA , definidas por las Ecuaciones ( 11 )

y

( 12 ) . Para ello , se asumen las características de los nodos y

y ( 12 ) . Para ello , se asumen las características de los nodos

y

los tubos recogidas en la Tabla 1 , donde también aparecen las condiciones de flujo

se asume que se dispone de medidas de nivel piezométrico en los depósitos ( ˜h1

y

˜h17 , nótese que la tilde refleja que se trata de medidas disponibles en el sistema

nótese que la tilde refleja que se trata de medidas disponibles en el sistema )

y

de demanda en el resto de los nodos ( ˜q2 , ˜q3 , … ,

en caso de pérdida de dos medidas de demanda , concretamente en los nodos 13

y

14 , ubicados en la zona mallada de la red . Escenario de partida

el método estocástico requiere cuantificar primeramente la incertidumbre inicial

y

condicionada asociada a este escenario . Comenzando por la incertidumbre inicial

ejemplo se asume que la desviación de los niveles piezométricos en los depósitos

y

de la demanda en los nodos de tránsito es nula ( tal y como se

condicionado será necesario considerar la existencia de sensores de nivel en los depósitos

y

medidores de caudal en los nodos de demanda . Asumiendo que la incertidumbre de la

las mismas distribuciones de incertidumbre aleatorias que en el escenario anterior

y

considerando la pérdida de la señal de los caudalímetros empleados para la medida

señal de los caudalímetros empleados para la medida de demanda en los nodos 13

y

14 , en este caso se obtiene SOIh=95 . 58 % , SOIq=96 . 66

, en este caso se obtiene SOIh=95 . 58 % , SOIq=96 . 66 %

y

SOWI=97 . 81 % , lo cual indica la pérdida de observabilidad de la red

obtenidos a nivel individual para cada una de las variables de la red ( flujo

y

nivel piezométrico ) se recogen en la Figura 2 . Esta imagen muestra que la observabilidad

correspondientes a los puntos en los que se mide el nivel piezométrico ( nodos 1

y

17 ) . Esto se justifica porque al dejar de disponer de las dos medidas

dos medidas de demanda en cuestión , el sistema no resulta compatible determinado

y

tan sólo se pueden conocer las variables de estado ( niveles piezométricos ) en

gracias al cálculo del índice de observabilidad , que complementa el resultado

y

da una idea del nivel de observabilidad de la red . Por otro lado , y

llegar a este escenario se ha producido la pérdida simultánea de las medidas ˜q13

y

˜q14 , resulta interesante analizar cómo contribuye la pérdida individual de cada

simplificado supone que la desviación en la medida de los niveles en los depósitos

y

en la medida de demanda en los nodos de tránsito es nula , no siendo

trabajo para el AO proporciona información sobre la observabilidad general de la red

y

sobre cada una de las variables de forma individual , permitiendo además la identificación

por tanto en su capacidad para analizar la observabilidad de la red de una forma rápida

y

además detectar las zonas donde el índice de observabilidad es menor , ya que en

en ellas será necesario añadir aparatos de medida para mejorar la observabilidad

y

por tanto los resultados de la EE . En este sentido , se debe resaltar que

rápida el interés que presenta la incorporación de un grupo de medidas en una zona

y

no otra del sistema , o incluso la recolocación de las medidas redundantes en regiones

rápido e intuitivo para la toma de decisiones en lo que respecta a la monitorización

y

explotación de redes de abastecimiento . CONCLUSIONES Este trabajo plantea un

escenarios : uno primero que asume que la incertidumbre de los niveles en los depósitos

y

las demandas en los nodos sigue distribuciones aleatorias normales , y un segundo

depósitos y las demandas en los nodos sigue distribuciones aleatorias normales ,

y

un segundo en el que esta estimación se ve modificada por la existencia de equipos

forma que este análisis matricial previo permita identificar las medidas redundantes

y

suprimirlas para el posterior análisis estocástico . Así , el método aquí presentado

estado hidráulico de la red en tiempo real para mejorar el funcionamiento del sistema

y

el servicio al usuario . La riada de Valencia de 1957 : reconstrucción hidrológica

servicio al usuario . La riada de Valencia de 1957 : reconstrucción hidrológica

y

sedimentológica y análisis comparativo con la situación actual INTRODUCCIÓN

. La riada de Valencia de 1957 : reconstrucción hidrológica y sedimentológica

y

análisis comparativo con la situación actual INTRODUCCIÓN La ciudad de Valencia

romana , testigo de los numerosos desbordamientos del río Turia . Desde entonces

y

hasta la actualidad , la ciudad se ha ido ampliando y ha tomado medidas para

. Desde entonces y hasta la actualidad , la ciudad se ha ido ampliando

y

ha tomado medidas para defenderse de las fuertes inundaciones . Entre las actuaciones

existente “ Fàbrica de Murs e Valls ” , encargada del mantenimiento de los desagües

y

de la muralla . Esta mejorada institución se creó con el fin de construir los

en el lecho del río : realineamiento de la corriente , refuerzo de las orillas

y

dragado del canal . Con estas prácticas , se pretendía rellenar y nivelar la parte

de las orillas y dragado del canal . Con estas prácticas , se pretendía rellenar

y

nivelar la parte interna de la obra y conseguir un mejor drenaje del agua de

Con estas prácticas , se pretendía rellenar y nivelar la parte interna de la obra

y

conseguir un mejor drenaje del agua de las crecidas ( Carmona , 1997 ) .

( Carmona , 1997 ) . Finalizadas las obras , el problema no queda resuelto

y

la ciudad de Valencia sigue sufriendo inundaciones constantes . La segunda gran

produce como consecuencia de la última riada del Turia , objeto de este estudio

y

comúnmente conocida como la “ Riuà ” de 1957 . La riada de octubre de

extraordinaria se produjo como consecuencia de un proceso convectivo de mesoescala ( Llasat

y

Puigcerver , 1994 ) que afectó fuertemente a la cuenca media y baja del río

mesoescala ( Llasat y Puigcerver , 1994 ) que afectó fuertemente a la cuenca media

y

baja del río Turia , con precipitaciones superiores a los 100 mm en 24 horas

la mañana del día 14 con un caudal máximo de 2700 m3 / s ,

y

una segunda onda de crecida que llegaba a las dos y media con un caudal

m3 / s , y una segunda onda de crecida que llegaba a las dos

y

media con un caudal máximo de 3700 m3 / s . Estos dos desbordamientos dejaron

dejaron la ciudad sumida en la catástrofe , con 81 víctimas , miles de damnificados

y

cuantiosos daños materiales ( Pérez , 2007 ) . Dada la gran conmoción ,

la Solución Sur , que significaba el desvío del río al sur de la ciudad

y

la construcción de una presa en Villamarchante . En 1972 , se termina la segunda

como sólidos , es decir , ha pretendido comprobar si con la información disponible

y

la inclusión del ciclo de sedimentos a un modelo hidrológico distribuido se es

actualidad un evento similar , teniendo en cuenta el cambio en los usos del suelo

y

la finalización de la construcción del embalse de Loriguilla y del Nuevo Cauce

los usos del suelo y la finalización de la construcción del embalse de Loriguilla

y

del Nuevo Cauce del Turia . La adición de este segundo objetivo ha buscado la

Su clima es mediterráneo semiárido con una precipitación media anual de 402 mm

y

una evapotranspiración potencial media anual de 1100 mm . En ella se encuentran

encuentran construidos cuatro embalses : Arquillo San Blas , Benagéber , Loriguilla

y

Buseo . De éstos , cabe destacar que el embalse de Loriguilla se finalizó el

como el fin último del modelo era la reconstrucción de la riada de 1957 ,

y

aguas arriba de Benagéber no se produjeron precipitaciones de carácter significativo

por presentar suelos bajos en contenido de arcillas ( valor medio del 23 % )

y

contenidos más altos de limos en la parte baja y de arena en la parte

medio del 23 % ) y contenidos más altos de limos en la parte baja

y

de arena en la parte media . Litológicamente , predominan las rocas detríticas

en la cuenca media , mientras que en la cuenca baja predominan las rocas calizas

y

cuaternarias ( Figura 2 ) . Respecto a los usos del suelo ( Soriano y

y cuaternarias ( Figura 2 ) . Respecto a los usos del suelo ( Soriano

y

Sanchis , 2014 ) , el cultivo de secano y de regadío , a excepción

usos del suelo ( Soriano y Sanchis , 2014 ) , el cultivo de secano

y

de regadío , a excepción del cultivo de cítricos que ha tomado más importancia

Información Hidrológica de la Confederación Hidrográfica del Júcar ( SAIH - CHJ )

y

el Centro de Estudios Hidrográficos del Centro de Estudios Experimentales ( CEH

recopiló información sobre la precipitación , así como información hidrológica

y

sedimentológica . Precipitación El periodo de precipitaciones se inició el día

sedimentológica . Precipitación El periodo de precipitaciones se inició el día 11

y

finalizó el día 14 . Aunque las precipitaciones de los días 11 y 12 no

día 11 y finalizó el día 14 . Aunque las precipitaciones de los días 11

y

12 no fueron de gran intensidad , éstas contribuyeron a aumentar la escorrentía

la mañana del día 14 en la cuenca media , en cambio , entre las 08

y

las 18 horas , se reprodujo el régimen tormentoso afectando a la cuenca baja .

precipitaciones de gran importancia , pues la variación de volumen no fue alta

y

el caudal medio de salida fue de 6 . 06 m3 / s .

Loriguilla ) se vio afectada por precipitaciones de gran intensidad durante la tarde

y

noche del día 13 , aunque se produjeron precipitaciones débiles desde la primera

hasta su cese a mediodía del 14 . En el caso de las ramblas Castellana

y

Escarihuela , a primera hora del día 13 ya se produjeron precipitaciones de poca

desde la primera hora del día 13 , alcanzando intensidades altas durante la tarde

y

noche del mismo día . Por último , en la margen derecha del Turia , las

día 13 . El régimen tormentoso en la parte alta se produjo durante la tarde

y

noche del día 13 , mientras que en la zona cercana a la costa ,

que en la zona cercana a la costa , éste se produjo durante la mañana

y

las primeras horas de la tarde del día 14 ( García y Carrasco , 1958 ;

produjo durante la mañana y las primeras horas de la tarde del día 14 ( García

y

Carrasco , 1958 ; Cánovas , 1958 ) . Además de esta información cualitativa

de 07 : 00 h a 07 : 00 h ) de los días 13

y

14 de octubre de 1957 , que se interpoló empleando el método de la inversa

inversa de la distancia al cuadrado ( utilizando las 3 estaciones más cercanas )

y

se obtuvo su distribución espacial , que puede observarse en la Figura 4 . Esta

a los encargados de las centrales hidroeléctricas que al encontrarse en servicio

y

ser sorprendidos por la avenida , se convirtieron en testigos , pudo obtenerse los

al embalse de Buseo , se supo que se encontraba a mitad de su capacidad

y

con los desagües de fondo cerrados , ya que no podían operarse debido a la

( 1 . 41 km2 de cauce , 11 . 20 km2 de zona urbana

y

9 . 51 km2 de zona no edificada ) , la altura media que alcanzaron

IMPLEMENTACIÓN DEL MODELO HIDROLÓGICO El modelo empleado fue el modelo hidrológico

y

del ciclo de sedimentos TETIS ( GIMHA , 2014 ) . El modelo hidrológico se

( GIMHA , 2014 ) es un modelo conceptual , con parámetros de base física

y

distribuido , permitiendo así considerar la heterogeneidad espacial de inputs ,

, permitiendo así considerar la heterogeneidad espacial de inputs , parámetros

y

variables de estado . El sub - modelo hidrológico está basado en una estructura

está basado en una estructura de tanques en cada celda , interconectados vertical

y

horizontalmente , representando los procesos de ladera y acuífero . El proceso de

interconectados vertical y horizontalmente , representando los procesos de ladera

y

acuífero . El proceso de acuífero está representado por un único tanque , mientras

nieve , interceptación , almacenamiento en superficie , almacenamiento estático

y

almacenamiento gravitacional . El resultado es la representación de la mayoría

evapotranspiración , infiltración , escorrentía directa , percolación , interflujo , flujo base

y

flujo subterráneo profundo . La escorrentía directa , el interflujo y el flujo

flujo base y flujo subterráneo profundo . La escorrentía directa , el interflujo

y

el flujo base se conectan a la red fluvial mediante la definición de dos áreas

mediante la definición de dos áreas umbrales . La red fluvial se divide en cárcavas

y

cauces . Las laderas ( escorrentía directa e interflujo ) , se conectan a la

( escorrentía directa e interflujo ) , se conectan a la red de cárcavas ,

y

sólo cuando existe flujo base , se conectan a la red de cauces . La

mapa correspondiente por un factor corrector , que es común para todas las celdas

y

diferente para cada parámetro . Por lo tanto , el número de variables a calibrar

el número de variables a calibrar se reduce a 9 ( ocho procesos de ladera

y

uno de propagación en la red fluvial ) . El sub - modelo sedimentológico

TETIS , ambos procesos son unidimensionales . Las tasas de producción , transporte

y

sedimentación están controladas por la disponibilidad de sedimentos en la cuenca

y

la capacidad de transporte de la corriente . Mientras que el transporte de materiales

separa las partículas de sedimentos en tres categorías de tamaño , arena , limo

y

arcilla , cada una de ellas con unas características físicas . De esta forma ,

se calculó el mapa de direcciones de drenaje , el de celdas drenantes acumuladas

y

el mapa de pendientes , necesarios para la correcta ejecución del modelo . Calculados

metodología propuesta por la Organización de las Naciones Unidas para la Alimentación

y

la Agricultura ( FAO , 2006 ) , y a partir de los valores adoptados

las Naciones Unidas para la Alimentación y la Agricultura ( FAO , 2006 ) ,

y

a partir de los valores adoptados en otros estudios desarrollados por el Grupo de

estudios desarrollados por el Grupo de Investigación de Modelación Hidrológica

y

Ambiental ( GIMHA ) . El parámetro de interceptación máxima para cada uno de los

de estudios previos realizados por el GIMHA ( Ramírez et al . , 2012 )

y

adaptados a los usos del suelo presentes en el área de estudio . La capacidad

la capa superior del suelo ( Figura 7b ) depende de la forma del terreno

y

del tipo de suelo . Para su cálculo se emplearon los mapas derivados de la

clasificó el suelo urbano en dos categorías : zonas totalmente urbanas ( ciudad )

y

zonas medias ( casas aisladas con jardín ) . En las primeras se tomó un

con jardín ) . En las primeras se tomó un 20 % del valor calculado

y

en las zonas medias , un 80 % . La capacidad de percolación del substrato

empleando el mapa de permeabilidad cualitativa del acuífero del Instituto Geológico

y

Minero de España ( IGME , 2013 ) . Esta permeabilidad cualitativa se transformó

en valores numéricos recurriendo a valores típicos de bibliografía ( Domenico

y

Schwartz , 1998 ) . En el caso de la conductividad hidráulica horizontal saturada

De igual forma ocurre con la conductividad hidráulica horizontal saturada del substrato

y

la capacidad de percolación en la zona profunda del suelo , pero en este caso

parámetros geomorfológicos se divide en la determinación de las áreas umbrales

y

la determinación de los coeficientes y exponentes empleados en la OCG . Los valores

la determinación de las áreas umbrales y la determinación de los coeficientes

y

exponentes empleados en la OCG . Los valores adoptados para las áreas umbrales fueron

adoptados para las áreas umbrales fueron 0 . 9 km2 en el caso del interflujo

y

15 km2 en el caso del flujo base . Para la OCG , se emplearon

modelo de sedimentos , fueron necesarios los mapas de contenido en arcilla , limo

y

arena ( Figura 2 ) , así como el factor de erosión del suelo ,

Figura 7e ) , el factor de cultivo , factor C ( Figura 7f )

y

el factor de prácticas de conservación , factor P . El factor K se calculó

al . ( 1994 ) . El factor C depende de los usos del suelo

y

se calculó adoptando los valores propuestos en Almorox et al . ( 1994 ) .

, 2012 ) . Por último , fue necesaria la información ( nivel , volumen ,

y

caudal de salida ) sobre los tres embalses situados en la subcuenca de estudio :

sobre los tres embalses situados en la subcuenca de estudio : Benagéber , Loriguilla

y

Buseo ( SAIH - CHJ y CEH - CEDEX ) ; así como la serie

en la subcuenca de estudio : Benagéber , Loriguilla y Buseo ( SAIH - CHJ

y

CEH - CEDEX ) ; así como la serie de caudales del aforo situado en

magnitud . Además , se empleó un segundo evento , el de los días 21

y

22 de septiembre de 2007 , debido a que el registro de caudal realizado por

2007 , debido a que el registro de caudal realizado por el SAIH - CHJ

y

por el CEH - CEDEX durante el evento de octubre de 2000 , fue diferente

CEH - CEDEX durante el evento de octubre de 2000 , fue diferente , tal

y

como se puede observar en la Figura 8a . Por tanto , dada la imposibilidad de

la Figura 8a . Por tanto , dada la imposibilidad de saber qué registro era correcto

y

la necesidad de emplear este evento por ser muy superior a los restantes , se

una calibración manual buscando una solución adecuada para los tres registros ,

y

dando mayor importancia a los registros del evento de 2000 por su magnitud .

los registros del evento de 2000 por su magnitud . Los índices de Nash

y

Sutcliffe ( Nash y Sutcliffe , 1970 ) para el evento de octubre de 2000

de 2000 por su magnitud . Los índices de Nash y Sutcliffe ( Nash

y

Sutcliffe , 1970 ) para el evento de octubre de 2000 fueron de 0 .

octubre de 2000 fueron de 0 . 6 para el registro del SAIH - CHJ

y

0 . 68 para el del CEH - DECEX . En cambio , el evento de

de septiembre de 2007 , no presentó un buen ajuste , su índice de Nash

y

Sutcliffe fue inferior a 0 . 5 , pero dados los buenos resultados obtenidos para

emplearon nueve eventos . Los resultados fueron satisfactorios , el índice de Nash

y

Sutcliffe tomó valores comprendidos entre 0 . 60 y 0 . 81 , salvo en

satisfactorios , el índice de Nash y Sutcliffe tomó valores comprendidos entre 0 . 60

y

0 . 81 , salvo en el caso del evento de Septiembre de 2010 en

a cabo con el periodo temporal 2000 - 2002 , obteniendo un índice de Nash

y

Sutcliffe de 0 . 81 . La validación se llevó a cabo con el

. La validación se llevó a cabo con el periodo temporal 2003 - 2008

y

se obtuvo un NSE de 0 . 64 ( Figura 10 ) .

un NSE de 0 . 64 ( Figura 10 ) . RECONSTRUCCIÓN HIDROLÓGICA

Y

SEDIMENTOLÓGICA Para poder realizar la reconstrucción de la riada , el primer

correctores obtenidos . Se reconstruyó a escala horaria la precipitación de los días 13

y

14 de octubre de 1957 registrada por el Servicio Meteorológico Nacional ( Tabla

además de esta información se emplearon como datos de control , los caudales máximos

y

horas de paso observados . Tras esto , dada la imposibilidad de calibrar el sub

embalse de Benagéber era conocido , pero , en el caso del embalse de Buseo

y

de la ataguía de Loriguilla no se conocía ni el hidrograma de salida ni la

, dato que había sido observado . La laminación del caudal total ( líquido

y

sólido ) se llevó a cabo mediante el método de Puls Modificado ( Ponce 1989

ello se empleó un modelo teóricamente basado , en condiciones de flujo turbulento

y

descarga variable . El modelo elegido fue el propuesto por la Environmental Protection

, vs es la velocidad de sedimentación de las partículas ( m / s )

y

vc es la velocidad crítica de sedimentación de las partículas ( m / s )

. 2 es un factor para tener en cuenta condiciones no ideales de sedimentación

y

A es el área del embalse ( m2 ) . De esta forma los

de atrapamiento para cada una de las tres categorías que emplea el modelo TETIS

y

para cada uno de los embalses son los que se muestran en la Tabla 3 .

se obtiene el volumen de sedimentos que quedará retenido en dicho paso de tiempo

y

que por lo tanto quedará depositado . El volumen restante queda suspendido y es

tiempo y que por lo tanto quedará depositado . El volumen restante queda suspendido

y

es susceptible de ser vertido . Conocido el caudal total de salida , éste deberá

por lo tanto , a partir de esta concentración se puede calcular el caudal líquido

y

sólido de salida . El volumen de sedimentos suspendidos que no ha sido vertido en

en el primer paso de tiempo , queda suspendido para el siguiente paso de tiempo

y

se sumará al volumen de sedimentos de entrada que no quede retenido . Reconstrucción

el hietograma horario de la tormenta en la ciudad ( Figura 11 ) . Tal

y

como se observa , la mayor parte de la precipitación se concentró en una hora

partir de este hietograma , la información cualitativa adicional sobre el inicio

y

el fin de las precipitaciones en la cuenca , así como la comparación parcial entre

cuenca , así como la comparación parcial entre los caudales máximos líquidos

y

sus horas de paso simulados y los caudales máximos y horas de paso observados en

comparación parcial entre los caudales máximos líquidos y sus horas de paso simulados

y

los caudales máximos y horas de paso observados en los puntos de control , fue

caudales máximos líquidos y sus horas de paso simulados y los caudales máximos

y

horas de paso observados en los puntos de control , fue posible realizar la discretización

precipitación era la que generaba el régimen tormentoso . Se consideró el 80 %

y

no el 95 % , debido a que en esta zona la precipitación se inicia

% , debido a que en esta zona la precipitación se inicia por la tarde

y

no a las 12 : 00 como es el caso de Valencia , por lo tanto el

cuenca , salvo en el caso de los pluviómetros situados en las ramblas Castellana

y

Escarihuela , en las que se produjo una situación diferente , el régimen tormentoso

en Valencia se consideró que cada bloque presentó una duración de dos horas ,

y

un medio de la precipitación . De igual forma que en el resto de pluviómetros ,

en las estaciones pluviométricas situadas más al oeste ( cuenca del río Tuéjar

y

alrededores ) se consideró que el 95 % de la precipitación diaria provocó el régimen

consideró que el 95 % de la precipitación diaria provocó el régimen tormentoso

y

en las estaciones de Chera y Bugarra se tomó el 99 % de la precipitación

precipitación diaria provocó el régimen tormentoso y en las estaciones de Chera

y

Bugarra se tomó el 99 % de la precipitación diaria . En el caso de

el 99 % de la precipitación diaria . En el caso de las ramblas Castellana

y

Escarihuela se consideró que el 87 % de la precipitación diaria contribuyó en

ciudad de Valencia . Sabiendo que la superficie inundada fue de 20 . 707 km2

y

que la altura media del sedimento fue de 25 cm , el volumen de sedimentos

se disponía de series de caudal sólido , el proceso de calibración fue manual

y

consistió en la obtención de unos factores correctores a escala diaria y horaria

manual y consistió en la obtención de unos factores correctores a escala diaria

y

horaria que obtuviesen el mismo volumen de sedimentos en la estación de aforo situada

de sedimentos en la estación de aforo situada en Manises al final del evento ,

y

que además , sabiendo que la capacidad del cauce era de 2300 m3 / s

y

realizando las hipótesis de concentración homogénea en toda la sección y depósito

/ s y realizando las hipótesis de concentración homogénea en toda la sección

y

depósito de todos los sedimentos desbordados , proporcionasen un volumen de sedimentos

volumen de 18 , 000 , 452 m3 en el caso de la escala diaria

y

un volumen de 18 , 000 , 449 m3 en la escala horaria , que

muestran los resultados en los diferentes puntos de simulación empleados ( Figura 13

y

Figura 14 ) . Igualmente , se obtuvo la variación de volumen y sedimentos depositados

Figura 13 y Figura 14 ) . Igualmente , se obtuvo la variación de volumen

y

sedimentos depositados , así como los caudales de entrada y salida , tanto de caudal

variación de volumen y sedimentos depositados , así como los caudales de entrada

y

salida , tanto de caudal líquido como sólido , en el embalse de Buseo y

y salida , tanto de caudal líquido como sólido , en el embalse de Buseo

y

en la ataguía de Loriguilla ( Figura 15 ) . Los resultados obtenidos demostraron

precipitación cabe destacar la importancia de la precipitación antecedente de los días 11

y

12 de octubre , que aunque por sí sola no tendría mayor importancia , contribuyó

Así como , la estructura espacio - temporal de la lluvia de los días 13

y

14 , incluyendo el desplazamiento noreste del epicentro de la tormenta . En el caso

la estación de aforo de Manises sufriendo únicamente desfase horario ( Bonache

y

Marco , 2014 ) y es la información más precisa de la que se dispuso

Confederación . De la comparación entre ambos hidrogramas se obtuvo unos errores de +8 %

y

– 1 . 5 % en el caso de los caudales máximos , y –

% y – 1 . 5 % en el caso de los caudales máximos ,

y

– 11 . 5 % en el caudal mínimo . Además , se puede

máximo alcanzado en este embalse hubiese sido inferior a su cota de coronación

y

por lo tanto no se hubiese producido vertido por coronación , lo que suponía una

inmediaciones del embalse . La erosión producida fue de 53 hm3 ( Tabla 4

y

Figura 16 ) , en cambio , el volumen de sedimentos que llegó a Manises fue

presenta pendientes tendidas que producen una disminución de la velocidad del flujo

y

por lo tanto de su capacidad de transporte , pero aguas arriba , en torno a las

. De los 18 hm3 de sedimentos , 11 . 25 hm3 eran de arcilla

y

6 . 15 hm3 de limo , lo que coincide con la descripción de los

capacidad de infiltración del suelo , la conductividad hidráulica del interflujo

y

el factor de cultivo . El proceso llevado a cabo , consistió en realizar

de Loriguilla , para comprobar el efecto del cambio en los usos del suelo ,

y

la segunda con la incorporación del embalse de Loriguilla . Como condiciones iniciales

de la segunda simulación , se inició con un volumen de 17 . 01 hm3 ,

y

se consideró que no sería necesario que realizase sueltas de caudal . Asimismo

la simulación de 1957 , ya que proviene principalmente de las ramblas Castellana

y

Escarihuela . En cambio , la primera onda se vería fuertemente reducida ( Figura

de almacenar todas las aportaciones que se producen aguas arriba del mismo , tal

y

como se consideró inicialmente . De hecho , los resultados mostraron que el nivel

importancia la estructura espacio - temporal de la precipitación de los días 13

y

14 , así como el desplazamiento noreste del epicentro de la tormenta , que provocó

el punto de vista de la peligrosidad , sí puede serlo en otras situaciones ,

y

en cualquier caso , lo es desde el punto de vista de los daños sociológicos y económicos

, y en cualquier caso , lo es desde el punto de vista de los daños sociológicos

y

económicos , pues la retirada de los sedimentos de la ciudad duró más de un

posible inundación generada con las precipitaciones de 1957 . Balances de agua

y

calor en la marisma de Doñana INTRODUCCIÓN Los recursos hídricos , incluyendo

Doñana INTRODUCCIÓN Los recursos hídricos , incluyendo las aguas superficiales

y

subterráneas , son un factor clave en el entorno de Doñana , ya que constituyen

constituyen la base estratégica del mantenimiento de los ecosistemas en general ,

y

en particular de los ambientes marismeños , los sistemas lagunares y los arroyos

en general , y en particular de los ambientes marismeños , los sistemas lagunares

y

los arroyos fluviales . El agua es un elemento esencial en la configuración de este

configuración de este territorio , tanto desde el punto de vista ecológico como económico

y

territorial ( CICCP , 1977 ; Junta de Andalucía , 2010 ) . Sin embargo ,

último siglo una serie de afecciones que han reducido considerablemente su extensión

y

han afectado a su funcionamiento natural y calidad hídrica . Tras el accidente minero

reducido considerablemente su extensión y han afectado a su funcionamiento natural

y

calidad hídrica . Tras el accidente minero de Aznalcóllar en 1998 , se intensificó

minero de Aznalcóllar en 1998 , se intensificó la recuperación de la cantidad

y

calidad de las aportaciones hídricas a la marisma , así como de su funcionalidad

modelo numérico bidimensional Iber , que permite conocer los campos de velocidades

y

profundidades de agua en la marisma ante diferentes escenarios de simulación ( Dolz

. , 2014 ) . También ha participado en el seguimiento de masas de agua

y

comunidades vegetales mediante teledetección ( Martí - Cardona et al . , 2010 ;

Martí - Cardona et al . , 2010 ; Marti - Cardona , 2014 )

y

en el estudio hidrometeorológico y de flujos de agua y calor en superficie y fondo

. , 2010 ; Marti - Cardona , 2014 ) y en el estudio hidrometeorológico

y

de flujos de agua y calor en superficie y fondo ( Prats et al .

- Cardona , 2014 ) y en el estudio hidrometeorológico y de flujos de agua

y

calor en superficie y fondo ( Prats et al . , 2011 , 2012 ;

) y en el estudio hidrometeorológico y de flujos de agua y calor en superficie

y

fondo ( Prats et al . , 2011 , 2012 ; Ramos Fuertes , 2012

, 2012 ; Ramos Fuertes , 2012 ) . La limnología tiene un importante

y

creciente papel en la toma de decisiones en relación a la gestión del agua embalsada

conocimiento limnológico de la marisma , estudiar el balance térmico de la misma

y

comprender los procesos de intercambio de calor que ocurren en ella resulta de interés

en los medios acuáticos , que afecta a un gran número de parámetros de calidad

y

a las reacciones químicas y biológicas que tienen lugar en ellos . En consecuencia

que afecta a un gran número de parámetros de calidad y a las reacciones químicas

y

biológicas que tienen lugar en ellos . En consecuencia , la evolución , distribución

y

ecología de los organismos acuáticos se encuentra profundamente afectada por la

Íntimamente relacionada con ella se encuentra el régimen de radiación solar en superficie

y

en el interior de la columna de agua , principal condicionante de la estructura térmica

columna de agua , principal condicionante de la estructura térmica de la masa de agua

y

determinante del desarrollo de los organismos acuáticos fotosintéticos . Por último

ámbito de Doñana es la evapotranspiración , flujo común a los balances de agua

y

calor que , como principal responsable del vaciado de la marisma , afecta a la

principal responsable del vaciado de la marisma , afecta a la disponibilidad de agua

y

al hidroperiodo o duración de la inundación . El trabajo aquí presentado tiene

la caracterización del balance térmico de la marisma de Doñana a escala horaria

y

estacional , haciendo especial énfasis en la evaporación y en el estudio del proceso

Doñana a escala horaria y estacional , haciendo especial énfasis en la evaporación

y

en el estudio del proceso de vaciado de la misma . Ello pretende contribuir al

conocimiento de los procesos físicos que determinan la disponibilidad de agua en la marisma

y

a la toma de decisiones en el manejo hídrico de la misma . MATERIAL

Y

MÉTODOS Área de estudio La marisma de Doñana es un humedal de reconocida

ecológica , principalmente como lugar de paso , cría e invernada de aves europeas

y

africanas . Situada en las provincias de Huelva y Sevilla , la actual marisma natural

cría e invernada de aves europeas y africanas . Situada en las provincias de Huelva

y

Sevilla , la actual marisma natural se extiende sobre aproximadamente 300 km2 en

Su conformación es el resultado de la colmatación del antiguo estuario del río

y

carece prácticamente de relieve , si bien las pequeñas variaciones de nivel del

1 , definen elementos microtopográficos como caños ( cursos internos de agua )

y

lucios ( depresiones ) que determinan la circulación y almacenamiento del agua en

cursos internos de agua ) y lucios ( depresiones ) que determinan la circulación

y

almacenamiento del agua en la marisma . Los procesos anuales de inundación y secado

circulación y almacenamiento del agua en la marisma . Los procesos anuales de inundación

y

secado se prolongan típicamente de mediados de otoño a mediados de verano ( Aragonés

torno a 160 hm3 anuales teniendo en cuenta el promedio histórico de precipitación )

y

la escorrentía superficial de arroyos tributarios por su límite noroeste ( Dolz

y

Velasco , 1990 ) . Los arroyos de mayor importancia son La Rocina ( 416

. Los arroyos de mayor importancia son La Rocina ( 416 km2 de cuenca )

y

El Partido ( 291 km2 de cuenca ) , para los cuales no se dispone

De manera orientativa , se asigna a las aguas dulces del arroyo de la Rocina

y

colindantes , con una cierta estabilidad en sus caudales a lo largo del año por

media anual entre 1 - 1 . 5 hm3 / año en ciclos de sequía

y

un máximo de 73 hm3 / año ( MMA , 2001 ) . La cubeta

/ año ( MMA , 2001 ) . La cubeta de la marisma es impermeable

y

la principal aportación freática se produce a través del caudal base de los arroyos

se producen en las zonas de contacto con sustratos arenosos en los límites norte

y

oeste de la marisma , y en los denominados “ ojos de la marisma ”

de contacto con sustratos arenosos en los límites norte y oeste de la marisma ,

y

en los denominados “ ojos de la marisma ” , pequeñas surgencias que atraviesan

actual se encuentra considerablemente mermada por la sobreexplotación del acuífero ,

y

la descarga directa neta en los límites de la marisma se estima del orden de

En el marco del presente estudio , la comunicación de la marisma con el estuario

y

el acuífero subyacente se considera nula a efectos de balances de agua y calor .

estuario y el acuífero subyacente se considera nula a efectos de balances de agua

y

calor . Los detalles de las hipótesis adoptadas sobre el comportamiento hídrico

comportamiento hídrico de la marisma se describen en la sección " Balance hídrico

y

curvas de agotamiento " . Instrumentación de medida El estudio del balance térmico

una red de estaciones de medida localizadas en áreas no vegetadas de la marisma

y

que ha estado operativa entre los años 2006 y 2011 ( Figura 1 ) .

áreas no vegetadas de la marisma y que ha estado operativa entre los años 2006

y

2011 ( Figura 1 ) . La red de estaciones cuenta con un total de

tres de temperatura del suelo a diferente profundidad ( D03 , D04 , D05 )

y

dos estaciones meteorológicas completas ( D04 , D05 ) . La información presentada

presentada en este trabajo consta de medidas de nivel de agua en todas las estaciones

y

de medidas térmicas y meteorológicas de la estación D04 . La Tabla 1 lista el

consta de medidas de nivel de agua en todas las estaciones y de medidas térmicas

y

meteorológicas de la estación D04 . La Tabla 1 lista el equipamiento completo de

) . Período de estudio Dado el marcado carácter pluvial de la marisma

y

la alta variabilidad propia del clima mediterráneo , los aportes de agua superficial

izquierda , que muestra medidas de profundidad de agua en diferentes puntos de la marisma

y

años hidrológicos . De los cinco años monitorizados , el año hidrológico 2006

, mantenimiento de niveles de agua ( diciembre de 2006 - febrero de 2007 )

y

vaciado ( marzo - julio de 2007 ) . El estudio del descenso del nivel

partir del mes de marzo de los años hidrológicos comprendidos entre 2006 / 2007

y

2010 / 2011 . El balance térmico superficial sobre la masa de agua se

donde S es el cambio en el contenido de calor por unidad de tiempo

y

unidad de área aplicado a la columna de agua en su total profundidad , RN

condensación ) , H el flujo de calor sensible intercambiado por contacto con la atmósfera

y

G el flujo de calor intercambiado con el sedimento subyacente , todos en W ·

de calor intercambiado con el sedimento subyacente , todos en W · m – 2

y

de signo positivo cuando suponen un aporte de calor para el agua . Se ignoran

un aporte de calor para el agua . Se ignoran en este esquema otras entradas

y

salidas de calor consideradas de menor importancia como son los flujos de calor debidos

topografía de la marisma , se presupone un lento movimiento de agua en la misma

y

no se considera la advección de calor en el plano horizontal . Ello no implica

calcula como : donde h representa el calado ( m ) , ∆Tw (

y

) es la variación de temperatura del agua ( °C o K ) ocurrida en

tiempo Δt ( s ) en una capa de grosor dy situada a una distancia

y

del suelo ( m ) , ρw es la densidad del agua ( 1000 kg

, ρw es la densidad del agua ( 1000 kg · m - 3 )

y

cw el calor específico del agua ( 4187 J · °C – 1 · kg

El cálculo del calor transmitido al suelo G con medidas discretas en el espacio i

y

el tiempo j toma la siguiente forma : siendo Uj el contenido de calor del

en sus extremos , ∆zi es el grosor de la capa i ( m )

y

∆t = 600 s . Se consideran tres capas de suelo , delimitadas por los

4 sensores de temperatura instalados a 0 , 0 . 25 , 0 . 75

y

1 m de profundidad . El intercambio de calor con la atmósfera ( H

función del ratio de Bowen β , dependiente de los gradientes observados de temperatura

y

humedad en superficie : siendo γ la constante psicrométrica ( 0 . 07 hPa ·

de vapor de saturación a la temperatura de la superficie de agua ( hPa )

y

ea la presión de vapor del aire ( hPa , ambas estimadas con medidas de

aire ( hPa , ambas estimadas con medidas de temperatura del agua , del aire

y

humedad ambiental ) . Finalmente , H y λE se definen como : Balance

de temperatura del agua , del aire y humedad ambiental ) . Finalmente , H

y

λE se definen como : Balance hídrico y curvas de agotamiento Con el

ambiental ) . Finalmente , H y λE se definen como : Balance hídrico

y

curvas de agotamiento Con el objetivo de evaluar la participación del balance térmico

estimaciones con la variación de calado observada . El balance de agua es muy simplificado

y

contempla únicamente las entradas por precipitación directa y las pérdidas por

muy simplificado y contempla únicamente las entradas por precipitación directa

y

las pérdidas por evaporación . Este balance supone que no hay contacto con el acuífero

justifica por la impermeabilidad del sedimento , compuesto principalmente por arcillas

y

limos ( Ruiz et al . , 2004 ) . Respecto a la segunda hipótesis , se

hipótesis , se puede suponer válida durante una gran parte del periodo de inundación ,

y

en particular durante el periodo de vaciado , a causa de la reducción de caudales

el periodo de vaciado , a causa de la reducción de caudales naturales de entrada

y

salida de agua por escorrentía superficial debida a la sobreexplotación de acuíferos

de agua por escorrentía superficial debida a la sobreexplotación de acuíferos

y

a la reconstrucción y prolongación del dique de la Montaña del Río ( Bayán y

escorrentía superficial debida a la sobreexplotación de acuíferos y a la reconstrucción

y

prolongación del dique de la Montaña del Río ( Bayán y Dolz , 1995 ;

y a la reconstrucción y prolongación del dique de la Montaña del Río ( Bayán

y

Dolz , 1995 ; Rodríguez - Ramírez y Clemente Salas , 2002 ) .

de la Montaña del Río ( Bayán y Dolz , 1995 ; Rodríguez - Ramírez

y

Clemente Salas , 2002 ) . Las observaciones de campo denotan una cierta concordancia

concordancia de los ritmos de vaciado de la marisma con independencia del lugar de medida

y

del año hidrológico , especialmente sobre la cota aproximada de 1 . 3 msnm (

de vaciado , las series de cotas superiores a 1 . 3 msnm se filtraron

y

agruparon en el eje de ordenadas . En primer lugar , se seleccionaron sólo aquellos

se seleccionaron sólo aquellos tramos de cota descendiente ( excluyendo avenidas

y

eventos de precipitación ) con duración superior a 30 días ( los huecos en las

, que pertenecen a un total de cinco puntos de medida , cinco años hidrológicos

y

contienen 1325 datos diarios de cota de agua . Dado que la magnitud de la

de las condiciones hídricas particulares de cada año hidrológico . RESULTADOS

Y

DISCUSIÓN Balance térmico diario La evolución media horaria de los componentes

es el principal destino del calor incorporado por radiación durante el día . RN

y

S alternan valores positivos diurnos ( ganancia de calor para el agua ) y negativos

RN y S alternan valores positivos diurnos ( ganancia de calor para el agua )

y

negativos nocturnos ( pérdida de calor ) . Durante las fases de llenado y mantenimiento

) y negativos nocturnos ( pérdida de calor ) . Durante las fases de llenado

y

mantenimiento de nivel , el intercambio de calor sensible y latente entre agua y

las fases de llenado y mantenimiento de nivel , el intercambio de calor sensible

y

latente entre agua y atmósfera ( H + λE ) supone una leve pérdida de

y mantenimiento de nivel , el intercambio de calor sensible y latente entre agua

y

atmósfera ( H + λE ) supone una leve pérdida de calor para el agua

siendo predominante a partir del mediodía por la intensificación típica del viento

y

del gradiente de presión de vapor entre las 12 h y las 18 h en

intensificación típica del viento y del gradiente de presión de vapor entre las 12 h

y

las 18 h en esa época . El calor intercambiado con el suelo inundado (

con respecto al análisis horario ( Figura 3 ) : dada la alternancia de ganancias

y

pérdidas de calor por parte de la columna de agua durante el día y la

ganancias y pérdidas de calor por parte de la columna de agua durante el día

y

la noche , respectivamente , el valor neto de S es poco relevante a escala

incorporando la intensificación de la radiación neta entre los meses de febrero

y

julio , consecuencia de la mayor radiación solar incidente . El intercambio de calor

radiación solar incidente . El intercambio de calor por conducción con la atmósfera H

y

con el suelo G son por lo general poco importantes . Entre los meses de noviembre

y

febrero , la pérdida de calor por evaporación supera a la entrada de calor por

la pérdida de calor por evaporación supera a la entrada de calor por radiación

y

por conducción . El resultado es una pérdida neta de calor y una disminución de

calor por radiación y por conducción . El resultado es una pérdida neta de calor

y

una disminución de la temperatura del agua en este periodo . A partir del mes

temperatura del agua en este periodo . A partir del mes de marzo la evaporación

y

la radiación neta tienden a igualarse , resultando en una pérdida de calor neta

acusada reducción del volumen de agua ( Figura 2 ) . Balance de agua

y

proceso de vaciado El balance de agua de la Figura 5 muestra los volúmenes mensuales

muestra los volúmenes mensuales de precipitación acumulada , evaporación acumulada

y

la variación de profundidad del agua observada durante el periodo de inundación

con la cota ) , información obtenida del modelo digital del terreno ( Cambronero

y

De Pourq , 2006 ; Ramos Fuertes , 2012 ) . Los meses más

Ramos Fuertes , 2012 ) . Los meses más lluviosos ( noviembre , enero

y

febrero ) presentan una discordancia importante entre el balance simplificado de

febrero ) presentan una discordancia importante entre el balance simplificado de agua

y

la variación real del calado , debido a que se producen entradas de agua por

de la marisma . El error absoluto del balance estimado es entre un 14 %

y

un 25 % de la variación real del calado en los meses de marzo a

% de la variación real del calado en los meses de marzo a junio ,

y

de un 36 % en diciembre . Las causas más probables del mismo son la

advección de caudal en el plano horizontal ( más probable en otoño - invierno )

y

la influencia de la transpiración vegetal en la fase de vaciado ( primavera - verano

disponible . La cota de la lámina de agua en diferentes puntos de medida

y

años hidrológicos presenta un ritmo de agotamiento similar con independencia de

et al . ( 2004 ) , con el límite superior de los principales lucios

y

caños de la marisma . Bajo esta cota , la marisma se encuentra en gran medida desconectada

la marisma . Bajo esta cota , la marisma se encuentra en gran medida desconectada

y

las diferentes masas de agua se comportan en su vaciado final de forma independiente

en su vaciado final de forma independiente , de acuerdo a la morfología del vaso

y

las condiciones hidrometeorológicas de cada periodo . Al analizar las series de

la intensificación de la radiación solar de marzo a julio ( Figura 4 ) ,

y

de la curva de capacidad de la marisma conectada . La pendiente de esta parábola

marisma de Doñana durante un ciclo típico de inundación a escalas temporales horaria

y

mensual . A escala horaria , destaca el papel del flujo de calor que absorbe

del flujo de calor que absorbe o libera la columna de agua como principal sumidero

y

fuente de la radiación neta en superficie . El intercambio térmico con la atmósfera

El intercambio térmico con la atmósfera es siempre favorable a la evaporación

y

se intensifica notablemente después del mediodía en la fase de vaciado de la marisma

poco relevante en el punto de estudio . A escala mensual , las ganancias

y

pérdidas de calor por parte del agua se compensan y resulta en un balance térmico

escala mensual , las ganancias y pérdidas de calor por parte del agua se compensan

y

resulta en un balance térmico reducido . El creciente aporte de calor a través

comparación de un balance simplificado de masa , basado en la precipitación medida

y

la evaporación estimada , comparado con la variación de calado observada . Utilizando

calado observada . Utilizando 18 series de vaciado de diferentes puntos de medida

y

años hidrológicos , se ha establecido un patrón de agotamiento de la lámina de

Análisis integral del impacto del Cambio Climático en los regímenes de agua , crecidas

y

sedimentos de una rambla mediterránea INTRODUCCIÓN El Cambio Climático es el

cambio en las características meteorológicas durante periodos de tiempo prolongados

y

que son el resultado de las modificaciones en el equilibrio del balance energético

modificaciones en el equilibrio del balance energético de la Tierra causado por procesos

y

agentes naturales y / o antropogénicos . Con cierta incertidumbre , se asume en

equilibrio del balance energético de la Tierra causado por procesos y agentes naturales

y

/ o antropogénicos . Con cierta incertidumbre , se asume en el mundo científico

en el periodo 1950 - 2010 , con efectos como : la fusión de hielo

y

nieve , el calentamiento de los océanos , la elevación del nivel del mar ,

de los océanos , la elevación del nivel del mar , alteraciones en la magnitud

y

frecuencia de extremos climáticos , … que ya han sido observados en la segunda

siglo XX ( Rodríguez et al . , 2014 ) . El Cambio Climático

y

sus efectos en el Ciclo Hidrológico es un tema muy importante en todo el planeta

mediterráneo español , dado que actualmente ya presentan un alto estrés hídrico

y

fuertes alteraciones antrópicas de su régimen . La utilización de modelos hidrológicos

preventivas adecuadas . El objetivo de este trabajo es mostrar una metodología integral

y

actual que permita estimar los posibles efectos del Cambio Climático en los ciclos

permita estimar los posibles efectos del Cambio Climático en los ciclos del agua

y

sedimentos de una cuenca representativa de la hidrología de las cuencas del arco

Viuda , utilizando el modelo hidrológico distribuido TETIS ( Francés et al . , 2007

y

Bussi et al . , 2014 ) y las proyecciones climáticas para la geografía española

( Francés et al . , 2007 y Bussi et al . , 2014 )

y

las proyecciones climáticas para la geografía española puestas a disposición

inglés ) en su quinto informe AR5 ( AEMET , 2017 ) . MATERIAL

Y

MÉTODOS Caso de estudio Como ya se ha indicado , el caso de

Mar Mediterráneo . Dentro de la misma se encuentran los embalses de Mª . Cristina

y

Alcora ( Figura 1 . b ) , el primero de ellos con problemas de aterramiento

y

probable infra - dimensionamiento de su aliviadero , por lo que es prioritario el

clima mediterráneo , caracterizado por ser una transición entre las zonas templadas

y

cálidas , con un régimen estacional de temperaturas parecido a los climas templados

estacional de temperaturas parecido a los climas templados , mientras por su aridez

y

escasas lluvias se asemeja a los climas tropicales secos . Tiene una litología común

donde se ubica , marcada por un dominio calcáreo , en buena parte karstificado ,

y

que en particular favorece la infiltración de agua en una buena parte de su red

TETIS es un modelo conceptual distribuido , de parámetros físicamente basados

y

que ha sido desarrollado para realizar simulación hidrológica y de ciclo de sedimentos

físicamente basados y que ha sido desarrollado para realizar simulación hidrológica

y

de ciclo de sedimentos en cuencas naturales con el objetivo de obtener de la mejor

procesos hidrológicos , incorporando la variabilidad espacial del ciclo hidrológico

y

de sedimentos con el menor número de parámetros posibles , cumpliendo por tanto

. Las entradas del modelo son series temporales de : precipitación , temperatura

y

evapotranspiración de referencia o potencial ( GIMHA , 2017 ) . La representación

( GIMHA , 2017 ) . La representación de los procesos físicos de ladera

y

acuífero en cada celda del sub - modelo hidrológico está basada en una estructura

modelo hidrológico está basada en una estructura de tanques interconectados vertical

y

horizontalmente , como se representa en las Figuras 2 y 3 . El acuífero es

interconectados vertical y horizontalmente , como se representa en las Figuras 2

y

3 . El acuífero es representado a través del tanque ( T4 ) . Los

) , superficie ( T2 ) ; agua capilar del suelo ( T1 ) ;

y

agua gravitacional del suelo ( T3 ) . Por otra parte , mientras la escorrentía directa

suelo ( T3 ) . Por otra parte , mientras la escorrentía directa , el interflujo

y

el flujo base se conectan con la red fluvial representada por otro tanque ( T5

MED , hasta alcanzar la red de drenaje formada por cárcavas ( flujo intermitente )

y

cauces ( flujo permanente ) . Por tanto existen tres sistemas diferenciados , definidos

CASC2D - SED , que simula los procesos de sedimentos en ladera en dos dimensiones

y

con una aproximación unidimensional en canales ( Rojas , 2002 ) . Para el transporte

Para el transporte en ladera , el sub - modelo utiliza la ecuación de Kilinc

y

Richardson modificada para estimar la capacidad de transporte ( Julien , 2010 ) .

está en función del caudal , la pendiente del terreno , el tipo de suelo

y

los usos de la tierra . El cálculo en cárcavas y cauce se lleva a

el tipo de suelo y los usos de la tierra . El cálculo en cárcavas

y

cauce se lleva a cabo usando la ecuación de Engelund y Hansen ( 1967 )

El cálculo en cárcavas y cauce se lleva a cabo usando la ecuación de Engelund

y

Hansen ( 1967 ) utilizando parámetros hidráulicos ( radio hidráulico , velocidad

1967 ) utilizando parámetros hidráulicos ( radio hidráulico , velocidad de flujo

y

ángulo de fricción ) y características de las partículas ( gravedad específica

parámetros hidráulicos ( radio hidráulico , velocidad de flujo y ángulo de fricción )

y

características de las partículas ( gravedad específica y diámetro ) , aplicando

ángulo de fricción ) y características de las partículas ( gravedad específica

y

diámetro ) , aplicando el concepto de potencia de corriente de Bagnold y el principio

específica y diámetro ) , aplicando el concepto de potencia de corriente de Bagnold

y

el principio de similitud , obteniendo la concentración de sedimentos por peso con

llevó a cabo utilizando la información meteorológica observada de SPAIN02 v . 04

y

los caudales de entrada al embalse de Mª Cristina , obtenidos de la base de

de Mª Cristina , obtenidos de la base de datos del CEDEX . La calibración

y

validación del modelo se realizó de manera cualitativa mediante la comparación

modelo se realizó de manera cualitativa mediante la comparación de los hidrogramas

y

cuantitativamente utilizando el coeficiente de eficiencia de Nash y Sutcliffe ( EF2

los hidrogramas y cuantitativamente utilizando el coeficiente de eficiencia de Nash

y

Sutcliffe ( EF2 ) . En este proceso fue necesario la introducción de pérdidas de

se correlacionaron las pérdidas con el área de cauce sobre las zonas kársticas

y

de tramos de cauce donde se identifica que existe conexión indirecta con el acuífero

ello el EF2 en calibración mejora de 0 . 905 a 0 . 928 ,

y

en validación de 0 . 739 a 0 . 924 , con respecto a una simulación previa

cuenca se vio dificultada debido a la falta de registros de la carga en suspensión

y

de lecho que circulan en ella . Sin embargo , se aprovechó información batimétrica

acumulado en el embalse no representa el total del volumen de sedimentos producidos

y

transportados por la cuenca hasta dicho punto , se utilizó el índice de eficiencia

transportada hasta el embalse . Esta información sirvió de referencia para modelizar

y

calibrar el sub - modelo sedimentológico de TETIS , obteniendo como resultado una

aproximación a la variación temporal de los procesos de producción , transporte

y

depósito de sedimentos dentro de la cuenca , principalmente en puntos de interés

informe AR5 . Las series temporales utilizadas son datos diarios de precipitación

y

temperaturas máxima y mínima , correspondientes a dos periodos : un primero , llamado

temporales utilizadas son datos diarios de precipitación y temperaturas máxima

y

mínima , correspondientes a dos periodos : un primero , llamado de Control ( 1961

que debe presentar características similares a los datos observados en el pasado ;

y

un segundo periodo que denominaremos Futuro ( 2006 - 2100 ) que trata de representar

( Figura 6 ) . Esto puede ser debido a la resolución de los MCG

y

/ o a las hipótesis de las técnicas de regionalización . Es por ello que fue

técnicas de regionalización . Es por ello que fue necesario descartar los peores

y

corregir el sesgo de las series de precipitación y temperatura del resto mediante

necesario descartar los peores y corregir el sesgo de las series de precipitación

y

temperatura del resto mediante regresiones lineales para cada estación del año

cada estación del año . Estas regresiones se estimaron durante el período común

y

fueron diferentes para cada modelo regional . Los métodos de corrección del sesgo

Los métodos de corrección del sesgo en la información proveniente de los MCG

y

sus regionalizaciones generalmente están basados en transformaciones estadísticas

al . , 2012 ) . Un ejemplo de la diferencia entre la información corregida

y

sin corregir se puede observar en la Figura 7 . Es de destacar la mala representación

presenta un régimen estacional de temperaturas similar al de los climas templados ,

y

la aridez y escasas precipitaciones características de los climas secos ; y a su

régimen estacional de temperaturas similar al de los climas templados , y la aridez

y

escasas precipitaciones características de los climas secos ; y a su litología

templados , y la aridez y escasas precipitaciones características de los climas secos ;

y

a su litología marcada por un dominio calcáreo que dificulta la conexión entre

litología marcada por un dominio calcáreo que dificulta la conexión entre río

y

acuífero y favorece la infiltración de agua durante su propagación en la red de

marcada por un dominio calcáreo que dificulta la conexión entre río y acuífero

y

favorece la infiltración de agua durante su propagación en la red de cauces .

correlacionaron dichas pérdidas con el área de cauce sobre las zonas kársticas

y

de tramos de cauce donde se identifica que existe conexión indirecta con el acuífero

hidrológico , lo que se vio reflejado en una mejor representación del hidrograma

y

en un mejor valor del coeficiente de eficiencia de Nash y Sutcliffe . La

representación del hidrograma y en un mejor valor del coeficiente de eficiencia de Nash

y

Sutcliffe . La implementación del sub - modelo sedimentológico se realizó utilizando

efecto de consolidación del sedimento utilizando la fórmula de Miller ( 1953 )

y

considerando la relación del volumen de sedimento transportado y acumulado mediante

Miller ( 1953 ) y considerando la relación del volumen de sedimento transportado

y

acumulado mediante el uso del índice de eficiencia de atrapamiento mediante la ecuación

Debido a errores en esta información , relacionados a la resolución de los MCG

y

/ o a las hipótesis de las técnicas de regionalización , fue necesario descartar

las hipótesis de las técnicas de regionalización , fue necesario descartar los peores

y

corregir las series de precipitación y temperatura del resto mediante regresiones

regionalización , fue necesario descartar los peores y corregir las series de precipitación

y

temperatura del resto mediante regresiones lineales , para obtener una mejor representación

evapotranspiración de referencia . Esta situación incrementará la evapotranspiración real

y

dará lugar a una reducción de los recursos disponibles . En lo que respecta a

erosivos con la producción de escorrentía directa , la reducción en la frecuencia

y

magnitud de las crecidas , significará también una disminución de los aportes

balsas para riego INTRODUCCIÓN A partir de los años 60 del siglo pasado

y

debido a la necesidad de aumentar la garantía de suministro en los cultivos , se

inexistencia de nuevas obras de regulación en los ríos , a la modernización de regadíos

y

a la utilización de aguas regeneradas así como desaladas , ha aumentado considerablemente

España es de 1230 , mientras que para balsas se estima un número entre 50000

y

80000 según González - Ortega ( 2010 ) . Estas obras se distribuyen a lo

obras se distribuyen a lo largo del país , encontrándose principalmente en el sur

y

suroeste ( principalmente en la Comunidad Valenciana , Región de Murcia y Andalucía

el sur y suroeste ( principalmente en la Comunidad Valenciana , Región de Murcia

y

Andalucía ) , constituyendo uno de los pilares básicos del sector agrario de regadío

seguridad es una exigencia más en el desarrollo de las sociedades modernas actuales

y

el nivel de seguridad exigido para las balsas de tierra , y en general para cualquier

sociedades modernas actuales y el nivel de seguridad exigido para las balsas de tierra ,

y

en general para cualquier infraestructura , es cada vez mayor . En este sentido ,

Especiales de Protección Civil ( PEPC ) . Estos planes hacen posible la coordinación

y

actuación de los distintos servicios y Administraciones implicadas . Las inundaciones

Estos planes hacen posible la coordinación y actuación de los distintos servicios

y

Administraciones implicadas . Las inundaciones consideradas por la DBPPC son todas

consideradas por la DBPPC son todas aquellas que representen un riesgo para la población

y

los bienes , produciendo daños en infraestructuras básicas o la interrupción de

operación incorrecta de obras de infraestructura hidráulica . Para la planificación

y

gestión de estos fenómenos , la DBPPC requiere delimitar las zonas que pudieran

infraestructuras ( balsas ) los necesarios estudios de zonificación territorial

y

estimación de los daños por rotura como herramienta para la gestión de futuras

inundadas , así como las posibles afecciones a personas , edificaciones , carreteras

y

/ o infraestructuras . Además , para mejorar la planificación y gestión de la

edificaciones , carreteras y / o infraestructuras . Además , para mejorar la planificación

y

gestión de la emergencia también es necesario conocer : la estimación de los tiempos

la onda de rotura , el tiempo para el cuál se obtienen los valores máximos

y

los tiempos de permanencia de la onda de rotura . Este tipo de estudios

tradicionalmente para las presas , existiendo gran experiencia al respecto ( Bruau García

y

García Navarro , 2001 ; Kumar et al . , 2017 ) . Las balsas

diferencias que suponen un tratamiento distinto en el análisis de los fenómenos de rotura

y

su posterior propagación . Sánchez - Romero ( 2014 ) indica que las principales

en cauces definidos . Por lo tanto , la ubicación de la rotura suele ser única

y

el flujo producido puede analizarse tanto de manera unidimensional como bidimensional

se realiza por zonas secas sin cauces definidos , que nunca han sido inundadas ,

y

en el que los modelos unidimensionales no son capaces de definir correctamente los

cálculo de estos caudales presenta cierta complejidad al ser un fenómeno natural

y

depender de multitud de variables . La falta de información en las series históricas

etc . . Las balsas reciben caudales significativamente menores , siendo conocidos

y

controlados . Ante cualquier eventualidad se puede actuar rápidamente sobre estos

diferencia entre considerar las distintas formas de colapso ( e . g . , tubificación

y

sobrevertido ) . A partir de los resultados obtenidos en las simulaciones hidráulicas

de los posibles daños , permitiendo una mejor interpretación de las afecciones

y

ayudando en la posible gestión en caso de emergencia . El presente artículo debe

, los cuáles deben ser acometidos desde la aprobación del RD _ 09 / 2008

y

que actualmente la mayoría de balsas para riego carecen de los mismos . MATERIAL

Y

MÉTODOS Proceso El proceso desarrollado ( Figura 1 ) describe las diferentes fases

Fase 4 ) propone el tratamiento de los resultados obtenidos , así como su análisis

y

síntesis para obtener como resultado final , los mapas de inundación generales

y

descriptivos para toda la zona , además de la información detallada de aquellos

las características más significativas ( volumen movilizable , altura interior

y

exterior , taludes y ancho de coronación ) , es necesaria la obtención de la

más significativas ( volumen movilizable , altura interior y exterior , taludes

y

ancho de coronación ) , es necesaria la obtención de la curva altura - área

de la zona , ya que puede determinar la elección del modelo numérico ( 1D

y

/ o 2D ) de propagación . El comportamiento hidrodinámico de carácter bidimensional

permite establecer las ubicaciones de las posibles roturas en el dique de la balsa

y

por lo tanto determinar el número de escenarios de rotura a analizar . Se deberán

carreteras , líneas eléctricas , etc . ) , siendo necesario la toma de fotografías

y

datos relativos a su ubicación . Además , en este tipo de estudios , es

produce de manera progresiva . En la actualidad existen diferentes metodologías

y

modelos numéricos que permiten predecir la forma y progresión de la brecha , así

existen diferentes metodologías y modelos numéricos que permiten predecir la forma

y

progresión de la brecha , así como determinar el hidrograma de rotura según los

, 2011 ) . Dentro de estos modelos destacan los paramétricos por su versatilidad

y

su fácil implementación en los modelos de propagación de la onda de rotura ( Greimann

paramétricos permiten al usuario la introducción de la forma final de la brecha

y

del tiempo de formación , así como otro tipo de parámetros que definen de una

debido a su capacidad de análisis de rotura de diques , presentando mejores opciones

y

técnicas que los demás modelos frente a problemas de estabilidad y convergencia

mejores opciones y técnicas que los demás modelos frente a problemas de estabilidad

y

convergencia numérica , además de tener una difusión gratuita . Esto ha permitido

Esto ha permitido que gran número de investigadores los hayan utilizado , testeado

y

modificado , siendo utilizados para la calibración de otros códigos numéricos

modificado , siendo utilizados para la calibración de otros códigos numéricos ( Bozkus

y

Kasap , 1998 ; Villanueva et al . , 1999 ; National Weather Service ,

de rotura ( Output 2 . 1 ) . La metodología desarrollada por Fread

y

Lewis ( 1988 ) e implementada en el modelo FLDWAV para el cálculo del hidrograma

caso balsa ) , el segundo proceso es el cálculo del caudal instantáneo de rotura

y

el último es el proceso de erosión y formación de la brecha . Respecto al

cálculo del caudal instantáneo de rotura y el último es el proceso de erosión

y

formación de la brecha . Respecto al balance de volúmenes el modelo tiene en cuenta

caudal entrante al sistema , la geometría del embalse ( curva altura - área )

y

el caudal instantáneo de rotura . El modelo calcula la altura de agua en el

caudales máximos por rotura ( tanto por sobrevertido como por tubificación , ) ,

y

aunque se tienen que tener en cuenta para iniciar el proceso de rotura por sobrevertido

brecha ( m ) , hb altura del fondo de la brecha ( m )

y

z talud de la brecha . El factor de sumergencia ks viene dado por la

éste . En esta transición se asume que el material entre el techo del orificio

y

la coronación de la balsa colapsa y es arrastrado , convirtiéndose el flujo en

que el material entre el techo del orificio y la coronación de la balsa colapsa

y

es arrastrado , convirtiéndose el flujo en un vertedero de lámina libre : Para

totalidad en un tiempo denominado tiempo de formación tf La anchura instantánea

y

la altura del fondo en la rotura vendrán dadas por donde bs es la

segundos ) , tf es el tiempo de formación de la brecha ( segundos )

y

ρ0 es el exponente de formación , que es interpretado como el grado de no

de no linealidad de crecimiento . Considerando ρ0=1 para roturas por sobrevertido

y

ρ0=4 para roturas por tubificación ( Fread , 1998 ; Sánchez - Romero , 2014

Esta metodología es similar a la utilizada por los modelos HEC - RAS , MIKE11

y

SMPDBK para el cálculo de los hidrogramas de rotura ( Greimann , 2004 ; Sánchez

, 2014 ) . Los valores para la anchura media , tiempo de formación

y

talud de brecha , son parámetros a definir por el usuario . Aunque existen muchas

1998 , 2004 ; Froehlich , 2008 ) , en España se utilizan las expresiones

y

recomendaciones propuestas por la “ Guía Técnica para la Clasificación de presas

Guía Técnica para la Clasificación de presas en función del riesgo potencial ”

y

“ Guía técnica para la elaboración de los Planes de Emergencia de presas ”

donde V es el volumen movilizable de la balsa en la rotura ( Hm3 )

y

hw es la altura de agua en la balsa en la rotura ( m )

diferencias para estos parámetros en rotura por sobrevertido o tubificación ( MacDonald

y

Langridge - Monopolis , 1984 ; Froehlich , 1995 ; Wahl , 1998 ) ,

y

entre presas ( gran volumen ) y balsas ( pequeños volúmenes ) ( Walder y

, 1995 ; Wahl , 1998 ) , y entre presas ( gran volumen )

y

balsas ( pequeños volúmenes ) ( Walder y O ’ Connor , 1997 ; Xu

y entre presas ( gran volumen ) y balsas ( pequeños volúmenes ) ( Walder

y

O ’ Connor , 1997 ; Xu and Zhang , 2009 ) , son utilizadas

sin embargo , en el presente se han generalizado los modelos bidimensionales ( 2D )

y

se empiezan a utilizar los tridimensionales ( 3D ) . En los modelos 1D el

de las secciones transversales cuando existen llanuras de inundación importantes

y

el cauce principal tiene cierta curvatura . Además la velocidad del agua no es constante

los modelos 2D - 3D , como un requerimiento menor en recursos de cálculo ,

y

por lo tanto menores tiempos de ejecución . Los programas de referencia 1 - D a

. Los programas de referencia 1 - D a nivel mundial han sido el FLDWAV

y

HEC - RAS , los dos de difusión libre , existiendo también otros modelos de

malla formada por una serie de celdas poligonales que representen la topografía

y

que se adapten al caso de estudio . Este proceso suele requerir de una mayor

, QUAD2D , RIVERFLOW2D , MIKE21 , SOBEK , RMA2 , HIVEL2D , INFOWORKS 2D

y

FLO 2D ) , en continuo desarrollo y mejora , que utilizan diferentes esquemas numéricos

, RMA2 , HIVEL2D , INFOWORKS 2D y FLO 2D ) , en continuo desarrollo

y

mejora , que utilizan diferentes esquemas numéricos para resolver los sistemas de

bidimensional de simulación de flujo turbulento en lámina libre en régimen no - permanente

y

de procesos medioambientales en hidráulica fluvial , siendo utilizado tanto para

por sobrevertido ( H2 ) , se supone la balsa a su nivel de coronación

y

un caudal de entrada máximo , considerando un fallo por sobrevertido por ineficacia

1 ) . Análisis de resultados ( Fase 4 ) Para la planificación

y

gestión en caso de rotura , es necesaria una descripción gráfica del fenómeno

intervienen en una posible emergencia , interpretar los posibles daños generados

y

poder determinar cuáles son las mejores estrategias para mitigarlos . Desde un punto

permite que , para cada uno de ellos , se obtengan la variación del calado

y

velocidad a lo largo del tiempo , lo que proporciona información sobre las posibles

se pueden combinar con datos georreferenciados , siendo posible una calificación

y

cuantificación de los posibles daños y afecciones ( Output 4 . 1c ) .

georreferenciados , siendo posible una calificación y cuantificación de los posibles daños

y

afecciones ( Output 4 . 1c ) . Sánchez - Romero ( 2014 )

proceso de rotura de una balsa , el vaciado se suele producir en minutos ,

y

en algunos casos en horas . Además no existe un caudal de aporte importante que

las inundaciones pueden generar daños sobre los elementos afectados , por calados

y

velocidades excesivas , son del orden de unos pocos minutos , e incluso inexistentes

afectada puede protegerse de la inundación , ya sea desplazándose a zonas seguras

y

/ o refugiándose en aquellas edificaciones que sean capaces de resistir la inundación

evaluar el grado de afecciones , fundamentalmente para las personas , vehículos

y

edificios , con el objetivo de poder priorizar las actuaciones en materia de autoprotección

peligrosidad para las personas está principalmente asociado con la velocidad , el calado

y

el producto de ambos ( Témez , 1992 ; Martínez - Gomariz et al .

aquella que cumpla con alguno de estos criterios v ≥ 1 m / s o

y

≥ 1 m o v · y ≥ 0 . 5 m2 / s ,

criterios v ≥ 1 m / s o y ≥ 1 m o v ·

y

≥ 0 . 5 m2 / s , siendo v velocidad , e y calado

· y ≥ 0 . 5 m2 / s , siendo v velocidad , e

y

calado , siendo el resto Zona no peligrosa . Según Nanía ( 1999 )

no ocasionarían ninguna pérdida humana , pero si cuantiosas pérdidas materiales

y

otros accidentes . Además , es importante evaluar también la afección a vehículos

importante evaluar también la afección a vehículos , para determinar su flotabilidad

y

su posible uso en caso de emergencia . La Agencia Catalana del Agua ( 2003

cumpla con alguno de estos criterios v ≥ 0 . 4 m / s ,

y

≥ 0 . 4 m y v · y ≥ 0 . 08 m2 /

v ≥ 0 . 4 m / s , y ≥ 0 . 4 m

y

v · y ≥ 0 . 08 m2 / s . Clausen y Clark

. 4 m / s , y ≥ 0 . 4 m y v ·

y

≥ 0 . 08 m2 / s . Clausen y Clark ( 1990 )

m y v · y ≥ 0 . 08 m2 / s . Clausen

y

Clark ( 1990 ) analizaron los daños en edificaciones producidos por la rotura de

los siguientes criterios : Zona de daños totales para v ≥ 2 m / s

y

v · y ≥ 7 m2 / s y Zona de daños parciales para v

: Zona de daños totales para v ≥ 2 m / s y v ·

y

≥ 7 m2 / s y Zona de daños parciales para v ≥ 2 m

v ≥ 2 m / s y v · y ≥ 7 m2 / s

y

Zona de daños parciales para v ≥ 2 m / s y v · y

m2 / s y Zona de daños parciales para v ≥ 2 m / s

y

v · y ≥ 3 m2 / s . Este mismo criterio es utilizado en

y Zona de daños parciales para v ≥ 2 m / s y v ·

y

≥ 3 m2 / s . Este mismo criterio es utilizado en la norma finlandesa

siguientes criterios de afección , que permiten a partir de los valores de velocidad

y

calado , determinar la posible afección tanto a personas , vehículos y edificios

velocidad y calado , determinar la posible afección tanto a personas , vehículos

y

edificios . Sánchez - Romero ( 2014 ) denomina y describe cada una de

a personas , vehículos y edificios . Sánchez - Romero ( 2014 ) denomina

y

describe cada una de estas zonas como : • Zona de afección leve ( v

• Zona de afección leve ( v ≤ 0 . 4 m / s e

y

≤ 0 . 4 m y v · y ≤ 0 . 08 m2 /

v ≤ 0 . 4 m / s e y ≤ 0 . 4 m

y

v · y ≤ 0 . 08 m2 / s ) , para esta zona

. 4 m / s e y ≤ 0 . 4 m y v ·

y

≤ 0 . 08 m2 / s ) , para esta zona el peligro para

de poco importancia . El daño producido sobre infraestructuras de comunicación

y

transporte son subsanables inmediatamente . • Zona de afección moderada ( v ≥

• Zona de afección moderada ( v ≥ 0 . 4 m / s o

y

≥ 0 . 4 m o v · y ≥ 0 . 08 m2 /

. 4 m / s o y ≥ 0 . 4 m o v ·

y

≥ 0 . 08 m2 / s ) . Para esta zona el peligro sobre

cortes en los sistemas de transporte , así como daños a estructuras de energía

y

comunicación . La inundación puede ser capaz de introducirse en las edificaciones

esta circunstancia para la planificación de posibles vías de escapatoria , avisos

y

corte de caminos y carreteras . • Zona de afección grave ( I )

para la planificación de posibles vías de escapatoria , avisos y corte de caminos

y

carreteras . • Zona de afección grave ( I ) ( v ≥ 1

Zona de afección grave ( I ) ( v ≥ 1 m / s o

y

≥ 1 m o v · y ≥ 0 . 5 m2 / s )

( v ≥ 1 m / s o y ≥ 1 m o v ·

y

≥ 0 . 5 m2 / s ) . Para esta zona existe afección grave

existe afección grave a vidas humanas . En los sistemas de transporte , energía

y

comunicación se pueden producir importantes daños , no pudiéndose reparar inmediatamente

con daño parcial en edificaciones ( II ) ( v ≥ 2 m / s

y

v · y ≥ 3 m2 / s ) . Mismas afecciones que para la

en edificaciones ( II ) ( v ≥ 2 m / s y v ·

y

≥ 3 m2 / s ) . Mismas afecciones que para la zona de afección

con daño total en edificaciones ( III ) ( v ≥ 2 m / s

y

v · y ≥ 7 m2 / s ) . Mismas afecciones que para la

en edificaciones ( III ) ( v ≥ 2 m / s y v ·

y

≥ 7 m2 / s ) . Mismas afecciones que para la zona de afección

s ) . Mismas afecciones que para la zona de afección grave ( I )

y

( II ) , siendo capaz de producir daños estructurales totales en las edificaciones

/ s ) en cada uno de los elementos de cálculo que componen ese corte

y

para cada tiempo de cálculo , se puede obtener el hidrograma de caudal ( m3

anteriormente se analizan de manera separada las hipótesis de rotura por tubificación

y

sobrevertido , definiéndose la relación entre caudales ( RQ ) como : Caso de

dique de las balsas analizadas consta de un ancho de coronación que oscila entre 4

y

5 metros , con taludes interiores entre 2H y 2 . 5H / 1V ,

de coronación que oscila entre 4 y 5 metros , con taludes interiores entre 2H

y

2 . 5H / 1V , taludes de terraplén entre 1 . 5H y 2H

2H y 2 . 5H / 1V , taludes de terraplén entre 1 . 5H

y

2H / 1V y taludes de desmonte entre 1H y 1 . 5H / 1V

5H / 1V , taludes de terraplén entre 1 . 5H y 2H / 1V

y

taludes de desmonte entre 1H y 1 . 5H / 1V . En la zona

terraplén entre 1 . 5H y 2H / 1V y taludes de desmonte entre 1H

y

1 . 5H / 1V . En la zona de terraplén está ubicada la caseta

prácticamente exclusiva a patologías asociadas a filtraciones a través del dique ,

y

concretamente las relacionadas con las vías preferentes inherentes a las estructuras

zonas . El número máximo de escenarios estudiado ha sido de 4 ( Balsas nº13

y

nº14 ) . Caso de estudio . Modelos de cálculo . En general ,

zonas de ladera con fuertes pendientes . Estas zonas confluyen a zonas cultivadas

y

/ o pequeñas vaguadas que terminan desembocando en el río Vinalopó , el cual constituye

desembocando en el río Vinalopó , el cual constituye el cauce principal de la cuenca

y

que presenta un cauce seco en la mayor parte de su recorrido . Para algunas

de su recorrido . Para algunas balsas el volumen de rotura termina almacenándose

y

drenándose lentamente o queda acumulado en una cuenca endorreica ( balsa nº10 )

. Se puede concluir que el cauce del río Vinalopó tiene capacidad para encauzar

y

desaguar la onda de rotura , estableciéndose las zonas de desembocadura al río

finalizado en las zonas de acumulación . Como condición inicial a toda la malla

y

para todos los casos se ha impuesto una situación de lecho seco . Para

un tamaño mínimo de malla de 1 metro en las zonas con topografía más compleja

y

con mayor densidad de puntos afectados . Los modelos realizados presentan una media

de elementos de cálculo de 1819047 , un máximo de 3363446 ( balsa nº12 )

y

un mínimo de 716942 ( balsa nº4 ) . En la Figura 5a y 5b ,

nº12 ) y un mínimo de 716942 ( balsa nº4 ) . En la Figura 5a

y

5b , se muestra ejemplo de la malla de cálculo conseguida para zonas donde era

realidad física existente . Para zonas puntuales donde roturas de motas , terraplenes

y

/ o márgenes pudieran aumentar el número de elementos afectados , se realizaron

afectados , se realizaron diferentes simulaciones modificando la geometría de la malla

y

analizando los posibles efectos que estas modificaciones pudieran tener sobre los

de variación para el coeficiente de Manning ha sido entre valores de 0 . 015

y

0 . 12 , estando los valores predominantemente en un rango entre 0 . 04

, ya que este proceso hidrológico disminuiría los volúmenes de propagación ,

y

por lo tanto reduciría la extensión tanto de las zonas afectadas como de las afecciones

objetivo de este tipo de análisis es la obtención de información para la gestión

y

planificación en caso de emergencia , lo que implica que los resultados finales

referirse a situaciones pésimas desde el punto de vista de las afecciones . RESULTADOS

Y

DISCUSIÓN En la Tabla 1 se muestran los resultados obtenidos para cada una de las hipótesis

consideradas en los casos de estudio . Se observa que los volúmenes movilizables

y

las alturas máximas de agua son similares para las hipótesis de rotura consideradas

de agua son similares para las hipótesis de rotura consideradas ( tubificación

y

sobrevertido ) , y por tanto los parámetros calculados para el tiempo de formación

similares para las hipótesis de rotura consideradas ( tubificación y sobrevertido ) ,

y

por tanto los parámetros calculados para el tiempo de formación y anchura media

sobrevertido ) , y por tanto los parámetros calculados para el tiempo de formación

y

anchura media de brecha son muy parecidos . Los caudales obtenidos para tubificación

. En la Figura 6a se muestran a modo de ejemplo , los hidrogramas de caudal

y

variación de la altura en el interior del vaso para la balsa nº1 , siendo

, siendo para el resto de balsas análogos a los mostrados . En la Figura 6b

y

Figura 6c se presenta la ratio de atenuación ( A ) y el tiempo de

En la Figura 6b y Figura 6c se presenta la ratio de atenuación ( A )

y

el tiempo de llegada de la onda de rotura , respectivamente , solo para los

tiempo de llegada de la onda de rotura , respectivamente , solo para los casos

y

escenarios donde ha sido posible realizar cortes en la malla sobre flujos principales

la Figura 6d , se muestran para estos casos la relación de caudal ( RQ )

y

la relación de tiempos de llegada ( Rt ) . En la Figura 6a se

el porcentaje de estos volúmenes respecto al volumen inicial oscila entre un 19

y

un 29 % , siendo el promedio del 24 % . Para sobrevertido oscila entre

% , siendo el promedio del 24 % . Para sobrevertido oscila entre un 50

y

un 62 % , siendo el promedio de un 56 % . En un proceso

las características de la rotura , de la topografía ( pendiente , obstrucciones )

y

de la rugosidad considerada ( coeficiente de Manning ) . En todos los casos analizados

cercanas a la balsa , la zona afectada es superior en la hipótesis por tubificación

y

para las zonas más alejadas la hipótesis de rotura por sobrevertido presenta mayores

distancia mínima recorrida por la onda de rotura ha sido de 1 . 25 km

y

la máxima de aproximadamente de 2 . 6 km , dependiendo esta distancia de las

permite definir sobre el terreno las zonas inundadas con cualquier grado de afección ,

y

determinar las estrategias más convenientes para el aviso a la población . En la

casos analizados , se observa cómo las zonas de mayores afecciones ( daño total

y

parcial a las edificaciones ) son muy reducidas y se emplazan junto a las balsas

mayores afecciones ( daño total y parcial a las edificaciones ) son muy reducidas

y

se emplazan junto a las balsas y en pequeñas vaguadas donde se concentra el flujo

parcial a las edificaciones ) son muy reducidas y se emplazan junto a las balsas

y

en pequeñas vaguadas donde se concentra el flujo . A medida que nos alejamos de

grado de afección va disminuyendo , apareciendo zonas con grados de afección leve

y

moderado . A partir de la ubicación de los posibles elementos afectados , es posible

de los posibles elementos afectados , es posible obtener la variación del calado

y

velocidad en cualquier punto mediante las herramientas que disponen los modelos ,

y

por tanto también la variación del grado de afección y su duración ( Figura 7c

disponen los modelos , y por tanto también la variación del grado de afección

y

su duración ( Figura 7c y Figura 7d ) . En el ejemplo mostrado ,

por tanto también la variación del grado de afección y su duración ( Figura 7c

y

Figura 7d ) . En el ejemplo mostrado , se observa que para el escenario

siguiente manera : 5 minutos de afección leve , 5 minutos de afección moderada

y

10 minutos de afección grave . CONCLUSIONES Los trabajos desarrollados han permitido

debido a la diferente naturaleza del comienzo del proceso de rotura en sobrevertido

y

tubificación , las alturas de agua son mayores por tubificación , obteniéndose

de rotura se produzca en el momento que se ha formado totalmente la brecha ,

y

que por lo tanto no existan diferencias significativas en este aspecto entre sobrevertido

por lo tanto no existan diferencias significativas en este aspecto entre sobrevertido

y

tubificación . Por otra parte , la rotura por sobrevertido desagua mayores volúmenes

volúmenes en los momentos iniciales , lo que implica menores tiempos de llegada

y

una menor atenuación a lo largo de la propagación de la onda de rotura .

rotura por sobrevertido comporta mayores afecciones en zonas alejadas a la balsa

y

menores tiempos de llegada a lo largo de toda la zona afectada . Aunque

estudios en balsas de tierra deberían establecer todos los escenarios posibles ,

y

en cada uno de estos escenarios analizar dos hipótesis de rotura ( tubificación

y

sobrevertido ) , determinando las circunstancias más desfavorables para la zona

sobrevertido ) , determinando las circunstancias más desfavorables para la zona afectada

y

por tanto las afecciones pésimas en los elementos existentes . La categorización

categorización de los diferentes grados de afección , permite la planificación

y

gestión ante una posible emergencia , pudiendo priorizar las actuaciones a llevar

ante una posible emergencia , pudiendo priorizar las actuaciones a llevar a cabo

y

ayudando a la toma de decisiones . La elaboración de mapas de inundación ,

inundación , donde se analicen adecuadamente las consecuencias en caso de rotura

y

con el suficiente detalle para que puedan adoptarse medidas que aminoren y / o anulen

rotura y con el suficiente detalle para que puedan adoptarse medidas que aminoren

y

/ o anulen los daños , en especial de vidas humanas , deben ser una

Algunos investigadores se han enfocado en el tratamiento que realizan los libros textos

y

su influencia en el aprendizaje [ 1 , 6 ] , otros han estudiado las

en el aprendizaje [ 1 , 6 ] , otros han estudiado las ideas previas

y

el razonamiento de los estudiantes [ 7 , 21 , 23 ] y otros han

ideas previas y el razonamiento de los estudiantes [ 7 , 21 , 23 ]

y

otros han sugerido que la discusión grupal podría ayudar a clarificar las preconcepciones

dificultades que tienen los alumnos para el aprendizaje de la fuerza de fricción

y

ponen de manifiesto la relevancia de investigar los procesos de enseñanza y aprendizaje

fricción y ponen de manifiesto la relevancia de investigar los procesos de enseñanza

y

aprendizaje sobre dicha fuerza . La complejidad de la actividad físico - matemática

complejidad de la actividad físico - matemática ( en términos de prácticas , contenidos

y

procesos , tanto físicos como matemáticos ) que realiza un grupo de estudiantes

fricción ya que permite tanto ( 1 ) explicar mejor las causas de los errores

y

dificultades en el aprendizaje de los alumnos , como ( 2 ) sugerir recomendaciones

supuesto de que la fuerza de fricción posee dos componentes indisociables , uno físico

y

otro matemático , que son esenciales para la comprensión adecuada de los fenómenos

fricción ” . La importancia del componente matemático nos ha llevado a utilizar

y

adaptar un marco teórico y metodológico que provienen de la Didáctica de la Matemática

del componente matemático nos ha llevado a utilizar y adaptar un marco teórico

y

metodológico que provienen de la Didáctica de la Matemática , el Enfoque Ontosemiótico

EOS , en adelante ) [ 11 ] . Este enfoque se está desarrollando

y

utilizando principalmente como referente teórico para la investigación en educación

principalmente como referente teórico para la investigación en educación matemática

y

, recientemente , algunos investigadores lo han aplicado a la investigación de los

investigadores lo han aplicado a la investigación de los procesos de enseñanza

y

aprendizaje de la Física [ 2 , 8 , 17 ] . En esta

prácticas físico - matemáticas , configuraciones epistémicas de objetos primarios

y

procesos activados en dichas prácticas . El término diagrama [ 27 ] es empleado

EOS , éstos pueden ser de tipo institucional si son construidos por un experto (

y

se corresponden con el llamado diagrama de cuerpo libre ) o personal si son construidos

se presenta una revisión bibliográfica de la investigación sobre la enseñanza

y

aprendizaje de la fuerza de fricción . En la tercera , se presenta el objetivo de

investigación , mientras que en la cuarta se describe la metodología utilizada

y

se explican algunos elementos teóricos provenientes del EOS , que se han utilizado

problema del corredor , en particular las nociones de configuración , epistémica

y

cognitiva , de objetos primarios . En la quinta sección se analizan y clasifican

epistémica y cognitiva , de objetos primarios . En la quinta sección se analizan

y

clasifican las respuestas de los participantes −se realiza una primera clasificación

situación macroscópica que ya está presente en los trabajos de Leonardo da Vinci

y

en el de los físicos Guillaume Amontons y Charles Augustin Coulomb [ 18 ] .

en los trabajos de Leonardo da Vinci y en el de los físicos Guillaume Amontons

y

Charles Augustin Coulomb [ 18 ] . Estos autores señalan , además , que lo

se sabe actualmente sobre la fricción no se ha incorporado en los libros de texto

y

que , por otra parte , estos no tienen en cuenta las ideas previas de los estudiantes

dificultades para identificar el punto de aplicación de la fuerza de rozamiento

y

tienen cierta confusión acerca de cuál es el cuerpo que ejerce la fuerza , no

fuerza , no especifican el sistema de referencia cuando describen algún movimiento

y

no están habituados además al empleo de la tercera ley de Newton cuando consideran

importante la fórmula de la fuerza de rozamiento estático que la segunda ley de Newton ,

y

no encuentran dificultad en negar el cumplimiento de ésta con tal de salvar aquélla

estudiantes universitarios sobre la fuerza de rozamiento en el movimiento circular ,

y

se señala la contradicción que hay entre la interpretación correcta ( dirigida

tendencia del movimiento , en el sentido de proceso , desarrollo o propensión )

y

la de muchos estudiantes que consideran que la fuerza de rozamiento siempre está

no son prototípicas donde se presenta la fuerza de rozamiento , en [ 14 ]

y

[ 21 ] reportan que muchos alumnos argumentan que la acción o resistencia de un

los estudiantes cuando se enfrentan a situaciones de rozamiento estático [ 23 ] ,

y

se han formulado además propuestas didácticas para la enseñanza de la fuerza de

aplicación del concepto de fuerza de rozamiento al análisis de fenómenos reales [ 21 ]

y

presentar a los estudiantes secuencias didácticas a través de experimentos en distintas

los textos describen leyes empíricas para la fuerza de rozamiento entre sólidos

y

clasifican la fuerza de rozamiento en dos tipos : en estática fe , entre superficies

de rozamiento en dos tipos : en estática fe , entre superficies en reposo ,

y

cinética fc , entre superficies en movimiento relativo . Con relación al sentido

. Las fuerzas de rozamiento se oponen automáticamente a este movimiento relativo

y

nunca contribuyen a él . Un ejemplo de los libros de textos analizados por

la fuerza de rozamiento : ( 1 ) es casi independiente del área de contacto

y

( 2 ) es proporcional a la magnitud de la fuerza normal ( esta última

( esta última es debida a las propiedades elásticas de los cuerpos en contacto

y

, por lo general , es igual al peso del cuerpo o a alguna de sus componentes

Por otra parte , el cociente entre la magnitud de la fuerza máxima de rozamiento estática

y

la magnitud de la fuerza normal se llama coeficiente de rozamiento estático µe

normal se llama coeficiente de rozamiento estático µe de las superficies implicadas

y

se establece que fe ≤ µeN , donde el signo de igualdad se obtiene cuando

cinética fc entre superficies secas no lubricadas también cumple con las mismas leyes

y

es independiente de la velocidad relativa entre los cuerpos . La relación entre

los cuerpos . La relación entre la magnitud de la fuerza de rozamiento cinético

y

la magnitud de la fuerza normal se llama coeficiente de rozamiento cinético µc

y

se establece que fc = µcN . En general , los valores de estos coeficientes dependen

valores de estos coeficientes dependen del material de las superficies en contacto

y

se tiene que µe >µc , por lo que la máxima fuerza de rozamiento estática

OBETIVO DE LA INVESTIGACIÓN En las explicaciones de un profesor en su clase

y

en los textos de Física , se emplean comúnmente representaciones visuales , modelos

actividad físico - matemática donde se usan conceptos , propiedades , procedimientos

y

argumentos , y se ponen en juego procesos diversos , entre otros los de significación

matemática donde se usan conceptos , propiedades , procedimientos y argumentos ,

y

se ponen en juego procesos diversos , entre otros los de significación y visualización

argumentos , y se ponen en juego procesos diversos , entre otros los de significación

y

visualización . Si se utiliza la metáfora del iceberg , la producción ostensiva

actividad físico - matemática ( caracterizada en términos de prácticas , contenidos

y

procesos , tanto físicos como matemáticos ) realizada para resolver un problema

matemáticos ) realizada para resolver un problema en el que se presentan fenómenos reales

y

en cuya solución los diagramas , en los que interviene la fuerza de rozamiento ,

dos conceptualizaciones [ 3 ] , como fuerza de reacción donde juega un rol motriz

y

como fuerza que se opone al movimiento , donde juega un papel resistivo . Los

) donde la conceptualización más adecuada es la de reacción . MARCO TEÓRICO

Y

METODOLÓGICO En esta investigación , de acuerdo con el EOS [ 10 - 11

con el propósito de resolver el problema propuesto , en términos de prácticas ,

y

configuraciones epistémicas y cognitivas de objetos primarios , así como de procesos

problema propuesto , en términos de prácticas , y configuraciones epistémicas

y

cognitivas de objetos primarios , así como de procesos activados en dichas prácticas

práctica que permite resolver el problema propuesto , el alumno necesita interpretar

y

realizar ciertas representaciones , las cuales son la parte ostensiva de una serie

, las cuales son la parte ostensiva de una serie de definiciones , proposiciones

y

procedimientos que intervienen en la elaboración de argumentos que permiten la solución

desencadena la práctica ( definiciones , proposiciones , procedimientos , representaciones

y

argumentos utilizados ) se le llama configuración de objetos primarios . Se distingue

se trata de la configuración que realizaría un experto ( perspectiva institucional )

y

configuraciones cognitivas cuando se trata de las que realiza cada alumno ( perspectiva

argumento , un proceso de argumentación , etc . Por tanto , además de prácticas

y

configuraciones , conviene tener en cuenta , para analizar la actividad físico -

, en la física escolar , se puede distinguir entre práctica , objetos primarios

y

procesos si consideramos a la actividad involucrada en la tarea de calcular la velocidad

y

la aceleración de un cuerpo que tiene En dicha tarea , los estudiantes llevan a

una secuencia de acciones ( práctica ) tales como ( 1 ) la lectura del enunciado

y

( 2 ) el cálculo de la primera y segunda derivada mediante la regla de la

( 1 ) la lectura del enunciado y ( 2 ) el cálculo de la primera

y

segunda derivada mediante la regla de la cadena para derivar funciones compuestas

problema del corredor : respuestas correctas , parcialmente correctas e incorrectas )

y

cualitativas ( tipo de solución del problema propuesto ) . En el estudio cualitativo

- 20 ] , que permite describir sistemáticamente las configuraciones epistémicas

y

cognitivas ( situación - problema , representaciones , conceptos , propiedades ,

situación - problema , representaciones , conceptos , propiedades , procedimientos

y

argumentos ) que intervienen en las prácticas que permiten la solución de la tarea

referencia ( solución experta ) , después la configuración cognitiva de los alumnos

y

por último se hacen comparaciones entre ambas . La metodología empleada tomó

que en el EOS se designan como el foco , el fin , la generalizabilidad

y

el nivel de la investigación . El foco fue epistémico ( configuraciones epistémicas

nivel de la investigación . El foco fue epistémico ( configuraciones epistémicas )

y

también cognitivo ( configuraciones cognitivas ) . El fin estuvo relacionado a la

contexto del corredor . Participantes Los participantes fueron 71 alumnos de entre 19

y

20 años de edad que cursaban el tercer semestre en la Facultad de Ciencias Químicas

al texto la Figura 1 , la cual muestra posiciones diferentes del corredor A , B

y

C , así como comentarios sobre ellas . La modificación del problema original fue

bloques . En el primer bloque se investiga el sentido de la fuerza de rozamiento

y

la resolución del problema , preguntas 1 y 16 respectivamente . El segundo ( preguntas

el sentido de la fuerza de rozamiento y la resolución del problema , preguntas 1

y

16 respectivamente . El segundo ( preguntas 2 a 4 ) indaga las fuerzas que actúan

) indaga las fuerzas que actúan cuando el corredor se encuentra en A , B

y

C . El tercero ( preguntas 5 a 7 ) averigua el efecto de la fuerza

7 ) averigua el efecto de la fuerza de fricción en el movimiento del corredor

y

la representación de las fuerzas en las posiciones A y C del corredor en la Figura

el movimiento del corredor y la representación de las fuerzas en las posiciones A

y

C del corredor en la Figura 1 . El cuarto ( preguntas 17 ) indaga la primera

empleo de un marco de referencia . En el quinto ( preguntas 8 , 9

y

15 ) se investiga la segunda ley de Newton mediante la condición de aceleración

mediante la condición de aceleración máxima . Mediante el sexto ( preguntas 12 , 13

y

14 ) se investiga la tercera ley de Newton al considerar la fuerza de fricción como

de fricción como fuerza de reacción . Y con el último bloque ( preguntas 10

y

11 ) se investiga la consideración de la propiedad P4 ( Tabla 1 ) mediante

consideración de la propiedad P4 ( Tabla 1 ) mediante la condición de resbalamiento

y

del valor máximo de la fuerza de fricción estática . Práctica , configuración

máximo de la fuerza de fricción estática . Práctica , configuración epistémica

y

procesos Para resolver la tarea el experto debe realizar una práctica físico -

práctica físico - matemática que consiste en la lectura del texto de la tarea

y

en la producción de un texto con las respuestas a las distintas cuestiones planteadas

producción de un texto con las respuestas a las distintas cuestiones planteadas ,

y

que será el resultado de una serie de acciones sujetas a reglas físico - matemáticas

movimiento pueda ser modelado como la combinación de un movimiento de traslación

y

otro de rotación . El resolutor experto realiza un tratamiento del problema similar

experto modela al corredor como una partícula situada en el contacto con la pista

y

representa mediante un diagrama , desde la perspectiva de un marco de referencia

un marco de referencia en reposo , a tres fuerzas ( peso , fuerza normal

y

rozamiento ) actuando sobre dicha partícula . En la Tabla 1 se presenta la configuración

que proponen [ 12 ] : a ) Comunicación de la forma , sus componentes

y

estructura , de objetos espaciales , o bien de objetos imaginados ( pensados o ideales

, de objetos espaciales , o bien de objetos imaginados ( pensados o ideales )

y

b ) Comunicación de la posición relativa de objetos en el espacio . En efecto

La tarea comunica información visual . Se explicitan los componentes ( corredor

y

suelo ) , forma ( pista plana ) , se sugiere movimiento ( mediante el

, se sugiere movimiento ( mediante el dibujo de varias posiciones del corredor )

y

se espera que el alumno realice un proceso de idealización y entienda que el corredor

posiciones del corredor ) y se espera que el alumno realice un proceso de idealización

y

entienda que el corredor se puede considerar como un punto ( objeto ideal ) .

posición relativa de objetos en el espacio ( por ejemplo , el corredor está arriba

y

el suelo abajo ) . Para resolver el problema , el alumno debe especificar un

, la mayoría de los conceptos , procedimientos , representaciones , propiedades

y

argumentos de esta configuración epistémica también se pueden considerar como

preguntas del cuestionario . Según el bloque 2 , las preguntas 2 , 3

y

4 , el 59 % , 39 % y 42 % de los estudiantes respectivamente

, las preguntas 2 , 3 y 4 , el 59 % , 39 %

y

42 % de los estudiantes respectivamente señalaron correctamente a las fuerzas presentes

respectivamente señalaron correctamente a las fuerzas presentes en las posiciones A , B

y

C del corredor ( Figura 1 ) . Sin embargo , en el bloque 1

( Figura 1 ) . Sin embargo , en el bloque 1 , preguntas 1

y

16 , el 17 % de los alumnos representó correctamente el sentido de la fuerza

de los alumnos representó correctamente el sentido de la fuerza de fricción en B

y

C , el resto lo dibujó en sentido opuesto al movimiento del corredor .

preguntas 6 el 46 % representó adecuadamente a las fuerzas en la posición A ,

y

la pregunta 7 verificó el resultado de las preguntas 1 y 16 . En

la posición A , y la pregunta 7 verificó el resultado de las preguntas 1

y

16 . En el bloque 4 , pregunta 17 , sólo el 52 %

pocos los estudiantes que consideraron la segunda ley de Newton . En las preguntas 8

y

15 , cerca del 30 % señaló una relación entre el empuje y el incremento

preguntas 8 y 15 , cerca del 30 % señaló una relación entre el empuje

y

el incremento de la aceleración , el resto sugirió que una fuerza de empuje constante

aceleración máxima . En relación con el bloque 6 , las respuestas a las preguntas 12

y

13 muestran que aproximadamente el 60 % de los estudiantes consideró necesario emplear

fuerza de reacción al empuje . Los resultados de los bloques 1 , 5

y

6 nos sugieren que un porcentaje relevante de los alumnos investigados considera

fuerza de rozamiento , como aquella que dificulta el movimiento de los cuerpos ,

y

no como fuerza de reacción , pues muchos consideran a dicha fuerza de magnitud constante

fuerza de reacción , pues muchos consideran a dicha fuerza de magnitud constante

y

de sentido opuesto el movimiento del cuerpo del corredor . Se trata de un

pregunta 10 señalan que el 38 % sugirió correctamente la condición de resbalamiento

y

en la pregunta 11 el 63 % identificó la relación entre la fuerza de fricción

y

el peso del corredor . Las cifras anteriores nos muestran la manera en que

de los estudiantes fueron clasificadas en cuatro grupos : ( 1 ) los que idealizan

y

resuelven correctamente , ( 2 ) los que idealizan y no resuelven correctamente ,

1 ) los que idealizan y resuelven correctamente , ( 2 ) los que idealizan

y

no resuelven correctamente , ( 3 ) los que no idealizan y resuelven correctamente

los que idealizan y no resuelven correctamente , ( 3 ) los que no idealizan

y

resuelven correctamente y ( 4 ) los que no idealizan y no resuelven correctamente

no resuelven correctamente , ( 3 ) los que no idealizan y resuelven correctamente

y

( 4 ) los que no idealizan y no resuelven correctamente , Figura 3 .

los que no idealizan y resuelven correctamente y ( 4 ) los que no idealizan

y

no resuelven correctamente , Figura 3 . El problema sólo fue resuelto correctamente

clasificación , a continuación se presentan los casos de los alumnos Luis , Elvira

y

Alberto para describir dos aspectos relacionados con la fuerza de fricción donde

estudiantes tuvieron más dificultades : ( i ) el sentido de la fuerza de rozamiento

y

( ii ) la maximización de la aceleración del corredor . En ambos casos se

preguntas 1 - 4 del cuestionario , los estudiantes asignaron correctamente la dirección

y

el sentido a la fuerza normal y al peso del corredor , pero no sucedió

estudiantes asignaron correctamente la dirección y el sentido a la fuerza normal

y

al peso del corredor , pero no sucedió lo mismo con el sentido de la

, el caso de Luis ; y el segundo que sí realiza el proceso de idealización

y

resuelve la tarea de manera adecuada , el caso de Elvira . El caso

alumno interpreta la fuerza de rozamiento como aquella que se opone al movimiento

y

toma en cuenta cuatro fuerzas : el peso , la fuerza normal , el empuje

cuenta cuatro fuerzas : el peso , la fuerza normal , el empuje del corredor

y

el rozamiento ( de sentido incorrecto ) . Así mismo , realiza una representación

haber idealizado , pues primero hace un proceso de esquematización del corredor

y

después un proceso de idealización que le lleva a “ concentrar ” al corredor

incorrecta de este alumno . En su representación diagramática ( parte central

y

derecha de la Figura 4 ) puede notarse que las fuerzas presentes en la situación

notarse que las fuerzas presentes en la situación son representadas mediante flechas (

y

letras que las acompañan ) . En ella , se pueden observar dos aspectos importantes

parte del alumno ( dibuja la fuerza de rozamiento opuesta al movimiento del corredor

y

cuatro fuerzas actuando sobre él ) , también emplea de manera inadecuada P2 al

atribuir la aceleración del corredor a la acción directa de la fuerza de empuje

y

no a la fuerza de rozamiento . Así mismo no emplea P3 y P4

empuje y no a la fuerza de rozamiento . Así mismo no emplea P3

y

P4 , ya que no considera la fuerza de rozamiento como reacción al empuje que

considera la fuerza de rozamiento como reacción al empuje que impulsa al corredor

y

no expresa la condición de resbalamiento . Finalmente , estos aspectos lo conducen

en el punto que representa el punto de contacto del pie con la pista )

y

no se equivoca al asignar el sentido a la fuerza de rozamiento . Se puede

“ E – F ” se refiere a que ambas fuerzas son iguales en magnitud

y

que no necesariamente está considerando que sean fuerzas vectoriales que se cancelen

describe el caso del alumno Alberto , el cual no realiza el proceso de idealización

y

no resuelve correctamente el problema . El caso de Alberto Alberto también pertenece

corredor se mantenía constante durante el movimiento ( pregunta 8 , Figura 2 ) ,

y

también al 63 % que señalaron la relación entre la fuerza de fricción y el

, y también al 63 % que señalaron la relación entre la fuerza de fricción

y

el peso del corredor , pregunta 11 . Por cuestiones de espacio se presenta su

( Figura 6 ) puede notarse la representación de algunas fuerzas mediante flechas

y

símbolos que las acompañan ( la fuerza normal N , el peso del corredor Fg

y

rozamiento FFric ) , sin embargo , en lugar de representar la fuerza de empuje

erróneo a la fuerza de fricción al suponerla contraria al movimiento del corredor

y

como aquella que dificulta su movimiento ( como en el caso del estudiante Alberto

cual utilizó para representar que la fuerza normal cancela al peso del corredor ,

y

concluir que el movimiento se lleva a cabo de manera horizontal al dibujar el eje

que la representación del corredor como un bloque tiene relación con el sentido

y

la magnitud constante que el estudiante atribuye a la fuerza de rozamiento . Esto

rozamiento como una fuerza que dificulta el movimiento del corredor ( sentido opuesto )

y

también el papel que juega en la expresión de la máxima aceleración “ A ”

que acelera al corredor como la resultante de la suma de la fuerza de empuje

y

la fuerza de rozamiento ( ambas constantes ) podría explicar también por qué el

corredor obtenga aceleración máxima . Dicha condición señala que la fuerza de empuje

y

la de fricción se incrementan y se cancelan mutuamente durante el movimiento hasta

Dicha condición señala que la fuerza de empuje y la de fricción se incrementan

y

se cancelan mutuamente durante el movimiento hasta que la fuerza de fricción es

cual , sólo algunos estudiantes señalaron una relación entre la fuerza de empuje

y

la aceleración . Cabe destacar que la mayoría de los estudiantes no sugirió una

mayoría de los estudiantes no sugirió una relación entre la fuerza de fricción

y

la velocidad . Una de las causas de sus conjeturas erróneas sobre ( 1 )

. Una de las causas de sus conjeturas erróneas sobre ( 1 ) el sentido

y

( 2 ) el incremento de la magnitud de la fuerza de rozamiento estática hasta

considerar ( a ) que la fuerza de rozamiento obstaculiza el movimiento del corredor

y

a ( b ) atribuir al corredor un movimiento con aceleración constante . Este resultado

aceleración constante . Este resultado es concordante con las investigaciones de [ 1 ]

y

[ 6 ] que señalan las dificultades que fomentan los libros de texto que utilizan

idealización , que lleva en un primer momento a representar al corredor de manera icónica

y

después como un punto al cual se le atribuye la masa del corredor y sobre

icónica y después como un punto al cual se le atribuye la masa del corredor

y

sobre el cual se aplican todas las fuerzas . Según los resultados del cuestionario

le asigna un sentido adecuado ) . El papel relevante de este proceso por una parte

y

el uso de representaciones prototípicas en forma de bloque en muchos libros , por

fuerza de rozamiento . También el hecho de idealizar a los cuerpos como bloques

y

no como partículas podría ser una explicación de los análisis de [ 5 ] ,

Newton o de la ley empírica de fricción ( P4 en la Tabla 1 , Tabla 2

y

Tabla 3 ) en la solución de la tarea del corredor , y en general de

Tabla 2 y Tabla 3 ) en la solución de la tarea del corredor ,

y

en general de cualquier otra tarea en el contexto de la mecánica , es lo que

podría explicar que los alumnos no terminen de dar sentido a la solución del problema

y

tengan dificultades para explicar su respuesta . Dichas dificultades se observan

8 - 15 del cuestionario . Por último , con base en las observaciones anteriores

y

con el hecho de que la mayoría de los autores de los libros de texto

5 de septiembre de 2017 ) Hibridación entre diferentes tipos de chiles

y

estimación de la heterosis para rendimiento y calidad de fruto Introducción

Hibridación entre diferentes tipos de chiles y estimación de la heterosis para rendimiento

y

calidad de fruto Introducción El chile ( Capsicum annuum L .

hortícolas a nivel mundial , su importancia radica en que es una hortaliza muy cotizada

y

consumida mundialmente ( La - borde - Cansino y Pozo - Campodónico , 1984 )

que es una hortaliza muy cotizada y consumida mundialmente ( La - borde - Cansino

y

Pozo - Campodónico , 1984 ) . En 2015 la producción mundial de chile fue

lugar en exportaciones a países como Estados Unidos , Japón , Canadá , Reino Unido

y

Alemania ( SAGARPA , 2015 ) . Sin embargo , cada día es más difícil cubrir

las crecientes demandas del mercado , por lo que se requieren genotipos innovadores

y

de alto po tencial genético que se adapten a condiciones de ambiente controlado

genético que se adapten a condiciones de ambiente controlado , garantizando alta calidad

y

rendimiento de fruto y el cultivo de hortalizas tiene cada vez más relevancia debido

condiciones de ambiente controlado , garantizando alta calidad y rendimiento de fruto

y

el cultivo de hortalizas tiene cada vez más relevancia debido a la necesidad de

hortalizas tiene cada vez más relevancia debido a la necesidad de diversificación

y

de mejorar la calidad de los productos alimenticios ( Rodríguez et al . , 2007

intensiva en invernadero , puesto que no existen genotipos de Jalapeño , Mirador

y

/ o Serrano de estas características y es mediante esta tecnología donde se obtienen

existen genotipos de Jalapeño , Mirador y / o Serrano de estas características

y

es mediante esta tecnología donde se obtienen mayores rendimientos y mejor calidad

características y es mediante esta tecnología donde se obtienen mayores rendimientos

y

mejor calidad de frutos . De esta forma , además se puede producir cuando no

, 2003 ) . México es el centro de origen de C . annuum

y

posee una amplia diversidad que oscila entre 50 tipos diferentes , entre los que

que oscila entre 50 tipos diferentes , entre los que destacan Jalapeño , Serrano

y

Pimiento . Se considera como una de las especies que exhibe mayor variabilidad de

presentando amplia diversidad de formas , tamaños , colores , rangos de maduración

y

grado de pungencia ( Valadez , 1998 ) , lo cual constituye un excelente recurso

constituye un excelente recurso para el mejoramiento genético ( Laborde - Cansino

y

PozoCampodónico , 1984 ; Bosland , 1994 ) . Actualmente el rendimiento promedio

mejoramiento genético puede ser útil en la obtención de genotipos de alto rendimiento

y

calidad de fruto , aprovechando la capacidad combinatoria y heterosis ( Pérez -

genotipos de alto rendimiento y calidad de fruto , aprovechando la capacidad combinatoria

y

heterosis ( Pérez - Grajales et al . , 2009 ) . La explotación

incrementar el rendimiento u otros caracteres económicos importantes ( Seneviratne

y

Kannangara , 2004 ) . Se considera que en Capsicum la heterosis es alta (

, 2004 ) . Se considera que en Capsicum la heterosis es alta ( Sousa

y

Maluf , 2003 ) , esto es debido a la amplia variedad de genotipos que

semicultivadas como silvestres , mismas que pueden aprovecharse para formar híbridos locales

y

nacionales . Además , actualmente un porcentaje importante de semilla hibrida que

cada vez más relevante , ya que no existen híbridos de chile Jalapeño , Mirador

y

/ o Serrano de alto rendimiento y alta calidad nutricional para producción en ambientes

existen híbridos de chile Jalapeño , Mirador y / o Serrano de alto rendimiento

y

alta calidad nutricional para producción en ambientes protegidos . Por ello se planteó

obtener híbridos intraespe cíficos , con potencial para su producción en invernadero

y

estimar la heterosis y heterobeltiosis en los híbridos formados . Material y métodos

cíficos , con potencial para su producción en invernadero y estimar la heterosis

y

heterobeltiosis en los híbridos formados . Material y métodos El presente trabajo

invernadero y estimar la heterosis y heterobeltiosis en los híbridos formados . Material

y

métodos El presente trabajo se realizó en la Universidad Autónoma Agraria Antonio

cuenta con extractores , calefactores , estación meteorológica , pared húmeda

y

control de temperatura ; registrando mínimas de 18°C y máximas de 36°C y una

meteorológica , pared húmeda y control de temperatura ; registrando mínimas de 18°C

y

máximas de 36°C y una humedad relativa promedio de 60 % . Formación de

húmeda y control de temperatura ; registrando mínimas de 18°C y máximas de 36°C

y

una humedad relativa promedio de 60 % . Formación de híbridos Los progenitores

de germinación Tur ba de Sphagnum Premier ( Angeles Millwork & Lumber Co USA )

y

perlita en una proporción 80 : 20 respectivamente . A los 45 días después de

después de la siembra se realizó el trasplante en camas con acolchado plástico negro

y

cintilla para riego localizado , las camas con 25 cm de altura y una separación

plástico negro y cintilla para riego localizado , las camas con 25 cm de altura

y

una separación de 1 , 60 m ; 30 cm entre plantas . El riego

Las polinizaciones se realizaron diariamente de las 7 : 00 a 10 : 00 h

y

cada flor polinizada fue cubierta , para evitar polinización no deseada , y etiquetada

00 h y cada flor polinizada fue cubierta , para evitar polinización no deseada ,

y

etiquetada in - dicando los progenitores utilizados . Cuando los frutos resultantes

frutos resultantes de los cruces alcanzaron su madurez fisiológica , se cortaron

y

almacenaron de 10 a 15 días a la sombra y a temperatura ambiente con el

madurez fisiológica , se cortaron y almacenaron de 10 a 15 días a la sombra

y

a temperatura ambiente con el objetivo de favorecer y obtener la madurez completa

a 15 días a la sombra y a temperatura ambiente con el objetivo de favorecer

y

obtener la madurez completa de las semillas , luego se procedió a la extracción

. Las semillas fueron secadas a la sombra a temperatura ambiente durante 10 días

y

posteriormente conservadas en sobres de papel estraza hasta su siembra . Establecimiento

conservadas en sobres de papel estraza hasta su siembra . Establecimiento en campo

y

manejo del cultivo , de híbridos y progenitores La semilla híbrida resultante

hasta su siembra . Establecimiento en campo y manejo del cultivo , de híbridos

y

progenitores La semilla híbrida resultante y las de sus progenitores se sembraron

manejo del cultivo , de híbridos y progenitores La semilla híbrida resultante

y

las de sus progenitores se sembraron en agosto de 2015 de la misma forma que

forma que en el primer ciclo . Cuando las plántulas tuvieron 12 cm de altura

y

50 días después de la siembra , fueron establecidas en camas elevadas de 25 cm

y

de 1 , 60 m de ancho , con acolchado plástico de color negro y

y de 1 , 60 m de ancho , con acolchado plástico de color negro

y

riego por goteo . Se establecieron a doble hilera en forma de tresbolillo con 40

azar con 3 repeticiones . Cada parcela experimental fue constituida por 10 plantas

y

como parcela útil 3 plantas con competencia completa ; se tomaron las plantas del

realizaron riegos cada tercer o cuarto día , dependiendo de las necesidades del cultivo ,

y

se realizaron aplicaciones de fertilizantes foliares para corregir deficiencias por

. Además se realizaron podas , tutoreos , fertilizaciones y para la prevención

y

control de mosca blanca , paratrioza ( pulgón saltador ) y trips , se usaron

y para la prevención y control de mosca blanca , paratrioza ( pulgón saltador )

y

trips , se usaron imidaclopid 30 , 7 % y metomilo 54 % . Para

paratrioza ( pulgón saltador ) y trips , se usaron imidaclopid 30 , 7 %

y

metomilo 54 % . Para prevenir las enfermedades se aplicó mancozeb + metalaxil +cloratolonil

estos productos se aplicaron mensualmente , para evitar la resistencia de plagas

y

enfermedades a los productos químicos . Mediciones de componentes de rendimiento

una balanza electrónica ( VELAB , Scientific , México ) con capacidad de 1 kg

y

una precisión de 0 , 001kg , se sumó el peso de fruto obtenido a

mediante el conteo de todos los frutos cosechados a lo largo del ciclo de producción

y

dividido entre el número de plantas cosechadas , mientras que el peso promedio de

de calidad del fruto fueron determinadas en el Laboratorio de Nutrición Vegetal

y

Cultivo de Tejidos del Departamento de Horticultura de la Universidad Autónoma Agraria

la absorbancia de las muestras analizadas , las cuales se leyeron por triplicado

y

el contenido se obtuvo mediante la siguiente fórmula : La cuantificación del contenido

de chile en la etapa de madurez fisiológica por el método descrito por Bennett

y

Kirby ( 1968 ) , en un espectrofotómetro Bio - 145025 BIO - MATE -

40 mg ml - 1 . Las lecturas se realizaron por triplicado para cada muestra

y

el contenido de capsaicina se expresó en unidades Scoville ( SHU ) . Estimación

Estimación de variables agronómicas Para la estimación de días a floración ( DAF )

y

días a cosecha ( DAC ) se contabilizaron los días transcurridos desde que se trasplantaron

genotipos hasta el día en que se presentó la primera florescencia de cada planta

y

hasta el día en que se cosecharon los primeros frutos , respectivamente . Mientras

, con una prueba de comparación de medias de Tukey ( P≤0 , 05 )

y

la estimación de heterosis se consideró como la relación entre el promedio del

promedio de sus progenitores co - rrespondientes , expresada en términos porcentuales

y

la heterobeltiosis ( o trasgresividad híbrida ) se estimó como la superioridad

superioridad del híbrido sobre el mejor progenitor también en escala porcentual ( Fonseca

y

Patterson , 1968 ) , como se indica a continuación : Resultados y discusión

( Fonseca y Patterson , 1968 ) , como se indica a continuación : Resultados

y

discusión Rendimiento y componentes del rendimiento Los análisis de varianza

, 1968 ) , como se indica a continuación : Resultados y discusión Rendimiento

y

componentes del rendimiento Los análisis de varianza aplicados a progenitores e

progenitores e híbridos exhibieron diferencias significativas ( P<0 , 01 ) entre híbridos

y

progenitores en RTF , NFP y PPF , lo que sugiere la existencia de variabilidad

diferencias significativas ( P<0 , 01 ) entre híbridos y progenitores en RTF , NFP

y

PPF , lo que sugiere la existencia de variabilidad genética de utilidad para desarrollar

híbrido P1xP7 con un rendimiento de 1980 g pl - 1 . Los híbridos P1xP5

y

P2xP6 fueron significativamente superiores a todos los progenitores es - tudiados

resultando genotipos prometedores para el desarrollo de variedades de chile Jalapeño

y

Mirador para su producción en invernadero , ya que el híbrido P1xP5 supero en 333

ya que el híbrido P1xP5 supero en 333 % la media de los progenitores hembras

y

en 159 % a los progenitores macho . Los resultados obtenidos permiten inferir que

significativa , ya que hubo la expresión de vigor híbrido en las progenies obtenidas

y

estudiadas . En este sentido Martínez - Sánchez et al . ( 2010 ) menciona

rendimiento es debido a que las plantas fueron estudiadas bajo condiciones de invernadero

y

estas condiciones favorecieron la expresión de los genes relacionados con el rendimiento

rendimiento de fruto . Duarte et al . ( 2012 ) logró la mayor producción

y

calidad de Jalapeño variedad Mitla , al usar gallinaza + 80N en invernade ro ,

significativamente superior a todos los progenitores machos , al progenitor hembra Mirador

y

al híbrido P1xP7 , sin embargo no fue significativamente diferente del resto de

superó en 104 % a la media de producción de fruto de los progenitores hembra

y

en 725 % a los progenitores macho , por lo tanto este híbrido se confirma

en invernadero . Tanto el híbrido P1xP5 como el P2xP6 coinciden con el mayor NFP

y

ADP , lo cual indicaría que el mayor desarrollo de la planta indujo mayor desarrollo

que el mayor desarrollo de la planta indujo mayor desarrollo de yemas reproductivas

y

por lo tanto mayor RTF . Para la variable PPF se encontró que los progenitores macho

el progenitor macho UANCn , superando significativamente a los progenitores hembra

y

a los híbridos , estos últimos tuvieron un peso promedio similar a los progenitores

a los progenitores hembra . Los PPF de progenitores hembra oscilan entre 9 , 2

y

19 , 6 g mientras que los progenitores macho fluctúan entre 86 , 4 y

y 19 , 6 g mientras que los progenitores macho fluctúan entre 86 , 4

y

105 , 6 g y los híbridos tienen pesos que van desde 25 , 2

mientras que los progenitores macho fluctúan entre 86 , 4 y 105 , 6 g

y

los híbridos tienen pesos que van desde 25 , 2 a 38 , 1 g

, como es el caso del contenido de ácido as - córbico , carotenoides totales

y

contenido de capsaicina . El análisis de varianza realizado a estas variables mostró

exhibieron una media de 186 , 8 mg / 0 , 1 Kg de fruto

y

solamente los híbridos P3xP4 y P3xP8 presentaron cantidades significativamente iguales

, 8 mg / 0 , 1 Kg de fruto y solamente los híbridos P3xP4

y

P3xP8 presentaron cantidades significativamente iguales a los progenitores macho

significativamente iguales a los progenitores macho ( Tabla 3 ) . Los progenitores hembras Mitla

y

Mirador presentaron cantidades de CAA significativamente inferiores a los progenitores

, 1 Kg de fruto , siendo el progenitor UANRd el que presentó mayor valor

y

fue significativamente diferente de los progenitores maternos que oscilaron entre

significativamente diferente de los progenitores maternos que oscilaron entre 55 , 8

y

101 , 2 mg / 0 , 1 Kg , mientras que los híbridos presentaron

mayor CAA . En relación al contenido de CT , cuatro de los progenitores macho

y

el progenitor hembra Mitla presentaron las mayores cantidades , con un rango de 6430

progenitores dentro de este rango fueron significativamente superiores al progenitor Mirador

y

a todos los híbridos estudiados . Es importante señalar que estos metabolitos actúan

encontraron seis que fueron significativamente iguales a los progenitores hembra

y

tres que fueron significativamente iguales a los progenitores macho . El contenido

hembra superaron en un 93 % en el contenido de CAPs a los progenitores macho

y

en 32 % a los híbridos . Lo observado anteriormente indica que en relación a

significativamente los costos a inicio de cosecha en la producción en invernadero

y

proporcionando cosechas más tempranas que pueden competir mejor en el mercado .

respecto a esta característica contrastando con lo observado por Montaño - Mata

y

Belisario - Ramos ( 2012 ) , quienes encontraron que el genotipo más precoz

tienen como principales ventajas adelantar cosechas así como evadir algunas plagas

y

/ o enfermedades , ya que algunas plagas y enfermedades son más severas en ciertas

cosechas así como evadir algunas plagas y / o enfermedades , ya que algunas plagas

y

enfermedades son más severas en ciertas épocas del año , logrando con ello menos

en ciertas épocas del año , logrando con ello menos aplicación de insecticidas

y

mayores cosechas de mayor calidad . En los DAC también se encontraron diferencias

el más tardío , aunque fue significativamente igual a cuatro progenitores más

y

tres híbridos , en un segundo grupo quedaron siete progenitores y ocho híbridos

progenitores más y tres híbridos , en un segundo grupo quedaron siete progenitores

y

ocho híbridos , que fueron significativamente iguales , mientras que en un tercer

el híbrido P2xP5 , aunque este fue significativamente igual a cinco progenitores

y

ocho híbridos . Los progenitores presentaron en promedio 133 , 9 días a cosecha

los 88 , 5 días después del trasplante se inició la cosecha en este híbrido

y

en sus progenitores , hembra y macho , a los 139 , 0 y 173

trasplante se inició la cosecha en este híbrido y en sus progenitores , hembra

y

macho , a los 139 , 0 y 173 , 7 días , respectivamente .

híbrido y en sus progenitores , hembra y macho , a los 139 , 0

y

173 , 7 días , respectivamente . En los progenitores se observó una diferencia

los progenitores se observó una diferencia de 57 , 0 días entre el más precoz

y

el más tardío , mientras que en los híbridos fue de 44 , 0 días

genotipos estudiados ( Tabla 4 ) , difiriendo con lo descrito por Montaño - Mata

y

Belisario - Ramos ( 2012 ) quines iniciaron la cosecha a los 75 , 0

que las temperaturas se mantuvieron bajas , retardando el crecimiento del cultivo ,

y

a que en este periodo la radiación solar también fue baja , en algunos días

también fue baja , en algunos días nula , lo que retrasó el crecimiento vegetativo

y

reproductivo de las plantas . En este sentido , se indica que temperaturas inferiores

mientras que las temperaturas óptimas para este cultivo en el día son 20 a 25°C

y

las nocturnas de 16 a 18°C ( Pilatti et al . , 2005 ) .

la primera cosecha del pimiento más precoz a los 91 , 0 días después del trasplante

y

el más tardío a los 117 , 0 días , valores inferiores a los encontrados

diferencias son consecuencia de la variabilidad presente entre las poblaciones de machos

y

hembras bajo estudio , además es posible inferir que en las poblaciones bajo estudio

, ( 2005 ) ; Wong - Romero et al . , ( 2006 )

y

Guerrero - Guerrero et al . , ( 2011 ) que identificaron híbridos en maíz

et al . , ( 2011 ) que identificaron híbridos en maíz con mayor altura

y

alto rendimiento , con alta variabilidad genética presente en los progenitores estudiados

condiciones de cultivo en invernadero existe un mayor control en los factores ambientales

y

manejo del cultivo que influye en el crecimiento y desarrollo . La comparación de

control en los factores ambientales y manejo del cultivo que influye en el crecimiento

y

desarrollo . La comparación de medias ( Tabla 4 ) nos muestra que todos los

todos los híbridos formados son superiores en altura a sus progenitores ( macho

y

hembra ) lo cual resultó conveniente puesto que se reflejó en el rendimiento de

se reflejó en el rendimiento de la planta , ya que había mayor área foliar

y

por lo tanto mayor acumulación de reservas que permitieron una alta productividad

% , lo que indica que los progenitores utilizados presentan divergencia genética

y

sobre todo complementariedad génica , que permite la expresión del vigor hibrido

NFP también se encontró un rango de heterosis que fue de 66 a 243 %

y

nuevamente fue el híbrido P2xP6 el que presentó el mayor valor ; sin embargo en

los frutos fueron más semejantes a los progenitores hembra . En la variable DAF

y

DAC también se encontraron valores elevados y negativos , lo cual es favorable ,

progenitores hembra . En la variable DAF y DAC también se encontraron valores elevados

y

negativos , lo cual es favorable , ya que en los híbridos formados se ganó

ADP el híbrido P2xP6 manifestó una heterosis de 128 % . En las variables CAA

y

CT se encontraron valores de heterosis negativa , lo cual indica que los valores

pimiento en cuanto a tamaño , sin embargo tiene alto RTF , NFP , ADP ,

y

CAPs , lo cual lo hace prometedor para el desarrollo de una variedad de chile

de la heterosis en un cruce sea de cuando menos 20 % ( Vasal

y

Córdova , 1996 ) . Sin embargo en el presente trabajo , la heterosis estimada en

, 1996 ) . Sin embargo en el presente trabajo , la heterosis estimada en RTF

y

NFP , está por encima del nivel mínimo recomendado , lo que indica que existe

puesto que la heterosis exhibida en sus cruzas simples depende de la aptitud rendidora

y

de la diversidad genética de los progenitores usados ( Moll et al . , 1962

, 1962 ) . Con la hibridación se busca mayor precocidad en la progenie

y

en este trabajo , los días a floración y días a cosecha de las primeras ramificaciones

busca mayor precocidad en la progenie y en este trabajo , los días a floración

y

días a cosecha de las primeras ramificaciones fueron indicadores de precocidad en

indicadores de precocidad en las poblaciones evaluadas . Los valores negativos en DAF

y

DAC son ventajosas ya que permiten tener cosechas más tempranas , datos similares

Heterobeltiosis La tabla 6 muestra que para RTF los híbridos P2xP6 , P1xP5

y

P1xP6 presentaron los valores más altos de heterobeltiosis , con valores de 178

presentaron los valores más altos de heterobeltiosis , con valores de 178 , 146

y

115 % , respectivamente , y estos mismos híbridos presentaron valores de 104 , 74

de heterobeltiosis , con valores de 178 , 146 y 115 % , respectivamente ,

y

estos mismos híbridos presentaron valores de 104 , 74 y 82 % de heterobeltiosis

115 % , respectivamente , y estos mismos híbridos presentaron valores de 104 , 74

y

82 % de heterobeltiosis para la variable NFP . Considerando estas dos variables ,

el desarrollo de variedades para producción en invernadero . En la variable DAF

y

DAC los valores más negativos indican que estos híbridos presentan las mayores

de cosechas más tempranas . Además , también se encontró heterobeltiosis alta

y

positiva en la variable ADP , lo cual permite tener un mejor aprovechamiento del

con la calidad de fruto no fue posible mejorar los contenidos de CAA , CT

y

CAPs , ya que todos los híbridos exhibieron valores negativos de heterobeltiosis

invernaderos de baja tecnología , al cruzar chiles picantes ( Jalapeño , Mirador

y

Serrano ) con pimientos , ya que en este trabajo se obtuvieron rendimientos que superaron

combinación adecuada de genes , permitió obtener híbridos con valores positivos

y

altos de heterosis y heterobeltiosis en rendimiento total de fruto y número de frutos

de genes , permitió obtener híbridos con valores positivos y altos de heterosis

y

heterobeltiosis en rendimiento total de fruto y número de frutos por planta . En

valores positivos y altos de heterosis y heterobeltiosis en rendimiento total de fruto

y

número de frutos por planta . En los híbridos desarrollados fue posible lograr

En los híbridos desarrollados fue posible lograr valores negativos , en heterosis

y

heterobeltiosis , en días a floración y días a cosecha , logrando tener genotipos

lograr valores negativos , en heterosis y heterobeltiosis , en días a floración

y

días a cosecha , logrando tener genotipos de cosecha temprana para producción en

presente investigación permiten inferir que es posible mejorar los rendimientos

y

calidad nutricional de frutos , con características de los progenitores hembras

preliminares obtenidos . Dado que se utilizaron líneas puras para formar los híbridos

y

estos fueron manejados en invernadero con control parcial de variables climáticas

cortos de comercialización territoriales : conceptualización , caracterización

y

aplicación para el caso de Córdoba , España Introducción Las dos últimas

seguridad alimentaria a nivel mundial . Primero , se ha rescatado el papel social

y

ambiental de la agricultura que ha despertado el interés por el desarrollo rural

agricultura que ha despertado el interés por el desarrollo rural territorial a nivel local

y

regional ( Marsden , 2017 ) . Segun do , ha surgido un consumo más

surgido un consumo más consciente en algunos territorios , en el que consumidores

y

productores , como ciudadanos activos , dan nuevas formas y valores a sus relaciones

en el que consumidores y productores , como ciudadanos activos , dan nuevas formas

y

valores a sus relaciones en el sistema agroalimentario , tomando gran relevancia

sus relaciones en el sistema agroalimentario , tomando gran relevancia el tiempo

y

el espacio en el que tienen lugar ( Goodman et al . , 2012 ;

2017 ) . Y , finalmente , han emergido cuestiones relacionadas con la sostenibilidad

y

sus conexiones con la “ nueva ecuación alimentaria ” caracterizada por la subida

degradación de los recursos naturales , los efectos del cambio climático ( Morgan

y

Sonnino , 2010 ) y la concentración de poder en las industrias agroalimentarias

recursos naturales , los efectos del cambio climático ( Morgan y Sonnino , 2010 )

y

la concentración de poder en las industrias agroalimentarias ( IPES - Food , 2017

. , 2017 ) . La reconfiguración de los vínculos entre las áreas rurales

y

urbanas , a través de movimientos alimentarios alternativos , gestados en oposición

alimentario sostenible , que satisface al mismo tiempo las demandas de ( in ) seguridad

y

sostenibilidad alimentaria ( Marsden y Morley , 2014 ) . Desde sus inicios ,

mismo tiempo las demandas de ( in ) seguridad y sostenibilidad alimentaria ( Marsden

y

Morley , 2014 ) . Desde sus inicios , estos nuevos movimientos fueron conceptualizados

la configuración ( o reconfiguración ) de un tipo de vínculos entre productores

y

consumidores diferentes a los establecidos bajo el régimen alimentario globalizado

. , 2016 ) . El constructo de RAA ha sido abordado desde diferentes enfoques

y

aproximaciones teóricas , siendo los Canales Cortos de Comercialización ( CCC )

más difundidas en la literatura científica . El concepto de CCC es polisémico

y

, desde su introducción en la literatura de los estudios rurales y agroalimentarios

es polisémico y , desde su introducción en la literatura de los estudios rurales

y

agroalimentarios , ha estado en constante cambio y evolución . Algunos académicos

literatura de los estudios rurales y agroalimentarios , ha estado en constante cambio

y

evolución . Algunos académicos diferencian dos perspectivas teóricas para abordarlo

perspectivas teóricas para abordarlo : la primera enfocada en el rol de los productores

y

la contribución potencial de este tipo de redes al desarrollo rural sostenible (

Goodman et al . , 2012 ; Renting et al . , 2012 ; Mundler

y

Laughrea , 2016 ) , y la segunda con un enfoque en el papel activo del

; Renting et al . , 2012 ; Mundler y Laughrea , 2016 ) ,

y

la segunda con un enfoque en el papel activo del consumidor , que a través del

enfoque en el papel activo del consumidor , que a través del consumo consciente

y

movimientos sociales intenta transformar las estructuras del sistema alimentario

perspectiva territorial , en el que sean considerados todos los elementos que componen

y

definen el territorio ; esto es , los recursos territoriales , los actores y los

componen y definen el territorio ; esto es , los recursos territoriales , los actores

y

los acuerdos institucionales ( Sánchez - Zamora , 2014 ) . En este contexto

presente artículo tiene dos objetivos : primero , avanzar en la conceptualización

y

caracterización de los Canales Cortos de Comercialización Territoriales ( CCCT

profundizando en sus orígenes , carac terísticas e implicaciones con el territorio ;

y

, segundo , analizar el caso concreto de los CCCT en la ciudad de Córdoba

Andalucía se ha producido un aumento en el número de establecimientos mayoristas

y

minoristas especializados en la venta de productos ecológicos ( MAPAMA , 2015 )

2012 ) asistiendo a una proliferación de diferentes formas de autoabastecimiento

y

CCCT . En 2016 , Córdoba se adhiere al Pacto de Política Alimentaria Urbana de

En 2016 , Córdoba se adhiere al Pacto de Política Alimentaria Urbana de Milán

y

, en 2017 , a la Red de Ciudades por la Agroecología . En ambas

ambas iniciativas se promueven sistemas agroalimentarios sostenibles , inclusivos

y

resilientes , donde los CCC se configuran como una línea prioritaria . Material

y

métodos Dada la naturaleza de este artículo , se ha establecido una metodología

objetivos propuestos . El proceso metodológico fue el siguiente : Conceptualización

y

caracterización de los CCCT Para conceptualizar y caracterizar los CCCT se ha procedido

siguiente : Conceptualización y caracterización de los CCCT Para conceptualizar

y

caracterizar los CCCT se ha procedido a revisar la literatura existente y a su posterior

conceptualizar y caracterizar los CCCT se ha procedido a revisar la literatura existente

y

a su posterior análisis crítico . Así el proceso metodológico se basó en tres

Selección de artículos Se hizo uso de la herramienta web de código abierto

y

gratuito Tree of science – ToS ( http : / / tos . manizales .

los artículos de acuerdo a tres indicadores : grado de entrada , intermediación

y

grado de salida . Los artículos con grado de entrada alto y salida cero ,

entrada , intermediación y grado de salida . Los artículos con grado de entrada alto

y

salida cero , se les denomina raíces ( artículos clásicos en la temática ) ;

artículos con un grado de intermediación alto tronco ( artículos que han desarrollado

y

/ o consolidado la temática ) ; y por último , y por cuestiones de visualización

alto tronco ( artículos que han desarrollado y / o consolidado la temática ) ;

y

por último , y por cuestiones de visualización , los artículos con un grado de

artículos que han desarrollado y / o consolidado la temática ) ; y por último ,

y

por cuestiones de visualización , los artículos con un grado de salida alto y un

, y por cuestiones de visualización , los artículos con un grado de salida alto

y

un grado de entrada cero , hojas ( artículos con perspectivas específicas o con

crítico A partir del análisis arrojado por ToS , se procede a la revisión

y

análisis crítico , a partir de los cuales se construye el concepto de CCCT .

parte de la información utilizada es primaria , recopilada entre noviembre de 2015

y

diciembre de 2016 en el marco geográfico de la ciudad de Córdoba . El proceso

CCCT a partir de las organizaciones vinculadas a una asociación de productores / as

y

tomando como base el inventario de CCC en Andalucía realizado por González et al

de la referenciación de las iniciativas de la asociación de productores / as ,

y

de otras organizaciones que acompañan los procesos de RAA y soberanía alimentaria

de productores / as , y de otras organizaciones que acompañan los procesos de RAA

y

soberanía alimentaria en la ciudad . Se obtuvo una base de iniciativas que implementan

ciudad . Se obtuvo una base de iniciativas que implementan como estrategia los CCCT

y

se seleccionaron las que comercializan los siguientes productos : frutas , verduras

siguientes productos : frutas , verduras , pan artesanal ecológico , vino ecológico

y

cerveza artesanal . Selección de los potenciales entrevistados Se recabó información

representantes legales o delegados de las iniciativas que participan en CCCT , entre julio

y

octubre de 2016 . Se utilizó un cuestionario semi - estructurado compuesto por dos

estructurado compuesto por dos partes : en la primera se describe la iniciativa ;

y

en la segunda se identifican las interacciones de las organizaciones . Se identificaron

actividades relacionas con temas de los CCCT ( no existe intercambio monetario ) ;

y

en el ámbito comercial , quién vende y quién compra productos en los CCCT (

( no existe intercambio monetario ) ; y en el ámbito comercial , quién vende

y

quién compra productos en los CCCT ( intercambio monetario ) . En cada interacción

intercambio monetario ) . En cada interacción se calificó la intensidad , entre 1

y

10 ( 1 = muy baja intensidad , 10 = muy alta intensidad de

deliberación , entre otros . La información se recogió en una libreta de campo ,

y

fue útil para complementar la información recopilada en las entrevistas . Análisis

de la información se siguieron técnicas de Análisis de Redes Sociales ( ARS )

y

se utilizó el software de código abierto y gratuito Gephi 0 . 9 . 1

Análisis de Redes Sociales ( ARS ) y se utilizó el software de código abierto

y

gratuito Gephi 0 . 9 . 1 ( https : / / gephi . org

conjunto de técnicas de análisis para el estudio de las relaciones entre actores

y

las estructuras sociales que surgen de éstas . Las redes están conformadas por

dos elementos básicos : los actores ( nodos ) representados por puntos en la red

y

las relaciones que se establecen ( aristas ) representados por líneas . Para su

análisis se hace uso del lenguaje matemático de la teoría de grafos , de matrices

y

de álgebra relacional ( Scott , 2017 ) . El ARS ha desarrollado una

( Scott , 2017 ) . El ARS ha desarrollado una serie de indicadores

y

estadísticos que permiten caracterizar las estructuras de las redes y la posición

indicadores y estadísticos que permiten caracterizar las estructuras de las redes

y

la posición de los nodos . La Tabla 1 muestra algunos indicadores importantes .

las que las relaciones no son tan fuertes , pero existen nodos conectores ( Buciega

y

Esparcia , 2013 ) . Resultados Conceptualización y caracterización de los Canales

nodos conectores ( Buciega y Esparcia , 2013 ) . Resultados Conceptualización

y

caracterización de los Canales Cortos de Comercialización Territoriales Concepto

propone analizar tres elementos clave : i ) origen ; ii ) perspectiva territorial ,

y

iii ) definición de los CCCT . i ) Origen Las RAA ,

los CCCT . i ) Origen Las RAA , originadas como reacciones individuales

y

colectivas a la crisis de civilización , y como propuestas alternativas a las contradicciones

originadas como reacciones individuales y colectivas a la crisis de civilización ,

y

como propuestas alternativas a las contradicciones de las cadenas alimentarias globalizadas

como los jardines comunitarios o los jardines productivos de carácter individual .

Y

la segunda caracterizada por el intercambio comercial de los alimentos , en las que

alimentos , en las que tienen lugar relaciones directas entre productores , consumidores

y

otros agentes del sistema alimentario . En esta última se enmarcan los CCCT con

distribución , como los mercados de productores , los sistemas de pedidos por internet ,

y

los esquemas de cesta , entre otros . ii ) Perspectiva territorial Las

entre otros . ii ) Perspectiva territorial Las nuevas relaciones entre productores

y

con - sumidores , gestadas en las RAA , se encuentran estrechamente ligadas al territorio

. , 2012 ) . Aquí , la noción de territorio sobrepasa la idea reduccionista

y

estática de un mero espacio físico dado , proveedor de recursos naturales y materiales

reduccionista y estática de un mero espacio físico dado , proveedor de recursos naturales

y

materiales , para ser concebida como un proceso de construcción social determinado

concebida como un proceso de construcción social determinado por los actores locales

y

sus interacciones para valorizar los recursos territoriales ( Sánchez - Zamora et

diferentes geografías resalta la singularidad de las narrativas que emergen en cada lugar

y

que dan forma a las estrategias alimentarias ( Sonnino , 2016 ) , limitadas o

propios del lugar donde se desarrollan . Cobra aquí importancia el enfoque dinámico

y

evolutivo del territorio , tanto por el espacio físico ( cambios en la disponibilidad

como por los modos de apropiación que hacen del espacio los actores ( Moranta Vidal

y

Pol Urrútia , 2005 ) que constantemente interactúan allí . Estas maneras de interacción

constantemente interactúan allí . Estas maneras de interacción social , convergentes

y

legitimadas , inciden en los procesos de organización territorial ( D ’ Aquino

inciden en los procesos de organización territorial ( D ’ Aquino , 2002 ) ,

y

en la configuración del capital social relacional . iii ) Hacia una definición

Así , al establecer el territorio como marco de desarrollo de las RAA ,

y

en particular para las formas que conllevan un intercambio comercial , surgen los

territorio alimentario específico , con objetivos relacionados con valores sostenibles ,

y

con una clara perspectiva de gobernanza agroalimentaria ( Fabbrizzi et al . , 2014

de distribución : uno a través de la venta directa ( sin intermediario ) ,

y

otro , con un único intermediario entre productor y consumidor . En la conceptualización

directa ( sin intermediario ) , y otro , con un único intermediario entre productor

y

consumidor . En la conceptualización aquí propuesta no se contemplan las formas

origen . Esto resalta la importancia de la proximidad territorial entre la producción

y

el consumo . Los CCCT contribuyen a la construcción social de mercados , entendiéndose

totalidad el proceso depende de la interacción humana ( Von Mises , 1998 ) ,

y

se constituye , por tanto , en un espacio de relaciones sociales , en el que

relaciones sociales , en el que se enfatiza el carácter dinámico de la interacción

y

de la comunicación , así como la convergencia de intereses más allá del mero

económico . Esta construcción social crea en cada territorio un mercado propio

y

particular , es decir , un mercado territorial o anidado ( van der Ploeg et

anglosajona ) , al subrayar su componente social como base para la acción económica

y

enfatizar la idea que las acciones económicas están enclavadas o embebidas en sistemas

en sistemas concretos de relaciones sociales ( Granovetter , 1985 ) . Finalmente ,

y

como consecuencia de todo lo expresado anteriormente , se ha elaborado la Figura

anteriormente , se ha elaborado la Figura 1 , en la que se recogen sintética

y

visualmente los elementos de análisis que , en nuestra opinión deben considerarse

considerarse para el caso de los CCCT . La Figura 1 , resalta la relación iterativa

y

bidireccional que tienen lugar entre las RAA y el territorio , al incidir los resultados

1 , resalta la relación iterativa y bidireccional que tienen lugar entre las RAA

y

el territorio , al incidir los resultados de cada elemento en la configuración de

otro elemento . Justamente , en estas interacciones que tienen lugar entre las RAA

y

el territorio es donde surgen los CCCT . Características de los CCCT Una

existencia de una base social robusta ; ii ) la participación de pequeños productores

y

minoristas ; iii ) la construcción confianza - calidad ; iv ) y la conexión

pequeños productores y minoristas ; iii ) la construcción confianza - calidad ; iv )

y

la conexión con economías transformadoras . Estos elementos serán detallados en

reúne . Así , esta construcción de narrativas conjuntas , genera cambios en el rol

y

el peso de los diferentes mecanismos de gobernanza en los sistemas agroalimentarios

peso de los diferentes mecanismos de gobernanza en los sistemas agroalimentarios ,

y

principalmente sugiere la necesidad de una reevaluación del rol de la sociedad civil

seguridad alimentaria ( Candel , 2014 ) , así como en la sostenibilidad ( Marsden

y

Morley , 2014 ) , al recombinar democracia , conocimiento y acción colectiva ( Sonnino

sostenibilidad ( Marsden y Morley , 2014 ) , al recombinar democracia , conocimiento

y

acción colectiva ( Sonnino et al . , 2016 ) . ii ) Participación

et al . , 2016 ) . ii ) Participación de pequeños productores

y

minoristas Vistos desde la producción , los CCCT pro - mueven la vinculación de

pequeños productores , tipificados en algunos países , como agricultura familiar

y

/ o campesina , en las que se resaltan las características de la multifuncionalidad

características de la multifuncionalidad de la agricultura ( o el territorio rural ) ( Marsden

y

Sonnino , 2008 ) , prevaleciendo la producción de policultivos y la pluriactividad

rural ) ( Marsden y Sonnino , 2008 ) , prevaleciendo la producción de policultivos

y

la pluriactividad ( van der Ploeg et al . , 2012 ) , así como

un movimiento agroecológico , en el que los alimentos llevan esta calificación ,

y

2 ) un movimiento campesino , en el que los alimentos son reconocidos como campesinos

campesinos , artesanales o naturales . Estas dos tendencias no compiten entre sí

y

en algunos territorios resulta imposible hacer la distinción entre ambas . Sin embargo

, como los europeos ( van der Ploeg et al . , 2012 ) ,

y

del movimiento campesino en países en vías de desarrollo , como los latinoamericanos

campesino en países en vías de desarrollo , como los latinoamericanos ( Ranaboldo

y

Arosio , 2014 ) . Van der Ploeg et al . ( 2012 ) señalan

siendo una de las principales caraterísticas la oferta de productos de alta calidad

y

de producción ecológica . Asimismo , indica que los consumidores involucrados pertenecen

proliferado los CCCT , se resalta la vinculación de economías familiares campesinas ,

y

la artículación de consumidores preocupados tanto por su propia seguridad alimentaria

proporción de los consumidores manifestó su interés por apoyar a las familias campesinas

y

la satisfacción por recordar sus vínculos con la ruralidad ( Parrado y Molina ,

campesinas y la satisfacción por recordar sus vínculos con la ruralidad ( Parrado

y

Molina , 2014 ) . Vistos desde la distribución , las formas de CCCT

intermediario , vinculan a pequeños minoristas locales , con ventas especializadas

y

diferenciadas de productos ecológicos , campesinos , naturales o artesanales . Algunos

artesanales . Algunos de estos establecimientos son emprendimientos individuales ,

y

otros de carácter colectivo como cooperativas de consumidores . iii ) Construcción

terceras partes , como mecanismo de entrada a ciertos canales de comercialización

y

como respaldo de credibilidad de sus productos . Sin embargo , en los CCCT no todos

se promueven nuevas formas de coordinación social como respaldo a sus prácticas

y

alimentos , tales como los Sistemas Participativos de Garantía o sellos de confianza

” del campesino para que el comprador pueda conocer a su proveedor , el lugar

y

el proceso de producción ( Parrado y Molina , 2014 ) . A través

pueda conocer a su proveedor , el lugar y el proceso de producción ( Parrado

y

Molina , 2014 ) . A través de estos mecanismos coordinados de calidad ,

) . A través de estos mecanismos coordinados de calidad , se busca rescatar

y

conservar la gastronomía propia de cada territorio , que abarca tanto el conjunto

como las costumbres alimentarias , las tradiciones , los procesos , las personas

y

los estilos de vida que se definen alrededor de la misma ( Di Clemente et

razas de animales , estacionalidad , etc . ) , procesos de elaboración ( recetas

y

procedimientos ) , utensilios de preparación y rituales sociales de consumo ( algunas

procesos de elaboración ( recetas y procedimientos ) , utensilios de preparación

y

rituales sociales de consumo ( algunas recetas únicamente se preparan en determinadas

) . En los CCCT se promueve la conservación de estas tradiciones , al recuperar

y

conservar semillas nativas , promover productos de elaboración artesanal y propiciar

recuperar y conservar semillas nativas , promover productos de elaboración artesanal

y

propiciar encuentros entre productores y consumidores , en los que comparten recomendaciones

promover productos de elaboración artesanal y propiciar encuentros entre productores

y

consumidores , en los que comparten recomendaciones de conservación y preparación

productores y consumidores , en los que comparten recomendaciones de conservación

y

preparación de los alimentos . iv ) Conexión con economías transformadoras Muchas

ecosistema de economías transformadoras , entre las que se encuentran : la social

y

solidaria , las colaborativas , las comunitarias , las feministas , la circular ,

y

el emprendimiento social , entre otras . Pero también vinculadas a economías no

diferente medida , buscan desplazar del centro del sistema económico el capital ,

y

colocar en su lugar el sostenimiento de la vida ( algunas con enfoque antropológico

antropológico ) , lo que permite transiciones como la economía basada en valores

y

la incorporación de raciocinios diferentes a los económicos en las elecciones de

de consumo . Así , se aboga por un consumo responsable , el comercio justo

y

la sostenibilidad . Estos son conceptos relevantes que involucran en muchos casos

desarrollando variadas formas de CCCT , adaptadas a las características del territorio ,

y

en constante proceso de evolución . En los siguientes apartados se describen cada

diferentes organizaciones vinculadas . Para ello i ) se identifican los actores , ii )

y

sus interacciones en dos ámbitos : comunitario y comercial ( ARS ) ; y iii

identifican los actores , ii ) y sus interacciones en dos ámbitos : comunitario

y

comercial ( ARS ) ; y iii ) se describe el tipo de CR que

) y sus interacciones en dos ámbitos : comunitario y comercial ( ARS ) ;

y

iii ) se describe el tipo de CR que se origina . i ) Actores

de instituciones de la sociedad civil ( p . ej . ONGs , medios de comunicación

y

diferentes formas de organización social ) . En la red comercial , la presencia

de las instituciones públicas orientada a la gestión de políticas desaparece ,

y

la de organismos de gestión del conocimiento disminuye notoriamente , prevaleciendo

entre ellos ( a través de ARS ) , tomando como referencia el ámbito comunitario

y

el comercial . La Tabla 3 recoge algunos indicadores importantes de la red que se

recoge algunos indicadores importantes de la red que se configura en cada ámbito

y

en la Figura 2 y la Figura 3 se Los datos de la Tabla 3 indican

indicadores importantes de la red que se configura en cada ámbito y en la Figura 2

y

la Figura 3 se Los datos de la Tabla 3 indican la participación de un

3 indican la participación de un mayor número de actores en la red comunitaria

y

por tanto de mayores interacciones con relación a la red comercial . La estructura

que se encuentran bases sólidas para la existencia de un importante CR ( Buciega

y

Esparcia , 2013 ) si se tiene en cuenta 1 ) que las redes son

que las redes son completas , es decir que no existe ningún actor desconectado ,

y

2 ) que al analizar el grado de centralidad o la posición de cada actor

concentrado en un actor en particular , como puede comprobarse visualmente en la Figura 2

y

Figura 3 . Para conocer los actores que tienen un rol importante en la red

comunitaria ( Figura 2 ) los actores - puente son los productores de alimentos frescos

y

procesados , los minoristas o mercado de productores y una organización de apoyo

productores de alimentos frescos y procesados , los minoristas o mercado de productores

y

una organización de apoyo ; mientras en la red comercial ( Figura 3 ) lo

mientras en la red comercial ( Figura 3 ) lo son los puntos de comercio

y

los productores de alimentos frescos y procesados . Estos actores al estar mejor

Figura 3 ) lo son los puntos de comercio y los productores de alimentos frescos

y

procesados . Estos actores al estar mejor integrados en la red , pueden tener mayor

integrados en la red , pueden tener mayor influencia en el proceso de planificación

y

gestión , así como intervenir en las interacciones y acceder a mayores recursos

proceso de planificación y gestión , así como intervenir en las interacciones

y

acceder a mayores recursos dentro de la red ( Lienert et al . , 2013

instituciones de carácter público , de gestión de políticas o del conocimiento ;

y

de carácter privado , con diversidad de organizaciones de la sociedad civil . Así

características , existencia de relaciones entre actores con características diferentes

y

existencia de intermediación , indican la configuración de un CR de tipo puente

) . El CR de tipo puente se refiere a relaciones menos estrechas e intensas

y

más diversas , que son importantes para que las redes avancen y evolucionen hacia

estrechas e intensas y más diversas , que son importantes para que las redes avancen

y

evolucionen hacia formas más cohesionadas ( Putman , 2001 ) . Las relaciones tipo

más cohesionadas ( Putman , 2001 ) . Las relaciones tipo puente son estratégicas

y

favorecen el acceso a otro tipo de recursos , conocimiento ( Buciega y Esparcia ,

estratégicas y favorecen el acceso a otro tipo de recursos , conocimiento ( Buciega

y

Esparcia , 2013 ) e innovaciones , que en caso de actuación individual de los

CCCT puede ser entendido como un “ recurso ” , que proporciona a los individuos

y

organizaciones vinculadas acceso a otros tipos de capitales , y allí se podría

a los individuos y organizaciones vinculadas acceso a otros tipos de capitales ,

y

allí se podría suscribir la idea de que detrás de la creación o mantenimiento

actores . No obstante , también puede ser interpretado como un “ instrumento ” ,

y

no tanto como un fin en sí mismo . Así , se resalta el potencial

influencia del marco institucional , que puede reforzarlo o debilitarlo ( Buciega

y

Esparcia , 2013 ) . De esta forma , el CR que se ha descrito

, el CR que se ha descrito no puede ser entendido como un elemento estable

y

permanente en Córdoba , sino como un atributo dinámico y cambiante en el tiempo

como un elemento estable y permanente en Córdoba , sino como un atributo dinámico

y

cambiante en el tiempo , y dependiente de la interacción de los diferentes actores

permanente en Córdoba , sino como un atributo dinámico y cambiante en el tiempo ,

y

dependiente de la interacción de los diferentes actores que se involucren . Participación

los diferentes actores que se involucren . Participación de pequeños productores

y

minoristas Los productores que se encuentran vinculados a los diferentes CCCT en

área de menos de tres ha . Mientas los productos procesados ocupan un área mayor

y

algunos de ellos no intervienen en la obtención de las materias primas . En cuanto a

algunos productores hacen mención específica de prácticas como la agroecología

y

la permacultura . La conservación y recuperación de semillas nativas son algunas

específica de prácticas como la agroecología y la permacultura . La conservación

y

recuperación de semillas nativas son algunas de las actividades que vienen desarrollando

en los CCCT en Córdoba . Así se han rescatado la zanahoria arcoíris ( morada

y

blanca ) , y algunas variedades de ajetes y nabos ( Tabla 4 ) .

Córdoba . Así se han rescatado la zanahoria arcoíris ( morada y blanca ) ,

y

algunas variedades de ajetes y nabos ( Tabla 4 ) . Por otro lado , la red

rescatado la zanahoria arcoíris ( morada y blanca ) , y algunas variedades de ajetes

y

nabos ( Tabla 4 ) . Por otro lado , la red comunitaria ( Figura 2 )

adoptan diferentes formas para llegar al consumidor final , algunas de venta directa

y

otras a través de un intermediario . Los productores se vinculan a los CCCT a

vinculan a los CCCT a título individual o a través de cooperativas de productores

y

consumidores , y , generalmente , participan en diferentes formas de CCCT paralelamente

a título individual o a través de cooperativas de productores y consumidores ,

y

, generalmente , participan en diferentes formas de CCCT paralelamente . La Tabla

la oferta de productos especializados , ya sean ecológicos o de carácter local

y

artesanal . En este grupo se ubican , emprendimientos individuales como las tiendas

grupo se ubican , emprendimientos individuales como las tiendas especializadas ,

y

emprendimientos colectivos como la cooperativa de productores y consumidores y el

especializadas , y emprendimientos colectivos como la cooperativa de productores

y

consumidores y el mercado social . Construcción confianza - calidad Como indican

, y emprendimientos colectivos como la cooperativa de productores y consumidores

y

el mercado social . Construcción confianza - calidad Como indican la Tabla 4

y

Tabla 5 , la gran mayoría de los CCCT comercializan productos respaldados por una

, ya que ésta se constituye en un requisito de entrada para algunos CCCT ,

y

para obtener ciertos beneficios otorgados por las entidades públicas , como pueden

) , en especial en aquelllos CCCT en los que participan productores certificados

y

no certificados ; esto puede aludir al carácter colectivo de este tipo de CCCT y

y no certificados ; esto puede aludir al carácter colectivo de este tipo de CCCT

y

que de cara a los consumidores se muestra como un solo colectivo social ; este

desarrollado un proceso completo de SPG , pero en sus estatutos de constitución

y

operación han establecido lineamientos acerca de las características de los productos

establecido lineamientos acerca de las características de los productos que compran

y

venden ; este es el caso del mercado social . Conexión con economías transformadoras

Los CCCT en Córdoba se han gestado en el marco de economías alternativas ,

y

bajo dichos principios realizan actividades de sensibilización y comercio . La Tabla

alternativas , y bajo dichos principios realizan actividades de sensibilización

y

comercio . La Tabla 6 resume la relación entre diferentes tipos de iniciativas de

economías transformadoras , siendo tres las más mencionadas : la economía social

y

solidaria , el comercio justo y la economía del bien común . Algunos de los

siendo tres las más mencionadas : la economía social y solidaria , el comercio justo

y

la economía del bien común . Algunos de los CCCT forman parte de organizaciones

de los CCCT forman parte de organizaciones supra - territoriales como federaciones

y

redes de carácter nacional , siendo las más representativas : Federación Andaluza

nacional , siendo las más representativas : Federación Andaluza de Consumidores

y

Productores Ecológicos ( FACPE ) , Red de Redes de Economía Alternativa y Solidaria

Consumidores y Productores Ecológicos ( FACPE ) , Red de Redes de Economía Alternativa

y

Solidaria ( que también existe a nivel provincial y de comunidad autónoma ) , Campo

Redes de Economía Alternativa y Solidaria ( que también existe a nivel provincial

y

de comunidad autónoma ) , Campo de Energía de Córdoba ( promotor de la Economía

, Campo de Energía de Córdoba ( promotor de la Economía del Bien Común )

y

Coordinadora Estatal de Comercio Justo ( CECJ ) . Discusión Las iniciativas de

Las iniciativas de CCCT responden a las características propias del territorio

y

al bagaje histórico que lo ha forjado . Así , aunque se identifiquen formatos de

uno de los componentes del territorio : recursos territoriales , tipo de actores

y

acuerdos institucionales . Estos componentes , inciden en la forma de gestión ,

componentes , inciden en la forma de gestión , los productos que se comercializan

y

los mecanismos de gobernanza que impulsan los CCCT . En este apartado , se exponen

CCCT . En este apartado , se exponen algunas comparaciones de los CCCT en Córdoba

y

otros territorios , que refuerzan este planteamiento . El Mercado de Productores

Córdoba es una iniciativa auto - gestionada por los productores , transformadores

y

una cooperativa de comercio justo , que en el momento de recolección de información

cestas en Córdoba , es una actividad auto - gestionada por los mismos productores ,

y

se viene desarrollando de manera similar en países como Reino Unido ( Kirwan , 2004

Kirwan , 2004 ) , Italia ( Fabbrizzi et al . , 2014 )

y

República Checa ( Zagata , 2012 ) . En este último país el ministerio de

intermediario encargado de acopiar la oferta de diferentes productores , armar las cestas

y

coordinar la entrega a los consumidores finales ; allí los productos son de producción

los consumidores finales ; allí los productos son de producción agroecológica

y

campesina , y no cuentan con una certificación de tercera parte ( Reina - Usuga

consumidores finales ; allí los productos son de producción agroecológica y campesina ,

y

no cuentan con una certificación de tercera parte ( Reina - Usuga et al .

importancia que han tenido las organizaciones de la sociedad civil en la promoción

y

el desarrollo de los CCCT , aunque el tipo de organización que los ha impulsado

, aunque el tipo de organización que los ha impulsado es de diferente naturaleza

y

está vinculada a las características propias de cada territorio . Se resalta la

resalta la importancia de este rol , ya que ello ha permitido identificar problemas

y

respuestas locales a cuestiones relacionadas con la seguridad alimentaria , así

desarrollar la capacidad entre diferentes agencias gubernamentales , sectores políticos

y

escalas de gobernanza ( Candel , 2014 ) . Conclusiones Aunque diferentes autores

. Conclusiones Aunque diferentes autores han sugerido que el lugar geográfico

y

el arraigo territorial tiene una dimensión relacional entre las unidades agrícolas

arraigo territorial tiene una dimensión relacional entre las unidades agrícolas

y

los CCC , no había sido abordada hasta ahora en la literatura académica una aproximación

unos recursos específicos , unos actores ( agentes territoriales ) que interactúan

y

crean acuerdos , e instituciones que transforman y desarrollan el mismo . Esto representa

territoriales ) que interactúan y crean acuerdos , e instituciones que transforman

y

desarrollan el mismo . Esto representa un carácter novedoso para esta investigación

poder . No obstante , existen actores - puente que conectan grupos de nodos aislados

y

destaca su papel como cohesionadores de la red . Se detecta una mayor participación

mayor participación de organizaciones en la red comunitaria que en la comercial ,

y

en ambas redes predominan las organizaciones de la sociedad civil . Los actores

ecológica o agroecológica , denominación bajo la cual comercializan sus productos ,

y

participan en varios CCCT simultáneamente ( participación multicanal ) . La construcción

principalmente en la certificación de producción ecológica ( por tercera parte ) ,

y

en algunos colectivos se observan procesos participativos de garantía . Por otra parte

vinculación de los CCCT a economías transformadoras , como la economía social

y

solidaria , economía del bien común y el comercio justo . Respecto a la metodología

transformadoras , como la economía social y solidaria , economía del bien común

y

el comercio justo . Respecto a la metodología , el Análisis de Redes Sociales

poco uso en el análisis de las RAA . Esta metodología presenta utilidad práctica

y

posibilidad de aplicación en otros estudios , lo que representa una aportación

como identificar la posible relación del CR con el tipo de actores involucrados

y

establecer criterios para analizar el CR en los CCCT ( por ejemplo , cuándo un

los cambios en los mecanismos de gobernanza del sistema alimentario territorial ,

y

de las políticas públicas . Características de los suelos y su fertilidad en

alimentario territorial , y de las políticas públicas . Características de los suelos

y

su fertilidad en relación con la incidencia de bitter pit en plantaciones de manzanas

aumentado en los últimos años debido principalmente a la creciente producción de China

y

de otros países del hemisferio Sur ( FAO , 2016 ) . En Europa se

, 2015 ) . La calidad de los frutos es el atributo de mayor interés

y

es un tema recurrente en la mayoría de los programas de mejora genética . En el

manzana ‘ Golden Delicious ’ se comercializa tras un periodo de frigoconservación

y

la minimización de pérdidas de manzana comercializable durante esta etapa es importante

importante para el sector . Las causas se asocian a la deshidratación , a fisiopatías

y

a podredumbres . Entre las fisiopatías causantes de estas mermas destaca el bitter

bitter pit . El bitter pit es un desorden fisiológico que se presenta en manzano

y

que se caracteriza por la aparición en el fruto de pequeñas manchas de color marrón

de mermas ( Ruralcat , 2017 ) . La fisiopatía se produce en el campo

y

, aunque la alteración puede ya presentarse en recolección , se desarrolla por

frio , cuando el fruto ha alcanzado un estado más avanzado de madurez . Mengel

y

Kirkby ( 1982 ) indicaron que el suministro subóptimo de Ca al fruto era la

años se conoce la estrecha correlación negativa entre aparición de bitter pit

y

el contenido de calcio ( Sharples , 1968 ) . Una elevada concentración de Ca

citar elevados ratios N / Ca en fruto ( Faust , 1989 ) . Trocme

y

Gras ( 1979 ) recogen trabajos de distintos autores que relacionan diversos índices

influyente en la aparición del bitter pit , correlacionando poca carga del árbol

y

calibres grandes con una alta incidencia del desorden ( Ben , 1995 ; Retamales y

y calibres grandes con una alta incidencia del desorden ( Ben , 1995 ; Retamales

y

Lepe , 2000 ) . La influencia del manejo del agua de riego tiene una

un 10 % ) limita el crecimiento radicular en manzano ( Ruhland , 1961 )

y

también porque el calcio se absorbe mayoritariamente en los ápices de las raíces

y

en las raíces de crecimiento secundario ; en consecuen - cia , cuando el contenido

limita la absorción de calcio . Desde la segunda vertiente , otros autores ( Ho

y

White , 2005 ) han demostrado que el estrés hídrico , la baja humedad ambiental

y

la salinidad inducen la deficiencia de calcio . La influencia de estas variables

contribuir a la manifestación del bitter pit . El efecto de las altas temperaturas

y

alta incidencia lumínica favorecen la rápida expansión del fruto , lo que diluye

rápida expansión del fruto , lo que diluye la concentración de calcio en el mismo

y

lo hace más sensible a las alteraciones fisiológicas ( De Freitas y Mitcham , 2012

en el mismo y lo hace más sensible a las alteraciones fisiológicas ( De Freitas

y

Mitcham , 2012 ) . Una mala polinización , que conduce a una disminución

) . Ante la diversidad de factores ( edáficos , agronómicos , meteorológicos

y

de manejo ) que influyen en la aparición de bitter pit ( Jemricet al .

calcio en la pulpa , generando procesos de degradación al transcurrir 10 - 15 días

y

obteniendo un efecto similar al del bitter pit ( Retamales y Valdes , 2001 )

10 - 15 días y obteniendo un efecto similar al del bitter pit ( Retamales

y

Valdes , 2001 ) ; la aceleración de la madurez con la aplicación de etileno

determinar las características edáficas que inciden en la aparición del bitter pit

y

que pueden estudiarse antes de la plantación , con la finalidad de que puedan ser

puedan ser utilizados en la decisión de conservar dichos frutos o no . Material

y

métodos Localización , climatología y parcelas estudiadas El área de estudio

conservar dichos frutos o no . Material y métodos Localización , climatología

y

parcelas estudiadas El área de estudio se centró en la denominada Plana de Lleida

zona son 383 mm ( Poal ) , 377 mm ( Vallfogona de Balaguer )

y

359 mm ( Raimat ) . La evapotranspiración para el mismo periodo , calculada según

( Allen et al . , 1998 ) , alcanza los 961 mm en Poal

y

953 mm en Raimat . La temperatura media de la zona es de 14 ,

la serie 1989 - 2016 , para las estaciones de Poal , Vallfogona de Balaguer

y

Raimat fue de 13 , 6ºC , 14 , 3ºC y 13 , 9ºC ,

, Vallfogona de Balaguer y Raimat fue de 13 , 6ºC , 14 , 3ºC

y

13 , 9ºC , respectivamente . La zona se distingue por su aridez para el

suelos desarrollados por el DARP ( Departamento de Agricultura , Ganadería , Pesca

y

Alimentación de la Generalitat de Cataluña ) en la zona de estudio y que ,

, Pesca y Alimentación de la Generalitat de Cataluña ) en la zona de estudio

y

que , posteriormente ( una parte de ellos ) , se han publicado dentro del

cabo de forma coordinada con los técnicos de las Agrupaciones de Defensa Vegetal

y

son una buena representación de los tipos de manejo de las plantaciones en la zona

que estuvieran representados los principales sistemas de riego ( por inundación

y

localizado ) además de los principales tipos de suelos de la zona . Se excluyeron

se estudiaron cincuenta plantaciones . Se ampliaron en el segundo año a sesenta

y

siete , en el tercer año se redujeron a sesenta y cinco y finalmente ,

el segundo año a sesenta y siete , en el tercer año se redujeron a sesenta

y

cinco y finalmente , en el cuarto año , se mantuvieron cincuenta y nueve parcelas

a sesenta y siete , en el tercer año se redujeron a sesenta y cinco

y

finalmente , en el cuarto año , se mantuvieron cincuenta y nueve parcelas . En

a sesenta y cinco y finalmente , en el cuarto año , se mantuvieron cincuenta

y

nueve parcelas . En dos campañas adicionales se validaron los modelos desarrollados

desarrollados ; para ello , en el quinto año se utilizaron 32 nuevas plantaciones

y

80 más el sexto . Estas últimas pertenecían a dos centrales hortofrutícolas (

manzano eran todas del grupo ‘ Golden ’ , que incluye ‘ Golden Delicious ’

y

algunos de sus mutantes . Estaban en plena producción e injertadas mayoritariamente

cubierta vegetal del suelo seguía la práctica común de la zona con la calle enyerbada

y

segada periódicamente , a la vez que se mantenía una banda en la línea de

profundidad enraizable , la clase de drenaje , la disposición de los horizontes

y

el material subyacente . Para diferenciar y posteriormente agrupar las parcelas se

drenaje , la disposición de los horizontes y el material subyacente . Para diferenciar

y

posteriormente agrupar las parcelas se utilizaron los tipos de suelo ( fases de series

( Danés et al . , 1991 ; Herrero et al . , 1993 )

y

en algunos casos criterios adicionales recogidos durante la prospección ( dificultad

freático muy elevado . Se descartaron 10 parcelas , 4 en la elaboración del modelo

y

6 en su validación con datos del sexto año . En la zona elegida

. En la zona elegida de cada plantación , al inicio de los experimentos

y

en el área asociada a la proyección vertical de la cubierta arbórea se obtuvieron

obtuvieron , mediante barrena , muestras compuestas de los primeros 30 cm de profundidad

y

de 30 a 60 cm , o cuando no era posible , hasta la profundidad

material vegetal , el sistema de riego , el sistema de formación , la edad

y

el vigor . El vigor de la plantación se determinó mediante la medida del diámetro

mediante la medida del diámetro del tronco a 20 centímetros sobre el punto de injerto

y

se empleó posteriormente para definir la eficiencia productiva . Se procedió a

anotar la producción , los abonados , los riegos , los tratamientos fitosanitarios

y

otras prácticas culturales realizadas por el agricultor ( datos no mostrados ) .

datos no mostrados ) . Cada año se realizó un análisis mineral de hojas

y

de frutos . Se muestrearon 100 hojas , situadas en el tercio central del brote

metodología descrita por Gagnard ( 1984 ) . En el momento de la cosecha ,

y

de un total de 20 frutos por plantación , se realizaron los análisis del fruto

20 frutos por plantación , se realizaron los análisis del fruto , de la piel

y

de la pulpa sin semillas . En todos los casos se determinaron : nitrógeno ,

determinaron : nitrógeno , fósforo , potasio , calcio , magnesio , boro , zinc

y

manganeso ( AENOR , 2008 ) . La lectura de los elementos se realizó mediante

en atmósfera controlada ( AC ) con el 2 - 3 % de oxígeno

y

un 3 % de dióxido de carbono . Tanto en los frutos que se mantuvieron

frigorífica se valoró la incidencia del bitter pit . En la fase de validación

y

durante la recolección se respetó el procedimiento del agricultor que deposita

de las características edáficas en función de los grupos de suelos establecidos

y

asociados a las distintas plantaciones estudiadas . La metodomodelo consistió en

Con esta selección de variables independientes se realizó un test de paralelismo

y

coincidencia con el objetivo de evaluar si las variables obtenidas en los diferentes

si las variables obtenidas en los diferentes ciclos de cultivo eran coincidentes

y

, por tanto , se podían tratar las variables independientes de los años de estudio

El modelo se elaboró para todas las parcelas ( todos los grupos de suelos )

y

también para una selección de parcelas pertenecientes a los cuatro grupos de suelos

desarrollados sobre materiales detríticos terrígenos finos ( DTF ) en fondos de valles

y

abanicos aluviales , profundos ( >80 cm ) , con alta capacidad de retención de

disponible ( CRAD2 , superior a 1500 m3 ha - 1 ) , bien drenados3

y

con buena estructura . No salinos Suelos desarrollados sobre materiales DTF , profundos

materiales DTF , profundos , con una CRAD superior a 1500 m3 ha - 1

y

con algún problema de drenaje . Presentan dificultades para el enraizamiento debido

drenaje . Presentan dificultades para el enraizamiento debido a las texturas finas

y

a la alta densidad aparente . Incluyen suelos no salinos a ligeramente salinos Suelos

de caliza fragmentada . Se encuentran en abanicos aluviales , terrazas fluviales

y

plataformas estructurales . Son poco profundos ( <0 , 4 m ) , con un

profundos ( <0 , 4 m ) , con un drenaje de bueno a rápido

y

baja CRAD ( <500 m3 ha - 1 ) . No salinos Suelos desarrollados

( <500 m3 ha - 1 ) . No salinos Suelos desarrollados sobre lutitas

y

materiales DTF , profundos y bien drenados . Han sido muy afectados por nivelaciones

1 ) . No salinos Suelos desarrollados sobre lutitas y materiales DTF , profundos

y

bien drenados . Han sido muy afectados por nivelaciones . No salinos o muy ligeramente

detríticos terrígenos con frecuentes elementos gruesos , situados en abanicos aluviales

y

terrazas fluviales , moderadamente profundos , con elevada CRAD , y bien drenados

abanicos aluviales y terrazas fluviales , moderadamente profundos , con elevada CRAD ,

y

bien drenados . No salinos Este grupo se corresponde con los suelos del grupo

existe , en la totalidad de la parcela o en parte , niveles freáticos elevados

y

peor drenaje 1 . Descripción en base a Herrero et al . ( 1993

, con valores de pH que oscilan entre 8 , 1 ± 0 , 2

y

8 , 3 ± 0 , 2 de 0 a 0 , 3 m o

profundidades el contenido de carbonato cálcico equivalente varía entre el 27 ± 6 %

y

el 28 ± 6 % . Los grupos de suelos analizados presentan diferencias significativas

suelos analizados presentan diferencias significativas en la profundidad enraizable ,

y

en los primeros 0 , 3 m de profundidad en los contenidos de arena y arcilla

, y en los primeros 0 , 3 m de profundidad en los contenidos de arena

y

arcilla , la salinidad , los contenidos de materia orgánica y en magnesio intercambiable

contenidos de arena y arcilla , la salinidad , los contenidos de materia orgánica

y

en magnesio intercambiable ( Tabla 2 ) . La arena predomina en suelos de origen

predomina en suelos de origen aluvial ( grupo 4 ) . El contenido de arcilla

y

limo es más elevado en los suelos del grupo 2 que además presentan dificultades

además presentan dificultades de drenaje , contrastando con los de los grupos 4

y

6 , que son suelos muy bien drenados debido a que el material subyacente son

de salinidad ( 0 , 48 dS m - 1 en la prueba previa )

y

de contenidos en magnesio de cambio por tener asociada una capa freática que suele

ser rica en este elemento . Los suelos cuyos materiales subyacentes son lutitas

y

arenitas ( grupos 3 y 5 ) son los que presentan menores contenidos de materia

elemento . Los suelos cuyos materiales subyacentes son lutitas y arenitas ( grupos 3

y

5 ) son los que presentan menores contenidos de materia orgánica ( aunque no difieren

materia orgánica se explica porque se encuentran en superficies con mayor pendiente

y

porque , en general , han sido objeto de movimientos de tierra ( nivelaciones y sistematizaciones

pendiente y porque , en general , han sido objeto de movimientos de tierra ( nivelaciones

y

sistematizaciones ) . Respecto a la fertilidad fósforo - potásica en la mayor parte

Las plantaciones con niveles bajos o muy bajos sólo representan un 6 , 3 %

y

un 13 , 4 % , respectivamente . El predominio de parcelas con contenidos elevados

, respectivamente . El predominio de parcelas con contenidos elevados de fósforo

y

potasio en el suelo se mantiene en las plantaciones empleadas para la validación

la validación ( Tabla 3 ) . Las producciones medias oscilaron entre 9 , 4

y

25 , 1 kg árbol - 1 y las eficiencias productivas entre 0 , 46

producciones medias oscilaron entre 9 , 4 y 25 , 1 kg árbol - 1

y

las eficiencias productivas entre 0 , 46 a 0 , 71 kg fruta cm -

- 1 variaron entre 25 - 35 kg N , 7 - 10 kg P2O5

y

70 - 80 kg K2O . Los grupos de suelos 4 y 6 presentan contenidos

10 kg P2O5 y 70 - 80 kg K2O . Los grupos de suelos 4

y

6 presentan contenidos de Ca en hoja significativamente más elevados pero similares

de Ca al fruto . Además , los valores de los cocientes N / Ca

y

K / Ca son inferiores a los del resto de grupos de suelos ( Tabla

los grupos de suelos en lo que se refiere al contenido de K , Ca

y

Mg en los frutos , pero se observa que los grupos 4 y 6 presentan

, Ca y Mg en los frutos , pero se observa que los grupos 4

y

6 presentan valores de los cocientes N / Ca y K / Ca significativamente inferiores

observa que los grupos 4 y 6 presentan valores de los cocientes N / Ca

y

K / Ca significativamente inferiores a los del resto de grupos ( Tabla 4 )

El cuarto año junto con el primero fueron los de mayor afectación ( Tabla 5 )

y

contrastan con el segundo y tercer años en los que los daños , tras el

el primero fueron los de mayor afectación ( Tabla 5 ) y contrastan con el segundo

y

tercer años en los que los daños , tras el almacenamiento a 20ºC durante 20

1 ) que son los más contrastados , en el periodo de cuatro años evaluado

y

tras 6 meses de conservación en AC , las plantaciones asociadas a los grupos de

meses de conservación en AC , las plantaciones asociadas a los grupos de suelos 1

y

2 presentaron un mayor porcentaje de frutos afectados por bitter pit que los suelos

drenaje ( grupo 2 ) indujeron una mayor afectación que los suelos poco profundos

y

con drenaje bueno o rápido ( grupo 4 ) . Estas diferencias concuerdan con los

4 de suelos ( Tabla 4 ) , mientras que la relación N / Ca

y

K / Ca es menor , excepto el valor de K / Ca para el

ha correlacionado negativamente ( Tabla 6 ) con el contenido de calcio en hoja ,

y

positivamente con las concentraciones de K y B foliar , aunque en alguno de los

con el contenido de calcio en hoja , y positivamente con las concentraciones de K

y

B foliar , aunque en alguno de los años estas correlaciones no han sido significativas

también se ha correlacionado positivamente con los valores de los cocientes K / Ca

y

( K+Mg ) / Ca en hoja y en fruto para todos los años individualmente

valores de los cocientes K / Ca y ( K+Mg ) / Ca en hoja

y

en fruto para todos los años individualmente , y para el conjunto de ellos (

K+Mg ) / Ca en hoja y en fruto para todos los años individualmente ,

y

para el conjunto de ellos ( Tabla 6 ) . Además , en fruto ,

variables que más influyen sobre la incidencia del bitter pit a salida de cámara

y

que no se ven afectadas por otras variables estudiadas , son la eficiencia productiva

por otras variables estudiadas , son la eficiencia productiva , el N en el fruto

y

la relación K / Ca en fruto . Se descartó introducir el boro foliar o

- 2 tronco ) , contenido de N ( en porcentaje sobre materia seca )

y

la relación K / Ca en el fruto . Se ha elaborado para el conjunto

Ca en el fruto . Se ha elaborado para el conjunto de todas las parcelas

y

también de manera separada para cada uno de los cuatro grupos de suelos con características

suelos 2 que son las de mayor incidencia de bitter pit ( Figura 2 )

y

el menor para el grupo de menor incidencia ( grupo 4 ) . En

( Tabla 8 ) . No obstante , cuando el modelo se aplicó a un año

y

tipo de suelo concreto no siempre resultó ser significativo . Los modelos obtenidos

estudios en situaciones reales de manejo en campo , con frutos de fincas comerciales ,

y

tal como han propuesto diversos autores ( Terblanche et al . , 1980 ; Tomala

Terblanche et al . , 1980 ; Tomala et al . , 1993 ; Aichner

y

Drahorad , 2003 ) resultan interesantes . Sin embargo , dada la diversidad de condiciones

resultan interesantes . Sin embargo , dada la diversidad de condiciones de manejo

y

la variabilidad agroclimática es difícil obtener relaciones estadísticas significativas

significativas , tanto en la elaboración del modelo como en su validación ( Tablas 7

y

8 ) . Un mayor número de frutos por árbol tiende a reducir la

crecimiento vegetativo tiene gran influencia en la presencia de fisiopatías ( Ferguson

y

Watkins , 1989 ; Prange et al . , 2011 ) . Las relaciones

. , 2011 ) . Las relaciones encontradas entre los contenidos minerales en fruto

y

el desarrollo de bitter pit también coinciden con otros autores ( Holland , 1980

Johnson , 1989 ; Marcelle , 1990 ; Fallahi et al . , 2010 )

y

especialmente con Sharples ( 1980 ) y Tomala et al . ( 1993 ) respecto a

; Fallahi et al . , 2010 ) y especialmente con Sharples ( 1980 )

y

Tomala et al . ( 1993 ) respecto a la concentración de calcio en fruto .

las parcelas según tipos de suelos ( grupos de suelo , Tabla 1 ) ,

y

el estudio de la incidencia del bitter pit en cada uno de ellos , pone de manifiesto

diferencial muy claro entre ellos ( Figura 2 ) , especialmente entre los grupos 2

y

4 . Los suelos del grupo 2 que presentan la mayor incidencia de bitter pit

bitter pit son profundos , de texturas moderadamente finas ( franco - ar - cillosas

y

franco - arcillo - limosas ) , sin elementos gruesos , con elevada capacidad de

, con elevada capacidad de retención de agua , con algunos problemas de drenaje

y

ligera salinidad . En cambio , en el grupo 4 donde existe la menor incidencia ,

, con baja capacidad de retención de agua , bien drenados o con drenaje rápido

y

no salinos . Las diferencias son evidentes , tanto en la incidencia de bitter pit

pit como en los parámetros considerados como más influyentes ( relaciones N / Ca

y

K / Ca ) . Serían estos factores ( drenaje , capacidad de retención de

/ Ca ) . Serían estos factores ( drenaje , capacidad de retención de agua

y

salinidad ) los determinantes en la aparición del bitter pit , factores que deberían

que deberían ser considerados a la hora de establecer una plantación de manzanos

y

, especialmente , cuando por razones de manejo se establecen variedades sensibles

injertado sobre “ EM9 ” , que favorecen la aparición de esta fisiopatía ( Trocme

y

Gras , 1979 ) . El grupo de suelos 4 , y en menor

( Trocme y Gras , 1979 ) . El grupo de suelos 4 ,

y

en menor medida el 6 , inducen un menor vigor de las plantaciones , lo

las plantaciones , lo que facilitaría un equilibrio más rápido de la plantación

y

una mayor eficiencia productiva . Este conocimiento tiene interés tanto a la hora

establecida , permite también desestimar lotes de fruta para largas conservaciones

y

puede plantear mejoras en el programa de nutrición ( K y P ) y manejo

para largas conservaciones y puede plantear mejoras en el programa de nutrición ( K

y

P ) y manejo del riego en futuras campañas . El modelo propuesto presenta

conservaciones y puede plantear mejoras en el programa de nutrición ( K y P )

y

manejo del riego en futuras campañas . El modelo propuesto presenta la dificultad

presenta la dificultad de que no permite predecir la totalidad de las situaciones

y

además , es necesaria la validación local con un número elevado de parcelas monitorizadas

re - tención de agua moderada , texturas medias , bien o rápidamente drenados ,

y

no salinos , presentan la menor incidencia de bitter pit . En el modelo de

incidencia de bitter pit propuesto , la eficiencia productiva , la relación K / Ca

y

el porcentaje de nitrógeno en el fruto son los factores más significativos . Cuando

determinación más elevado se obtiene en los suelos con mayor incidencia ( grupo 2 )

y

el menor en los de menor incidencia ( grupo 4 ) . No obstante , el

salida de AC próxima al 12 % ( grupo 2 , suelos con texturas finas

y

dificultades de drenaje ) a ser inferior al 3 % ( grupo 4 , suelos

consume de las almendras ( Prunus amygdalus Batsch ) es la semilla , denominada pepita

y

considerada como un alimento de un valor nutritivo y medicinal muy elevado . Las

semilla , denominada pepita y considerada como un alimento de un valor nutritivo

y

medicinal muy elevado . Las almendras se han consumido desde tiempo inmemorial ,

sólo para el consumo humano directo , sino también en las industrias farmacéutica

y

cosmética . Ya en tiempos antiguos se ofrecían a los invitados a todas las bodas

bodas cinco peladillas como símbolos de salud , fertilidad , buena suerte , riqueza

y

felicidad . En sentido contrario , las almendras amargas se han asociado con la mala

. En sentido contrario , las almendras amargas se han asociado con la mala suerte

y

con la muerte ( Verdú et al . , 2017 ) . Por otro lado

registradas en los sistemas médicos más antiguos , como en los antiguos griegos

y

persas , así como en las medicinas tradicionales china y ayurvédica ( Albala ,

en los antiguos griegos y persas , así como en las medicinas tradicionales china

y

ayurvédica ( Albala , 2009 ) . El consumo de almendras ha aumentado en los

hecho que responde a la evolución de su consideración desde la de un aperitivo

y

un elemento de diversos postres , especialmente los navideños , hacia la de un alimento

( Lovejoy et al . , 2002 ; Scott et al . , 2003 )

y

las enfermedades cardiovasculares ( Spiller et al . , 1992 ; Fulgoni et al .

. , 2003 ) . Las almendras se encuentran en un gran número de productos

y

de aplicaciones gastronómicas y se consideran un ingrediente básico de la dieta

almendras se encuentran en un gran número de productos y de aplicaciones gastronómicas

y

se consideran un ingrediente básico de la dieta mediterránea . Se pueden consumir

. Permiten su transformación con el fin de producir turrón , mazapán , helados

y

una gran cantidad de productos de pastelería . Las almendras dan un gusto delicado

y

sabroso y una textura especial a cualquier producto por la calidad de su aceite (

cantidad de productos de pastelería . Las almendras dan un gusto delicado y sabroso

y

una textura especial a cualquier producto por la calidad de su aceite ( Verdú et

importante de calorías , no contiene colesterol sino esteroles de origen vegetal

y

además una elevada proporción de ácidos grasos mono y poli - insaturados frente

esteroles de origen vegetal y además una elevada proporción de ácidos grasos mono

y

poli - insaturados frente a los saturados ( Sabaté y Hook , 1996 ) .

proporción de ácidos grasos mono y poli - insaturados frente a los saturados ( Sabaté

y

Hook , 1996 ) . A pesar de su elevado valor energético , el cuerpo

humano lo compensa con la elevada sensación de saciedad de las almendras ( Hollis

y

Mattes , 2007 ) . Además , la absorción de la energía de las almendras

. Además , la absorción de la energía de las almendras no es muy eficiente

y

su consumo crónico puede aumentar el gasto de energía de reposo ( Cassady et al

fluctuaciones de la glucosa de la sangre ( Jenkins et al . , 2006 )

y

la relación del colesterol TL : HDL y de los triglicéridos ( Foster et al

Jenkins et al . , 2006 ) y la relación del colesterol TL : HDL

y

de los triglicéridos ( Foster et al . , 2012 ) . Considerando todos estos

deben presentar una calidad elevada para satisfacer las exigencias de la industria

y

las preferencias de los consumidores ( Socias i Company et al . , 2008 )

composición química es fundamental al considerar las diferentes aplicaciones industriales

y

la gran diversidad de productos derivados que se pueden obtener . Teniendo en cuenta

porcentaje de ácido oleico en este aceite , la relación entre los ácidos oleico

y

linoleico ( O / L ) , y la concentración en tocoferoles y otras sustancias

, la relación entre los ácidos oleico y linoleico ( O / L ) ,

y

la concentración en tocoferoles y otras sustancias con carácter antioxidante como

ácidos oleico y linoleico ( O / L ) , y la concentración en tocoferoles

y

otras sustancias con carácter antioxidante como los polifenoles ( Socias i Company

antioxidante como los polifenoles ( Socias i Company et al . , 2008 ; Kodad

y

Socias i Company , 2008 ) . La tendencia de las almendras a enranciarse durante

Company , 2008 ) . La tendencia de las almendras a enranciarse durante su almacenaje

y

transporte es una pérdida de calidad relacionada con la oxidación de los ácidos

( Kester et al . , 1993 ) . Por lo tanto , la estabilidad del aceite

y

la composición en ácidos grasos , esencialmente la relación O / L , son unos

punto de vista nutritivo , así como el contenido en antioxidantes como la vitamina E

y

los polifenoles , que previenen la oxidación de los ácidos grasos ( Socias i Company

( Reische et al . , 1998 ) , con diferentes homólogos dependiendo del número

y

la posición de los grupos metilo , presentes en las plantas en cuatro formas ,

en las plantas en cuatro formas , α - , β - , γ -

y

δ - tocoferol . Estos componentes intervienen en diferentes funciones fisiológicas

Estos componentes intervienen en diferentes funciones fisiológicas , biológicas

y

bioquímicas , sobre todo por su acción antioxidante , pero también como estabilizadores

su acción antioxidante , pero también como estabilizadores de la membrana ( Azzi

y

Stocker , 2000 ) . Su fun ción bioquímica principal parece ser la de proteger

proteger los ácidos grasos poli - insaturados contra la peroxidación ( Kamal - Eldin

y

Appelqvist , 1996 ) , además de tener un papel pro tector en los sistemas

los sistemas biológicos , con propiedades hipo - colesterólicas , anti - cáncer

y

neuro - protectivas ( Sen et al . , 2007 ) . La concentración en

almendra un papel muy importante en la protección de los lípidos contra la oxidación

y

, por lo tanto , su posibilidad de almacenaje ( García - Pascual et al

( Kodad , 2017 ) , lo que puede determinar la utilización de las almendras

y

la predicción de su vida útil de almacenaje . Al considerar la variación en el

contenido superior , pueden tener más capacidad de almacenaje que los otros ( Fourie

y

Basson , 1989 ) . Por lo tanto , el objetivo de este trabajo ha sido la

cultivares mallorquines de almendro en cuanto a su composición en ácidos grasos

y

en tocoferoles con el fin de establecer su valoración cualitativa , así como un

producto más o menos diferencial entre los otros cultivares de almendro . Material

y

métodos Material vegetal Se estudió un total de 23 cultivares mallorquines de almendro

establecida en la Granja Experimental de la Ciutat de Mallorca durante los años 1950

y

1960 , además de introducciones posteriores llevadas a cabo por el técnico Joan

de eliminar el efecto de variabilidad debido a las condiciones ambientales ( Kodad

y

Socias i Company , 2008 ) . Las almendras se recogieron durante dos años consecutivos

2010 - 2011 ) una vez maduras , cuando el mesocarpo ya estaba totalmente seco

y

abierto y la abscisión del pedúnculo era fácil ( Felipe , 1977 ) . Cada

2011 ) una vez maduras , cuando el mesocarpo ya estaba totalmente seco y abierto

y

la abscisión del pedúnculo era fácil ( Felipe , 1977 ) . Cada año se

almendras de los tres árboles diferentes de cada cultivar . Determinación del aceite

y

de los ácidos grasos El contenido en aceite y la proporción de los diferentes ácidos

cultivar . Determinación del aceite y de los ácidos grasos El contenido en aceite

y

la proporción de los diferentes ácidos grasos se determinó después de pelar les

grasos se determinó después de pelar les almendras por escaldado , dejarlas secar

y

molerlas en un molinillo de cocina hasta obtener una harina de grano fino . El

éter de petróleo como solvente ( rango de ebullición de 40 - 60º ) ,

y

manteniendo la fuente de calor a 135ºC . La extracción duró 2 horas , ya

a 135ºC . La extracción duró 2 horas , ya que experiencias previas ( Kodad

y

Socias i Company , 2008 ) habían demostrado que con este tiempo la extracción era

expresó como la diferencia entre el peso seco ( PS ) de la muestra antes

y

después del proceso de extracción ( % de PS ) . El porcentaje relativo de

( ISO , 2000 ) . Estos EMAGs se separaron en un cromatógrafo HP 6890

y

después se detectaron con un detector ionizante de llama equipado con una columna

capilar ( HP - Innowax × 0 , 25 mm i . d . )

y

un film de 0 , 25 μm de grosor ( Agilent Technologies , Waldbronn ,

una velocidad de flujo de 1 mL min - 1 . La temperatura del inyector

y

el detector fueron de 220 y 275ºC respectivamente . La temperatura inicial de la

mL min - 1 . La temperatura del inyector y el detector fueron de 220

y

275ºC respectivamente . La temperatura inicial de la columna fue de 100ºC durante

1 , de 150 a 200ºC a un aumento de 3º min - 1 ,

y

desde 200 a un aumento de 13º min - 1 , manteniéndose a 240º durante

, España ) . Ca da año se calcularon las medias de las tres determinaciones

y

los valores de los diferentes ácidos grasos se expresaron como porcentajes sobre

alto rendimiento ( CLAR o HPLC ) , equipada con una bomba de doble pistón

y

de auto - muestreo . Se utilizó una columna Luna 100RP - 18 de ×

μL del extracto propanólico se inyectaron a la HPLC . Las cantidades de δ -

y

γ - tocoferol se midieron con un detector de fluorescencia FLD SFM25 ( Kontron )

SFM25 ( Kontron ) , calibrado a longitud de onda de excitación de 295 nm

y

de emisión de 325 nm . La cantidad de α - tocoferol se midió con

de los cultivares mallorquines se compararon con los de otras regiones españolas

y

de países extranjeros con los datos de Kodad y Socias i Company ( 2008 )

los de otras regiones españolas y de países extranjeros con los datos de Kodad

y

Socias i Company ( 2008 ) y de Kodad et al . ( 2011 y

de países extranjeros con los datos de Kodad y Socias i Company ( 2008 )

y

de Kodad et al . ( 2011 y 2014 ) . Resultados Variabilidad Todos los

y Socias i Company ( 2008 ) y de Kodad et al . ( 2011

y

2014 ) . Resultados Variabilidad Todos los componentes estudiados mostraron una gran

diferentes cultivares mallorquines de almendro osciló entre 47 , 40 % de ‘ Guarim ’

y

61 , 62 % de ‘ Andreu ’ , con un promedio para el conjunto

procedencia , que muestran una variación global del contenido en aceite entre 20

y

68 % del PS , incluyendo Argentina , California , Egipto , España , Francia

Grecia , India , Irán , Italia , Marruecos , Portugal , Serbia , Túnez

y

Turquía ( Kodad , 2017 ) . Ácidos grasos En relación a la composición del aceite

, el gadoleico ( C20 : 1 ) con un 0 , 08 % ,

y

el behénico ( C22 : 0 ) con un 0 , 01 % . El

cultivares mallorquines analizados osciló entre el 58 , 65 % en ‘ Vinagrilla ’

y

el 78 , 44 % en ‘ Corona ’ , con una media del 69

2 n6 ) que osciló entre el 12 , 01 % en ‘ Corona ’

y

el 29 , 40 % en ‘ Vinagrilla ’ , y una media del 21

en ‘ Corona ’ y el 29 , 40 % en ‘ Vinagrilla ’ ,

y

una media del 21 , 14 % , confirmando la correlación negati va entre el

ácidos grasos ( Yada et al . , 2011 ) al encontrarse los valores máximos

y

mínimos en los mismos cultivares . Ya que la relación entre los ácidos oleico

y

linoleico ( O / L ) es esencial para el mantenimiento de la estabilidad del

valor se considera un criterio importante para evaluar la calidad de las almendras ,

y

muestra los valores extremos en los mismos cultivares : 6 , 53 en ‘ Corona

los valores extremos en los mismos cultivares : 6 , 53 en ‘ Corona ’

y

1 , 99 en ‘ Vinagrilla ’ , con una media de 3 , 27

de 2 , 38 % para el ácido esteárico ( C18 : 0 ) ;

y

de 0 , 38 % en ‘ Torres ’ a 0 , 61 % en

analizadas osciló entre 258 , 5 mg kg - 1 aceite para ‘ Alzina ’

y

500 , 2 mg kg - 1 aceite para ‘ Bord de Santa Maria ’

kg - 1 aceite . Estos dos mismos cultivares presentan los mismos valores mínimos

y

máximos para los otros homólogos de tocoferol : 0 , 02 mg kg - 1

de tocoferol : 0 , 02 mg kg - 1 aceite para ‘ Alzina ’

y

1 , 13 mg kg - 1 aceite para ‘ Bord de Santa Maria ’

tocoferoles osciló entre 260 , 41 mg kg - 1 aceite para ‘ Alzina ’

y

522 , 13 para ‘ Bord de Santa Maria ’ , teniendo en cuenta el

, el porcentaje de ácido oleico en el aceite , la relación O / L

y

la cantidad de anti - oxidantes ( Socias i Company et al . , 2008

así como la extracción de aceite , especialmente para las industrias farmacéutica

y

cosmética . Almendras con un porcentaje bajo de aceite facilitan la molienda al

con un porcentaje bajo de aceite facilitan la molienda al no formar tantos grumos ,

y

así se obtiene más fácilmente la harina de almendra o la leche de almendra ,

para los californianos , 57 % a 63 , 9 % para los ucranianos ,

y

63 % a 66 % para los indios , con unas medias similares ( Tabla

Tabla 5 ) . Cabe remarcar que la mayoría de cultivares de Francia , Italia

y

los Estados Unidos tienen un contenido en aceite superior al 60 % , mientras que

superior al 60 % , como ‘ Andreu ’ , ‘ Bord del Raiguer ’

y

‘ Corona ’ , podrían servir para producir turrón blando , aunque en otras regiones

Por otro lado , cultivares con los porcentajes menores de aceite , como ‘ Guarim ’

y

‘ Vinagrilla ’ , serían más apropiados para la obtención de harina o de leche

% ) , España ( del 67 , 42 al 74 , 92 % )

y

Portugal ( del 67 , 66 al 74 , 10 % ) ( Kodad et

resultando en la producción de diversos productos que afectan a la calidad ( Sung

y

Jeng , 1994 ) , incluido el gusto rancio ( Harris et al . ,

grasos insaturados , la luz , el oxígeno , los iones metálicos , la temperatura

y

las enzimas ( Gou et al . , 2000 ; Zacheo et al . ,

relación de ácidos grasos insaturados / ácidos grasos saturados ( Serra Bonvehí

y

Ventura Coll , 1993 ; Pershern et al . , 1995 ) . Aunque en general

) . En todo caso hay algún cultivar que puede destacar en un componente concreto

y

que sea de interés para poder establecer un uso determinado de sus almendras . Sí

. Sí puede considerarse que en general el contenido en aceite , en ácido oleico

y

en tocoferoles de los cultivares mallorquines es inferior a la media de las otras

-5	-4	-3	-2	-1	Keyword	+1	+2	+3	+4	+5

Concordances from the Ingeniería corpus