Application of landsat images for quantifying the energy balance under conditions of land use changes in the semi-arid region of Brazil

In the Nilo Coelho irrigation scheme, Brazil, the natural vegetation has been replaced by irrigated agriculture, bringing importance for the quantification of the effects on the energy exchanges between the mixed vegetated surfaces and the lower atmosphere. Landsat satellite images and agro-meteorological stations from 1992 to 2011 were used together, for modelling these exchanges. Surface albedo (α0), NDVI and surface temperature (T0) were the basic remote sensing retrieving parameters necessary to calculate the latent heat flux (λE) and the surface resistance to evapotranspiration (rs) on a large scale. The daily net radiation (Rn) was obtained from α0, air temperature (Ta) and short-wave transmissivity (τsw) throughout the slob equation, allowing the quantification of the daily sensible heat flux (H) by residual in the energy balance equation. With a threshold value for rs, it was possible to separate the energy fluxes from crops and natural vegetation. The averaged fractions of Rn partitioned as H and λE, were in average 39 and 67%, respectively. It was observed an increase of the energy used for the evapotranspiration process inside irrigated areas from 51% in 1992 to 80% in 2011, with the ratio λE/Rn presenting an increase of 3 % per year. The tools and models applied in the current research, can subsidize the monitoring of the coupled climate and land use changes effects in irrigation perimeters, being valuable when aiming the sustainability of the irrigated agriculture in the future, avoiding conflicts among different water users. © 2012 SPIE.

Balanço de energia de fluxo de CO2: cerrado da floresta estadual de Assis

Based on climatological data and energy flow, this paper analyzes the behavior observed in microclimatic an important remnant of cerrado in the São Paulo State. The seasonal climate is well marked, in the study area, with two climatic periods (one dry and one wet). The vegetation presents physiology as a function of water availability and the local thermal conditions and can be seen marked changes in the landscape due to the more or less presence of solar radiation. It’s important to understand the ecosystems behave in the context of global change. It has gained the attention of many researchers in the world.

Influence of salinity on the interannual heat storage trends in the Atlantic estimated from altimeters and Pilot Research Moored Array in the Tropical Atlantic data

Changes in the oceanic heat storage (HS) can reveal important evidences of climate variability related to ocean heat fluxes. Specifically, long-term variations in HS are a powerful indicator of climate change as HS represents the balance between the net surface energy flux and the poleward heat transported by the ocean currents. HS is estimated from sea surface height anomaly measured from the altimeters TOPEX/Poseidon and Jason 1 from 1993 to 2006. To characterize and validate the altimeter-based HS in the Atlantic, we used the data from the Pilot Research Moored Array in the Tropical Atlantic (PIRATA) array. Correlations and rms differences are used as statistical figures of merit to compare the HS estimates. The correlations range from 0.50 to 0.87 in the buoys located at the equator and at the southern part of the array. In that region the rms differences range between 0.40 and 0.51 x 10(9) Jm(-2). These results are encouraging and indicate that the altimeter has the precision necessary to capture the interannual trends in HS in the Atlantic. Albeit relatively small, salinity changes can also have an effect on the sea surface height anomaly. To account for this effect, NCEP/GODAS reanalysis data are used to estimate the haline contraction. To understand which dynamical processes are involved in the HS variability, the total signal is decomposed into nonpropagating basin-scale and seasonal (HS(l)) planetary waves, mesoscale eddies, and small-scale residual components. In general, HS(l) is the dominant signal in the tropical region. Results show a warming trend of HS(l) in the past 13 years almost all over the Atlantic basin with the most prominent slopes found at high latitudes. Positive interannual trends are found in the halosteric component at high latitudes of the South Atlantic and near the Labrador Sea. This could be an indication that the salinity anomaly increased in the upper layers during this period. The dynamics of the South Atlantic subtropical gyre could also be subject to low-frequency changes caused by a trend in the halosteric component on each side of the South Atlantic Current.

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.

Bio-physical interactions and feedbacks in a global climate model

This PhD thesis addresses the topic of large-scale interactions between climate and marine biogeochemistry. To this end, centennial simulations are performed under present and projected future climate conditions with a coupled ocean-atmosphere model containing a complex marine biogeochemistry model. The role of marine biogeochemistry in the climate system is first investigated. Phytoplankton solar radiation absorption in the upper ocean enhances sea surface temperatures and upper ocean stratification. The associated increase in ocean latent heat losses raises atmospheric temperatures and water vapor. Atmospheric circulation is modified at tropical and extratropical latitudes with impacts on precipitation, incoming solar radiation, and ocean circulation which cause upper-ocean heat content to decrease at tropical latitudes and to increase at middle latitudes. Marine biogeochemistry is tightly related to physical climate variability, which may vary in response to internal natural dynamics or to external forcing such as anthropogenic carbon emissions. Wind changes associated with the North Atlantic Oscillation (NAO), the dominant mode of climate variability in the North Atlantic, affect ocean properties by means of momentum, heat, and freshwater fluxes. Changes in upper ocean temperature and mixing impact the spatial structure and seasonality of North Atlantic phytoplankton through light and nutrient limitations. These changes affect the capability of the North Atlantic Ocean of absorbing atmospheric CO2 and of fixing it inside sinking particulate organic matter. Low-frequency NAO phases determine a delayed response of ocean circulation, temperature and salinity, which in turn affects stratification and marine biogeochemistry. In 20th and 21st century simulations natural wind fluctuations in the North Pacific, related to the two dominant modes of atmospheric variability, affect the spatial structure and the magnitude of the phytoplankton spring bloom through changes in upper-ocean temperature and mixing. The impacts of human-induced emissions in the 21st century are generally larger than natural climate fluctuations, with the phytoplankton spring bloom starting one month earlier than in the 20th century and with ~50% lower magnitude. This PhD thesis advances the knowledge of bio-physical interactions within the global climate, highlighting the intrinsic coupling between physical climate and biosphere, and providing a framework on which future studies of Earth System change can be built on.

Retrospective analysis of a nonforecasted rain-on-snow flood in the Alps – a matter of model limitations or unpredictable nature?

A rain-on-snow flood occurred in the Bernese Alps, Switzerland, on 10 October 2011, and caused significant damage. As the flood peak was unpredicted by the flood forecast system, questions were raised concerning the causes and the predictability of the event. Here, we aimed to reconstruct the anatomy of this rain-on-snow flood in the Lötschen Valley (160 km2) by analyzing meteorological data from the synoptic to the local scale and by reproducing the flood peak with the hydrological model WaSiM-ETH (Water Flow and Balance Simulation Model). This in order to gain process understanding and to evaluate the predictability. The atmospheric drivers of this rain-on-snow flood were (i) sustained snowfall followed by (ii) the passage of an atmospheric river bringing warm and moist air towards the Alps. As a result, intensive rainfall (average of 100 mm day-1) was accompanied by a temperature increase that shifted the 0° line from 1500 to 3200 m a.s.l. (meters above sea level) in 24 h with a maximum increase of 9 K in 9 h. The south-facing slope of the valley received significantly more precipitation than the north-facing slope, leading to flooding only in tributaries along the south-facing slope. We hypothesized that the reason for this very local rainfall distribution was a cavity circulation combined with a seeder-feeder-cloud system enhancing local rainfall and snowmelt along the south-facing slope. By applying and considerably recalibrating the standard hydrological model setup, we proved that both latent and sensible heat fluxes were needed to reconstruct the snow cover dynamic, and that locally high-precipitation sums (160 mm in 12 h) were required to produce the estimated flood peak. However, to reproduce the rapid runoff responses during the event, we conceptually represent likely lateral flow dynamics within the snow cover causing the model to react "oversensitively" to meltwater. Driving the optimized model with COSMO (Consortium for Small-scale Modeling)-2 forecast data, we still failed to simulate the flood because COSMO-2 forecast data underestimated both the local precipitation peak and the temperature increase. Thus we conclude that this rain-on-snow flood was, in general, predictable, but requires a special hydrological model setup and extensive and locally precise meteorological input data. Although, this data quality may not be achieved with forecast data, an additional model with a specific rain-on-snow configuration can provide useful information when rain-on-snow events are likely to occur.

Stable oxygen isotope record and lithogenic composition of ODP Hole 117-722B

The upper 38 m of Hole 722B sediments (Owen Ridge, northwest Arabian Sea) was sampled at 20 cm intervals and used to develop records of lithogenic percent, mass accumulation rate, and grain size spanning the past 1 m.y. Over this interval, the lithogenic component of Owen Ridge sediments can be used to infer variability in the strength of Arabian Sea summer monsoon winds (median grain size) and the aridity of surrounding dust source-areas (mass accumulation rate; MAR in g/cm**2/k.y). The lithogenic MAR has strong 100, 41, and 23 k.y. cyclicities and is forced primarily by changes in source-area aridity associated with glacial-interglacial cycles. The lithogenic grain size, on the other hand, exhibits higher frequency variability (23 k.y.) and is forced by the strength of summer monsoon winds which, in turn, are forced by the effective sensible heating of the Indian-Asian landmass and by the availability of latent heat from the Southern Hemisphere Indian Ocean. These forcing mechanisms combine to produce a wind-strength record which has no strong relationship to glacial-interglacial cycles. Discussion of the mechanisms responsible for production of primary Milankovitch cyclicities in lithogenic records from the Owen Ridge is presented elsewhere (Clemens and Prell, 1990, doi:10.1029/PA005i002p00109). Here we examine the 1 m.y. record from Hole 722B focusing on different aspects of the lithogenic components including an abrupt change in the monsoon wind-strength record at 500 k.y., core-to-core reproducibility, comparison with magnetic susceptibility, coherency with a wind-strength record from the Pacific Ocean, and combination frequencies in the wind-strength record. The Hole 722B lithogenic grain-size record shows an abrupt change at 500 k.y. possibly indicating decreased monsoon wind-strength over the interval from 500 k.y. to present. The grain-size decrease appears to be coincident with a loss of spectral power near the 41 k.y. periodicity. However, the grain-size decrease is not paralleled in the Globigerina bulloides upwelling record, an independent record of summer monsoon wind-strength (Prell, this volume). These observations leave us with competing hypotheses possibly involving: (1) a decrease in the sensitivity of monsoon windstrength to obliquity forcing, (2) decoupling of the grain size and G. bulloides records via a decoupling of the nutrient supply from wind-driven upwelling, and/or (3) a change in dust source-area or the patterns of dust transporting winds. Comparison of the lithogenic grain size and weight percent records from Hole 722B with those from a nearby core shows that the major and most minor events are well replicated. These close matches establish our confidence in the lithogenic extraction techniques and measurements. Further, reproducibility on a core-to-core scale indicates that the eolian depositional signal is regionally strong, coherent, and well preserved. The lithogenic weight percent and magnetic susceptibility are extremely well correlated in both the time and frequency domains. From this we infer that the magnetically susceptible component of Owen Ridge sediments is of terrestrial origin and transported to the Owen Ridge via summer monsoon winds. Because of the high correlation with the lithogenic percent record, the magnetic susceptibility record can be cast in terms of lithogenic MAR and used as a high resolution proxy for continental aridity. In addition to primary Milankovitch periodicities, the Hole 722B grain-size record exhibits periodicity at 52 k.y. and at 29 k.y. Both periodicities are also found in the grain-size record from piston core RC11-210 in the equatorial Pacific Ocean. Comparison of the two grain-size records shows significant coherence and zero phase relationships over both the 52 and 29 k.y. periodicities suggesting that the strengths of the Indian Ocean monsoon and the Pacific southeasterly trade winds share common forcing mechanisms. Two possible origins for the 52 and 29 k.y. periodicities in the Hole 722B wind-strength record are (1) direct Milankovitch forcing (54 and 29 k.y. components of obliquity) and (2) combination periodicities resulting from nonlinear interactions within the climate system. We find that the 52 and 29 k.y. periodicities show stronger coherency with crossproducts of eccentricity and obliquity (29 k.y.) and precession and obliquity (52 k.y.) than with direct obliquity forcing. Our working hypothesis attributes these periodicities to nonlinear interaction between external insolation forcing and internal climatic feedback mechanisms involving an interdependence of continental snow/ice-mass (albedo) and the hydrological cycle (latent heat availability).

Evapotranspiration of grapevine trained to a gable trellis system under netting and black plastic mulching

The evapotranspiration (ET c) of a table grape vineyard (Vitis vinifera, cv. Red Globe) trained to a gable trellis under netting and black plastic mulching was determined under semiarid conditions in the central Ebro River Valley during 2007 and 2008. The netting was made of high-density polyethylene (pores of 12 mm2) and was placed just above the ground canopy about 2.2 m above soil surface. Black plastic mulching was used to minimize soil evaporation. The surface renewal method was used to obtain values of sensible heat flux (H) from high-frequency temperature readings. Later, latent heat flux (LE) values were obtained by solving the energy balance equation. For the May–October period, seasonal ET c was about 843 mm in 2007 and 787 mm in 2008. The experimental weekly crop coefficients (K cexp) fluctuated between 0.64 and 1.2. These values represent crop coefficients adjusted to take into account the reduction in ET c caused by the netting and the black plastic mulching. Average K cexp values during mid- and end-season stages were 0.79 and 0.98, respectively. End-season K cexp was higher due to combination of factors related to the precipitation and low ET o conditions that are typical in this region during fall. Estimated crop coefficients using the Allen et al. (1998) approach adjusting for the effects of the netting and black plastic mulching (K cFAO) showed a good agreement with the experimental K cexp values.

Evapotranspiration and crop coefficients of rice (Oryza sativa L.) under sprinkler irrigation in a semiarid climate determined by the surface renewal method

The evapotranspiration (ETc) of sprinkler-irrigated rice was determined for the semiarid conditions of NE Spain during 2001, 2002 and 2003. The surface renewal method, after calibration against the eddy covariance method, was used to obtain values of sensible heat flux (H) from high-frequency temperature readings. Latent heat flux values were obtained by solving the energy balance equation. Finally, lysimeter measurements were used to validate the evapotranspiration values obtained with the surface renewal method. Seasonal rice evapotranspiration was about 750–800 mm. Average daily ETc for mid-season (from 90 to 130 days after sowing) was 5.1, 4.5 and 6.1 mm day−1 for 2001, 2002 and 2003, respectively. The experimental weekly crop coefficients fluctuated in the range of 0.83–1.20 for 2001, 0.81–1.03 for 2002 and 0.84–1.15 for 2003. The total growing season was about 150–160 days. In average, the crop coefficients for the initial (Kcini), mid-season (Kcmid) and late-season stages (Kcend) were 0.92, 1.06 and 1.03, respectively, the length of these stages being about 55, 45 and 25 days, respectively.

Design and analysis of micro-channel heat-exchanger embedded in Low Temperature Co-fire Ceramic (LTCC)

Increased device density, switching speeds of integrated circuits and decrease in package size is placing new demands for high power thermal-management. The convectional method of forced air cooling with passive heat sink can handle heat fluxes up-to 3-5W/cm2; however current microprocessors are operating at levels of 100W/cm2, This demands the usage of novel thermal-management systems. In this work, water-cooling systems with active heat sink are embedded in the substrate. The research involved fabricating LTCC substrates of various configurations - an open-duct substrate, the second with thermal vias and the third with thermal vias and free-standing metal columns and metal foil. Thermal testing was performed experimentally and these results are compared with CFD results. An overall thermal resistance for the base substrate is demonstrated to be 3.4oC/W-cm2. Addition of thermal vias reduces the effective resistance of the system by 7times and further addition of free standing columns reduced it by 20times.

Global surface net-radiation at 5 km from MODIS Terra

Reliable and fine resolution estimates of surface net-radiation are required for estimating latent and sensible heat fluxes between the land surface and the atmosphere. However, currently, fine resolution estimates of net-radiation are not available and consequently it is challenging to develop multi-year estimates of evapotranspiration at scales that can capture land surface heterogeneity and are relevant for policy and decision-making. We developed and evaluated a global net-radiation product at 5 km and 8-day resolution by combining mutually consistent atmosphere and land data from the Moderate Resolution Imaging Spectroradiometer (MODIS) on board Terra. Comparison with net-radiation measurements from 154 globally distributed sites (414 site-years) from the FLUXNET and Surface Radiation budget network (SURFRAD) showed that the net-radiation product agreed well with measurements across seasons and climate types in the extratropics (Wilmott’s index ranged from 0.74 for boreal to 0.63 for Mediterranean sites). Mean absolute deviation between the MODIS and measured net-radiation ranged from 38.0 ± 1.8 W∙m−2 in boreal to 72.0 ± 4.1 W∙m−2 in the tropical climates. The mean bias was small and constituted only 11%, 0.7%, 8.4%, 4.2%, 13.3%, and 5.4% of the mean absolute error in daytime net-radiation in boreal, Mediterranean, temperate-continental, temperate, semi-arid, and tropical climate, respectively. To assess the accuracy of the broader spatiotemporal patterns, we upscaled error-quantified MODIS net-radiation and compared it with the net-radiation estimates from the coarse spatial (1° × 1°) but high temporal resolution gridded net-radiation product from the Clouds and Earth’s Radiant Energy System (CERES). Our estimates agreed closely with the net-radiation estimates from the CERES. Difference between the two was less than 10 W•m−2 in 94% of the total land area. MODIS net-radiation product will be a valuable resource for the science community studying turbulent fluxes and energy budget at the Earth’s surface.

A STUDY OF HEAT TRANSFER AND FLOW CHARACTERISTICS OF RISING TAYLOR BUBBLES

Practical application of flow boiling to ground- and space-based thermal management systems hinges on the ability to predict the system’s heat removal capabilities under expected operating conditions. Research in this field has shown that the heat transfer coefficient within two-phase heat exchangers can be largely dependent on the experienced flow regime. This finding has inspired an effort to develop mechanistic heat transfer models for each flow pattern which are likely to outperform traditional empirical correlations. As a contribution to the effort, this work aimed to identify the heat transfer mechanisms for the slug flow regime through analysis of individual Taylor bubbles. An experimental apparatus was developed to inject single vapor Taylor bubbles into co-currently flowing liquid HFE 7100. The heat transfer was measured as the bubble rose through a 6 mm inner diameter heated tube using an infrared thermography technique. High-speed flow visualization was obtained and the bubble film thickness measured in an adiabatic section. Experiments were conducted at various liquid mass fluxes (43-200 kg/m2s) and gravity levels (0.01g-1.8g) to characterize the effect of bubble drift velocity on the heat transfer mechanisms. Variable gravity testing was conducted during a NASA parabolic flight campaign. Results from the experiments showed that the drift velocity strongly affects the hydrodynamics and heat transfer of single elongated bubbles. At low gravity levels, bubbles exhibited shapes characteristic of capillary flows and the heat transfer enhancement due to the bubble was dominated by conduction through the thin film. At moderate to high gravity, traditional Taylor bubbles provided small values of enhancement within the film, but large peaks in the wake heat transfer occurred due to turbulent vortices induced by the film plunging into the trailing liquid slug. Characteristics of the wake heat transfer profiles were analyzed and related to the predicted velocity field. Results were compared and shown to agree with numerical simulations of colleagues from EPFL, Switzerland. In addition, a preliminary study was completed on the effect of a Taylor bubble passing through nucleate flow boiling, showing that the thinning thermal boundary layer within the film suppressed nucleation, thereby decreasing the heat transfer coefficient.

Échanges d’énergie et d’eau des écosystèmes nordiques dans un contexte de changement climatique

Le réchauffement climatique affecte fortement les régions nordiques du Canada où le dégel du pergélisol discontinu à sa limite sud est accompagné du mouvement de la limite des arbres vers le nord en zone de pergélisol continu. Ces altérations faites aux paysages de la Taïga des Plaines sont le point de départ de plusieurs rétroactions puisque les changements apportés aux caractéristiques de la surface (au niveau de l’albédo, l’humidité du sol et la rugosité de la surface) vont à leur tour entraîner des modifications biophysiques et éventuellement influencer l’augmentation ou la diminution subséquente des températures et de l’humidité de l’air. Seulement, il y a un nombre important de facteurs d’influence qu’il est difficile de projeter toutes les boucles rétroactives qui surviendront avec les présents changements climatiques en régions nordiques. Dans le but de caractériser les échanges d’eau et d’énergie entre la surface et l’atmosphère de trois sites des Territoires du Nord-Ouest subissant les conséquences de l’augmentation des températures de l’air, la méthode micro-météorologique de covariance des turbulences fut utilisée en 2013 aux sites de Scotty Creek (forêt boréale et tourbière nordique en zone de pergélisol sporadique-discontinu), de Havikpak Creek (forêt boréale nordique en zone de pergélisol continu) et de Trail Valley Creek (toundra arctique en zone de pergélisol continu). En identifiant les procédés biotiques et abiotiques (ex. intensité lumineuse, disponibilité en eau, etc.) d’évapotranspiration aux trois sites, les contrôles par l’eau et l’énergie furent caractérisés et permirent ainsi de projeter une augmentation de la limitation en eau, mais surtout en énergie du site de Trail Valley Creek. La répartition de l’énergie projetée est semblable à celle de Havikpak Creek, avec une augmentation de la proportion du flux de chaleur sensible au détriment de celui latent suite aux modifications des caractéristiques de la surface (albédo, rugosité et humidité du sol). L’augmentation relative du flux d’énergie sensible laisse présager une boucle rétroactive positive de l’augmentation des températures de l’air à ce site. Ensuite, en comparant des données modelées de la hauteur de la couche limite planétaire et des données provenant de profils atmosphériques d’Environnement Canada entre les trois sites, les changements de hauteur de cette couche atmosphérique furent aussi projetés. Trail Valley Creek pourrait connaître une hausse de la hauteur de sa couche limite planétaire avec le temps alors que Scotty Creek connaîtrait une diminution de celle-ci. Ces changements au niveau des couches atmosphériques liés à la répartition des flux d’énergie dans les écosystèmes se répercuteraient alors sur le climat régional de façon difficile à déterminer pour l’instant. Les changements apportés désignent une boucle rétroactive positive des températures de l’air à Trail Valley Creek et l’inverse à Scotty Creek. Les deux axes d’analyse arrivent donc aux mêmes conclusions et soulignent aussi l’importance de l’influence mutuelle entre le climat et les caractéristiques spécifiques des écosystèmes à la surface.