Evaluation of a Photosynthesis-Based Canopy Resistance Formulation in the Noah Land-Surface Model

Accurately representing complex land-surface processes balancing complexity and realism remains one challenge that the weather modelling community is facing nowadays. In this study, a photosynthesis-based Gas-exchange Evapotranspiration Model (GEM) is integrated into the Noah land-surface model replacing the traditional Jarvis scheme for estimating the canopy resistance and transpiration. Using 18-month simulations from the High Resolution Land Data Assimilation System (HRLDAS), the impact of the photosynthesis-based approach on the simulated canopy resistance, surface heat fluxes, soil moisture, and soil temperature over different vegetation types is evaluated using data from the Atmospheric Radiation Measurement (ARM) site, Oklahoma Mesonet, 2002 International H₂O Project (IHOP_2002), and three Ameriflux sites. Incorporation of GEM into Noah improves the surface energy fluxes as well as the associated diurnal cycle of soil moisture and soil temperature during both wet and dry periods. An analysis of midday, average canopy resistance shows similar day-to-day trends in the model fields as seen in observed patterns. Bias and standard deviation analyses for soil temperature and surface fluxes show that GEM responds somewhat better than the Jarvis scheme, mainly because the Jarvis approach relies on a parametrised minimum canopy resistance and meteorological variables such as air temperature and incident radiation. The analyses suggest that adding a photosynthesis-based transpiration scheme such as GEM improves the ability of the land-data assimilation system to simulate evaporation and transpiration under a range of soil and vegetation conditions.     

Regional climate modelling over complex terrain: an evaluation study of COSMO-CLM hindcast model runs for the Greater Alpine Region

In this study the results of the regional climate model COSMO-CLM (CCLM) covering the Greater Alpine Region (GAR, 4°–19°W and 43°–49°N) were evaluated against observational data. The simulation was carried out as a hindcast run driven by ERA-40 reanalysis data for the period 1961–2000. The spatial resolution of the model data presented is approx. 10 km per grid point. For the evaluation purposes a variety of observational datasets were used: CRU TS 2.1, E-OBS, GPCC4 and HISTALP. Simple statistics such as mean biases, correlations, trends and annual cycles of temperature and precipitation for different sub-regions were applied to verify the model performance. Furthermore, the altitude dependence of these statistical measures has been taken into account. Compared to the CRU and E-OBS datasets CCLM shows an annual mean cold bias of ?0.6 and ?0.7 °C, respectively. Seasonal precipitation sums are generally overestimated by +8 to +23 % depending on the observational dataset with large variations in space and season. Bias and correlation show a dependency on altitude especially in the winter and summer seasons. Temperature trends in CCLM contradict the signals from observations, showing negative trends in summer and autumn which are in contrast to CRU and E-OBS.     

Climate Suitability for Stable Malaria Transmission in Zimbabwe Under Different Climate Change Scenarios

Climate is one factor that determines the potential range of malaria. As such, climate change may work with or against efforts to bring malaria under control. We developed a model of future climate suitability for stable Plasmodium falciparum malaria transmission in Zimbabwe. Current climate suitability for stable malaria transmission was based on the MARA/ARMA model of climatic constraints on the survival and development of the Anopheles vector and the Plasmodium falciparum malaria parasite. We explored potential future geographic distributions of malaria using 16 projections of climate in 2100. The results suggest that, assuming no future human-imposed constraints on malaria transmission, changes in temperature and precipitation could alter the geographic distribution of malaria in Zimbabwe, with previously unsuitable areas of dense human population becoming suitable for transmission. Among all scenarios, the highlands become more suitable for transmission, while the lowveld and areas with low precipitation show varying degrees of change, depending on climate sensitivity and greenhouse gas emission stabilization scenarios, and depending on the general circulation model used. The methods employed can be used within or across other African countries.     

Multi-species chemical data assimilation with the Danish Eulerian hemispheric model: system description and verification

Satellite retrievals of atmospheric composition provide a wealth of data on a global scale. These complement results from atmospheric chemistry-transport models (CTMs), and can be combined using data assimilation. We present two assimilation schemes coupled to the Danish Eulerian Hemispheric Model (DEHM), a three-dimensional, off-line CTM with full photochemistry: a variant on the ensemble Kalman filter and the three-dimensional variational scheme. The aim of this paper is to describe the two schemes and present an initial assessment of their impacts on model skill. Retrievals of multiple atmospheric trace gases are assimilated, namely: NO₂ tropospheric column densities, CH₄ total column densities, and partial column concentrations of O₃, CO and CH₄; these data are retrieved from four satellite sensors. Data for each species are assimilated independently of one another, and other species are only adjusted indirectly via the model’s chemistry and dynamics. Assimilation results are compared with measurements from surface monitoring stations and other satellite retrievals, and preliminary validation results are presented.Reference simulations (without assimilation) grossly underestimate surface CO concentrations, and both assimilation schemes eliminate this large and systematic model bias. The assimilation improves the spatial correlation of modelled CO with surface observations, and improves the spatial correlation between forecasts and retrievals for CO, NO₂ and O₃. Results for CH₄ show a loss of skill due to a mismatch in model bias between two assimilated CH₄ data-sets. Finally, we discuss differences in methodology and results between this paper and a recent study on multi-species chemical data assimilation. Joint optimisation of initial conditions and emission rates offers a promising direction for improving modelled boundary-layer concentrations.     

Middle Paleolithic occupation on a Marine Isotope Stage 5 lakeshore in the Nefud Desert,Saudi Arabia

Major hydrological variations associated with glacial and interglacial climates in North Africa and the Levant have been related to Middle Paleolithic occupations and dispersals, but suitable archaeological sites to explore such relationships are rare on the Arabian Peninsula. Here we report the discovery of Middle Paleolithic assemblages in the Nefud Desert of northern Arabia associated with stratified deposits dated to 75,000 years ago. The site is located in close proximity to a substantial relict lake and indicates that Middle Paleolithic hominins penetrated deeply into the Arabian Peninsula to inhabit landscapes vegetated by grasses and some trees. Our discovery supports the hypothesis of range expansion by Middle Paleolithic populations into Arabia during the final humid phase of Marine Isotope Stage 5, when environmental conditions were still favorable.     

Influence of a coarse interlayer on seawater intrusion and contaminant migration in coastal aquifers

Vertical 2D slice laboratory experiments were carried out in homogenous and layered sand tanks to elucidate the effects of a highly permeable (coarse‐grained sand) interlayer on seawater intrusion and transport of contaminants to a coastal sea. Tidal fluctuations produced oscillations in the seawater–freshwater transition zone, fluctuations of the contaminant infiltration rate and a zigzag contaminant plume outline. The seawater wedge became discontinuous at the (vertical) edges of the interlayer because of increased lateral movement of the seawater–freshwater interface within the interlayer. The contaminant plume formed a tail within the interlayer depending on the tidal stage, and similar to the wedge, its movement was accentuated. A simple analytical model that neglected vertical flow reliably predicted steady‐state seawater intrusion into the coastal aquifer. Numerical modeling was used to gain insight into the groundwater hydrodynamics and contaminant migration. The numerical results confirmed the experimental findings, i.e. that a highly permeable interlayer can provide a rapid transit path for contaminants to reach the seaward boundary and that the interlayer amplifies the effects of tidal fluctuations, resulting in wider transition zones for the seawater wedge and contaminant plume. Numerical simulations further showed that, with increasing interlayer hydraulic conductivity, the maximum seawater intrusion distance inside the interlayer increases approximately linearly. For the fixed‐head contaminant injection condition used, the model showed that contaminant infiltration increases approximately logarithmically with increasing interlayer hydraulic conductivity (other factors held fixed). Copyright © 2013 John Wiley & Sons, Ltd.     

Characterization of dissolved organic matter across the Subtropical Convergence off the South Island,New Zealand

Dissolved organic matter (DOM) was investigated along a gradient across the Subtropical Convergence (STC) off the South Island, New Zealand. Ultrahigh resolution electrospray ionization Fourier transform ion cyclotron resonance mass spectrometry (FT-ICR-MS), excitation emission matrix fluorescence (EEM) spectroscopy, and molecular lignin analysis techniques were used to study this DOM. The analysis revealed a group of compounds found only in the coastal DOM samples, which were also characterized by an elevated terrestrial DOM fluorescence pattern and elevated lignin content if compared to the offshore samples. This group exhibited a high degree of carbon unsaturation, as evident from high double bond equivalence minus oxygen values (DBE-O > 9) and maximum fluorescence intensity. Sulfur-containing molecular formulae for summer DOM samples were much more abundant across the entire transect of the STC compared to winter and exhibited distinctly different hydrogen:carbon and oxygen:carbon elemental ratios, suggesting a possible correlation between biological activity and sulfur compounds. The molecular formulae common to all STC samples were compared with those calculated for DOM extracted from freshwater collected from a stream discharging into Doubtful Sound (West Coast, New Zealand). ESI-FT-ICR-MS measurements undertaken in negative electrospray ionization (ESI?) mode indicated that 30% of the molecular formulae were present in both types of DOM, while in positive ESI mode (ESI+) over 90% of the formulae were present in all samples. Hence, a significant portion of the molecular formulae assigned to the solid-phase extractable DOM pool (SPE-DOM) appear to be identical in oceanic and freshwater samples.     

Diatom-inferred depth models in 8 Canadian boreal lakes: inferred changes in the benthic:planktonic depth boundary and implications for assessment of past droughts

Assessment of past drought in boreal regions, a region predicted to be strongly affected by climate warming, can provide insights into future availability of water. However, limited instrumental data and paleoclimatic data are available for this assessment. To address this lack of data in the boreal region of northwest Ontario, a regional study of lakes in the Winnipeg River Drainage Basin was initiated. Diatom-inferred (D-I) depth models were developed based on surface samples collected along a depth gradient within 8 small boreal lakes. Weighted-averaging and modern analog approaches provided robust within-lake depth models for each of the study lakes, with bootstrapped r² values ranging from 0.90 to 0.98, and root-mean-squared-errors of prediction (RMSEP) between 1.1 and 2.5 m. Large differences in the estimated depth optima for three representative, but common diatom species across our 8 study lakes suggested that within-lake calibration datasets are more appropriate for inferring past drought based in depth models than a regional multi-lake calibration dataset, and that light and related variables are controlling factors governing the maximum depth of benthic taxa. A down-core application of the D-I depth models on a near-shore core from Meekin Lake, retrieved near the present-day ecotone between the benthic and planktonic diatom assemblages indicated highly similar trends in inferred depth (r = 0.96). The models have significant correlations with other metrics of changes in depth including diatom species richness (r = 0.74–0.78) and evenness (r = 0.76–0.8), thereby allowing a check on the strength and direction of the depth inferences down-core. Near-shore cores located near the benthic:planktonic transition is a sensitive region that can provide estimates of past droughts in lakes where such inferences have been difficult to estimate.     

Effect of bed surface roughness on longitudinal dispersion in artificial open channels

In several empirical and modelling studies on river hydraulics, dispersion was negatively correlated to surface roughness. In this study, it was aimed to investigate the influence of surface roughness on longitudinal dispersion under controlled conditions. In artificial flow channels with a length of 104 m, tracer experiments with variations in channel bed material were performed. By use of measured tracer breakthrough curves, average flow velocity, mean longitudinal dispersion, and mean longitudinal dispersivity were calculated. Longitudinal dispersion coefficients ranged from 0·018 m² s^?1 in channels with smooth bed surface up to 0·209 m² s^?1 in channels with coarse gravel as bed material. Longitudinal dispersion was linearly related to mean flow velocity. Accordingly, longitudinal dispersivities ranged between 0·152 ± 0·017 m in channels with smooth bed surface and 0·584 ± 0·015 m in identical channels with a coarse gravel substrate. Grain size and surface roughness of the channel bed were found to correlate positively to longitudinal dispersion. This finding contradicts several existing relations between surface roughness and dispersion. Future studies should include further variation in surface roughness to derive a better‐founded empirical equation forecasting longitudinal dispersion from surface roughness. Copyright © 2011 John Wiley & Sons, Ltd.