Internal waves,baroclinic energy fluxes and mixing at the European shelf edge

The energy flux in internal waves generated at the Celtic Sea shelf break was estimated by (i) applying perturbation theory to a week-long dataset from a mooring at 200 m depth, and (ii) using a 2D non-hydrostatic circulation model over the shelf break. The dataset consisted of high resolution time-series of currents and vertical stratification together with two 25-h sets of vertical profiles of the dissipation of turbulent kinetic energy. The observations indicated an average energy flux of 139 W m⁻¹, travelling along the shelf break towards the northwest. The average energy flux across the shelf break at the mooring was only 8 W m⁻¹. However, the waves propagating onshelf transported up to 200 W m⁻¹, but they were only present 51% of the time. A comparison between the divergence of the baroclinic energy flux and observed dissipation within the seasonal thermocline at the mooring showed that the dissipation was at least one order of magnitude larger. Results from a 2D model along a transect perpendicular to the shelf break showed a time-averaged onshelf energy flux of 153–425 W m⁻¹, depending on the magnitude of the barotropic forcing. A divergence zone of the energy flux was found a few kilometre offshore of the location of the observations in the model results, and fluxes on the order of several kW m⁻¹ were present in the deep waters further offshelf from the divergence zone. The modelled fluxes exhibited qualitative agreements with the phase and hourly onshelf magnitudes of the observed energy fluxes. Both the observations and the model results show an intermittent onshelf energy flux of 100–200 W m⁻¹, but these waves could only propagate ∼20–30 km onshore before dissipating. This conclusion was supported by a 25-h dataset sampled some 180 km onto the shelf, where a weak wave energy flux was found going towards the shelf break. We therefore conclude that shelf break generated internal waves are unlikely to be the main source of energy for mixing on the inner part of the shelf.     

Utility of an Automated Thermal-Based Approach for Monitoring Evapotranspiration

A very simple remote sensing-based model for water use monitoring is presented. The model acronym DATTUTDUT (Deriving Atmosphere Turbulent Transport Useful To Dummies Using Temperature) is a Dutch word which loosely translates as “it’s unbelievable that it works”. DATTUTDUT is fully automated and only requires a surface temperature map, making it simple to use and providing a rapid estimate of spatially-distributed fluxes. The algorithm is first tested over a range of environmental and land-cover conditions using data from four short-term field experiments and then evaluated over a growing season in an agricultural region. Flux model output is in satisfactory agreement with observations and established remote sensing-based models, except under dry and partial canopy cover conditions. This suggests that DATTUTDUT has utility in identifying relative water use and as an operational tool providing initial estimates of ET anomalies in data-poor regions that would be confirmed using more robust modeling techniques.     

Computing cross-isotherm volume transports from ocean temperature observations and surface heat fluxes,with application to the Barents Sea inflow

A method for determining the cross-isotherm ocean transport from surface heat flux and ocean temperature data is derived. By computing the volume flux through the isotherm that extend from 19°E, 74°N to the eastern part of the Kola Peninsula, the flow through the western entrance of the Barents Sea south of 74°N is estimated. Using three different surface heat flux datasets, the inflow is found to range from 2.9 to 4.5 Sv in winter (October–March) and from 0.4 to 1.4 Sv in summer (April–September; 1 Sv=10⁶ m³ s⁻¹). The seasonal variations are stronger than indicated by results from direct current measurements, probably because the seasonal cycle of the surface heat fluxes is overestimated along the considered isotherm. The annual mean inflow ranges from 1.9 to 2.2 Sv during a cold period (1986–1988), and from 2.4 to 3.0 Sv during a warm period (1990–1992), close to reported observations.     

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.     

Modeling studies of the upper ocean response to a tropical cyclone

A coupled ocean and boundary layer flux numerical modeling system is used to study the upper ocean response to surface heat and momentum fluxes associated with a major hurricane, namely, Hurricane Dennis (July 2005) in the Gulf of Mexico. A suite of experiments is run using this modeling system, constructed by coupling a Navy Coastal Ocean Model simulation of the Gulf of Mexico to an atmospheric flux model. The modeling system is forced by wind fields produced from satellite scatterometer and atmospheric model wind data, and by numerical weather prediction air temperature data. The experiments are initialized from a data assimilative hindcast model run and then forced by surface fluxes with no assimilation for the time during which Hurricane Dennis impacted the region. Four experiments are run to aid in the analysis: one is forced by heat and momentum fluxes, one by only momentum fluxes, one by only heat fluxes, and one with no surface forcing. An equation describing the change in the upper ocean hurricane heat potential due to the storm is developed. Analysis of the model results show that surface heat fluxes are primarily responsible for widespread reduction (0.5°–1.5°C) of sea surface temperature over the inner West Florida Shelf 100–300 km away from the storm center. Momentum fluxes are responsible for stronger surface cooling (2°C) near the center of the storm. The upper ocean heat loss near the storm center of more than 200 MJ/m² is primarily due to the vertical flux of thermal energy between the surface layer and deep ocean. Heat loss to the atmosphere during the storm’s passage is approximately 100–150 MJ/m². The upper ocean cooling is enhanced where the preexisting mixed layer is shallow, e.g., within a cyclonic circulation feature, although the heat flux to the atmosphere in these locations is markedly reduced.     

Early spring turbulent mixing in an ice-covered Arctic fjord during transition to melting

Observations are presented of currents, hydrography and turbulence in a jet-type tidally forced fjord in Svalbard. The fjord was ice covered at the time of the experiment in early spring 2004. Turbulence measurements were conducted by both moored instruments within the uppermost 5 m below the ice and a microstructure profiler covering 3–60 m at 75 m depth. Tidal choking at the mouth of the fjord induces a tidal jet advecting relatively warmer water past the measurement site and dominating the variability in hydrography. While there was no strong correlation with the observed hydrography or mixing and the phase of the semidiurnal tidal cycle, the mean structure in dissipation of turbulent kinetic energy, work done under the ice and the mixing in the water column correlated with the current when conditionally sampled for tidal jet events. Observed levels of dissipation of turbulent kinetic energy per unit mass, 1.1×10⁻⁷ W kg⁻¹, and eddy diffusivity, 7.3×10⁻⁴ m² s⁻¹, were comparable to direct measurements at other coastal sites and shelves with rough topography and strong forcing. During spring tides, an average upward heat flux of 5 W m⁻² in the under-ice boundary layer was observed. Instantaneous (1 h averaged) large heat flux events were correlated with periods of large inflow, hence elevated heat fluxes were associated with the tidal jet and its heat content. Vertical heat fluxes are derived from shear-probe measurements by employing a novel model for eddy diffusivity [Shih et al., 2005. Parameterization of turbulent fluxes and scales using homogeneous sheared stably stratified turbulence simulations. Journal of Fluid Mechanics 525, 193–214]. When compared to the direct heat flux measurements using the eddy correlation method at 5 m below the ice, the upper 4–6 m averaged heat flux estimates from the microstructure profiler agreed with the direct measurements to within 10%. During the experiment water column was stably, but weakly, stratified. Destabilizing buoyancy fluxes recorded close to the ice were absent at 5 m below the ice, and overall, turbulence production was dominated by shear. A scaling for dissipation employing production by both stress and buoyancy [Lombardo and Gregg, 1989. Similarity scaling of viscous and thermal dissipation in a convecting boundary layer. Journal of Geophysical Research 94, 6273–6284] was found to be appropriate for the under-ice boundary layer.     

The NAFE’06 data set: Towards soil moisture retrieval at intermediate resolution

The National Airborne Field Experiment 2006 (NAFE’06) was conducted during a three week period of November 2006 in the Murrumbidgee River catchment, located in southeastern Australia. One objective of NAFE’06 was to explore the suitability of the area for SMOS (Soil Moisture and Ocean Salinity) calibration/validation and develop downscaling and assimilation techniques for when SMOS does come on line. Airborne L-band brightness temperature was mapped at 1 km resolution 11 times (every 1–3 days) over a 40 by 55 km area in the Yanco region and 3 times over a 40 by 50 km area that includes Kyeamba Creek catchment. Moreover, multi-resolution, multi-angle and multi-spectral airborne data including surface temperature, surface reflectance (green, read and near infrared), lidar data and aerial photos were acquired over selected areas to develop downscaling algorithms and test multi-angle and multi-spectral retrieval approaches. The near-surface soil moisture was measured extensively on the ground in eight sampling areas concurrently with aircraft flights, and the soil moisture profile was continuously monitored at 41 sites. Preliminary analyses indicate that (i) the uncertainty of a single ground measurement was typically less than 5% vol. (ii) the spatial variability of ground measurements at 1 km resolution was up to 10% vol. and (iii) the validation of 1 km resolution L-band data is facilitated by selecting pixels with a spatial soil moisture variability lower than the point-scale uncertainty. The sensitivity of passive microwave and thermal data is also compared at 1 km resolution to illustrate the multi-spectral synergy for soil moisture monitoring at improved accuracy and resolution. The data described in this paper are available at www.nafe.unimelb.edu.au.     

Calibration of an integrated land surface process and radiobrightness (LSP/R) model during summertime

In this study, a soil vegetation and atmosphere transfer (SVAT) model was linked with a microwave emission model to simulate microwave signatures for different terrain during summertime, when the energy and moisture fluxes at the land surface are strong. The integrated model, land surface process/radiobrightness (LSP/R), was forced with weather and initial conditions observed during a field experiment. It simulated the fluxes and brightness temperatures for bare soil and brome grass in the Northern Great Plains. The model estimates of soil temperature and moisture profiles and terrain brightness temperatures were compared with the observed values. Overall, the LSP model provides realistic estimates of soil moisture and temperature profiles to be used with a microwave model. The maximum mean differences and standard deviations between the modeled and the observed temperatures (canopy and soil) were 2.6 K and 6.8 K, respectively; those for the volumetric soil moisture were 0.9% and 1.5%, respectively. Brightness temperatures at 19 GHz matched well with the observations for bare soil, when a rough surface model was incorporated indicating reduced dielectric sensitivity to soil moisture by surface roughness. The brightness temperatures of the brome grass matched well with the observations indicating that a simple emission model was sufficient to simulate accurate brightness temperatures for grass typical of that region and surface roughness was not a significant issue for grass-covered soil at 19 GHz. Such integrated SVAT-microwave models allow for direct assimilation of microwave observations and can also be used to understand sensitivity of microwave signatures to changes in weather forcings and soil conditions for different terrain types.     

Hydrodynamics and sediment dynamics in the German Bight. A focus on observations and numerical modelling in the East Frisian Wadden Sea

This work deals with analysis of hydrographic observations and results of numerical simulations. The data base includes acoustic Doppler current profilers (ADCP) observations, continuous measurements on data stations and satellite data originating from the medium resolution imaging spectrometer (MERIS) onboard the European Space Agency (ESA) satellite ENVISAT with a spatial resolution of 300 m. Numerical simulations use nested models with horizontal resolutions ranging from 1 km in the German Bight to 200 m in the East Frisian Wadden Sea coupled with a suspended matter transport model. Modern satellite observations have now a comparable horizontal resolution with high-resolution numerical model of the entire area of the East Frisian Wadden Sea allowing to describe and validate new and so far unknown patterns of sediment distribution. The two data sets are consistent and reveal an oscillatory behaviour of sediment pools to the north of the back-barrier basins and clear propagation patterns of tidally driven suspended particulate matter outflow into the North Sea. The good agreement between observations and simulations is convincing evidence that the model simulates the basic dynamics and sediment transport processes, which motivates its further use in hindcasting, as well as in the initial steps towards forecasting circulation and sediment dynamics in the coastal zone.     

Determination of averaging period parameter and its effects analysis for eddy covariance measurements

It is more and more popular to estimate the exchange of water vapor, heat and CO2fluxes between the land surface and the atmosphere using the eddy covariance technique. To get believable fluxes, it is necessary to correct the observations based on the different surface conditions and to determine relevant techinical parameters. The raw 10 Hz eddy covariance data observed in the Yucheng and Changbai Mountains stations were recalculated by various averaging periods (from 1 to 720 min) respectively, and the recalculated results were compared with the results calculated by the averaging period of 30 mins. Meanwhile, the distinctions of fluxes calculated by different averaging periods were analyzed. The continuous 15 days observations over wheat fields in the Yucheng station were mainly analyzed. The results are shown that: (i) In the Yucheng station, compared with the observations by 30 min, when the averaging period changes from 10 to 60 min, the variations of the eddy-covariance estimates of fluxes were less than 2%; when the averaging period changes less than 10 min, the estimate of fluxes reduced obviously with the reduction of the averaging period (the max relative error was -12%); and when the averaging period exceeds 120 min, the eddy covariance estimates of fluxes will be increased and become unsteady (the max relative error is over 10%); (ii) the eddy covariance estimates of fluxes over wheat field in the Yucheng station suggusted that it is much better to take 10 min as an averaging period in studying diurnal change of fluxes, and take 30min for a long-term flux observation; and (iii) normalized ratio was put forward to determine the range of averaging period of eddy covariance measurements. By comparing the observations over farmlands and those over forests, it is indicated that the increase of eddy covariance estimates over tall forest was more than that over short vegetation when the averaging period increased.     

Effect of remote clouds on surface UV irradiance

Clouds affect local surface UV irradiance, even if the horizontal distance from the radiation observation site amounts to several kilometers. In order to investigate this effect, which we call remote clouds effect, a 3-dimensional radiative transfer model is applied. Assuming the atmosphere is subdivided into a quadratic based sector and its surrounding, we quantify the influence of changing cloud coverage within this surrounding from 0% to 100% on surface UV irradiance at the sector center. To work out this remote clouds influence as a function of sector base size, we made some calculations for different sizes between 10 km × 10 km and 100 km × 100 km. It appears that in the case of small sectors (base size 20 km × 20 km) the remote clouds effect is highly variable: Depending on cloud structure, solar zenith angle and wavelength, the surface UV irradiance may be enhanced up to 15% as well as reduced by more than 50%. In contrast, for larger sectors it is always the case that enhancements become smaller by 5% if sector base size exceeds 60 km × 60 km. However, these values are upper estimates of the remote cloud effects and they are found only for special cloud structures. Since these structures might occur but cannot be regarded as typical, different satellite observed cloud formations (horizontal resolution about 1 km × 1 km) have also been investigated. For these more common cloud distributions we find remote cloud effects to be distinctly smaller than the corresponding upper estimates, e.g., for a sector with base size of 25 km × 25 km the surface UV irradiance error due to ignoring the actual remote clouds and replacing their influence with periodic horizontal boundary conditions is less than 3%, whereas the upper estimate of remote clouds effect would suggest an error close to 10%.     

Simulating the effects of geologic heterogeneity and transient boundary conditions on streambed temperatures — Implications for temperature-based water flux calculations

Analytical solutions to the one-dimensional heat transport equation for steady-state conditions can provide simple means to quantify groundwater surface water exchange. The errors in exchange flux calculations that are introduced when the underlying assumptions of homogeneous sediments and constant temperature boundary conditions are violated were systematically evaluated in a simulation study. Temperatures in heterogeneous sediments were simulated using a numerical model. Heterogeneity in the sediments was represented by discrete, binary geologic units. High contrasts between the hydraulic conductivities (K) of the geologic units were found to lead to large errors, while the influence of the structural arrangement of the units was smaller. The effects of transient temperature boundary conditions were investigated using an analytical equation. Errors introduced by transient boundary conditions were small for Darcy-velocities > 0.1 m d^− 1 in the period near maximum and minimum annual surface water temperatures. For smaller fluxes, however, errors can be large. Assuming steady-state conditions and vertical flow in homogeneous sediments is acceptable at certain times of the year and for medium to high exchange fluxes, but pronounced geologic heterogeneity can lead to large errors.     

Estimation of spatially distributed surface energy fluxes using remotely‐sensed data for agricultural fields

Land surface energy fluxes are required in many environmental studies, including hydrology, agronomy and meteorology. Surface energy balance models simulate microscale energy exchange processes between the ground surface and the atmospheric layer near ground level. Spatial variability of energy fluxes limits point measurements to be used for larger areas. Remote sensing provides the basis for spatial mapping of energy fluxes. Remote‐sensing‐based surface energy flux‐mapping was conducted using seven Landsat images from 1997 to 2002 at four contiguous crop fields located in Polk County, northwestern Minnesota. Spatially distributed surface energy fluxes were estimated and mapped at 30 m pixel level from Landsat Thematic Mapper and Enhanced Thematic Mapper images and weather information. Net radiation was determined using the surface energy balance algorithm for land (SEBAL) procedure. Applying the two‐source energy balance (TSEB) model, the surface temperature and the latent and sensible heat fluxes were partitioned into vegetation and soil components and estimated at the pixel level. Yield data for wheat and soybean from 1997 to 2002 were mapped and compared with latent heat (evapotranspiration) for four of the fields at pixel level. The spatial distribution and the relation of latent heat flux and Bowen ratio (ratio of sensible heat to latent heat) to crop yield were studied. The root‐mean‐square error and the mean absolute percentage of error between the observed and predicted energy fluxes were between 7 and 22 W m⁻² and 12 and 24% respectively. Results show that latent heat flux and Bowen ratio were correlated (positive and negative) to the yield data. Wheat and soybean yields were predicted using latent heat flux with mean R² = 0·67 and 0·70 respectively, average residual means of −4·2 bushels/acre and 0·11 bushels/acre respectively, and average residual standard deviations of 16·2 bushels/acre and 16·6 bushels/acre respectively (1 bushel/acre ≈ 0·087 m³ ha⁻¹). The flux estimation procedure from the SEBAL‐TSEB model was useful and applicable to agricultural fields. Copyright © 2005 John Wiley & Sons, Ltd.     

Simulation of evapotranspiration and carbon dioxide flux in the wheat‐maize rotation croplands of the North China Plain using the Simple Biosphere Model

The North China Plain, which is critical for food production in China, is encountering serious water shortage due to heavy agricultural water requirement. The accurate modelling of carbon dioxide flux and evapotranspiration (ET) in croplands is thus essential for yield prediction and water resources management. The land surface model is powerful in simulating energy and carbon dioxide fluxes between land and atmosphere. Some key processes in the Simple Biosphere Model (Version 2, SiB2) were parameterized based on the observations. The simulated fluxes were tested against the eddy covariance flux measurements over two typical winter wheat/maize double cropping fields. A simple diagnostic parameterisation of soil respiration, not included in SiB2, was added and calibrated using the observations to model the carbon budget. The Ball‐Berry stomatal conductance model was calibrated using observed leaf gas exchange rate, showing that the original SiB2 could significantly underpredict the ET in the wheat field. SiB2 significantly underpredicted soil resistance at the Weishan site, leading to overpredict the ET. Overall, there was a close agreement between the simulated and observed latent heat fluxes and net CO₂ exchange using the re‐parameterized SiB2. These findings are important when the model is used for the regional simulation in the North China Plain. Copyright © 2011 John Wiley & Sons, Ltd.     

Analysis and modeling of spatial correlation structure in small-scale rainfall in Central Oklahoma

Spatial correlation structure in small-scale rainfall is analyzed based on a dense cluster of raingauges in Central Oklahoma. This cluster, called the EVAC PicoNet, consists of 53 gauges installed in 25 measurement stations covering an area of about 3 km by 3 km. Two raingauges are placed in 24 stations and five in the central station. Three aspects of the estimated spatial correlation functions are discussed: dependence on time-scale ranging from 1 min to 24 h, inter-storm variability, and dependence on rainfall intensity. The results show a regular dependence of the correlogram parameters on the averaging time-scale, large differences of the correlograms in the individual storms, and the dominance of storms with high spatial variability on the average large sample characteristics. The authors also demonstrate and discuss the ambiguities in correlation estimates conditioned on rainfall intensities. The findings of this study have implications for raingauge network design, rainfall modeling, and conclusive evaluation of radar and satellite estimates of rainfall.     

Volatile emissions and gas geochemistry of Hot Spring Basin,Yellowstone National Park,USA

We characterize and quantify volatile emissions at Hot Spring Basin (HSB), a large acid-sulfate region that lies just outside the northeastern edge of the 640 ka Yellowstone Caldera. Relative to other thermal areas in Yellowstone, HSB gases are rich in He and H₂, and mildly enriched in CH₄ and H₂S. Gas compositions are consistent with boiling directly off a deep geothermal liquid at depth as it migrates toward the surface. This fluid, and the gases evolved from it, carries geochemical signatures of magmatic volatiles and water–rock reactions with multiple crustal sources, including limestones or quartz-rich sediments with low K/U (or ^40?Ar/^4?He). Variations in gas chemistry across the region reflect reservoir heterogeneity and variable degrees of boiling. Gas-geothermometer temperatures approach 300 °C and suggest that the reservoir feeding HSB is one of the hottest at Yellowstone. Diffuse CO₂ flux in the western basin of HSB, as measured by accumulation-chamber methods, is similar in magnitude to other acid-sulfate areas of Yellowstone and is well correlated to shallow soil temperatures. The extrapolation of diffuse CO₂ fluxes across all the thermal/altered area suggests that 410 ± 140 t d^− 1 CO₂ are emitted at HSB (vent emissions not included). Diffuse fluxes of H₂S were measured in Yellowstone for the first time and likely exceed 2.4 t d^− 1 at HSB. Comparing estimates of the total estimated diffuse H₂S emission to the amount of sulfur as SO₄²⁻ in streams indicates ~ 50% of the original H₂S in the gas emission is lost into shallow groundwater, precipitated as native sulfur, or vented through fumaroles. We estimate the heat output of HSB as ~ 140–370 MW using CO₂ as a tracer for steam condensate, but not including the contribution from fumaroles and hydrothermal vents. Overall, the diffuse heat and volatile fluxes of HSB are as great as some active volcanoes, but they are a small fraction (1–3% for CO₂, 2–8% for heat) of that estimated for the entire Yellowstone system.     

Statistical analysis for thermal data in the Japanese Islands

We calculated statistical average of thermal data to speculate regional thermal structure of the forearc area of the Japanese Islands. The three thermal statistical averages show a difference of a high thermal regime in the western part of forearc inner zone and a low in the Kanto forearc outer zone. The Kanto zone marks 18 K km⁻¹ for mean geothermal gradient, 44 mW m⁻² for mean heat flow, while the western inner zone shows 27 K km⁻¹ for mean geothermal gradient, 63 mW m⁻² for mean heat flow. The geothermal gradients of the Nobi Plain and the Osaka Plain in the western inner zone are 29 and 36 K km⁻¹, respectively, while the value of the Kanto Plain in the Kanto zone is 21 K km⁻¹. Taking account of the effect of accumulation of sediments, we see the difference in the thermal regime between the plains and conclude that the difference is significant. Heat flux in the crust depends on the volume of granite rich in radioactive elements. There are few granitic rocks in the Kanto zone, while granitic rocks are dominant in the western inner zone. The heat flow of 20 mW m⁻² is attributed to the granitic rocks of about 8 km in thickness. There are two oceanic plate subductions of the Pacific plate and the Philippine Sea plate under the Kanto zone, while only the Philippine Sea plate has been subducting under the western inner zone. The model simulation based on thermal and subduction model shows a heat flow ranging 50-60 mW m⁻² in the southwest Japan forarc area and a low value of about 20 mW m⁻² in the northeast Japan forearc area. The heat flux from the cooling oceanic lithosphere depends on the age of plate. The Shikoku Basin, a part of the Philippine Sea plate, off the western inner zone is 15-30 Ma, while the Pacific plate off the Kanto zone is 122-132 Ma. Theoretically, heat flux values of 15 and 50 Ma oceanic plates range 60-120 mW m⁻² and those of 122-132 Ma could be about 10 mW m⁻². If the heat flux contribution from the Philippine Sea plate under the Kanto zone is smaller than the plate under the western inner zone, there could be a thermal regime difference in order of several tens of mW m⁻². Conclusively, the cause of the difference of heat flux could be the uneven granitic rocks distribution and/or the difference of heat flux between the two subducting plate.     

CLM4-LISSS模型对太湖多时间尺度水热通量模拟性能的评估

湖泊模型为数值天气预报模型提供热量通量、水汽通量和动量通量等下边界条件,但是不同时间尺度上湖泊水热通量变化的控制因子不同,因此有必要对湖泊模型进行多时间尺度上的离线评估.本文利用2012-2016年太湖中尺度通量网避风港站的气象资料和辐射数据驱动CLM4-LISSS模型(Community Land Model version 4-Lake,Ice,Snow and Sediment Simulator),并与涡度相关观测(Eddy Covariance,EC)结果进行对比,以年平均潜热通量模拟结果最佳为目标调整了模式中的消光系数、粗糙度长度方案,研究了该模型从半小时到年尺度上对湖表温度和水热通量的模拟性能.结果表明:模型对湖表温度的模拟在各时间尺度上均比较理想,但是模拟的日较差较小;从半小时到年尺度上潜热通量的变化趋势都能被很好地模拟出来,但在季节尺度上,潜热通量的模拟出现了秋冬季偏高、春夏季偏低的情况,季节变化模拟不准确.湖表温度和潜热通量模拟偏差的原因可能是消光系数的参数化方案.相比之下,感热通量尽管年际变化趋势的模拟值与观测值一致,但是从半小时到年尺度均被高估.特别地,冷锋过境期间,模型能较好地模拟出潜热通量和感热通量的变化趋势,但对于高风速条件下的感热通量模拟效果不佳.本文的研究结果能为湖泊模式的应用与发展提供有用信息.     

ParBreZo: A parallel,unstructured grid,Godunov-type,shallow-water code for high-resolution flood inundation modeling at the regional scale

Topographic data are increasingly available at high resolutions (<10 m) over large spatial extents to support detailed flood inundation modeling and loss estimation analyses required for flood risk management. This paper describes ParBreZo, the parallel implementation of a two-dimensional, Godunov-type, shallow-water code, to address the computational demand of high-resolution flood modeling at the regional scale (10²–10⁴ km²). A systematic approach to unstructured grid partitioning (domain decomposition) is presented, and the Single Process Multiple Data (SPMD) paradigm of distributed-memory parallelism is implemented so the code can be executed on computer clusters with distributed memory, shared memory, or some combination of the two (now common with multi-core architectures). In a fully-wetted, load-balanced test problem, the code scales very well with a parallel efficiency of close to 100% on up to 512 processes (maximum tested). A weighted grid partitioning is used to partially address the load balancing challenge posed by partially wetted domains germane to flooding applications, where the flood extent varies over time, while the partitioning remains static. An urban dam-break flood test problem shows that weighted partitions achieve a parallel efficiency exceeding 70% using up to 48 processes. This corresponds to a 97% reduction in execution time so results are obtained in a matter of minutes, which is attractive for routine engineering analyses. A hurricane storm surge test problem shows that a 10 m resolution, 12 h inundation forecast for a 40 km length of coastline can be completed in under 2 h using 512 processors. Hence, if coupled to a hurricane forecast system capable of resolving storm surge, inundation forecasts could be made at 10 m resolution with at least a 10 h lead time.