Numerical Investigation on the Heat Transfer Enhancement Using Micro/Nano Phase-Change Particulate Flow

The introduction of phase change material fluid and nanofluid in micro-channel heat sink design can significantly increase the cooling capacity of the heat sink because of the unique features of these two kinds of fluids. To better assist the design of a high performance micro-channel heat sink using phase change fluid and nanofluid, the heat transfer enhancement mechanism behind the flow with such fluids must be completely understood. A detailed parametric study is conducted to further investigate the heat transfer enhancement of the phase change material particle suspension flow, by using the two-phase non-thermal-equilibrium model developed by Hao and Tao (2004). The parametric study is conducted under normal conditions with Reynolds numbers of Re=600-900 and phase change material particle concentrations ¡Ü0.25 , as well as extreme conditions of very low Reynolds numbers (Re < 50) and high phase change material particle concentration (0.5-0.7) slurry flow. By using the two newly-defined parameters, named effectiveness factor and performance index, respectively, it is found that there exists an optimal relation between the channel design parameters, particle volume fraction, Reynolds number, and the wall heat flux. The influence of the particle volume fraction, particle size, and the particle viscosity, to the phase change material suspension flow, are investigated and discussed. The model was validated by available experimental data. The conclusions will assist designers in making their decisions that relate to the design or selection of a micro-pump suitable for micro or mini scale heat transfer devices. To understand the heat transfer enhancement mechanism of the nanofluid flow from the particle level, the lattice Boltzmann method is used because of its mesoscopic feature and its many numerical advantages. By using a two-component lattice Boltzmann model, the heat transfer enhancement of the nanofluid is analyzed, through incorporating the different forces acting on the nanoparticles to the two-component lattice Boltzmann model. It is found that the nanofluid has better heat transfer enhancement at low Reynolds numbers, and the Brownian motion effect of the nanoparticles will be weakened by the increase of flow speed.

Monitoring of active layer dynamics at two transects on Svalbard using multi-channel ground-penetrating radar

Multi-channel ground-penetrating radar is used to investigate the late-summer evolution of the thaw depth and the average soil water content of the thawed active layer at a high-arctic continuous permafrost site on Svalbard, Norway. Between mid of August and mid of September 2008, five surveys have been conducted over transect lengths of 130 and 175 m each. The maximum thaw depths range from 1.6 m to 2.0 m, so that they are among the deepest thaw depths recorded for Svalbard so far. The thaw depths increase by approximately 0.2 m between mid of August and beginning of September and subsequently remain constant until mid of September. The thaw rates are approximately constant over the entire length of the transects within the measurement accuracy of about 5 to 10 cm. The average volumetric soil water content of the thawed soil varies between 0.18 and 0.27 along the investigated transects. While the measurements do not show significant changes in soil water content over the first four weeks of the study, strong precipitation causes an increase in average soil water content of up to 0.04 during the last week. These values are in good agreement with evapotranspiration and precipitation rates measured in the vicinity of the the study site. While we cannot provide conclusive reasons for the detected spatial variability of the thaw depth at the study site, our measurements show that thaw depth and average soil water content are not directly correlated. The study demonstrates the potential of multi-channel ground-penetrating radar for mapping thaw depth in permafrost areas. The novel non-invasive technique is particularly useful when the thaw depth exceeds 1.5 m, so that it is hardly accessible by manual probing. In addition, multi-channel ground-penetrating radar holds potential for mapping the latent heat content of the active layer and for estimating weekly to monthly averages of the ground heat flux during the thaw period.

(Table 1) Age determination of sediment cores CH07-98-22 and KNR140-2-59

Western subtropical North Atlantic oceanic and atmospheric circulations connect tropical and subpolar climates. Variations in these circulations can generate regional climate anomalies that are not reflected in Northern Hemisphere averages. Assessing the significance of anthropogenic climate change at regional scales requires proxy records that allow recent trends to be interpreted in the context of long-term regional variability. We present reconstructions of Gulf Stream sea surface temperature (SST) and hydrographic variability during the past two millennia based on the magnesium/calcium ratio and oxygen isotopic composition of planktic foraminifera preserved in two western subtropical North Atlantic sediment cores. Reconstructed SST suggests low-frequency variability of ~1°C during an interval that includes the Medieval Climate Anomaly (MCA) and the Little Ice Age (LIA). A warm interval near 1250 A.D. is distinct from regional and hemispheric temperature, possibly reflecting regional variations in ocean-atmosphere heat flux associated with changes in atmospheric circulation (e.g., the North Atlantic Oscillation) or the Atlantic Meridional Overturning Circulation. Seawater d18O, which is marked by a fresher MCA and a more saline LIA, covaries with meridional migrations of the Atlantic Intertropical Convergence Zone. The northward advection of tropical salinity anomalies by mean surface currents provides a plausible mechanism linking Carolina Slope and tropical Atlantic hydrology.

Gas hydrate investigation at the Amsterdam mud volcano

We investigated gas hydrate in situ inventories as well as the composition and principal transport mechanisms of fluids expelled at the Amsterdam mud volcano (AMV; 2,025 m water depth) in the Eastern Mediterranean Sea. Pressure coring (the only technique preventing hydrates from decomposition during recovery) was used for the quantification of light hydrocarbons in near-surface deposits. The cores (up to 2.5 m in length) were retrieved with an autoclave piston corer, and served for analyses of gas quantities and compositions, and pore-water chemistry. For comparison, gravity cores from sites at the summit and beyond the AMV were analyzed. A prevalence of thermogenic light hydrocarbons was inferred from average C1/C2+ ratios <35 and d13C-CH4 values of -50.6 per mil. Gas venting from the seafloor indicated methane oversaturation, and volumetric gas-sediment ratios of up to 17.0 in pressure cores taken from the center demonstrated hydrate presence at the time of sampling. Relative enrichments in ethane, propane, and iso-butane in gas released from pressure cores, and from an intact hydrate piece compared to venting gas suggest incipient crystallization of hydrate structure II (sII). Nonetheless, the co-existence of sI hydrate can not be excluded from our dataset. Hydrates fill up to 16.7% of pore volume within the sediment interval between the base of the sulfate zone and the maximum sampling depth at the summit. The concave-down shapes of pore-water concentration profiles recorded in the center indicate the influence of upward-directed advection of low-salinity fluids/fluidized mud. Furthermore, the SO42- and Ba2+ pore-water profiles in the central part of the AMV demonstrate that sulfate reduction driven by the anaerobic oxidation of methane is complete at depths between 30 cm and 70 cm below seafloor. Our results indicate that methane oversaturation, high hydrostatic pressure, and elevated pore-water activity caused by low salinity promote fixing of considerable proportions of light hydrocarbons in shallow hydrates even at the summit of the AMV, and possibly also of other MVs in the region. Depending on their crystallographic structure, however, hydrates will already decompose and release hydrocarbon masses if sediment temperatures exceed ca. 19.3°C and 21.0°C, respectively. Based on observations from other mud volcanoes, the common occurrence of such temperatures induced by heat flux from below into the immediate subsurface appears likely for the AMV.

Geochemistry of basalts from the Caribbean Large Igneous Province

Oceanic flood basalts are poorly understood, short-term expressions of highly increased heat flux and mass flow within the convecting mantle. The uniqueness of the Caribbean Large Igneous Province (CLIP, 92-74 Ma) with respect to other Cretaceous oceanic plateaus is its extensive sub-aerial exposures, providing an excellent basis to investigate the temporal and compositional relationships within a starting plume head. We present major element, trace element and initial Sr-Nd-Pb isotope composition of 40 extrusive rocks from the Caribbean Plateau, including onland sections in Costa Rica, Colombia and Curaçao as well as DSDP Sites in the Central Caribbean. Even though the lavas were erupted over an area of ~3*10**6 km**2, the majority have strikingly uniform incompatible element patterns (La/Yb=0.96+/-0.16, n=64 out of 79 samples, 2sigma) and initial Nd-Pb isotopic compositions (e.g. 143Nd/144Ndin=0.51291+/-3, epsilon-Nd i=7.3+/-0.6, 206Pb/204Pbin=18.86+/-0.12, n=54 out of 66, 2sigma). Lavas with endmember compositions have only been sampled at the DSDP Sites, Gorgona Island (Colombia) and the 65-60 Ma accreted Quepos and Osa igneous complexes (Costa Rica) of the subsequent hotspot track. Despite the relatively uniform composition of most lavas, linear correlations exist between isotope ratios and between isotope and highly incompatible trace element ratios. The Sr-Nd-Pb isotope and trace element signatures of the chemically enriched lavas are compatible with derivation from recycled oceanic crust, while the depleted lavas are derived from a highly residual source. This source could represent either oceanic lithospheric mantle left after ocean crust formation or gabbros with interlayered ultramafic cumulates of the lower oceanic crust. High 3He/4He in olivines of enriched picrites at Quepos are ~12 times higher than the atmospheric ratio suggesting that the enriched component may have once resided in the lower mantle. Evaluation of the Sm-Nd and U-Pb isotope systematics on isochron diagrams suggests that the age of separation of enriched and depleted components from the depleted MORB source mantle could have been <=500 Ma before CLIP formation and interpreted to reflect the recycling time of the CLIP source. Mantle plume heads may provide a mechanism for transporting large volumes of possibly young recycled oceanic lithosphere residing in the lower mantle back into the shallow MORB source mantle.

Combined Conduction-Convection-Radiation Heat Transfer of Slip Flow inside a Micro-channel Filled with a Cellular Porous Material

Combined conduction–convection–radiation heat transfer is investigated numerically in a micro-channel filled with a saturated cellular porous medium, with the channel walls held at a constant heat flux. Invoking the velocity slip and temperature jump, the thermal behaviour of the porous–fluid system are studied by considering hydrodynamically fully developed flow and applying the Darcy–Brinkman flow model. One energy equation model based on the local thermal equilibrium condition is adopted to evaluate the temperature field within the porous medium. Combined conduction and radiation heat transfer is treated as an effective conduction process with a temperature-dependent effective thermal conductivity. Results are reported in terms of the average Nusselt number and dimensionless temperature distribution, as a function of velocity slip coefficient, temperature jump coefficient, porous medium shape parameter and radiation parameters. Results show that increasing the radiation parameter (Tr)(Tr) and the temperature jump coefficient flattens the dimensionless temperature profile. The Nusselt numbers are more sensitive to the variation in the temperature jump coefficient rather than to the velocity slip coefficient. Such that for high porous medium shape parameter, the Nusselt number is found to be independent of velocity slip. Furthermore, it is found that as the temperature jump coefficient increases, the Nusselt number decrease. In addition, for high temperature jump coefficients, the Nusselt number is found to be insensitive to the radiation parameters and porous medium shape parameter. It is also concluded that compared with the conventional macro-channels, wherein using a porous material enhances the rate of heat transfer (up to about 40 % compared to the clear channel), insertion of a porous material inside a micro-channel in slip regime does not effectively enhance the rate of heat transfer that is about 2 %.

Effect of boiling on deposition of metallic colloids

In this study the relationship between heterogeneous nucleate boiling surfaces and deposition of suspended metallic colloidal particles, popularly known as crud or corrosion products in process industries, on those heterogeneous sites is investigated. Various researchers have reported that hematite is a major constituent of crud which makes it the primary material of interest; however the models developed in this work are irrespective of material choice. Qualitative hypotheses on the deposition process under boiling as proposed by previous researchers have been tested, which fail to provide explanations for several physical mechanisms observed and analyzed. In this study a quantitative model of deposition rate has been developed on the basis of bubble dynamics and colloid-surface interaction potential. Boiling from a heating surface aids in aggregation of the metallic particulates viz. nano-particles, crud particulate, etc. suspended in a liquid, which helps in transporting them to heating surfaces. Consequently, clusters of particles deposit onto the heating surfaces due to various interactive forces, resulting in formation of porous or impervious layers. The deposit layer grows or recedes depending upon variations in interparticle and surface forces, fluid shear, fluid chemistry, etc. This deposit layer in turn affects the rate of bubble generation, formation of porous chimneys, critical heat flux (CHF) of surfaces, activation and deactivation of nucleation sites on the heating surfaces. Several problems are posed due to the effect of boiling on colloidal deposition, which range from research initiatives involving nano-fluids as a heat transfer medium to industrial applications such as light water nuclear reactors. In this study, it is attempted to integrate colloid and surface science with vapor bubble dynamics, boiling heat transfer and evaporation rate. Pool boiling experiments with dilute metallic colloids have been conducted to investigate several parameters impacting the system. The experimental data available in the literature is obtained by flow experiments, which do not help in correlating boiling mechanism with the deposition amount or structure. With the help of experimental evidences and analysis, previously proposed hypothesis for particle transport to the contact line due to hydrophobicity has been challenged. The experimental observations suggest that deposition occurs around the bubble surface contact line and extends underneath area of the bubble microlayer as well. During the evaporation the concentration gradient of a non-volatile species is created, which induces osmotic pressure. The osmotic pressure developed inside the microlayer draws more particles inside the microlayer region or towards contact line. The colloidal escape time is slower than the evaporation time, which leads to the aggregation of particles in the evaporating micro-layer. These aggregated particles deposit onto or are removed from the heating surface, depending upon their total interaction potential. Interaction potential has been computed with the help of surface charge and van der Waals potential for the materials in aqueous solutions. Based upon the interaction-force boundary layer thickness, which is governed by debye radius (or ionic concentration and pH), a simplified quantitative model for the attachment kinetics is proposed. This attachment kinetics model gives reasonable results in predicting attachment rate against data reported by previous researchers. The attachment kinetics study has been done for different pH levels and particle sizes for hematite particles. Quantification of colloidal transport under boiling scenarios is done with the help of overall average evaporation rates because generally waiting times for bubbles at the same position is much larger than growth times. In other words, from a larger measurable scale perspective, frequency of bubbles dictates the rate of collection of particles rather than evaporation rate during micro-layer evaporation of one bubble. The combination of attachment kinetics and colloidal transport kinetics has been used to make a consolidated model for prediction of the amount of deposition and is validated with the help of high fidelity experimental data. In an attempt to understand and explain boiling characteristics, high speed visualization of bubble dynamics from a single artificial large cavity and multiple naturally occurring cavities is conducted. A bubble growth and departure dynamics model is developed for artificial active sites and is validated with the experimental data. The variation of bubble departure diameter with wall temperature is analyzed with experimental results and shows coherence with earlier studies. However, deposit traces after boiling experiments show that bubble contact diameter is essential to predict bubble departure dynamics, which has been ignored previously by various researchers. The relationship between porosity of colloid deposits and bubbles under the influence of Jakob number, sub-cooling and particle size has been developed. This also can be further utilized in variational wettability of the surface. Designing porous surfaces can having vast range of applications varying from high wettability, such as high critical heat flux boilers, to low wettability, such as efficient condensers.

A Generalized Methodology to Characterize Composite Materials for Pyrolysis Models

The predictive capabilities of computational fire models have improved in recent years such that models have become an integral part of many research efforts. Models improve the understanding of the fire risk of materials and may decrease the number of expensive experiments required to assess the fire hazard of a specific material or designed space. A critical component of a predictive fire model is the pyrolysis sub-model that provides a mathematical representation of the rate of gaseous fuel production from condensed phase fuels given a heat flux incident to the material surface. The modern, comprehensive pyrolysis sub-models that are common today require the definition of many model parameters to accurately represent the physical description of materials that are ubiquitous in the built environment. Coupled with the increase in the number of parameters required to accurately represent the pyrolysis of materials is the increasing prevalence in the built environment of engineered composite materials that have never been measured or modeled. The motivation behind this project is to develop a systematic, generalized methodology to determine the requisite parameters to generate pyrolysis models with predictive capabilities for layered composite materials that are common in industrial and commercial applications. This methodology has been applied to four common composites in this work that exhibit a range of material structures and component materials. The methodology utilizes a multi-scale experimental approach in which each test is designed to isolate and determine a specific subset of the parameters required to define a material in the model. Data collected in simultaneous thermogravimetry and differential scanning calorimetry experiments were analyzed to determine the reaction kinetics, thermodynamic properties, and energetics of decomposition for each component of the composite. Data collected in microscale combustion calorimetry experiments were analyzed to determine the heats of complete combustion of the volatiles produced in each reaction. Inverse analyses were conducted on sample temperature data collected in bench-scale tests to determine the thermal transport parameters of each component through degradation. Simulations of quasi-one-dimensional bench-scale gasification tests generated from the resultant models using the ThermaKin modeling environment were compared to experimental data to independently validate the models.

Superconducting Transition-Edge Sensor Physics

Despite record-setting performance demonstrated by superconducting Transition Edge Sensors (TESs) and growing utilization of the technology, a theoretical model of the physics governing TES devices superconducting phase transition has proven elusive. Earlier attempts to describe TESs assumed them to be uniform superconductors. Sadleir et al. 2010 shows that TESs are weak links and that the superconducting order parameter strength has significant spatial variation. Measurements are presented of the temperature T and magnetic field B dependence of the critical current Ic measured over 7 orders of magnitude on square Mo/Au bilayers ranging in length from 8 to 290 microns. We find our measurements have a natural explanation in terms of a spatially varying order parameter that is enhanced in proximity to the higher transition temperature superconducting leads (the longitudinal proximity effect) and suppressed in proximity to the added normal metal structures (the lateral inverse proximity effect). These in-plane proximity effects and scaling relations are observed over unprecedentedly long lengths (in excess of 1000 times the mean free path) and explained in terms of a Ginzburg-Landau model. Our low temperature Ic(B) measurements are found to agree with a general derivation of a superconducting strip with an edge or geometric barrier to vortex entry and we also derive two conditions that lead to Ic rectification. At high temperatures the Ic(B) exhibits distinct Josephson effect behavior over long length scales and following functional dependences not previously reported. We also investigate how film stress changes the transition, explain some transition features in terms of a nonequilibrium superconductivity effect, and show that our measurements of the resistive transition are not consistent with a percolating resistor network model.

AIR SIDE HEAT TRANSFER ENHANCEMENT IN HEAT EXCHANGERS UTILIZING INNOVATIVE DESIGNS AND THE ADDITIVE MANUFACTURING TECHNIQUE

Over the last decade, rapid development of additive manufacturing techniques has allowed the fabrication of innovative and complex designs. One field that can benefit from such technology is heat exchanger fabrication, as heat exchanger design has become more and more complex due to the demand for higher performance particularly on the air side of the heat exchanger. By employing the additive manufacturing, a heat exchanger design was successfully realized, which otherwise would have been very difficult to fabricate using conventional fabrication technologies. In this dissertation, additive manufacturing technique was implemented to fabricate an advanced design which focused on a combination of heat transfer surface and fluid distribution system. Although the application selected in this dissertation is focused on power plant dry cooling applications, the results of this study can directly and indirectly benefit other sectors as well, as the air-side is often the limiting side for in liquid or single phase cooling applications. Two heat exchanger designs were studied. One was an advanced metallic heat exchanger based on manifold-microchannel technology and the other was a polymer heat exchanger based on utilization of prime surface technology. Polymer heat exchangers offer several advantages over metals such as antifouling, anticorrosion, lightweight and often less expensive than comparable metallic heat exchangers. A numerical modeling and optimization were performed to calculate a design that yield an optimum performance. The optimization results show that significant performance enhancement is noted compared to the conventional heat exchangers like wavy fins and plain plate fins. Thereafter, both heat exchangers were scaled down and fabricated using additive manufacturing and experimentally tested. The manifold-micro channel design demonstrated that despite some fabrication inaccuracies, compared to a conventional wavy-fin surface, 15% - 50% increase in heat transfer coefficient was possible for the same pressure drop value. In addition, if the fabrication inaccuracy can be eliminated, an even larger performance enhancement is predicted. Since metal based additive manufacturing is still in the developmental stage, it is anticipated that with further refinement of the manufacturing process in future designs, the fabrication accuracy can be improved. For the polymer heat exchanger, by fabricating a very thin wall heat exchanger (150μm), the wall thermal resistance, which usually becomes the limiting side for polymer heat exchanger, was calculated to account for only up to 3% of the total thermal resistance. A comparison of air-side heat transfer coefficient of the polymer heat exchanger with some of the commercially available plain plate fin surface heat exchangers show that polymer heat exchanger performance is equal or superior to plain plate fin surfaces. This shows the promising potential for polymer heat exchangers to compete with conventional metallic heat exchangers when an additive manufacturing-enabled fabrication is utilized. Major contributions of this study are as follows: (1) For the first time demonstrated the potential of additive manufacturing in metal printing of heat exchangers that benefit from a sophisticated design to yield a performance substantially above the respective conventional systems. Such heat exchangers cannot be fabricated with the conventional fabrication techniques. (2) For the first time demonstrated the potential of additive manufacturing to produce polymer heat exchangers that by design minimize the role of thermal conductivity and deliver a thermal performance equal or better that their respective metallic heat exchangers. In addition of other advantages of polymer over metal like antifouling, anticorrosion, and lightweight. Details of the work are documented in respective chapters of this thesis.

Heat Transfer Measurements in a Supersonic Film Flow

This thesis presents measurements of wall heat flux and flow structure in a canonical film cooling configuration with Mach 2.3 core flow in which the coolant is injected parallel to the wall through a two-dimensional louver. Four operating conditions are investigated: no film (i.e. flow over a rearward-facing step), subsonic film, pressure-matched film, and supersonic film. The overall objective is to provide a set of experimental data with well characterized boundary conditions that can be used for code validation. The results are compared to RANS and LES simulations which overpredict heat transfer in the subsonic film cases and underpredict heat transfer in supersonic cases after film breakdown. The thesis also describes a number of improvements that were made to the experimental facility including new Schlieren optics, a better film heater, more data at more locations, and a verification of the heat flux measurement hardware and data reduction methods.

É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.

É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.

Lake–atmosphere interactions at Alqueva reservoir: a case study in the summer of 2014

The study of lake–atmosphere interactions was the main purpose of a 2014 summer experiment at Alqueva reservoir in Portugal. Near-surface fluxes of momentum, heat and mass [water vapour (H2O) and carbon dioxide (CO2)] were obtained with the new Campbell Scientific’s IRGASON Integrated Open-Path CO2/H2O Gas Analyser and 3D Sonic Anemometer between 2 June and 2 October. On average, the reservoir was releasing energy in the form of sensible and latent heat flux during the study period. At the end of the 75 d, the total evaporation was estimated as 490.26 mm. A high correlation was found between the latent heat flux and the wind speed (R = 0.97). The temperature gradient between air and water was positive between 12 and 21 UTC, causing a negative sensible heat flux, and negative during the rest of the day, triggering a positive sensible heat flux. The reservoir acted as a sink of atmospheric CO2 with an average rate of −0.026 mg m−2 s−1. However, at a daily scale we found an unexpected uptake between 0 and 9 UTC and almost null flux between 13 and 19 UTC. Potential reasons for this result are further discussed. The net radiation was recorded for the same period and water column heat storage was estimated using water temperature profiles. The energy balance closure for the analysed period was 81%. In-water solar spectral downwelling irradiance profiles were measured with a new device allowing measurements independent of the solar zenith angle, which enabled the computation of the attenuation coefficient of light in the water column. The average attenuation coefficient for the photosynthetically active radiation spectral region varied from 0.849 ± 0.025 m−1 on 30 July to 1.459 ± 0.007 m−1 on 25 September.