Erick R. Rivera
Current Organizations / Institutes
The organizations/institutes from the scholar's latest article
Previous Organizations / Institutes
The organizations / institutes from the scholar's other articles in the past ten years
Research Subjects
Erick R. Rivera has a total of 53 co-authors. Here are the top co-authors by number of shared publications:
- Jorge A. Amador
Universidad de Costa Rica
8 shared publications - Eric Jose Alfaro Martínez
Escuela de Física, Centro de Investigaciones Geofísicas y Centro de Investigación en Ciencias del Mar y Limnología, Universidad de Costa Rica, San José, Costa Rica
3 shared publications - Tito Maldonado
School of Physics, University of Costa Rica, San Jose 11501, Costa Rica
3 shared publications - Ana María Durán-Quesada
Centro de Investigación en Contaminación Ambiental (CICA), Universidad de Costa Rica, Escuela de Física, San José, Costa Rica
2 shared publications B. Calderón
Centro de Investigaciones Geofísicas, Universidad de Costa Rica, 11501-2060 San José, Costa Rica
2 shared publicationsChristopher L. Castro
Univ Arizona, Dept Hydrol & Atmospher Sci, Tucson, AZ USA
2 shared publications
Latest Publications
- Open Access
- 19 March 2026
- Published by MDPI AG in Meteorology
- Vol. 5 (1) , 6
On the Interaction of Tropical Easterly Waves and the Caribbean Low-Level Jet Using Observed, ERA5 and WWLLN Data over the Intra-Americas Seas During OTREC 2019
Propagating easterly waves (EW) are analyzed here, within the dynamical environment of the Caribbean Low-Level Jet (CLLJ) using radiosondes from the Organization of Tropical East Pacific Convection (OTREC) field campaign, ERA5 reanalysis, and lightning from the World Wide Lightning Location Network (WWLLN) over – N, – W during 21 August–30 September 2019. Radiosondes resolve the vertical structure of the waves at San Andrés (Colombia), Limón and Santa Cruz–Guanacaste (Costa Rica), while ERA5 provides spatial–temporal continuity and vertically integrated diagnostics—namely, the vertically integrated moisture flux divergence (VIMFD) and the vertically integrated geopotential flux divergence (VIGFD). Lightning from WWLLN and precipitation from ERA5 and the Integrated Multi-satellite Retrievals for the Global Precipitation Measurement mission (GPM IMERG) offer independent convective proxies to track disturbances. Mean profiles from radiosondes and ERA5 show strong agreement at Limón and Guanacaste and some differences at San Andrés, yet all datasets capture coherent, phase-locked anomalies in zonal wind, meridional wind, temperature, humidity, vertical velocity and vorticity used to diagnose EW–CLLJ interactions. VIMFD, VIGFD, lightning and precipitation exhibit westward-propagating cores that align with the above anomalies, indicating that organized convection is coupled to the disturbances, whereas the mean state preconditions the environment to enable wave-induced upward motion. A robust vertical adjustment of the CLLJ is documented: the core shifts from near 925 hPa over the Caribbean Sea to about 700 hPa over the Eastern Tropical Pacific ( hPa). This feature is reproduced by a 30-year ERA5 climatology, consistent with jet-exit forcing and enhanced boundary-layer coupling over land. Conditions favorable for barotropic instability using the Rayleigh–Kuo criterion, were present over most of the period. A qualitative barotropic conversion proxy, computed from the eddy momentum covariance , shows positive values in the lower troposphere at Guanacaste and in the layer 850–700 hPa at San Andrés, suggesting mean-to-eddy momentum transfer, whereas the signal at Limón is weaker. Together, these results provide a physically consistent view of EW–CLLJ interactions across the IAS; therefore, a schematic of those mechanisms is proposed here. The results highlight the need for high-resolution modeling and full energy-budget analyses.- 19 March 2026
- Published by MDPI AG in Meteorology
- Vol. 5 (1) , 6
Open AccessOn the Interaction of Tropical Easterly Waves and the Caribbean Low-Level Jet Using Observed, ERA5 and WWLLN Data over the Intra-Americas Seas During OTREC 2019
Propagating easterly waves (EW) are analyzed here, within the dynamical environment of the Caribbean Low-Level Jet (CLLJ) using radiosondes from the Organization of Tropical East Pacific Convection (OTREC) field campaign, ERA5 reanalysis, and lightning from the World Wide Lightning Location Network (WWLLN) over – N, – W during 21 August–30 September 2019. Radiosondes resolve the vertical structure of the waves at San Andrés (Colombia), Limón and Santa Cruz–Guanacaste (Costa Rica), while ERA5 provides spatial–temporal continuity and vertically integrated diagnostics—namely, the vertically integrated moisture flux divergence (VIMFD) and the vertically integrated geopotential flux divergence (VIGFD). Lightning from WWLLN and precipitation from ERA5 and the Integrated Multi-satellite Retrievals for the Global Precipitation Measurement mission (GPM IMERG) offer independent convective proxies to track disturbances. Mean profiles from radiosondes and ERA5 show strong agreement at Limón and Guanacaste and some differences at San Andrés, yet all datasets capture coherent, phase-locked anomalies in zonal wind, meridional wind, temperature, humidity, vertical velocity and vorticity used to diagnose EW–CLLJ interactions. VIMFD, VIGFD, lightning and precipitation exhibit westward-propagating cores that align with the above anomalies, indicating that organized convection is coupled to the disturbances, whereas the mean state preconditions the environment to enable wave-induced upward motion. A robust vertical adjustment of the CLLJ is documented: the core shifts from near 925 hPa over the Caribbean Sea to about 700 hPa over the Eastern Tropical Pacific ( hPa). This feature is reproduced by a 30-year ERA5 climatology, consistent with jet-exit forcing and enhanced boundary-layer coupling over land. Conditions favorable for barotropic instability using the Rayleigh–Kuo criterion, were present over most of the period. A qualitative barotropic conversion proxy, computed from the eddy momentum covariance , shows positive values in the lower troposphere at Guanacaste and in the layer 850–700 hPa at San Andrés, suggesting mean-to-eddy momentum transfer, whereas the signal at Limón is weaker. Together, these results provide a physically consistent view of EW–CLLJ interactions across the IAS; therefore, a schematic of those mechanisms is proposed here. The results highlight the need for high-resolution modeling and full energy-budget analyses. - Open Access
- 30 April 2025
- Published by Springer Nature in Climatic Change
- Vol. 178 (5) , 1-21
Detection and attribution of trends of meteorological extremes in Central America
- 30 April 2025
- Published by Springer Nature in Climatic Change
- Vol. 178 (5) , 1-21
Open AccessDetection and attribution of trends of meteorological extremes in Central America
- Open Access
- 22 December 2022
- Published by Inter-Research Science Center in Climate Research
- Vol. 89, 61-83
Sensitivity of precipitation and atmospheric low-level circulation patterns to domain size and choice of parameterization schemes in RegCM4.4 over Central America
The sensitivity of regional climate model simulations to domain size and position is becoming increasingly important for generating reliable climate scenarios. In this study, the Central America CORDEX domain (CCA) at 50 km horizontal resolution with a relaxation zone of 10° around the boundaries (CCA+) was taken as the basis to increase domain size in the RegCM4.4 model. The low-level circulation and precipitation patterns over Central America, the western Caribbean Sea, and the eastern tropical Pacific do not show strong sensitivity to domain size changes for an area increment of 18, 32, and 52% with respect to the size of CCA+. The physical configuration in RegCM4.4 has a greater impact on the representation of relevant climate characteristics and atmospheric processes. Simulated 925 hPa winds over the Caribbean low-level jet (CLLJ) region show unrealistic winter and summer intensities, especially in July. Similar issues were found in other studies for the CCA domain, using different models and physical configurations. A reduction of the critical Richardson number in the selected planetary boundary layer scheme resulted in little change in the strength of the summer component of the CLLJ. The model simulations do not completely capture key regional precipitation patterns such as the mid-summer drought, due to limitations in adequately representing low-level circulation. However, simulations using the Grell cumulus parameterization perform relatively better than those using a mixed scheme (Grell over land-Emanuel over ocean).
- 22 December 2022
- Published by Inter-Research Science Center in Climate Research
- Vol. 89, 61-83
Open AccessSensitivity of precipitation and atmospheric low-level circulation patterns to domain size and choice of parameterization schemes in RegCM4.4 over Central America
The sensitivity of regional climate model simulations to domain size and position is becoming increasingly important for generating reliable climate scenarios. In this study, the Central America CORDEX domain (CCA) at 50 km horizontal resolution with a relaxation zone of 10° around the boundaries (CCA+) was taken as the basis to increase domain size in the RegCM4.4 model. The low-level circulation and precipitation patterns over Central America, the western Caribbean Sea, and the eastern tropical Pacific do not show strong sensitivity to domain size changes for an area increment of 18, 32, and 52% with respect to the size of CCA+. The physical configuration in RegCM4.4 has a greater impact on the representation of relevant climate characteristics and atmospheric processes. Simulated 925 hPa winds over the Caribbean low-level jet (CLLJ) region show unrealistic winter and summer intensities, especially in July. Similar issues were found in other studies for the CCA domain, using different models and physical configurations. A reduction of the critical Richardson number in the selected planetary boundary layer scheme resulted in little change in the strength of the summer component of the CLLJ. The model simulations do not completely capture key regional precipitation patterns such as the mid-summer drought, due to limitations in adequately representing low-level circulation. However, simulations using the Grell cumulus parameterization perform relatively better than those using a mixed scheme (Grell over land-Emanuel over ocean).
- Open Access
- 1 July 2022
- Published by Universidad Nacional de Costa Rica in Universidad en Diálogo: Revista de Extensión
- Vol. 12 (2) , 31-45
Aspectos científicos y sociales del diseño, construcción, instalación y operación de una estación meteorológica para la comunidad cabécar de Pú Ksogë, Talamanca, Limón
A mediados del 2018 un grupo de investigadores e investigadoras del Centro de Investigaciones Geofísicas y de la Escuela de Física de la Universidad de Costa Rica iniciaron un proyecto para compartir conocimientos sobre las ciencias del tiempo y el clima con la comunidad cabécar de Pú Ksogë (9.591°N, 82.978°W, 114.4 m s. n. m.), Talamanca, Limón. La iniciativa fue apoyada por líderes y autoridades educativas de esa comunidad indígena a finales de ese año. Con el aporte laboral local y la autorización respectiva, se instaló y dejó en operación una estación meteorológica automática (Siwá Etka) en terrenos de la escuela comunal a principios de 2019. Se describen aquí detalles de este proyecto y se muestran los primeros datos atmosféricos observados en la historia de esa comunidad. La falta de un adecuado control comunal y el desinterés de las autoridades de la escuela en el cuido del instrumental provocaron problemas que no permitieron darle continuidad al proyecto y la estación tuvo que ser desinstalada en octubre de 2019. La experiencia tuvo resultados positivos, pero otros aspectos como la conectividad a celular e internet, en ocasiones de muy baja calidad o inexistente, también incidieron en la decisión de descontinuar la iniciativa. Los contactos realizados con una asociación de desarrollo integral de la región podrían ser una alternativa para darle continuidad al proyecto.
- 1 July 2022
- Published by Universidad Nacional de Costa Rica in Universidad en Diálogo: Revista de Extensión
- Vol. 12 (2) , 31-45
Open AccessAspectos científicos y sociales del diseño, construcción, instalación y operación de una estación meteorológica para la comunidad cabécar de Pú Ksogë, Talamanca, Limón
A mediados del 2018 un grupo de investigadores e investigadoras del Centro de Investigaciones Geofísicas y de la Escuela de Física de la Universidad de Costa Rica iniciaron un proyecto para compartir conocimientos sobre las ciencias del tiempo y el clima con la comunidad cabécar de Pú Ksogë (9.591°N, 82.978°W, 114.4 m s. n. m.), Talamanca, Limón. La iniciativa fue apoyada por líderes y autoridades educativas de esa comunidad indígena a finales de ese año. Con el aporte laboral local y la autorización respectiva, se instaló y dejó en operación una estación meteorológica automática (Siwá Etka) en terrenos de la escuela comunal a principios de 2019. Se describen aquí detalles de este proyecto y se muestran los primeros datos atmosféricos observados en la historia de esa comunidad. La falta de un adecuado control comunal y el desinterés de las autoridades de la escuela en el cuido del instrumental provocaron problemas que no permitieron darle continuidad al proyecto y la estación tuvo que ser desinstalada en octubre de 2019. La experiencia tuvo resultados positivos, pero otros aspectos como la conectividad a celular e internet, en ocasiones de muy baja calidad o inexistente, también incidieron en la decisión de descontinuar la iniciativa. Los contactos realizados con una asociación de desarrollo integral de la región podrían ser una alternativa para darle continuidad al proyecto.
- Open Access
- 9 December 2020
- Published by MDPI AG in Atmosphere
- Vol. 11 (12) , 1333
A Sea Breeze Study during Ticosonde-NAME 2004 in the Central Pacific of Costa Rica: Observations and Numerical Modeling
Surface and upper air observations and MM5v3 simulations examined the structure and inland penetration of sea breeze (SB) along the Grande de Tárcoles river basin (GTRB), central Pacific, Costa Rica, for two different intensity regimes of the Caribbean Low-Level Jet (CLLJ). Data comprise the period of 1 July to 16 September 2004 from Ticosonde-North American Monsoon Experiment, and a local University of Costa Rica-National Meteorological Institute field campaign. Maximum precipitation occurs between 14:00–17:00 LST, showing a time lag of 2 to 3 h after the temperature maximum, suggesting that local diurnal heating is key to convection. July–August precipitation exhibited a rainfall decrease along GTRB due to the SB dynamical processes interaction with a strong CLLJ. The SB maximum inland incursion was 24 km, with no evidence of its penetration into the Central Valley. The MM5v3 simulations for two convective and boundary layer (BL) schemes captured some SB structure features along the GTRB. Comparison of model results with observed data shows deficiencies in the model representation of the surface flow near coastal regions. Differences may be the result of time lag model’s poor responses to actual early morning BL sea–land temperature gradients. MM5v3 configurations used in this study resulted in biased wind speed simulations.
- 9 December 2020
- Published by MDPI AG in Atmosphere
- Vol. 11 (12) , 1333
Open AccessA Sea Breeze Study during Ticosonde-NAME 2004 in the Central Pacific of Costa Rica: Observations and Numerical Modeling
Surface and upper air observations and MM5v3 simulations examined the structure and inland penetration of sea breeze (SB) along the Grande de Tárcoles river basin (GTRB), central Pacific, Costa Rica, for two different intensity regimes of the Caribbean Low-Level Jet (CLLJ). Data comprise the period of 1 July to 16 September 2004 from Ticosonde-North American Monsoon Experiment, and a local University of Costa Rica-National Meteorological Institute field campaign. Maximum precipitation occurs between 14:00–17:00 LST, showing a time lag of 2 to 3 h after the temperature maximum, suggesting that local diurnal heating is key to convection. July–August precipitation exhibited a rainfall decrease along GTRB due to the SB dynamical processes interaction with a strong CLLJ. The SB maximum inland incursion was 24 km, with no evidence of its penetration into the Central Valley. The MM5v3 simulations for two convective and boundary layer (BL) schemes captured some SB structure features along the GTRB. Comparison of model results with observed data shows deficiencies in the model representation of the surface flow near coastal regions. Differences may be the result of time lag model’s poor responses to actual early morning BL sea–land temperature gradients. MM5v3 configurations used in this study resulted in biased wind speed simulations.
- Open Access
- 4 December 2020
- Published by MDPI AG in Atmosphere
- Vol. 11 (12) , 1317
Examination of WRF-ARW Experiments Using Different Planetary Boundary Layer Parameterizations to Study the Rapid Intensification and Trajectory of Hurricane Otto (2016)
Hurricane Otto (2016) was characterised by remarkable meteorological features of relevance for the scientific community and society. Scientifically, among the most important attributes of Otto is that it underwent a rapid intensification (RI) process. For society, this cyclone severely impacted Costa Rica and Nicaragua, leaving enormous economic losses and many fatalities. In this study, a set of three numerical simulations are performed to examine the skill of model estimations in reproducing RI and trajectory of Hurricane Otto by comparing the results of a global model to a regional model using three different planetary boundary layer parameterizations (PBL). The objective is to set the basis for future studies that analyse the physical reasons why a particular simulation (associated with a certain model setup) performs better than others in terms of reproducing RI and trajectory. We use the regional model Weather Research and Forecasting—Advanced Research WRF (WRF-ARW) with boundary and initial conditions provided by the Global Forecast System (GFS) analysis (horizontal resolution of 0.5 degrees). The PBL used are the Medium Range Forecast, the Mellor-Yamada-Janjic (MYJ), and the Yonsei University (YSU) parameterizations. The regional model is run in three static domains with horizontal grid spacing of 27, 9 and 3 km, the latter covering the spacial extent of Otto during the simulation period. WRF-ARW results improve the GFS forecast, in almost every aspect evaluated in this study, particularly, the simulated trajectories in WRF-ARW show a better representation of the cyclone path and movement compared to GFS. Even though the MYJ experiment was the only one that exhibited an abrupt 24-h change in the storm’s surface wind, close to the 25-knot threshold, the YSU scheme presented the fastest intensification, closest to reality.
- 4 December 2020
- Published by MDPI AG in Atmosphere
- Vol. 11 (12) , 1317
Open AccessExamination of WRF-ARW Experiments Using Different Planetary Boundary Layer Parameterizations to Study the Rapid Intensification and Trajectory of Hurricane Otto (2016)
Hurricane Otto (2016) was characterised by remarkable meteorological features of relevance for the scientific community and society. Scientifically, among the most important attributes of Otto is that it underwent a rapid intensification (RI) process. For society, this cyclone severely impacted Costa Rica and Nicaragua, leaving enormous economic losses and many fatalities. In this study, a set of three numerical simulations are performed to examine the skill of model estimations in reproducing RI and trajectory of Hurricane Otto by comparing the results of a global model to a regional model using three different planetary boundary layer parameterizations (PBL). The objective is to set the basis for future studies that analyse the physical reasons why a particular simulation (associated with a certain model setup) performs better than others in terms of reproducing RI and trajectory. We use the regional model Weather Research and Forecasting—Advanced Research WRF (WRF-ARW) with boundary and initial conditions provided by the Global Forecast System (GFS) analysis (horizontal resolution of 0.5 degrees). The PBL used are the Medium Range Forecast, the Mellor-Yamada-Janjic (MYJ), and the Yonsei University (YSU) parameterizations. The regional model is run in three static domains with horizontal grid spacing of 27, 9 and 3 km, the latter covering the spacial extent of Otto during the simulation period. WRF-ARW results improve the GFS forecast, in almost every aspect evaluated in this study, particularly, the simulated trajectories in WRF-ARW show a better representation of the cyclone path and movement compared to GFS. Even though the MYJ experiment was the only one that exhibited an abrupt 24-h change in the storm’s surface wind, close to the 25-knot threshold, the YSU scheme presented the fastest intensification, closest to reality.
- Open Access
- 1 April 2018
- Published by Instituto de Ciencias de la Atmosfera y Cambio Climatico in Atmósfera
- Vol. 31 (2) , 185-188
Putting into action the REGCM4.6 regional climate model for the study of climate change, variability and modeling over Central America and Mexico
Putting into action the REGCM4.6 regional climate model for the study of climate change, variability and modeling over Central America and Mexico
- 1 April 2018
- Published by Instituto de Ciencias de la Atmosfera y Cambio Climatico in Atmósfera
- Vol. 31 (2) , 185-188
Open AccessPutting into action the REGCM4.6 regional climate model for the study of climate change, variability and modeling over Central America and Mexico
Putting into action the REGCM4.6 regional climate model for the study of climate change, variability and modeling over Central America and Mexico
- 1 June 2017
- Published by Elsevier in Journal of Hydrology
- Vol. 549, 194-207
Framework for incorporating climate change on flood magnitude and frequency analysis in the upper Santa Cruz River
- 1 June 2017
- Published by Elsevier in Journal of Hydrology
- Vol. 549, 194-207
Framework for incorporating climate change on flood magnitude and frequency analysis in the upper Santa Cruz River
- Open Access
- 2 March 2016
- Published by Universidad de Costa Rica in Revista de Biologia Tropical
- Vol. 64 (1) , 23-57
The easternmost tropical Pacific. Part II: Seasonal and intraseasonal modes of atmospheric variability
The easternmost tropical Pacific. Part II: Seasonal and intraseasonal modes of atmospheric variability
- 2 March 2016
- Published by Universidad de Costa Rica in Revista de Biologia Tropical
- Vol. 64 (1) , 23-57
Open AccessThe easternmost tropical Pacific. Part II: Seasonal and intraseasonal modes of atmospheric variability
The easternmost tropical Pacific. Part II: Seasonal and intraseasonal modes of atmospheric variability
- Open Access
- 2 March 2016
- Published by Universidad de Costa Rica in Revista de Biologia Tropical
- Vol. 64 (1) , 1-22
The easternmost tropical Pacific. Part I: A climate review
The tropics are characterized by a variety of atmospheric and oceanic systems dominated by multi-scale interaction processes. This is Part I of a two-part review study on climate and climate variability of the Eastern Tropical Pacific (ETP). Here, the mean fields of incoming short wave radiation, surface energy fluxes, sea surface temperature, sea level pressure, salinity, wind patterns at 10m height, wind stress curl, precipitation, and evaporation are analyzed and synthetized using available data including that from the last decade. Given the strong interaction between the ETP and Caribbean Sea-Atlantic Ocean, mean field discussions are presented for two different domains, a relatively large domain from 24° S - 36° N and between 129° W - 17° W, and a smaller domain embedding only the Caribbean Sea and the easternmost part of the ETP. Most results on the climate of these two regions can be used as the base line for climate change studies. Interannual variability of tropical cyclones is also investigated over the domain of the latter two basins in the smaller domain. The study is complemented with a short review of some low frequency modes, such as, El Niño-Southern Oscillation, the Pacific Decadal Oscillation, and the Atlantic Multi-decadal Oscillation. Such low-frequency modes are known to modulate regional systems, including tropical cyclone frequency.
- 2 March 2016
- Published by Universidad de Costa Rica in Revista de Biologia Tropical
- Vol. 64 (1) , 1-22
Open AccessThe easternmost tropical Pacific. Part I: A climate review
The tropics are characterized by a variety of atmospheric and oceanic systems dominated by multi-scale interaction processes. This is Part I of a two-part review study on climate and climate variability of the Eastern Tropical Pacific (ETP). Here, the mean fields of incoming short wave radiation, surface energy fluxes, sea surface temperature, sea level pressure, salinity, wind patterns at 10m height, wind stress curl, precipitation, and evaporation are analyzed and synthetized using available data including that from the last decade. Given the strong interaction between the ETP and Caribbean Sea-Atlantic Ocean, mean field discussions are presented for two different domains, a relatively large domain from 24° S - 36° N and between 129° W - 17° W, and a smaller domain embedding only the Caribbean Sea and the easternmost part of the ETP. Most results on the climate of these two regions can be used as the base line for climate change studies. Interannual variability of tropical cyclones is also investigated over the domain of the latter two basins in the smaller domain. The study is complemented with a short review of some low frequency modes, such as, El Niño-Southern Oscillation, the Pacific Decadal Oscillation, and the Atlantic Multi-decadal Oscillation. Such low-frequency modes are known to modulate regional systems, including tropical cyclone frequency.