Total Net Radiation Research Articles

MONTPEL: A multi-component Penman-Monteith energy balance model

Mechanistic modelling is gradually replacing empiricism in crop models, focusing on leaf-level physiological processes. This shift necessitates simulating crop surface temperature at infra-canopy sub-daily scales but many crop models still rely on empirical formulations for canopy temperature estimation, typically on a daily basis. We developed MONTPEL, a multi-component Penman-Monteith model that allows simulating the crop energy balance with flexible canopy representations (“BigLeaf” vs. “Layered”, “Lumped” vs. “Sunlit-Shaded”) and accounts for atmospheric stability conditions. We analyzed the model behavior, sensitivity and accuracy, using measurements from four wheat (Triticum aestivum L.) experiments conducted under varying pedoclimatic and water stress conditions. Measurements included hourly energy balance terms (total net radiation, soil heat flux, sensible and latent energy fluxes), hourly temperature of the canopy surface or of leaves at different depths inside the canopy, and sunlit and shaded leaf temperatures around solar noon at different dates. MONTPEL reproduced the measured energy balance terms with a root mean square error (RMSE) between 21 and 87 Wm-2 and a coefficient of determination (R²) exceeding 0.65. The model's accuracy in simulating canopy temperature, with RMSE ≤ 2.2 °C and R² ≥ 0.92, remained consistent regardless of measurement scale. Adjusting the aerodynamic resistance for atmospheric stability minimized simulated canopy temperature errors, notably in semi-arid conditions. Crop latent energy flux and temperature were most sensitive to the maximal stomatal conductance (gs,max) parameter. However, using a single gs,max value across the simulated experiments yielded satisfactory results, suggesting a weak sensitivity to the temporal and site-to-site variability of gs,max. Distinguishing sunlit from shaded canopy fractions systematically resulted in lower latent energy fluxes compared to “Lumped” canopy representation results. Analysis identified limitations in the multi-component approach, particularly an unrealistic uniform temperature shift across leaf layers when soil surface temperature changes.

Effect of Drought and Pluvial Climates on the Production and Stability of Different Types of Peanut Cultivars in Guangdong, China

The production and breeding of peanuts was restricted by the frequently extreme climatic conditions in Guangdong province, China. To understand the influence of drought and pluvial climates on peanut traits and yield, a phenotypic investigation of seventy peanut cultivars was conducted from 2018 to 2022; comprehensive field meteorological data collection, and typical drought (2021) and pluvial (2022) climates were recorded. The results revealed that the cultivars achieved the highest single plant pod weight (SPPW) and single plant seed weight (SPSW) of 61.03 g and 45.84 g, respectively, in drought conditions, followed by the control, and finally the pluvial. The SPPW, SPSW and eight agronomy traits exhibited significant differences across the different climatic conditions. Correlation analysis revealed the yield traits and key yield-related traits were positively or negatively correlated with soil water content (SWC), total global radiation (TGR), total precipitation (TP) and total net radiation (TNR). The intermediate and Spanish type cultivars were more stable and productive than the other botanical types of cultivars, commercial varieties exhibited better performance than landraces, and seven cultivars were identified with good production potential, under drought and pluvial conditions. Our study showed that pluvial climate was detrimental to peanut yield, and the SPPW and SPSW were significantly influenced by climates with genotype differences.

Forest impacts on snow accumulation and melt in a semi-arid mountain environment

Snowmelt is complex under heterogeneous forest cover due to spatially variable snow surface energy and mass balances and snow accumulation. Forest canopies influence the under-canopy snowpack net total radiation energy balance by enhancing longwave radiation, shading the surface from shortwave radiation, in addition to intercepting snow, and protecting the snow surface from the wind. Despite the importance of predicting snowmelt timing for water resources, there are limited observations of snowmelt timing in heterogeneous forest cover across the Intermountain West. This research seeks to evaluate the processes that control snowmelt timing and magnitude at two paired forested and open sites in semi-arid southern Idaho, USA. Snow accumulation, snowmelt, and snow energy balance components were measured at a marginal snowpack and seasonal snowpack location in the forest, sparse vegetation, forest edge, and open environments. At both locations, the snow disappeared either later in the forest or relatively uniformly in the open and forest. At the upper elevation location, a later peak in maximum snow depth resulted in more variable snow disappearance timing between the open and forest sites with later snow disappearance in the forest. Snow disappearance timing at the marginal snowpack location was controlled by the magnitude and duration of a late season storm increasing snow depth variability and reducing the shortwave radiation energy input. Here, a shorter duration spring storm resulted in more uniform snowmelt in the forest and open. At both locations, the low-density forests shaded the snow surface into the melt period slowing the melt rate in the forest. However, the forest site had less cold content to overcome before melting started, partially canceling out the forest shading effect. Our results highlight the regional similarities and differences of snow surface energy balance controls on the timing and duration of snowmelt.

Thermal Performance Analysis of Hydronic Ceiling Radiant Panel Heaters Under Different Parameters

Abstract A numerical analysis of a ceiling-type radiant panel heater system was performed to examine the heating performance under different parameters, using the FloEFD code. Three-dimensional models of the room and radiant panel heater were created, and the effects of the Reynolds number, water inlet temperature, pipe diameter and pipe runs on the heating performance of the system were examined in detail. The effects of these parameters on the total heat load, the net radiation rate, and the average surface temperature on the sheet and insulation material have been presented. The total heat load and net radiation rate obtained from the system increase with the increase in the Reynolds number. In addition, a rise in the water inlet temperature increases the heat output of the system. An increase of approximately 500 W was observed in the total heat output as the pipe diameter increased. It was observed, too, that the heat output increased with an increase in pipe runs, although above a certain value the heat output became almost constant. The results of this study could offer information to engineers and manufacturers on the design and use of hydronic radiant systems.

Study of energy balance components at different growth stages of green gram grown in new alluvial zone of West Bengal

An experiment was conducted in the experimental farm of Bidhan Chandra KrishiViswavidyalaya, Nadia, West Bengal to study the radiation pattern and its balance over green gram (Vignaradiata var. Samrat). The BREB method was used to determine the sensible heat flux and latent energy. The net radiation was measured through net radiometer and the ground heat flux was measured using Fourier's law. Both the diurnal and seasonal variation of net radiation were studied. Similarly, the energy balance components were studied regularly for different crop growth stages as well as on diurnal basis. It is observed that the net radiation varies from 6.32 Wm-2 to 606.43 Wm-2. The latent heat flux constitutes more than 50% of the net radiation for all growth stages as depicted by energy balance partitioning. The sensible heat flux is partitioned into 10% to 20% of total net radiation throughout the growth stages of green gram, which is the lowest in magnitude among all three energy fluxes. The relationship between Bowen ratio and Vapour pressure deficit (VPD), Bowen ratio and Canopy air temperature difference (CATD) was studied. It was found that Bowen ratio is negatively correlated with VPD but positively correlated with CATD. This study enables to monitor ET pattern through latent heat flux and microclimatic characteristics through sensible and ground heat flux.

Closure in Surface Flux Estimation by Energy Balance Model: Comparison of Priestly-Taylor and Penman-Monteith Computations for a Tropical Site in Ibadan

Increasing demand to further understand complexity in surface energy flux partitioning necessitates the adaptation of numerous estimation methods to fit the site of observation. This is useful for reducing the uncertainty in physically measurable parameters especially those in tropical regions with high human interference in the atmospheric boundary layer. In this study, we used computations from two methods - the Priestley-Taylor (PT) and the Penman-Monteith (PM), based on the Energy Balance model to ascertain closure performance in the surface flux estimations. The study was carried out at the Nigerian Meteorological Experiment III site (7.38 o N and 3.93 o E, 224.2m) located in Ibadan, Southwest Nigeria. Thirty days of a year (2006) dataset were examined using the Bowen ratio (BR) energy balance model to validate the PT and PM methods. The systems were examined across daily and diurnal cycles to better understand the differences in energy partitioning. Results showed that both systems generally favored latent heat flux compared to sensible heat flux perhaps due to above-normal rainfall during the period. The PM method performed better than the PT method with a period average for the sensible heat and latent heat fluxes as 32.05 Wm -2 and 67.66 Wm -2 respectively, accounting for 29.22% and 61.39% of the total net radiation. The PT method underestimates the sensible heat flux by as much as 19.70 Wm -2 compared to the PM method, with a period average of 12.36 Wm -2 representing 11.26% of total net radiation. The PM method also gives a period average Bowen ratio estimate of 0.55, consistent with the standard range for grasslands. The study suggests that the performance of the PM method is related to its response to heat and water vapor transfer over humid regions and would contribute to further research on land-surface interactions over the tropics. Finally, we propose that the measurement of available energy, net radiation, and ground heat flux should be separated for different collocated systems in order to reduce the forcing of closure and aid in proper partitioning of the fluxes. Keywords: surface energy flux, energy balance model, Priestly-Taylor, Penman-Monteith, West Africa, latent heat, sensible heat, NIMEX_3 DOI: 10.7176/JEES/11-5-05 Publication date: May 31 st 2021

Variability of Surface Radiation Budget Components Over the U.S. From 1996 to 2019—Has Brightening Ceased?

AbstractThe record of downwelling solar irradiance and other surface radiation budget components for the U.S. has been extended through 2019 using SURFRAD Network data. Brightening of surface solar irradiance of +7.36 Wm−2/decade occurred from 1996 through 2012. In 2013, surface solar radiation sharply decreased to the long‐term mean (representing 1996–2019) and remained near that level through 2017. Successive decreases in 2018 and 2019 yielded a dimming trend of −3.90 Wm−2/decade after 2012, but with a high uncertainty owing to the observed variability and brief period covered. Individually, all stations but Penn State showed brightening trends consistent with the network average, and surface solar irradiance decreased at all stations after 2012. Total surface net radiation showed similar tendencies but the reversal from increasing to decreasing was more gradual because of the response of surface net longwave to the changing solar input. Aerosol optical depth decreased continuously throughout the tenure of the network but accounted for only 3% of the variability of surface solar irradiance, while cloud fraction explained 62%. The mean cloud fraction was 2.4% greater during the dimming period than the brightening period but showed no trends due to high interannual variability. However, annual anomalies of direct‐normal solar radiation, which relate to sun duration and clouds, generally increased to 2012 and then decreased thereafter. Collectively, these results indicate that changing cloud cover was the primary source of brightening and dimming over the U.S. from 1996 to 2019.

Estimation of the distribution of the total net radiative flux from satellite and automatic weather station data in the Upper Blue Nile basin, Ethiopia

Net radiation is an important factor in studies of land–atmosphere processes, water resource management, and global climate change. This is particularly true for the Upper Blue Nile (UBN) basin, where significant parts of the basin are dry and evapotranspiration (ET) is a major mechanism for water loss. However, net radiation has not yet been appropriately parameterized in the basin. In this study, we estimated the instantaneous distribution of the net radiation flux in the basin using data from the Moderate Resolution Imaging Spectroradiometer (MODIS) sensor onboard the Terra satellite and Automatic Weather Station (AWS) data. Downward shortwave radiation and air temperature usually vary with topography, so we applied residual kriging spatial interpolation techniques to convert AWS data for point locations into gridded surface data. Simulated net radiation outputs were validated through comparison with independent field measurements. Validation results show that our method successfully reproduced the downward shortwave, upward shortwave, and net radiation fluxes. Using AWS data and residual kriging spatial interpolation techniques makes our results robust and comparable to previous works that used satellite data at a finer spatial resolution than MODIS. The estimated net shortwave, longwave, and total radiation fluxes were in close agreement with ground truth measurements, with mean bias (MB) values of − 14.84, 5.7, and 20.53 W m−2 and root mean square error (RMSE) values 83.43, 32.54, and 78.07 W m−2, respectively. The method presented here has potential applications in research focused on energy balance, ET estimation, and weather prediction for regions with similar physiographic features to those of the Nile basin.

Evaluation of Noah-MP Land-Model Uncertainties over Sparsely Vegetated Sites on the Tibet Plateau

The water budget and energy exchange over the Tibetan Plateau (TP) region play an important role on the Asian monsoon. However, it is not well presented in the current land surface models (LSMs). In this study, uncertainties in the Noah with multiparameterization (Noah-MP) LSM are assessed through physics ensemble simulations in three sparsely vegetated sites located in the central TP. The impact of soil organic matter on energy flux and water cycles, along with the influence of uncertainties in precipitation are explored using observations at those sites during the third Tibetan Plateau Experiment from 1August2014 to31July2015. The greatest uncertainties are in the subprocesses of the canopy resistance, soil moisture limiting factors for evaporation, runoff (RNF) and ground water, and surface-layer parameterization. These uncertain subprocesses do not change across the different precipitation datasets. More precipitation can increase the annual total net radiation (Rn), latent heat flux (LH) and RNF, but decrease sensible heat flux (SH). Soil organic matter enlarges the annual total LH by ~26% but lessens the annual total Rn, SH, and RNF by ~7%, 7%, and 39%, respectively. Its effect on the LH and RNF at the Nagqu site, which has a sand soil texture type, is greater than that at the other two sites with sandy loam. This study highlights the importance of precipitation uncertainties and the effect of soil organic matter on the Noah-MP land-model simulations. It provides a guidance to improve the Noah-MP LSM further and hence the land-atmosphere interactions simulated by weather and climate models over the TP region.

Evaluation of the Solar Radiation in a Seismic Zone

The solar radiation affects the life on Earth and it is important in meteorology, climatology, solar energy, agriculture, hydrology and seismology, too. The Sun warms the earth which reflects a part of the energy. An effect of tectonic stress increases the ground temperature that is radiated into the atmosphere. We study the possibility to use the variation of reflected energy in correlation with seismicity, radon concentration, CO2 emission, ULF radio waves, telluric currents, air ionization, temperature in borehole and acoustic noise. In this case, the solar radiation is a precursor factor. This research is important to study the effects of climate change. Analyzing the solar energy budget we find information about the atmosphere: aerosols, ozone, sulfur dioxide, carbon dioxide, nitrogen dioxide and pollution. Our multidisciplinary network (AeroSolSys) monitors Vrancea, a Romanian area characterized by deep earthquakes. For this study, we use two stations, Plostina and Vrancioaia, where we have a net radiometer and a pyranometer. The first sensor includes two pyranometers (up and low) and two pyrgeometers (up and low), a Pt-100 and a Thermistor. The spectral range covers both the Solar Radiation, 0.3 to 3 micrometers, and the Far Infrared Radiation, 4.5 to 42 micrometers. In this case, we measure and determine the temperature of sensor necessary for pyrgeometers, pyranometer short-wave radiation up and low EyU - EyL, pyrgeometers long-wave radiation up and low EgU - EgL, Albedo, net solar radiation NSR, net far infrared radiation NFIR, net total radiation NTR, sky temperature SkyT and, land surface temperature LST. In Vrancioaia, we have a pyranometer with range 0.3 – 2.8 micrometers and maxim 1600 W/m2. Two video cameras monitor the sky and the measurements are correlated with meteorological equipment installed in each location. We use information about solar activity from NOAA satellites and Copernicus CAMS. The humidity of ground is important in energy budget and it is estimated from meteorological data and telluric currents. The interval of analysis is 3 hours, between 2016.01.01 – 2018.05.31. It is difficult to notice the effect of tectonic stress since the daily and seasonal variations of solar radiation are high but global warming is obvious. Solar radiation monitoring helps us to understand other phenomena like radon emission or air ionization.

Refining surface net radiation estimates in arid and semi-arid climates of Iran

Although the downwelling fluxes exhibit space-time scales of dependency on characteristic of atmospheric variations, especially clouds, the upward fluxes and, hence the net radiation, depends on the variation of surface properties, particularly surface skin temperature and albedo. Evapotranspiration at the land surface depends on the properties of that surface and is determined primarily by the net surface radiation, mostly absorbed solar radiation. Thus, relatively high spatial resolution net radiation data are needed for evapotranspiration studies. Moreover, in more arid environments, the diurnal variations of surface (air and skin) temperature can be large so relatively high (sub-daily) time resolution net radiation is also needed. There are a variety of radiation and surface property products available but they differ in accuracy, space-time resolution and information content. This situation motivated the current study to evaluate multiple sources of information to obtain the best net radiation estimate with the highest space-time resolution from ISCCP FD dataset. This study investigates the accuracy of the ISCCP FD and AIRS surface air and skin temperatures, as well as the ISCCP FD and MODIS surface albedos and aerosol optical depths as the leading source of uncertainty in ISCCP FD dataset. The surface air temperatures, 10-cm soil temperatures and surface solar insolation from a number of surface sites are used to judge the best combinations of data products, especially on clear days. The corresponding surface skin temperatures in ISCCP FD, although they are known to be biased somewhat high, disagreed more with AIRS measurements because of the mismatch of spatial resolutions. The effect of spatial resolution on the comparisons was confirmed using the even higher resolution MODIS surface skin temperature values. The agreement of ISCCP FD surface solar insolation with surface measurements is good (within 2.4–9.1%), but the use of MODIS aerosol optical depths as an alternative was checked and found to not improve the agreement. The MODIS surface albedos differed from the ISCCP FD values by no more than 0.02–0.07, but because these differences are mostly at longer wavelengths, they did not change the net solar radiation very much. Therefore to obtain the best estimate of surface net radiation with the best combination of spatial and temporal resolution, we developed a method to adjust the ISCCP FD surface longwave fluxes using the AIRS surface air and skin temperatures to obtain the higher spatial resolution of the latter (45 km), while retaining the 3-h time intervals of the former. Overall, the refinements reduced the ISCCP FD longwave flux magnitudes by about 25.5–42.1 W/m2 RMS (maximum difference −27.5 W/m2 for incoming longwave radiation and −59 W/m2 for outgoing longwave radiation) with the largest differences occurring at 9:00 and 12:00 UTC near local noon. Combining the ISCCP FD net shortwave radiation data and the AIRS-modified net longwave radiation data changed the total net radiation for summertime by 4.64 to 61.5 W/m2 and for wintertime by 1.06 to 41.88 W/m2 (about 11.1–39.2% of the daily mean).

Quantifying deforestation and forest degradation with thermal response

Deforestation and forest degradation cause the deterioration of resources and ecosystem services. However, there are still no operational indicators to measure forest status, especially for forest degradation. In the present study, we analysed the thermal response number (TRN, calculated by daily total net radiation divided by daily temperature range) of 163 sites including mature forest, disturbed forest, planted forest, shrubland, grassland, savanna vegetation and cropland. TRN generally increased with latitude, however the regression of TRN against latitude differed among vegetation types. Mature forests are superior as thermal buffers, and had significantly higher TRN than disturbed and planted forests. There was a clear boundary between TRN of forest and non-forest vegetation (i.e. grassland and savanna) with the exception of shrubland, whose TRN overlapped with that of forest vegetation. We propose to use the TRN of local mature forest as the optimal TRN (TRNopt). A forest with lower than 75% of TRNopt was identified as subjected to significant disturbance, and forests with 66% of TRNopt was the threshold for deforestation within the absolute latitude from 30° to 55°. Our results emphasized the irreplaceable thermal buffer capacity of mature forest. TRN can be used for early warning of deforestation and degradation risk. It is therefore a valuable tool in the effort to protect forests and prevent deforestation.

Actual and Reference Evapotranspiration in a Cornfield in the Zhangye Oasis, Northwestern China

Evapotranspiration (ET) is an important component of the surface energy balance and water cycle, especially in arid and semiarid regions. The characteristics of the actual evapotranspiration (ETa), which was calculated using the eddy covariance method, and the reference evapotranspiration (ET0), which was estimated using the Food and Agriculture Organisation (FAO) Penman–Monteith method, were analysed. This work focussed on the seasonal variations in evapotranspiration and crop coefficient (Kc) above the heterogeneous canopy of an arid oasis ecosystem in a cornfield of the Zhangye oasis in northwestern China. The results showed that in 2008, the total net radiation (Rn) was 2457.73 MJ∙m−2 and that the rainfall was 117 mm. The average wind velocity, air temperature, and specific humidity, which were observed 2 m above the ground surface, were 1.23 m∙s−1, 7.07 °C, and 3.66 g∙kg−1, respectively. The total ETa and ET0 were 654.69 mm and 1039.92 mm, respectively; thus, the ET0 was higher than the ETa. The difference between the ET0 and ETa was high in summer and autumn, and low in winter and spring. The ETa was greatly influenced by irrigation events, whereas the ET0 was not influenced by irrigation. The ETa and ET0 were both greatly influenced by meteorological elements. The Kc values were less than 0.5 outside of the maize-growing stage and greater than 0.5 during the entire maize-growing stage (from 20 April to 22 September 2008). The Kc values were 0.63, 0.75, 0.78, 0.76, 0.61 and 0.71 at the seedling, shooting, heading, filling, and maturity stages and the entire growth stage, respectively.

Eta model simulations using two radiation schemes in clear-sky conditions

This work evaluates the performance of two radiation parameterization schemes in 30-day clear-sky runs of the Eta model over a region in Southeast Brazil. Two versions of the Eta model are compared: a version using the radiation scheme developed by the Geophysical Fluid Dynamics Laboratory (GFDL) and a recently developed version using the Rapid Radiative Transfer Model for GCM (RRTMG). These simulations are compared against CMSAF satellite data and surface station data. The simulation using RRTMG produced downward surface shortwave radiation fluxes closer to observations and reduced the systematic positive bias of the Eta simulation using the GFDL scheme. The 2-m temperature negative bias found in the Eta-GFDL simulations is reduced in the Eta-RRTMG simulations, which results from a larger net total radiation in the Eta-RRTMG simulations. The new version has better accuracy than the Eta using the GFDL scheme for most of the evaluated variables, particularly for clear-sky conditions.

Spatial variation in annual actual evapotranspiration of terrestrial ecosystems in China: Results from eddy covariance measurements

Understanding the spatial variation in annual actual evapotranspiration (AET) and its influencing factors is crucial for a better understanding of hydrological processes and water resources management. By synthesizing ecosystem-level observations of eddy-covariance flux sites in China (a total of 61 sites), we constructed the most complete AET dataset in China up to now. Based on this dataset, we quantified the statistic characteristics of AET and water budgets (defined as the ratio of AET to annual mean precipitation (MAP), AET/MAP) of terrestrial ecosystems in China. Results showed that AET differed significantly among both different vegetation types and climate types in China, with overall mean AET of 534.7±232.8 mm yr-1. AET/MAP also differed significantly among different climate types, but there were no distinct differences in AET/MAP values across vegetation types, with mean AET/MAP of 0.82±0.28 for non-irrigated ecosystems. We further investigated how the main climatic factors and vegetation attributes control the spatial variation in AET. Our findings revealed that the spatial variation of AET in China was closely correlated with the geographical patterns of climate and vegetation, in which the effects of total annual net radiation (R n), MAP and mean annual air temperature (MAT) were dominant. Thus, we proposed an empirical equation to describe the spatial patterns of AET in China, which could explain about 84% of the spatial variation in AET of terrestrial ecosystems in China. Based on the constructed dataset, we also evaluated the uncertainties of five published global evapotranspiration products in simulating site-specific AET in China. Results showed that large biases in site-specific AET values existed for all five global evapotranspiration products, which indicated that it is necessary to involve more observation data of China in their parameterization or validation, while our AET dataset would provide a data source for it.

Performance of the Goddard multiscale modeling framework with Goddard ice microphysical schemes

AbstractThe multiscale modeling framework (MMF), which replaces traditional cloud parameterizations with cloud‐resolving models (CRMs) within a host atmospheric general circulation model (GCM), has become a new approach for climate modeling. The embedded CRMs make it possible to apply CRM‐based cloud microphysics directly within a GCM. However, most such schemes have never been tested in a global environment for long‐term climate simulation. The benefits of using an MMF to evaluate rigorously and improve microphysics schemes are here demonstrated. Four one‐moment microphysical schemes are implemented into the Goddard MMF and their results validated against three CloudSat/CALIPSO cloud ice products and other satellite data. The new four‐class (cloud ice, snow, graupel, and frozen drops/hail) ice scheme produces a better overall spatial distribution of cloud ice amount, total cloud fractions, net radiation, and total cloud radiative forcing than earlier three‐class ice schemes, with biases within the observational uncertainties. Sensitivity experiments are conducted to examine the impact of recently upgraded microphysical processes on global hydrometeor distributions. Five processes dominate the global distributions of cloud ice and snow amount in long‐term simulations: (1) allowing for ice supersaturation in the saturation adjustment, (2) three additional correction terms in the depositional growth of cloud ice to snow, (3) accounting for cloud ice fall speeds, (4) limiting cloud ice particle size, and (5) new size‐mapping schemes for snow and graupel. Despite the cloud microphysics improvements, systematic errors associated with subgrid processes, cyclic lateral boundaries in the embedded CRMs, and momentum transport remain and will require future improvement.

Factors controlling Nitrous Oxide emission from a spruce forest ecosystem on drained organic soil, derived using the CoupModel

High Nitrous Oxide (N2O) emissions have been identified in hemiboreal forests in association with draining organic soils. However, the specific controlling factors that regulate the emissions remain unclear. To examine the importance of different factors affecting N2O emissions in a spruce forest on drained organic soil, a process-based model, CoupModel, was calibrated using the generalized likelihood uncertainty estimation (GLUE) method. The calibration also aims to estimate parameter density distributions, the covariance matrix of estimated parameters and the correlation between parameters and variables information, useful when applying the model on other peat soil sites and for further model improvements. The calibrated model reproduced most of the high resolution data (total net radiation, soil temperature, groundwater level, net ecosystem exchange, etc.) very well, as well as cumulative measured N2O emissions (simulated 8.7±1.1kgN2Oha−1year−1 (n=97); measured 8.7±2.7kgN2Oha−1year−1 (n=6)), but did not capture every measured peak. Parameter uncertainties were reduced after calibration, in which 16 out of 20 parameters changed from uniform distributions into normal distributions or log normal distributions. Four parameters describing bypass water flow, oxygen diffusion and soil freezing changed significantly after calibration. Inter-connections and correlations between many calibrated parameters and variables reflect the complex and interrelated nature of pedosphere, biosphere and atmosphere interactions. This also highlights the need to calibrate a number of parameters simultaneously. Model sensitivity analysis indicated that N2O emissions during growing seasons are controlled by competition between plants and microbes for nitrogen, while during the winter season snow melt periods are important. Our results also indicate that N2O is mainly produced in the capillary fringe close to the groundwater table by denitrification in the anaerobic zone. We conclude that, in afforested drained peatlands, the plants and groundwater level have important influences on soil N availability, ultimately controlling N2O emissions.

Soil heat flux calculation for sunlit and shaded surfaces under row crops: 1. Model development and sensitivity analysis

Soil heat flux at the surface (G0) is strongly influenced by whether the soil is shaded or sunlit, and therefore can have large spatial variability for incomplete vegetation cover, such as across the interrows of row crops. Most practical soil–plant–atmosphere energy balance models calculate G0 as a function of either total (RN) or soil net radiation (RN,S). Even though RN,S includes sunlit and shaded conditions, this is seldom considered, even at spatial scales of a few m. In order to improve the utility of surface energy balance models designed for row crops at relatively small spatial scales, a method was developed to calculate G0 as a function of shaded, partially sunlit, or fully sunlit RN,S. Calculation of RN,S was derived using a geometric approach, and G0 was derived by the calorimetric method using measurements of soil temperature and volumetric soil water content under upland cotton (Gossypium hirsutum L.) over a wide range of canopy cover conditions. Calorimetric G0 and calculated RN,S were related by assuming their normalized values were equal at 24h time steps and in phase. The method required only a single empirical parameter, which related the 24h minimum G0 (G0,MIN) to the 24h maximum RN,S (RN,S,MAX) as G0,MIN=a×RN,S,MAX, and a=−0.31 was found by simple linear regression (p<0.01). Model sensitivity (SM) of calculated G0 was calculated for sparse, medium, and full canopy cover; nighttime, shaded, partially sunlit, and fully sunlit surface conditions; and north–south and east–west row orientations, where input values of agronomic, shortwave, and longwave input variables were varied ±25% of their base values. The method was most sensitive (1.0<SM<36) to canopy width, canopy height, leaf area index, row spacing, canopy and soil emittances, and canopy and soil temperatures for medium to full canopy cover. Also, there was generally greater sensitivity for shaded and partially sunlit surfaces compared with sunlit and nighttime surfaces, and north–south rows compared with east–west rows. However, the method had little sensitivity (SM usually <0.50) to input variables used to calculate the shortwave components of RN,S.

Changing Trends in Meteorological Elements and Reference Evapotranspiration in a Mega City: A Case Study in Shenzhen City, China

Shenzhen city was a farmland region before 1978, and it then developed to a mega city in China. This type of change in city can greatly affect the climatic conditions. In this study, the daily, monthly, and annual climatic variables and the reference crop evapotranspiration (ET0) for Shenzhen from 1954 to 2012 were computed using the FAO Penman-Monteith equation (PM), and these parameters were analyzed to study the temporal trends of ET0and meteorological factors. The trends and the time points of abrupt changes of ET0and meteorological factors were tested using Mann-Kendall methods. Results show that, in the past 59 years, the annual ET0first decreased from 1954 to 1978, then increased from 1979 to 1990, and now varied slightly after 1990. The mean air temperature rose gradually, and the relative humidity decreased as a whole. These trends finally resulted in an increasing trend in vapor pressure deficit (VPD). The wind speed showed a slightly decreasing trend. Both the annual total sunshine duration and net radiation showed trends of rapid decline. ET0change is sensitive to the hours of sunshine and VPD. The significant increase in ET0after 1979 was mainly due to the increased air temperature and decreased relative humidity.

Analysis of the microclimate and CO2 flux characteristics in arid desert wetland ecosystem in summer: A case study in Xihu desert wetland, Dunhuang, China

Xihu desert wetland is an important and unusual environment in China or even in the world. However, until now, little research has been focused on the microclimate and CO2 flux characteristics in this area. This paper reports the characteristics of daily variations of microclimate and CO2 flux in the Dunhuang Xihu desert wetland, based on data observed in the desert wetland during a period of continuous fine weather in summer 2012. Results indicate that the characteristics of the micrometeorology were significantly affected by the land–lake breeze during the study period, and updrafts were prevalent in this region. The friction wind speed and the vertical velocity were much greater than those in the Maqu grasslands. The energy budget was strongly imbalanced: the latent heat flux was significantly higher than the sensible flux. The daily mean values of total solar radiation and net radiation were larger than those in Maqu grasslands and Jinta oasis. There was a temperature inversion and inverse humidity gradient in the atmospheric surface layer at night. The desert wetland ecosystem was a carbon sink during the whole of the observation period, and the maximum rate of carbon absorption usually occurred at about 11:00 hr each day in this region.