The energy and water balance of a Eucalyptus plantation in southeast Brazil

The eddy covariance method was used to measure energy and water balance of a plantation of Eucalyptus (grandis x urophylla) hybrids over a 2 year period. The average daily evaporation rates were 5.4 (+/- 2.0) mm day(-1) in summer, but fell to 1.2 (+/- 0.3) mm day(-1) in winter. In contrast, the sensible heat flux was relatively low in summer but dominated the energy balance in winter. Evaporation accounted for 80% and 26% of the available energy, in summer and winter respectively. The annual evaporation was 82% (1124 mm) and 96% (1235 mm) of the annual rainfall recorded during the first and second year, respectively. Daily average canopy and aerodynamic conductance to water vapour were in the summer 51.9 (+/- 38.4) mm s(-1) 84.1 (+/- 25.6) mm s(-1), respectively; and in the winter 6.0 (+/- 10.5) mm s(-1) and 111.6 (+/- 24.6) mm s(-1), respectively. (C) 2010 Elsevier B.V. All rights reserved.

Atmosphere and hydrological controls of the evapotranspiration over a floodplain forest in the Bananal Island region, Amazonia

This article discusses seasonal and interannual variations of the evapotranspiration (ET) rates in Bananal Island floodplain, Brazil. Measurements included ET and sensible heat flux using the eddy covariance method, atmospheric forcings (net radiation, Rn, vapor pressure deficit, VPD, wind speed and air temperature), soil moisture profiles, groundwater level and flood height, taken from November 2003 to December 2006. For the hydrological years (October-September) of 2003/2004, 2004/2005 and 2005/2006, the accumulated precipitation was 1692, 1471, 1914 mm and the accumulated ET was 1361, 1318 and 1317 mm, respectively. Seasonal analyses indicated that ET decreased in the dry season (average 3.7 mm day(-1)), despite the simultaneous increase in Rn, air temperature and VPD. The increase of ET in the wet season and particularly in the flood period (average 4.1 mm day(-1)) showed that the free water surface evaporation strongly influenced the energy exchange. Soil moisture, which was substantially depleted during the dry season, and adaptative vegetation mechanisms such as leaf senescence contributed to limit the dry season ET. Strong drainage within permeable sandy soils helped to explain the soil moisture depletion. These results suggest that the Bananal flooding area shows a different pattern in relation to the upland Amazon forests, being more similar to the savanna strictu senso areas in central Brazil. For example, seasonal ET variation was not in phase with Rn; the wet season ET was higher than the dry season ET; and the system stored only a tiny memory of the flooding period, being sensitive to extended drought periods.

The land-atmosphere water flux in the tropics

Tropical vegetation is a major source of global land surface evapotranspiration, and can thus play a major role in global hydrological cycles and global atmospheric circulation. Accurate prediction of tropical evapotranspiration is critical to our understanding of these processes under changing climate. We examined the controls on evapotranspiration in tropical vegetation at 21 pan-tropical eddy covariance sites, conducted a comprehensive and systematic evaluation of 13 evapotranspiration models at these sites, and assessed the ability to scale up model estimates of evapotranspiration for the test region of Amazonia. Net radiation was the strongest determinant of evapotranspiration (mean evaporative fraction was 0.72) and explained 87% of the variance in monthly evapotranspiration across the sites. Vapor pressure deficit was the strongest residual predictor (14%), followed by normalized difference vegetation index (9%), precipitation (6%) and wind speed (4%). The radiation-based evapotranspiration models performed best overall for three reasons: (1) the vegetation was largely decoupled from atmospheric turbulent transfer (calculated from X decoupling factor), especially at the wetter sites; (2) the resistance-based models were hindered by difficulty in consistently characterizing canopy (and stomatal) resistance in the highly diverse vegetation; (3) the temperature-based models inadequately captured the variability in tropical evapotranspiration. We evaluated the potential to predict regional evapotranspiration for one test region: Amazonia. We estimated an Amazonia-wide evapotranspiration of 1370 mm yr(-1), but this value is dependent on assumptions about energy balance closure for the tropical eddy covariance sites; a lower value (1096 mm yr(-1)) is considered in discussion on the use of flux data to validate and interpolate models.

Aerosol properties, in-canopy gradients, turbulent fluxes and VOC concentrations at a pristine forest site in Amazonia

Aerosol physical and chemical properties were measured in a forest site in central Amazonia (Cuieiras reservation, 2.61S; 60.21W) during the dry season of 2004 (Aug-Oct). Aerosol light scattering and absorption, mass concentration, elemental composition and size distributions were measured at three tower levels (Ground: 2 m; Canopy: 28 m, and Top: 40 m). For the first time, simultaneous eddy covariance fluxes of fine mode particles and volatile organic compounds (VOC) were measured above the Amazonian forest canopy. Aerosol fluxes were measured by eddy covariance using a Condensation Particle Counter (CPC) and a sonic anemometer. VOC fluxes were measured by disjunct eddy covariance using a Proton Transfer Reaction Mass Spectrometer (PTR-MS). At nighttime, a strong vertical gradient of phosphorus and potassium in the aerosol coarse mode was observed, with higher concentrations at Ground level. This suggests a source of primary biogenic particles below the canopy. Equivalent black carbon measurements indicate the presence of light-absorbing aerosols from biogenic origin. Aerosol number size distributions typically consisted of superimposed Aitken (76 nm) and accumulation modes (144 nm), without clear events of new particle formation. Isoprene and monoterpene fluxes reached respectively 7.4 and 0.82 mg m(-2) s(-1) around noon. An average fine particle flux of 0.05 +/- 0.10 10(6) m(-2) s(-1) was calculated, denoting an equilibrium between emission and deposition fluxes of fine mode particles at daytime. No significant correlations were found between VOC and fine mode aerosol concentrations or fluxes. (C) 2009 Elsevier Ltd. All rights reserved.

Fluxes of CO(2) above a plantation of Eucalyptus in southeast Brazil

Eddy-covariance measurements of net ecosystem exchange of CO(2) (NEE) and estimates of gross ecosystem productivity (GEP) and ecosystem respiration (R(E)) were obtained in a 2-4 year old Eucalyptus plantation during two years with very different winter rainfall In the first (drier) year the annual NEE GEP and RE were lower than the sums in the second (normal) year and conversely the total respiratory costs of assimilated carbon were higher in the dry year than in the normal year Although the net primary production (NPP) in the first year was 23% lower than that of the second year the decrease in the carbon use efficiency (CUE = NPP/GEP) was 11% and autotrophic respiration utilized more resources in the first dry year than in the second normal year The time variations in NEE were followed by NPP because in these young Eucalyptus plantations NEE is very largely dominated by NPP and heterotrophic respiration plays only a relatively minor role During the dry season a pronounced hysteresis was observed in the relationship between NEE and photosynthetically active radiation and NEE fluxes were inversely proportional to humidity saturation deficit values greater than 0 8 kPa Nighttime fluxes of CO(2) during calm conditions when the friction velocity (u) was below the threshold (0 25 ms(-1)) were estimated based on a Q(10) temperature-dependence relationship adjusted separately for different classes of soil moisture content which regulated the temperature sensitivity of ecosystem respiration (C) 2010 Elsevier B V All rights reserved

Experimental and theoretical investigation of first hyperpolarizability in aminophenols

This work reports theoretical and experimental studies on the first hyperpolarizability (beta) of aminophenols, evaluating the influence of the NH(2) group position relative to the OH group on the hyperpolarizability. A new extension of hyper-Rayleigh scattering technique using picosecond pulse trains was employed to obtain the experimental absolute values of (beta). The theoretical static beta(0) values were calculated using AMI method implemented in the AMPAC program. The theoretical and experimental data show a clear dependence between beta and the relative position of the electron donor (D) and acceptor (A) groups, presenting the 2-aminophenol the higher values. Moreover, calculations show excellent qualitative agreement between theoretical and experimental data, which are improved when the simulations considering the solvated molecule in a combination of discrete solvent molecules interacting with the solute and the application of continuous dielectric model. Besides, the study indicates that the experimental hyperpolarizabilities seem to be a property of the solute-solvation shell system. These facts have affirmed that the theoretical approach employed can be successfully used to foresee the variation in beta due to modifications in the D/A position. Moreover, a theoretical study of the ground state absorption is performed and compared with experimental data. (C) 2008 Elsevier B.V. All rights reserved.

Reliability of a method for evaluating porosity in denture base resins

Background: The method of porosity analysis by water absorption has been carried out by the storage of the specimens in pure water, but it does not exclude the potential plasticising effect of the water generating unreal values of porosity. Objective: The present study evaluated the reliability of this method of porosity analysis in polymethylmethacrylate denture base resins by the determination of the most satisfactory solution for storage (S), where the plasticising effect was excluded. Materials and methods: Two specimen shapes (rectangular and maxillary denture base) and two denture base resins, water bath-polymerised (Classico) and microwave-polymerised (Acron MC) were used. Saturated anhydrous calcium chloride solutions (25%, 50%, 75%) and distilled water were used for specimen storage. Sorption isotherms were used to determine S. Porosity factor (PF) and diffusion coefficient (D) were calculated within S and for the groups stored in distilled water. anova and Tukey tests were performed to identify significant differences in PF results and Kruskal-Wallis test and Dunn multiple comparison post hoc test, for D results (alpha = 0.05). Results: For Acron MC denture base shape, FP results were 0.24% (S 50%) and 1.37% (distilled water); for rectangular shape FP was 0.35% (S 75%) and 0.19% (distilled water). For Classico denture base shape, FP results were 0.54% (S 75%) and 1.21% (distilled water); for rectangular shape FP was 0.7% (S 50%) and 1.32% (distilled water). FP results were similar in S and distilled water only for Acron MC rectangular shape (p > 0.05). D results in distilled water were statistically higher than S for all groups. Conclusions: The results of the study suggest that an adequate solution for storing specimens must be used to measure porosity by water absorption, based on excluding the plasticising effect.

A METHOD FOR THE STUDY OF ACCRETION DISK EMISSION IN CATACLYSMIC VARIABLES. I. THE MODEL

We have developed a spectrum synthesis method for modeling the ultraviolet (UV) emission from the accretion disk from cataclysmic variables (CVs). The disk is separated into concentric rings, with an internal structure from the Wade & Hubeny disk-atmosphere models. For each ring, a wind atmosphere is calculated in the comoving frame with a vertical velocity structure obtained from a solution of the Euler equation. Using simple assumptions, regarding rotation and the wind streamlines, these one-dimensional models are combined into a single 2.5-dimensional model for which we compute synthetic spectra. We find that the resulting line and continuum behavior as a function of the orbital inclination is consistent with the observations, and verify that the accretion rate affects the wind temperature, leading to corresponding trends in the intensity of UV lines. In general, we also find that the primary mass has a strong effect on the P Cygni absorption profiles, the synthetic emission line profiles are strongly sensitive to the wind temperature structure, and an increase in the mass-loss rate enhances the resonance line intensities. Synthetic spectra were compared with UV data for two high orbital inclination nova-like CVs-RW Tri and V347 Pup. We needed to include disk regions with arbitrary enhanced mass loss to reproduce reasonably well widths and line profiles. This fact and a lack of flux in some high ionization lines may be the signature of the presence of density-enhanced regions in the wind, or alternatively, may result from inadequacies in some of our simplifying assumptions.

An automatic correction method for the heel effect in digitized mammography images

The most significant radiation field nonuniformity is the well-known Heel effect. This nonuniform beam effect has a negative influence on the results of computer-aided diagnosis of mammograms, which is frequently used for early cancer detection. This paper presents a method to correct all pixels in the mammography image according to the excess or lack on radiation to which these have been submitted as a result of the this effect. The current simulation method calculates the intensities at all points of the image plane. In the simulated image, the percentage of radiation received by all the points takes the center of the field as reference. In the digitized mammography, the percentages of the optical density of all the pixels of the analyzed image are also calculated. The Heel effect causes a Gaussian distribution around the anode-cathode axis and a logarithmic distribution parallel to this axis. Those characteristic distributions are used to determine the center of the radiation field as well as the cathode-anode axis, allowing for the automatic determination of the correlation between these two sets of data. The measurements obtained with our proposed method differs on average by 2.49 mm in the direction perpendicular to the anode-cathode axis and 2.02 mm parallel to the anode-cathode axis of commercial equipment. The method eliminates around 94% of the Heel effect in the radiological image and the objects will reflect their x-ray absorption. To evaluate this method, experimental data was taken from known objects, but could also be done with clinical and digital images.

NMR Chemical Shielding and Spin-Spin Coupling Constants of Liquid NH(3): A Systematic Investigation using the Sequential QM/MM Method

The NMR spin coupling parameters, (1)J(N,H) and (2)J(H,H), and the chemical shielding, sigma((15)N), of liquid ammonia are studied from a combined and sequential QM/MM methodology. Monte Carlo simulations are performed to generate statistically uncorrelated configurations that are submitted to density functional theory calculations. Two different Lennard-Jones potentials are used in the liquid simulations. Electronic polarization is included in these two potentials via an iterative procedure with and without geometry relaxation, and the influence on the calculated properties are analyzed. B3LYP/aug-cc-pVTZ-J calculations were used to compute the V(N,H) constants in the interval of -67.8 to -63.9 Hz, depending on the theoretical model used. These can be compared with the experimental results of -61.6 Hz. For the (2)J(H,H) coupling the theoretical results vary between -10.6 to -13.01 Hz. The indirect experimental result derived from partially deuterated liquid is -11.1 Hz. Inclusion of explicit hydrogen bonded molecules gives a small but important contribution. The vapor-to-liquid shifts are also considered. This shift is calculated to be negligible for (1)J(N,H) in agreement with experiment. This is rationalized as a cancellation of the geometry relaxation and pure solvent effects. For the chemical shielding, U(15 N) Calculations at the B3LYP/aug-pcS-3 show that the vapor-to-liquid chemical shift requires the explicit use of solvent molecules. Considering only one ammonia molecule in an electrostatic embedding gives a wrong sign for the chemical shift that is corrected only with the use of explicit additional molecules. The best result calculated for the vapor to liquid chemical shift Delta sigma((15)N) is -25.2 ppm, in good agreement with the experimental value of -22.6 ppm.

Dating archeological ceramics from the Valley of Vitor, Arequipa by the TL method

The age of some ancient pottery from the Valley of Vitor in the region of Arequipa, Peru, is determined by the thermoluminescence (TL) method. For dating, a 325 degrees C TL peak was used and irradiation with -dose from 5 to 50Gy was carried out for the additive method, and from 0.4 to 5Gy for the regeneration method. For these dose values, the TL intensity is observed to grow linearly, obtaining an accumulated dose of 1.62 +/- 0.09Gy and 1.36 +/- 0.03Gy for the additive and regeneration methods, respectively. The age (A) of the sample was calculated by the two methods, being A=867 +/- 195 years after Christ (AC) for the additive method and A=1050 +/- 157 years AC for the regeneration method. Both results are within 800-1200 years AC, which is the period of the Wari culture.

A powder X-ray diffraction method for detection of polyprenylated benzophenones in plant extracts associated with HPLC for quantitative analysis

A robust, direct, rapid and non-destructive X-ray diffraction crystallography method to detect the polyprenylated benzophenones 7-epi-clusianone (1) and guttiferone A (2) in extracts from Garcinia brasiliensis is presented. Powder samples of benzophenones 1 and 2, dried hexane extracts from G. brasiliensis seeds and fruit`s pericarp, and the dried ethanolic extract from G. brasiliensis seeds were unambiguously characterized by powder X-ray diffractometry. The calculated X-ray diffraction peaks from crystal structures of analytes 1 and 2, previously determined by single-crystal X-ray diffraction technique, were overlaid to those of the experimental powder diffractograms, providing a practical identification of these compounds in the analyzed material and confirming the pure contents of the powder samples. Using the X-ray diffraction crystallography method, the studied polyprenylated benzophenones were selectively and simultaneously detected in the extracts which were mounted directly on sample holder. In addition, reference materials of the analytes were not required for analyses since the crystal structures of the compounds are known. High performance liquid chromatography analyses also were comparatively carried out to quantify the analytes in the same plant extracts showing to be in agreement with X-ray diffraction crystallography method. (C) 2010 Elsevier B.V. All rights reserved.

Efficient solutions to robust, semi-implicit discretizations of the immersed boundary method

The immersed boundary method is a versatile tool for the investigation of flow-structure interaction. In a large number of applications, the immersed boundaries or structures are very stiff and strong tangential forces on these interfaces induce a well-known, severe time-step restriction for explicit discretizations. This excessive stability constraint can be removed with fully implicit or suitable semi-implicit schemes but at a seemingly prohibitive computational cost. While economical alternatives have been proposed recently for some special cases, there is a practical need for a computationally efficient approach that can be applied more broadly. In this context, we revisit a robust semi-implicit discretization introduced by Peskin in the late 1970s which has received renewed attention recently. This discretization, in which the spreading and interpolation operators are lagged. leads to a linear system of equations for the inter-face configuration at the future time, when the interfacial force is linear. However, this linear system is large and dense and thus it is challenging to streamline its solution. Moreover, while the same linear system or one of similar structure could potentially be used in Newton-type iterations, nonlinear and highly stiff immersed structures pose additional challenges to iterative methods. In this work, we address these problems and propose cost-effective computational strategies for solving Peskin`s lagged-operators type of discretization. We do this by first constructing a sufficiently accurate approximation to the system`s matrix and we obtain a rigorous estimate for this approximation. This matrix is expeditiously computed by using a combination of pre-calculated values and interpolation. The availability of a matrix allows for more efficient matrix-vector products and facilitates the design of effective iterative schemes. We propose efficient iterative approaches to deal with both linear and nonlinear interfacial forces and simple or complex immersed structures with tethered or untethered points. One of these iterative approaches employs a splitting in which we first solve a linear problem for the interfacial force and then we use a nonlinear iteration to find the interface configuration corresponding to this force. We demonstrate that the proposed approach is several orders of magnitude more efficient than the standard explicit method. In addition to considering the standard elliptical drop test case, we show both the robustness and efficacy of the proposed methodology with a 2D model of a heart valve. (C) 2009 Elsevier Inc. All rights reserved.

Calculating Flash Point Numbers from Molecular Structure: An Improved Method for Predicting the Flash Points of Acyclic Alkanes

We report a novel method for calculating flash points of acyclic alkanes from flash point numbers, N(FP), which can be calculated from experimental or calculated boiling point numbers (Y(BP)) with the equation N(FP) = 1.020Y(BP) - 1.083 Flash points (FP) are then determined from the relationship FP(K) = 23.369N(FP)(2/3) + 20.010N(FP)(1/3) + 31.901 For it data set of 102 linear and branched alkanes, the correlation of literature and predicted flash points has R(2) = 0.985 and an average absolute deviation of 3.38 K. N(FP) values can also be estimated directly from molecular structure to produce an even closer correspondence of literature and predicted FP values. Furthermore, N(FP) values provide a new method to evaluate the reliability of literature flash point data.

Simple Method to Evaluate and to Predict Flash Points of Organic Compounds

Flash points (T(FP)) of organic compounds are calculated from their flash point numbers, N(FP), with the relationship T(FP) = 23.369N(FP)(2/3) + 20.010N(FP)(1/3) + 31.901. In turn, the N(FP) values can be predicted from boiling point numbers (Y(BP)) and functional group counts with the equation N(FP) = 0.974Y(BP) + Sigma(i)n(i)G(i) + 0.095 where G(i) is a functional group-specific contribution to the value of N(FP) and n(i) is the number of such functional groups in the structure. For a data set consisting of 1000 diverse organic compounds, the average absolute deviation between reported and predicted flash points was less than 2.5 K.