Radiacin Solar

CAPITULO 2

RADIACION SOLAR
2.1.

Generalidades

La fuente de energa de todos los procesos fsicos y biolgicos que ocurren sobre la superficie terrestre
proviene del Sol. Monteith (1958) mencion que la agricultura era una explotacin de la energa solar,
hecha posible por un adecuado suministro de agua y nutrientes para mantener el desarrollo de las
plantas.
De acuerdo con la Termodinmica se conoce que la energa se transforma pero no se pierde y el
ejemplo tal vez, ms significativo en la agricultura es la fotosntesis, cuya ecuacin qumica en su
forma ms simple es:
CO2 + H2O + energa (CH2O) + O2 - 112,000 caloras
en la que la energa proviene del sol.
Best (1962) clasific a los efectos de la radiacin solar sobre las plantas verdes en dos tipos de
procesos, de la siguiente forma:
I. Procesos Foto-energticos: como la Fotosntesis
II. Procesos Foto-estimulantes:
A. Procesos de Movimiento: como el Fototropismo
B. Procesos Formativos: como Elongacin de Estomas, Expansin de las
Hojas, Formacin de Pigmentos, etc.
Segn Brown y Escombe, citados por Demolon (1966), s la energa incidente fuera descompuesta
porcentualmente, su distribucin seria la siguiente:
Energa Consumida por la Fotosntesis
Energa Consumida por la Transpiracin
Energa Transmitida por las Hojas
Energa Perdida por Irradiacin

0.66 %
48.39 %
31.40 %
19.55 %

De lo anterior se puede inferir la significacin del estudio de la Radiacin Solar.

2.2.

Caractersticas de la Radiacin Solar

La radiacin solar recibida en la superficie de la tierra, es el origen de casi todos los fenmenos
meteorolgicos y de sus variaciones en el curso de los das y de los aos.
La radiacin es un proceso fsico por medio del cual se transmite energa en forma de ondas
electromagnticas.

Radiacin Solar

La radiacin solar que llega del lmite superior de la atmsfera est formada por rayos de diferentes
longitudes de onda, principalmente por:
a) Rayos Luminosos, que son visibles y su longitud de onda vara entre 0.36 a 0.76 micrones
(1 mm = 1000 micrones);
b) Rayos Ultravioleta o Qumicos, no son visibles y su longitud de onda es pequea, menor
de 0.36 micrones; y
c) Rayos Infrarrojos o Calorficos, tampoco son visibles y su longitud de onda es mayor a
0.76 micrones.
Los rayos ultravioleta causan efectos detrimentales o germicidas sobre las plantas, los rayos luminosos
intervienen en muchos procesos metablicos, principalmente en el fototropismo (0.36 a 0.49) y la
fotosntesis (0.49 a 0.76) y sobre los rayos infrarrojos no se conocen efectos especficos sobre las
plantas, aparentemente son absorbidos y transformados en calor sin interferir con los procesos
bioqumicos.
2.3.

Constante Solar

La Constante Solar se define como el flujo de energa proveniente del sol, que incide en una
superficie perpendicular a la direccin de propagacin de la radiacin solar, fuera de la atmsfera
terrestre y a una distancia media de la tierra al sol.
La tierra gira alrededor del sol en una rbita elptica y el sol esta ubicado en un foco, como se ilustra en
la Figura 2.1.

152 x

147 x 106
km

Figura 2.1. Distancia de la Tierra al Sol.

La distancia ms cercana de la tierra al sol es de 147 millones de Km y ocurre aproximadamente el da

3 de enero, conocida como PERIHELIO y como AFELIO a la distancia ms lejana, que ocurre el 4 de

Radiacin Solar

julio y tiene un valor de 152 millones de Km. La distancia promedio que ms se usa es la de 149.6
millones de km.
Es posible estimar la distancia de la tierra al sol, para cualquier da del ao, a travs de diferentes
frmulas, de las ms comunes, se han escogido a las siguientes:
dm

360
1 0.0167 cos
n
365

dm

360
n
365

1 0.033 cos

En ambas los ngulos son medidos en grados a diferencia de:

dm

n
365

1 0.033 cos

cuyo ngulo se expresa en radianes.

En las tres frmulas dm representa a la distancia media de la tierra al sol, 149.6 millones de km, d es la
distancia real y n es el da del ao codificado de 1 (para enero 1) a 365 (para diciembre 31). Esta
codificacin tambin es conocida como da juliano. Para establecer el valor de n, existen ecuaciones
como las reportadas en el Cuadro 2.1 o a travs de valores directos como se muestran en el Cuadro 2.2.
Si se tuviera el inters de calcular la distancia de la tierra al sol para el da 4 de julio, n tendra el valor
del da 185 del ao, 181 + 4, y si se empleara a la ltima de las frmulas, el procedimiento sera:
dm

d
dm

1 0.033 cos
n
365

2 3.1416

185
365

1 0.033 cos

1 0.033 cos 3.1846

Radiacin Solar

Cuadro 2.1. Ecuaciones para convertir el da del mes en el

nmero del da del ao
Mes

n para el i-simo
da del mes

Enero
Febrero
Marzo
Abril
Mayo
Junio
Julio
Agosto
Septiembre
Octubre
Noviembre
Diciembre

dm

i
31 + i
59 + i
90 + i
120 + i
151 + i
181 + i
212 + i
243 + i
273 + i
304 + i
334 + i

0.9670

dm
0.9670 0.9834
d

y
d

dm
149.6

152.1
0.9834 0.9834

entonces:
d = 152.1 millones de km
Con lo anterior es posible inferir que si la Constante Solar se calcul con la distancia media, dm,
2
dm

entonces, la d influir en aumentar o disminuir su valor.

Durante mucho tiempo el valor de la constante solar fue estimado como 1.94 Langleys/min (1.94
cal/cm2/min), sobre la base de las mediciones del Instituto Smithsoniano. Posteriormente Johnson
(1954) estableci el valor de 2.0 cal/cm2/min (2 Langleys/min o 2 Ly/min), el mayor valor result por
la consideracin de la parte ultravioleta del espectro. Otro dato es el reportado por Sellers (1965), con
un valor de 1.974 Ly/min y el ms reciente por Driessen y Konijn (1992) de 1.9328 Ly/min.
Actualmente las unidades que se emplean para reportar a la radiacin son J/m2/s; para transformar
Ly/min en J/m2/s, se multiplica a los Ly/min por 700. As 1.9328 Ly/min multiplicado por 700 resulta
1353 J/m2/s.

Radiacin Solar

Cuadro 2.2. Calendario Juliano

Da

ENE

FEB

MAR

ABR

MAY

JUN

JUL

AGO

SEP

OCT

NOV

DIC

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31

32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
**
**
**

60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90

91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
**

121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151

152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
**

182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212

213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243

244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
**

274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304

305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
**

335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365

2.4.

Variacin de la Radiacin

La superficie terrestre recibe en promedio 300 cal/cm2/da y fcilmente podemos inferir a travs de la
constante solar que la energa recibida en la alta atmsfera es mucho mayor. Lo anterior nos indica la
existencia de prdidas de energa al atravesar a la atmsfera, las cuales pueden atribuirse a fenmenos
como absorcin, reflexin y dispersin.
La absorcin de la radiacin solar puede ser explicada por la Ley de Boliguer, que indica que la
intensidad calorfica de una radiacin que atraviesa un medio transparente decrece en progresin
geomtrica cuando la masa atravesada crece en progresin aritmtica, por ejemplo (De Fina y Ravelo
1973):
Masas atravesadas
Intensidad de Radiacin

1, 2,
3,
4,.........,n
1/3, 1/9, 1/27, 1/81,... ,1/3n

La reflexin vara considerablemente para las diferentes superficies, como se reporta a continuacin:

Radiacin Solar

Bosque
Ocano
Campo arado
Arena seca
Pradera
Nubes
Nieve fresca

10 %
12 %
15 %
20 %
25 %
75 %
80 %

Estos valores son promedios.

La dispersin es un fenmeno similar a la reflexin, diferencindose de sta, en que la radiacin
modifica sus caracteres al ser devuelta o desviada.
Ms especficamente, ya sobre la tierra, la cantidad diaria de calor recibida por centmetro cuadrado de
suelo horizontal, vara segn la poca del ao. Ello es debido a que:
a) Los rayos solares llegan a la superficie terrestre con distinta inclinacin, segn
la poca del ao.
b) La duracin del da es diferente segn las pocas del ao.
c) Los das son ms largos cuando los rayos solares son ms perpendiculares.
Estos tres fenmenos son debidos a que el eje de rotacin de la tierra forma un ngulo (de 23 27) con
la perpendicular, como se muestra en la Figura 2.2.
El Cenit, que es el da ms largo, para los trpicos, corresponde a los solsticios, de verano (21 de
junio) para el Trpico de Cncer y de invierno (21 de diciembre) para el Trpico de Capricornio. Los
equinoccios, de primavera (21 de marzo) y de otoo (23 de septiembre) son Cenit del Ecuador.
Grficamente podemos expresar lo anterior en las Figuras 2.3 y 2.4.
Al analizar tres puntos sobre la superficie terrestre en funcin de la energa solar (R) a travs del
tiempo, se obtendran las grficas de la Figura 2.5.
En el Cuadro 2.3 se reporta la duracin astronmica de la insolacin (N) que est en funcin de la
latitud.
La duracin astronmica de la insolacin (N) se considera como el tiempo en que habra luz solar en
un da promedio, suponiendo que la Tierra fuera plana. En cierta medida, como se analizar
posteriormente, la insolacin se relaciona con la radiacin solar y al graficar la N a la latitud de 10 y a
la latitud de 40 notaremos que a mayor latitud existe una poca del ao en la que se tiene un exceso de
energa (Figura 2.6).

10

Radiacin Solar

= 23o 27
Eje de
Rotacin

90o
Clima Polar

Clima Templado

66o 33 N
23o 27 N
0o
Clima Trrido
23o 27 S
Clima Templado
66o 33 S

Clima Trrido

Clima Polar
90o

Figura 2.2. Clima Solar

Solsticio de
Invierno

Solsticio de
Verano

Figura 2.3. Cenit de los Trpicos

Figura 2.4. Cenit de los Trpicos

11

Radiacin Solar

Figura 2.5. Variacin de la Radiacin Solar a diferentes Latitudes.

Cuadro 2.3. Duracin Astronmica de la Insolacin (N) Promedio Mensual, en Horas y Dcimos de Hora.
LAT
N

9.6

10.7

11.9

13.3

14.4

15.0

14.7

13.7

12.5

11.2

10.0

9.8

35
33
32
31
30

10.1
10.2
10.2
10.3
10.4

11.0
11.0
11.1
11.1
11.1

11.9
11.9
12.0
12.0
12.0

13.1
13.0
13.0
12.9
12.9

14.0
13.8
13.8
13.7
13.6

14.5
14.3
14.2
14.1
14.0

14.3
14.1
14.1
14.0
13.9

13.5
13.4
13.3
13.3
13.2

12.4
12.4
12.4
12.4
12.4

11.3
11.4
11.4
11.5
11.5

10.3
10.4
10.5
10.5
10.6

9.8
10.0
10.0
10.1
10.2

29
28
27
26
25

10.5
10.5
10.6
10.6
10.7

11.1
11.2
11.2
11.3
11.3

12.0
12.0
12.0
12.0
12.0

12.9
12.8
12.8
12.7
12.7

13.5
13.5
13.4
13.4
13.3

13.9
13.9
13.8
13.8
13.7

13.8
13.8
13.7
13.6
13.5

13.2
13.2
13.1
13.0
13.0

12.4
12.4
12.3
12.3
12.3

11.5
11.5
11.6
11.6
11.6

10.7
10.7
10.8
10.8
10.9

10.3
10.3
10.4
10.5
10.6

24
23
22
21
20

10.8
10.8
10.9
10.9
11.0

11.3
11.4
11.4
11.5
11.5

12.0
12.0
12.0
12.0
12.0

12.7
12.7
12.6
12.6
12.6

13.3
13.2
13.2
13.1
13.1

13.6
13.6
13.4
13.4
13.3

13.4
13.4
13.3
13.3
13.2

13.0
12.9
12.9
12.8
12.8

12.3
12.3
12.3
12.3
12.3

11.6
11.6
11.7
11.7
11.7

11.0
11.0
11.1
11.1
11.2

10.7
10.7
10.8
10.8
10.9

19
18
17
16
15
14

11.1
11.1
11.2
11.2
11.3
11.4

11.5
11.5
11.6
11.6
11.6
11.6

12.0
12.0
12.0
12.0
12.0
12.0

12.6
12.6
12.5
12.5
12.5
12.5

13.0
13.0
12.9
12.9
12.8
12.8

13.2
13.2
13.1
13.1
13.0
12.9

13.1
13.1
13.0
13.0
12.9
12.8

12.8
12.7
12.7
12.6
12.6
12.6

12.3
12.3
12.2
12.2
12.2
12.2

11.7
11.7
11.8
11.8
11.8
11.8

11.2
11.3
11.3
11.4
11.4
11.4

11.0
11.0
11.1
11.1
11.2
11.3

10

11.6

11.8

12.0

12.3

12.6

12.7

12.6

12.4

12.1

11.8

11.6

11.5

40

12

Radiacin Solar

Figura 2.6. Duracin Astronmica de la Insolacin (N) a travs del Ao para dos Latitudes.

Dicho exceso de energa afecta la produccin de cultivos, principalmente gramneas. Un ejemplo es el

trabajo de Chang (1981), en el que se analiz la produccin promedio de maz de tres aos, 19751977, para 98 pases. En la Figura 2.7 se muestra la grfica del rendimiento del maz contra la latitud y

Figura 2.7. Rendimientos de Maz como una funcin de la Latitud. (-----) Rendimiento Promedio para
cada 10 de latitud ( Chang, 1981)

13

Radiacin Solar

en el Cuadro 2.4 un resumen de los datos. En sus conclusiones, Chang indica que la diferencia entre
los rendimientos de las zonas templadas respecto a las tropicales es de 4:1 y que con mejoramiento
tecnolgico tal relacin puede disminuir a 2.5:1, como lo sugiere el examen de los rendimientos
mximos experimentales.
Cuadro 2.4. Rendimientos Promedio de Maz y Duracin del Da durante el Ciclo por Fajas Latitudinales
de 10 (Chang, 1981).
FAJA
LATITUDINAL

No. de
Pases

5N - 5S
5 - 15
15 - 25
25 - 35
35 - 45
45 - 55

10
27
23
15
12
11

Rendimiento
Promedio de
Maz (Kg./ha)
998
1018
1190
1788
3920
4714

Duracin del Da
Minutos
733
755
784
830
883
940

Horas
12.2
12.6
13.1
13.8
14.7
15.7

Sobre los rendimientos mximos experimentales para maz, Chang (1981) reporta para Asia; 6,500
kg./ha para Filipinas (10-20 N), 10,000 kg./ha para Pakistn (28 N) y 16,000 kg./ha para Corea del
Sur (36 N), indicando que resultados similares se han encontrado en Amrica. El rcord mundial de
24,100 kg./ha de Maz fue obtenido a una latitud de 42 N en Kalamazoo, Michigan.
2.5.

Fotoperiodo

Durante los equinoccios, 21 de marzo y 21 de septiembre, la duracin de los das y de las noches en
todas las latitudes es aproximadamente igual. Sin embargo durante el solsticio de verano, 21 de junio,
la duracin del da en el ecuador es de 12 horas, a 40 de latitud norte es de 15, a los 60 de latitud norte
es de 19 y en el polo norte de 24 horas.
Esta diferencia en la duracin del da y la noche, con respecto a la latitud, es un factor importante en la
distribucin de las plantas (Allard, 1948 citado por Wilsie, 1966). La respuesta de las plantas a la
duracin relativa del da y la noche es conocida como fotoperiodismo. Las plantas que se desarrollan y
reproducen normalmente cuando el fotoperodo es mayor a un mnimo critico se les denomina de da
largo 1y las que se desarrollan con un fotoperodo menor a un mximo critico reciben el nombre de da
corto.
Los primeros estudios que dieron origen a estos conceptos fueron los realizados por Garner y Allard en
1920 (citados por Wilsie, 1966), sobre tabaco y soya. Con la variedad de tabaco Maryland Mammoth,
se demostr, que cuando se cultivaba al aire libre en verano en Washington, D.C., solo creca
vegetativamente, mientras que, durante el invierno con das cortos floreca espectacularmente, la soya
mostr un comportamiento semejante. Ambas fueron consideradas como plantas de da corto.

14

Radiacin Solar

En investigaciones posteriores se propuso a un fotoperodo crtico, comprendido entre 12 y 14 horas,

permitiendo la clasificacin de las especies en dos grupos, de da largo y de da corto. Tambin se
encontr que otras especies no eran afectadas por el fotoperodo, denominndolas fotoneutras.
En el Cuadro 2.5 se muestra la respuesta al fotoperodo de algunas especies cultivadas. El fotoperodo
se ha usado como sinnimo de la duracin astronmica de la insolacin (N) o de la duracin del da
(DUR o DL) por diferentes autores. Sus valores pueden obtenerse del Cuadro 2.3 o calcularse a travs
de frmulas.
Cuadro 2.5. Respuesta al Fotoperodo durante la Floracin de algunas Especies Cultivadas (Doorenbos y
Kassam, 1979)
Cultivo
Alfalfa
Algodn
Arroz
Cacahuate
Caa de
azcar
Cebolla
Col
Ctricos
Frijol
Girasol
Guisante

Tipo de Fotoperodo
para inducir la
Floracin
Neutros
Das cortos/ Neutros
Das cortos/ Neutros
Neutros

Cultivo
Maz
Papa
Pimentero
Pltano
Remolacha
azucarera
Sanda
Sorgo
Soya
Tabaco
Tomate
Trigo

Das cortos/ Neutros

Das Largos/ Neutros
Das Largos
Neutros
Das cortos/ Neutros
Das cortos/ Neutros
Neutros

Tipo de Fotoperodo
para inducir la
Floracin
Neutros/Das cortos
Das Largos/ Neutros
Das Cortos/ Neutros
Neutros
Das Largos
Neutros
Das cortos/ Neutros
Das cortos/ Neutros
Das cortos/ Neutros
Neutros
Neutros /Das Largos

De las diferentes frmulas que existen mencionaremos cuatro, la propuesta por De Witt, et al. (1978);
por Hardy (1990), por Charles-Edwards, et al. (1986) y por Stuff y Dale (1973), por ser contrastantes.

Procedimiento de De Witt, et al.

15

Radiacin Solar

La duracin del da (DL) es una funcin del da del ao y la latitud del lugar (De Witt, et al. 1978). La
frmula que emplea es:
DL = 12 * (PI + 2 * asen (SSCC))/PI
con:
SSCC = SSEN / CCOS
SSEN = sen(LAT * RAD) * sen(DEC * RAD)
CCOS = cos(LAT * RAD) * cos(DEC * RAD)
DEC = - 23.45 cos (2*PI*(DAY + 10)/365)
donde:
RAD = factor de conversin (grados a radianes;
RAD = PI/180);
LAT = Latitud, en grados;
DEC = Declinacin solar, en grados;
DAY = Da Juliano en el Hemisferio Norte
Procedimiento de Hardy
La formula que emplea es:
DL = (24/) cos-1 (- tan * tan )
donde:
DL = Duracin del da en horas;
= Latitud en radianes y
= declinacin solar en radianes
de acuerdo con Iqbal (1983) el clculo de basado en series de Fourier es:
= 0.006918 - 0.399912 cos(da) + 0.070257 sen(da) - 0.006758 cos(2*da)
+ 0.000907 sen(2*da) - 0.002697 cos(3*da) + 0.00148 sen(3*da)
donde:
da = (2/365)(Dj - 1)
con:
Dj = da juliano.

Procedimiento de Charles-Edwards, et al.

16

Radiacin Solar

Este procedimiento es diferente a los anteriores porque calcula la duracin del da (h) en segundos, con
la frmula:
h = 4.36 x 104 + hD * sen [2 (t + z)/365]
donde:
t = Da Juliano (Dj);
z = 283 para el hemisferio norte;
hD = Variacin estacional, calculada como:
hD = exp(7.42 + 0.045 LAT) =

(7.42 +0.045 LAT)

para 10 < LAT < 55 .

Finalmente la duracin del da (DL), se puede transformar en horas con:

DL = h / 3600
Los Das Julianos que se mencionan en los mtodos de De Witt, como DAY, de Hardy, como Dj, y de
Charles-Edwards, como t, se obtienen numerando los das del ao de 1 (para el 1 de enero) a 365 (para
el 31 de diciembre), como se mostr en el Cuadro 2.2.
Mtodo de Stuff y Dale
Este mtodo emplea a la frmula:
DL = 12.14 + 3.34* tan(LA)* cos(0.0172*CD - 1.95)
donde:
DL = Duracin del da, en horas;
LA = Latitud, en Radianes
CD = Da climtico, que se calcula como:
CD = ENTERO(30.6* mes + da del mes - 91.3)

que se interpreta como el nmero entero ms cercano al valor obtenido dentro del parntesis. con:
mes = nmero del mes en estudio (enero = 1, febrero = 2,
marzo = 3, ................), con una restriccin;
S mes < 3, entonces, mes = mes + 12
Ejemplos de Clculo:
Supongamos que se desea calcular la duracin del da 3 de Septiembre y 12 de Febrero para el
observatorio de Chapingo, Mxico, con los cuatro mtodos indicados:
a) De Witt
Lo primero que se requiere en este procedimiento es el conocer el da juliano:
Para el 3 de Septiembre DAY = 246 con el cual se procede a calcular la declinacin solar,

17

Radiacin Solar

DEC = - 23.45 * cos[2*PI*(DAY + 10) / 365]

= - 23.45 * cos[2 * 3.1416(246 +10)/ 365]
= - 23.45 * cos(4.4068) = - 23.45 (- 0.3008)
= 7.0542
DEC * RAD = 7.0542* PI/180 = 7.0542*3.1416/180 = 0.1231
con la latitud:
LAT = 19 30 = 19 + (30 / 60) = 19.4833
y
LAT * RAD = 19.50 * PI/180 = 0.3403
con la latitud y la declinacin se procede a calcular CCOS, SSEN y SSCC:
CCOS = cos(LAT*RAD)*cos(DEC*RAD)
= cos(0.3403)*cos(0.1231) = 0.9427 * 0.9924
= 0.9355
SSEN = sen(LAT*RAD)*sen(DEC*RAD)
= sen(0.3403)*sen(0.1231) = 0.3338 * 0.1228
= 0.0410
SSCC = SSEN / CCOS = 0.0410 / 0.9355 = 0.0438
y finalmente:
DL = 12 * [PI + 2 * asen(SSCC)]/ PI
es conveniente aclarar que asen = arc sen, por lo que:
DL = 12 * [PI + 2 * arc sen(0.0438)]/ PI
= 12.3344
esto es:
DL = 12.33
Aplicando el mismo procedimiento para el da 12 de Febrero se tiene:
DAY = 43
DEC = - 14.3487,
DEC*RAD = - 0.2504,
LAT*RAD = 0.3403,
CCOS = 0.9132,
SSEN = - 0.0827,
SSCC = - 0.0906
DL = 11.31

b) Hardy
Para las fechas indicadas se tiene:

18

Radiacin Solar

FECHA
3 de Septiembre
Dj
246
da =(2/365)(Dj-1)
4.2175

0.006918
0.006918
-0.399912 cos(da)
0.189939
+0.070257 sen(da)
-0.061827
-0.006758 cos(2*da)
0.003709
+0.000907 sen(2*da)
0.000758
-0.002697 cos(3*da)
-0.002687
+0.00148 sen(3*da)
0.000127
=

12 de Febrero
43
0.7230
0.006918
-0.299865
0.046484
-0.000841
0.000899
0.001518
0.001223

0.1369

-0.2437

tan

0.1378

-0.2486

0.3403

0.3403

tan

0.3541

0.3541

1.5220

1.4827

cos-1 (-tan tan )

DL=(24/ )cos-1 (-tan tan )

12.37

11.33

c) Charles-Edwards
Los pasos necesarios son:
FECHA
t = Dj
hD = exp(7.42 + 0.045 LAT)
LAT
hD

3 de Septiembre
246

12 de Febrero
43

19.50

19.50

4010.81

4010.81

h = 4.36 x 104+ hD sen[2 (t + 283)/ 365]

h
44855.8147
DL = h / 3600
DL

41105. 0849

12.46

11.42

d) Stuff y Dale
La frmula a aplicar es:

19

Radiacin Solar

DL = 12.14 + 3.34* tan(LA)* cos(0.0172*CD - 1.95)

siendo CD la diferencia mayor con respecto a los otros mtodos, ya que:
CD = Da climtico, que se calcula como:
CD = ENTERO(30.6* mes + da del mes - 91.3)

Para septiembre 3, mes = 9 y da del mes = 3, por lo que:

CD = ENTERO (30.6* 9 + 3 - 91.3) = ENTERO (187.1) = 187
y

para febrero 12, como mes = 2 < 3, entonces mes = 2 + 12 = 14 y da del mes = 12, por lo que:
CD = ENTERO (30.6*14 + 12 - 91.3) = ENTERO (349.1) = 349

para ambas fechas la tan (LA) = 0.3541, por lo que, la frmula:

DL = 12.14 + 3.34* tan(LA)* cos(0.0172*CD - 1.95)
queda
DL = 12.14 + 3.34* 0.3541* cos(0.0172*CD - 1.95)
faltando por sustituir al da climtico, como sigue:
para septiembre 3:
DL = 12.14 + 3.34* 0.3541* cos(0.0172*187 - 1.95)
= 12.14 + 3.34* 0.3541* cos(1.2664) = 12.49
y para febrero 12:
DL = 12.14 + 3.34* 0.3541* cos(0.0172*349 - 1.95)
= 12.14 + 3.34* 0.3541* cos(4.0528) = 11.42
Una comparacin de los resultados de los cuatro mtodos se presenta en el Cuadro 2.6, en el cual
puede apreciarse que se obtienen diferentes valores, siendo el procedimiento de De Witt el que da los
ms bajos y el de Stuff y Dale los ms altos, la diferencia entre ellos es menor a 10 minutos.
Cuadro 2.6. Comparacin de cuatro mtodos para estimar la Duracin del da en Chapingo, Mxico para el
da 3 de Septiembre y el da 12 de Febrero.
Fecha
3 de Septiembre
12 de Febrero

De Witt
12.33
11.31

Hardy
12.37
11.33

Charles-Edwards
12.46
11.42

Stuff y Dale
12.49
11.42

Como podr notarse con los procedimientos descritos y aplicados es posible calcular DL, N o DUR,
para cualquier da del ao a diferencia de los valores reportados en tablas como en el Cuadro 2.3, que
corresponden al da 15 de cada mes.
Tambin a partir de la duracin del da se ha propuesto calcular la hora de salida y entrada del
sol, como sigue:

20

Radiacin Solar

DL
2
SALIDA =
DL
12
2
ENTRADA =
12

As para el da 3 de septiembre con una DL = 12.33 (Mtodo de De Witt) tendramos:

SALIDA = 12 (12.33/2) = 12 6.165 = 5.84 ( 5 hr 50 min) y
ENTRADA = 12 + (12.33/2) = 12 + 6.165 = 18.17 (18hr 10 min).
A veces se reportan en los textos tablas, con la Duracin Astronmica de la Insolacin promedio
mensual, como la reportada en el Cuadro 2.3; el trmino promedio mensual en realidad corresponde al
da 15 de cada mes, es decir, para el Observatorio Chapingo, tendramos los clculos y resultados que
se muestran en el Cuadro 2.7.
Cuadro 2.7. Duracin Astronmica de la Insolacin promedio mensual,con dos mtodos DeWitt y
Charles-Edwards para el Observatorio Chapingo, Estado de Mxico.
Mes
ENE 15
FEB 15
MAR 15
ABR 15
MAY 15
JUN 15
JUL 15
AGO 15
SEP 15
OCT 15
NOV 15
DIC 15

2.6.

Dj
15
46
74
105
135
166
196
227
258
288
319
349

DeWitt
10.94
11.36
11.86
12.44
12.92
13.17
13.07
12.67
12.11
11.55
11.06
10.83

Charles-Edwards
11.09
11.46
11.95
12.54
12.99
13.22
13.14
12.79
12.24
11.68
11.22
11.01

Estimacin de la Radiacin Solar

La radiacin solar recibida sobre la superficie terrestre consiste de dos partes: la radiacin solar directa
y la radiacin difusa en el cielo. La suma de estos dos componentes se les conoce como Radiacin
Global (Rg).
Toda la superficie terrestre recibe anualmente 860 x 10 18 Kcal de radiacin global. Este valor
representa un milln de veces el total de energa producida por la industria elctrica mundial (Schulze,
1970).
Griffiths (1978) indicaba que alrededor del mundo existan aproximadamente unas 600 estaciones
meteorolgicas en las que se medan sobre una base continua a la radiacin. Ubicadas principalmente

21

Radiacin Solar

en Estados Unidos, algunos pases de Europa y en la Repblica de Sudfrica. Sin embargo, esta
situacin en la actualidad ha cambiado y ya se dispone de un gran nmero de estaciones que realizan
esta medicin, por ejemplo para Mxico se cuenta con 94 estaciones automticas, que miden de hecho
cada minuto a la radiacin solar (CNA, 2005).
En el pasado por la dificultad de obtener datos de radiacin medidos en forma directa, diferentes
autores correlacionaban a la radiacin global recibida en un lugar determinado con la duracin de la
insolacin o con el grado de nubosidad del cielo.
En los siguientes prrafos se describen los mtodos empleados en Mxico para estimar a la radiacin.
a) Mtodo de ngstrom y Prescott
Uno de los primeros mtodos en los que se busc la relacin Radiacin Global - Insolacin fue el de
Angstrm dado a conocer en 1924, quien propuso la frmula siguiente:

Rg
n
a b

Rgo
N
donde:
Rg, es la radiacin global,
Rgo, es la radiacin global de un da muy despejado,
n, es el nmero real de horas de insolacin,
N, es la duracin astronmica de la insolacin, y
a y b, coeficientes de regresin, originalmente, a + b = 1.
Posteriormente, Prescott (1940), propuso una modificacin a la frmula de Angstrm, en la que se
usan los valores de Angot (RA) o radiacin terica, es decir, la que recibira la superficie terrestre en
ausencia de atmsfera. Expresndose la frmula modificada de la siguiente manera:

Rg
n
a b

RA
N
adems los coeficientes a y b al sumarse podran tener un valor diferente de 1. Esta frmula ha sido la
ms ampliamente utilizada.
En el Cuadro 2.8 se reportan los valores de Angot (a partir de Torres, 1983). Al igual que para la
duracin del da, existen frmulas para su clculo, como la propuesta por FAO (1992):
RA = 37.586*d*[Hs * sen() * sen() + cos() * cos() * sen(Hs)]
donde:
Cuadro 2.8. Valores Mensuales de la Radiacin Solar Terica, RA, en cal/cm2/da (Torres, 1983)
LAT

ENE

FEB

MAR

ABR

MAY

JUN

22

JUL

AGO

SEP

OCT

NOV

DIC

Radiacin Solar
(N)
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33

697
685
672
660
647
635
622
609
596
582
569
555
542
527
514
500
486
472
458

772
763
753
743
734
723
710
703
692
681
670
659
647
636
623
611
599
587
574

846
841
836
830
824
818
811
805
798
791
783
775
768
756
751
742
734
725
715

896
896
896
895
894
893
892
891
889
887
884
882
879
871
872
869
865
861
856

909
913
917
920
924
927
930
932
935
937
939
941
942
943
944
945
945
946
946

907
913
919
924
930
934
940
944
949
953
960
961
964
971
971
973
976
978
981

909
914
919
923
928
932
936
940
943
947
950
952
955
957
959
961
963
964
965

904
905
907
908
909
909
910
910
910
909
909
908
906
905
903
901
899
897
894

868
865
861
858
853
849
844
839
834
829
823
817
811
807
798
791
784
776
769

797
789
781
772
763
754
745
736
726
716
706
696
685
678
664
652
641
630
618

715
704
692
680
668
656
643
634
618
605
592
579
566
553
540
525
512
498
484

673
660
647
634
621
608
594
581
567
553
539
525
511
495
482
468
454
439
424

RA = Radiacin Extraterrestre, MJ/m2*da;

= Latitud. radianes;
= Declinacin Solar, radianes;
calculada por:
2

0.4093 sen
Dj 1.405
360

con Dj = Da Juliano;
Hs = ngulo solar a la hora de salida, radianes
Hs = arc cos[- tan()* tan()]
d = Distancia relativa de la tierra al sol,
calculada por:

d 1 0.033 cos
Dj
365

con Dj = Da Juliano

Es conveniente mencionar que la distancia relativa d representa a (dm/d)2 referida con anterioridad al
tratar a la constante solar.
Por considerarlo de inters, en el Cuadro 2.9 se presenta una tabla para la conversin de las unidades
Cuadro 2.9. Tabla de Conversiones de Unidades de Energa.
Para
De

convertir
A

23

Multiplique
por

Radiacin Solar

Cal/cm2/da
Cal/cm2/da
Cal/cm2/da
mm/da
mm/da
mm/da
MJ/m2.da
MJ/m2.da
MJ/m2.da
W/m2
W/m2
W/m2

mm/da
MJ/m2.da
W/m2
Cal/cm2/da
MJ/m2.da
W/m2
Cal/cm2/da
mm/da
W/m2
Cal/cm2/da
mm/da
MJ/m2.da

0.0171
0.0419
0.4846
58.508
2.4497
28.352
23.884
0.4082
11.574
2.0636
0.0353
0.0864

de energa solar ms empleadas.

Adems, se han obtenido modelos relacionando la latitud con los valores de RA, como se muestran en
el Cuadro 2.10, con los cuales se facilita la interpolacin.
Cuadro 2.10. Modelos Mensuales par Estimar RA, en cal/cm2/da, en funcin de la Latitud en Grados
(Ortiz, 2007)
MES
MODELO
R2
ENERO
RA = 872.72 - 10.511 (L) - 0.0619 (L)2
0.9999
FEBRERO
RA = 899.602 - 6.59394 (L) - 0.0956439 (L)2
0.9999
MARZO
RA = 911.156 - 1.74545 (L) - 0.119318 (L)2
0.9999
ABRIL
RA = 878.348 + 4.11212 (L) - 0.132576 (L)2
0.9994
MAYO
RA = 827.229 + 8.13636 (L) - 0.125947 (L)2
0.9992
2
JUNIO
RA = 795.826 + 9.75455 (L) - 0.116477 (L)
0.9998
JULIO
RA = 806.453 + 9.00909 (L) - 0.123106 (L)2
0.9997
AGOSTO
RA = 852.912 + 5.50303 (L) - 0.130682 (L)2
0.9976
SEPTIEMBRE
RA = 892.229 + 0.136364 (L) - 0.125947 (L)2
0.9999
OCTUBRE
RA = 898.73 - 5.78485 (L) - 0.09375 (L)2
0.9999
NOVIEMBRE
RA = 874.644 - 10.0242 (L) - 0.0625 (L)2
0.9999
2
DICVIEMBRE
RA = 864.364 - 11.7909 (L) 0.0454545 (L)
0.9999

El problema en el uso de la frmula de Angstrm, modificada por Prescott radicaba en la definicin de

los valores de a y b. Para Europa, Penman (1948-1963) propuso los valores de 0.18 y 0.55,
respectivamente; mientras que Glover y MacCulloch (1958), a partir de sus estudios en las zonas
montaosas de frica, consideran que:
a = 0.29 cos y b = 0.52
donde es la latitud del lugar, siendo uno de los procedimientos ms empleados y difundidos. A pesar
de ello es conveniente indicar que existen en la literatura una gran diversidad de valores que se les
puede asignar a dichos coeficientes.
b) Mtodo de Hargreaves y Samani

24

Radiacin Solar

El mtodo de Hargreaves y Samani fue propuesto por sus autores en 1982, el cul calcula a la
radiacin con la siguiente frmula:
Rg = KT*RA*D1/2
donde:
Rg = Radiacin global, cal/cm2/da;
RA = Radiacin Astronmica, cal/cm2/da;
D = Oscilacin Trmica = Tmax Tmin;
KT = Coeficiente de calibracin, con:
KT = 0.075 (S/D)1/2 y
S = Porcentaje de Insolacin = (n/N)x 100
c) Mtodo de Allen
Mtodo propuesto en 1997, se indica que se parti del mtodo anterior y la frmula que emplea es:

Rg
RA

Kr TM Tm

0.5

donde:
Rg = Radiacin global, cal/cm2/da;
RA= Radiacin Astronmica, cal/cm2/da;
TM = Temperatura media mxima;
Tm = Temperatura media mnima;
Kr = Coeficiente de calibracin, con:
Kr = Kra (P/Po)0.5
Para Kra se han propuesto dos valores: 0.20 para las regiones costeras y 0.17 para las regiones del
interior, y; P = Presin atmosfrica del sitio, en kPa, la cul se puede calcular a partir de:
P = 101 0.0116 ELEV + 5.44 x 10-7(ELEV)2
donde:
ELEV = Elevacin sobre el nivel del mar, en m, y
2.7.

Po = Presin a nivel del mar = 101.3 kPa.

Estimacin de la Radiacin Solar en Mxico

En Mxico se reportan para el perodo 1981 2000 datos del Total de Horas de Insolacin Mensuales
para 73 Observatorios, con lo que es posible estimar a la radiacin global con los primeros dos
mtodos.

25

Radiacin Solar

Mtodo de Angstrm-Prescott
En el Cuadro 2.11 se reportan los datos necesarios, para cuatro meses y los resultados de los clculos
para obtener la radiacin global para el Observatorio de Tacubaya, D.F. con el mtodo de Angstrm y
Prescott.
Cuadro 2.11. Clculo de la Radiacin Global con la Frmula de ngstrom - Prescott para Cuatro Meses con
Datos del perodo 1981 2000 para el Observatorio de Tacubaya,D.F.
DATO
Total de Horas de Insolacin
No. de Das del Mes
n
N (DL o DUR)
n/N
Rg/RA
RA (cal/ cm2/da)
Rg (cal/ cm2/da)

ENERO
240
31
7.74
11.09
0.70
0.64
645
411

ABRIL
232
30
7.73
12.53
0.62
0.59
908
540

JULIO
176
31
5.68
13.14
0.43
0.50
935
466

OCTUBRE
194
31
6.26
11.68
0.54
0.55
751
415

Dado que el Observatorio de Tacubaya, D.F, se localiza a una latitud de 19 24 03, al utilizar la
frmula:

Rg
n
a b

RA
N
lo primero que se puede calcular es el valor de la constante a, con la frmula de Glover y MacCulloch
(1958), como sigue:
a = 0.29 cos = 0.29 cos(19 24 03) = 0.29*cos (19.40361) = 0.2735
por lo que para Tacubaya, especficamente se tendra como frmula a:

Rg
n
0.27 0.52

RA
N
Para obtener el valor de n, primero debemos entender que representa el valor promedio de un da de
las horas de insolacin, en cada mes. Razn por la cual en el Cuadro 2.11, se presentan en primer
trmino dos hileras que corresponden al Total de Horas de Insolacin mensuales como se reportan en
los Observatorios Normales (SMN, 2007) y al Nmero de das del mes, generndose n a travs del
cociente de las dos anteriores, es decir:
n = (Total de Horas de Insolacin) / (No. de das del mes)
As para Enero, se tiene:

26

Radiacin Solar

n = 240 / 31 = 7.74
El valor de N o DL, puede obtenerse por interpolacin de los valores reportados en una tabla, como los
del Cuadro 2.3, para la latitud del lugar o bien, se puede aplicar uno de los procedimientos descritos
(De Witt, Hardy, Charles-Edwards o Stuff y Dale), en el ejemplo se reporta a N empleando el mtodo
de Charles - Edwards. Para generar el cociente n/N.
En la frmula especfica de Tacubaya se sustituyen los valores de n/N y se obtienen los valores de
Rg/RA.
La hilera de RA se obtiene a partir de interpolacin de la tabla reportada en el Cuadro 2.8 o por la
frmula de FAO (1992) o de los modelos del Cuadro 2.10. Que multiplicados por (Rg/R A) generan los
valores de Rg.
Mtodo de Angstrm Prescott modificado
Se utiliza la misma frmula y los mismos coeficientes especficos para el Observatorio de inters, en
nuestro caso Tacubaya. La modificacin consiste en la forma como se estima a la insolacin media
mensual (n), la cul se basa en la frmula propuesta por Vargas y Tejeda (1996) a partir de la
nubosidad, como sigue:

x 0.5 x2

n N15 0.31 0.48 3
x

2
3
1
donde:
n = nubosidad media mensual, en horas;
N15 = Duracin Astronmica de la Insolacin para el da 15
de cada mes, en horas;
x1 = No. de das del mes de nublados a cerrados;
x2 = No. de das del mes medio nublados;
x3 = No. de das del mes despejados.
Los autores del modelo indican que trabajaron con 672 datos de Mxico, que obtuvieron una
correlacin de 0.90 y que es vlido para las latitudes de 14o < LAT < 33o.
Dado que el modelo tambin puede expresarse como:
x 0.5 x2
n

0.31 0.48 3
N15
x1 x2 x3

Se decidi sustituir el valor de n/N en la frmula de Angstrm por n/N15 de Vargas y Tejeda.
En el Cuadro 2.12 se reportan los datos, para los mismos cuatro meses y los resultados de los clculos
para obtener la radiacin global para el Observatorio de Tacubaya, con el mtodo de Angstrm y
Prescott modificado con el modelo de Vargas y Tejeda.

27

Radiacin Solar

Cuadro 2.12. Clculo de la Radiacin Global con el Mtodo de ngstrom - Prescott modificado con el modelo
de Vargas y Tejeda para Cuatro Meses con Datos de las Normales de 1980 para el Observatorio
de Chapingo, Mxico.
DATO
x3: No. de das despejados
x2. No. de das medio nublados
x1: No. de das de nublados a
cerrados
n/N
Rg/RA
RA (cal/ cm2/da)
Rg (cal/ cm2/da)

ENERO
11.2
15.6

ABRIL
4.6
18.4

JULIO
0.8
8.5

OCTUBRE
3.3
12.2

4.2
0.60
0.59
645
379

7.0
0.53
0.55
908
499

21.7
0.39
0.48
935
444

15.5
0.46
0.51
751
383

Nuevamente se utiliza a la frmula especfica para Tacubaya:

Rg
n
0.27 0.52

RA
N
Calculando el valor de n/N con los datos de nubosidad, que comnmente se reportan al final de los
observatorios, bajo la denominacin de Fenmenos Especiales. Una vez obtenido dicho cociente se
aplica el mismo procedimiento que en el mtodo anterior.
Es evidente que al cambiar de mtodo tambin cambien los valores.
Mtodo de Hargreaves y Samani
En el Cuadro 2.13 se reportan los datos para cuatro meses y los resultados de los clculos para obtener
la radiacin global para el Observatorio de Tacubaya, con el mtodo de Hargreaves y Samani.
Como se indic la radiacin global con el mtodo de Hargreaves y Samani se calcula con la frmula:

Cuadro 2.13. Clculo de la Radiacin Global con el Mtodo de Hargreaves y Samani para Cuatro Meses con
Datos del perodo 1981 2000 para el Observatorio de Tacubaya,D.F.
DATO
D: Oscilacin Trmica
S: Porcentaje de insolacin
KT: Coeficiente de Calibracin
RA (cal/ cm2/da)
Rg (cal/ cm2/da)

ENERO
14.3
70
0.1657
645
404

ABRIL
14.5
62
0.1547
908
535

28

JULIO
11.3
43
0.1467
935
461

OCTUBRE
11.7
54
0.1605
751
412

Radiacin Solar

Rg = KT *RA *D1/2
La oscilacin trmica (D), se puede obtener directamente de la informacin de la estacin o calcularse
a travs de la diferencia entre Tmax Tmin , que son los valores normales de las temperaturas mxima y
mnima mensuales respectivamente.
El porcentaje de insolacin (S), es el cociente n/N multiplicado por 100, por lo que resulta necesario
disponer de datos de insolacin real (n) y calcular la insolacin astronmica (N), con alguno de los
mtodos ya indicados.
El coeficiente de calibracin KT se calcula con:
KT = 0.075 (S/D)1/2
La RA como se mencion con anterioridad se obtiene a partir de interpolacin con los valores de la
tabla reportada en el Cuadro 2.8 o por la frmula de la FAO (1992) o de los modelos del Cuadro 2.10.
Finalmente se multiplican los valores de la raz cuadrada de D, la RA y del coeficiente de calibracin
KT para obtener los valores de Rg.
Mtodo de Hargreaves y Samani modificado con la Nubosidad
La diferencia con el mtodo anterior radica en el clculo de S, empleando la frmula:
S = (n/N15)*100
Donde n/N15 se genera aplicando el modelo de Vargas y Tejeda (1996). En el Cuadro 2.14 se muestran
los resultados con tal modificacin.
Mtodo de Hargreaves y Samani Simplificado
Al reconocer que la disponibilidad de los datos meteorolgicos en nuestro pas es una limitante, y en
particular al no contar con datos de insolacin y de nubosidad, para un buen nmero de estaciones, un
procedimiento que resulta til es el propuesto por Ortiz y Pjaro (1993), que se basa en el trabajo de
Hargreaves y Samani (1982), es decir, parte tambin de la frmula:
Cuadro 2.14. Clculo de la Radiacin Global con el Mtodo de Hargreaves y Samanimodificado con el modelo
de Vargas y Tejeda, para Cuatro Meses con Datos del perodo 1981 - 2000 para el Observatorio
de Tacubaya, D. F.
DATO
D: Oscilacin Trmica
S*: Porcentaje de insolacin
KT: Coeficiente de Calibracin
RA (cal/ cm2/da)
Rg (cal/ cm2/da)

ENERO
14.3
60
0.1542
645
376

ABRIL
14.5
53
0.1435
908
496

JULIO
11.3
39
0.1390
935
437

* Calculado como S = (n/N15)*100 y (n/N15) = 0.31 + 0.48 (x3 + 0.5 x2)/(x1 + x2 + x3)

29

OCTUBRE
11.7
46
0.1480
751
380

Radiacin Solar

Rg = KT RA D1/2
Desde un inicio se consider que el problema para su aplicacin en nuestro medio, era el valor de K T,
una serie de modelos fueron estimados a partir de la oscilacin trmica y la precipitacin, como
sustituto de la nubosidad, en el mbito mensual como se muestran en el Cuadro 2.15, estos modelos se
generaron con los datos de 40 Observatorios y de las normales del perodo 1941 1970 (SMN, 1976).
Adems, Ortiz y Pjaro (1993), en su trabajo indicaron los intervalos para la oscilacin trmica y la
precipitacin, donde los modelos son vlidos, los cuales se reportan en el Cuadro 2.16.
Recomendando el procedimiento siguiente:
1. Comparar la informacin mensual de oscilacin trmica y precipitacin con el Cuadro 2.16, s
se ubican dentro de los intervalos sealados, aplique los modelos del Cuadro 2.15, para
estimar a KT.
Cuadro 2.15. Modelos para Estimar Valores Mensuales de K T a partir de la Oscilacin Trmica (D) y la
Precipitacin (P).
MES
ENERO
FEBRERO
MARZO
ABRIL
MAYO
JUNIO
JULIO
AGOSTO
SEPTIEMBRE
OCTUBRE
NOVIEMBRE
DICIEMBRE

MODELO
KT = 0.349276 - 0.020433 (D) + 0.000511 (D2) - 0.001063 (P)
KT = 0.34827 - 0.018933 (D) + 0.000441 (D2) - 0.001361 (P)
KT = 0.335301 - 0.018357 (D) + 0.000438 (D2) - 0.001174 (P)
KT = 0.329488 - 0.020221 (D) + 0.000530 (D2) - 0.000457 (P)
KT = 0.367345 - 0.026878 (D) + 0.000796 (D2) - 0.000166 (P)
KT = 0.400671 - 0.032887 (D) + 0.001049 (D2) - 0.000102 (P)
KT = 0.376354 - 0.029467 (D) + 0.000917 (D2) - 0.000072 (P)
KT = 0.390603 - 0.030780 (D) + 0.000953 (D2) - 0.000064 (P)
KT = 0.383426 - 0.032185 (D) + 0.001068 (D2) - 0.000071 (P)
KT = 0.389938 - 0.030368 (D) + 0.000957 (D2) - 0.000121 (P)
KT = 0.381204 - 0.025179 (D) + 0.000686 (D2) - 0.000573 (P)
KT = 0.349787 - 0.021413 (D) + 0.000548 (D2) - 0.000893 (P)

R2
0.8133
0.8604
0.8665
0.7756
0.8006
0.8767
0.8706
0.8518
0.8106
0.7921
0.7902
0.8051

Cuadro 2.16. Intervalos de la Oscilacin Trmica (D) y de la Precipitacin (P) Mensuales para
los Modelos de KT.
Mes
ENERO
FEBRERO
MARZO
ABRIL
MAYO
JUNIO
JULIO
AGOSTO

Intervalo para D
5.9 - 21.6
6.1 - 22.4
5.7 - 22.9
5.5 - 21.7
5.5 - 18.9
4.4 - 16.9
4.6 - 16.1
5.0 - 15.5

30

Intervalo para P
0.5 - 43.2
0.9 - 46.1
0.9 - 54.7
0.5 - 103.2
0.2 - 293.0
1.9 - 325.8
38.1 - 421.5
43.7 - 353.5

Radiacin Solar

SEPTIEMBRE
OCTUBRE
NOVIEMBRE
DICIEMBRE

5.0 - 15.4
5.0 - 17.1
5.3 - 19.9
5.7 - 20.8

27.2 - 488.5
14.6 - 402.3
5.1 - 75.5
0.0 - 46.4

2. Obtener los valores de RA para la estacin de inters, por cualquiera de los procedimientos
previamente descritos.
3. Generar los valores de Rg, a partir de la frmula de Hargreaves y Samani.
Al aplicar el procedimiento a los datos del Observatorio de Tacubaya, se genera el Cuadro 2.17.
Cuadro 2.17. Clculo de la Radiacin Global con el Mtodo de Hargreaves y Samani simplificado por Ortiz y
Pjaro para Cuatro Meses con Datos del perodo 1981 2000 para el Observatorio de
Tacubaya,D.F.
DATO
D: Oscilacin Trmica
P: Precipitacin
KT: Coeficiente de Calibracin
RA (cal/ cm2/da)
Rg (cal/ cm2/da)

ENERO
14.3
7.6
0.1533
645
374

ABRIL
14.5
22.5
0.1374
908
475

JULIO
11.3
189.5
0.1468
935
461

OCTUBRE
11.7
72.4
0.1569
751
403

Mtodo de Allen
Otro procedimiento que en la actualidad tiene grandes posibilidades de aplicacin en Mxico, es el
mtodo de Allen (1997), que como se mencion emplea la frmula:

Rg
0.5
Kr TM Tm
RA
De hecho la TM Tm es equivalente a la oscilacin trmica D. En realidad este mtodo constituye otra
forma modificada del procedimiento de Hargreaves y Samani, para el clculo de Rg.
Al aplicar el procedimiento al Observatorio de Tacubaya, por ser una regin del interior su Kra = 0.17
y como su ELEV = 2309 m, genera una P = 77.416 kPa, por lo que tiene una
Kr = 0.17 (77.416/101.3)0.5 = 0.1311
En el Cuadro 2.18 se muestran los resultados. Que como era de esperarse varan con respecto al
mtodo anterior.
Cuadro 2.18. Clculo de la Radiacin Global con el Mtodo de Allen para Cuatro Meses con Datos del perodo
1981 2000 para el Observatorio de Tacubaya,D.F.
DATO

ENERO

ABRIL

31

JULIO

OCTUBRE

Radiacin Solar

D = (TM Tm)
Rg/RA
RA (cal/ cm2/da)
Rg (cal/ cm2/da)

14.3
0.50
645
320

14.5
0.50
908
454

11.3
0.44
935
412

11.7
0.45
751
337

Rg a partir de Datos de Estaciones Meteorolgicas

Datos de estaciones meteorolgicas del pas, se pueden obtener a partir de las Normales del
perodo 1971 2000, recurriendo al portal del Servicio Meteorolgico Nacional, smn.cna.mx
(SMN, 2007) o a travs del sistema ERIC III (IMTA, 2006).
La diferencia entre ambas fuentes es que con el sistema ERIC III es posible establecer datos de
nubosidad. Por ejemplo, para el observatorio Chapingo, Estado de Mxico, la estadstica histrica
para el % de das nublados arroja los datos del Cuadro 2.19.
Cuadro 2.19. Estadstica Histrica para el % de Das Nublados generada con el Sistema
ERIC III para el Observatorio de Chapingo, Mxico.
Dato

ENE

FEB

MAR

ABR

MAY

JUN

JUL

AGO

SEP

OCT

NOV

DIC

Das
Aos
MNub
Desv
Nubl
Desv

1457
47
47.3
25.6
9.2
11.2

1329
47
50.7
23.9
6.1
7.5

1457
47
52.5
26.9
6.5
8.7

1410
47
61.1
28.4
10.4
15.5

1488
48
61.0
29.6
14.2
17.9

1440
48
56.9
26.6
29.5
24.5

1457
47
59.1
27.1
34.7
25.1

1488
48
61.4
29.4
30.3
27.4

1440
48
58.2
25.2
32.1
23.6

1488
48
55.7
24.9
18.1
16.1

1440
48
56.3
24.0
9.1
11.3

1488
48
54.2
23.6
8.7
10.5

Con los datos del Cuadro 2.19 se genera la informacin del Cuadro 2.20, relativa al nmero de
das Nublados, Medio Nublados y Despejados para Chapingo. Para ello, se consider al nmero
total de das del mes y los das Despejados se calcularon por diferencia.
Con los datos del Cuadro 2.20 se puede aplicar el procedimiento de ngstrom Prescott
modificado.

Cuadro 2.20. Nmero de Das Nublados, Medio Nublados y Despejados para el

Observatorio de Chapingo, Mxico.
Dato

Das
Nublados
Das
MNublados
Das
Despejados

ENE

FEB

MAR

ABR

MAY

JUN

JUL

AGO

SEP

OCT

NOV

DIC

2.9

1.7

2.0

3.1

4.4

8.9

10.8

9.4

9.6

5.6

2.7

2.7

14.7

14.2

16.3

18.3

18.9

17.1

18.3

19.0

17.5

17.3

16.9

16.9

13.5

12.1

12.7

8.6

7.7

4.1

1.9

2.6

2.9

8.1

10.4

10.4

32

Radiacin Solar

Adems, con las dos fuentes se pueden obtener datos de Temperaturas Mxima y Mnima y
Precipitacin, para con ellos aplicar los procedimientos de Hargreaves y Samani y el mtodo de
Allen.
2.8.

Distribucin de la Radiacin Solar en Mxico

Se han realizado varios intentos para elaborar mapas mundiales de la distribucin de la radiacin solar.
Uno de los ltimos es el elaborado por Landsberg (1961) quien utiliz ms de 300 estaciones, su mapa
es muy generalizado y no considera variaciones mensuales. En este trabajo para Mxico se reportan
tres isolneas, la de 180 Kcal/cm2/ao al Noroeste del pas, la de 160 Kcal/cm 2/ao que pasa por los
limites de las Baja Californias, sigue por el Ocano Pacfico y se introduce al territorio al sur de
Nayarit, baja a una latitud inferior a 20 y otra isolnea tambin de 160 Kcal/cm 2/ao que rodea al
estado de Chiapas.
Jen-hu Chang (l968) indica que los mapas de radiacin son guas muy tiles para evaluar las
potencialidades agrcolas de diferentes regiones.
De Vries (1963) es de la opinin que, desde un punto de vista biolgico, seria preferible utilizar a la
radiacin solar como un factor climtico fundamental en lugar de la temperatura del aire en la
clasificacin de microambientes.
En la Figura 2.8 se muestra un mapa de la radiacin global media anual para la Repblica
Mexicana elaborada por estudiantes de la Universidad Autnoma Chapingo en el ao de 1982.
2.9.

Utilizacin de la Radiacin por los Cultivos

En el campo, la eficiencia de la utilizacin de la radiacin solar por los cultivos es baja y esto es debido
a dos razones: (1) la superficie del suelo no est completamente cubierta por el cultivo, desgastndose
una gran proporcin de radiacin, y (2) existen deficiencias variables en agua y nutrientes minerales,
daos de plagas y enfermedades y temperaturas desfavorables. Otro factor que contribuye a la aparente
baja eficiencia es la exclusin de informacin sobre el material radical, el cual puede ser ms del 30%
del total de materia seca en datos de cosecha (Schuyrman y Makkink, l955).

33

Radiacin Solar

Figura 2.8. Radiacin Global Media Anual para la Repblica Mexicana

La eficiencia de la utilizacin de la radiacin puede ser calculada por comparacin entre el valor
calrico de la materia orgnica producida por unidad de rea cultivada y la radiacin incidente en la
misma rea durante el mismo perodo. Algunos datos de la eficiencia de la utilizacin de la radiacin
se presentan en el Cuadro 2.21. Estos datos indican que ordinariamente un cultivo convierte menos del
1% de la radiacin solar. Spoehr (1956) estim que bajo condiciones ptimas, el maz puede convertir
cerca del 1.5% de la radiacin incidente en materia orgnica, durante un perodo de crecimiento de
cuatro meses.
En Hawaii, para Caa de Azcar se ha generado un modelo que relaciona a las Toneladas de Azcar
por acre (TSA) con la Radiacin global (Rg), como sigue:
TSA = - 17.77 + 0.055 Rg
Ese modelo tiene un coeficiente de correlacin de 0.80, que de acuerdo a las especificaciones para
considerar un buen modelo con datos de nuestro pas no sera aceptable (R2 > 0.75). Sin embargo, dada
la alta confiabilidad de la informacin generada en ese lugar, se usa para predecir los rendimientos de
azcar.

34

Radiacin Solar

Cuadro 2.21. Produccin de Materia Seca y Eficiencia de la Utilizacin de la Radiacin para varios Cultivos
Agrcolas.
Perodo Vegetativo

Produccin
(ton/ha)

Eficiencia
(%)

Abril - Agosto
Noviembre - Marzo
Mayo - Octubre
Mayo - Octubre
Mayo - Octubre
Junio - Septiembre
22 meses

9.60
10.45
16.00
16.00
6.86
15.52
129.48

0.50
0.52
0.90
0.90
0.39
1.05
1.43

CULTIVO

Papa
Trigo de Invierno
Remolacha Azucarera
Remolacha Forrajera
Zanahorias
Maz
Caa de Azcar

Otro uso que tiene el dato de Rg es el ndice de Productividad del IRRI (International Rice Research
Institute), que adopt el modelo de Yoshida y Parao (citado por Murata, 1975), para predecir los
rendimientos mximos de arroz. Dicho modelo es:
Y = Rg (278 - 7.07 t) x 0.86 x 18.1 x 10-5
donde:
Y = Rendimiento Mximo Experimental, Ton/ha;
Rg = Radiacin global, cal/cm2/da, durante los 25 das antes de la floracin;
t = Temperatura en C, durante los 25 das antes de la floracin;
0.86 = Porcentaje promedio de llenado de grano;
18.1 = Peso promedio de 1000 granos;
10-5 = Factor de correccin.
que es considerado como una modificacin del modelo propuesto por Murata (1964):
Y = Rg {1.20 - 0.021(t - 21.5)2}
donde:
Y = rendimiento de grano en Kg./1000 m2 ;
Rg = Radiacin solar (cal/cm2 /da) en Agosto y Septiembre y;
t = Temperatura media diaria (C), en Agosto y Septiembre.
Este tipo de modelos para caa de azcar y para arroz pueden ser de inters en Mxico y sera muy
recomendable su investigacin con datos de nuestro medio.
Dentro de la metodologa de Zonas Agroecolgicas propuesta por la FAO (1978) para la estimacin
de rendimientos de cultivos, se menciona la existencia de la Radiacin Fotosintticamente Activa
(RFA), cuyos valores han sido establecidos por De Witt para diferentes latitudes y das perfectamente
claros (Ac), que se reportan en el Cuadro 2.22. La RFA en un da totalmente cerrado es el 20% de

35

Radiacin Solar

Cuadro 2.22. Radiacin Fotosintticamente Activa en das Muy Claros (Ac), en cal/cm 2/da y Tasa de
Fotosntesis Bruta diaria en Cultivos cerrados en das Muy Claros (bc) y en das Cerrados
(bo) en Kg./ha/da, para una Fotosntesis Mxima (fm) de 20 Kg. de CH2O/ha/h (De Witt).
LAT
N
0 Ac
bc
bo

ENE

FEB

MAR

ABR

MAY

JUN

JUL

AGO

SEP

OCT

NOV

DIC

343
413
219

360
424
226

369
429
230

364
426
228

349
417
221

337
410
216

342
413
218

357
422
225

368
429
230

365
427
228

349
418
222

337
410
216

10 Ac
bc
bo

299
376
197

332
401
212

359
422
225

375
437
234

377
440
236

374
440
235

375
440
236

377
439
235

369
431
230

345
411
218

311
385
203

291
370
193

20 Ac
bc
bo

249
334
170

293
371
193

337
407
215

375
439
235

394
460
246

400
468
250

399
465
249

386
451
242

357
425
226

313
387
203

264
348
178

238
325
164

30 Ac
bc
bo

191
281
137

245
333
168

303
385
200

363
437
232

400
471
251

417
489
261

411
483
258

384
456
243

333
412
216

270
356
182

210
299
148

179
269
130

40 Ac
bc
bo

131
218
99

190
283
137

260
353
178

339
327
223

396
480
253

422
506
268

413
497
263

369
455
239

298
390
200

220
314
155

151
241
112

118
204
91

Ac y la RFA es el 50% de la Radiacin Global (Rg). Entonces, la fraccin del da-tiempo de desarrollo
cuando el cielo est cerrado (F), es:

Ac 0.5Rg
0.8 Ac

La produccin mxima de biomasa bruta por cultivo caracterizado por un ndice de rea Foliar
mximo, puede calcularse para das perfectamente claros y para das totalmente cerrados (bc y bo
respectivamente). La tasa de produccin mxima de biomasa bruta (bgm), es entonces calculada como
sigue:
bgm = F x bo + (1 - F) x bc
Los valores de bc y bo para diferentes latitudes, calculadas por De Witt, se reportan en el Cuadro 2.22,
con la aclaracin de que estos datos corresponden a una Fotosntesis Mxima (fm) de 20 Kg./ha/h. Las
modificaciones que se tienen que efectuar cuando se tiene un cultivo con una fm mayor o menor de 20
se explicaran en el tema sobre estimacin de rendimientos.
Tambin se han generado modelos para estimar los valores de Ac, bc y bo a partir de la R A, como se
reportan en el Cuadro 2.23.

36

Radiacin Solar

Cuadro 2.23. Modelos para estimar a Ac, Radiacin fotositticamente activa en das despejados, en
cal/cm2/da, bc y bo, Tasas de Fotosntesis bruta diaria en das despejados y das cerrados,
en Kg/ha/da, como funcin de la RA en cal/cm2/da.
MODELO

R2

Ac = 0.4566 (RA) 34.815

0.9921

Bc = 0.4051 (RA) + 80.334

0.9759

Bo = 0.2469 (RA) + 15.277

0.9896

Por la importancia que tiene la Radiacin Fotosintticamente Activa (RFA), a continuacin se

menciona una frmula para su clculo, propuesta por Driessen y Konijn (1992):
La Radiacin Fotosintticamente Activa (RFA) en el extremo superior de la atmsfera se calcula por:
RFA = 0.5 * [CS * (1 + 0.033 * cos(2 * PI * DAY / 365))] * RDN
con:
RDN = SSEN + 24 * CCOS * (1 - (SSCC)2 )0.5 / (DL * PI)
donde:
RFA = Radiacin Fotosintticamente Activa, fuera de la atmsfera, J/m2/s;
CS = Constante Solar (CS = 1353 J/m2/s),
RDN = Fraccin de la CS a la Latitud LAT y el da DAY.
Se recomienda al lector revisar dentro de este captulo las definiciones de DL, SSCC, SSEN, CCOS,
LAT y DAY.
Ejemplo: Supongamos que para Chapingo, Mxico, deseamos estimar la tasa de produccin mxima
de biomasa bruta para el cultivo de maz cuyo perodo vegetativo esta entre Mayo y Septiembre. Los
datos necesarios son la Radiacin global (Rg), que puede ser estimada a partir de los mtodos
sealados, particularmente en este caso, se empleo el mtodo de Allen; los otros datos, Radiacin
Fotosintticamente Activa (Ac), la Tasa de Fotosntesis en das muy claros (bc) y en das cerrados
(bo), se generan a travs de interpolaciones del Cuadro 2.22, considerando la latitud de Chapingo (19
30) y los meses de inters. Resultando la tabla, reportada en el Cuadro 2.24.
Clculo de F:
F = (Ac - 0.5 Rg)/ 0.8 Ac = (384 - 0.5 x 494)/ 0.8x384 = 137 / 307.2 = 0.45
Ntese que los clculos se realizan con los valores promedio.
Cuadro 2.24. Datos necesarios para calcular la produccin de biomasa bruta durante en ciclo del cultivo de
Maz en Chapingo, Mxico.

37

Radiacin Solar

Dato
Rg
Ac
bc
bo

Mayo
537
390
456
244

Junio
508
396
464
248

Julio
490
394
462
247

Agosto
484
386
450
242

Septiembre
450
354
423
224

Promedio
494
384
451
241

Clculo de bgm:
bgm = F x bo + (1 - F) x bc = 0.45 x 241 + (1-0.45) x 451
= 108.45 + 248.05 = 356.5 = 357 kg/ha/d
Es decir, se puede producir en promedio, entre Mayo y Septiembre, una biomasa bruta mxima de 357
kg/ha diariamente y como se explicar posteriormente, parte de esta biomasa se pierde por la
respiracin y una parte proporcional se convierte en el producto aprovechable por el hombre.

38