Net radiation dynamics: performance of 20 daily net radiation models as related to model structure and intricacy in two climates.

We compared daily net radiation (Rn) estimates from 19 methods with the ASCE-EWRI Rn estimates in two climates: Clay Center, Nebraska (sub-humid) and Davis, California (semi-arid) for the calendar year. The performances of all 20 methods, including the ASCE-EWRI Rn method, were then evaluated against Rn data measured over a non-stressed maize canopy during two growing seasons in 2005 and 2006 at Clay Center. Methods differ in terms of inputs, structure, and equation intricacy. Most methods differ in estimating the cloudiness factor, emissivity (e), and calculating net longwave radiation (Rnl). All methods use albedo (a) of 0.23 for a reference grass/alfalfa surface. When comparing the performance of all 20 Rn methods with measured Rn, we hypothesized that the a values for grass/alfalfa and non-stressed maize canopy were similar enough to only cause minor differences in Rn and grass- and alfalfa-reference evapotranspiration (ETo and ETr) estimates. The measured seasonal average a for the maize canopy was 0.19 in both years. Using a = 0.19 instead of a = 0.23 resulted in 6% overestimation of Rn. Using a = 0.19 instead of a = 0.23 for ETo and ETr estimations, the 6% difference in Rn translated to only 4% and 3% differences in ETo and ETr, respectively, supporting the validity of our hypothesis. Most methods had good correlations with the ASCE-EWRI Rn (r2 > 0.95). The root mean square difference (RMSD) was less than 2 MJ m-2 d-1 between 12 methods and the ASCE-EWRI Rn at Clay Center and between 14 methods and the ASCE-EWRI Rn at Davis. The performance of some methods showed variations between the two climates. In general, r2 values were higher for the semi-arid climate than for the sub-humid climate. Methods that use dynamic e as a function of mean air temperature performed better in both climates than those that calculate e using actual vapor pressure. The ASCE-EWRI-estimated Rn values had one of the best agreements with the measured Rn (r2 = 0.93, RMSD = 1.44 MJ m-2 d-1), and estimates were within 7% of the measured Rn. The Rn estimates from six methods, including the ASCE-EWRI, were not significantly different from measured Rn. Most methods underestimated measured Rn by 6% to 23%. Some of the differences between measured and estimated Rn were attributed to the poor estimation of Rnl. We conducted sensitivity analyses to evaluate the effect of Rnl on Rn, ETo, and ETr. The Rnl effect on Rn was linear and strong, but its effect on ETo and ETr was subsidiary. Results suggest that the Rn data measured over green vegetation (e.g., irrigated maize canopy) can be an alternative Rn data source for ET estimations when measured Rn data over the reference surface are not available. In the absence of measured Rn, another alternative would be using one of the Rn models that we analyzed when all the input variables are not available to solve the ASCE-EWRI Rn equation. Our results can be used to provide practical information on which method to select based on data availability for reliable estimates of daily Rn in climates similar to Clay Center and Davis.

Optimal Nonlinear Control and Estimation for a Reusable Launch Vehicle During Reentry Phase

In this paper a nonlinear optimal controller has been designed for aerodynamic control during the reentry phase of the Reusable Launch Vehicle (RLV). The controller has been designed based on a recently developed technique Optimal Dynamic Inversion (ODI). For full state feedback the controller has required full information about the system states. In this work an Extended Kalman filter (EKF) is developed to estimate the states. The vehicle (RLV) has been has been consider as a nonlinear Six-Degree-Of-Freedom (6-DOF) model. The simulation results shows that EKF gives a very good estimation of the states and it is working well with ODI. The resultant trajectories are very similar to those obtained by perfect state feedback using ODI only.

Construction and use of linear regression models for processor performance analysis

Processor architects have a challenging task of evaluating a large design space consisting of several interacting parameters and optimizations. In order to assist architects in making crucial design decisions, we build linear regression models that relate Processor performance to micro-architecture parameters, using simulation based experiments. We obtain good approximate models using an iterative process in which Akaike's information criteria is used to extract a good linear model from a small set of simulations, and limited further simulation is guided by the model using D-optimal experimental designs. The iterative process is repeated until desired error bounds are achieved. We used this procedure to establish the relationship of the CPI performance response to 26 key micro-architectural parameters using a detailed cycle-by-cycle superscalar processor simulator The resulting models provide a significance ordering on all micro-architectural parameters and their interactions, and explain the performance variations of micro-architectural techniques.

Pipelined ring algorithm for matrix multiplication on transputer networks: performance analysis and estimation

A parallel matrix multiplication algorithm is presented, and studies of its performance and estimation are discussed. The algorithm is implemented on a network of transputers connected in a ring topology. An efficient scheme for partitioning the input matrices is introduced which enables overlapping computation with communication. This makes the algorithm achieve near-ideal speed-up for reasonably large matrices. Analytical expressions for the execution time of the algorithm have been derived by analysing its computation and communication characteristics. These expressions are validated by comparing the theoretical results of the performance with the experimental values obtained on a four-transputer network for both square and irregular matrices. The analytical model is also used to estimate the performance of the algorithm for a varying number of transputers and varying problem sizes. Although the algorithm is implemented on transputers, the methodology and the partitioning scheme presented in this paper are quite general and can be implemented on other processors which have the capability of overlapping computation with communication. The equations for performance prediction can also be extended to other multiprocessor systems.

Identification and estimation of size from the beaks of 18 species of cephalopods from the Pacific Ocean

A method of identifying the beaks and estimating body weight and mantle length of 18 species of cephalopods from the Pacific Ocean is presented. Twenty specimens were selected from each of the following cephalopod species: Symplectoteuthis oualaniensis, Dosidicus gigas, Ommastrephes bartramii, S. luminosa, Todarodes pacificus, Nototodarus hawaiiensis, Ornithoteuthis volalilis, Hyaloteuthis pelagica, Onychoteuthis banksii, Pterygioteuthis giardi, Abraliopsis affinis, A. felis, Liocranchia reinhardti, Leachia danae, Histioteuthis heteropsis, H. dofleini, Gonalus onyx, and Loligo opalescens. Dimensions measured on the upper and lower beak are converted to ratios and compared individually among the species using an analysis of variance procedure with Tukey's omega and Duncan's multiple range tests. Significant differences (P =0.05) observed among the species' beak ratio means and structural characteristics are used to construct artificial keys for the upper and lower beaks of the 18 species. Upper and lower beak dimensions are used as independent variables in a linear regression model with mantle length and body weight (log transformed). (PDF file contains 56 pages.)

New directions in sparse sampling and estimation for underdetermined systems

A central objective in signal processing is to infer meaningful information from a set of measurements or data. While most signal models have an overdetermined structure (the number of unknowns less than the number of equations), traditionally very few statistical estimation problems have considered a data model which is underdetermined (number of unknowns more than the number of equations). However, in recent times, an explosion of theoretical and computational methods have been developed primarily to study underdetermined systems by imposing sparsity on the unknown variables. This is motivated by the observation that inspite of the huge volume of data that arises in sensor networks, genomics, imaging, particle physics, web search etc., their information content is often much smaller compared to the number of raw measurements. This has given rise to the possibility of reducing the number of measurements by down sampling the data, which automatically gives rise to underdetermined systems.

In this thesis, we provide new directions for estimation in an underdetermined system, both for a class of parameter estimation problems and also for the problem of sparse recovery in compressive sensing. There are two main contributions of the thesis: design of new sampling and statistical estimation algorithms for array processing, and development of improved guarantees for sparse reconstruction by introducing a statistical framework to the recovery problem.

We consider underdetermined observation models in array processing where the number of unknown sources simultaneously received by the array can be considerably larger than the number of physical sensors. We study new sparse spatial sampling schemes (array geometries) as well as propose new recovery algorithms that can exploit priors on the unknown signals and unambiguously identify all the sources. The proposed sampling structure is generic enough to be extended to multiple dimensions as well as to exploit different kinds of priors in the model such as correlation, higher order moments, etc.

Recognizing the role of correlation priors and suitable sampling schemes for underdetermined estimation in array processing, we introduce a correlation aware framework for recovering sparse support in compressive sensing. We show that it is possible to strictly increase the size of the recoverable sparse support using this framework provided the measurement matrix is suitably designed. The proposed nested and coprime arrays are shown to be appropriate candidates in this regard. We also provide new guarantees for convex and greedy formulations of the support recovery problem and demonstrate that it is possible to strictly improve upon existing guarantees.

This new paradigm of underdetermined estimation that explicitly establishes the fundamental interplay between sampling, statistical priors and the underlying sparsity, leads to exciting future research directions in a variety of application areas, and also gives rise to new questions that can lead to stand-alone theoretical results in their own right.

Beyond the blur: Construction and characterization of the first autonomous AO system and An AO survey of magnetar proper motions

Adaptive optics (AO) corrects distortions created by atmospheric turbulence and delivers diffraction-limited images on ground-based telescopes. The vastly improved spatial resolution and sensitivity has been utilized for studying everything from the magnetic fields of sunspots upto the internal dynamics of high-redshift galaxies. This thesis about AO science from small and large telescopes is divided into two parts: Robo-AO and magnetar kinematics.

In the first part, I discuss the construction and performance of the world’s first fully autonomous visible light AO system, Robo-AO, at the Palomar 60-inch telescope. Robo-AO operates extremely efficiently with an overhead < 50s, typically observing about 22 targets every hour. We have performed large AO programs observing a total of over 7,500 targets since May 2012. In the visible band, the images have a Strehl ratio of about 10% and achieve a contrast of upto 6 magnitudes at a separation of 1′′. The full-width at half maximum achieved is 110–130 milli-arcsecond. I describe how Robo-AO is used to constrain the evolutionary models of low-mass pre-main-sequence stars by measuring resolved spectral energy distributions of stellar multiples in the visible band, more than doubling the current sample. I conclude this part with a discussion of possible future improvements to the Robo-AO system.

In the second part, I describe a study of magnetar kinematics using high-resolution near-infrared (NIR) AO imaging from the 10-meter Keck II telescope. Measuring the proper motions of five magnetars with a precision of upto 0.7 milli-arcsecond/yr, we have more than tripled the previously known sample of magnetar proper motions and proved that magnetar kinematics are equivalent to those of radio pulsars. We conclusively showed that SGR 1900+14 and SGR 1806-20 were ejected from the stellar clusters with which they were traditionally associated. The inferred kinematic ages of these two magnetars are 6±1.8 kyr and 650±300 yr respectively. These ages are a factor of three to four times greater than their respective characteristic ages. The calculated braking index is close to unity as compared to three for the vacuum dipole model and 2.5-2.8 as measured for young pulsars. I conclude this section by describing a search for NIR counterparts of new magnetars and a future promise of polarimetric investigation of a magnetars’ NIR emission mechanism.

Articulated wooden blades for a horizontal-axis wind turbine. Design, construction and testing.

A pair of blades were constructed following a Tapered Chord, Zero Twist pattern after Anderson. The construction uses the Wood Epoxy Saturation Technique, with a solid Beech main spar and leading edge joined together with laminated veneers of beech forming a D-section; the trailing edge is formed from millimetre ply skins, foam filled to resist compressive loads. This construction leads to an extremely light, flexible blade, with the centres of gravity and torsion well forward, giving good stability. Each blade has three built-in strain gauges, alowing flapwise bending to be measured. Stiffness, and natural frequencies, were measured, to input to a numerical computer model to calculate blade deformation during operation, and to determine stability boundaries of the blade. Preliminary aerodynamic performance measurements are presented and close agreement is found with theory.

Use of contexts in language model interpolation and adaptation

Language models (LMs) are often constructed by building multiple individual component models that are combined using context independent interpolation weights. By tuning these weights, using either perplexity or discriminative approaches, it is possible to adapt LMs to a particular task. This paper investigates the use of context dependent weighting in both interpolation and test-time adaptation of language models. Depending on the previous word contexts, a discrete history weighting function is used to adjust the contribution from each component model. As this dramatically increases the number of parameters to estimate, robust weight estimation schemes are required. Several approaches are described in this paper. The first approach is based on MAP estimation where interpolation weights of lower order contexts are used as smoothing priors. The second approach uses training data to ensure robust estimation of LM interpolation weights. This can also serve as a smoothing prior for MAP adaptation. A normalized perplexity metric is proposed to handle the bias of the standard perplexity criterion to corpus size. A range of schemes to combine weight information obtained from training data and test data hypotheses are also proposed to improve robustness during context dependent LM adaptation. In addition, a minimum Bayes' risk (MBR) based discriminative training scheme is also proposed. An efficient weighted finite state transducer (WFST) decoding algorithm for context dependent interpolation is also presented. The proposed technique was evaluated using a state-of-the-art Mandarin Chinese broadcast speech transcription task. Character error rate (CER) reductions up to 7.3 relative were obtained as well as consistent perplexity improvements. © 2012 Elsevier Ltd. All rights reserved.

Automatic construction and natural-language description of nonparametric regression models

Copyright © 2014, Association for the Advancement of Artificial Intelligence (www.aaai.org). All rights reserved. This paper presents the beginnings of an automatic statistician, focusing on regression problems. Our system explores an open-ended space of statistical models to discover a good explanation of a data set, and then produces a detailed report with figures and natural- language text. Our approach treats unknown regression functions non- parametrically using Gaussian processes, which has two important consequences. First, Gaussian processes can model functions in terms of high-level properties (e.g. smoothness, trends, periodicity, changepoints). Taken together with the compositional structure of our language of models this allows us to automatically describe functions in simple terms. Second, the use of flexible nonparametric models and a rich language for composing them in an open-ended manner also results in state- of-the-art extrapolation performance evaluated over 13 real time series data sets from various domains.

Construction and evaluation of DNA vaccines encoding Edwardsiella tarda antigens

Edwardsiella tarda is an opportunistic pathogen that can infect humans, animal, and fish. Two E. tarda antigens, Eta6 and FliC, which are homologues to an ecotin precursor and the FliC flagellin, respectively, were identified by in vivo-induced antigen technology from a pathogenic E. tarda strain isolated from diseased fish. When used as a subunit vaccine, purified recombinant Eta6 was moderately protective against lethal challenge of E. tarda in a Japanese flounder model, whereas purified recombinant FliC showed no apparent immunciprotectivity. Similarly, DNA vaccines based on eta6 and fliC in the form of plasmids pEta6 and pFliC induced, respectively, moderate and marginal protection against E. tarda infection. To improve the vaccine efficacy of eta6, a chimeric DNA vaccine, pCE6, was constructed, which encodes Eta6 fused in-frame to FliC. pCE6 was found to induce significantly higher level of protection than pEta6. Likewise, another chimeric DNA vaccine, pCE18, which expresses FliC fused to a previously identified E. tarda antigen Et18, elicited significantly stronger protective immunity than the DNA vaccine based on et18 alone. Fish immunized with pEta6 and pCE6 produced specific serum antibodies and exhibited significantly enhanced expression of the genes encoding elements that are involved in both innate and adaptive immune responses. Furthermore, the induction magnitudes of most of these genes were significantly higher in pCE6-vaccinated fish than in pEta6-vaccinated fish. (C) 2009 Elsevier Ltd. All rights reserved.

Construction and analysis of an experimental Streptococcus iniae DNA vaccine

Streptococcus iniae is a severe aquaculture pathogen that can also infect humans and animal. A putative secretory antigen, Slat 0, was identified from a pathogenic S. iniae strain by in vivo-induced antigen technology. Using turbot as an animal model, the immunoprotective effect of Sia10 was examined as a DNA vaccine in the form of plasmid pSia10, which expresses sia10 under the cytomegalovirus immediate-early promoter. In fish vaccinated with pSia10, transcription of sia10 was detected in muscle, liver, spleen, and kidney at 7, 14, 21, 28, 35, 42, and 49 days post-vaccination. In addition, production of Sia10 protein was also detected in the muscle tissues of pSia10-vaccinated fish. Fish vaccinated with pSia10 exhibited a relative percent survival (RPS) of 73.9% and 92.3%, respectively, when challenged with high and low doses (producing a cumulative mortality of 92% and 52%, respectively, in the control groups) of S. iniae. Immunological and transcriptional analyses showed that vaccination with pSia10(i) induced much stronger chemiluminescence response and significantly higher levels of nitric oxide production and acid phosphatase activity in head kidney macrophages; (ii) caused the production of specific serum antibodies, which afforded apparent immunoprotection when transferred passively into naive fish; and (iii) upregulated the expression of the genes encoding proteins that are possibly involved in both innate and adaptive immune responses. Taken together, these results indicated that pSia10 is an effective vaccine candidate and may be used in the control of S. iniae infection in aquaculture. (C) 2010 Elsevier Ltd. All rights reserved.

Analysis of the sea-level variability along the Chinese coast and estimation of the impact of a CO2-perturbed atmospheric circulation

The relationship between monthly sea-level data measured at stations located along the Chinese coast and concurrent large-scale atmospheric forcing in the period 1960-1990 is examined. It is found that sea-level varies quite coherently along the whole coast, despite the geographical extension of the station set. A canonical correlation analysis between sea-level and sea-level pressure (SLP) indicates that a great part of the sea-level variability can be explained by the action of the wind stress on the ocean surface. The relationship between sea-level and sea-level pressure is analyzed separately for the summer and winter half-years. In winter, one factor affecting sea-level variability at all stations is the SLP contrast between the continent and the Pacific Ocean, hence the intensity of the winter Monsoon circulation. Another factor that affects coherently all stations is the intensity of the zonal circulation at mid-latitudes. In the summer half year, on the other hand, the influence of SLP on sea-level is spatially less coherent: the stations in the Yellow Sea are affected by a more localized circulation anomaly pattern, whereas the rest of the stations is more directly connected to the intensity of the zonal circulation. Based on this analysis, statistical models (different for summer and winter) to hindcast coastal sealevel anomalies from the large-scale SLP field are formulated. These models have been tested by fitting their internal parameters in a test period and reproducing reasonably the sea-level evolution in an independent period. These statistical models are also used to estimate the contribution of the changes of the atmospheric circulation on sea-level along the Chinese coast in an altered climate. For this purpose the ouput of 150 year-long experiment with the coupled ocean-atmosphere model ECHAM1-LSG has been analyzed, in which the atmospheric concentration of greenhouse gases was continuously increased from 1940 until 2090, according to the Scenario A projection of the Intergovermental Panel on Climate Change. In this experiment the meridional (zonal) circulation relevant for sea-level tends to become weaker (stronger) in the winter half year and stronger (weaker) in summer. The estimated contribution of this atmospheric circulation changes to coastal sea-level is of the order of a few centimeters at the end of the integration, being in winter negative in the Yellow Sea and positive in the China Sea with opposite signs in the summer half-year.

A von Bertalanffy Based Model for the Estimation of Oyster (Crassostrea virginica) Growth on Restored Oyster Reefs in Chesapeake Bay

A model to estimate the mean monthly growth of Crassostrea virginica oysters in Chesapeake Bay was developed. This model is based on the classic von Bertalanffy growth function, however the growth constant is changed every monthly timestep in response to short term changes in temperature and salinity. Using a dynamically varying growth constant allows the model to capture seasonal oscillations in growth, and growth responses to changing environmental conditions that previous applications of the von Bertalanffy model do not capture. This model is further expanded to include an estimation of Perkinsus marinus impacts on growth rates as well as estimations of ecosystem services provided by a restored oyster bar over time. The model was validated by comparing growth estimates from the model to oyster shell height observations from a variety of restoration sites in the upper Chesapeake Bay. Without using the P. marinus impact on growth, the model consistently overestimates mean oyster growth. However, when P. marinus effects are included in the model, the model estimates match the observed mean shell height closely for at least the first 3 years of growth. The estimates of ecosystem services suggested by this model imply that even with high levels of mortality on an oyster reef, the ecosystem services provided by that reef can still be maintained by growth for several years. Because larger oyster filter more water than smaller ones, larger oysters contribute more to the filtration and nutrient removal ecosystem services of the reef. Therefore a reef with an abundance of larger oysters will provide better filtration and nutrient removal. This implies that if an oyster restoration project is trying to improve water quality through oyster filtration, it is important to maintain the larger older oysters on the reef.