Use of Land Surface Temperature to Estimate Surface Energy Fluxes: Contributions of Wilfried Brutsaert and Collaborators

Land surface temperature (LST) plays a key role in governing the land surface energy budget, and measurements or estimates of LST are an integral part of many land surface models and methods to estimate land surface sensible heat (H) and latent heat fluxes. In particular, the LST anchors the potential temperature profile in Monin-Obukhov similarity theory, from which H can be derived. Brutsaert has made important contributions to our understanding the nature of surface temperature measurements as well as the practical but theoretically sound use of LST in this framework. His work has coincided with the wide-spread availability of remotely sensed LST measurements. Use of remotely sensed LST estimates inevitably involves complicating factors, such as: varying spatial and temporal scales in measurements, theory, and models; spatial variability of LST and H; the relationship between measurements of LST and the temperature felt by the atmosphere; and the need to correct satellite-based radiometric LST measurements for the radiative effects of the atmosphere. This paper reviews the progress made in research in these areas by tracing and commenting on Brutsaert's contributions.

Performance characteristics and design recommendations for biomass-burning stoves using earthen construction materials

This report provides an analysis of the thermal performance and emissions characteristics of improved biomass stoves constructed using earthen materials. Commonly referred to as mud stoves, this type of improved stove incorporates high clay content soil with an organic binder in the construction of its combustion chamber and body. When large quantities of the mud material are used to construct the stove body, the stove does not offer significant improvements in fuel economy or air quality relative to traditional open fire cooking. This is partly because a significant amount of heat is absorbed by the mass of the stove reducing combustion efficiency and heat transfer to the cook pot. An analysis of the thermal and mechanical properties of stove materials was also performed. A material mixture containing a one‐to‐one ratio by volume of high content clay soil and straw was found to have thermal properties comparable to fired ceramics used in more advanced improved stove designs. Feedback from mud stove users in Mauritania and Mali, West Africa was also collected during implementation. Suggestions for stove design improvements were developed based on this information and the data collected in the performance, emissions, and material properties analysis. Design suggestions include reducing stove height to accommodate user cooking preferences and limiting overall stove mass to reduce heat loss to the stove body.

Open top chambers and infrared lamps : a comparison of heating efficacy and CO₂/CH₄ dynamics in a Lake Superior coastal peatland

Experimental warming provides a method to determine how an ecosystem will respond to increased temperatures. Northern peatland ecosystems, sensitive to changing climates, provide an excellent setting for experimental warming. Storing great quantities of carbon, northern peatlands play a critical role in regulating global temperatures. Two of the most common methods of experimental warming include open top chambers (OTCs) and infrared (IR) lamps. These warming systems have been used in many ecosystems throughout the world, yet their efficacy to create a warmer environment is variable and has not been widely studied. To date, there has not been a direct, experimentally controlled comparison of OTCs and IR lamps. As a result, a factorial study was implemented to compare the warming efficacy of OTCs and IR lamps and to examine the resulting carbon dioxide (CO2) and methane (CH4) flux rates in a Lake Superior peatland. IR lamps warmed the ecosystem on average by 1-2 #°C, with the majority of warming occurring during nighttime hours. OTC's did not provide any long-term warming above control plots, which is contrary to similar OTC studies at high latitudes. By investigating diurnal heating patterns and micrometeorological variables, we were able to conclude that OTCs were not achieving strong daytime heating peaks and were often cooler than control plots during nighttime hours. Temperate day-length, cloudy and humid conditions, and latent heat loss were factors that inhibited OTC warming. There were no changes in CO2 flux between warming treatments in lawn plots. Gross ecosystem production was significantly greater in IR lamp-hummock plots, while ecosystem respiration was not affected. CH4 flux was not significantly affected by warming treatment. Minimal daytime heating differences, high ambient temperatures, decay resistant substrate, as well as other factors suppressed significant gas flux responses from warming treatments.

COMBUSTION INITIATION BY ELECTRICAL-DISCHARGE-INDUCED PLASMA IN LEAN AND DILUTE METHANE-AIR MIXTURE: EXPERIMENTAL AND MODELING INVESTIGATION

This dissertation represents experimental and numerical investigations of combustion initiation trigged by electrical-discharge-induced plasma within lean and dilute methane air mixture. This research topic is of interest due to its potential to further promote the understanding and prediction of spark ignition quality in high efficiency gasoline engines, which operate with lean and dilute fuel-air mixture. It is specified in this dissertation that the plasma to flame transition is the key process during the spark ignition event, yet it is also the most complicated and least understood procedure. Therefore the investigation is focused on the overlapped periods when plasma and flame both exists in the system. Experimental study is divided into two parts. Experiments in Part I focuses on the flame kernel resulting from the electrical discharge. A number of external factors are found to affect the growth of the flame kernel, resulting in complex correlations between discharge and flame kernel. Heat loss from the flame kernel to code ambient is found to be a dominant factor that quenches the flame kernel. Another experimental focus is on the plasma channel. Electrical discharges into gases induce intense and highly transient plasma. Detailed observation of the size and contents of the discharge-induced plasma channel is performed. Given the complex correlation and the multi-discipline physical/chemical processes involved in the plasma-flame transition, the modeling principle is taken to reproduce detailed transitions numerically with minimum analytical assumptions. Detailed measurement obtained from experimental work facilitates the more accurate description of initial reaction conditions. The novel and unique spark source considering both energy and species deposition is defined in a justified manner, which is the key feature of this Ignition by Plasma (IBP) model. The results of numerical simulation are intuitive and the potential of numerical simulation to better resolve the complex spark ignition mechanism is presented. Meanwhile, imperfections of the IBP model and numerical simulation have been specified and will address future attentions.

Two Stable Equilibria of the Atlantic Subpolar Gyre

The cyclonic circulation of the Atlantic subpolar gyre is a key mechanism for North Atlantic climate variability on a wide range of time scales. It is generally accepted that it is driven by both cyclonic winds and buoyancy forcing, yet the individual importance and dynamical interactions of the two contributions remain unclear. The authors propose a simplified four-box model representing the convective basin of the Labrador Sea and its shallow and deep boundary current system, the western subpolar gyre. Convective heat loss drives a baroclinic flow of relatively light water around the dense center. Eddy salt flux from the boundary current to the center increases with a stronger circulation, favors the formation of dense waters, and thereby sustains a strong baroclinic flow, approximately 10%–25% of the total. In contrast, when the baroclinic flow is not active, surface waters may be too fresh to convect, and a buoyancy-driven circulation cannot develop. This situation corresponds to a second stable circulation mode. A hysteresis is found for variations in surface freshwater flux and the salinity of the near-surface boundary current. An analytical solution is presented and analyzed.

Understanding the performance of the FLake model over two African Great Lakes

The ability of the one-dimensional lake model FLake to represent the mixolimnion temperatures for tropical conditions was tested for three locations in East Africa: Lake Kivu and Lake Tanganyika's northern and southern basins. Meteorological observations from surrounding automatic weather stations were corrected and used to drive FLake, whereas a comprehensive set of water temperature profiles served to evaluate the model at each site. Careful forcing data correction and model configuration made it possible to reproduce the observed mixed layer seasonality at Lake Kivu and Lake Tanganyika (northern and southern basins), with correct representation of both the mixed layer depth and water temperatures. At Lake Kivu, mixolimnion temperatures predicted by FLake were found to be sensitive both to minimal variations in the external parameters and to small changes in the meteorological driving data, in particular wind velocity. In each case, small modifications may lead to a regime switch, from the correctly represented seasonal mixed layer deepening to either completely mixed or permanently stratified conditions from similar to 10 m downwards. In contrast, model temperatures were found to be robust close to the surface, with acceptable predictions of near-surface water temperatures even when the seasonal mixing regime is not reproduced. FLake can thus be a suitable tool to parameterise tropical lake water surface temperatures within atmospheric prediction models. Finally, FLake was used to attribute the seasonal mixing cycle at Lake Kivu to variations in the near-surface meteorological conditions. It was found that the annual mixing down to 60m during the main dry season is primarily due to enhanced lake evaporation and secondarily to the decreased incoming long wave radiation, both causing a significant heat loss from the lake surface and associated mixolimnion cooling.

Dependence of Eemian Greenland temperature reconstructions on the ice sheet topography

The influence of a reduced Greenland Ice Sheet (GrIS) on Greenland's surface climate during the Eemian interglacial is studied using a set of simulations with different GrIS realizations performed with a comprehensive climate model. We find a distinct impact of changes in the GrIS topography on Greenland's surface air temperatures (SAT) even when correcting for changes in surface elevation, which influences SAT through the lapse rate effect. The resulting lapse-rate-corrected SAT anomalies are thermodynamically driven by changes in the local surface energy balance rather than dynamically caused through anomalous advection of warm/cold air masses. The large-scale circulation is indeed very stable among all sensitivity experiments and the Northern Hemisphere (NH) flow pattern does not depend on Greenland's topography in the Eemian. In contrast, Greenland's surface energy balance is clearly influenced by changes in the GrIS topography and this impact is seasonally diverse. In winter, the variable reacting strongest to changes in the topography is the sensible heat flux (SHF). The reason is its dependence on surface winds, which themselves are controlled to a large extent by the shape of the GrIS. Hence, regions where a receding GrIS causes higher surface wind velocities also experience anomalous warming through SHF. Vice-versa, regions that become flat and ice-free are characterized by low wind speeds, low SHF, and anomalous low winter temperatures. In summer, we find surface warming induced by a decrease in surface albedo in deglaciated areas and regions which experience surface melting. The Eemian temperature records derived from Greenland proxies, thus, likely include a temperature signal arising from changes in the GrIS topography. For the Eemian ice found in the NEEM core, our model suggests that up to 3.1 °C of the annual mean Eemian warming can be attributed to these topography-related processes and hence is not necessarily linked to large-scale climate variations.

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.

Daily polynya area and ice production in the North Water polynya region

The North Water (NOW) Polynya is a regularly-forming area of open-water and thin-ice, located between northwestern Greenland and Ellesmere Island (Canada) at the northern tip of Baffin Bay. Due to its large spatial extent, it is of high importance for a variety of physical and biological processes, especially in wintertime. Here, we present a long-term remote sensing study for the winter seasons 1978/1979 to 2014/2015. Polynya characteristics are inferred from (1) sea ice concentrations and brightness temperatures from passive microwave satellite sensors (Advanced Microwave Scanning Radiometer (AMSR-E and AMSR2), Scanning Multichannel Microwave Radiometer (SMMR), Special Sensor Microwave Imager/Sounder (SSM/I-SSMIS)) and (2) thin-ice thickness distributions, which are calculated using MODIS ice-surface temperatures and European Center for Medium-Range Weather Forecasts (ECMWF) atmospheric reanalysis data in a 1D thermodynamic energy-balance model. Daily ice production rates are retrieved for each winter season from 2002/2003 to 2014/2015, assuming that all heat loss at the ice surface is balanced by ice growth. Two different cloud-cover correction schemes are applied on daily polynya area and ice production values to account for cloud gaps in the MODIS composites. Our results indicate that the NOW polynya experienced significant seasonal changes over the last three decades considering the overall frequency of polynya occurrences, as well as their spatial extent. In the 1980s, there were prolonged periods of a more or less closed ice cover in northern Baffin Bay in winter. This changed towards an average opening on more than 85% of the days between November and March during the last decade. Noticeably, the sea ice cover in the NOW polynya region shows signs of a later-appearing fall freeze-up, starting in the late 1990s. Different methods to obtain daily polynya area using passive microwave AMSR-E/AMSR2 data and SSM/I-SSMIS data were applied. A comparison with MODIS data (thin-ice thickness < 20 cm) shows that the wintertime polynya area estimates derived by MODIS are about 30 to 40% higher than those derived using the polynya signature simulation method (PSSM) with AMSR-E data. In turn, the difference in polynya area between PSSM and a sea ice concentration (SIC) threshold of 70% is fairly low (approximately 10%) when applied to AMSR-E data. For the coarse-resolution SSM/I-SSMIS data, this difference is much larger, particularly in November and December. Instead of a sea ice concentration threshold, the PSSM method should be used for SSM/I-SSMIS data. Depending on the type of cloud-cover correction, the calculated ice production based on MODIS data reaches an average value of 264.4 ± 65.1 km**3 to 275.7 ± 67.4 km**3 (2002/2003 to 2014/2015) and shows a high interannual variability. Our achieved long-term results underline the major importance of the NOW polynya considering its influence on Arctic ice production and associated atmosphere/ocean processes.

Turbulent flux and meteorological measurements during ARCTEX-2006 campaign

This dataset present result from the DFG- funded Arctic-Turbulence-Experiment (ARCTEX-2006) performed by the University of Bayreuth on the island of Svalbard, Norway, during the winter/spring transition 2006. From May 5 to May 19, 2006 turbulent flux and meteorological measurements were performed on the monitoring field near Ny-Ålesund, at 78°55'24'' N, 11°55'15'' E Kongsfjord, Svalbard (Spitsbergen), Norway. The ARCTEX-2006 campaign site was located about 200 m southeast of the settlement on flat snow covered tundra, 11 m to 14 m above sea level. The permanent sites used for this study consisted of the 10 m meteorological tower of the Alfred Wegener Institute for Polar- and Marine Research (AWI), the international standardized radiation measurement site of the Baseline Surface Radiation Network (BSRN), the radiosonde launch site and the AWI tethered balloon launch sites. The temporary sites - set up by the University of Bayreuth - were a 6 m meteorological gradient tower, an eddy-flux measurement complex (EF), and a laser-scintillometer section (SLS). A quality assessment and data correction was applied to detect and eliminate specific measurement errors common at a high arctic landscape. In addition, the quality checked sensible heat flux measurements are compared with bulk aerodynamic formulas that are widely used in atmosphere-ocean/land-ice models for polar regions as described in Ebert and Curry (1993, doi:10.1029/93JC00656) and Launiainen and Cheng (1995). These parameterization approaches easily allow estimation of the turbulent surface fluxes from routine meteorological measurements. The data show: - the role of the intermittency of the turbulent atmospheric fluctuation of momentum and scalars, - the existence of a disturbed vertical temperature profile (sharp inversion layer) close to the surface, - the relevance of possible free convection events for the snow or ice melt in the Arctic spring at Svalbard, and - the relevance of meso-scale atmospheric circulation pattern and air-mass advection for the near-surface turbulent heat exchange in the Arctic spring at Svalbard. Recommendations and improvements regarding the interpretation of eddy-flux and laser-scintillometer data as well as the arrangement of the instrumentation under polar distinct exchange conditions and (extreme) weather situations could be derived.

Strontium chemistry of anhydrite from DSDP/ODP Hole 504B

A unique record of the chemical evolution of seawater during hydrothermal recharge into oceanic crust is preserved by anhydrite from the volcanic sequences and sheeted dike complex in ODP Hole 504B. Chemical and isotopic analyses 87Sr/86Sr, delta18O, delta34S of anhydrite constrain the changing composition of fluids due to reaction with basalt. There is a general trend of decreasing 87Sr/86Sr of anhydrite, corresponding to the minor incorporation of basaltic strontium with depth in the volcanic rocks. 87Sr/86Sr ratios decrease rapidly with depth in the dikes to values identical to host basalt (0.7029). Sr/Ca ratios (<0.1 mmol/mol) suggest that recharge fluids have very low Sr concentrations and fluids evolve by first precipitating Sr-bearing phases before extensive exchange of Sr with the host basalt. There is a background trend of decreasing sulfate delta18O with depth from +12-13? in the lower volcanics to +7? in the lower sheeted dikes recording an increase in recharge fluid temperature from c. 150° to c. 250°C, and confirming the presence of sulfate in hydrothermal fluids at elevated temperatures. From the amount of anhydrite recovered from Hole 504B and the amount of seawater sulfur that has been reduced to sulfide, a minimum seawater recharge flux can be calculated. This value is 4-25 times lower than estimates of high-temperature fluid fluxes based on either thermal constraints or global chemical budgets and suggests that there is significant deficit of seawater-derived sulfur in the oceanic crust. Only a minor proportion of the seawater that percolates into the crust near the axis is heated to high temperatures and exits as black smoker-type fluids. A significant proportion of the axial heat loss must be advected at 200-250°C by sulfate-bearing hydrothermal solutions that egress diffusely from the crust. These fluids penetrate into the dikes and exchange both heat and chemical tracers without the extensive clogging of porosity by anhydrite precipitation, which would halt hydrothermal circulation for any reasonable fluid flux. The heating of the major proportion of hydrothermal fluids to only moderate temperatures (c. 250°C) reconciles estimates of hydrothermal fluxes derived from thermal models and global geochemical budgets. The flux of hydrothermal sulfate would be of a magnitude similar to the riverine input, and oxygen-isotopic exchange at 200-250°C between dissolved sulfate and recharge fluids during hydrothermal circulation provides a mechanism to continuously buffer seawater sulfate oxygen to the light isotopic composition observed.

Composition of sediment pore water and isotopes in carbonate and benthic foraminifera at DSDP Holes 68-501, 69-505 and Site 69-504

Two sites on the southern flank of the Costa Rica Rift were drilled on DSDP Legs 68 and 69, one on crust 3.9 m.y. old and the other on crust 5.9 m.y. old. The basement of the younger site is effectively cooled by the circulation of seawater. The basement of the older site has been sealed by sediment, and an interval in the uppermost 560 meters of basement recently reheated to temperatures of 60 to 120°C. Although the thickness of the sediments at the two sites is similar (150-240 m versus 270 m), the much rougher basement topography at the younger Site 505 produces occasional basement outcrops, through which 80 to 90% of the total heat loss apparently occurs by advection of warm seawater. This seawater has been heated only slightly, however; the temperature at the base of the sediments is only 9°C. Changes in its composition due to reaction with the basement basalts are negligible, as indicated by profiles of sediment pore water chemistry. Bacterial sulfate reduction in the sediments produces a decrease in SO4 (and Ca) and an increase in alkalinity (and Sr and NH3) as depth increases to an intermediate level, but at deeper levels these trends reverse, and all of these species plus Mg, K, Na, and chlorinity approach seawater values near basement. Si, however, is higher, and Li may be lower. At the older site, Site 501/504, where heat loss is entirely by conduction, the temperature at the sediment/basement contact is 59°C. Sediment pore water chemistry is heavily affected by reaction with the basaltic basement, as indicated by large decreases in d18O, Mg, alkalinity, Na, and K and an increase in Ca with increasing depth. The size of the changes in d18O, Mg, alkalinity, Ca, Sr, and SO4 varies laterally over 500 meters, indicating lateral gradients in pore water chemistry that are nearly as large as the vertical gradients. The lateral gradients are believed to result from similar lateral gradients in the composition of the basement formation water, which propagate upward through the sediments by diffusion. A model of the d18O profile suggests that the basement at Site 501/504 was sealed off from advection about 1 m.y. ago, so that reaction rates began to dominate the basement pore water chemistry. A limestone-chert diagenetic front began to move upward through the lower sediments less than 200,000 yr. ago.

(Table 2) Oxygen, carbon, and strontium isotopic composition of carbonate-brucite material from the Lost City hydrothermal field

The isotopic (dD, d18O, d13C, and 87Sr/86Sr) and geochemical characteristics of hydrothermal solutions from the Mid-Atlantic Ridge and the material of brucite-carbonate chimneys at the Lost City hydrothermal field at 30°N, MAR, were examined to assay the role of the major factors controlling the genesis of the fluid and hydrothermal chimneys of the Lost City field. The values of dD and d18O in fluid samples indicates that solutions at the Lost City field were produced during the serpentinization of basement ultramafic rocks at temperatures higher than 200°C and at relatively low fluid/rock ratios (<1). The active role of serpentinization processes in the genesis of the Lost City fluid also follows from the results of the electron-microscopic studying of the material of hydrothermal chimneys at this field. The isotopic (d18O, d13C, and 87Sr/86Sr) and geochemical (Sr/Ca and REE) signatures indicate that, before its submarine discharging at the Lost City field, the fluid filtered through already cold altered outer zones of the Atlantis Massif and cooled via conductive heat loss. During this stage, the fluid could partly dissolve previously deposited carbonates in veins cutting serpentinite at the upper levels of the Atlantis Massif and the carbonate cement of sedimentary breccias underlying the hydrothermal chimneys. Because of this, the age of modern hydrothermal activity at the Lost City field can be much younger than 25 ka.

Stable isotope record of foraminifera from the subarctic Pacific

Benthic (Uvigerina spp., Cibicidoides spp., Gyroidinoides spp.) and planktonic (N. pachyderma sinistral, G. bulloides) stable isotope records from three core sites in the central Gulf of Alaska are used to infer mixed-layer and deepwater properties of the late glacial Subarctic Pacific. Glacial-interglacial amplitudes of the planktonic delta18O records are 1.1-1.3 per mil, less than half the amplitude observed at core sites at similar latitudes in the North Atlantic; these data imply that a strong, negative deltaw anomaly existed in the glacial Subarctic mixed layer during the summer, which points to a much stronger low-salinity anomaly than exists today. If true, the upper water column in the North Pacific would have been statically more stable than today, thus suppressing convection even more efficiently. This scenario is further supported by vertical (i.e., planktic versus benthic) delta18O and delta13C gradients of >1 per mil, which suggest that a thermohaline link between Pacific deep waters and the Subarctic Pacific mixed layer did not exist during the late glacial. Epibenthic delta13C in the Subarctic Pacific is more negative than at tropical-subtropical Pacific sites but similar to that recorded at Southern Ocean sites, suggesting ventilation of the deep central Pacific from mid-latitude sources, e.g., from the Sea of Japan and Sea of Okhotsk. Still, convection to intermediate depths could have occurred in the Subarctic during the winter months when heat loss to the atmosphere, sea ice formation, and wind-driven upwelling of saline deep waters would have been most intense. This would be beyond the grasp of our planktonic records which only document mixed-layer temperature-salinity fields extant during the warmer seasons. Also we do not have benthic isotope records from true intermediate water depths of the Subarctic Pacific.