Daily pan-Arctic sea-ice lead maps for 2003-2015, with links to maps in netCDF format

The presence of sea-ice leads represents a key feature of the Arctic sea ice cover. Leads promote the flux of sensible and latent heat from the ocean to the cold winter atmosphere and are thereby crucial for air-sea-ice-ocean interactions. We here apply a binary segmentation procedure to identify leads from MODIS thermal infrared imagery on a daily time scale. The method separates identified leads into two uncertainty categories, with the high uncertainty being attributed to artifacts that arise from warm signatures of unrecognized clouds. Based on the obtained lead detections, we compute quasi-daily pan-Arctic lead maps for the months of January to April, 2003-2015. Our results highlight the marginal ice zone in the Fram Strait and Barents Sea as the primary region for lead activity. The spatial distribution of the average pan-Arctic lead frequencies reveals, moreover, distinct patterns of predominant fracture zones in the Beaufort Sea and along the shelf-breaks, mainly in the Siberian sector of the Arctic Ocean as well as the well-known polynya and fast-ice locations. Additionally, a substantial inter-annual variability of lead occurrences in the Arctic is indicated.

An assessment of the latent and sensible heat flux on the simulated regional climate over Southwestern South Atlantic Ocean

A Regional Climate Model (RegCM3) 10-year (1990-1999) simulation over southwestern South Atlantic Ocean (SAO) is evaluated to assess the mean climatology and the simulation errors of turbulent fluxes over the sea. Moreover, the relationship between these fluxes and the rainfall over some cyclogenetic areas is also analyzed. The RegCM3 results are validated using some reanalyses datasets (ERA40, R2, GPCP and WHOI). The summer and winter spatial patterns of latent and sensible heat fluxes simulated by the RegCM3 are in agreement with the reanalyses (WHOI, R2 and ERA40). They show large latent heat fluxes exchange in the subtropical SAO and at higher latitudes in the warm waters of Brazil Current. In particular, the magnitude of RegCM3 latent heat fluxes is similar to the WHOI, which is probably related to two factors: (a) small specific humidity bias, and (b) the RegCM3 flux algorithm. In contrast, the RegCM3 presents large overestimation of sensible heat flux, though it simulates well their spatial pattern. This simulation error is associated with the RegCM3 underestimation of the 2-m air temperature. In southwestern SAO, in three known cyclogenetic areas, the reanalyses and the RegCM3 show the existence of different physical mechanisms that control the annual cycles of latent/sensible heating and rainfall. It is shown that over the eastern coast of Uruguay (35A degrees-43A degrees S) and the southeastern coast of Argentina (44A degrees-52A degrees S) the sea-air moisture and heat exchange play an important role to control the annual cycle of precipitation. This does not happen on the south/southeastern coast of Brazil.

Multi-scale sensible heat fluxes in the urban environment from large aperture scintillometry and eddy covariance

Sensible heat fluxes (QH) are determined using scintillometry and eddy covariance over a suburban area. Two large aperture scintillometers provide spatially integrated fluxes across path lengths of 2.8 km and 5.5 km over Swindon, UK. The shorter scintillometer path spans newly built residential areas and has an approximate source area of 2-4 km2, whilst the long path extends from the rural outskirts to the town centre and has a source area of around 5-10 km2. These large-scale heat fluxes are compared with local-scale eddy covariance measurements. Clear seasonal trends are revealed by the long duration of this dataset and variability in monthly QH is related to the meteorological conditions. At shorter time scales the response of QH to solar radiation often gives rise to close agreement between the measurements, but during times of rapidly changing cloud cover spatial differences in the net radiation (Q*) coincide with greater differences between heat fluxes. For clear days QH lags Q*, thus the ratio of QH to Q* increases throughout the day. In summer the observed energy partitioning is related to the vegetation fraction through use of a footprint model. The results demonstrate the value of scintillometry for integrating surface heterogeneity and offer improved understanding of the influence of anthropogenic materials on surface-atmosphere interactions.

Heat fluxes and surface radiation measurements from Samoylov Island, Lena Delta, Siberia, October 2007 to March 2008

In this study, we present the winter time surface energy balance at a polygonal tundra site in northern Siberia based on independent measurements of the net radiation, the sensible heat flux and the ground heat flux from two winter seasons. The latent heat flux is inferred from measurements of the atmospheric turbulence characteristics and a model approach. The long-wave radiation is found to be the dominant factor in the surface energy balance. The radiative losses are balanced to about 60 % by the ground heat flux and almost 40 % by the sensible heat fluxes, whereas the contribution of the latent heat flux is small. The main controlling factors of the surface energy budget are the snow cover, the cloudiness and the soil temperature gradient. Large spatial differences in the surface energy balance are observed between tundra soils and a small pond. The ground heat flux released at a freezing pond is by a factor of two higher compared to the freezing soil, whereas large differences in net radiation between the pond and soil are only observed at the end of the winter period. Differences in the surface energy balance between the two winter seasons are found to be related to differences in snow depth and cloud cover which strongly affect the temperature evolution and the freeze-up at the investigated pond.

Determination of turbulent heat fluxes using a large aperture scintillometer over undulating mixed agricultural terrain

Scintillometry is an established technique for determining large areal average sensible heat fluxes. The scintillometer measurement is related to sensible heat flux via Monin–Obukhov similarity theory, which was developed for ideal homogeneous land surfaces. In this study it is shown that judicious application of scintillometry over heterogeneous mixed agriculture on undulating topography yields valid results when compared to eddy covariance (EC). A large aperture scintillometer (LAS) over a 2.4 km path was compared with four EC stations measuring sensible (H) and latent (LvE) heat fluxes over different vegetation (cereals and grass) which when aggregated were representative of the LAS source area. The partitioning of available energy into H and LvE varied strongly for different vegetation types, with H varying by a factor of three between senesced winter wheat and grass pasture. The LAS derived H agrees (one-to-one within the experimental uncertainty) with H aggregated from EC with a high coefficient of determination of 0.94. Chronological analysis shows individual fields may have a varying contribution to the areal average sensible heat flux on short (weekly) time scales due to phenological development and changing soil moisture conditions. Using spatially aggregated measurements of net radiation and soil heat flux with H from the LAS, the areal averaged latent heat flux (LvELAS) was calculated as the residual of the surface energy balance. The regression of LvELAS against aggregated LvE from the EC stations has a slope of 0.94, close to ideal, and demonstrates that this is an accurate method for the landscape-scale estimation of evaporation over heterogeneous complex topography.

Multi-season eddy covariance observations of energy, water and carbon fluxes over a suburban area in Swindon, UK

Eddy covariance measurements of the turbulent sensible heat, latent heat and carbon dioxide fluxes for 12 months (2011–2012) are reported for the first time for a suburban area in the UK. The results from Swindon are comparable to suburban studies of similar surface cover elsewhere but reveal large seasonal variability. Energy partitioning favours turbulent sensible heat during summer (midday Bowen ratio 1.4–1.6) and latent heat in winter (0.05–0.7). A significant proportion of energy is stored (and released) by the urban fabric and the estimated anthropogenic heat flux is small but non-negligible (0.5–0.9 MJ m−2 day−1). The sensible heat flux is negative at night and for much of winter daytimes, reflecting the suburban nature of the site (44% vegetation) and relatively low built fraction (16%). Latent heat fluxes appear to be water limited during a dry spring in both 2011 and 2012, when the response of the surface to moisture availability can be seen on a daily timescale. Energy and other factors are more relevant controls at other times; at night the wind speed is important. On average, surface conductance follows a smooth, asymmetrical diurnal course peaking at around 6–9 mm s−1, but values are larger and highly variable in wet conditions. The combination of natural (vegetative) and anthropogenic (emission) processes is most evident in the temporal variation of the carbon flux: significant photosynthetic uptake is seen during summer, whilst traffic and building emissions explain peak release in winter (9.5 g C m−2 day−1). The area is a net source of CO2 annually. Analysis by wind direction highlights the role of urban vegetation in promoting evapotranspiration and offsetting CO2 emissions, especially when contrasted against peak traffic emissions from sectors with more roads. Given the extent of suburban land use, these results have important implications for understanding urban energy, water and carbon dynamics.

An assessment of air-sea heat fluxes from ocean and coupled reanalyses

Sixteen monthly air–sea heat flux products from global ocean/coupled reanalyses are compared over 1993–2009 as part of the Ocean Reanalysis Intercomparison Project (ORA-IP). Objectives include assessing the global heat closure, the consistency of temporal variability, comparison with other flux products, and documenting errors against in situ flux measurements at a number of OceanSITES moorings. The ensemble of 16 ORA-IP flux estimates has a global positive bias over 1993–2009 of 4.2 ± 1.1 W m−2. Residual heat gain (i.e., surface flux + assimilation increments) is reduced to a small positive imbalance (typically, +1–2 W m−2). This compensation between surface fluxes and assimilation increments is concentrated in the upper 100 m. Implied steady meridional heat transports also improve by including assimilation sources, except near the equator. The ensemble spread in surface heat fluxes is dominated by turbulent fluxes (>40 W m−2 over the western boundary currents). The mean seasonal cycle is highly consistent, with variability between products mostly <10 W m−2. The interannual variability has consistent signal-to-noise ratio (~2) throughout the equatorial Pacific, reflecting ENSO variability. Comparisons at tropical buoy sites (10°S–15°N) over 2007–2009 showed too little ocean heat gain (i.e., flux into the ocean) in ORA-IP (up to 1/3 smaller than buoy measurements) primarily due to latent heat flux errors in ORA-IP. Comparisons with the Stratus buoy (20°S, 85°W) over a longer period, 2001–2009, also show the ORA-IP ensemble has 16 W m−2 smaller net heat gain, nearly all of which is due to too much latent cooling caused by differences in surface winds imposed in ORA-IP.

Hurricane disturbance and recovery of energy balance, CO2 fluxes and canopy structure in a mangrove forest of the Florida Everglades

Eddy covariance (EC) estimates of carbon dioxide (CO2) fluxes and energy balance are examined to investigate the functional responses of a mature mangrove forest to a disturbance generated by Hurricane Wilma on October 24, 2005 in the Florida Everglades. At the EC site, high winds from the hurricane caused nearly 100% defoliation in the upper canopy and widespread tree mortality. Soil temperatures down to -50 cm increased, and air temperature lapse rates within the forest canopy switched from statically stable to statically unstable conditions following the disturbance. Unstable conditions allowed more efficient transport of water vapor and CO2 from the surface up to the upper canopy layer. Significant increases in latent heat fluxes (LE) and nighttime net ecosystem exchange (NEE) were also observed and sensible heat fluxes (H) as a proportion of net radiation decreased significantly in response to the disturbance. Many of these impacts persisted through much of the study period through 2009. However, local albedo and MODIS (Moderate Resolution Imaging Spectro-radiometer) data (the Enhanced Vegetation Index) indicated a substantial proportion of active leaf area recovered before the EC measurements began 1 year after the storm. Observed changes in the vertical distribution and the degree of clumping in newly emerged leaves may have affected the energy balance. Direct comparisons of daytime NEE values from before the storm and after our measurements resumed did not show substantial or consistent differences that could be attributed to the disturbance. Regression analyses on seasonal time scales were required to differentiate the storm's impact on monthly average daytime NEE from the changes caused by interannual variability in other environmental drivers. The effects of the storm were apparent on annual time scales, and CO2 uptake remained approximately 250 g C m-2 yr-1 lower in 2009 compared to the average annual values measured in 2004-2005. Dry season CO2 uptake was relatively more affected by the disturbance than wet season values. Complex leaf regeneration dynamics on damaged trees during ecosystem recovery are hypothesized to lead to the variable dry versus wet season impacts on daytime NEE. In contrast, nighttime CO2 release (i.e., nighttime respiration) was consistently and significantly greater, possibly as a result of the enhanced decomposition of litter and coarse woody debris generated by the storm, and this effect was most apparent in the wet seasons compared to the dry seasons. The largest pre- and post-storm differences in NEE coincided roughly with the delayed peak in cumulative mortality of stems in 2007-2008. Across the hurricane-impacted region, cumulative tree mortality rates were also closely correlated with declines in peat surface elevation. Mangrove forest-atmosphere interactions are interpreted with respect to the damage and recovery of stand dynamics and soil accretion processes following the hurricane.

The role of anomalous SST and surface fluxes over the Southeastern North Atlantic in the explosive development of windstorm Xynthia

In late February 2010 the extraordinary windstorm Xynthia crossed over Southwestern and Central Europe and caused severe damage, affecting particularly the Spanish and French Atlantic coasts. The storm was embedded in uncommon large-scale atmospheric and boundary conditions prior to and during its development, namely enhanced sea surface temperatures (SST) within the low-level entrainment zone of air masses, an unusual southerly position of the polar jet stream, and a remarkable split jet structure in the upper troposphere. To analyse the processes that led to the rapid intensification of this exceptional storm originating close to the subtropics (30°N), the sensitivity of the cyclone intensification to latent heat release is determined using the regional climate model COSMO-CLM forced with ERA-Interim data. A control simulation with observed SST shows that moist and warm air masses originating from the subtropical North Atlantic were involved in the cyclogenesis process and led to the formation of a vertical tower with high values of potential vorticity (PV). Sensitivity studies with reduced SST or increased laminar boundary roughness for heat led to reduced surface latent heat fluxes. This induced both a weaker and partly retarded development of the cyclone and a weakening of the PV-tower together with reduced diabatic heating rates, particularly at lower and mid levels. We infer that diabatic processes played a crucial role during the phase of rapid deepening of Xynthia and thus to its intensity over the Southeastern North Atlantic. We suggest that windstorms like Xynthia may occur more frequently under future climate conditions due to the warming SSTs and potentially enhanced latent heat release, thus increasing the windstorm risk for Southwestern Europe.

An event of stratospheric air intrusion and its associated secondary surface cyclogenesis over the South Atlantic Ocean

This work presents an analysis of a lowermost stratospheric air intrusion event over the coast of Brazil, which may have been responsible for a secondary surface cyclogenesis over the southwestern Atlantic Ocean. The surface cyclone initiated at 0600 UTC 17 April 1999 in a cold air mass in the rear of a cold front after a primary cyclone developed over the same region. The analysis of the secondary cyclone revealed the presence of lowermost stratospheric air intrusion characterized by anomalous potential vorticity (PV), dry air, and high concentration of ozone in atmospheric column. The system developed on the eastern side of an upper level core of PV anomaly, which induced a cyclonic wind circulation at lower levels and favored the onset of the secondary cyclone. In midlevels (500 hPa), the cutoff low development contributed to reduce the propagation speed of the wave pattern. This feature seemed to (1) allow the low-level cold/dry air to heat/moisten associated with sensible and latent fluxes transferred from the ocean to the atmosphere, which intensified a baroclinic zone parallel to the coast, and (2) contribute to the long duration of the system. The present analysis indicates that this secondary cyclone development could be the result of the coupling between the PV anomaly in the upper levels and low-level air-sea interaction.

Latent heat loss and sweat gland histology of male goats in an equatorial semi-arid environment

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

Computational simulations of latent heat thermal energy storage systems - with innovative and first-principles based simulation for the underlying unsteady melting (and solidification) processes

This thesis develops an effective modeling and simulation procedure for a specific thermal energy storage system commonly used and recommended for various applications (such as an auxiliary energy storage system for solar heating based Rankine cycle power plant). This thermal energy storage system transfers heat from a hot fluid (termed as heat transfer fluid - HTF) flowing in a tube to the surrounding phase change material (PCM). Through unsteady melting or freezing process, the PCM absorbs or releases thermal energy in the form of latent heat. Both scientific and engineering information is obtained by the proposed first-principle based modeling and simulation procedure. On the scientific side, the approach accurately tracks the moving melt-front (modeled as a sharp liquid-solid interface) and provides all necessary information about the time-varying heat-flow rates, temperature profiles, stored thermal energy, etc. On the engineering side, the proposed approach is unique in its ability to accurately solve – both individually and collectively – all the conjugate unsteady heat transfer problems for each of the components of the thermal storage system. This yields critical system level information on the various time-varying effectiveness and efficiency parameters for the thermal storage system.

Experimental investigation of regular fluids and nanofluids during flow boiling in a single microchannel at different heat fluxes and mass fluxes

The dissipation of high heat flux from integrated circuit chips and the maintenance of acceptable junction temperatures in high powered electronics require advanced cooling technologies. One such technology is two-phase cooling in microchannels under confined flow boiling conditions. In macroscale flow boiling bubbles will nucleate on the channel walls, grow, and depart from the surface. In microscale flow boiling bubbles can fill the channel diameter before the liquid drag force has a chance to sweep them off the channel wall. As a confined bubble elongates in a microchannel, it traps thin liquid films between the heated wall and the vapor core that are subject to large temperature gradients. The thin films evaporate rapidly, sometimes faster than the incoming mass flux can replenish bulk fluid in the microchannel. When the local vapor pressure spike exceeds the inlet pressure, it forces the upstream interface to travel back into the inlet plenum and create flow boiling instabilities. Flow boiling instabilities reduce the temperature at which critical heat flux occurs and create channel dryout. Dryout causes high surface temperatures that can destroy the electronic circuits that use two-phase micro heat exchangers for cooling. Flow boiling instability is characterized by periodic oscillation of flow regimes which induce oscillations in fluid temperature, wall temperatures, pressure drop, and mass flux. When nanofluids are used in flow boiling, the nanoparticles become deposited on the heated surface and change its thermal conductivity, roughness, capillarity, wettability, and nucleation site density. It also affects heat transfer by changing bubble departure diameter, bubble departure frequency, and the evaporation of the micro and macrolayer beneath the growing bubbles. Flow boiling was investigated in this study using degassed, deionized water, and 0.001 vol% aluminum oxide nanofluids in a single rectangular brass microchannel with a hydraulic diameter of 229 µm for one inlet fluid temperature of 63°C and two constant flow rates of 0.41 ml/min and 0.82 ml/min. The power input was adjusted for two average surface temperatures of 103°C and 119°C at each flow rate. High speed images were taken periodically for water and nanofluid flow boiling after durations of 25, 75, and 125 minutes from the start of flow. The change in regime timing revealed the effect of nanoparticle suspension and deposition on the Onset of Nucelate Boiling (ONB) and the Onset of Bubble Elongation (OBE). Cycle duration and bubble frequencies are reported for different nanofluid flow boiling durations. The addition of nanoparticles was found to stabilize bubble nucleation and growth and limit the recession rate of the upstream and downstream interfaces, mitigating the spreading of dry spots and elongating the thin film regions to increase thin film evaporation.

(Table 1) Sea ice and snow characteristics and heat fluxes observed during R/V Aurora Australis cruise to East Antarctica in September/October 2007

The properties of snow on East Antarctic sea ice off Wilkes Land were examined during the Sea Ice Physics and Ecosystem Experiment (SIPEX) in late winter of 2007, focusing on the interaction with sea ice. This observation includes 11 transect lines for the measurement of ice thickness, freeboard, and snow depth, 50 snow pits on 13 ice floes, and diurnal variation of surface heat flux on three ice floes. The detailed profiling of topography along the transects and the d18O, salinity, and density datasets of snow made it possible to examine the snow-sea-ice interaction quantitatively for the first time in this area. In general, the snow displayed significant heterogeneity in types, thickness (mean: 0.14 +- 0.13 m), and density (325 +- 38 kg/m**3), as reported in other East Antarctic regions. High salinity was confined to the lowest 0.1 m. Salinity and d18O data within this layer revealed that saline water originated from the surface brine of sea ice in 20% of the total sites and from seawater in 80%. From the vertical profiles of snow density, bulk thermal conductivity of snow was estimated as 0.15 W/K/m on average, only half of the value used for numerical sea-ice models. Although the upward heat flux within snow estimated with this value was significantly lower than that within ice, it turned out that a higher value of thermal conductivity (0.3 to 0.4 W/K/m) is preferable for estimating ice growth amount in current numerical models. Diurnal measurements showed that upward conductive heat flux within the snow and net long-wave radiation at the surface seem to play important roles in the formation of snow ice from slush. The detailed surface topography allowed us to compare the air-ice drag coefficients of ice and snow surfaces under neutral conditions, and to examine the possibility of the retrieval of ice thickness distribution from satellite remote sensing. It was found that overall snow cover works to enhance the surface roughness of sea ice rather than moderate it, and increases the drag coefficient by about 10%. As for thickness retrieval, mean ice thickness had a higher correlation with ice surface roughness than mean freeboard or surface elevation, which indicates the potential usefulness of satellite L-band SAR in estimating the ice thickness distribution in the seasonal sea-ice zone.