Economic shocks and crime: evidence from the crash of ponzi schemes

In November 2008, Colombian authorities dismantled a network of Ponzi schemes, making hundreds of thousands of investors lose tens of millions of dollars throughout the country. Using original data on the geographical incidence of the Ponzi schemes, this paper estimates the impact of their break down on crime. We find that the crash of Ponzi schemes differentially exacerbated crime in affected districts. Confirming the intuition of the standard economic model of crime, this effect is only present in places with relatively weak judicial and law enforcement institutions, and with little access to consumption smoothing mechanisms such as microcredit. In addition, we show that, with the exception of economically-motivated felonies such as robbery, violent crime is not affected by the negative shock.

Construcción y validación de un instrumento diagnóstico para la representación cognitiva de la Orientación Psicopedagógica en Galicia. 'Construction and validation of a diagnostic instrument for the cognitive representation of the Guidance in Galicia'.

Resumen tomado de la publicación. Con el apoyo económico del departamento MIDE de la UNED

Numerical model of solid phase transformations governed by nucleation and growth: microstructure development during isothermal crystallization

A simple numerical model which calculates the kinetics of crystallization involving randomly distributed nucleation and isotropic growth is presented. The model can be applied to different thermal histories and no restrictions are imposed on the time and the temperature dependences of the nucleation and growth rates. We also develop an algorithm which evaluates the corresponding emerging grain-size distribution. The algorithm is easy to implement and particularly flexible, making it possible to simulate several experimental conditions. Its simplicity and minimal computer requirements allow high accuracy for two- and three-dimensional growth simulations. The algorithm is applied to explore the grain morphology development during isothermal treatments for several nucleation regimes. In particular, thermal nucleation, preexisting nuclei, and the combination of both nucleation mechanisms are analyzed. For the first two cases, the universal grain-size distribution is obtained. The high accuracy of the model is stated from its comparison to analytical predictions. Finally, the validity of the Kolmogorov-Johnson-Mehl-Avrami model SSSR, is verified for all the cases studied

A Validated Methodology for Genetic Identification of Tuna Species (Genus Thunnus)

Tuna species of the genus Thunnus, such as the bluefin tunas, are some of the most important and yet most endangered trade fish in the world. Identification of these species in traded forms, however, may be difficult depending on the presentation of the products, which may hamper conservation efforts on trade control. In this paper, we validated a genetic methodology that can fully distinguish between the eight Thunnus species from any kind of processed tissue. Methodology: After testing several genetic markers, a complete discrimination of the eight tuna species was achieved using Forensically Informative Nucleotide Sequencing based primarily on the sequence variability of the hypervariable genetic marker mitochondrial DNA control region (mtDNA CR), followed, in some specific cases, by a second validation by a nuclear marker rDNA first internal transcribed spacer (ITS1). This methodology was able to distinguish all tuna species, including those belonging to the subgenus Neothunnus that are very closely related, and in consequence can not be differentiated with other genetic markers of lower variability. This methodology also took into consideration the presence of introgression that has been reported in past studies between T. thynnus, T. orientalis and T. alalunga. Finally, we applied the methodology to cross-check the species identity of 26 processed tuna samples. Conclusions: Using the combination of two genetic markers, one mitochondrial and another nuclear, allows a full discrimination between all eight tuna species. Unexpectedly, the genetic marker traditionally used for DNA barcoding, cytochrome oxidase 1, could not differentiate all species, thus its use as a genetic marker for tuna species identification is questioned

Empowerment: model i instrument de mesura basat en les variables primàries o efectives PIEQ (Primary Items Empowerment Questionnaire)

In this thesis a review of the state of the art on empowerment is found. A proposal is put forward for a model of the effective variables of empowerment, a new theoretical model that provides a distinct classification of individual variables. The model consists of four metavariables called reciprocal, unidirectional, shared and reflexive. Finally, a measurement tool for measuring and representing these variables graphically is described, including its process of elaboration through a preliminary pilot study that served to refine and improve the questionnaire. Both the model and the instrument proposed aim to determine the state of the primary variables involved in determining the predisposition and potential of a working group towards empowerment. With the measurement tool both organisational, departmental, working group and/or individual data can be obtained. This instrument can help companies and organisations to discover the limitations existing within a working group and act in consequence.

La investigació sobre l'èxit i el fracàs escolar des la la perspectiva dels factors de risc. Implicacions per a la recerca i la pràctica educatives

La situació problemàtica de la qual partim en aquesta tesi és la constatació de l'existència d'unes dinàmiques escolars negatives -expressades amb males notes reiterades- difícils de modificar, que determinats infants inicien i desenvolupen al llarg de la seva escolarització i que els condueix a una situació de fracàs escolar. Les males notes són el senyal que alerta pares i educadors de la presència de problemes escolars en els alumnes i constitueixen la causa explícita que motiva la cerca de solucions. Sovint es busquen solucions fora de l'escola, la qual cosa suggereix que, malgrat els esforços realitzats, l'escola i el sistema educatiu han tingut dificultats per a ajudar els infants a millorar globalment els seus resultats escolars. D'una banda, l'escola troba dificultats per a identificar el mes aviat possible quins infants seran mes susceptibles de desenvolupar unes dinàmiques escolars negatives que els puguin conduir a una situació de fracàs escolar. D'altra banda, també hi ha dificultats per trobar i aplicar estratègies preventives d'intervenció educativa a l'aula, que resultin adequades per a prevenir el desenvolupament de dinàmiques escolars negatives en alguns infants. Partint de la situació problemàtica descrita, la finalitat de la tesi es obtenir informació teòrica , empírica sobre algunes variables que puguin resultar rellevants per a articular, des de l'aula escolar, intervencions educatives destinades a prevenir el desenvolupament de les dinàmiques escolars negatives. Des de la perspectiva de la prevenció, la rellevància de les variables hauria d'establir-se en funció de la seva utilitat per a: A) Identificar des de l'aula escolar situacions de més perill -de més risc- de desenvolupar aquestes dinàmiques negatives i, consegüentment, arribar a la situació de fracàs escolar . Això implica, per tant, que les variables han de permetre la identificació abans que la situació de fracàs escolar es produeixi . B) Intervenir educativament des de l'escola; per tant, cal que siguin variables sobre les quals l'escola pugui incidir. La modificabilitat de les variables ha de permetre que es puguin emprendre accions educatives, des de la mateixa aula escolar. Variables que resultin suficientment importants de cara a disminuir el perill o evitar el desenvolupament de dinàmiques negatives. Es a dir, que la seva modificació contribueixi a evitar que els infants arribin a la situació de fracàs escolar. Per assolir aquesta finalitat es realitzen un segut de passes en funció de les quals s'ha estructurat la tesi en dues parts: Un marc teòric i un estudi de casos. EI Marc teòric té dos objectius: 1. Definir la situació problemàtica. En el capítol primer del marc teòric de la tesi, s'exposen les dimensions d'aquesta situació problemàtica. La revisió bibliogràfica entorn del tema del fracàs escolar ens ajuda a emmarcar la qüestió de les dinàmiques escolars en el fenomen complex del fracàs i l'èxit escolar, i del rendiment. Aquestes aportacions teòriques juntament amb les aportacions de la recerca educativa en relació a l'estabilitat del rendiment al llarg dels cursos són la base per a definir la situació problemàtica. 2. Delimitar, des d'un punt de vista teòric, algunes variables rellevants per a la prevenció del desenvolupament de dinàmiques escolars negatives, conduents a la situació de fracàs escolar. La primera passa per a dur a terme aquesta delimitació teòrica, que es presenta al segon capítol, ha estat revisar les aportacions de les investigacions sobre variables que incideixen en el rendiment escolar, les quals s'han analitzat en funció de la seva rellevància per a la prevenció. Aquesta revisió ha permès constatar un seguit de problemes de caire terminològic, metodològic i sobre la repercussió d'aquesta recerca en la practica educativa que afecten directament la utilitat de les aportacions d'aquestes investigacions de cara a identificar variables rellevants per a la prevenció. De l'anàlisi dels resultats d'aquestes recerques es desprèn que: a) Hi ha moltes variables associades al rendiment escolar, algunes difícilment modificables mitjançant la intervenció educativa escolar. b) EI fet que una variable estigui associada al rendiment no implica que sigui rellevant per a la prevenció c) S'obté poca informació sobre variables que contribueixin a disminuir la probabilitat de fracàs escolar. En base a aquests resultats es constata que cal buscar una perspectiva d'anàlisi de les variables mes adequades a l'enfocament preventiu, perspectiva que ha d'orientar-se a l'estudi del risc de fracàs escolar. La segona passa que es duu a terme per arribar a una delimitació teòrica de les variables rellevants es l'estudi del concepte de risc i d'altres conceptes relacionats: signe de risc, marcador de risc, factor de risc, factor protector, població en risc, infant en situació de risc, així com la revisió d'estudis i recerques que s'han plantejat en aquesta línia. Aquest treball ha permès: a) Clarificar aquests conceptes i aplicar-los en l'àmbit educatiu, en referència al problema del fracàs escolar, estructurant un marc teòric en funció del qual plantejar una anàlisi de les variables associades al rendiment escolar. b) Concretar un model per a l'anàlisi, des de l'aula escolar, de les variables que incideixen en el risc de fracàs escolar. Aquest model, que es situa en la perspectiva de l'aula escolar i que pren en consideració el paper actiu que els estudiants tenen respecte al seu aprenentatge, consta de tres components: unes variables, les seves relacions i la funció que exerceixen en relació al risc de fracàs escolar. La conclusió a la qual s'arriba es que, des d'un punt de vista teòric, una variable serà rellevant per a articular intervencions educatives preventives des de l'escola si constitueix o bé un factor de risc, o be un factor protector o compensador del risc. La delimitació empírica de variables rellevants per a la prevenció del fracàs escolar -que constitueix l'objectiu de la segona part de la tesi- es duu a terme mitjançant un estudi de casos que es deriva del plantejament teòric elaborat. L'objectiu d'aquest estudi es identificar variables que han pogut constituir factors protectors en joves que es troben en situació de risc per circumstàncies sociofamiliars. Se seleccionen tres noies i dos nois que, malgrat trobar-se en situació de risc, han assolit un cert nivell d'èxit escolar. Basant-nos en el model s'ha recollit informació sobre característiques actitudinals dels estudiants, del seu procés d'autoaprenentatge i de l'ambient d'aprenentatge. S'utilitza un disseny qualitatiu d'estudi de casos, utilitzant entrevistes amb profunditat per a recollir informació, la qual s'analitza mitjançant tècniques d'anàlisi de continguts. L'estudi de cada un dels cinc casos i la seva posterior comparació ha permès identificar algunes variables que poden haver constituït factors protectors del risc de fracàs escolar. Entre elles podem citar 1. La consciència de la pròpia situació complexa i desfavorable que han viscut o estan vivint. 2. Tenir un projecte vital a mig o llarg termini, en el qual els estudis són concebuts com a una via per assolir-lo. 3. Ser autoresponsables dels aprenentatges. 4. Haver identificat models a seguir en altres persones. Del treball realitzat tant des d'una perspectiva teòrica com empírica i les conclusions a les quals s'ha arribat se'n desprenen implicacions per a la practica educativa, per a la recerca i per a la formació professional dels educadors i educadores socials. Pel que fa a les implicacions per a la practica educativa, es proposa el model com a base per a la identificació de situacions de risc i per al disseny d'intervencions educatives amb l'objectiu de prevenir el fracàs escolar. En aquest sentit, i en funció dels resultats obtinguts a l'estudi de casos, es proposen unes línies d'intervenció preventiva en casos de risc similars als que han estat objecte d'estudi, línies que poden prendre en consideració tant els/les mestres com els educadors/es. Pel que fa a les implicacions per a la investigació educativa, es deriven quatre línies de recerca: investigació sobre factors de risc, investigació sobre factors protectors, investigació sobre el potencial preventiu d'intervencions educatives dissenyades en base a factors de risc i factors protectors, l'investigació sobre com potenciar des de diferents àmbits (escola i vida quotidiana) els fadors protectors. Quant a les implicacions per a la formació professional dels educadors/es socials, els resultats de l'estudi de casos com a possibles aspectes a treballar per part dels educadors/es impliquen un treball en el qual aquests professionals han d'haver rebut formació sobre: 1. La relació educativa com a recurs per a la intervenció educativa professional. 2. La necessitat d'un treball coordinat interdisciplinari com a estratègia de treball professional. 3. EI coneixement de programes coherents i estratègies d'intervenció sobre factors protectors. 4. L'elaboració de programes educatius, de manera que els educadors/es puguin adaptar les intervencions a les necessitats educatives dels subjectes. 5. La intervenció educativa en famílies.

Identification of Putative Wintering Areas and Ecological Determinants of Population Dynamics of Common House-Martin (Delichon urbicum) and Common Swift (Apus apus) Breeding in Northern Italy

To identify the causes of population decline in migratory birds, researchers must determine the relative influence of environmental changes on population dynamics while the birds are on breeding grounds, wintering grounds, and en route between the two. This is problematic when the wintering areas of specific populations are unknown. Here, we first identified the putative wintering areas of Common House-Martin (Delichon urbicum) and Common Swift (Apus apus) populations breeding in northern Italy as those areas, within the wintering ranges of these species, where the winter Normalized Difference Vegetation Index (NDVI), which may affect winter survival, best predicted annual variation in population indices observed in the breeding grounds in 1992–2009. In these analyses, we controlled for the potentially confounding effects of rainfall in the breeding grounds during the previous year, which may affect reproductive success; the North Atlantic Oscillation Index (NAO), which may account for climatic conditions faced by birds during migration; and the linear and squared term of year, which account for nonlinear population trends. The areas thus identified ranged from Guinea to Nigeria for the Common House-Martin, and were located in southern Ghana for the Common Swift. We then regressed annual population indices on mean NDVI values in the putative wintering areas and on the other variables, and used Bayesian model averaging (BMA) and hierarchical partitioning (HP) of variance to assess their relative contribution to population dynamics. We re-ran all the analyses using NDVI values at different spatial scales, and consistently found that our population of Common House-Martin was primarily affected by spring rainfall (43%–47.7% explained variance) and NDVI (24%–26.9%), while the Common Swift population was primarily affected by the NDVI (22.7%–34.8%). Although these results must be further validated, currently they are the only hypotheses about the wintering grounds of the Italian populations of these species, as no Common House-Martin and Common Swift ringed in Italy have been recovered in their wintering ranges.

Spherical nonparametric estimators applied to the UGAMP model integration for AMIP

The aim of this paper is essentially twofold: first, to describe the use of spherical nonparametric estimators for determining statistical diagnostic fields from ensembles of feature tracks on a global domain, and second, to report the application of these techniques to data derived from a modern general circulation model. New spherical kernel functions are introduced that are more efficiently computed than the traditional exponential kernels. The data-driven techniques of cross-validation to determine the amount elf smoothing objectively, and adaptive smoothing to vary the smoothing locally, are also considered. Also introduced are techniques for combining seasonal statistical distributions to produce longer-term statistical distributions. Although all calculations are performed globally, only the results for the Northern Hemisphere winter (December, January, February) and Southern Hemisphere winter (June, July, August) cyclonic activity are presented, discussed, and compared with previous studies. Overall, results for the two hemispheric winters are in good agreement with previous studies, both for model-based studies and observational studies.

Data assimilation for a coastal area morphodynamic model: Morecambe Bay

This paper investigates the use of data assimilation in coastal area morphodynamic modelling using Morecambe Bay as a study site. A simple model of the bay has been enhanced with a data assimilation scheme to better predict large-scale changes in bathymetry observed in the bay over a 3-year period. The 2DH decoupled morphodynamic model developed for the work is described, as is the optimal interpolation scheme used to assimilate waterline observations into the model run. Each waterline was acquired from a SAR satellite image and is essentially a contour of the bathymetry at some level within the inter-tidal zone of the bay. For model parameters calibrated against validation observations, model performance is good, even without data assimilation. However the use of data assimilation successfully compensates for a particular failing of the model, and helps to keep the model bathymetry on track. It also improves the ability of the model to predict future bathymetry. Although the benefits of data assimilation are demonstrated using waterline observations, any observations of morphology could potentially be used. These results suggest that data assimilation should be considered for use in future coastal area morphodynamic models.

An integrated model of soil-canopy spectral radiances, photosynthesis, fluorescence, temperature and energy balance

This paper presents the model SCOPE (Soil Canopy Observation, Photochemistry and Energy fluxes), which is a vertical (1-D) integrated radiative transfer and energy balance model. The model links visible to thermal infrared radiance spectra (0.4 to 50 μm) as observed above the canopy to the fluxes of water, heat and carbon dioxide, as a function of vegetation structure, and the vertical profiles of temperature. Output of the model is the spectrum of outgoing radiation in the viewing direction and the turbulent heat fluxes, photosynthesis and chlorophyll fluorescence. A special routine is dedicated to the calculation of photosynthesis rate and chlorophyll fluorescence at the leaf level as a function of net radiation and leaf temperature. The fluorescence contributions from individual leaves are integrated over the canopy layer to calculate top-of-canopy fluorescence. The calculation of radiative transfer and the energy balance is fully integrated, allowing for feedback between leaf temperatures, leaf chlorophyll fluorescence and radiative fluxes. Leaf temperatures are calculated on the basis of energy balance closure. Model simulations were evaluated against observations reported in the literature and against data collected during field campaigns. These evaluations showed that SCOPE is able to reproduce realistic radiance spectra, directional radiance and energy balance fluxes. The model may be applied for the design of algorithms for the retrieval of evapotranspiration from optical and thermal earth observation data, for validation of existing methods to monitor vegetation functioning, to help interpret canopy fluorescence measurements, and to study the relationships between synoptic observations with diurnally integrated quantities. The model has been implemented in Matlab and has a modular design, thus allowing for great flexibility and scalability.

Evaluating the century C model using long-term fertilizer trials in the Indo-Gangetic Plains, India

The GEFSOC Project developed a system for estimating soil carbon (C) stocks and changes at the national and sub-national scale. As part of the development of the system, the Century ecosystem model was evaluated for its ability to simulate soil organic C (SOC) changes in environmental conditions in the Indo-Gangetic Plains, India (IGP). Two long-term fertilizer trials (LTFT), with all necessary parameters needed to run Century, were used for this purpose: a jute (Corchorus capsularis L.), rice (Oryza sativa L.) and wheat (Triticum aestivum L.) trial at Barrackpore, West Bengal, and a rice-wheat trial at Ludhiana, Punjab. The trials represent two contrasting climates of the IGP, viz. semi-arid, dry with mean annual rainfall (MAR) of < 800 mm and humid with > 1600 turn. Both trials involved several different treatments with different organic and inorganic fertilizer inputs. In general, the model tended to overestimate treatment effects by approximately 15%. At the semi-arid site, modelled data simulated actual data reasonably well for all treatments, with the control and chemical N + farm yard manure showing the best agreement (RMSE = 7). At the humid site, Century performed less well. This could have been due to a range of factors including site history. During the study, Century was calibrated to simulate crop yields for the two sites considered using data from across the Indian IGP. However, further adjustments may improve model performance at these sites and others in the IGP. The availability of more longterm experimental data sets (especially those involving flooded lowland rice and triple cropping systems from the IGP) for testing and validation is critical to the application of the model's predictive capabilities for this area of the Indian sub-continent. (C) 2007 Elsevier B.V. All rights reserved.

The Phosphorus Indicators Tool: a simple model of diffuse P loss from agricultural land to water

The Phosphorus Indicators Tool provides a catchment-scale estimation of diffuse phosphorus (P) loss from agricultural land to surface waters using the most appropriate indicators of P loss. The Tool provides a framework that may be applied across the UK to estimate P loss, which is sensitive not only to land use and management but also to environmental factors such as climate, soil type and topography. The model complexity incorporated in the P Indicators Tool has been adapted to the level of detail in the available data and the need to reflect the impact of changes in agriculture. Currently, the Tool runs on an annual timestep and at a 1 km(2) grid scale. We demonstrate that the P Indicators Tool works in principle and that its modular structure provides a means of accounting for P loss from one layer to the next, and ultimately to receiving waters. Trial runs of the Tool suggest that modelled P delivery to water approximates measured water quality records. The transparency of the structure of the P Indicators Tool means that identification of poorly performing coefficients is possible, and further refinements of the Tool can be made to ensure it is better calibrated and subsequently validated against empirical data, as it becomes available.

Characterization and identification of mixed-metal phosphates in soils: the application, of Raman spectroscopy

The development of protocols for the identification of metal phosphates in phosphate-treated, metal-contaminated soils is a necessary yet problematical step in the validation of remediation schemes involving immobilization of metals as phosphate phases. The potential for Raman spectroscopy to be applied to the identification of these phosphates in soils has yet to be fully explored. With this in mind, a range of synthetic mixed-metal hydroxylapatites has been characterized and added to soils at known concentrations for analysis using both bulk X-ray powder diffraction (XRD) and Raman spectroscopy. Mixed-metal hydroxylapatites in the binary series Ca-Cd, Ca-Pb, Ca-Sr and Cd-Pb synthesized in the presence of acetate and carbonate ions, were characterized using a range of analytical techniques including XRD, analytical scanning electron microscopy (SEM), infrared spectroscopy (IR), inductively coupled plasma-atomic emission spectrometry (ICP-AES) and Raman spectroscopy. Only the Ca-Cd series displays complete solid solution, although under the synthesis conditions of this study the Cd-5(PO4)(3)OH end member could not be synthesized as a pure phase. Within the Ca-Cd series the cell parameters, IR active modes and Raman active bands vary linearly as a function of Cd content. X-ray diffraction and extended X-ray absorption fine structure spectroscopy (EXAFS) suggest that the Cd is distributed across both the Ca(1) and Ca(2) sites, even at low Cd concentrations. In order to explore the likely detection limits for mixed-metal phosphates in soils for XRD and Raman spectroscopy, soils doped with mixed-metal hydroxylapatites at concentrations of 5, 1 and 0.5 wt.% were then studied. X-ray diffraction could not confirm unambiguously the presence or identity of mixed-metal phosphates in soils at concentrations below 5 wt.%. Raman spectroscopy proved a far more sensitive method for the identification of mixed-metal hydroxylapatites in soils, which could positively identify the presence of such phases in soils at all the dopant concentrations used in this study. Moreover, Raman spectroscopy could also provide an accurate assessment of the degree of chemical substitution in the hydroxylapatites even when present in soils at concentrations as low as 0.1%.

Spatial analysis of the error in a model of soil nitrogen

Models of the dynamics of nitrogen in soil (soil-N) can be used to aid the fertilizer management of a crop. The predictions of soil-N models can be validated by comparison with observed data. Validation generally involves calculating non-spatial statistics of the observations and predictions, such as their means, their mean squared-difference, and their correlation. However, when the model predictions are spatially distributed across a landscape the model requires validation with spatial statistics. There are three reasons for this: (i) the model may be more or less successful at reproducing the variance of the observations at different spatial scales; (ii) the correlation of the predictions with the observations may be different at different spatial scales; (iii) the spatial pattern of model error may be informative. In this study we used a model, parameterized with spatially variable input information about the soil, to predict the mineral-N content of soil in an arable field, and compared the results with observed data. We validated the performance of the N model spatially with a linear mixed model of the observations and model predictions, estimated by residual maximum likelihood. This novel approach allowed us to describe the joint variation of the observations and predictions as: (i) independent random variation that occurred at a fine spatial scale; (ii) correlated random variation that occurred at a coarse spatial scale; (iii) systematic variation associated with a spatial trend. The linear mixed model revealed that, in general, the performance of the N model changed depending on the spatial scale of interest. At the scales associated with random variation, the N model underestimated the variance of the observations, and the predictions were correlated poorly with the observations. At the scale of the trend, the predictions and observations shared a common surface. The spatial pattern of the error of the N model suggested that the observations were affected by the local soil condition, but this was not accounted for by the N model. In summary, the N model would be well-suited to field-scale management of soil nitrogen, but suited poorly to management at finer spatial scales. This information was not apparent with a non-spatial validation. (c),2007 Elsevier B.V. All rights reserved.

Integrated Nitrogen CAtchment model (INCA) applied to a tropical catchment in the Atlantic Forest, Sao Paulo, Brazil

Stream-water flows and in-stream nitrate and ammonium concentrations in a small (36.7 ha) Atlantic Forest catchment were simulated using the Integrated Nitrogen in CAtchments (INCA) model version 1.9.4. The catchment, at Cunha, is in the Serra do Mar State Park, SE Brazil and is nearly pristine because the nearest major conurbations, Sao Paulo and Rio, are some 450 km distant. However, intensive farming may increase nitrogen (N) deposition and there are growing pressures for urbanisation. The mean-monthly discharges and NO3-N concentration dynamics were simulated adequately for the calibration and validation periods with (simulated) loss rates of 6.55 kg.ha(-1) yr(-1) for NO3-N and 3.85 kg.ha(-1) yr(-1) for NH4-N. To investigate the effects of elevated levels of N deposition in the future, various scenarios for atmospheric deposition were simulated; the highest value corresponded to that in a highly polluted area of Atlantic Forest in Sao Paulo City. It was found that doubling the atmospheric deposition generated a 25% increase in the N leaching rate, while at levels approaching the highly polluted Sao Paulo deposition rate, five times higher than the current rate, leaching increased by 240%, which would create highly eutrophic conditions, detrimental to downstream water quality. The results indicate that the INCA model can be useful for estimating N concentration and fluxes for different atmospheric deposition rates and hydrological conditions.