Numerical Modelling of the Observed Diurnal Cycle of Indian Monsoon Rainfall

We have analysed the diurnal cycle of rainfall over the Indian region (10S-35N, 60E-100E) using both satellite and in-situ data, and found many interesting features associated with this fundamental, yet under-explored, mode of variability. Since there is a distinct and strong diurnal mode of variability associated with the Indian summer monsoon rainfall, we evaluate the ability of the Weather Research and Forecasting Model (WRF) to simulate the observed diurnal rainfall characteristics. The model (at 54km grid-spacing) is integrated for the month of July, 2006, since this period was particularly favourable for the study of diurnal cycle. We first evaluate the sensitivity of the model to the prescribed sea surface temperature (SST), by using two different SST datasets, namely, Final Analyses (FNL) and Real-time Global (RTG). It was found that with RTG SST the rainfall simulation over central India (CI) was significantly better than that with FNL. On the other hand, over the Bay of Bengal (BoB), rainfall simulated with FNL was marginally better than with RTG. However, the overall performance of RTG SST was found to be better than FNL, and hence it was used for further model simulations. Next, we investigated the role of the convective parameterization scheme on the simulation of diurnal cycle of rainfall. We found that the Kain-Fritsch (KF) scheme performs significantly better than Betts-Miller-Janjić (BMJ) and Grell-Devenyi schemes. We also studied the impact of other physical parameterizations, namely, microphysics, boundary layer, land surface, and the radiation parameterization, on the simulation of diurnal cycle of rainfall, and identified the “best” model configuration. We used this configuration of the “best” model to perform a sensitivity study on the role of various convective components used in the KF scheme. In particular, we studied the role of convective downdrafts, convective timescale, and feedback fraction, on the simulated diurnal cycle of rainfall. The “best” model simulations, in general, show a good agreement with observations. Specifically, (i) Over CI, the simulated diurnal rainfall peak is at 1430 IST, in comparison to the observed 1430-1730 IST peak; (ii) Over Western Ghats and Burmese mountains, the model simulates a diurnal rainfall peak at 1430 IST, as opposed to the observed peak of 1430-1730 IST; (iii) Over Sumatra, both model and observations show a diurnal peak at 1730 IST; (iv) The observed southward propagating diurnal rainfall bands over BoB are weakly simulated by WRF. Besides the diurnal cycle of rainfall, the mean spatial pattern of total rainfall and its partitioning between the convective and stratiform components, are also well simulated. The “best” model configuration was used to conduct two nested simulations with one-way, three-level nesting (54-18-6km) over CI and BoB. While, the 54km and 18km simulations were conducted for the whole of July, 2006, the 6km simulation was carried out for the period 18 - 24 July, 2006. The results of our coarse- and fine-scale numerical simulations of the diurnal cycle of monsoon rainfall will be discussed.

On a Novel Geometric Representation of Rotation

This paper presents a novel method of representing rotation and its application to representing the ranges of motion of coupled joints in the human body, using planar maps. The present work focuses on the viability of this representation for situations that relied on maps on a unit sphere. Maps on a unit sphere have been used in diverse applications such as Gauss map, visibility maps, axis-angle and Euler-angle representations of rotation etc. Computations on a spherical surface are difficult and computationally expensive; all the above applications suffer from problems associated with singularities at the poles. There are methods to represent the ranges of motion of such joints using two-dimensional spherical polygons. The present work proposes to use multiple planar domain “cube” instead of a single spherical domain, to achieve the above objective. The parameterization on the planar domains is easy to obtain and convert to spherical coordinates. Further, there is no localized and extreme distortion of the parameter space and it gives robustness to the computations. The representation has been compared with the spherical representation in terms of computational ease and issues related to singularities. Methods have been proposed to represent joint range of motion and coupled degrees of freedom for various joints in digital human models (such as shoulder, wrist and fingers). A novel method has been proposed to represent twist in addition to the existing swing-swivel representation.

Wave dispersion in a cellular composite with modulated microstructure

In this paper we report a modeling technique and analysis of wave dispersion in a cellular composite laminate with spatially modulated microstructure, which can be modeled by parameterization and homogenization in an appropriate length scale. Higher order beam theory is applied and the system of wave equations are derived. Homogenization of these equations are carried out in the scale of wavelength and frequency of the individual wave modes. Smaller scale scattering below the order of cell size are filtered out in the present approach. The longitudinal dispersion relations for different values of a modulation parameter are analyzed which indicates the existence of stop and pass band patterns. Dispersion relations for flexural-shear case are also analyzed which indicates a tendency toward forming the stop and pass bands for increasing values of a shear stiffness modulation parameter. The effect the phase angle (θ) of the incident wave indicates the existence more number of alternative stop bands and pass bands for θ = 45°.

Topology Optimization of Micromachined Structures with the Manufacturing Constraints of KOH Etching of (110) Silicon

The focus of this paper is on designing useful compliant micro-mechanisms of high-aspect-ratio which can be microfabricated by the cost-effective wet etching of (110) orientation silicon (Si) wafers. Wet etching of (110) Si imposes constraints on the geometry of the realized mechanisms because it allows only etch-through in the form of slots parallel to the wafer's flat with a certain minimum length. In this paper, we incorporate this constraint in the topology optimization and obtain compliant designs that meet the specifications on the desired motion for given input forces. Using this design technique and wet etching, we show that we can realize high-aspect-ratio compliant micro-mechanisms. For a (110) Si wafer of 250 µm thickness, the minimum length of the etch opening to get a slot is found to be 866 µm. The minimum achievable width of the slot is limited by the resolution of the lithography process and this can be a very small value. This is studied by conducting trials with different mask layouts on a (110) Si wafer. These constraints are taken care of by using a suitable design parameterization rather than by imposing the constraints explicitly. Topology optimization, as is well known, gives designs using only the essential design specifications. In this work, we show that our technique also gives manufacturable mechanism designs along with lithography mask layouts. Some designs obtained are transferred to lithography masks and mechanisms are fabricated on (110) Si wafers.

Topology Optimization of Micromachined Structures with Surface Micromachining Manufacturing Constraints

Topology optimization methods have been shown to have extensive application in the design of microsystems. However, their utility in practical situations is restricted to predominantly planar configurations due to the limitations of most microfabrication techniques in realizing structures with arbitrary topologies in the direction perpendicular to the substrate. This study addresses the problem of synthesizing optimal topologies in the out-of-plane direction while obeying the constraints imposed by surface micromachining. A new formulation that achieves this by defining a design space that implicitly obeys the manufacturing constraints with a continuous design parameterization is presented in this paper. This is in contrast to including manufacturing cost in the objective function or constraints. The resulting solutions of the new formulation obtained with gradient-based optimization directly provide the photolithographic mask layouts. Two examples that illustrate the approach for the case of stiff structures are included.

On the role of cloud adjustment time scale in simulating precipitation with Relaxed Arakawa-Schubert convection scheme

The precipitation by Relaxed Arakawa-Schubert cumulus parameterization in a General Circulation Model (GCM) is sensitive to the choice of relaxation parameter or specified cloud adjustment time scale. In the present study, we examine sensitivity of simulated precipitation to the choice of cloud adjustment time scale (tau(adj)) over different parts of the tropics using National Center for Environmental Prediction (NCEP) Seasonal Forecast Model (SFM) during June-September. The results show that a single specified value of tau(adj) performs best only over a particular region and different values are preferred over different parts of the world. To find a relation between tau(adj) and cloud depth (convective activity) we choose six regions over the tropics. Based on the observed relation between outgoing long-wave radiation and tau(adj), we propose a linear cloud-type dependent relaxation parameter to be used in the model. The simulations over most parts of the tropics show improved results due to this newly formulated cloud-type dependent relaxation parameter.

Simulating the daily discharge of the Mandovi River, west coast of India

A hydrological modelling framework was assembled to simulate the daily discharge of the Mandovi River on the Indian west coast. Approximately 90% of the west-coast rainfall, and therefore discharge, occurs during the summer monsoon (June-September), with a peak during July-August. The modelling framework consisted of a digital elevation model (DEM) called GLOBE, a hydrological routing algorithm, the Terrestrial Hydrological Model with Biogeochemistry (THMB), an algorithm to map the rainfall recorded by sparse rain-gauges to the model grid, and a modified Soil Conservation Service Curve Number (SCS-CN) method. A series of discharge simulations (with and without the SCS method) was carried out. The best simulation was obtained after incorporating spatio-temporal variability in the SCS parameters, which was achieved by an objective division of the season into five regimes: the lean season, monsoon onset, peak monsoon, end-monsoon, and post-monsoon. A novel attempt was made to incorporate objectively the different regimes encountered before, during and after the Indian monsoon, into a hydrological modelling framework. The strength of our method lies in the low demand it makes on hydrological data. Apart from information on the average soil type in a region, the entire parameterization is built on the basis of the rainfall that is used to force the model. That the model does not need to be calibrated separately for each river is important, because most of the Indian west-coast basins are ungauged. Hence, even though the model has been validated only for the Mandovi basin, its potential region of application is considerable. In the context of the Prediction in Ungauged Basins (PUB) framework, the potential of the proposed approach is significant, because the discharge of these (ungauged) rivers into the eastern Arabian Sea is not small, making them an important element of the local climate system.

The Treewidth of MDS and Reed-Muller Codes

The constraint complexity of a graphical realization of a linear code is the maximum dimension of the local constraint codes in the realization. The treewidth of a linear code is the least constraint complexity of any of its cycle-free graphical realizations. This notion provides a useful parameterization of the maximum-likelihood decoding complexity for linear codes. In this paper, we show the surprising fact that for maximum distance separable codes and Reed-Muller codes, treewidth equals trelliswidth, which, for a code, is defined to be the least constraint complexity (or branch complexity) of any of its trellis realizations. From this, we obtain exact expressions for the treewidth of these codes, which constitute the only known explicit expressions for the treewidth of algebraic codes.

A study of vertical cloud structure of the Indian summer monsoon using CloudSat data

Precise specification of the vertical distribution of cloud optical properties is important to reduce the uncertainty in quantifying the radiative impacts of clouds. The new global observations of vertical profiles of clouds from the CloudSat mission provide opportunities to describe cloud structures and to improve parameterization of clouds in the weather and climate prediction models. In this study, four years (2007-2010) of observations of vertical structure of clouds from the CloudSat cloud profiling radar have been used to document the mean vertical structure of clouds associated with the Indian summer monsoon (ISM) and its intra-seasonal variability. Active and break monsoon spells associated with the intra-seasonal variability of ISM have been identified by an objective criterion. For the present analysis, we considered CloudSat derived column integrated cloud liquid and ice water, and vertically profiles of cloud liquid and ice water content. Over the South Asian monsoon region, deep convective clouds with large vertical extent (up to 14 km) and large values of cloud water and ice content are observed over the north Bay of Bengal. Deep clouds with large ice water content are also observed over north Arabian Sea and adjoining northwest India, along the west coast of India and the south equatorial Indian Ocean. The active monsoon spells are characterized by enhanced deep convection over the Bay of Bengal, west coast of India and northeast Arabian Sea and suppressed convection over the equatorial Indian Ocean. Over the Bay of Bengal, cloud liquid water content and ice water content is enhanced by similar to 90 and similar to 200 % respectively during the active spells. An interesting feature associated with the active spell is the vertical tilting structure of positive CLWC and CIWC anomalies over the Arabian Sea and the Bay of Bengal, which suggests a pre-conditioning process for the northward propagation of the boreal summer intra-seasonal variability. It is also observed that during the break spells, clouds are not completely suppressed over central India. Instead, clouds with smaller vertical extent (3-5 km) are observed due to the presence of a heat low type of circulation. The present results will be useful for validating the vertical structure of clouds in weather and climate prediction models.

Performance evaluation of chemistry transport models over India

Using continuous and near-real time measurements of the mass concentrations of black carbon (BC) aerosols near the surface, for a period of 1 year (from January to December 2006) from a network of eight observatories spread over different environments of India, a space-time synthesis is generated. The strong seasonal variations observed, with a winter high and summer low, are attributed to the combined effects of changes in synoptic air mass types, modulated strongly by the atmospheric boundary layer dynamics. Spatial distribution shows much higher BC concentration over the Indo-Gangetic Plain (IGP) than the peninsular Indian stations. These were examined against the simulations using two chemical transport models, GOCART (Goddard Global Ozone Chemistry Aerosol Radiation and Transport) and CHIMERE for the first time over Indian region. Both the model simulations significantly deviated from the measurements at all the stations; more so during the winter and pre-monsoon seasons and over mega cities. However, the CHIMERE model simulations show better agreement compared with the measurements. Notwithstanding this, both the models captured the temporal variations; at seasonal and subseasonal timescales and the natural variabilities (intra-seasonal oscillations) fairly well, especially at the off-equatorial stations. It is hypothesized that an improvement in the atmospheric boundary layer (ABL) parameterization scheme for tropical environment might lead to better results with GOCART.

Dynamics of the Indian-Ocean oxygen minimum zones

In the Indian Ocean, mid-depth oxygen minimum zones (OMZs) occur in the Arabian Sea and the Bay of Bengal. The lower part of the Arabian-Sea OMZ (ASOMZ; below 400 m) intensifies northward across the basin; in contrast, its upper part (above 400 m) is located in the central/eastern basin, well east of the most productive regions along the western boundary. The Bay-of-Bengal OMZ (BBOMZ), although strong, is weaker than the ASOMZ. To investigate the processes that maintain the Indian-Ocean OMZs, we obtain a suite of solutions to a coupled biological/physical model. Its physical component is a variable-density, 6 1/2-layer model, in which each layer corresponds to a distinct dynamical regime or water-mass type. Its biological component has six compartments: nutrients, phytoplankton, zooplankton, two size classes of detritus, and oxygen. Because the model grid is non-eddy resolving (0.5 degrees), the biological model also includes a parameterization of enhanced mixing based on the eddy kinetic energy derived from satellite observations. To explore further the impact of local processes on OMZs, we also obtain analytic solutions to a one-dimensional, simplified version of the biological model. Our control run is able to simulate basic features of the oxygen, nutrient, and phytoplankton fields throughout the Indian Ocean. The model OMZs result from a balance, or lack thereof, between a sink of oxygen by remineralization and subsurface oxygen sources due primarily to northward spreading of oxygenated water from the Southern Hemisphere, with a contribution from Persian-Gulf water in the northern Arabian Sea. The northward intensification of the lower ASOMZ results mostly from horizontal mixing since advection is weak in its depth range. The eastward shift of the upper ASOMZ is due primarily to enhanced advection and vertical eddy mixing in the western Arabian Sea, which spread oxygenated waters both horizontally and vertically. Advection carries small detritus from the western boundary into the central/eastern Arabian Sea, where it provides an additional source of remineralization that drives the ASOMZ to suboxic levels. The model BBOMZ is weaker than the ASOMZ because the Bay lacks a remote source of detritus from the western boundary. Although detritus has a prominent annual cycle, the model OMZs do not because there is not enough time for significant remineralization to occur.

Role of the cloud adjustment time scale in simulation of the interannual variability of Indian summer monsoon

The simulation of precipitation in a general circulation model relying on relaxed mass flux cumulus parameterization scheme is sensitive to cloud adjustment time scale (CATS). In this study, the frequency of the dominant intra-seasonal mode and interannual variability of Indian summer monsoon rainfall (ISMR) simulated by an atmospheric general circulation model is shown to be sensitive to the CATS. It has been shown that a longer CATS of about 5 h simulates the spatial distribution of the ISMR better. El Nio Southern Oscillation-ISMR relationship is also sensitive to CATS. The equatorial Indian Ocean rainfall and ISMR coupling is sensitive to CATS. Our study suggests that a careful choice of CATS is necessary for adequate simulation of spatial pattern as well as interannual variation of Indian summer monsoon precipitation.

Freeform Skeletal Shape Optimization of Compliant Mechanisms

Compliant mechanisms are elastic continua used to transmit or transform force and motion mechanically. The topology optimization methods developed for compliant mechanisms also give the shape for a chosen parameterization of the design domain with a fixed mesh. However, in these methods, the shapes of the flexible segments in the resulting optimal solutions are restricted either by the type or the resolution of the design parameterization. This limitation is overcome in this paper by focusing on optimizing the skeletal shape of the compliant segments in a given topology. It is accomplished by identifying such segments in the topology and representing them using Bezier curves. The vertices of the Bezier control polygon are used to parameterize the shape-design space. Uniform parameter steps of the Bezier curves naturally enable adaptive finite element discretization of the segments as their shapes change. Practical constraints such as avoiding intersections with other segments, self-intersections, and restrictions on the available space and material, are incorporated into the formulation. A multi-criteria function from our prior work is used as the objective. Analytical sensitivity analysis for the objective and constraints is presented and is used in the numerical optimization. Examples are included to illustrate the shape optimization method.

Energy Sharing for Multiple Sensor Nodes With Finite Buffers

We consider the problem of finding optimal energy sharing policies that maximize the network performance of a system comprising of multiple sensor nodes and a single energy harvesting (EH) source. Sensor nodes periodically sense the random field and generate data, which is stored in the corresponding data queues. The EH source harnesses energy from ambient energy sources and the generated energy is stored in an energy buffer. Sensor nodes receive energy for data transmission from the EH source. The EH source has to efficiently share the stored energy among the nodes to minimize the long-run average delay in data transmission. We formulate the problem of energy sharing between the nodes in the framework of average cost infinite-horizon Markov decision processes (MDPs). We develop efficient energy sharing algorithms, namely Q-learning algorithm with exploration mechanisms based on the epsilon-greedy method as well as upper confidence bound (UCB). We extend these algorithms by incorporating state and action space aggregation to tackle state-action space explosion in the MDP. We also develop a cross entropy based method that incorporates policy parameterization to find near optimal energy sharing policies. Through simulations, we show that our algorithms yield energy sharing policies that outperform the heuristic greedy method.

Isotopic fractionation in proteins as a measure of hydrogen bond length

If a deuterated molecule containing strong intramolecular hydrogen bonds is placed in a hydrogenated solvent, it may preferentially exchange deuterium for hydrogen. This preference is due to the difference between the vibrational zero-point energy for hydrogen and deuterium. It is found that the associated fractionation factor (I) is correlated with the strength of the intramolecular hydrogen bonds. This correlation has been used to determine the length of the H-bonds (donor-acceptor separation) in a diverse range of enzymes and has been argued to support the existence of short low-barrier H-bonds. Starting with a potential energy surface based on a simple diabatic state model for H-bonds, we calculate (I) as a function of the proton donor-acceptor distance R. For numerical results, we use a parameterization of the model for symmetric 0-H. ``.0 bonds R. H. McKenzie, Chem. Phys. Lett. 535, 196 (2012)]. We consider the relative contributions of the 0-H stretch vibration, O-H bend vibrations (both in plane and out of plane), tunneling splitting effects at finite temperature, and the secondary geometric isotope effect. We compare our total (I) as a function of R with NMR experimental results for enzymes, and in particular with an earlier model parametrization (D(R), used previously to determine bond lengths. (C) 2015 AIP Publishing LLC.