Understanding the rapid summer warming and changes in temperature extremes since the mid-1990s over Western Europe

Analysis of observations indicates that there was a rapid increase in summer (June-August, JJA) mean surface air temperature (SAT) since the mid-1990s over Western Europe. Accompanying this rapid warming are significant increases in summer mean daily maximum temperature, daily minimum temperature, annual hottest day temperature and warmest night temperature, and an increase in frequency of summer days and tropical nights, while the change in the diurnal temperature range (DTR) is small. This study focuses on understanding causes of the rapid summer warming and associated temperature extreme changes. A set of experiments using the atmospheric component of the state-of-the-art HadGEM3 global climate model have been carried out to quantify relative roles of changes in sea surface temperature (SST)/sea ice extent (SIE), anthropogenic greenhouse gases (GHGs), and anthropogenic aerosols (AAer). Results indicate that the model forced by changes in all forcings reproduces many of the observed changes since the mid-1990s over Western Europe. Changes in SST/SIE explain 62.2% ± 13.0% of the area averaged seasonal mean warming signal over Western Europe, with the remaining 37.8% ± 13.6% of the warming explained by the direct impact of changes in GHGs and AAer. Results further indicate that the direct impact of the reduction of AAer precursor emissions over Europe, mainly through aerosol-radiation interaction with additional contributions from aerosol-cloud interaction and coupled atmosphere-land surface feedbacks, is a key factor for increases in annual hottest day temperature and in frequency of summer days. It explains 45.5% ± 17.6% and 40.9% ± 18.4% of area averaged signals for these temperature extremes. The direct impact of the reduction of AAer precursor emissions over Europe acts to increase DTR locally, but the change in DTR is countered by the direct impact of GHGs forcing. In the next few decades, greenhouse gas concentrations will continue to rise and AAer precursor emissions over Europe and North America will continue to decline. Our results suggest that the changes in summer seasonal mean SAT and temperature extremes over Western Europe since the mid-1990s are most likely to be sustained or amplified in the near term, unless other factors intervene.

Response surface models to predict broiler performance and applications for economic analysis

A study was carried out to elaborate response surface models using broiler performance data recovered from literature in order to predict performance and elaborate economic analyses. Nineteen studies published between 1995 and 2005 were retrieved using the systematic literature review method. Weight gain and feed conversion data were collected from eight studies that fulfilled the pre-established inclusion criteria, and a response surface model was adjusted using crude protein, environmental temperature, and age as independent variables. The models produced for weight gain (r² = 0.93) and feed conversion (r² = 0.85) were accurate, precise, and not biased. Protein levels, environmental temperature and age showed linear and quadratic effects on weight gain and feed conversion. There was no interaction between protein level and environmental temperature. Age and crude protein showed interaction for weight gain and feed conversion, whereas interaction between age and temperature was detected only for weight gain. It was possible to perform economic analyses to determine maximum profit as a function of the variables that were included in the model. It was concluded that the response surface models are effective to predict the performance of broiler chickens and allow the elaboration of economic analyses to optimize profit.

Semi-automatic modeling of buildings from digital surface models

Semi-automatic building detection and extraction is a topic of growing interest due to its potential application in such areas as cadastral information systems, cartographic revision, and GIS. One of the existing strategies for building extraction is to use a digital surface model (DSM) represented by a cloud of known points on a visible surface, and comprising features such as trees or buildings. Conventional surface modeling using stereo-matching techniques has its drawbacks, the most obvious being the effect of building height on perspective, shadows, and occlusions. The laser scanner, a recently developed technological tool, can collect accurate DSMs with high spatial frequency. This paper presents a methodology for semi-automatic modeling of buildings which combines a region-growing algorithm with line-detection methods applied over the DSM.

Imagens multitemporais ALOS/PALSAR baseada em orientação a objeto na discriminação da cobertura da terra

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Inverse Modeling of GPS Multipath for Snow Depth Estimation-Part I: Formulation and Simulations

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Métrica para definição de pontos estáveis utilizando-se da Análise de Robustez no posicionamento por ponto preciso

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES)

Radiation balance at the surface in the city of So Paulo, Brazil: diurnal and seasonal variations

The main goal of this work is to describe the diurnal and seasonal variations of the radiation balance components at the surface in the city of So Paulo based on observations carried out during 2004. Monthly average hourly values indicate that the amplitudes of the diurnal cycles of net radiation (Q*), downwelling and upwelling shortwave radiation (SW(DW), SW(UP)), and longwave radiations (LW(DW), LW(UP)) in February were, respectively, 37%, 14%, 19%, 11%, and 5% larger than they were in August. The monthly average daily values indicate a variation of 60% for Q*, with a minimum in June and a maximum in December; 45% for SW(DW), with a minimum in May and a maximum in September; 50% for SW(UP), with a minimum in June and a maximum in September; 13% for LW(DW), with a minimum in July and a maximum in January; and 9% for LW(UP), with a minimum in July and a maximum in February. It was verified that the atmospheric broadband transmissivity varied from 0.36 to 0.57; the effective albedo of the surface varied from 0.08 to 0.10; and the atmospheric effective emissivity varied from 0.79 to 0.92. The surface effective emissivity remained approximately constant and equal to 0.96. The albedo and surface effective emissivity for So Paulo agreed with those reported for urban areas in Europe and North America cities. This indicates that material and geometric effects on albedo and surface emissivity in So Paulo are similar to ones observed in typical middle latitudes cities. On the other hand, it was found that So Paulo city induces an urban heat island with daytime maximum intensity varying from 2.6A degrees C in July (16:00 LT) to 5.5A degrees C in September (15:00 LT). The analysis of the radiometric properties carried out here indicate that this daytime maximum is a primary response to the seasonal variation of daily values of net solar radiation at the surface.

A MODIS-Based Energy Balance to Estimate Evapotranspiration for Clear-Sky Days in Brazilian Tropical Savannas

Evapotranspiration (ET) plays an important role in global climate dynamics and in primary production of terrestrial ecosystems; it represents the mass and energy transfer from the land to atmosphere. Limitations to measuring ET at large scales using ground-based methods have motivated the development of satellite remote sensing techniques. The purpose of this work is to evaluate the accuracy of the SEBAL algorithm for estimating surface turbulent heat fluxes at regional scale, using 28 images from MODIS. SEBAL estimates are compared with eddy-covariance (EC) measurements and results from the hydrological model MGB-IPH. SEBAL instantaneous estimates of latent heat flux (LE) yielded r(2) = 0.64 and r(2) = 0.62 over sugarcane croplands and savannas when compared against in situ EC estimates. At the same sites, daily aggregated estimates of LE were r(2) = 0.76 and r(2) = 0.66, respectively. Energy balance closure showed that turbulent fluxes over sugarcane croplands were underestimated by 7% and 9% over savannas. Average daily ET from SEBAL is in close agreement with estimates from the hydrological model for an overlay of 38,100 km(2) (r(2) = 0.88). Inputs to which the algorithm is most sensitive are vegetation index (NDVI), gradient of temperature (dT) to compute sensible heat flux (H) and net radiation (Re). It was verified that SEBAL has a tendency to overestimate results both at local and regional scales probably because of low sensitivity to soil moisture and water stress. Nevertheless the results confirm the potential of the SEBAL algorithm, when used with MODIS images for estimating instantaneous LE and daily ET from large areas.

A GIS based model for solving the planar Huff problem considering different demand distributions and forbidden regions

[EN] In this paper, we have used Geographical Information Systems (GIS) to solve the planar Huff problem considering different demand distributions and forbidden regions. Most of the papers connected with the competitive location problems consider that the demand is aggregated in a finite set of points. In other few cases, the models suppose that the demand is distributed along the feasible region according to a functional form, mainly a uniform distribution. In this case, in addition to the discrete and uniform demand distributions we have considered that the demand is represented by a population surface model, that is, a raster map where each pixel has associated a value corresponding to the population living in the area that it covers...

Experimental and numerical analyses about the efficiency of flow through devices for the sediment controll in urban runoff

As land is developed, the impervious surfaces that are created increase the amount of runoff during rainfall events, disrupting the natural hydrologic cycle, with an increment in volume of runoff and in pollutant loadings. Pollutants deposited or derived from an activity on the land surface will likely end up in stormwater runoff in some concentration, such as nutrients, sediment, heavy metals, hydrocarbons, gasoline additives, pathogens, deicers, herbicides and pesticides. Several of these pollutants are particulate-bound, so it appears clear that sediment removal can provide significant water-quality improvements and it appears to be important the knowledge of the ability of stromwater treatment devices to retain particulate matter. For this reason three different units which remove sediments have been tested through laboratory. In particular a roadside gully pot has been tested under steady hydraulic conditions, varying the characteristics of the influent solids (diameter, particle size distribution and specific gravity). The efficiency in terms of particles retained has been evaluated as a function of influent flow rate and particles characteristics; results have been compared to efficiency evaluated applying an overflow rate model. Furthermore the role of particles settling velocity in efficiency determination has been investigated. After the experimental runs on the gully pot, a standard full-scale model of an hydrodynamic separator (HS) has been tested under unsteady influent flow rate condition, and constant solid concentration at the input. The results presented in this study illustrate that particle separation efficiency of the unit is predominately influenced by operating flow rate, which strongly affects the particles and hydraulic residence time of the system. The efficiency data have been compared to results obtained from a modified overflow rate model; moreover the residence time distribution has been experimentally determined through tracer analyses for several steady flow rates. Finally three testing experiments have been performed for two different configurations of a full-scale model of a clarifier (linear and crenulated) under unsteady influent flow rate condition, and constant solid concentration at the input. The results illustrate that particle separation efficiency of the unit is predominately influenced by the configuration of the unit itself. Turbidity measures have been used to compare turbidity with the suspended sediments concentration, in order to find a correlation between these two values, which can allow to have a measure of the sediments concentration simply installing a turbidity probe.

Interactions between the hydrological cycle and atmospheric chemistry

The land-atmosphere exchange of atmospheric trace gases is sensitive to meteorological conditions and climate change. It contributes in turn to the atmospheric radiative forcing through its effects on tropospheric chemistry. The interactions between the hydrological cycle and atmospheric processes are intricate and often involve different levels of feedbacks. The Earth system model EMAC is used in this thesis to assess the direct role of the land surface components of the terrestrial hydrological cycle in the emissions, deposition and transport of key trace gases that control tropospheric chemistry. It is also used to examine its indirect role in changing the tropospheric chemical composition through the feedbacks between the atmospheric and the terrestrial branches of the hydrological cycle. Selected features of the hydrological cycle in EMAC are evaluated using observations from different data sources. The interactions between precipitation and the water vapor column, from the atmospheric branch of the hydrological cycle, and evapotranspiration, from its terrestrial branch, are assessed specially for tropical regions. The impacts of changes in the land surface hydrology on surface exchanges and the oxidizing chemistry of the atmosphere are assessed through two sensitivity simulations. In the first, a new parametrization for rainfall interception in the densely vegetated areas in the tropics is implemented, and its effects are assessed. The second study involves the application of a soil moisture forcing that replaces the model calculated soil moisture. Both experiments have a large impact on the local hydrological cycle, dry deposition of soluble and insoluble gases, emissions of isoprene through changes in surface temperature and the Planetary Boundary Layer height. Additionally the soil moisture forcing causes changes in local vertical transport and large-scale circulation. The changes in trace gas exchanges affect the oxidation capacity of the atmosphere through changes in OH, O$_3$, NO$_x$ concentrations.

Soil biogenic emissions of nitric oxide from arid and semi-arid ecosystems

The biogenic production of NO in the soil accounts for between 10% and 40% of the global total. A large degree of the uncertainty in the estimation of the biogenic emissions stems from a shortage of measurements in arid regions, which comprise 40% of the earth’s land surface area. This study examined the emission of NO from three ecosystems in southern Africa which cover an aridity gradient from semi-arid savannas in South Africa to the hyper-arid Namib Desert in Namibia. A laboratory method was used to determine the release of NO as a function of the soil moisture and the soil temperature. Various methods were used to up-scale the net potential NO emissions determined in the laboratory to the vegetation patch, landscape or regional level. The importance of landscape, vegetation and climatic characteristics is emphasized. The first study occurred in a semi-arid savanna region in South Africa, where soils were sampled from 4 landscape positions in the Kruger National Park. The maximum NO emission occurred at soil moisture contents of 10%-20% water filled pore space (WFPS). The highest net potential NO emissions came from the low lying landscape positions, which have the largest nitrogen (N) stocks and the largest input of N. Net potential NO fluxes obtained in the laboratory were converted in field fluxes for the period 2003-2005, for the four landscape positions, using soil moisture and temperature data obtained in situ at the Kruger National Park Flux Tower Site. The NO emissions ranged from 1.5-8.5 kg ha-1 a-1. The field fluxes were up-scaled to a regional basis using geographic information system (GIS) based techniques, this indicated that the highest NO emissions occurred from the Midslope positions due to their large geographical extent in the research area. Total emissions ranged from 20x103 kg in 2004 to 34x103 kg in 2003 for the 56000 ha Skukuza land type. The second study occurred in an arid savanna ecosystem in the Kalahari, Botswana. In this study I collected soils from four differing vegetation patch types including: Pan, Annual Grassland, Perennial Grassland and Bush Encroached patches. The maximum net potential NO fluxes ranged from 0.27 ng m-2 s-1 in the Pan patches to 2.95 ng m-2 s-1 in the Perennial Grassland patches. The net potential NO emissions were up-scaled for the year December 2005-November 2006. This was done using 1) the net potential NO emissions determined in the laboratory, 2) the vegetation patch distribution obtained from LANDSAT NDVI measurements 3) estimated soil moisture contents obtained from ENVISAT ASAR measurements and 4) soil surface temperature measurements using MODIS 8 day land surface temperature measurements. This up-scaling procedure gave NO fluxes which ranged from 1.8 g ha-1 month-1 in the winter months (June and July) to 323 g ha-1 month-1 in the summer months (January-March). Differences occurred between the vegetation patches where the highest NO fluxes occurred in the Perennial Grassland patches and the lowest in the Pan patches. Over the course of the year the mean up-scaled NO emission for the studied region was 0.54 kg ha-1 a-1 and accounts for a loss of approximately 7.4% of the estimated N input to the region. The third study occurred in the hyper-arid Namib Desert in Namibia. Soils were sampled from three ecosystems; Dunes, Gravel Plains and the Riparian zone of the Kuiseb River. The net potential NO flux measured in the laboratory was used to estimate the NO flux for the Namib Desert for 2006 using modelled soil moisture and temperature data from the European Centre for Medium Range Weather Forecasts (ECMWF) operational model on a 36km x 35km spatial resolution. The maximum net potential NO production occurred at low soil moisture contents (<10%WFPS) and the optimal temperature was 25°C in the Dune and Riparian ecosystems and 35°C in the Gravel Plain Ecosystems. The maximum net potential NO fluxes ranged from 3.0 ng m-2 s-1 in the Riparian ecosystem to 6.2 ng m-2 s-1 in the Gravel Plains ecosystem. Up-scaling the net potential NO flux gave NO fluxes of up to 0.062 kg ha-1 a-1 in the Dune ecosystem and 0.544 kg h-1 a-1 in the Gravel Plain ecosystem. From these studies it is shown that NO is emitted ubiquitously from terrestrial ecosystems, as such the NO emission potential from deserts and scrublands should be taken into account in the global NO models. The emission of NO is influenced by various factors such as landscape, vegetation and climate. This study looks at the potential emissions from certain arid and semi-arid environments in southern Africa and other parts of the world and discusses some of the important factors controlling the emission of NO from the soil.

Global modeling of pollutant transport and deposition from anthropogenic emission hot spots on global scale

Urban centers significantly contribute to anthropogenic air pollution, although they cover only a minor fraction of the Earth's land surface. Since the worldwide degree of urbanization is steadily increasing, the anthropogenic contribution to air pollution from urban centers is expected to become more substantial in future air quality assessments. The main objective of this thesis was to obtain a more profound insight in the dispersion and the deposition of aerosol particles from 46 individual major population centers (MPCs) as well as the regional and global influence on the atmospheric distribution of several aerosol types. For the first time, this was assessed in one model framework, for which the global model EMAC was applied with different representations of aerosol particles. First, in an approach with passive tracers and a setup in which the results depend only on the source location and the size and the solubility of the tracers, several metrics and a regional climate classification were used to quantify the major outflow pathways, both vertically and horizontally, and to compare the balance between pollution export away from and pollution build-up around the source points. Then in a more comprehensive approach, the anthropogenic emissions of key trace species were changed at the MPC locations to determine the cumulative impact of the MPC emissions on the atmospheric aerosol burdens of black carbon, particulate organic matter, sulfate, and nitrate. Ten different mono-modal passive aerosol tracers were continuously released at the same constant rate at each emission point. The results clearly showed that on average about five times more mass is advected quasi-horizontally at low levels than exported into the upper troposphere. The strength of the low-level export is mainly determined by the location of the source, while the vertical transport is mainly governed by the lifting potential and the solubility of the tracers. Similar to insoluble gas phase tracers, the low-level export of aerosol tracers is strongest at middle and high latitudes, while the regions of strongest vertical export differ between aerosol (temperate winter dry) and gas phase (tropics) tracers. The emitted mass fraction that is kept around MPCs is largest in regions where aerosol tracers have short lifetimes; this mass is also critical for assessing the impact on humans. However, the number of people who live in a strongly polluted region around urban centers depends more on the population density than on the size of the area which is affected by strong air pollution. Another major result was that fine aerosol particles (diameters smaller than 2.5 micrometer) from MPCs undergo substantial long-range transport, with about half of the emitted mass being deposited beyond 1000 km away from the source. In contrast to this diluted remote deposition, there are areas around the MPCs which experience high deposition rates, especially in regions which are frequently affected by heavy precipitation or are situated in poorly ventilated locations. Moreover, most MPC aerosol emissions are removed over land surfaces. In particular, forests experience more deposition from MPC pollutants than other land ecosystems. In addition, it was found that the generic treatment of aerosols has no substantial influence on the major conclusions drawn in this thesis. Moreover, in the more comprehensive approach, it was found that emissions of black carbon, particulate organic matter, sulfur dioxide, and nitrogen oxides from MPCs influence the atmospheric burden of various aerosol types very differently, with impacts generally being larger for secondary species, sulfate and nitrate, than for primary species, black carbon and particulate organic matter. While the changes in the burdens of sulfate, black carbon, and particulate organic matter show an almost linear response for changes in the emission strength, the formation of nitrate was found to be contingent upon many more factors, e.g., the abundance of sulfuric acid, than only upon the strength of the nitrogen oxide emissions. The generic tracer experiments were further extended to conduct the first risk assessment to obtain the cumulative risk of contamination from multiple nuclear reactor accidents on the global scale. For this, many factors had to be taken into account: the probability of major accidents, the cumulative deposition field of the radionuclide cesium-137, and a threshold value that defines contamination. By collecting the necessary data and after accounting for uncertainties, it was found that the risk is highest in western Europe, the eastern US, and in Japan, where on average contamination by major accidents is expected about every 50 years.

Single fluoroscopic image based 2D/3D reconstruction of a scaled, patient-specific vertebral model

Purpose Accurate three-dimensional (3D) models of lumbar vertebrae can enable image-based 3D kinematic analysis. The common approach to derive 3D models is by direct segmentation of CT or MRI datasets. However, these have the disadvantages that they are expensive, timeconsuming and/or induce high-radiation doses to the patient. In this study, we present a technique to automatically reconstruct a scaled 3D lumbar vertebral model from a single two-dimensional (2D) lateral fluoroscopic image. Methods Our technique is based on a hybrid 2D/3D deformable registration strategy combining a landmark-to-ray registration with a statistical shape model-based 2D/3D reconstruction scheme. Fig. 1 shows different stages of the reconstruction process. Four cadaveric lumbar spine segments (total twelve lumbar vertebrae) were used to validate the technique. To evaluate the reconstruction accuracy, the surface models reconstructed from the lateral fluoroscopic images were compared to the associated ground truth data derived from a 3D CT-scan reconstruction technique. For each case, a surface-based matching was first used to recover the scale and the rigid transformation between the reconstructed surface model Results Our technique could successfully reconstruct 3D surface models of all twelve vertebrae. After recovering the scale and the rigid transformation between the reconstructed surface models and the ground truth models, the average error of the 2D/3D surface model reconstruction over the twelve lumbar vertebrae was found to be 1.0 mm. The errors of reconstructing surface models of all twelve vertebrae are shown in Fig. 2. It was found that the mean errors of the reconstructed surface models in comparison to their associated ground truths after iterative scaled rigid registrations ranged from 0.7 mm to 1.3 mm and the rootmean squared (RMS) errors ranged from 1.0 mm to 1.7 mm. The average mean reconstruction error was found to be 1.0 mm. Conclusion An accurate, scaled 3D reconstruction of the lumbar vertebra can be obtained from a single lateral fluoroscopic image using a statistical shape model based 2D/3D reconstruction technique. Future work will focus on applying the reconstructed model for 3D kinematic analysis of lumbar vertebrae, an extension of our previously-reported imagebased kinematic analysis. The developed method also has potential applications in surgical planning and navigation.

Validation of statistical shape model based reconstruction of the proximal femur - A morphology study

Seventeen bones (sixteen cadaveric bones and one plastic bone) were used to validate a method for reconstructing a surface model of the proximal femur from 2D X-ray radiographs and a statistical shape model that was constructed from thirty training surface models. Unlike previously introduced validation studies, where surface-based distance errors were used to evaluate the reconstruction accuracy, here we propose to use errors measured based on clinically relevant morphometric parameters. For this purpose, a program was developed to robustly extract those morphometric parameters from the thirty training surface models (training population), from the seventeen surface models reconstructed from X-ray radiographs, and from the seventeen ground truth surface models obtained either by a CT-scan reconstruction method or by a laser-scan reconstruction method. A statistical analysis was then performed to classify the seventeen test bones into two categories: normal cases and outliers. This classification step depends on the measured parameters of the particular test bone. In case all parameters of a test bone were covered by the training population's parameter ranges, this bone is classified as normal bone, otherwise as outlier bone. Our experimental results showed that statistically there was no significant difference between the morphometric parameters extracted from the reconstructed surface models of the normal cases and those extracted from the reconstructed surface models of the outliers. Therefore, our statistical shape model based reconstruction technique can be used to reconstruct not only the surface model of a normal bone but also that of an outlier bone.