The role of convective parameterization in the simulation of a cyclone over the South Atlantic

Numerical simulations are carried out to examine the role of the Kuo and Kain-Fritsch (KF) cumulus parameterization schemes and dry dynamics on a cyclone development, in a weak baroclinic atmosphere, over subtropical South Atlantic Ocean. The initial phase of the cyclone development is investigated with a coarse horizontal mesh (75 km) and when the cyclone reaches the mature stage two different horizontal resolutions are used (75 and 25 km). The best performance simulation for the cyclone initial phase occurs when the Kuo convective scheme is applied, and this may be attributed to a greater diabatic warming in the troposphere. On the other hand, the dry simulation is not capable of simulating the correct location and intensity of the cyclone in its initial phase. During the mature phase, a cyclone over deepening occurs in the Kuo scheme experiment associated with larger latent heat release in a deep vertical column. The presence of downdraft currents in the KF scheme, which acts to cool and dry the lower levels, is essential to stabilize the atmosphere and to reproduce the nearest observation cyclone deepening rate. The largest cyclone deepening is found in the Kuo scheme high resolution experiment. This suggests that the KF convective scheme is less sensitive to the horizontal grid resolution. It was also revealed that the diabatic processes are crucial to simulate the observed features of this marine cyclone over subtropical region.

Precipitation diurnal cycle and summer climatology assessment over South America: An evaluation of Regional Climate Model version 3 simulations

Regional Climate Model version 3 (RegCM3) simulations of 17 summers (1988-2004) over part of South America south of 5 degrees S were evaluated to identify model systematic errors. Model results were compared to different rainfall data sets (Climate Research Unit (CRU), Climate Prediction Center (CPC), Global Precipitation Climatology Project (GPCP), and National Centers for Environmental Prediction (NCEP) reanalysis), including the five summers mean (1998-2002) precipitation diurnal cycle observed by the Tropical Rainfall Measuring Mission (TRMM)-Precipitation Radar (PR). In spite of regional differences, the RegCM3 simulates the main observed aspects of summer climatology associated with the precipitation (northwest-southeast band of South Atlantic Convergence Zone (SACZ)) and air temperature (warmer air in the central part of the continent and colder in eastern Brazil and the Andes Mountains). At a regional scale, the main RegCM3 failures are the underestimation of the precipitation in the northern branch of the SACZ and some unrealistic intense precipitation around the Andes Mountains. However, the RegCM3 seasonal precipitation is closer to the fine-scale analyses (CPC, CRU, and TRMM-PR) than is the NCEP reanalysis, which presents an incorrect north-south orientation of SACZ and an overestimation of its intensity. The precipitation diurnal cycle observed by TRMM-PR shows pronounced contrasts between Tropics and Extratropics and land and ocean, where most of these features are simulated by RegCM3. The major similarities between the simulation and observation, especially the diurnal cycle phase, are found over the continental tropical and subtropical SACZ regions, which present afternoon maximum (1500-1800 UTC) and morning minimum (0900-1200 UTC). More specifically, over the core of SACZ, the phase and amplitude of the simulated precipitation diurnal cycle are very close to the TRMM-PR observations. Although there are amplitude differences, the RegCM3 simulates the observed nighttime rainfall in the eastern Andes Mountains, over the Atlantic Ocean, and also over northern Argentina. The main simulation deficiencies are found in the Atlantic Ocean and near the Andes Mountains. Over the Atlantic Ocean the convective scheme is not triggered; thus the rainfall arises from the grid-scale scheme and therefore differs from the TRMM-PR. Near the Andes, intense (nighttime and daytime) simulated precipitation could be a response of an incorrect circulation and topographic uplift. Finally, it is important to note that unlike most reported bias of global models, RegCM3 does not trigger the moist convection just after sunrise over the southern part of the Amazon.

The Coupled Aerosol and Tracer Transport model to the Brazilian developments on the Regional Atmospheric Modeling System (CATT-BRAMS) - Part 1: Model description and evaluation

We introduce the Coupled Aerosol and Tracer Transport model to the Brazilian developments on the Regional Atmospheric Modeling System (CATT-BRAMS). CATT-BRAMS is an on-line transport model fully consistent with the simulated atmospheric dynamics. Emission sources from biomass burning and urban-industrial-vehicular activities for trace gases and from biomass burning aerosol particles are obtained from several published datasets and remote sensing information. The tracer and aerosol mass concentration prognostics include the effects of sub-grid scale turbulence in the planetary boundary layer, convective transport by shallow and deep moist convection, wet and dry deposition, and plume rise associated with vegetation fires in addition to the grid scale transport. The radiation parameterization takes into account the interaction between the simulated biomass burning aerosol particles and short and long wave radiation. The atmospheric model BRAMS is based on the Regional Atmospheric Modeling System (RAMS), with several improvements associated with cumulus convection representation, soil moisture initialization and surface scheme tuned for the tropics, among others. In this paper the CATT-BRAMS model is used to simulate carbon monoxide and particulate material (PM(2.5)) surface fluxes and atmospheric transport during the 2002 LBA field campaigns, conducted during the transition from the dry to wet season in the southwest Amazon Basin. Model evaluation is addressed with comparisons between model results and near surface, radiosondes and airborne measurements performed during the field campaign, as well as remote sensing derived products. We show the matching of emissions strengths to observed carbon monoxide in the LBA campaign. A relatively good comparison to the MOPITT data, in spite of the fact that MOPITT a priori assumptions imply several difficulties, is also obtained.

Robust relations between CCN and the vertical evolution of cloud drop size distribution in deep convective clouds

In-situ measurements in convective clouds (up to the freezing level) over the Amazon basin show that smoke from deforestation fires prevents clouds from precipitating until they acquire a vertical development of at least 4 km, compared to only 1-2 km in clean clouds. The average cloud depth required for the onset of warm rain increased by similar to 350 m for each additional 100 cloud condensation nuclei per cm(3) at a super-saturation of 0.5% (CCN0.5%). In polluted clouds, the diameter of modal liquid water content grows much slower with cloud depth (at least by a factor of similar to 2), due to the large number of droplets that compete for available water and to the suppressed coalescence processes. Contrary to what other studies have suggested, we did not observe this effect to reach saturation at 3000 or more accumulation mode particles per cm(3). The CCN0.5% concentration was found to be a very good predictor for the cloud depth required for the onset of warm precipitation and other microphysical factors, leaving only a secondary role for the updraft velocities in determining the cloud drop size distributions. The effective radius of the cloud droplets (r(e)) was found to be a quite robust parameter for a given environment and cloud depth, showing only a small effect of partial droplet evaporation from the cloud's mixing with its drier environment. This supports one of the basic assumptions of satellite analysis of cloud microphysical processes: the ability to look at different cloud top heights in the same region and regard their r(e) as if they had been measured inside one well developed cloud. The dependence of r(e) on the adiabatic fraction decreased higher in the clouds, especially for cleaner conditions, and disappeared at r(e)>=similar to 10 mu m. We propose that droplet coalescence, which is at its peak when warm rain is formed in the cloud at r(e)=similar to 10 mu m, continues to be significant during the cloud's mixing with the entrained air, cancelling out the decrease in r(e) due to evaporation.

Hybrid modeling of convective laminar flow in a permeable tube associated with the cross-flow process

The confined flows in tubes with permeable surfaces arc associated to tangential filtration processes (microfiltration or ultrafiltration). The complexity of the phenomena do not allow for the development of exact analytical solutions, however, approximate solutions are of great interest for the calculation of the transmembrane outflow and estimate of the concentration, polarization phenomenon. In the present work, the generalized integral transform technique (GITT) was employed in solving the laminar and permanent flow in permeable tubes of Newtonian and incompressible fluid. The mathematical formulation employed the parabolic differential equation of chemical species conservation (convective-diffusive equation). The velocity profiles for the entrance region flow, which are found in the connective terms of the equation, were assessed by solutions obtained from literature. The velocity at the permeable wall was considered uniform, with the concentration at the tube wall regarded as variable with an axial position. A computational methodology using global error control was applied to determine the concentration in the wall and concentration boundary layer thickness. The results obtained for the local transmembrane flux and the concentration boundary layer thickness were compared against others in literature. (C) 2007 Elsevier B.V. All rights reserved.

Parameterization and Evaluation of Predictions of DSSAT/CANEGRO for Brazilian Sugarcane

The DSSAT/CANEGRO model was parameterized and its predictions evaluated using data from five sugarcane (Sacchetrum spp.) experiments conducted in southern Brazil. The data used are from two of the most important Brazilian cultivars. Some parameters whose values were either directly measured or considered to be well known were not adjusted. Ten of the 20 parameters were optimized using a Generalized Likelihood Uncertainty Estimation (GLUE) algorithm using the leave-one-out cross-validation technique. Model predictions were evaluated using measured data of leaf area index (LA!), stalk and aerial dry mass, sucrose content, and soil water content, using bias, root mean squared error (RMSE), modeling efficiency (Eff), correlation coefficient, and agreement index. The Decision Support System for Agrotechnology Transfer (DSSAT)/CANEGRO model simulated the sugarcane crop in southern Brazil well, using the parameterization reported here. The soil water content predictions were better for rainfed (mean RMSE = 0.122mm) than for irrigated treatment (mean RMSE = 0.214mm). Predictions were best for aerial dry mass (Eff = 0.850), followed by stalk dry mass (Eff = 0.765) and then sucrose mass (Eff = 0.170). Number of green leaves showed the worst fit (Eff = -2.300). The cross-validation technique permits using multiple datasets that would have limited use if used independently because of the heterogeneity of measures and measurement strategies.

Climatic simulations of the eastern Andes low-level jet and its dependency on convective parameterizations

The South American low level jet (SALLJ) of the Eastern Andes is investigated with Regional Climate Model version 3 (RegCM3) simulations during the 2002-2003 austral summer using two convective parameterizations (Grell and Emanuel). The simulated SALLJ is compared with the special observations of SALLJEX (SALLJ Experiment). Both the Grell and Emanuel schemes adequately simulate the low level flow over South America. However, there are some intensity differences. Due to the larger (smaller) convective activity, the Emanuel (Grell) scheme simulates more intense (weaker) low level wind than analysis in the tropics and subtropics. The objectives criteria of Sugahara (SJ) and Bonner (BJ) were used for LLJ identification. When applied to the observations, both criteria suggest a larger frequency of the SALLJ in Santa Cruz, followed by Mariscal, Trinidad and Asuncin. In Mariscal and Asuncin, the diurnal cycle indicates that SJ occurs mainly at 12 UTCs (morning), while the BJ criterion presents the SALLJ as more homogenously distributed. The concentration into two of the four-times-a-day observations does not allow conclusions about the diurnal cycle in Santa Cruz and Trinidad. The simulated wind profiles result in a lower than observed frequency of SALLJ using both the SJ and BJ criteria, with fewer events obtained with the BJ. Due to the stronger simulated winds, the Emanuel scheme produces an equal or greater relative frequency of SALLJ than the Grell scheme. However, the Grell scheme using the SJ criterion simulates the SALLJ diurnal cycle closer to the observed one. Although some discrepancies between observed and simulated mean vertical profiles of the horizontal wind are noted, there is large agreement between the composites of the vertical structure of the SALLJ, especially when the SJ criterion is used with the Grell scheme. On an intraseasonal scale, a larger southward displacement of SALLJ in February and December when compared with January has been noted. The Grell and Emanuel schemes simulated this observed oscillation in the low-level flow. However, the spatial pattern and intensity of rainfall and circulation anomalies simulated by the Grell scheme are closer to the analyses than those obtained with the Emanuel scheme.

Resonant Wave Interactions in the Presence of a Diurnally Varying Heat Source

Resonant interactions among equatorial waves in the presence of a diurnally varying heat source are studied in the context of the diabatic version of the equatorial beta-plane primitive equations for a motionless, hydrostatic, horizontally homogeneous and stably stratified background atmosphere. The heat source is assumed to be periodic in time and of small amplitude [i.e., O(epsilon)] and is prescribed to roughly represent the typical heating associated with deep convection in the tropical atmosphere. In this context, using the asymptotic method of multiple time scales, the free linear Rossby, Kelvin, mixed Rossby-gravity, and inertio-gravity waves, as well as their vertical structures, are obtained as leading-order solutions. These waves are shown to interact resonantly in a triad configuration at the O(e) approximation, and the dynamics of these interactions have been studied in the presence of the forcing. It is shown that for the planetary-scale wave resonant triads composed of two first baroclinic equatorially trapped waves and one barotropic Rossby mode, the spectrum of the thermal forcing is such that only one of the triad components is resonant with the heat source. As a result, to illustrate the role of the diurnal forcing in these interactions in a simplified fashion, two kinds of triads have been analyzed. The first one refers to triads composed of a k = 0 first baroclinic geostrophic mode, which is resonant with the stationary component of the diurnal heat source, and two dispersive modes, namely, a mixed Rossby-gravity wave and a barotropic Rossby mode. The other class corresponds to triads composed of two first baroclinic inertio-gravity waves in which the highest-frequency wave resonates with a transient harmonic of the forcing. The integration of the asymptotic reduced equations for these selected resonant triads shows that the stationary component of the diurnal heat source acts as an ""accelerator"" for the energy exchanges between the two dispersive waves through the excitation of the catalyst geostrophic mode. On the other hand, since in the second class of triads the mode that resonates with the forcing is the most energetically active member because of the energy constraints imposed by the triad dynamics, the results show that the convective forcing in this case is responsible for a longer time scale modulation in the resonant interactions, generating a period doubling in the energy exchanges. The results suggest that the diurnal variation of tropical convection might play an important role in generating low-frequency fluctuations in the atmospheric circulation through resonant nonlinear interactions.

Impact of biomass burning aerosols on precipitation in the Amazon: A modeling case study

A study of the potential role of aerosols in modifying clouds and precipitation is presented using a numerical atmospheric model. Measurements of cloud condensation nuclei (CCN) and cloud size distribution properties taken in the southwestern Amazon region during the transition from dry to wet seasons were used as guidelines to define the microphysical parameters for the simulations. Numerical simulations were carried out using the Brazilian Development on Regional Atmospheric Modeling System, and the results presented considerable sensitivity to changes in these parameters. High CCN concentrations, typical of polluted days, were found to result in increases or decreases in total precipitation, depending on the level of pollution used as a reference, showing a complexity that parallels the aerosol-precipitation interaction. Our results show that on the grids evaluated, higher CCN concentrations reduced low-to-moderate rainfall rates and increased high rainfall rates. The principal consequence of the increased pollution was a change from a warm to a cold rain process, which affected the maximum and overall mean accumulated precipitation. Under polluted conditions, cloud cover diminished, allowing greater amounts of solar radiation to reach the surface. Aerosol absorption of radiation in the lower layers of the atmosphere delayed convective evolution but produced higher maximum rainfall rates due to increased instability. In addition, the intensity of the surface sensible heat flux, as well as that of the latent heat flux, was reduced by the lower temperature difference between surface and air, producing greater energy stores at the surface.

Bayesian inference for a skew-normal IRT model under the centred parameterization

Item response theory (IRT) comprises a set of statistical models which are useful in many fields, especially when there is interest in studying latent variables. These latent variables are directly considered in the Item Response Models (IRM) and they are usually called latent traits. A usual assumption for parameter estimation of the IRM, considering one group of examinees, is to assume that the latent traits are random variables which follow a standard normal distribution. However, many works suggest that this assumption does not apply in many cases. Furthermore, when this assumption does not hold, the parameter estimates tend to be biased and misleading inference can be obtained. Therefore, it is important to model the distribution of the latent traits properly. In this paper we present an alternative latent traits modeling based on the so-called skew-normal distribution; see Genton (2004). We used the centred parameterization, which was proposed by Azzalini (1985). This approach ensures the model identifiability as pointed out by Azevedo et al. (2009b). Also, a Metropolis Hastings within Gibbs sampling (MHWGS) algorithm was built for parameter estimation by using an augmented data approach. A simulation study was performed in order to assess the parameter recovery in the proposed model and the estimation method, and the effect of the asymmetry level of the latent traits distribution on the parameter estimation. Also, a comparison of our approach with other estimation methods (which consider the assumption of symmetric normality for the latent traits distribution) was considered. The results indicated that our proposed algorithm recovers properly all parameters. Specifically, the greater the asymmetry level, the better the performance of our approach compared with other approaches, mainly in the presence of small sample sizes (number of examinees). Furthermore, we analyzed a real data set which presents indication of asymmetry concerning the latent traits distribution. The results obtained by using our approach confirmed the presence of strong negative asymmetry of the latent traits distribution. (C) 2010 Elsevier B.V. All rights reserved.

Bias correction in a multivariate normal regression model with general parameterization

This paper derives the second-order biases Of maximum likelihood estimates from a multivariate normal model where the mean vector and the covariance matrix have parameters in common. We show that the second order bias can always be obtained by means of ordinary weighted least-squares regressions. We conduct simulation studies which indicate that the bias correction scheme yields nearly unbiased estimators. (C) 2009 Elsevier B.V. All rights reserved.

Controle de temperatura em intervenção cirúrgica abdominal convencional: comparação entre os métodos de aquecimento por condução e condução associada à convecção

JUSTIFICATIVA E OBJETIVOS: Hipotermia intra-operatória é complicação frequente, favorecida por operação abdominal. A eficácia da associação dos métodos de aquecimento por condução e convecção na prevenção de hipotermia e seus efeitos no período de recuperação pós-operatória foram os objetivos deste estudo. MÉTODO: Quarenta e três pacientes de ambos os sexos de 18 a 88 anos de idade, submetidos à laparotomia xifopúbica sob anestesia geral e monitorização da temperatura esofágica, foram distribuídos de modo aleatório em dois grupos de aquecimento: COND (n = 24), com colchão de circulação de água a 37°C no dorso e COND + CONV (n = 19), com a mesma condição associada à manta de ar aquecido a 42°C sobre o tórax e membros superiores. Analisados peso, sexo, idade, duração da operação e anestesia, temperaturas na indução anestésica (Mi), horas consecutiva (M1, M2), final da operação (Mfo) e anestesia (Mfa), entrada (Me-REC) e saída (Ms-REC) da recuperação pós-anestésica (SRPA), além das incidências de tremores e queixas de frio no pós-operatório. RESULTADOS: Os grupos foram semelhantes em todas as variáveis analisadas, exceto nas temperaturas em M2, M3, M4, Mfo e Mfa. O grupo COND reduziu a temperatura a partir da segunda hora da indução anestésica, mas o grupo COND + CONV só na quarta hora. Em COND, observou-se hipotermia na entrada e saída da SRPA. CONCLUSÕES: Associar métodos de aquecimento retardou a instalação e diminui a intensidade da hipotermia intra-operatória, mas não reduziu a incidência das queixas de frio e tremores.

Convecção linearmente organizada na área de Petrolina, semi-árido do Nordeste do Brasil: aspectos em meso e grande escala

É apresentado um estudo sobre sistemas convectivos linearmente organizados e observados por um radar meteorológico banda-C na região semi-árida do Nordeste do Brasil. São analisados três dias (27 a 29) de março de 1985, com ênfase na investigação do papel desempenhado por fatores locais e de grande escala no desenvolvimento dos sistemas. No cenário de grande escala, a área de cobertura do radar foi influenciada por um cavado de ar superior austral no dia 27 e por um vórtice ciclônico de altos níveis no dia 29. A convergência de umidade próxima à superfície favoreceu a atividade convectiva nos dias 27 e 29, enquanto que divergência de umidade próxima à superfície inibiu a atividade convectiva no dia 28. No cenário de mesoescala, foi observado que o aquecimento diurno é um fator importante para a formação de células convectivas, somando-se a ele o papel determinante da orografia na localização dos ecos. De maneira geral, as imagens de radar mostram os sistemas convectivos linearmente organizados em áreas elevadas e núcleos convectivos intensos envolvidos por uma área de precipitação estratiforme. Os resultados indicam que convergência do fluxo de umidade em grande escala e aquecimento radiativo, são fatores determinantes na evolução e desenvolvimento dos ecos na área de estudo.

Modelagem dos processos de remoção sulfato e dióxido de enxofre presente no particulado em diferentes localidades da região metropolitana de São Paulo

This article deals with the scavenging processes modeling of the particulate sulfate and the gas sulfur dioxide, emphasizing the synoptic conditions at different sampling sites in order to verify the domination of the in-cloud or below-cloud scavenging processes in the Metropolitan Area of São Paulo (RMSP). Three sampling sites were chosen: GV (Granja Viana) at RMSP surroundings, IAG-USP and Mackenzie (RMSP center). Basing on synoptic conditions, it was chosen a group of events where the numerical modeling, a simple scavenging model, was used. These synoptic conditions were usually convective cloud storms, which are usual at RMSP. The results show that the in-cloud processes were dominant (80%) for sulfate/sulfur dioxide scavenging processes, with below-cloud process indicating around 20% of the total. Clearly convective events, with total rainfall higher than 20 mm, are better modeled than the stratiform events, with correlation coefficient of 0.92. There is also a clear association with events presenting higher rainfall amount and the ratio between modeled and observed data set with correlation coefficient of 0.63. Additionally, the suburb sampling site, GV, as expected due to the pollution source distance, presents in general smaller amount of rainwater sulfate (modeled and observed) than the center sampling site, Mackenzie, where the characterization event explains partially the rainfall concentration differences.

Análise de jatos de baixos níveis associados aum sistema convectivo de mesoescala na américa do sul: um estudo de caso

The Levei Low Jet (LLJ) observed in the Porto Alegre metropolitan region, Rio Grande do Sul State, Brazil, was analyzed using 1989-2003 at 00:00 and 12:00 UTC upper-air observations. The LLJ classification criteria proposed by Bonner (1968) and modified by Whiteman et aI. (1997) were applied to determine the LLJ occurrence. Afterwards was selected a LLJ event, that was one of the most intense observed in the summer (01/27/2002 at 12:00 UTC), during the study period. ln this study were used as tools: atmospheric soundings, GOES-8 satellite images, and wind, temperature and specific humidity fields from GLOBAL, ETA and BRAMS models. Based on the numerical analysis was possible to verify that the three models overestimated the specific humidity and potential temperature values, at LLJ time occurrence. The wind speed was underestimated by the models. It was observed in the study region, at 12:00 UTC (LLJ detected hour in the Porto Alegre region), by three models, warm and wet air from north, generating conditions to Mesoscale Convective System (MCS) formation and intensification.