Agronoma Trop. 62(1 - 4): 137-144.

2012

Comparacin de mtodos de estimacin

de la radiacin solar en Maracay, Venezuela
Comparison of methods to estimate the solar radiation
in Maracay, Venezuela
Vannessa G. Aylln C.1
1

Ingeniero Agrnomo. Tesis de Pregrado. Universidad Central de Venezuela (UCV). Facultad de Agronoma. Maracay 2101,
estado Aragua. Venezuela. Correo electrnico: [email protected]

RESUMEN

ABSTRACT

La radiacin solar global es el parmetro de entrada

para numerosas aplicaciones de ndole agronmico,
sin embargo, los datos estn poco disponibles para
el usuario, por esto, surge la necesidad de utilizar
mtodos empricos para estimarla que manejen
variables de fcil consecucin disponibles en todas
las estaciones. Utilizando datos diarios de la Estacin
Meteorolgica del CENIAP-Maracay, serial 2604,
perteneciente a la red agrometeorolgica del Instituto
Nacional de Investigaciones Agrcolas, estado Aragua,
Venezuela, se realiz un estudio comparativo en cinco
mtodos de estimacin: Angstrom-Prescott (1940);
Glover y McCulloch (1958); Bristow y Campbell (1984);
Hargreaves et al. (1985); De Jong y Stewart (1993),
con la finalidad de seleccionar el que estime mejor
la radiacin solar global para la zona en estudio. Las
comparaciones se ejecutaron a travs de la prueba
estadstica de prueba ji cuadrado, encontrndose
que, todos los mtodos estn en concordancia con
los valores de radiacin medidos en esta estacin.
Sin embargo, el mtodo que present la menor
sobreestimacin, fue De Jong y Stewart.

The solar radiation is the input parameter for numerous

applications, however, hardly available to the user, this
arises the need to use empirical methods to estimate
it which use easy achievement variables available
in all meteorological stations. Using daily data from
the meteorological station of Ceniap Maracay serial
2604, belonging to the agrometeorological network of
INIA in Maracay, Venezuela, a study was conducted
to compare 5 estimation methods: Angstrom-Prescott
(1940), Glover and McCulloch (1958), Bristow and
Campbell (1984), Hargreaves et al. (1985) and De
Jong and Stewart (1993); in order to select which of
these it considers better estimator of solar radiation for
the study area. Comparisons have been implemented
through the statistical test chi square, finding that,
statistically, all methods are consistent with the values
of radiation measured for this station. However, the
method presented by the minor overestimation, was
found to be De Jong and Stewart.

Key words: global solar radiation, AngstromPrescott, Glover and McCulloch, Bristow and
Campbell, Hargreaves et al., De Jong and Stewart.

Palabras clave: radiacin solar global, AngstromPrescott, Glover y McCulloch, Bristow y Campbell,
Hargreaves, De Jong y Stewart.

Recibido: 26/04/2013

Aprobado: 05/11/2014

137

Vol. 62 (1 - 4)

AGRONOMA TROPICAL

2012

en las estaciones climatolgicas, con el objeto

de comparar los datos estimados de radiacin
global con los medidos, y as determinar cul de
los mtodos es el que mejor la estima para la
zona en estudio.

INTRODUCCIN
El conocimiento de la radiacin solar global es
esencial para muchas aplicaciones, tales como
el uso de la energa solar como fuente alterna,
modelos de simulacin de crecimiento de los
cultivos, en la arquitectura en el diseo de
sistemas, estimacin de la evapotranspiracin,
entre otras. Sin embargo, la disponibilidad
de datos de radiacin solar global es escasa,
haciendo que el uso de tcnicas numricas
sea una alternativa esencial. Con tales tcnicas
indirectas, otros datos meteorolgicos son
matemticamente explorados con el fin de
estimar la cantidad de radiacin solar global que
llega a la tierra (Tymvios et al., 2004).

MATERIALES Y MTODOS
Para el estudio se utiliz una serie de datos
climatolgicos comprendidos desde el ao 1971
hasta 2010 de la estacin del Centro Nacional
de Investigaciones Agropecuarias (CENIAP)
Maracay serial 2604, perteneciente a la red
agrometeorolgica del Instituto Nacional de
Investigaciones Agrcolas (INIA), ubicada en
Maracay, estado Aragua, latitud: 10 17' N,
longitud: 67 37 W, altitud: 455 m s. n. m.

Aunque se han desarrollado varios modelos

empricos para la estimacin de la radiacin
global, por regla general, los datos que se
requieren para su aplicacin son bsicamente
los mismos. Los parmetros meteorolgicos
que normalmente se utilizan, son: la humedad
relativa, la temperatura, la altitud, la nubosidad,
as como la insolacin relativa.

Para realizar el procesamiento de datos, estos

se organizaron en forma diaria tal como se
muestra en el Cuadro 1.
Posteriormente se evalu la calidad de los datos
a travs de la siguiente metodologa:
Datos Faltantes: se eliminaron aquellos aos en
los que existan ms de 10 datos diarios faltantes
para cualquier elemento del clima, siguiendo los
lineamientos establecidos por la Organizacin
Meteorolgica Mundial (Lpez, 2009).

Los modelos que ms aceptacin han tenido

son aquellos que consideran nicamente la
insolacin relativa o bien la cantidad de nubes,
ya que, son simples expresiones de regresin
del proceso real de transferencia de radiacin
en la atmsfera (Estrada y Almanza, 2005).

Deteccin de datos atpicos: para radiacin

global, se compar cada serie con los datos
correspondientes de radiacin extraterrestre,
teniendo en cuenta que debe encontrarse entre
el 10% y 85% (Atlas de Radiacin Solar, 2002).

En este estudio se estim la radiacin global a

travs de diferentes metodologas que emplean
otras variables climticas medidas comnmente

Cuadro 1. Disposicin de los elementos climticos para el mes de enero.

Ao

Da

Precipitacin
(mm)

1971

1971

31

2010
2010

Radiacin
(MJ m-2 d-1)

Insolacin
(h)

Tmx
(C)

Tmn
(C)

16,8269

6,0

30,2

17,0

18,2082

8,9

31,8

14,8

18,5012

3,7

33,6

19,6

31

21,4732

8,7

33,8

14,6

138

Aylln - Comparacin de mtodos de estimacin de la radiacin solar...

Segn Castro (2008) el rango propuesto para las

observaciones de temperatura del aire es de 5 a
40 C para zonas bajas tropicales. En el caso del
porcentaje de humedad relativa, considerando
los valores promedios mensuales para la zona a
estudiar, se consideran valores atpicos aquellos
que sean menores a 30%. Para insolacin se
consideran valores atpicos aquellos que sean
mayores a la duracin del da astronmico (N)2.

Donde: Ra: radiacin extraterrestre (MJ m-2 da-1);

Gsc: constante solar = 0,0820 (MJ m-2 da-1); dr: inverso de la distancia relativa entre la Tierra y el Sol;
s: ngulo de la puesta de sol (radianes); : latitud (radianes); : declinacin solar (radianes).

Homogeneidad de las series climatolgicas

dr = 1 + 0,033 cos

Para verificar la homogeneidad de las series

climatolgicas se utiliz la prueba de rachas,
segn lo propone Prez et al. (2005), se corri la
prueba a travs del programa estadstico Statistix
versin 8.0 con un nivel de significacin del 5%,
se emplearon los promedios mensuales, ya que
podra existir cierta dependencia o relacin en
los valores diarios de las variables en estudio3.

Donde J es el nmero del da del ao (1 a 365)

El inverso de la distancia relativa entre la Tierra

y el Sol (dr), y la declinacin () estn dados por
las siguientes ecuaciones:
2

( )
365

= 0,409 sen

( ) - 1,39

365

El ngulo de radiacin a la hora de la puesta del

sol, s, viene dado por:
s = arc cos [- tan () tan ()]
Mtodo de Glover y McCulloch (1958)
Rs

Para la estimacin de la radiacin se utilizaron

cinco mtodos:

Ra

Mtodo Angstrom-Prescott (1940)

= 0,29 cos(lat) + 0,55 (n/N)

Donde Rs, Ra, n, N, n/N, lat., son similares a

los indicados previamente. Lat. es la latitud
radianes. El mtodo es vlido para latitudes de
0 a 60.
Mtodo de Bristow y Campbell (1984)

Donde: Rs: radiacin global (MJ m da );

Ra: radiacin extraterrestre (MJ m-2 da-1); n:
insolacin (h); N: duracin mxima posible de
la insolacin (h); n/N: duracin relativa de la
insolacin (adimensional); a y b: coeficientes
empricos que se pueden estimar para una zona
en particular.
-2

-1

Rs
Ra

Donde Rs y Ra, son similares a lo indicado

previamente, mientras que A, B y C son
coeficientes empricos que se pueden estimar
para un sitio en particular a partir de datos
medidos de radiacin solar y D: rango diurno de
la temperatura del aire, se calcula:

Se fijaron los valores de los coeficientes a y b

como 0,2 y 0,5, respectivamente, tal como lo
recomiendan Prez et al. (2005) cuando no se
dispone de datos de Rs y Ra.

D = Tmx -

La radiacin extraterrestre se calcul de la

siguiente manera (Allen et al., 1998):
Ra =

2
3

24(60)

= A(1 - exp [-B DC])

Tmn (j) + Tmn (j + 1)

2

Donde: Tmx: temperatura mxima diaria (C);

Tmn(j): temperatura mnima en un da (C);
Tmn(j+1): temperatura mnima diaria al da
siguiente.

Gscdr[ssen()sen()+cos()cos()sen(s)]

Ingeniera Agrnomo Mara Jos Mndez, comunicacin personal - [email protected]

Ingeniera Agrnomo Genady Bracho, comunicacin personal - [email protected]

139

Vol. 62 (1 - 4)

AGRONOMA TROPICAL

Para el presente estudio se establecieron los

coeficientes de acuerdo con las directrices
utilizadas por dichos autores, donde A y C se
mantuvieron constantes con valores de 0,7 y
2,4, respectivamente y B variable con el objeto
de estudiar la estacionalidad de los datos. El
valor de B est directamente relacionado con la
amplitud trmica diaria (T) y se calcula de la
siguiente manera:

RESULTADOS Y DISCUSIN
Con respecto al anlisis de calidad de los datos,
despus de la eliminacin de datos atpicos se
cumpli con la longitud adecuada de las series
(30 aos) para asegurar la confiabilidad de los
resultados (Prez et al., 2005).
La no uniformidad se present en la precipitacin
para el mes de marzo, la radiacin global para
los meses enero, febrero, marzo, abril, mayo y
diciembre; y por ltimo la temperatura mnima
en el mes de octubre. Para el resto de los
elementos climticos en estudio si se obtuvo
homogeneidad.

B = 0,036 exp (-0,154 T)

Donde:
T = Tmx - Tmn

Mtodo de Hargreaves et al. (1985)

Rs
Ra

Segn Prez et al. (2005), los valores de

probabilidad por debajo del nivel de significacin
indican que los trminos de la serie tienen una
alta dependencia positiva y por lo tanto, se
sospecha que la misma tiene una tendencia a
aumentar o disminuir.

= KT T

Donde: KT: coeficiente de calibracin que puede

ser estimado para una zona en particular; T:
diferencia entre Tmx y Tmn diarias.

Referente a las estimaciones de los coeficientes

del modelo de De Jong y Stewart (1993), el
Cuadro 2 muestra los valores obtenidos y los
estadsticos calculados.

El coeficiente de calibracin (KT) que se utiliz

es de 0,134, segn seala Lpez (2009).
Mtodo de De Jong y Stewart (1993)
Rs = Ra[a(T)b (1 + cP + dP2)]

Los bajos valores del coeficiente de determinacin (R2) indican una deficiencia en el
modelo, sin embargo, al observar los valores
de los cuadrados medios del error (RMSE),
resultaron de baja magnitud, lo que indicara
poca discrepancia entre los valores observados
y los medidos.

Donde Rs, Ra, T son smilares a lo indicado

anteriormente y a, b, c, d son coeficientes
empricos que varan segn la poca del ao.
Los coeficientes se estimaron a travs de ajustes
no lineales, por medio del programa estadstico
JMP versin 6.0, luego para verificar la
bondad del ajuste y siguiendo lo propuesto por
los autores, se realiz el anlisis de regresin
para la variable Y (Rs estimada/Ra) y la variable
X (Rs medida/Ra) obtenindose los valores de
los coeficientes de determinacin (R2) y errores
cuadrticos medios (RMSE) para cada mes.

Es importante destacar que decidir la bondad

del modelo solamente basado en el valor de R2
podra llegar a establecer conclusiones erradas.
Comparacin de las estimaciones
Segn lo observado en la Figura 1, las estimaciones hechas por el mtodo de De Jong y Stewart
(1993) son las ms precisas, exceptuando
en el mes de julio, mientras que para el resto
de los mtodos se observa que Bristow y
Campbell (1984) y Glover y McCulloch (1958)
sobreestiman la radiacin global, en cambio que
Angstrom-Prescott (1940) y Hargreaves et al.
(1985) tienden a subestimar la radiacin.

Por ltimo, se realiz la prueba de ji cuadrado

para comparar las estimaciones de cada mtodo
(E) con los valores reales de radiacin solar
global (O). El estadstico de prueba se basa en
las diferencias entre los Oi y Ei y se define como:
k

2 =

i=1

2012

(Oi - Ei)2
Ei

140

Aylln - Comparacin de mtodos de estimacin de la radiacin solar...

Cuadro 2. Coeficientes mensuales y estadsticos para el modelo de De Jong y Stewart.

Mes

R2

RMSE

Enero

0,544

0,013

-0,032

0,00172

0,008

0,104

Febrero

0,287

0,244

-0,068

0,00335

0,087

0,106

Marzo

0,209

0,352

-0,043

0,00137

0,177

0,098

Abril

0,140

0,483

-0,008

0,00010

0,316

0,105

Mayo

0,165

0,425

-0,005

0,00005

0,233

0,107

Junio

0,230

0,313

-0,006

0,00005

0,161

0,108

Julio

0,251

0,284

-0,006

0,00005

0,135

0,104

Agosto

0,213

0,346

-0,006

0,00007

0,148

0,488

Septiembre

0,302

0,216

-0,005

0,00008

0,062

0,114

Octubre

0,300

0,210

-0,005

0,00005

0,064

0,107

Noviembre

0,386

0,118

-0,008

0,00020

0,029

0,113

Diciembre

0,361

0,153

-0,022

0,00051

0,047

0,105

Figura 1. Radiacin global medida vs. estimada.

Como se observa, durante el ao 1987 los

valores estimados por el mtodo de AngstromPrescott (1940) resultaron ms precisos en los
meses secos, inclusive en aquellos casos donde
existen valores bajos de radiacin, no obstante,
el mtodo requiere datos medidos de radiacin
global para estimar los coeficientes a y b de
la zona en particular. Por otra parte, el mtodo
de Glover y McCulloch (1958), sobreestima la

Este resultado podra deberse al efecto de solapamiento causado por las medias mensuales,
para verificar dicho efecto, se compararon las
estimaciones diarias usando el ao 1987, por
ser el nico ao en el que no existen datos
faltantes para los 365 das. En contraste, se presentan en las Figuras 2 y 3 las estimaciones
diarias del mes de enero (mes seco) con las del
mes de agosto (mes hmedo), respectivamente.
141

Vol. 62 (1 - 4)

AGRONOMA TROPICAL

2012

una desventaja con respecto a los otros mtodos

que requieren de datos medidos en estaciones
ms sencillas. En el caso del mtodo de Bristow
y Campbell (1984) su principal debilidad se
debe a que, son asumidos los valores de los
coeficientes A y C del estudio original, por lo
que no se tienen valores caractersticos de
la zona en estudio, generando errores en las
estimaciones. El mtodo de Hargreaves et al.
(1985) requiere de datos medidos de radiacin
solar global para estimar el coeficiente KT de la
zona en particular.

radiacin tanto en estimaciones diarias como

en las mensuales, en contraste con el mtodo
de Angstrom-Prescott (1940) y Hargreaves
et al. (1985) no son necesarios datos de
radiacin global para la estimacin y su principal
desventaja es que solamente es vlido para
latitudes entre 0 y 60.
Tanto en el mtodo Angstrom-Prescott (1940)
como el de Glover y McCulloch (1958) son
necesarios datos de insolacin, que solamente
se miden en estaciones tipo C1, lo que implica

Figura 2. Radiacin global medida vs. estimada para el mes de enero.

Figura 3. Radiacin global medida vs. estimada para el mes de agosto.

142

Aylln - Comparacin de mtodos de estimacin de la radiacin solar...

En comparacin con los resultados obtenidos

en la Figura 1, donde el mejor estimador es el
modelo de De Jong y Stewart (1993), se observa
que en los datos diarios esto se debe ms a
la uniformidad de las estimaciones, que a la
precisin de las mismas, puesto que los valores
estimados diarios por encima y por debajo
de los medidos, se compensan al calcular los
promedios mensuales.

la sobreestim para todos los casos como se

mencion anteriormente, contrario a esto, se
tiene que el mtodo de De Jong y Stewart (1993)
es el que mejor estim la radiacin, obtenindose
el valor mnimo de ji cuadrado para el mes de
enero y el mximo en el mes de mayo.
Considerando los resultados obtenidos del
anlisis de regresin, donde los valores de los
coeficientes de determinacin mensuales son
relativamente bajos al igual que los cuadrados
medios del error, el modelo de De Jong y Stewart
(1993) es superior a los otros por la calidad de las
estimaciones. La principal desventaja es que la
obtencin de los coeficientes se realiza a travs
de metodologas que pueden resultar complejas
en usuarios que no estn familiarizados con
modelos no lineales.

Segn se observa en el Cuadro 3, los resultados

obtenidos en la prueba de ji cuadrado, para todos
los mtodos de estimacin los valores estimados
estn en concordancia con los valores medidos,
sin embargo existen diferencias entre los valores
de ji cuadrado para los distintos mtodos. El
mtodo de Bristow y Campbell (1984) fue el
que peor estim la radiacin global en cuanto

Cuadro 3. Valores de ji-Cuadrado para cada mtodo.

AngstromPrescott

Glover y
McCulloch

Bristow y
Campbell

Hargreaves et al.

De Jong y Stewart

Enero

0,71

0,86

0,85

0,85

0,63

Febrero

0,68

0,87

0,94

0,81

0,69

Marzo

0,81

0,95

1,06

0,73

0,66

Abril

1,36

1,02

2,64

0,94

0,94

Mayo

1,18

0,90

2,37

1,07

1,02

Junio

1,60

0,91

3,79

1,13

0,92

Julio

1,34

0,83

1,65

1,11

0,84

Agosto

1,25

0,90

1,78

1,08

0,91

Septiembre

1,43

0,94

1,87

1,36

0,98

Octubre

1,10

0,92

1,56

1,02

0,82

Noviembre

1,10

1,01

1,32

1,12

0,82

Diciembre

0,96

0,91

0,99

0,83

0,64

Media

1,13

0,92

1,73

1,00

0,82

143

Vol. 62 (1 - 4)

AGRONOMA TROPICAL

2012

Cientficos de la Selva y Manejo de Informacin.

Compilado por Enrique Castro. Revisado por:
D. Brenes, W. Lpez, O. Madrigal, B. Matarrita,
J. Serrano y Z. Zahawi. Sarapiqu, Costa Rica.

CONCLUSIONES
Evaluar la calidad de los datos resulta fundamental en cualquier anlisis climatolgico,
puesto que de ello depender la bondad de los
resultados.

De Jong, R. and D. W. Stewart. 1993. Estimating

global solar radiation from common meteorological observations in western Canada.
Canadian Journal of Plant Science (Canad).
73(2):509-518.

El uso futuro de un mtodo en particular, depender del alcance que se persiga, si se desea
trabajar con datos mensuales, el mtodo de De
Jong y Stewart (1993) proporcionar la mejor
herramienta para la estimacin, mientras que si
se trabajara con datos diarios, es recomendable
el mtodo de Angstrom-Prescott (1940). Por
otra parte, se deber evaluar si de acuerdo al
propsito de la investigacin conviene ms un
mtodo que sobreestime o subestime la radiacin
solar global, porque de ello depender su impacto
en otras variables.

Estrada, V. y R. Almanza. 2005. Irradiaciones global,

directa y difusa en superficies horizontales
e inclinadas, as como irradiacin directa
normal, en la Repblica Mexicana. Universidad
Nacional Autnoma de Mxico, Mxico. 64 p.
Lpez J. 2009. Efecto y estimacin de datos
diarios faltantes de radiacin global sobre
la calidad de la serie mensual. Trabajo de
grado Ingeniero Agrnomo. Maracay, Ven.
Universidad Central de Venezuela. Facultad
de Agronoma. 57 p.

Los coeficientes obtenidos en los distintos

modelos, exceptuando el de Angstrom-Prescott
son particulares de la zona en estudio, por lo que
no es recomendable utilizarlos en regiones con
caractersticas dismiles.

Hargreaves, G. L., G. H. Hargreaves and J. Riley.

1985. Irrigation water requirement for the
Senegal River Basin. Journal of Irrigation and
Drainage Engineering. 111(3):265-75.

LITERATURA CITADA
Allen, R., L. Pereira, D. Raes and M. Smith.
1998. Crop evapotranspiration. Guidelines for
computing crop water requirements. FAO-56.
Roma. 301 p.

Prez, M., M. Puche y G. Bracho. 2005.

Herramientas para el control de calidad de
datos climticos y su anlisis con fines de uso
agropecuarios. Revista Digital CENIAP HOY.
Revista tcnica de las ciencias del agro editada
por el CENIAP del INIA. Seccin Preimpresos.

Atlas de Radiacin Solar. 2002. Instituto de

hidrologa, meteorologa y estudios ambientales IDEAM. Unidad de planeacin minero
energtica UPME. Repblica de Colombia.
Apndice. D.

Prescott, J. 1940. Evaporation from a Water

Surface in Relation to Solar Radiation.
Quarterly Journal of the Royal Meteorological
Society. 64(1):114-118.
Tymvios, F. S., C. P. Jacovides, S. C. Michaelides
and C. Scouteli. 2005. Comparative study
of Angstroms and Artificial Neural Networks
methodologies in estimating global solar
radiation. Solar Energy (Cyprus). 78(6):752-762.

Bristow, K. and G. Campbell. 1984. On the relationship between incoming solar radiation and
daily maximum and minimum temperature.
Agricultural and Forest Meteorology (EE.
UU.).31(2):159-166.
Castro, F. 2008. Organizacin para Estudios
Tropicales. Manual de procedimientos para
las estaciones meteorolgicas. Departamentos

144