Methodology for construction, calibration and validation of a national hydrological model for Denmark

An integrated groundwater/surface water hydrological model with a 1 km² grid has been constructed for Denmark covering 43,000 km². The model is composed of a relatively simple root zone component for estimating the net precipitation, a comprehensive three-dimensional groundwater component for estimating recharge to and hydraulic heads in different geological layers, and a river component for streamflow routing and calculating stream–aquifer interaction. The model was constructed on the basis of the MIKE SHE code and by utilising comprehensive national databases on geology, soil, topography, river systems, climate and hydrology. The present paper describes the modelling process for the 7330 km² island of Sjælland with emphasis on the problems experienced in combining the classical paradigms of groundwater modelling, such as inverse modelling of steady-state conditions, and catchment modelling, focussing on dynamic conditions and discharge simulation. Three model versions with different assumptions on input data and parameter values were required until the performance of the final, according to pre-defined accuracy criteria, model was evaluated as being satisfactory. The paper highlights the methodological issues related to establishment of performance criteria, parameterisation and assessment of parameter values from field data, calibration and validation test schemes. Most of the parameter values were assessed directly from field data, while about 10 ‘free’ parameters were subject to calibration using a combination of inverse steady-state groundwater modelling and manual trial-and-error dynamic groundwater/surface water modelling. Emphasising the importance of tests against independent data, the validation schemes included combinations of split-sample tests (another period) and proxy-basin tests (another area).     

Scaling effects on modeled surface energy-balance components using the NOAH-OSU land surface model

As surface exchange processes are highly non-linear and heterogeneous in space and time, it is important to know the appropriate scale for the reasonable prediction of these exchange processes. For example, the explicit representation of surface variability has been vital in predicting mesoscale weather events such as late-afternoon thunderstorms initiated by latent heat exchanges in mid-latitude regions of the continental United States. This study was undertaken to examine the effects of different spatial scales of input data on modeled fluxes, so as to better understand the resolution needed for accurate modeling. A statistical procedure was followed to select two cells from the Southern Great Plains 1997 hydrology experiment region, each 20 km×20 km, representing the most homogeneous and the most heterogeneous surface conditions (based on soil and vegetation) within the study region. The NOAH-OSU (Oregon State University) Land Surface Model (LSM) was employed to estimate surface energy fluxes. Three scales of study (200 m, 2 and 20 km) were considered in order to investigate the impacts of the aggregation of input data, especially soil and vegetation inputs, on the model output. Model results of net radiation and latent, sensible and ground heat fluxes were compared for the three scales. For the heterogeneous area, the model output at the 20-km resolution showed some differences when compared with the 200-m and 2-km resolutions. This was more pronounced in latent heat (12% decrease), sensible heat (22% increase), and ground heat flux (44% increase) estimation than in net radiation. The scaling effects were much less for the relatively homogeneous land area with 5% increase in sensible heat and 4% decrease in ground heat flux estimation. All of the model outputs for the 2- and 20-km resolutions were in close agreement. The results suggested that, for this study region, soils and vegetation input resolution of about 2 km should be chosen for realistic modeling of surface exchange processes. This resolution was sufficient to capture the effects of sub-grid scale heterogeneity, while avoiding the data and computational difficulties associated with higher spatial resolutions.     

On the circulation in the East Frisian Wadden Sea: numerical modeling and data analysis

 In this paper we use a combination of numerical modeling and data analysis to gain a better understanding of the major characteristics of the circulation in the East Frisian Wadden Sea. In particular, we concentrate on the asymmetry of the tidal wave and its modulation in the coastal area, which results in a complex pattern of responses to the sea-level forcing from the North Sea. The numerical simulations are based on the 3-D primitive equation General Estuarine Transport Model (GETM) with a horizontal resolution of 200 m and terrain-following vertical coordinates. The model is forced at its open boundaries with sea-level data from an operational model for the German Bight (German Hydrographic Office). The validation data for our model simulations include time series of tidal gauge data and surface currents measured at a pile in the back-barrier basin of the Island Langeoog, as well as several ADCP transects in the Accumer Ee tidal inlet. Circulation and turbulence characteristics are investigated for typical situations driven by spring and neap tides, and the analysis is focused on dominating temporal and spatial patterns. By investigating the response of five back-barrier basins with rather different morphologies to external forcing, an attempt is made to elucidate the dominating physical balances controlling the circulation in the individual sub-basins. It is demonstrated that the friction at the seabed tends to slow down the tidal signal in the shallow water. This leads to the establishment of flood dominance in the shallow sea north of the barrier islands. South of the islands, where the water volume of the channels at low tide is smaller than the tidal prism, the asymmetry of the tidal signal is shifted towards ebb dominance, a feature which is particularly pronounced at spring tide. At the northern open boundary, the tidal wave propagating from west to east generates a sea-level difference of ∼1 m along the boundary, and thereby triggers vigorous alongshore currents. The frictional control in the model is located in the inlets, as well as along the northern boundary. The correlation between velocity and turbulent kinetic energy tends to the establishment of a net southward transport, giving theoretical support to the observed accumulation of sediments on the intertidal flats. Weak turbulence along the northern shores of the barrier islands and the small magnitude of the residual currents there promote accumulation of suspended matter in these areas, although wave action will generally counteract this effect. Received: 29 May 2002 / Accepted: 26 September 2002 Responsible Editor: Jean-Marie Beckers Acknowledgements We are indebted to S. Dick for providing the data from the operational model of BSH and to B. Flemming for the useful discussions. The topography data and Fig. 1 have been prepared in cooperation with F. Meyer. Figure 2 has been prepared by G. Brink-Spalink. We also thank for the comments from an anonymous reviewer which helped to improve our paper.     

观测有效性的度量方法

基于Hilbert空间理论 ,通过对观测方程的系数阵实施格拉姆_施密特 (Gram_Schmidt,G_S)正交化过程 ,分析观测的两个特性 ,从一个新的角度论证了参数估计的必要观测问题 ,以及多余观测与必要观测的关系 ;用欧氏范数‖·‖2 和Frobenius‖·‖ F 范数作为基本度量工具对观测提供参数估计的信息量进行度量分析 ,研究观测的有效性和参数估计的可靠性。采用观测空间分析方法研究了测绘学科中必要观测、多余观测等基本概念 ;提出了观测信息量、观测有效性的概念和相应的度量方法 ;从观测的空间结构出发分析参数估计系统的可靠性 ,从一个侧面扩展和丰富了测绘学中的可靠性理论 ;丰富了研究测绘学中一些基本问题的分析方法 ,也为信息学中的数据挖掘问题提供了一种分析思路     

一种新的多通道GPS共视资料的处理算法

提出了一种按仰角加权的多通道GPS共视资料处理算法，充分利用了所有观测数据，同时抑制了多径等因素影响，从而可提高共视比对精度。对日本邮政省通讯研究所(CRL)和国家授时中心(NTSC)2002年5月～2003年3月的多通道共视比对的处理结果表明，新的处理算法结果更接近真实钟差。     

The climatic impacts of land surface change and carbon management, and the implications for climate-change mitigation policy

Strategies to mitigate anthropogenic climate change recognize that carbon sequestration in the terrestrial biosphere can reduce the build-up of carbon dioxide in the Earth’s atmosphere. However, climate mitigation policies do not generally incorporate the effects of these changes in the land surface on the surface albedo, the fluxes of sensible and latent heat to the atmosphere, and the distribution of energy within the climate system. Changes in these components of the surface energy budget can affect the local, regional, and global climate. Given the goal of mitigating climate change, it is important to consider all of the effects of changes in terrestrial vegetation and to work toward a better understanding of the full climate system. Acknowledging the importance of land surface change as a component of climate change makes it more challenging to create a system of credits and debits wherein emission or sequestration of carbon in the biosphere is equated with emission of carbon from fossil fuels. Recognition of the complexity of human-caused changes in climate does not, however, weaken the importance of actions that would seek to minimize our disturbance of the Earth’s environmental system and that would reduce societal and ecological vulnerability to environmental change and variability.     

Comparison of the Langrangian Footprint

We compare flux and concentration footprint estimates of athree-dimensional Lagrangian stochastic dispersion modelapplying backward trajectories with the results of ananalytical footprint model by Kormann and Meixner.The comparison is performed for varying stability regimesof the surface layer as well as for different measurementheights. In general, excellent correspondence is found.     

全球观测系统实地观测的内容和要求

全球性观测系统在检测，监测和预测地球系统的变化中发挥着越工重要的作用。地球系统的复杂性，多尺度性，非线性性，突变性，非平衡性等要求各种全球性观测系统之间必须加强联系与协调。评述了全球性观测系统联系与协调的发展趋势，并依据《全球性观测系统实地观测的内容和体系纲要》、介绍了全球性观测系统实地观测的内容和规范要求。     

煤层倾角与覆岩变形破裂分带

地下采煤中三带(冒落带、裂隙带、弯曲带)高度的判断是预测采空引起地表变形破坏程度的依据。在诸多的影响因素中, 煤层倾角是控制三带高度的最主要因素。建立在实测统计基础上的经验范围值用于三带高度的判断带有一定的人为性。急倾斜煤层中冒落带和裂隙带高度随煤层倾角变化的规律已被褐示[7].通过弹塑性岩石材料的非线性有限元模拟, 本文提出了利用应力重分布图判断中、缓倾角煤层采空区覆岩三带高度的方法, 并应用于实际工程。     

Comparison of Aircraft and Ground-Based FluxMeasurements during OASIS95

Aircraft and ground-based measurements made during the1995 Australian OASIS field campaign are compared. The aircraft data were recorded during low-level flightsat 6 m above ground level and grid flights at altitudes of between 15 and 65 m, allin unstable atmospheric conditions. The low-level flights revealed an inadequate temperaturesensor response time, a correction for which was determined from subsequent work ina wind tunnel. Aircraft and ground-based measurements of mean wind speed, wind directionand air temperature agree to within 0.2 m s^-1, 4° and 0.9 °C respectively.Comparisons between aircraft and ground-based observations of the standarddeviations of vertical velocity, horizontal wind speed, air temperature and specifichumidity have slopes of 0.96, 0.97, 0.92 and 0.99 respectively but the observed scatter isroughly twice the random error expected due to the averaging length of the aircraft data andthe averaging period of the ground-based data. For the low-level flights, the ground-basedand aircraft measurements of sensible and latent heat flux show mean differences of 27 and-25 W m^-2 respectively, which reduce to 11 and -4 W m^-2 respectivelywhen analysis of aircraft data is limited to areas immediately adjacent to the fluxtowers. For the flights at 15 to 65 m above ground level, the mean differences between theground-based and aircraft measurements of sensible and latent heat flux are -22 and-1 W m^-2 respectively and these change to -1 and -7 W m^-2 respectively oncethe effect of surface heterogeneity is included. Aircraft and ground-based measurementsof net radiation agree to within 6% at one ground-based site but differ by 20% at a second.Aircraft measurements of friction velocity at 6 m above the ground agree well withground-based data, but those from flights between 15 and 65 m above ground level do not.This is because at these heights the aircraft measurements provide the local shear stress,not the surface shear stress. Overall, the level of agreement allows confidence in the aircraftdata provided due care is taken of instrument response times and differences in thesurfaces sampled by aircraft and ground-based systems.