A model study of oceanic mechanisms affecting Equatorial Pacific sea surface temperature during the 1997-98 El Niño

In this study, the processes affecting sea surface temperature variability over the 1992–98 period, encompassing the very strong 1997–98 El Niño event, are analyzed. A tropical Pacific Ocean general circulation model, forced by a combination of weekly ERS1–2 and TAO wind stresses, and climatological heat and freshwater fluxes, is first validated against observations. The model reproduces the main features of the tropical Pacific mean state, despite a weaker than observed thermal stratification, a 0.1 m s−1 too strong (weak) South Equatorial Current (North Equatorial Countercurrent), and a slight underestimate of the Equatorial Undercurrent. Good agreement is found between the model dynamic height and TOPEX/Poseidon sea level variability, with correlation/rms differences of 0.80/4.7 cm on average in the 10°N–10°S band. The model sea surface temperature variability is a bit weak, but reproduces the main features of interannual variability during the 1992–98 period. The model compares well with the TAO current variability at the equator, with correlation/rms differences of 0.81/0.23 m s−1 for surface currents. The model therefore reproduces well the observed interannual variability, with wind stress as the only interannually varying forcing. This good agreement with observations provides confidence in the comprehensive three-dimensional circulation and thermal structure of the model. A close examination of mixed layer heat balance is thus undertaken, contrasting the mean seasonal cycle of the 1993–96 period and the 1997–98 El Niño. In the eastern Pacific, cooling by exchanges with the subsurface (vertical advection, mixing, and entrainment), the atmospheric forcing, and the eddies (mainly the tropical instability waves) are the three main contributors to the heat budget. In the central–western Pacific, the zonal advection by low-frequency currents becomes the main contributor. Westerly wind bursts (in December 1996 and March and June 1997) were found to play a decisive role in the onset of the 1997–98 El Niño. They contributed to the early warming in the eastern Pacific because the downwelling Kelvin waves that they excited diminished subsurface cooling there. But it is mainly through eastward advection of the warm pool that they generated temperature anomalies in the central Pacific. The end of El Niño can be linked to the large-scale easterly anomalies that developed in the western Pacific and spread eastward, from the end of 1997 onward. In the far-western Pacific, because of the shallower than normal thermocline, these easterlies cooled the SST by vertical processes. In the central Pacific, easterlies pushed the warm pool back to the west. In the east, they led to a shallower thermocline, which ultimately allowed subsurface cooling to resume and to quickly cool the surface layer.

Effects of biosolids from tomato processing on soil fertility and wheat growth in a Greek alfisol

The effects of biosolids from tomato processing on soil properties and wheat growth were investigated in an Alfisol from central Greece. Biosolids were mixed with soil from the surface (Ap) or subsurface (Bt) horizon in plastic containers at rates of 1%, 5%, and 10% by dry weight (d.w.; equivalent to 10, 50, and 100 Mg ha–1). Biosolid treatments were compared to an NH4Cl application (50 mg N kg–1) and an untreated control in (1) a 102 d incubation experiment at 28°C to determine biosolid nitrification potential and (2) a 45 d outdoor experiment to evaluate effects on soil fertility and wheat growth. Mineralization of biosolids in the incubation experiment resulted in accumulation of nitrate-N and indicated that biosolids were able to supply N that was in excess of crop needs in treatments of 5% and 10%. After 45 d of wheat growth, available soil nutrients (N, P) and P uptake by wheat were distinctly lower in the Bt than in the Ap horizon. However, soil pH, electrical conductivity, organic matter, total N, nitrate-N, extractable P, and exchangeable K increased with increasing rate of biosolid application in both soils. These were followed by corresponding increases in wheat nutrient uptake and biomass production, thus demonstrating the importance of this organic material for sustaining production in soils of low immediate fertility. Compared to the NH4Cl treatment (50 kg N ha–1 equivalent), biosolid application rates of 5% and 10% had higher available soil nutrients, similar or higher nutrient uptake and higher wheat biomass. But only an application of 10% biosolids provided sufficient N levels for wheat in the surface soil, and even higher applications were required for providing sufficient N and P in the Bt horizon.

A model study of factors influencing projected changes in regional sea level over the 21st century

In addition to projected increases in global mean sea level over the 21st century, model simulations suggest there will also be changes in the regional distribution of sea level relative to the global mean. There is a considerable spread in the projected patterns of these changes by current models, as shown by the recent Intergovernmental Panel on Climate Change (IPCC) Fourth Assessment (AR4). This spread has not reduced from that given by the Third Assessment models. Comparison with projections by ensembles of models based on a single structure supports an earlier suggestion that models of similar formulation give more similar patterns of sea level change. Analysing an AR4 ensemble of model projections under a business-as-usual scenario shows that steric changes (associated with subsurface ocean density changes) largely dominate the sea level pattern changes. The relative importance of subsurface temperature or salinity changes in contributing to this differs from region to region and, to an extent, from model-to-model. In general, thermosteric changes give the spatial variations in the Southern Ocean, halosteric changes dominate in the Arctic and strong compensation between thermosteric and halosteric changes characterises the Atlantic. The magnitude of sea level and component changes in the Atlantic appear to be linked to the amount of Atlantic meridional overturning circulation (MOC) weakening. When the MOC weakening is substantial, the Atlantic thermosteric patterns of change arise from a dominant role of ocean advective heat flux changes.

Hyperresolution global land surface modeling: meeting a grand challenge for monitoring Earth's terrestrial water

Monitoring Earth's terrestrial water conditions is critically important to many hydrological applications such as global food production; assessing water resources sustainability; and flood, drought, and climate change prediction. These needs have motivated the development of pilot monitoring and prediction systems for terrestrial hydrologic and vegetative states, but to date only at the rather coarse spatial resolutions (∼10–100 km) over continental to global domains. Adequately addressing critical water cycle science questions and applications requires systems that are implemented globally at much higher resolutions, on the order of 1 km, resolutions referred to as hyperresolution in the context of global land surface models. This opinion paper sets forth the needs and benefits for a system that would monitor and predict the Earth's terrestrial water, energy, and biogeochemical cycles. We discuss six major challenges in developing a system: improved representation of surface‐subsurface interactions due to fine‐scale topography and vegetation; improved representation of land‐atmospheric interactions and resulting spatial information on soil moisture and evapotranspiration; inclusion of water quality as part of the biogeochemical cycle; representation of human impacts from water management; utilizing massively parallel computer systems and recent computational advances in solving hyperresolution models that will have up to 109 unknowns; and developing the required in situ and remote sensing global data sets. We deem the development of a global hyperresolution model for monitoring the terrestrial water, energy, and biogeochemical cycles a “grand challenge” to the community, and we call upon the international hydrologic community and the hydrological science support infrastructure to endorse the effort.

Coupled atmosphere-biophysics-hydrology models for environmental modeling

The formulation and implementation of LEAF-2, the Land Ecosystem–Atmosphere Feedback model, which comprises the representation of land–surface processes in the Regional Atmospheric Modeling System (RAMS), is described. LEAF-2 is a prognostic model for the temperature and water content of soil, snow cover, vegetation, and canopy air, and includes turbulent and radiative exchanges between these components and with the atmosphere. Subdivision of a RAMS surface grid cell into multiple areas of distinct land-use types is allowed, with each subgrid area, or patch, containing its own LEAF-2 model, and each patch interacts with the overlying atmospheric column with a weight proportional to its fractional area in the grid cell. A description is also given of TOPMODEL, a land hydrology model that represents surface and subsurface downslope lateral transport of groundwater. Details of the incorporation of a modified form of TOPMODEL into LEAF-2 are presented. Sensitivity tests of the coupled system are presented that demonstrate the potential importance of the patch representation and of lateral water transport in idealized model simulations. Independent studies that have applied LEAF-2 and verified its performance against observational data are cited. Linkage of RAMS and TOPMODEL through LEAF-2 creates a modeling system that can be used to explore the coupled atmosphere–biophysical–hydrologic response to altered climate forcing at local watershed and regional basin scales.

Contribution of N Species and P Fractions to Stream Water Quality in Agricultural Catchments

The contribution from agricultural catchments to stream nitrogen and phosphorus concentrations was assessed by evaluation of the chemical composition of these nutrients in agricultural runoff for both surface and subsurface flow pathways. A range of land uses (grazed and ungrazed grassland, cereals, roots) in intensive agricultural systems was studied at scales from hillslope plots (0.5m2) to large catchment (>300km2). By fractionating the total nutrient load it was possible to establish that most of the phosphorus was transported in the unreactive (particulate and organic) fraction via surface runoff. This was true regardless of the scale of measurement. The form of the nitrogen load varied with land use and grazing intensity. High loads of dissolved inorganic nitrogen (with >90% transported as NH4-N) were recorded in surface runoff from heavily grazed land. In subsurface flow from small (2km2) subcatchments and in larger (>300 km2) catchments, organic nitrogen was found to be an important secondary constituent of the total nitrogen load, comprising 40% of the total annual load.

Modelling the impact of land use change on water quality in agricultural catchments

Export coefficient modelling was used to model the impact of agriculture on nitrogen and phosphorus loading on the surface waters of two contrasting agricultural catchments. The model was originally developed for the Windrush catchment where the highly reactive Jurassic limestone aquifer underlying the catchment is well connected to the surface drainage network, allowing the system to be modelled using uniform export coefficients for each nutrient source in the catchment, regardless of proximity to the surface drainage network. In the Slapton catchment, the hydrological path-ways are dominated by surface and lateral shallow subsurface flow, requiring modification of the export coefficient model to incorporate a distance-decay component in the export coefficients. The modified model was calibrated against observed total nitrogen and total phosphorus loads delivered to Slapton Ley from inflowing streams in its catchment. Sensitivity analysis was conducted to isolate the key controls on nutrient export in the modified model. The model was validated against long-term records of water quality, and was found to be accurate in its predictions and sensitive to both temporal and spatial changes in agricultural practice in the catchment. The model was then used to forecast the potential reduction in nutrient loading on Slapton Ley associated with a range of catchment management strategies. The best practicable environmental option (BPEO) was found to be spatial redistribution of high nutrient export risk sources to areas of the catchment with the greatest intrinsic nutrient retention capacity.

The long-term stability of a possible aqueous ammonium sulfate ocean inside Titan

We model the thermal evolution of a subsurface ocean of aqueous ammonium sulfate inside Titan using a parameterized convection scheme. The cooling and crystallization of such an ocean depends on its heat flux balance, and is governed by the pressure-dependent melting temperatures at the top and bottom of the ocean. Using recent observations and previous experimental data, we present a nominal model which predicts the thickness of the ocean throughout the evolution of Titan; after 4.5 Ga we expect an aqueous ammonium sulfate ocean 56 km thick, overlain by a thick (176 km) heterogeneous crust of methane clathrate, ice I and ammonium sulfate. Underplating of the crust by ice I will give rise to compositional diapirs that are capable of rising through the crust and providing a mechanism for cryovolcanism at the surface. We have conducted a parameter space survey to account for possible variations in the nominal model, and find that for a wide range of plausible conditions, an ocean of aqueous ammonium sulfate can survive to the present day, which is consistent with the recent observations of Titan's spin state from Cassini radar data [Lorenz, R.D., Stiles, B.W., Kirk, R.L., Allison, M.D., del Marmo, P.P., Iess, L., Lunine, J.I., Ostro, S.J., Hensley, S., 2008. Science 319, 1649–1651].

Have aerosols caused the observed Atlantic multidecadal variability?

Identifying the prime drivers of the twentieth-century multidecadal variability in the Atlantic Ocean is crucial for predicting how the Atlantic will evolve in the coming decades and the resulting broad impacts on weather and precipitation patterns around the globe. Recently, Booth et al. showed that the Hadley Centre Global Environmental Model, version 2, Earth system configuration (HadGEM2-ES) closely reproduces the observed multidecadal variations of area-averaged North Atlantic sea surface temperature in the twentieth century. The multidecadal variations simulated in HadGEM2-ES are primarily driven by aerosol indirect effects that modify net surface shortwave radiation. On the basis of these results, Booth et al. concluded that aerosols are a prime driver of twentieth-century North Atlantic climate variability. However, here it is shown that there are major discrepancies between the HadGEM2-ES simulations and observations in the North Atlantic upper-ocean heat content, in the spatial pattern of multidecadal SST changes within and outside the North Atlantic, and in the subpolar North Atlantic sea surface salinity. These discrepancies may be strongly influenced by, and indeed in large part caused by, aerosol effects. It is also shown that the aerosol effects simulated in HadGEM2-ES cannot account for the observed anticorrelation between detrended multidecadal surface and subsurface temperature variations in the tropical North Atlantic. These discrepancies cast considerable doubt on the claim that aerosol forcing drives the bulk of this multidecadal variability.

Spatial interaction among nontoxic phytoplankton, toxic phytoplankton and zooplankton: emergence in space and time

In homogeneous environments, by overturning the possibility of competitive exclusion among phytoplankton species, and by regulating the dynamics of overall plankton population, toxin-producing phytoplankton (TPP) potentially help in maintaining plankton diversity—a result shown recently. Here, I explore the competitive effects of TPP on phytoplankton and zooplankton species undergoing spatial movements in the subsurface water. The spatial interactions among the species are represented in the form of reaction-diffusion equations. Suitable parametric conditions under which Turing patterns may or may not evolve are investigated. Spatiotemporal distributions of species biomass are simulated using the diffusivity assumptions realistic for natural planktonic systems. The study demonstrates that spatial movements of planktonic systems in the presence of TPP generate and maintain inhomogeneous biomass distribution of competing phytoplankton, as well as grazer zooplankton, thereby ensuring the persistence of multiple species in space and time. The overall results may potentially explain the sustainability of biodiversity and the spatiotemporal emergence of phytoplankton and zooplankton species under the influence of TPP combined with their physical movement in the subsurface water.

Influence of UV radiation on chlorophyll, and antioxidant enzymes of wetland plants in different types of constructed wetland

A surface- and vertical subsurface-flow-constructed wetland were designed to study the response of chlorophyll and antioxidant enzymes to elevated UV radiation in three types of wetland plants (Canna indica, Phragmites austrail, and Typha augustifolia). Results showed that (1) chlorophyll content of C. indica, P. austrail, and T. augustifolia in the constructed wetland was significantly lower where UV radiation was increased by 10 and 20 % above ambient solar level than in treatment with ambient solar UV radiation (p < 0.05). (2) The malondialdehyde (MDA) content, guaiacol peroxidase (POD), superoxide dismutase (SOD), and catalase (CAT) activities of wetland plants increased with elevated UV radiation intensity. (3) The increased rate of MDA, SOD, POD, and CAT activities of C. indica, P. australis, and T. angustifolia by elevated UV radiation of 10 % was higher in vertical subsurface-flow-constructed wetland than in surface-flow-constructed wetland. The sensitivity of MDA, SOD, POD, and CAT activities of C. indica, P. austrail, and T. augustifolia to the elevated UV radiation was lower in surface-flow-constructed wetland than in the vertical subsurface-flow-constructed wetland, which was related to a reduction in UV radiation intensity through the dissolved organic carbon and suspended matter in the water. C. indica had the highest SOD and POD activities, which implied it is more sensitive to enhanced UV radiation. Therefore, different wetland plants had different antioxidant enzymes by elevated UV radiation, which were more sensitive in vertical subsurface-flow-constructed wetland than in surface-flow-constructed wetland.

Multi‐method global sensitivity analysis (MMGSA) for modelling floodplain hydrological processes

When studying hydrological processes with a numerical model, global sensitivity analysis (GSA) is essential if one is to understand the impact of model parameters and model formulation on results. However, different definitions of sensitivity can lead to a difference in the ranking of importance of the different model factors. Here we combine a fuzzy performance function with different methods of calculating global sensitivity to perform a multi-method global sensitivity analysis (MMGSA). We use an application of a finite element subsurface flow model (ESTEL-2D) on a flood inundation event on a floodplain of the River Severn to illustrate this new methodology. We demonstrate the utility of the method for model understanding and show how the prediction of state variables, such as Darcian velocity vectors, can be affected by such a MMGSA. This paper is a first attempt to use GSA with a numerically intensive hydrological model.

Ocean heat uptake processes: a model intercomparison

We compare the quasi-equilibrium heat balances, as well as their responses to 4×CO2 perturbation, among three global climate models with the aim to identify and explain inter-model differences in ocean heat uptake (OHU) processes. We find that, in quasi-equilibrium, convective and mixed layer processes, as well as eddy-related processes, cause cooling of the subsurface ocean. The cooling is balanced by warming caused by advective and diapycnally diffusive processes. We also find that in the CO2-perturbed climates the largest contribution to OHU comes from changes in vertical mixing processes and the mean circulation, particularly in the extra-tropics, caused both by changes in wind forcing, and by changes in high-latitude buoyancy forcing. There is a substantial warming in the tropics, a significant part of which occurs because of changes in horizontal advection in extra-tropics. Diapycnal diffusion makes only a weak contribution to the OHU, mainly in the tropics, due to increased stratification. There are important qualitative differences in the contribution of eddy-induced advection and isopycnal diffusion to the OHU among the models. The former is related to the different values of the coefficients used in the corresponding scheme. The latter is related to the different tapering formulations of the isopycnal diffusion scheme. These differences affect the OHU in the deep ocean, which is substantial in two of the models, the dominant region of deep warming being the Southern Ocean. However, most of the OHU takes place above 2000 m, and the three models are quantitatively similar in their global OHU efficiency and its breakdown among processes and as a function of latitude.

Implication of methodological uncertainties for mid-Holocene sea surface temperature reconstructions

We present and examine a multi-sensor global compilation of mid-Holocene (MH) sea surface temperatures (SST), based on Mg/Ca and alkenone palaeothermometry and reconstructions obtained using planktonic foraminifera and organic-walled dinoflagellate cyst census counts. We assess the uncertainties originating from using different methodologies and evaluate the potential of MH SST reconstructions as a benchmark for climate-model simulations. The comparison between different analytical approaches (time frame, baseline climate) shows the choice of time window for the MH has a negligible effect on the reconstructed SST pattern, but the choice of baseline climate affects both the magnitude and spatial pattern of the reconstructed SSTs. Comparison of the SST reconstructions made using different sensors shows significant discrepancies at a regional scale, with uncertainties often exceeding the reconstructed SST anomaly. Apparent patterns in SST may largely be a reflection of the use of different sensors in different regions. Overall, the uncertainties associated with the SST reconstructions are generally larger than the MH anomalies. Thus, the SST data currently available cannot serve as a target for benchmarking model simulations. Further evaluations of potential subsurface and/or seasonal artifacts that may contribute to obscure the MH SST reconstructions are urgently needed to provide reliable benchmarks for model evaluations.

Multi-method global sensitivity analysis (MMGSA) for modelling floodplain hydrological processes

When studying hydrological processes with a numerical model, global sensitivity analysis (GSA) is essential if one is to understand the impact of model parameters and model formulation on results. However, different definitions of sensitivity can lead to a difference in the ranking of importance of the different model factors. Here we combine a fuzzy performance function with different methods of calculating global sensitivity to perform a multi-method global sensitivity analysis (MMGSA). We use an application of a finite element subsurface flow model (ESTEL-2D) on a flood inundation event on a floodplain of the River Severn to illustrate this new methodology. We demonstrate the utility of the method for model understanding and show how the prediction of state variables, such as Darcian velocity vectors, can be affected by such a MMGSA. This paper is a first attempt to use GSA with a numerically intensive hydrological model