Abstract. Land surface models rarely incorporate the terrestrial phosphorus cycle and its interactions with the carbon cycle, despite the extensive scientific debate about the importance of nitrogen and phosphorus supply for future land carbon uptake. We describe a representation of the terrestrial phosphorus cycle for the ORCHIDEE land surface model, and evaluate it with data from nutrient manipulation experiments along a soil formation chronosequence in Hawaii.ORCHIDEE accounts for the influence of the nutritional state of vegetation on tissue nutrient concentrations, photosynthesis, plant growth, biomass allocation, biochemical (phosphatase-mediated) mineralization, and biological nitrogen fixation. Changes in the nutrient content (quality) of litter affect the carbon use efficiency of decomposition and in return the nutrient availability to vegetation. The model explicitly accounts for root zone depletion of phosphorus as a function of root phosphorus uptake and phosphorus transport from the soil to the root surface.The model captures the observed differences in the foliage stoichiometry of vegetation between an early (300-year) and a late (4.1 Myr) stage of soil development. The contrasting sensitivities of net primary productivity to the addition of either nitrogen, phosphorus, or both among sites are in general reproduced by the model. As observed, the model simulates a preferential stimulation of leaf level productivity when nitrogen stress is alleviated, while leaf level productivity and leaf area index are stimulated equally when phosphorus stress is alleviated. The nutrient use efficiencies in the model are lower than observed primarily due to biases in the nutrient content and turnover of woody biomass.We conclude that ORCHIDEE is able to reproduce the shift from nitrogen to phosphorus limited net primary productivity along the soil development chronosequence, as well as the contrasting responses of net primary productivity to nutrient addition.

No abstract

Phosphorus (P) is an essential element for plant growth. Low P availability in soils is likely to limit crop yields in many parts of the world, but this effect has never been quantified at the global scale by process‐based models. Here we attempt to estimate P limitation in three major cereals worldwide for the year 2000 by combining information on soil P distribution in croplands and a generic crop model, while accounting for the nature of soil‐plant P transport. As a global average, the diffusion‐limited soil P supply meets the crop's P demand corresponding to the climatic yield potential, due to the legacy soil P in highly fertilized areas. However, when focusing on the spatial distribution of P supply versus demand, we found strong limitation in regions like North and South America, Africa, and Eastern Europe. Averaged over grid cells where P supply is lower than demand, the global yield gap due to soil P is estimated at 22, 55, and 26% in winter wheat, maize, and rice. Assuming that a fraction (20%) of the annual P applied in fertilizers is directly available to the plant, the global P yield gap lowers by only 5–10%, underlying the importance of the existing soil P supply in sustaining crop yields. The study offers a base for exploring P limitation in crops worldwide but with certain limitations remaining. These could be better accounted for by describing the agricultural P cycle with a fully coupled and mechanistic soil‐crop model.

No abstract

Abstract. Land surface models rarely incorporate the terrestrial phosphorus cycle and its interactions with the carbon cycle, despite the extensive scientific debate about the importance of nitrogen and phosphorus supply for future land carbon uptake. We describe a representation of the terrestrial phosphorus cycle for the land surface model ORCHIDEE, and evaluate it with data from nutrient manipulation experiments along a soil formation chronosequence in Hawaii. ORCHIDEE accounts for influence of nutritional state of vegetation on tissue nutrient concentrations, photosynthesis, plant growth, biomass allocation, biochemical (phosphatase-mediated) mineralization and biological nitrogen fixation. Changes in nutrient content (quality) of litter affect the carbon use efficiency of decomposition and in return the nutrient availability to vegetation. The model explicitly accounts for root zone depletion of phosphorus as a function of root phosphorus uptake and phosphorus transport from soil to the root surface. The model captures the observed differences in the foliage stoichiometry of vegetation between an early (300yr) and a late stage (4.1 Myr) of soil development. The contrasting sensitivities of net primary productivity to the addition of either nitrogen, phosphorus or both among sites are in general reproduced by the model. As observed, the model simulates a preferential stimulation of leaf level productivity when nitrogen stress is alleviated, while leaf level productivity and leaf area index are stimulated equally when phosphorus stress is alleviated. The nutrient use efficiencies in the model are lower as observed primarily due to biases in the nutrient content and turnover of woody biomass. We conclude that ORCHIDEE is able to reproduce the shift from nitrogen to phosphorus limited net primary productivity along the soil development chronosequence, as well as the contrasting responses of net primary productivity to nutrient addition.

The impact of mesoscale structures on the occurrence of anomalous propagation (AP) conditions for radio waves, including ducts, superrefractive, and subrefractive conditions, was studied. The chosen meteorological situations are the bora wind and the sporadic sea/land breeze (SB/LB) during three selected cases over a large portion of the northern Adriatic. For this purpose, we used available radio soundings and numerical mesoscale model simulations (of real cases and their sensitivity tests) at a horizontal resolution of 1.5 km and 81 vertical levels. The model simulated the occurrences of AP conditions satisfactorily, although their intensities and frequency were underestimated at times. Certain difficulties appeared in reproducing the vertical profile of the modified refractive index, which is mainly dependent on the accuracy of the modeled humidity. The spatial distributions of summer AP conditions reveal that the surface layer above the sea (roughly between 30 and 100 m asl) is often covered by superrefractive conditions and ducts. The SB is highly associated with the formations of AP conditions: (i) in the first 100 m asl, where trapping and superrefractive conditions form because of the advection of cold and moist air, and (ii) inside the transition layer between the SB body and the elevated return flow in the form of subrefractive conditions. When deep convection occurs, all three types of AP conditions are caused by the downdraft beneath the cumulonimbus cloud base in its mature phase that creates smaller but marked pools of cold and dry air. The bora wind usually creates a pattern of AP conditions associated with the hydraulic jump and influences distribution of AP conditions over the sea surface.

No abstract

Abstract. The northeastern (NE) Adriatic in the northernMediterranean is the area with (i) the highest frequency of thunderstorms in Croatia, and (ii) frequent appearances of sea breeze (SB) along the coast. This study investigates the impact of the combined large-scale wind (associated with particular synoptic conditions) and the SB on the moist convection development over the NE Adriatic. The four selected cases were (i) chosen on the basis of a daytime moist convection; (ii) supplemented by one of the dominant large-scale winds with seaward (NE, NW) and landward (SW, SE) directions and (iii) simulated by WRF numerical model.The near-surface wind patterns consisted of SBs along the coastline, generated a narrow eastward-moving convergence zone (CZ) along the area if the large-scale wind was less than 9 m s −1 (below 500 hPa). Apart from the low-level CZ, the advection of large-scale wind influenced the lifetime and movement of the initial Cb cells. While the local front collision with the NE wind advection caused the thunderstorm to propagate southward, the CZ and fronts interaction determined the afternoon northwestward storm movement against the NW large-scale wind. Due to particular synoptic background, the thunderstorm event in SE case was the shortest with only a minor impact on the SB. While the origins and locations of storm cells were completely controlled by the low-level CZ and the upward advection of low-level moisture at the SB front, the most typical convective case with SW warm-wet wind only partially supported the SB-Cb interaction.