Evaluation of boundary-layer type in a weather forecast model utilising long-term Doppler lidar observations

Many studies evaluating model boundary-layer schemes focus either on near-surface parameters or on short-term observational campaigns. This reflects the observational datasets that are widely available for use in model evaluation. In this paper we show how surface and long-term Doppler lidar observations, combined in a way to match model representation of the boundary layer as closely as possible, can be used to evaluate the skill of boundary-layer forecasts. We use a 2-year observational dataset from a rural site in the UK to evaluate a climatology of boundary layer type forecast by the UK Met Office Unified Model. In addition, we demonstrate the use of a binary skill score (Symmetric Extremal Dependence Index) to investigate the dependence of forecast skill on season, horizontal resolution and forecast leadtime. A clear diurnal and seasonal cycle can be seen in the climatology of both the model and observations, with the main discrepancies being the model overpredicting cumulus capped and decoupled stratocumulus capped boundary-layers and underpredicting well mixed boundary-layers. Using the SEDI skill score the model is most skillful at predicting the surface stability. The skill of the model in predicting cumulus capped and stratocumulus capped stable boundary layer forecasts is low but greater than a 24 hr persistence forecast. In contrast, the prediction of decoupled boundary-layers and boundary-layers with multiple cloud layers is lower than persistence. This process based evaluation approach has the potential to be applied to other boundary-layer parameterisation schemes with similar decision structures.

A new boundary-based morphological model

Mathematical morphology addresses the problem of describing shapes in an n-dimensional space using the concepts of set theory. A series of standardized morphological operations are defined, and they are applied to the shapes to transform them using another shape called the structuring element. In an industrial environment, the process of manufacturing a piece is based on the manipulation of a primitive object via contact with a tool that transforms the object progressively to obtain the desired design. The analogy with the morphological operation of erosion is obvious. Nevertheless, few references about the relation between the morphological operations and the process of design and manufacturing can be found. The non-deterministic nature of classic mathematical morphology makes it very difficult to adapt their basic operations to the dynamics of concepts such as the ordered trajectory. A new geometric model is presented, inspired by the classic morphological paradigm, which can define objects and apply morphological operations that transform these objects. The model specializes in classic morphological operations, providing them with the determinism inherent in dynamic processes that require an order of application, as is the case for designing and manufacturing objects in professional computer-aided design and manufacturing (CAD/CAM) environments. The operators are boundary-based so that only the points in the frontier are handled. As a consequence, the process is more efficient and more suitable for use in CAD/CAM systems.

Commentary: Using the convection-dispersion model and transit time density functions in the analysis of organ distribution kinetics

The convection-dispersion model and its extended form have been used to describe solute disposition in organs and to predict hepatic availabilities. A range of empirical transit-time density functions has also been used for a similar purpose. The use of the dispersion model with mixed boundary conditions and transit-time density functions has been queried recently by Hisaka and Sugiyanaa in this journal. We suggest that, consistent with soil science and chemical engineering literature, the mixed boundary conditions are appropriate providing concentrations are defined in terms of flux to ensure continuity at the boundaries and mass balance. It is suggested that the use of the inverse Gaussian or other functions as empirical transit-time densities is independent of any boundary condition consideration. The mixed boundary condition solutions of the convection-dispersion model are the easiest to use when linear kinetics applies. In contrast, the closed conditions are easier to apply in a numerical analysis of nonlinear disposition of solutes in organs. We therefore argue that the use of hepatic elimination models should be based on pragmatic considerations, giving emphasis to using the simplest or easiest solution that will give a sufficiently accurate prediction of hepatic pharmacokinetics for a particular application. (C) 2000 Wiley-Liss Inc. and the American Pharmaceutical Association J Pharm Sci 89:1579-1586, 2000.

The simulation of the opposing fluxes of latent heat and CO2 over various land-use types: coupling a gas exchange model to a mesoscale atmospheric model

A mesoscale meteorological model (FOOT3DK) is coupled with a gas exchange model to simulate surface fluxes of CO2 and H2O under field conditions. The gas exchange model consists of a C3 single leaf photosynthesis sub-model and an extended big leaf (sun/shade) sub-model that divides the canopy into sunlit and shaded fractions. Simulated CO2 fluxes of the stand-alone version of the gas exchange model correspond well to eddy-covariance measurements at a test site in a rural area in the west of Germany. The coupled FOOT3DK/gas exchange model is validated for the diurnal cycle at singular grid points, and delivers realistic fluxes with respect to their order of magnitude and to the general daily course. Compared to the Jarvis-based big leaf scheme, simulations of latent heat fluxes with a photosynthesis-based scheme for stomatal conductance are more realistic. As expected, flux averages are strongly influenced by the underlying land cover. While the simulated net ecosystem exchange is highly correlated with leaf area index, this correlation is much weaker for the latent heat flux. Photosynthetic CO2 uptake is associated with transpirational water loss via the stomata, and the resulting opposing surface fluxes of CO2 and H2O are reproduced with the model approach. Over vegetated surfaces it is shown that the coupling of a photosynthesis-based gas exchange model with the land-surface scheme of a mesoscale model results in more realistic simulated latent heat fluxes.

Implementation of non-local boundary layer schemes in the Regional Atmospheric Modeling System and its impact on simulated mesoscale circulations

This paper proposes the implementation of different non-local Planetary Boundary Layer schemes within the Regional Atmospheric Modeling System (RAMS) model. The two selected PBL parameterizations are the Medium-Range Forecast (MRF) PBL and its updated version, known as the Yonsei University (YSU) PBL. YSU is a first-order scheme that uses non-local eddy diffusivity coefficients to compute turbulent fluxes. It is based on the MRF, and improves it with an explicit treatment of the entrainment. With the aim of evaluating the RAMS results for these PBL parameterizations, a series of numerical simulations have been performed and contrasted with the results obtained using the Mellor and Yamada (MY) scheme, also widely used, and the standard PBL scheme in the RAMS model. The numerical study carried out here is focused on mesoscale circulation events during the summer, as these meteorological situations dominate this season of the year in the Western Mediterranean coast. In addition, the sensitivity of these PBL parameterizations to the initial soil moisture content is also evaluated. The results show a warmer and moister PBL for the YSU scheme compared to both MRF and MY. The model presents as well a tendency to overestimate the observed temperature and to underestimate the observed humidity, considering all PBL schemes and a low initial soil moisture content. In addition, the bias between the model and the observations is significantly reduced moistening the initial soil moisture of the corresponding run. Thus, varying this parameter has a positive effect and improves the simulated results in relation to the observations. However, there is still a significant overestimation of the wind speed over flatter terrain, independently of the PBL scheme and the initial soil moisture used, even though a different degree of accuracy is reproduced by RAMS taking into account the different sensitivity tests.

Parametric CAD model based shape optimization using adjoint functions

Adjoint methods have proven to be an efficient way of calculating the gradient of an objective function with respect to a shape parameter for optimisation, with a computational cost nearly independent of the number of the design variables [1]. The approach in this paper links the adjoint surface sensitivities (gradient of objective function with respect to the surface movement) with the parametric design velocities (movement of the surface due to a CAD parameter perturbation) in order to compute the gradient of the objective function with respect to CAD variables.
For a successful implementation of shape optimization strategies in practical industrial cases, the choice of design variables or parameterisation scheme used for the model to be optimized plays a vital role. Where the goal is to base the optimization on a CAD model the choices are to use a NURBS geometry generated from CAD modelling software, where the position of the NURBS control points are the optimisation variables [2] or to use the feature based CAD model with all of the construction history to preserve the design intent [3]. The main advantage of using the feature based model is that the optimized model produced can be directly used for the downstream applications including manufacturing and process planning.
This paper presents an approach for optimization based on the feature based CAD model, which uses CAD parameters defining the features in the model geometry as the design variables. In order to capture the CAD surface movement with respect to the change in design variable, the “Parametric Design Velocity” is calculated, which is defined as the movement of the CAD model boundary in the normal direction due to a change in the parameter value.
The approach presented here for calculating the design velocities represents an advancement in terms of capability and robustness of that described by Robinson et al. [3]. The process can be easily integrated to most industrial optimisation workflows and is immune to the topology and labelling issues highlighted by other CAD based optimisation processes. It considers every continuous (“real value”) parameter type as an optimisation variable, and it can be adapted to work with any CAD modelling software, as long as it has an API which provides access to the values of the parameters which control the model shape and allows the model geometry to be exported. To calculate the movement of the boundary the methodology employs finite differences on the shape of the 3D CAD models before and after the parameter perturbation. The implementation procedure includes calculating the geometrical movement along a normal direction between two discrete representations of the original and perturbed geometry respectively. Parametric design velocities can then be directly linked with adjoint surface sensitivities to extract the gradients to use in a gradient-based optimization algorithm.
The optimisation of a flow optimisation problem is presented, in which the power dissipation of the flow in an automotive air duct is to be reduced by changing the parameters of the CAD geometry created in CATIA V5. The flow sensitivities are computed with the continuous adjoint method for a laminar and turbulent flow [4] and are combined with the parametric design velocities to compute the cost function gradients. A line-search algorithm is then used to update the design variables and proceed further with optimisation process.

Low Level Jet stream of Asian Summer Monsoon and its Variability

The main objective of the of present study are to study the intraseasonal variability of LLJ and its relation with convective heating of the atmosphere, to establish whether LLJ splits into two branches over the Arabian sea as widely believed, the role of horizonatal wind shear of LLJ in the episodes of intense rainfall events observed over the west coast of India, to perform atmospheric modeling work to test whether small (meso) scale vortices form during intense rainfall events along the west coast; and to study the relation between LLJ and monsoon depression genesis. The results of a study on the evolution of Low Level Jetstream (LLJ) prior to the formation of monsoon depressions are presented. A synoptic model of the temporal evolution of monsoon depression has been produced. There is a systematic temporal evolution of the field of deep convection strength and position of the LLJ axis leading to the genesis of monsoon depression. One of the significant outcomes of the present thesis is that the LLJ plays an important role in the intraseasonal and the interannual variability of Indian monsoon activity. Convection and rainfall are dependent mainly on the cyclonic vorticity in the boundary layer associated with LLJ. Monsoon depression genesis and the episodes of very heavy rainfall along the west coast of India are closely related to the cyclonic shear of the LLJ in the boundary layer and the associated deep convection. Case studies by a mesoscale numerical model (MM5) have shown that the heavy rainfall episodes along the west coast of India are associated with generation of mesoscale cyclonic vortices in the boundary layer.

Modeling Of Work-Related Musculoskeletal Disorders In Chemical Process Industry In India

In this modern complex world, stress at work is found to be increasingly a common feature in day to day life. For the same reason, job stress is one of the active areas in occupational health and safety research for over last four decades and is continuing to attract researchers in academia and industry. Job stress in process industries is of concern due to its influence on process safety, and worker‘s safety and health. Safety in process (chemical and nuclear material) industry is of paramount importance, especially in a thickly populated country like India. Stress at job is the main vector in inducing work related musculoskeletal disorders which in turn can affect the worker health and safety in process industries. In view of the above, the process industries should try to minimize the job stress in workers to ensure a safe and healthy working climate for the industry and the worker. This research is mainly aimed at assessing the influence of job stress in inducing work related musculoskeletal disorders in chemical process industries in India

A numerical investigation of the atmosphere-ocean thermal contrast over the coastal upwelling region of Cabo Frio, Brazil

An one-dimensional atmospheric second order closure model, coupled to an oceanic mixed layer model, is used to investigate the short term variation of the atmospheric and oceanic boundary layers in the coastal upwelling area of Cabo Frio, Brazil (23 degrees S, 42 degrees 08`W). The numerical simulations were carried out to evaluate the impact caused by the thermal contrast between atmosphere and ocean on the vertical extent and other properties of both atmospheric and oceanic boundary layers. The numerical simulations were designed taking as reference the observations carried out during the passage of a cold front that disrupted the upwelling regime in Cabo Frio in July of 1992. The simulations indicated that in 10 hours the mechanical mixing, sustained by a constant background flow of 10 in s(-1), increases the atmospheric boundary layer in 214 in when the atmosphere is initially 2 K warmer than the ocean (positive thermal contrast observed during upwelling regime). For an atmosphere initially -2 K colder than the ocean (negative thermal contrast observed during passage of the cold front), the incipient thermal convection intensifies the mechanical mixing increasing the vertical extent of the atmospheric boundary layer in 360 in. The vertical evolution of the atmospheric boundary layer is consistent with the observations carried out in Cabo Frio during upwelling condition. When the upwelling is disrupted, the discrepancy between the simulated and observed atmospheric boundary layer heights in Cabo Frio during July of 1992 increases considerably. During the period of 10 hours, the simulated oceanic mixed layer deepens 2 in and 5.4 in for positive and negative thermal contrasts of 2 K and -2 K, respectively. In the latter case, the larger vertical extent of the oceanic mixed layer is due to the presence of thermal convection in the atmospheric boundary layer, which in turn is associated to the absence of upwelling caused by the passage of cold fronts in Cabo Frio.

Hydrogeologic Setting and Ground-Water Flow Simulations of the Eastern High Plains Regional Study Area, Nebraska

The transport of anthropogenic and natural contaminants to public-supply wells was evaluated in a part of the High Plains aquifer near York, Nebraska, as part of the U.S. Geological Survey National Water-Quality Assessment Program. The aquifer in the Eastern High Plains regional study area is composed of Quaternary alluvial deposits typical of the High Plains aquifer in eastern Nebraska and Kansas, is an important water source for agricultural irrigation and public water supply, and is susceptible and vulnerable to contamination. A six-layer, steady-state ground-water flow model of the High Plains aquifer near York, Nebraska, was constructed and calibrated to average conditions for the time period from 1997 to 2001. The calibrated model and advective particle-tracking simulations were used to compute areas contributing recharge and travel times from recharge areas to selected public-supply wells. Model results indicate recharge from agricultural irrigation return flow and precipitation (about 89 percent of inflow) provides most of the ground-water inflow, whereas the majority of ground-water discharge is to pumping wells (about 78 percent of outflow). Particle-tracking results indicate areas contributing recharge to public-supply wells extend northwest because of the natural ground-water gradient from the northwest to the southeast across the study area. Particle-tracking simulations indicate most ground-water travel times from areas contributing recharge range from 20 to more than 100 years but that some ground water, especially that in the lower confined unit, originates at the upgradient model boundary instead of at the water table in the study area and has travel times of thousands of years.

Impact of the Manaus urban plume on trace gas mixing ratios near the surface in the Amazon Basin: Implications for the NO-NO2-O-3 photostationary state and peroxy radical levels

We measured the mixing ratios of NO, NO2, O-3, and volatile organic carbon as well as the aerosol light-scattering coefficient on a boat platform cruising on rivers downwind of the city of Manaus (Amazonas State, Brazil) in July 2001 (Large-Scale Biosphere-Atmosphere Experiment in Amazonia-Cooperative LBA Airborne Regional Experiment-2001). The dispersion and impact of the Manaus plume was investigated by a combined analysis of ground-based (boat platform) and airborne trace gas and aerosol measurements as well as by meteorological measurements complemented by dispersion calculations (Hybrid Single-Particle Lagrangian Integrated Trajectory model). For the cases with the least anthropogenic influence (including a location in a so far unexplored region similar to 150 km west of Manaus on the Rio Manacapuru), the aerosol scattering coefficient, sigma(s), was below 11 Mm(-1), NOx mixing ratios remained below 0.6 ppb, daytime O-3 mixing ratios were mostly below 20 ppb and maximal isoprene mixing ratios were about 3 ppb in the afternoon. The photostationary state (PSS) was not established for these cases, as indicated by values of the Leighton ratio, Phi, well above unity. Due to the influence of river breeze systems and other thermally driven mesoscale circulations, a change of the synoptic wind direction from east-northeast to south-southeast in the afternoon often caused a substantial increase of ss and trace gas mixing ratios (about threefold for sigma(s), fivefold for NOx, and twofold for O-3), which was associated with the arrival of the Manaus pollution plume at the boat location. The ratio F reached unity within its uncertainty range at NOx mixing ratios of about 3 ppb, indicating "steady-state" conditions in cases when radiation variations, dry deposition, emissions, and reactions mostly involving peroxy radicals (XO2) played a minor role. The median midday/afternoon XO2 mixing ratios estimated using the PSS method range from 90 to 120 parts per trillion (ppt) for the remote cases (sigma(s) < 11 Mm(-1) and NOx < 0.6 ppb), while for the polluted cases our estimates are 15 to 60 ppt. These values are within the range of XO2 estimated by an atmospheric chemistry box model (Chemistry As A Box model Application-Module Efficiently Calculating the Chemistry of the Atmosphere (CAABA/MECCA)-3.0).

Paleoproductivity reconstruction for the eastern equatorial Pacific

Productivity at six core locations in the eastern equatorial Pacific (EEP) was reconstructed with a benthic foraminiferal transfer function. The core records show strong regionality, especially where affected by Peru margin upwelling of deeper Equatorial Undercurrent Water (EUC) (originally coming from the subantarctic). This "Peru margin" record differs from that seen along the equator where divergence leads to shallow upwelling, and it is generally inverse to that seen in cores outside the areas of equatorial upwelling. Principal components analysis shows that the main productivity pattern correlates well to the global oxygen isotope record and has lowest values during isotope stages 2 and 4. In addition to this, equatorial cores show a higher frequency pattern of variation which becomes much more pronounced during MIS 3 and 2. The reconstructions based on benthic foraminifera were tested against those from nonaccumulation rate based inorganic chemical proxies of export production. These were found to correlate well in the region influenced by Peru upwelling, and also to share common features for sites along the equator. All the nonaccumulation rate based paleotracers are consistent with one another and differ from accumulation rate derived proxies. The differences between the two classes of paleotracers may result from uncertainties in calculating actual biogenic fluxes since 230Th-normalized results conform more to those we obtained. Analysis of planktonic carbon isotope values for the EEP, and their comparison to the record of the Pacific subantarctic, indicates that the subantarctic contribution to the EUC was reduced during MIS 3 and 2.

Paleomagnetism at DSDP Sites 86-576 and 86-578

The upper sections of Deep Sea Drilling Project Sites 576 (32°21.4'N, 164°16.5'E) and 578 (33°55.6'N, 151°37.7'E) both have stable detrital remanence that can be correlated with the standard reversal stratigraphy. Site 576 contains all reversals above the base of the Gilbert Epoch (5 m.y.) at about 25 m, whereas Site 578 contains a remarkable section of about 60 reversals extending to Anomaly 5B (15 m.y.) at about 150 m sub-bottom depth. In both cases, the paleomagnetic stratigraphy breaks down when accumulation rates drop below 2 m/m.y. At both sites, authigenic manganiferous clays deposited from 70 to 16 m.y. ago accumulated at about 0.4 m/m.y. Similarly, at both sites, the Pleistocene pulse of eolian debris increased accumulation rates by about 6 m/m.y.**2. From 16 to 2 m.y. ago, however, sediment accumulated at Site 578 about five times as rapidly as at Site 576, apparently because of augmented input to the western site by bottom currents.