Spatial characterization of land surface energy fluxes and uncertainty estimation at the Salar de Atacama,Northern Chile

We use Advanced Spaceborne Thermal Emission and Reflection radiometer (ASTER) data to estimate spatial energy flux and evaporation distributions at the Salar de Atacama, a playa in Northern Chile. Our approach incorporates ASTER surface kinetic temperature, emissivity, and reflectance data, ground-based meteorological measurements, and empirical parameters. Energy flux distributions are estimated using either spatially constant or spatially distributed values of model parameters, with spatially distributed parameters assigned separately to each land cover category in an image classification. We test the sensitivity of energy budget calculations to state variable and parameter values by conducting Monte Carlo simulations for regions with ground energy budget measurements. Results show that assigning spatially distributed model parameters via land cover classifications yields significant improvements to ground and sensible heat flux predictions. Latent heat fluxes cannot, however, be predicted with sufficient accuracy to allow estimation of area-integrated evaporative moisture loss at this low-evaporation playa.     

A threshold in soil formation at Earth’s arid-hyperarid transition

The soils of the Atacama Desert in northern Chile have long been known to contain large quantities of unusual salts, yet the processes that form these soils are not yet fully understood. We examined the morphology and geochemistry of soils on post-Miocene fans and stream terraces along a south-to-north (27° to 24° S) rainfall transect that spans the arid to hyperarid transition (21 to ∼2 mm rain y⁻¹). Landform ages are ? 2 My based on cosmogenic radionuclide concentrations in surface boulders, and Ar isotopes in interbedded volcanic ash deposits near the driest site indicate a maximum age of 2.1 My. A chemical mass balance analysis that explicitly accounts for atmospheric additions was used to quantify net changes in mass and volume as a function of rainfall. In the arid (21 mm rain y⁻¹) soil, total mass loss to weathering of silicate alluvium and dust (−1030 kg m⁻²) is offset by net addition of salts (+170 kg m⁻²). The most hyperarid soil has accumulated 830 kg m⁻² of atmospheric salts (including 260 kg sulfate m⁻² and 90 kg chloride m⁻²), resulting in unusually high volumetric expansion (120%) for a soil of this age. The composition of both airborne particles and atmospheric deposition in passive traps indicates that the geochemistry of the driest soil reflects accumulated atmospheric influxes coupled with limited in-soil chemical transformation and loss. Long-term rates of atmospheric solute addition were derived from the ion inventories in the driest soil, divided by the landform age, and compared to measured contemporary rates. With decreasing rainfall, the soil salt inventories increase, and the retained salts are both more soluble and present at shallower depths. All soils generally exhibit vertical variation in their chemistry, suggesting slow and stochastic downward water movement, and greater climate variability over the past 2 My than is reflected in recent (∼100 y) rainfall averages. The geochemistry of these soils shows that the transition from arid to hyperarid rainfall levels marks a fundamental geochemical threshold: in wetter soils, the rate and character of chemical weathering results in net mass loss and associated volumetric collapse after 10⁵ to 10⁶ years, while continuous accumulation of atmospheric solutes in hyperarid soils over similar timescales results in dramatic volumetric expansion. The specific geochemistry of hyperarid soils is a function of atmospheric sources, and is expected to vary accordingly at other hyperarid sites. This work identifies key processes in hyperarid soil formation that are likely to be independent of location, and suggests that analogous processes may occur on Mars.     

Controls on nutrient and phytoplankton dynamics during normal flow and storm runoff conditions, southern Kaneohe Bay, Hawaii

Fluvial effects on nutrient and phytoplankton dynamics were evaluated in southern Kaneohe Bay, Oahu, Hawaii. Fluvial inputs occurred as small, steady baseflows interrupted by intense pulses of storm runoff. Baseflow river inputs only affected restricted areas around stream mouths, but the five storm events sampled during this study produced transient runoff plumes of much greater spatial extent. Nutrient loading via runoff generally led to an increase of the phytoplankton biomass and gross primary productivity in southern Kaneohe Bay, but the rapid depletion of nutrients resulted in a decline of the algal populations in the relatively short time of days. Under baseline conditions, water column primary productivity in southern Kaneohe Bay is normally nitrogen limited. Following storm events, the high ratio of dissolved inorganic nitrogen to dissolved inorganic phosphorus (DIN:DIP, 25–29) fluxes of runoff nutrients drove bay waters towards phosphorus limitation. A depletion of phosphate relative to DIN in surface waters was observed following all storm events. Due to high flushing rates, recovery times of bay waters from storm perturbations ranged from 3 to 8 d and appeared to be correlated with tidal range. Storm inputs have a significant effect on the water column ecosystem and biogeochemistry in southern Kaneohe Bay. The perturbations were only transient events and the system rapidly recovered to prestorm conditions.     

Evaporation and land surface energy budget at the Salar de Atacama, Northern Chile

Playa systems are driven by evaporation processes, yet the mechanisms by which evaporation occurs through playa salt crusts are still poorly understood. In this study we examine playa evaporation as it relates to land surface energy fluxes, salt crust characteristics, groundwater and climate at the Salar de Atacama, a 3000 km² playa in northern Chile containing a uniquely broad range of salt crust types. Land surface energy budget measurements were taken at eight representative sites on this playa during winter (August 2001) and summer (January 2002) seasons. Measured values of net all-wave radiation were highest at vegetated and rough halite crust sites and lowest over smooth, highly reflective salt crusts. Over most of the Salar de Atacama, net radiation was dissipated by means of soil and sensible heat fluxes. Dry salt crusts tended to heat and cool very quickly, whereas soil heating and cooling occurred more gradually at wetter vegetated sites. Sensible heating was strongly linked to wind patterns, with highest sensible heat fluxes occurring on summer days with strong afternoon winds. Very little energy available at the land surface was used to evaporate water. Eddy covariance measurements could only constrain evaporation rates to within 0.1 mm d⁻¹, and some measured evaporation rates were less than this margin of uncertainty. Evaporation rates ranged from 0.1 to 1.1 mm d⁻¹ in smooth salt crusts around the margin of the salar and from 0.4 to 2.8 mm d⁻¹ in vegetated areas. No evaporation was detected from the rugged halite salt crust that covers the interior of the salar, though the depth to groundwater is less than 1 m in this area. These crusts therefore represent a previously unrecorded end member condition in which the salt crusts form a practically impermeable barrier to evaporation.     

A numerical study on the influence of geometry on the formation of sandbanks

Applying a two-dimensional, non-linear hydrodynamic numerical model in combination with a semiempirical equation for bedload sediment transport, the influence of geometry on the formation of sandbanks is investigated. In the first experiment, the formation of sandbanks in an ideal rectangular basin, resembling the Taylor’s problem, was calculated. Sandbanks occur in a small area at the closed boundary. Similar experiments were carried out for a range of wavelengths of the incident Kelvin wave. The results reveal that large wavelengths favor the generation of sandbanks. In subsequent calculations, the basin was modified by introducing new geographical features like bays and peninsulas. The numerical experiments show that geometry is a fundamental factor to determine the position where groups of sandbanks appear. The results suggest that in regions where the Kelvin wave is diffracted, the formation of sandbanks occurs. An experiment, in which we applied an ideal geometric configuration representing that of the Southern Bight of the North Sea, generated sandbank patterns resembling those observed in the region.     

The Madden–Julian oscillation in ECHAM4 coupled and uncoupled general circulation models

The Madden-Julian oscillation (MJO) dominates tropical variability on timescales of 30–70 days. During the boreal winter/spring, it is manifested as an eastward propagating disturbance, with a strong convective signature over the eastern hemisphere. The space–time structure of the MJO is analyzed using simulations with the ECHAM4 atmospheric general circulation model run with observed monthly mean sea-surface temperatures (SSTs), and coupled to three different ocean models. The coherence of the eastward propagation of MJO convection is sensitive to the ocean model to which ECHAM4 is coupled. For ECHAM4/OPYC and ECHO-G, models for which ~100 years of daily data is available, Monte Carlo sampling indicates that their metrics of eastward propagation are different at the 1% significance level. The flux-adjusted coupled simulations, ECHAM4/OPYC and ECHO-G, maintain a more realistic mean-state, and have a more realistic MJO simulation than the nonadjusted scale interaction experiment (SINTEX) coupled runs. The SINTEX model exhibits a cold bias in Indian Ocean and tropical West Pacific Ocean sea-surface temperature of ~0.5°C. This cold bias affects the distribution of time-mean convection over the tropical eastern hemisphere. Furthermore, the eastward propagation of MJO convection in this model is not as coherent as in the two models that used flux adjustment or when compared to an integration of ECHAM4 with prescribed observed SST. This result suggests that simulating a realistic basic state is at least as important as air–sea interaction for organizing the MJO. While all of the coupled models simulate the warm (cold) SST anomalies that precede (succeed) the MJO convection, the interaction of the components of the net surface heat flux that lead to these anomalies are different over the Indian Ocean. The ECHAM4/OPYC model in which the atmospheric model is run at a horizontal resolution of T42, has eastward propagating zonal wind anomalies and latent heat flux anomalies. However, the integrations with ECHO-G and SINTEX, which used T30 atmospheres, produce westward propagation of the latent heat flux anomalies, contrary to reanalysis. It is suggested that the differing ability of the models to represent the near-surface westerlies over the Indian Ocean is related to the different horizontal resolutions of the atmospheric model employed.     

A Full Waveform Test of the Southern California Velocity Model by the Reciprocity Method

— We apply the reciprocity method (Eisner and Clayton, 2001a) to compare the full waveform synthetic seismograms with a large number of observed seismograms. The reciprocity method used in the finite-difference modeling allows for the use of high quality data observed from the earthquakes distributed over the wide range of azimuths and depths. We have developed a methodology to facilitate the comparison between data and synthetics using a set of attributes to characterize the seismograms. These attributes are maximum amplitude, time delay and coda decay of the magnitude of the displacement vector. For the Southern California Velocity Model, Version 1 (Magistrale et al., 1996), we have found misfits between data and synthetics for paths traveling outside of the sedimentary basins and the western part of the Los Angeles and San Fernando basins.     

Terminal Pleistocene/Early Holocene Environmental Change at the Sunshine Locality, North-Central Nevada, U.S.A.

Sedimentological, faunal, and archaeological investigations at the Sunshine Locality, Long Valley, Nevada reveal a history of human adaptation and environmental change at the last glacial–interglacial transition in North America's north-central Great Basin. The locality contains a suite of lacustrine, alluvial, and eolian deposits associated with fluvially reworked faunal remains and Paleoindian artifacts. Radiocarbon-dated stratigraphy indicates a history of receding pluvial lake levels followed by alluvial downcutting and subsequent valley filling with marsh-like conditions at the end of the Pleistocene. A period of alluvial deposition and shallow water tables (9,800 to 11,000 ¹⁴C yr B.P.) correlates to the Younger Dryas. Subsequent drier conditions and reduced surface runoff mark the early Holocene; sand dunes replace wetlands by 8,000 ¹⁴C yr B.P. The stratigraphy at Sunshine is similar to sites located 400 km south and supports regional climatic synchroneity in the central and southern Great Basin during the terminal Pleistocene/early Holocene. Given regional climate change and recurrent geomorphic settings comparable to Sunshine, we believe that there is a high potential for buried Paleoindian features in primary association with extinct fauna elsewhere in the region yet to be discovered due to limited stratigraphic exposure and consequent low visibility.     

Validation of ADEOS-II GLI ocean color products using in-situ observations

The Global Imager (GLI) aboard the Advanced Earth Observing Satellite-II (ADEOS-II) made global observations from 2 April 2003 to 24 October 2003. In cooperation with several institutes and scientists, we obtained quality controlled match-ups between GLI products and in-situ data, 116 for chlorophyll-a concentration (CHLA), 249 for normalized water-leaving radiance (nLw) at 443 nm, and 201 for aerosol optical thickness at 865 nm (Tau_865) and Angstrom exponent between 520 and 865 nm (Angstrom). We evaluated the GLI ocean color products and investigated the causes of errors using the match-ups. The median absolute percentage differences (MedPD) between GLI and in-situ data were 14.1–35.7% for nLws at 380–565 nm, 52.5–74.8% nLws at 625–680 nm, 47.6% for Tau_865, 46.2% for Angstrom, and 46.6% for CHLA, values that are comparable to the ocean-color products of other sensors. We found that some errors in GLI products are correlated with observational conditions; nLw values were underestimated when nLw at 680 nm was high, CHLA was underestimated in absorptive aerosol conditions, and Tau_865 was overestimated in sunglint regions. The error correlations indicate that we need to improve the retrievals of the optical properties of absorptive aerosols and seawater and sea surface reflection for further applications, including coastal monitoring and the combined use of products from multiple sensors.     

Interactive Videos Enhance Learning about Socio-Ecological Systems

Two forms of interactive video were assessed in an online course focused on conservation. The hypothesis was that interactive video enhances student perceptions about learning and improves mental models of social-ecological systems. Results showed that students reported greater learning and attitudes toward the subject following interactive video. Moreover, metrics of network complexity applied to students' concept maps of socio-ecological systems increased longitudinally through the course highlighting greater awareness of socio-ecological linkages. We conclude that interactive video can increase awareness of interdisciplinary connections in socio-ecological systems and should be more widely employed in general education geography courses.