(pdf)
SciELO 
SciELO 
Permalink
INTRODUCCIÓN
El estudio de la energía solar ha ido evolucionando conforme los años. No obstante, hace poco más de 50 años comenzó el interés por la radiación solar en México (Hernández et al., 1991). El término radiación solar es comúnmente utilizado para referirse a la energía proveniente del sol en forma de radiación electromagnética, sin embargo, existen otras definiciones que componen la nomenclatura de la energía solar. En otras palabras, el concepto de radiación solar es empleado erróneamente para referirse a términos como irradiación e irradiancia. Por un lado, irradiación es la energía que, en forma de radiación, se integra durante un intervalo de tiempo en una superficie, generalmente sus unidades son
Las primeras aproximaciones superficiales desarrolladas para México fueron mapas de isohelias por interpolación a mano libre (Hernández et al., 1991), donde se ilustra la distribución espacial desigual de la irradiación, es decir, que existen zonas con mayor incidencia de energía solar que otras (Estrada-Cajigal, 2008). Estas variaciones son causadas por factores geográficos, como la orografía, latitud, altitud y paisaje, a su vez atmosféricos, como la presencia de nubes, aerosoles, entre otros (Hernández et al., 1991).
El factor de nubosidad se ha relacionado ampliamente con la dispersión e incidencia de la irradiación global (Calbó et al., 2016). Por lo que, a partir de las propiedades y comportamiento de las nubes presentes, surgió la idea de implementar satélites meteorológicos geoestacionarios para estimar los valores de la energía solar (Estrada-Cajigal, 2008). En México las observaciones superficiales de nubosidad se realizan mediante el conteo de días despejados, medio nublados o nublados cerrados en los observatorios atmosféricos (Tejeda-Martínez y Gómez-Azpeitia, 2015). Dichos datos de cobertura nubosa en el cielo son necesarios para perfeccionar los métodos de estimación, empero, con el tiempo se han dejado de registrar continuamente. Además, Díaz-Félix (2012) menciona que estas observaciones son factibles sólo si se toman máximo a 100 km de distancia de los puntos de medición de irradiancia.
Con la recopilación de ambos datos meteorológicos se elaboran numerosos productos (mapas de distribución, gráficos de energía disponible, entre otros) que se han integrado a investigaciones y recopilaciones como el Atlas Solar de la República Mexicana (1991) y el Prontuario Solar (2015). Donde más allá de temas de relevancia científica e histórica se describen softwares y bases de datos de irradiación solar, determinando que el comportamiento solar en superficie está constantemente ligado a medir sus componentes, ya que su incidencia no es la misma a lo largo del país (Tejeda-Martínez y Gómez-Azpeitia, 2015). Pese a esto, en el país resultan insuficientes los puntos de observación con utilidad en el desarrollo de proyectos científicos, energéticos, entre otros. La interpretación de los valores por satélites meteorológicos geoestacionarios continuó siendo una herramienta significativa para sustituir las mediciones in situ, ya que algunas organizaciones de monitoreo atmosférico se encargaron de crear plataformas web para visualizar y descargar los registros de cada satélite a nivel de la superficie. Sumado a ello, existen empresas y organismos que cuentan con información privada del recurso solar mexicano para fines comerciales que limitan ampliamente la utilización de datos de calidad.
Ante esta situación, actualmente las bases de datos por reanálisis son el método más empleado gracias a su accesibilidad y descarga gratuita. Los diferentes tipos de reanálisis integran principalmente observaciones de sitios de todo el mundo, en un conjunto de datos completo y consistente, mediante la asimilación simultáneamente restringida por leyes físicas a partir de modelos numéricos de predicción del tiempo (Habte et al., 2017). Entre los principales reanálisis se encuentran ERA5, The National Centers for Environmental Prediction con the National Center for Atmospheric Research (NCEP/NCAR), The Modern-Era Retrospective analysis for Research and Applications Version 2 (MERRA 2) y Prediction of Worldwide Energy Resources (POWER). En conclusión, existen tres diferentes métodos para la generación de datos solarimétricos en México. Estos son la disposición de una red de instrumentos en superficie, información de satélites geoestacionarios meteorológicos y productos de modelos de estimación por reanálisis. Se encuentran dos problemas primordiales al ejecutar los conjuntos generados por sensores remotos y reanálisis: en el primero se mantiene un margen de error y calidad desconocido, en el segundo destaca su falta de disponibilidad temporal-geográfica del país.
Considerando al sol como la mayor fuente inagotable de energía, la implementación de estos datos está lejos de acabar (Petrarca, 2021). Es necesario realizar un estudio para conocer las bases públicas más competentes de estimaciones solares y de nubosidad para el caso de México, puesto que, actualmente no existen trabajos de investigación que generen dicha información.
Según lo expuesto, el objetivo de este estudio es evaluar las bases de datos disponibles de irradiación e irradiancia. De manera que se determinará el error de estimación y adaptabilidad de cada una para seleccionar aquellas que alcancen los mejores resultados replicando el comportamiento de cada parámetro solar. Para el caso de la nubosidad, se complementará la información solarimétrica con las estimaciones satelitales del área total de nubes en los puntos de evaluación. Mostrando como el comportamiento espaciotemporal de la energía solar en superficie es influenciado por la presencia o ausencia de nubosidad.
El presente trabajo pretende apoyar con nuevas técnicas de adquisición de bases de datos solarimétricos, para su posterior aplicación climatológica, agrícola, energética, entre otras según el usuario en cuestión.
METODOLOGÍA
Recopilación de las bases de datos
Por razones de vigencia y continuidad en los conjuntos encontrados, se establecieron las siguientes cinco características específicas para verificar la confiabilidad de la fuente y la calidad de los datos salientes:
Manejo de datos solarimétricos y nubosidad
Las mediciones horarias de irradiancia
El valor estándar utilizado para determinar la cantidad anual de irradiación fue determinado con los promedios mensuales anteriores, se calculó la irradiación promedio anual mediante la ecuación 2:
Para el caso de la nubosidad, el parámetro ocupado fue fracción total de nubosidad, conocida con otros términos como “cantidad de nubes” o “cobertura nubosa”, ya que, de acuerdo con Enríquez-Alonso (2016) es la propiedad que describe esta condición en general. La fracción total de nubosidad se expresa en porcentaje (%) de 0 a 100, para este caso se integran los valores en un promedio mensual. En la Tabla I se estableció el porcentaje equivalente al número de octas para los términos “despejado”, “parcialmente nublado”, “nublado” y “completamente nublado” según la Organización Meteorológica Mundial (OMM).
Tabla I Clasificación del porcentaje de fracción de nubosidad.
| Porcentaje | Categoría |
| 0-30% | Despejado |
| 31%-60% | Parcialmente nublado |
| 61%-90% | Nublado |
| 91%-100% | Completamente nublado |
Fuente: elaboración propia (transformación de número de octas a porcentaje) a partir de datos de la Organización Meteorológica Mundial (OMM), http://worldweather.wmo.int/oktas.htm.
Fuentes de datos en superficie
Se efectuó una búsqueda exhaustiva en los sitios web oficiales de las instituciones nacionales que realizan mediciones horarias de irradiancia y que cumplen con el PROY-NMX-AA-166. Cabe recalcar, que mientras se realizó esta evaluación la página oficial del SMN sufrió una intromisión en su servidor, debido a esto, se inhabilitó la proporción de cualquier tipo de datos provenientes de la red de Estaciones Meteorológicas Automáticas (EMA).
Finalmente se trabajó con las bases de datos en superficie de las siguientes dependencias:
Métricos estadísticos de comparación
Con la finalidad de realizar la evaluación entre las series de datos y valores estimados, se implementaron los siguientes cuatro métodos estadísticos donde Mo representa el valor estimado, Me el valor real y n el número total de datos.
1) La raíz del error cuadrático medio (RMSD)
El RMSD es la raíz cuadrada de la diferencia de los valores estimados y los valores medidos in situ (ecuación 3). Especialmente, se utiliza para evaluar la precisión de modelos de estimación. De igual modo, se caracteriza por ser una medida sensible a cualquier valor extremo. Es decir, cada cambio instantáneo se muestra de manera explícita en las series de tiempo evaluadas (Porrini, 2017).
2) Error de sesgo medio (MBE)
El MBE captura las desviaciones promedio entre dos conjuntos de datos (ecuación 4). Los valores cercanos a 0 son los mejores, los valores negativos indican subestimación y los valores positivos indican sobreestimación (Turkish State Meteorological Service, s.f.).
3) Coeficiente de correlación de Pearson (r)
El Coeficiente de Correlación de Pearson (r) mide la covariación entre al menos dos variables (ecuación 5). Dicho de otro modo, si los valores de una característica cambian consistentemente conforme a los de la otra. De manera que, el coeficiente de correlación se empleó para cuantificar la covarianza entre los valores de las observaciones in situ y los valores estimados (Hernández-Lalinde et al., 2018).
4) Error porcentual absoluto medio (MAPE)
El MAPE es un indicador de desempeño de datos para establecer la precisión de valores estimados con respecto a valores observados en la superficie (ecuación 6). A menudo, se usa para evaluar la fiabilidad absoluta de los modelos de predicción, dando como resultado en qué porcentaje la información estimada es precisa (Crawford-Vidal et al., 2020).
RESULTADOS
Bases de datos
Se sometieron a revisión más de 25 ensambles de datos globales, pero solamente 15 de ellos cumplieron con los criterios anteriores. En la Tabla II se muestra el nombre completo de cada base de datos de sensores remotos o reanálisis y la abreviatura otorgada para fines prácticos.
Tabla II Bases de datos seleccionadas para la evaluación.
| Abreviatura | Nombre completo |
| POWER | Prediction of Worldwide Energy Resources |
| NSRDB | National Solar Radiation Database |
| CERES_EBAF | Clouds and the Earth’s Radiant Energy System_Energy Balanced and Filled |
| CERES_SYN | Clouds and the Earth’s Radiant Energy System_ SYNoptic TOA and surface fluxes and clouds |
| ERA 5 | European Centre for Medium-Range Weather Forecasts (ECMWF) Reanalysis Version 5 |
| NCEP (NCAR) | The National Centers for Environmental Prediction (NCEP) with the National Center for Atmospheric Research (NCAR) |
| NCEP (DOE) | The National Centers for Environmental Prediction (NCEP) with the National Center for Atmospheric Research (NCAR) Version 2 |
| MERRA 2 | The Modern-Era Retrospective analysis for Research and Applications Version 2 |
| CM SAF | The Satellite Application Facility on Climate Monitoring |
| FLDAS | Famine Early Warning Systems Network (FEWS NET) Land Data Assimilation System |
| CERES_CLD | Clouds and the Earth’s Radiant Energy System_ CldTypHist Algorithm |
| CLIMCAPS | Community Long-term Infrared Microwave Coupled Product System |
| MODIS | Moderate-Resolution Imaging Spectroradiometer |
| ISCCP | International Satellite Cloud Climatology Project |
| PATMOS-x | The Pathfinder Atmospheres Extended (PATMOS-x) Cloud Properties Climate Data Record (CDR) |
Fuente: elaboración propia a partir de la recopilación realizada en este estudio.
Cabe señalar que, con respecto al sitio Clouds and the Earth’s Radiant Energy System (CERES), se utilizaron tres subconjuntos como bases independientes, puesto que, la composición de cada algoritmo es diferente. De las 15 bases seleccionadas, 10 proporcionan irradiancia e irradiación promedio y 12 calculan la fracción total del área nubosa promedio diario y mensual (Tabla III).
Tabla III Clasificación del tipo de parámetro proporcionado por cada base de datos.
| Parámetros | Base de Datos | ||||||||||||||
| POWER | NSRDB | CERES_EBAF | CERES_SYN | ERA 5 | NCEP
(NCAR) |
NCEP
(DOE) |
MERRA 2 | CM SAF | FLDAS | CERES_CLD | CLIMCAPS | MODIS | ISCCP | PATMOS-x | |
| Irradiancia
e irradiación |
x | x | x | x | x | x | x | x | x | x | |||||
| Nubosidad | x | x | x | x | x | x | x | x | x | x | x | x | |||
Fuente: elaboración propia a partir de la información obtenida de cada base de datos de sensores remotos y reanálisis presentada en la Tabla II.
Puntos de referencia
La consulta y descarga de mediciones horarias de irradiancia se efectuó mediante los sitios oficiales de las dependencias universitarias y de investigación mencionadas. Para empezar, se llevó a cabo un análisis exploratorio donde se estableció que cada sitio de medición debía cumplir, mínimo, con el 75 % del total de datos hasta 2023. Por tal razón, fueron descartados 30 puntos de evaluación: tres observatorios atmosféricos, 12 estaciones mareográficas y 15 estaciones solarimétricas.
Las series temporales fueron analizadas para determinar el periodo de evaluación, donde la mayoría mostró variaciones iniciales debido al momento de la instalación del instrumental correspondiente. De manera que, la homogenización de las mediciones comenzó el 1 de junio de 2014, como resultado, el periodo determinado fue de junio 2014 a junio 2023. La evaluación se llevó a cabo a través de 20 puntos de referencia con observaciones de irradiancia y datos de irradiación a lo largo de todo México (Figura I). De los cuales 13 son observatorios atmosféricos pertenecientes a la RUOA, tres estaciones del SSM, tres estaciones de la Red mareográfica de CICESE y la estación del SAMMO (Tabla IV).
Fuente: elaboración propia en QGIS a partir de los puntos de referencia presentados en la Tabla IV.
Figura I Estaciones y observatorios utilizados en la evaluación.
Tabla IV Observatorios atmosféricos, estaciones mareográficas y solarimétricas.
| CLAVE | ESTACIÓN/OBSERVATORIO | LATITUD | LONGITUD | ALTITUD
(msnm) |
| AGSC | Observatorio Atmosférico Aguascalientes, Aguascalientes | 21.9157 | -102.319 | 1868 |
| ALTZ | Observatorio Atmosférico Altzomoni Amecameca de Juárez, Edo. de México | 19.1187 | -98.6552 | 3985 |
| BLAN | Estación Mareográfica Bahía de los Ángeles, B.C.N. | 28.57 | -113.33 | 0 |
| CDJ | Estación solarimétrica Ciudad Juárez, Chihuahua | 31.743 | -106.432 | 1127 |
| CHAM | Observatorio Atmosférico Chamela | 19.4983 | -105.0446 | 91 |
| COE | Estación solarimétrica Coneo, Michoacán | 19.813 | -101.695 | 1990 |
| ERNO | Observatorio Atmosférico Hermosillo, Sonora | 29.0814 | -110.9706 | 154 |
| JQRO | Observatorio Atmosférico Juriquilla, Querétaro | 20.703 | -100.4473 | 1945 |
| LPAZ | Estación Mareográfica La Paz, B.C.S | 24.1 | -110.21 | 0 |
| LTUX | Observatorio Atmosférico Las Tuxtlas, Veracruz | 18.5853 | -95.0752 | 180 |
| MAZT | Observatorio Atmosférico Mazatlán, Sinaloa | 23.1836 | -106.425 | 20 |
| MEDA | Observatorio Atmosférico Mérida, Yucatán | 20.9838 | -89.6452 | 25 |
| MORE | Observatorio Atmosférico Morelia, Michoacán | 19.6493 | -101.2221 | 1936 |
| SAMMO | Estación Mareográfica SAMMO | 20.8682 | -86.8683 | 5 |
| SISL | Observatorio Atmosférico Temixco, Morelos | 21.1645 | -90.0484 | 8 |
| SNFN | Estación Mareográfica San Felipe, B.C.N. | 30.59816 | -114.49325 | 0 |
| SSLO | Observatorio Atmosférico Saltillo, Coahuila | 25.3532 | -101.0332 | 1786 |
| TMIX | Observatorio Atmosférico Temixco, Morelos | 18.8405 | -99.2362 | 1253 |
| UNAM | Observatorio Atmosférico UNAM, Ciudad de México | 19.3262 | -99.1761 | 2280 |
| ZAC | Estación solarimétrica Zacatecas, Zacatecas | 22.772 | -102.644 | 2316 |
Fuente: elaboración propia con los puntos de referencia ilustrados en la Figura I.
Métricos estadísticos
De manera general se presentan los resultados en promedio del MBE y MAPE (Tabla V) en los 20 sitios de evaluación. Dando como primera aproximación que nueve de diez bases subestiman la irradiancia, CERES_EBAF y MERRA2 señalan el menor y mayor porcentaje de error respectivamente (MAPE). Asimismo, las diez subestiman la irradiación en menor magnitud a comparación del parámetro anterior y nuevamente CERES_EBAF con 5.8 % mantiene el porcentaje de error mínimo y MERRA2 con 11.1 % el máximo.
Tabla V Resultados promedio de MBE y MAPE de irradiancia (Mj/m2) e irradiación (kWatt-hora/m2 día).
| PARÁMETRO | IRRADIANCIA | IRRADIACIÓN | ||
| MÉTRICA | MBE | MAPE | MBE | MAPE |
| POWER | 0.2 | 7.7% | 0.12 | 6.9% |
| NSRDB | 0.8 | 8.8% | 0.26 | 7.5% |
| CERES_EBAF | 0.04 | 7.5% | 0.08 | 5.8% |
| CERES_SYN | 0.1 | 8.8% | 0.1 | 6.9% |
| ERA5 | 0.3 | 9.1% | 0.17 | 7.4% |
| NCEP (NCAR) | 1.5 | 12.5% | 0.53 | 10.8% |
| NCEP (DOE) | -0.4 | 10.9% | 0.04 | 8.7% |
| MERRA2 | 0.4 | 12.8% | 0.19 | 11.1% |
| CM SAF | 1 | 9.6% | 0.35 | 7.7% |
| FLDAS | 0.3 | 12.7% | 0.14 | 10.8% |
Fuente: elaboración propia con los promedios resultantes de MBE y MAPE en los puntos de referencia de la Tabla IV.
De manera general, se expone que las bases de datos evaluadas estiman con mayor precisión la irradiación que la irradiancia.
Evaluación por regiones
Con el fin de evaluar la variación de un parámetro meteorológico en términos de su distribución espaciotemporal, y de acuerdo con Diaz-Félix (2012), se sugiere dividir en subregiones o subperíodos de tiempo (meses o estaciones del año). Ya que, se presentan diferencias considerables entre los diversos tipos de clima y la orografía del país. Estos factores, como ya se mencionó, afectan la distribución de la irradiancia, a la par, se deben reflejar en las estimaciones suministradas por las bases de datos. A fin de generar un desarrollo más amplio, se seleccionaron siete de los 20 sitios sometidos a esta evaluación para cada zona principal del país. En la Figura II se ilustra la distribución geográfica de la estación mareográfica y los seis observatorios atmosféricos.
Fuente: elaboración en QGIS propia a partir de los puntos de referencia determinados.
Figura II Estación y observatorios utilizados para el análisis de regiones.
En el apartado de ANEXOS se presentan las tablas completas del RMSD, MBE, r y MAPE de la irradiancia mensual y la irradiación anual. Tras dicha evaluación estadística, se observó que los valores estimados con mejores resultados fueron en el observatorio de Hermosillo, Sonora (ERNO); en contrate los errores más grandes se obtuvieron en el observatorio de Las Tuxtlas, Veracruz (LTUX). Con base en el criterio anterior, se muestran los gráficos del comportamiento mensual de cada parámetro suministrados en ambos lugares.
En la Figura III, en ERNO se muestra que seis bases de datos replican con precisión los valores reales. A excepción de NCEP (NCAR) que sobrestima los datos durante la mayoría de los meses, MERRA2 y FLDAS que muestran un comportamiento desfasado.
Fuente: elaboración propia mediante los datos mensuales
Figura III Irradiancia mensual observada y estimada para a) ERNO y b) LTUX. promedio de la irradiancia a lo largo de un año.
En LTUX las diferencias entre los valores reales y estimados son ampliamente notorias, a excepción de los meses de otoño (septiembre, octubre y noviembre) donde las variaciones máximas fueron de
En la Figura IV se observan valores promedio anuales entre 5 a 6
Fuente: elaboración propia mediante la transformación de datos de irradiancia mensual a irradiación promedio anual en el periodo de evaluación.
Figura IV Irradiación promedio anual para a) ERNO y b) LTUX.
En este segundo apartado se ilustran los gráficos correspondientes al comportamiento de la fracción total de nubosidad promedio mensual, del periodo establecido, estimada por las doce bases datos.
En la Figura V la característica principal, en ambos casos, es que la mayoría de las estimaciones se comportan de manera similar. Aunque no se contó con datos de referencia en superficie, es notorio que los comportamientos son congruentes entre sí, lo que indica que la mayoría de las bases replican la fracción de nubosidad de manera eficiente. Los valores más altos se encuentran en julio, agosto y septiembre en ambos puntos. Lo anterior está ampliamente relacionado a la presencia de la dinámica atmosférica de cada región, según la época del año.
Fuente: elaboración propia donde mediante los porcentajes de fracción de nubosidad promedio a lo largo de un año.
Figura V Fracción de nubosidad promedio mensual para a) ERNO y b) LTUX.
La Figura VI muestra la cobertura nubosa presente con base en la variación de la irradiación en superficie. Es decir, la climatología conocida sobre el débil desarrollo de nubes en la zona noroeste y norte del país provoca un alto nivel entrante de energía solar. Ambos casos concuerdan con la zona del país en la que se ubican. Por un lado, ERNO, en el noroeste del territorio nacional, predominan meses despejados a parcialmente nublados y, LTUX, en el este, prevalecen los meses parcialmente nublados a nublados. En este mismo sentido, la distribución espacial de las nubes es difícil de estimar actualmente debido a los múltiples criterios que intervienen en ello. Para esta clasificación, los datos obtenidos por cada base de datos fueron relativamente similares, aunque se destaca un mayor número de coincidencias en ERNO, debido a su ubicación geográfica sobre los desiertos de México, donde mayormente existe una menor cantidad de nubosidad; contrariamente a LTUX, donde hay una mayor discrepancia por su cercanía al Golfo de México, y, por ende, una mayor dinámica para la formación de nubes.
Fuente: elaboración propia mediante la clasificación establecida de la OMM.
Figura VI Meses totalmente nublados, nublados, parcialmente nublados y despejados en promedio para a) ERNO y b) LTUX.
Adicionalmente, se plantea que el grado de nubosidad afecta la asimilación simultánea de los diferentes tipos de reanálisis y sensores remotos, reduciendo principalmente la eficiencia en que estos integran los datos observacionales, propiedades físicas y atmosféricas de sitios de todo el mundo. En el caso de LTUX, con mayor cobertura nubosa, se obtuvo una discrepancia mayormente alta de los valores solarimétricos estimados, en contraste, las métricas estadísticas en ERNO presentaron mayor exactitud y niveles inferiores al 60 % en fracción total de nubosidad.
Determinación de las bases de datos “ideales”
Después de analizar los resultados obtenidos por las bases de datos solarimétricos, se expone que POWER fue la mejor evaluada, debido a su menor margen de error en las estimaciones de la irradiancia, además sus valores de irradiación fueron muy similares al promedio de las mediciones. NSRDB obtuvo mejor desempeño en el cálculo de irradiación promedio anual, no obstante, sus actualizaciones permanecen disponibles hasta el año 2023. Esto limita de manera significativa la implementación de esta base de datos para futuros proyectos. CERES_EBAF fue el conjunto de irradiancia más preciso establecido por los métodos estadísticos y para el caso de la irradiación su desempeño fue muy competente. CERES_SYN obtuvo un alto rendimiento en las estimaciones de irradiación, llegando a destacar en las zonas centrales del país. ERA5 demostró un porcentaje de error pequeño, sin embargo, no superó a las otras bases de datos. NCEP (NCAR) subestimó la irradiación y fracción de nubosidad anual en gran magnitud, lo que ilustró una deficiencia para estimar las variables de interés. NCEP (DOE) siendo la versión mejorada del producto anterior, suministró un mayor grado de exactitud en los valores de irradiancia e irradiación, a su vez se encontraron similitudes en la estimación de la nubosidad. Asimismo, CM SAF demostró una implementación óptima para este caso, ya que, los valores salientes fueron muy parecidos a las mediciones de irradiancia en la mayoría de los puntos de evaluación. MERRA2 tuvo un desempeño deficiente, debido a que se adaptó mejor a ciertas áreas del país que a otras, por lo que, no es recomendado para su aplicación nacional. Por último, FLDAS fue la base con el peor rendimiento solarimétrico, careciendo de adaptabilidad en el territorio mexicano.
Tomando en cuenta lo anterior, se proponen las bases POWER y CERES_EBAF como las sobresalientes en la generación de datos solarimétricos.
Por un lado, POWER es un sitio web con una interfaz amigable que, a través de un mapa interactivo de respuesta rápida le permite al usuario descargar, en diferentes formatos, las estimaciones de irradiancia global, directa, difusa, entre otros. Se encuentra en constante mantenimiento al ser una base compuesta por varios modelos y datos de sistemas de percepción remota pertenecientes a la NASA. Cuenta con una resolución espacial de 0.5°× 0.5° (celda de 55.56 × 55.56 km). Aunque subestima ligeramente las estimaciones de irradiancia e irradiación, consiguió correlaciones fuertes (0.66 a 0.99) y generó del 1.7 % al 15.15 % de error absoluto para irradiancia y 0.8 % a 12.35 % para irradiación. Además, mostró mayor adaptabilidad espacial en el análisis por regiones.
La segunda opción CERES_EBAF es una página web con diseño burdo y simple, donde se debe crear una cuenta para tener acceso gratuito a los datos de flujos solares en superficie (irradiancia global, de onda corta y de onda larga en la atmósfera superior, entre otros). A diferencia de POWER, el formato de descarga disponible únicamente es netCDF3 y la resolución espacial es ligeramente inferior, con una cuadrícula regional de 1°× 1° (celda de 111.1 × 111.1 km). Los datos de irradiancia destacaron con diferencias máximas de 3.2
Para el caso del área total de nubosidad, MODIS fue de los resultados más similares para estimarla, en especial, en las zonas norte y noroeste. PATMOS-x implementó conjuntos definidos que incluyen parámetros más específicos de la cobertura nubosa del cielo, mostrando mayor definición temporal y espacial para cada punto. CERES_CLD es de los datos de nubosidad más precisos, al destacar que fueron similares a los productos de CERES_EBAF, CERES_SYN, ERA5, CM SAF y CLIMCAPS. Por último, ISCCP modeló el factor de nubosidad en promedio parcialmente aceptable, obteniendo un comportamiento subestimado a comparación de otras bases.
Las doce bases que estiman el porcentaje de fracción total de nubosidad se compararon entre sí, replicando el comportamiento atmosférico exitosamente. Cabe destacar que, a falta de observaciones es difícil determinar concretamente los porcentajes de nubosidad. A pesar de ello, se destacan aquellos productos que tenían un nivel superior de resolución espaciotemporal para generar porcentajes con mayor exactitud. Tales como CERES_CLD, MODIS y PATMOS- x.
DISCUSIÓN
La cantidad de información solarimétrica es un factor importante para realizar evaluaciones de datos modelados y estimados. Es bien conocido que, al tener un mayor número de datos se tendrán resultados más satisfactorios y aceptables estadísticamente. Asimismo, el periodo de tiempo no ha sido lo suficientemente continuo. En numerosos estudios se han ocupado rangos de años diferentes, como en el caso de dos evaluaciones realizadas al NSRDB. En dichos ejemplos, Sengupta et al. (2015) tomaron en cuenta siete años, mientras que Habte et al. (2017) ocuparon 17 años, por la antigüedad de los registros en superficie. Y más aún, otros autores han realizado productos tomando en cuenta menos datos, como en el caso de Tejeda-Martínez y Gómez-Azpeitia (2015), Díaz-Félix (2012) y Miranda-Jiménez (2016). En esta evaluación se consideró prudente ocupar los datos disponibles por minuto de los últimos nueve años, con el fin de obtener hallazgos significativos y precisos en la implementación de las bases de datos. Por otro lado, los 20 puntos de referencia ocupados son insuficientes para evaluar en todo el país, no obstante, como primer acercamiento este estudio plantea ser lo suficientemente aceptable para generar información utilizable. Las bases de irradiación e irradiancia estimada se han convertido en una importante herramienta, gracias a su aplicabilidad en los últimos años. Los estudios realizados por Tejeda-Martínez y Gómez-Azpeitia (2015), Díaz-Félix (2012) y Petrarca (2021) apoyan la implementación de la irradiación solar en México como fuente de energía renovable, debido a las características inagotables y constantes del sol. Del mismo modo, con dichos valores solarimétricos disponibles se pueden realizar proyectos científicos nacionales de gran importancia. Autores como Silvero et al. (2022) han caracterizado las tendencias a largo plazo de la irradiancia solar en superficie utilizando datos directamente de productos de reanálisis.
CONCLUSIONES
Se cumplió con el objetivo de evaluar la capacidad de cada base de datos de irradiación, irradiancia y nubosidad a lo largo del país. Se planteó que los registros de irradiancia global en superficie provienen de una cantidad limitada de sensores, provenientes de los observatorios/estaciones meteorológicas, solarimétricas y mareográficas, lo que afectó la aplicación de las métricas estadísticas. Actualmente, los productos siguen siendo limitados para el caso del país, a pesar de ello, las bases de datos evaluadas se adaptaron eficientemente a México como método de estimación en puntos donde no existen instrumentos para realizar mediciones.
Los valores de irradiancia e irradiación promedio, extraídos de las diez bases evaluadas, tuvieron una aproximación cercana al comportamiento de las mediciones en superficie. Destacando POWER y CERES_EBAF como las bases de datos mejor evaluadas para estimar los valores de los parámetros solares, principalmente por la resolución espacio-temporal y procedimiento de descarga accesible. Se estima que la base de datos POWER es la más recomendable para estimar los parámetros de irradiancia e irradiación para cualquier área de México.
Aunque este trabajo está más enfocado en los valores solarimétricos, será necesario en un futuro incluir con datos de calidad la fracción total de nubosidad para los puntos de evaluación, se debe generar y mejorar estos datos, estimaciones puntuales de nubosidad brindarán mayor conocimiento sobre el comportamiento espaciotemporal de la energía solar en las diferentes regiones del país.
De este modo, al comparar los mejores resultados estimados de irradiancia e irradiación con las estimaciones del parámetro nuboso se concluye que, la inexactitud de los valores estimados fue mayor en los puntos de medición donde se presenta mayor porcentaje de cobertura nubosa mensual. De manera que es importante tomar en cuenta este primer acercamiento para validar las afectaciones potenciales que tiene las bases de datos de sensores remotos y reanálisis bajo condiciones nubladas a completamente nubladas.
Las 15 fuentes presentan similitudes con potencial de ser utilizadas en México, sin embargo, se sugiere analizar una resolución temporal más específica, tal como la diaria u horaria. Asimismo, se recomienda utilizar parámetros de fracción total de nubosidad con la clasificación de nubes altas, medias y bajas debido a los tipos de nubes en capas individuales tienen efectos directos de la irradiancia entrante.
En resumen, el presente estudio temporal y espacial de la solarimetría en México demuestra un alto potencial de incidencia anual del orden de 4.5 a 6
