SciELO - Scientific Electronic Library Online

 
vol.29 special issueClimate estimates for Eastern Amazon with OLAM modelRadar-observed spatial and temporal rainfall variability near the Tapajós-Amazon confluence author indexsubject indexarticles search
Home Pagealphabetic serial listing  

Services on Demand

Journal

Article

Indicators

Related links

Share


Revista Brasileira de Meteorologia

Print version ISSN 0102-7786

Rev. bras. meteorol. vol.29 no.spe São Paulo Dec. 2014

http://dx.doi.org/10.1590/0102-778620130026 

ARTIGOS

 

Weather forecasting for Eastern Amazon with OLAM model

 

Previsão para o leste da Amazônia com o modelo OLAM

 

 

Renato Ramos da SilvaI; Adilson Wagner GandúII; Julia Clarinda CohenIII; Paulo KuhnIV; Maria Aurora MotaIII

IUniversidade Federal de Santa Catarina, Centro de Ciências Físicas e Matemáticas, Departamento de Física, Florianópolis, SC, Brasil.renato.ramos@ufsc.br
IIFundação Cearense de Meteorologia e Recursos Hídricos, Fortaleza, CE, Brasil.adwgandu@gmail.com
IIIUniversidade Federal do Pará, Departamento de Meteorologia, Belém, PA, Brasil. jcpcohen@gmail.com, aurora@ufpa.br
IVUniversidade Federal do Amazonas, Faculdade de Ciências Agrárias, Manaus, AM, Brasil.paulokuhn.1@gmail.com

 

 


ABSTRACT

The OLAM model has as its characteristics the advantage to represent simultaneously the global and regional meteorological phenomena using the application of a grid refinement scheme. During the REMAM project the model was applied for a few case studies to evaluate its performance on numerical weather prediction for the eastern Amazon region. Case studies were performed for the twelve months of the year of 2009. The model results for those numerical experiments were compared with the observed data for the region of study. Precipitation data analysis showed that OLAM is able to represent the average mean accumulated precipitation and the seasonal features of the events occurrence, but can't predict the local total amount of precipitation. However, individual evaluation for a few cases had shown that OLAM was able to represent the dynamics and forecast a few days in advance the development of coastal meteorological systems such as the squall lines that are one of the most important precipitating systems of the Amazon.

Keywords: Weather forecasting, OLAM, Amazon


RESUMO

O modelo OLAM tem como característica a vantagem de representar simultaneamente os fenômenos meteorológicos de escala global e regional através de um esquema de refinamento de grades. Durante o projeto REMAM, o modelo foi aplicado para alguns estudos de caso com objetivo de avaliar o desempenho do modelo na previsão numérica de tempo para a região leste da Amazônia. Estudos de caso foram feitos para os doze meses do ano de 2009. Os resultados do modelo para estes casos foram comparados com dados observados na região de estudo. A análise dos dados de precipitação mostrou que o modelo consegue representar a distribuição média da precipitação acumulada e os aspectos da sazonalidade da ocorrência dos eventos, mas não consegue prever individualmente a acumulação de precipitação local. No entanto, avaliação individual de alguns casos mostrou que o modelo OLAM conseguiu representar dinamicamente e prever, com alguns dias de antecedência, o desenvolvimento de fenômenos meteorológicos costeiros como as linhas de instabilidade, que são um dos mais importantes sistemas precipitantes da Amazônia.

Palavras-chave: Previsão numérica de tempo, OLAM, Amazônia


 

 

1. INTRODUCTION

Rainfall prediction at the tropical regions is very difficult because these regions are characterized by high downward solar radiation and atmospheric moisture producing strong turbulent fluxes and convection that are highly non-linear phenomena (Cotton et al, 2011).

Remote effects are constantly modifying the tropical weather. For instance, the seasonal migration of the Inter-tropical Convergence Zone (ITCZ) sets an important large scale variability that modifies the local meteorological conditions. In addition, the Madden-Julian Oscillation (MJO, Madden and Julian, 1972) can have important effects on the local precipitation (Souza and Ambrizi, 2006). Further climate variability may occur driven by the Pacific and Atlantic sea surface temperature conditions. For example, during the El Niño conditions, when the west Pacific is anomalously warm, precipitations usually are lower over the eastern Amazon (Grimm, 2003; Marengo, 2004). In addition, data analysis shows that a warm tropical Atlantic can be associated with weaker wind shear and less precipitation over the region (Nobre and Shukla, 1996).

Furthermore, the surface heterogeneity has an important effect on the local atmospheric circulations and rainfall distribution. The eastern Amazon is bounded by the Atlantic Ocean; the region is crossed by large rivers; it has a landscape covered by heterogeneous vegetation; and it is covered by topographic hills. These surface heterogeneities can produce horizontal temperature gradients that develop local atmospheric circulations affecting the convection and precipitation. Examples of these circulations are the river breeze (Silva Dias et al., 2004), the effects of topography on the convection development (Silva Dias et al., 2002), the cloud streets formation (Ramos da Silva et al., 2011); and the sea breeze and its influence on the squall lines formation (Cohen et al., 1995).

In general, the global models don't have enough spatial resolution to represent the mesoscale meteorological systems. Otherwise, the limited area models have better spatial resolution, but are slave of atmospheric boundary conditions from the global models. To fully represent the local mesoescale features and to account for the influence of remote and large scale systems simultaneously, a new modeling approach has to be applied.

The new Ocean-Land-Atmosphere Model (OLAM) has this capability. This model can represent the global atmospheric circulation using a hexagonal (or triangular) coarser icosahedral grid and simultaneously have an improved regional resolution through a grid refinement (Walko and Avissar, 2008a, b; Ramos da Silva et al., 2009). Recent studies with OLAM for the South America and the Amazon had shown good results on representing the regional climate variables (Medivgy et al., 2008, 2011).

During the project REMAM (Rede de Monitoramento e Pesquisa de Fenômenos Meteorológicos Extremos na Amazônia) a series of modeling experiments were performed with OLAM to evaluate its forecasting skill for the eastern Amazon.

The main goal of this study was to evaluate the OLAM model performance for several environmental conditions. First, the region of study, the model characteristics, and the numerical experiments are briefly described. Then, the results for all the case studies are evaluated trough a comparison with observed data. Finally, the model characteristics and performance are discussed.

 

2. METHODOLOGY

2.1 Region of study

The eastern Amazon region is located near the equator at the east of the South America (Figure 1). The Amazon River crosses the region and reaches the Atlantic Ocean at the east of the domain. The major land cover is rainforest, but urbanization and deforestation is modifying the landscape of the region. Topography of the region can reach about 700 meters over the north and south of the domain (Figure 1b). Dynamically, the atmospheric region is dominated by the trade winds that transport moist air from the Atlantic Ocean into the continental areas. The interface between land and the ocean sets local breeze circulations and convection producing squall lines that can propagate westward over the Amazon basin (Cohen et al., 1995). These squall lines are the major systems that produce the precipitation distribution over the Amazon basin (Cohen et al., 1995).

2.2 OLAM model

The OLAM model is an evolution of the mesoscale Regional Atmospheric Modeling System (RAMS, Cotton et al., 2003) but improves in several aspects such as with the possibility of a global grid that uses an icosahedral structure (Figure 1). Higher resolution can be achieved using a grid refinement for the region of interest. Therefore, the model is design with a unique grid and avoids problems that occur on the usual grid nesting schemes. The dynamical core for the grid structure can have a triangular or hexagonal grid mesh. The prognostic equations are solved using the finite volume technique and conserves momentum, mass and energy (Walko and Avissar, 2008a, b).

The major physical parametrizations adopted by OLAM are the same used by the RAMS model including the land-vegetation-atmosphere interaction (Walko et al., 2000a), the cloud microphysics scheme (Walko et al., 1995; Meyers et al., 1997; Walko et al., 2000b), the radiation transfer (Chen and Cotton, 1983; Harrington and Olson, 2001), and cumulus parametrization (Kain and Fritsch, 1990; Grell and Dévény, 2002). The RAMS mesoscale model was used in several applications for the South America studies such as on the deforestation impacts (Gandu et al., 2004; Ramos da Silva and Avissar, 2006; Ramos da Silva et al., 2008); lightning formation (Gonçalves et al., 2012), convection transport (Freitas et al., 2000), squall lines formation (Bender and Dias de Freitas, 2013), and weather forecasting (Freitas et al., 2005; Ramos da Silva et al., 2007; Longo et al., 2013).

2.3 Numerical design

In this study, the OLAM model was set with a global grid having a refined spatial resolution for the eastern Amazon (Figure 1). The global coarser grid was design having about 240 km of spacing between the hexagonal grid cells (Figure 1a). A five degree of refinement was set with the innermost grids cells having space on the order of 120 km, 60 km, 30 km and 15 km, respectively. The innermost high resolution domain was centered at 04 degree latitude South and 48 Longitude West, at the eastern Amazon (Figure 1b).

The atmospheric initial conditions for horizontal winds, temperature, geopotential height and moisture were obtained from the NCEP/NCAR (National Centers for Environment Prediction/National Center for Atmospheric Research) reanalysis project version 02 (Kanamitsu et al., 2002). Since OLAM is not a limited area model, nudging was not necessary. For OLAM, the coarser global grid produces the boundary conditions for the innermost model domain. Weekly sea surface temperature (SST) was obtained from NOAA (National Oceanic and Atmospheric Administration) as described by Reynolds et al. (2002). The SST fields were constantly updated during the atmospheric integration. Global data of topography, soil texture, and vegetation types, were used as a surface boundary condition.

To better evaluate the model, a series of case studies were performed for the year of 2009. The model was integrated for a total of seven days on twelve numerical experiments, one for each month for the year of 2009 (Table 1). These twelve cases were chosen to represent broad meteorological conditions including wet and dry periods of the year.

 

 

Precipitation data obtained from NCEP reanalysis and Tropical Rainfall Measuring Mission (TRMM) satellite estimates were used to evaluate the model results. The NCEP global precipitation data merges satellite and gauge measurement on a 2.5º × 2.5º latitude–longitude grid (Kanamitsu et al., 2002). TRMM precipitation data merges polar and geostationary satellites to provide a 0.25 x 0.25 latitude-longitude resolution dataset (Huffman et al., 2007). Meteorological images from the GOES satellite were used to identify important local precipitating systems. Also, precipitation data from the Instituto Nacional de Meteorologia (INMET) local weather stations were used to evaluate the local model performance.

 

3. RESULTS AND DISCUSSION

3.1 Precipitation forecast for eastern Amazon

The results of the 7-day accumulated precipitation simulations for each of the twelve case studies are presented at Figure 2. The model results and TRMM estimates show higher precipitation accumulation for the cases that occurred on the first half of the year (Figure 2 and Figure 3). The higher precipitation at this time of the year is caused by the presence of the ITCZ that is located over these low-latitude positions. Thus, winds and moisture converge over the region causing high frequency of convective cells. The precipitation distribution had influence not only from the presence of the ITCZ but also from local circulations caused by the local heterogeneity such as the sea-breeze and local topography. As compared with the TRMM precipitation estimates it is noted that the model was not able to fully simulate the high rainfall accumulation observed over the north of the domain including the Atlantic Ocean area. Statistical spatial correlation between the daily accumulated precipitation from the model and from the TRMM estimates shows low correlation between the simulation and the observation (Table 1). The best spatial correlations were obtained for the cases that occurred in March and May, which corresponds to a transition period when the ITCZ is moving north. However, for most of the cases, the model was able to represent the major precipitation accumulation distribution. Furthermore, the model was able to simulate reasonably well the domain averaged rainfall accumulation (Figure 3).

Figure 3 shows the forecasted mean precipitation accumulated over the domain. The results show that the model is able to estimate the 7-day accumulation of precipitation as compared with the observations from NCEP and TRMM. Between December and June there is more precipitation because the proximity of the ITCZ and the higher transport of moisture from the Atlantic Ocean. Between July and November, model and observations show very low precipitation as the ITCZ is located at far north latitudes.

To evaluate the overall domain daily rainfall rate a comparison between the model simulations and TRMM estimates are presented at Figure 4. The precipitation distribution from the model and TRMM show a similar number frequency of grid points and precipitation bins. Therefore, although the model can't predict the local total accumulation, it can estimate reasonable well the domain average and rainfall rate distribution.

To evaluate the local model performance on the precipitation forecast, the simulated results were compared with data from local standard weather stations. Figure 5 shows the locations of thirteen weather stations from the INMET network. Data comparison between the model precipitation accumulation for the seven days and the observations from the local stations are presented at Table 2. Except for a few cases, the results show that local simulated convection is not well correlated with the recorded data from the weather stations.

3.2 Individual case studies

To evaluate the model performance on representing the dynamics of the mesoscale local features, comparison between model results and satellite images were conducted for various selected cases. Satellite images analysis show that the interface between the continent and the ocean in combination with topography were responsible for several convective cells and squall lines formation.

Figure 6 represents the model forecast for the late afternoon of the 3rd day of forecast as compared with the GOES satellite image for 03 February 2009 20:00 UTC. A comparison between the model results and satellite image show that OLAM was able to predict with three days in advance the formation of the high convective clouds observed near the coast and at the south of the domain where there is a more steep topography (Figure 2). The total spatial mean precipitation forecasted for that day was 13.4 mm, while the TRMM estimates for the same area was 13.3 mm. These results show that OLAM had a very good predictability skill for that day.

Figure 7 shows the model forecast valid for the late afternoon of the 7th day of integration (i.e. 164 hours) as compared with the GOES satellite image for 08 March 2009 at 20:00 UTC. A comparison between the model results and the satellite image shows that OLAM was able to predict with seven days in advance the formation of the squall line at the coast and the convection over the south of the domain. The total spatial mean precipitation forecasted for that day was 8.7 mm, while the TRMM estimates for the same area was 4.5 mm. These results show that OLAM was able to predict with great time in advance the precipitation for that particular day.

Figure 8 shows the results from another case in which, the convective cells formed along with the coast and over the south of the domain as pictured by the satellite image from 07 April 2009 at 20:00 UTC. The model results show a 6-day prediction of these cells with great accordance with the satellite image. The mean accumulated precipitation for that day was 3.2 mm, while the TRMM data mean precipitation estimate for the same area was 6.2 mm. Again, these results show good model predictability skill.

Figure 9 shows the model forecast valid for the afternoon of 05 August 2009. During this period the ITCZ is far north and the satellite image showed that only a few convective cells formed over the continent. Furthermore, the winds are from southeast and are slightly modified by the sea breeze circulation. The forecasted mean accumulated precipitation for that day was 0.6 mm, while the TRMM mean estimates was 0.7 mm. Therefore, OLAM was able to forecast the observed system five days in advance with great performance.

 

4. CONCLUSIONS

A series of case studies were performed with the OLAM model to evaluate its forecast skill for the eastern Amazon region. The model had only an initial atmospheric condition and was integrated for a seven day run covering twelve cases, one for each month of the year 2009. The results showed that OLAM was not able to predict well the local precipitation accumulation, but was able to represent the dynamics of the major convective coastal systems such as the squall lines and the topographical induced convective cells. The results show that it is still a challenging problem to forecast well the convective cells that form in this region. Convection is itself a very nonlinear phenomenon. Indeed, a recent precipitation data measurement on several locations at the city of Belem had shown a very high variability on measurements obtained on stations very close to each other. The results suggest that new approaches such as ensemble techniques will have to be applied to improve the weather forecast for this region. In addition, better initial conditions are likely to help to improve de predictions.

In this study we used as initial atmospheric conditions the NCEP/NCAR reanalysis data that has low spatial resolution. A better initial condition based on data assimilation of meteorological observations from local stations, radar and satellite retrievals would likely improve the model predictability. This new approach is part of the next steps of this study. The Data Assimilation Research Testbed (DART, Anderson et al., 2009) is now being tested and implemented in conjunction with the OLAM and we expect to have and improved meteorological forecasting system for the region in the near future.

 

5. ACKNOWLEDGMENTS

The authors thanks to "Ministério de Ciência e Tecnologia" (MCT) and to "Financiadora de Estudos e Projetos" (FINEP), for the support to REMAM Project, and Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq) for further grants.

 

6. BIBLIOGRAPHIC REFERENCES

ANDERSON J. et al. The Data Assimilation Research Testbed: A Community Facility. Bulletin of the American Meteorological Society, v. 90, p. 1283–1296, 2009.         [ Links ]

BENDER, A.; DIAS DE FREITAS, E. Evaluation of BRAMS' Turbulence Schemes during a Squall Line Occurrence in São Paulo, Brazil. American Journal of Environmental Engineerin, v. 3, p. 1-7, 2013.         [ Links ]

CHEN, C.; COTTON, W. R., A One-dimensional Simulation of the Stratocumulus-capped Mixed Layer, Boundary-Layer Meteorology, v. 25, p. 289-321, 1983.         [ Links ]

COHEN, J. C. P.; SILVA DIAS. M. A. F; NOBRE, C. A.. Environmental conditions associated with Amazonian squall lines: a case study. Monthly Weather Review, v. 123, n. 11, p. 3163-3174, 1995.         [ Links ]

COTTON, W., BRYAN G., HEEVER S. Storm and cloud Dynamics, Academic Press, San Diego, 2011.         [ Links ]

COTTON W. R. et al. RAMS: Current Status and future directions. Meteorology and Atmospheric Physic, v. 82, p. 5–29, 2003.         [ Links ]

FREITAS, S. R. et al. A convective kinematic trajectory technique for low resolution atmospheric models. Journal of Geophysical Research. v. 105, n. D19, p. 24,375-24,386, 2000.         [ Links ]

FREITAS, E. D. et al.. A simple photochemical module implemented in RAMS for tropospheric ozone concentration forecast in the metropolitan area of Sao Paulo, Brazil: Coupling and validation. Atmospheric Environment, v. 39, p. 6352-6361, 2005.         [ Links ]

GANDU, A. W. et al.. .Simulations of deforestation in Eastern Amazon using a higher-resolution model. Theoretical and Applied Climatology, v. 78, p. 123-135, 2004.         [ Links ]

GRELL, G.A.; DÉVÉNYI, D.. A generalized approach to parameterizing convection combining ensemble and data assimilation techniques. Geophysical Research Letter, v. 29, n. 14, p. 1693, doi:10.1029/2002GL015311, 2002.         [ Links ]

GRIMM, A. M. The El Nino impact on the summer monsoon in Brazil: Regional processes versus remote influences. Journal of Climate, v. 16, p. 263–280, 2003.         [ Links ]

GONÇALVES, F. L. T. et al. Effect of bacterial ice nuclei on the frequency and intensity of lightning activity inferred by the BRAMS model. Atmospheric Chemistry and Physics, v. 12, p. 5677-5689, doi:10.5194/acp-12-5677-2012, 2012.         [ Links ]

HARRINGTON, J. Y.; OLSSON P. Q., A method for the Parameterization of Cloud Optical Properties in Bulk and Bin Microphysical Models: Implications for Arctic Cloudy Boundary Layers. Atmospheric Research., v. 57, p. 51-80, 2001.         [ Links ]

HUFFMAN, G. J. et al.. The TRMM Multi-satellite Precipitation Analysis: Quasi-global, multi-year, combined-sensor precipitation estimates at fine scale. Journal of Hydrometeorology, v. 8, n. 1, p. 38-55, 2007.         [ Links ]

KAIN, J. S.; FRITSCH, J. M. A One-Dimensional Entraining/Detraining Plume Model and Its Application in Convective Parameterization. Journal of the Atmospheric Sciences, v. 47, p. 2784–2802, 1990.         [ Links ]

KANAMITSU et al., NCEP-DEO AMIP-II Reanalysis (R-2). Bulletin of the American Meteorological Society, v. 83, p. 1631-1643, 2002.         [ Links ]

LONGO, K. M. et al. The chemistry CATT BRAMS model (CCATT BRAMS 4.5): a regional atmospheric model system for integrated air quality and weather forecasting and research. Geosci. Model Develop. Discussions, v. 6, p. 1173-1222, 2013.         [ Links ]

MADDEN, R. A.; JULIAN, P. R. Description of global-scale circulation cells in the tropics with a 40-50 days period. Journal of the Atmospheric Sciences, v. 29, p.1109-1123, 1972.         [ Links ]

MARENGO, J. A. Interdecadal variability and trends of rainfall across the Amazon basin. Theoretical and Applied Climatology, v. 78, p. 79-96, 2004.         [ Links ]

MEDVIGY D., WALKO R.; AVISSAR, R.. Effects of deforestation on spatiotemporal distributions of precipitation in South America, Journal of Climate, v. 24, p. 2147–2163, doi: 10.1175/2010JCLI3882.1, 2011.         [ Links ]

MEDVIGY D., WALKO R.; AVISSAR, R., Modeling interannual variability of the Amazon hydroclimate, Geophysical Research Letter, v. 35, p. L15817, doi:10.1029/2008GL034941, 2008.         [ Links ]

NOBRE, P.; SHUKLA, J. Variations of sea surface temperature, wind stress, and rainfall over the tropical Atlantic and South America. Journal of Climate, v. 9, p. 2464-2479, 1996.         [ Links ]

MEYERS, M. P. et al.. New RAMS cloud microphysics parameterization .2. The two-moment scheme. Atmospheric Research, v. 45, p. 3-39, 1997.         [ Links ]

RAMOS DA SILVA, R.; AVISSAR R. The hydrometeorology of a deforested region of the Amazon. Journal of Hydrometeorology, v. 7, n. 1028-1042, 2006.         [ Links ]

RAMOS DA SILVA, R. et al.. Cloud streets and land water interactions in the Amazon. Biogeochemistry (Dordrecht), p. 1-11, 2011.         [ Links ]

RAMOS DA SILVA, R.; SILVA DIAS, P. L.; MOREIRA, D. S.; SOUZA, E. B. Modelo OLAM – (Ocean Land Atmosphere Model) Descrição, Aplicações e Perspectivas, Revista Brasileira de Meteorologia, v. 24, n. 2, 144-157, 2009.         [ Links ]

RAMOS DA SILVA, R. et al. Progressos na Detecção e Previsão de Eventos Meteorológicos Extremos na Amazônia Oriental, Boletim da Sociedade Brasileira de Meteoroligia, SBMET, 2007.         [ Links ]

RAMOS DA SILVA, R., D. WERTH; AVISSAR R. Regional Impacts of Future Land-Cover Changes on the Amazon Basin Wet-Season Climate. Journal of Climate, v. 21, n. 6, 1153-1170, 2008.         [ Links ]

REYNOLDS, R. W. et al., An improved in situ and satellite SST analysis for climate. Journal of Climate, v. 15, 1609-1625, 2002.         [ Links ]

SILVA-DIAS, M. A. F. et al. Cloud and rain processes in a biosphere atmosphere interaction context in the Amazon Region. Journal of Geophysical Research, v. 107, n. D20, 8072, 2002.         [ Links ]

SILVA DIAS, M. A. F. et al. River breeze circulation in Eastern Amazon: observations and modeling results. Theoretical and Applied Climatology, v. 78, 111-121, 2004.         [ Links ]

SOUZA E.; AMBRIZZI T. Modulation of the intraseasonal rainfall over tropical Brazil by the Madden Julian Oscillation. International Journal of Climatology, v. 26, n. 13, p. 1759-1776, 2006.         [ Links ]

WALKO, R. L.; AVISSAR R. The Ocean–Land–Atmosphere Model (OLAM). Part I: Shallow-water testes. Monthly Weather Review, v. 136, p. 4033-4044, 2008a.         [ Links ]

WALKO, R.; AVISSAR R. The Ocean-Land-Atmosphere Model (OLAM). Part II: Formulation and Tests of the Non-Hydrostatic Dynamic Core, Monthly Weather Review, v. 136, p. 4045-4062, 2008b.         [ Links ]

WALKO, R. L. et al. Coupled atmosphere-biophysics-hydrology models for environmental modeling. Journal of Applied Meteorolog, v. 39, p. 931-944, 2000a.         [ Links ]

WALKO, R. L. et al. Efficient computation of vapor and heat diffusion between hydrometeors in a numerical model. Atmospheric Research, v. 53, p.171–183, 2000b.         [ Links ]

WALKO, R. L. et al. New RAMS cloud microphysics parameterization .1. The single-moment scheme. Atmospheric Research, v. 38, p. 29-62, 1995.         [ Links ]

 

 

Received August 2013
Accepted April 2014

Creative Commons License All the contents of this journal, except where otherwise noted, is licensed under a Creative Commons Attribution License