Spatio-Temporal Scanning and Statistical Test of the Accelerating Moment Release (AMR) Model Using Australian Earthquake Data

The Accelerating Moment Release (AMR) preceding earthquakes with magnitude above 5 in Australia that occurred during the last 20 years was analyzed to test the Critical Point Hypothesis. Twelve earthquakes in the catalog were chosen based on a criterion for the number of nearby events. Results show that seven sequences with numerous events recorded leading up to the main earthquake exhibited accelerating moment release. Two occurred near in time and space to other earthquakes preceded by AM R. The remaining three sequences had very few events in the catalog so the lack of AMR detected in the analysis may be related to catalog incompleteness. Spatio-temporal scanning of AMR parameters shows that 80% of the areas in which AMR occurred experienced large events. In areas of similar background seismicity with no large events, 10 out of 12 cases exhibit no AMR, and two others are false alarms where AMR was observed but no large event followed. The relationship between AMR and Load-Unload Response Ratio (LURR) was studied. Both methods predict similar critical region sizes, however, the critical point time using AMR is slightly earlier than the time of the critical point LURR anomaly.

Scientific Data Application Profile Scoping Study

This study was undertaken by UKOLN on behalf of the Joint Information Systems Committee (JISC) in the period April to September 2008. Application profiles are metadata schemata which consist of data elements drawn from one or more namespaces, optimized for a particular local application. They offer a way for particular communities to base the interoperability specifications they create and use for their digital material on established open standards. This offers the potential for digital materials to be accessed, used and curated effectively both within and beyond the communities in which they were created. The JISC recognized the need to undertake a scoping study to investigate metadata application profile requirements for scientific data in relation to digital repositories, and specifically concerning descriptive metadata to support resource discovery and other functions such as preservation. This followed on from the development of the Scholarly Works Application Profile (SWAP) undertaken within the JISC Digital Repositories Programme and led by Andy Powell (Eduserv Foundation) and Julie Allinson (RRT UKOLN) on behalf of the JISC. Aims and Objectives 1.To assess whether a single metadata AP for research data, or a small number thereof, would improve resource discovery or discovery-to-delivery in any useful or significant way. 2.If so, then to:a.assess whether the development of such AP(s) is practical and if so, how much effort it would take; b.scope a community uptake strategy that is likely to be successful, identifying the main barriers and key stakeholders. 3.Otherwise, to investigate how best to improve cross-discipline, cross-community discovery-to-delivery for research data, and make recommendations to the JISC and others as appropriate. Approach The Study used a broad conception of what constitutes scientific data, namely data gathered, collated, structured and analysed using a recognizably scientific method, with a bias towards quantitative methods. The approach taken was to map out the landscape of existing data centres, repositories and associated projects, and conduct a survey of the discovery-to-delivery metadata they use or have defined, alongside any insights they have gained from working with this metadata. This was followed up by a series of unstructured interviews, discussing use cases for a Scientific Data Application Profile, and how widely a single profile might be applied. On the latter point, matters of granularity, the experimental/measurement contrast, the quantitative/qualitative contrast, the raw/derived data contrast, and the homogeneous/heterogeneous data collection contrast were discussed. The Study report was loosely structured according to the Singapore Framework for Dublin Core Application Profiles, and in turn considered: the possible use cases for a Scientific Data Application Profile; existing domain models that could either be used or adapted for use within such a profile; and a comparison existing metadata profiles and standards to identify candidate elements for inclusion in the description set profile for scientific data. The report also considered how the application profile might be implemented, its relationship to other application profiles, the alternatives to constructing a Scientific Data Application Profile, the development effort required, and what could be done to encourage uptake in the community. The conclusions of the Study were validated through a reference group of stakeholders.

KE Workshop report: Report Knowledge Exchange workshop Making data count - Research Data and Research Assessment

The possibilities of digital research have altered the production, publication and use of research results. Academic research practice and culture are changing or have already been transformed, but to a large degree the system of academic recognition has not yet adapted to the practices and possibilities of digital research. This applies especially to research data, which are increasingly produced, managed, published and archived, but play hardly a role yet in practices of research assessment. The aim of the workshop was to bring experts and stakeholders from research institutions, universities, scholarly societies and funding agencies together in order to review, discuss and build on possibilities to implement the culture of sharing and to integrate publication of data into research assessment procedures. The report 'The Value of Research Data - Metrics for datasets from a cultural and technical point of view' was presented and discussed. Some of the key finding were that data sharing should be considered normal research practice, in fact not sharing should be considered malpractice. Research funders and universities should support and encourage data sharing. There are a number of important aspects to consider when making data count in research and evaluation procedures. Metrics are a necessary tool in monitoring the sharing of data sets. However, data metrics are at present not very well developed and there is not yet enough experience in what these metrics actually mean. It is important to implement the culture of sharing through codes of conducts in the scientific communities. For further key findings please read the report.

Seafood safety: economics of Hazard Analysis and Critical Control Point (HACCP) programmes

Seventy percent of the world's catch of fish and fishery products is consumed as food. Fish and shellfish products represent 15.6 percent of animal protein supply and 5.6 percent of total protein supply on a worldwide basis. Developing countries account for almost 50 percent of global fish exports. Seafood-borne disease or illness outbreaks affect consumers both physically and financially, and create regulatory problems for both importing and exporting countries. Seafood safety as a commodity cannot be purchased in the marketplace and government intervenes to regulate the safety and quality of seafood. Theoretical issues and data limitations create problems in estimating what consumers will pay for seafood safety and quality. The costs and benefits of seafood safety must be considered at all levels, including the fishers, fish farmers, input suppliers to fishing, processing and trade, seafood processors, seafood distributors, consumers and government. Hazard Analysis Critical Control Point (HACCP) programmes are being implemented on a worldwide basis for seafood. Studies have been completed to estimate the cost of HACCP in various shrimp, fish and shellfish plants in the United States, and are underway for some seafood plants in the United Kingdom, Canada and Africa. Major developments within the last two decades have created a set of complex trading situations for seafood. Current events indicate that seafood safety and quality can be used as non-tariff barriers to free trade. Research priorities necessary to estimate the economic value and impacts of achieving safer seafood are outlined at the consumer, seafood production and processing, trade and government levels. An extensive list of references on the economics of seafood safety and quality is presented. (PDF contains 56 pages; captured from html.)

Turbulent collision of inertial particles: point-particle based, hybrid simulations and beyond

Point-particle based direct numerical simulation (PPDNS) has been a productive research tool for studying both single-particle and particle-pair statistics of inertial particles suspended in a turbulent carrier flow. Here we focus on its use in addressing particle-pair statistics relevant to the quantification of turbulent collision rate of inertial particles. PPDNS is particularly useful as the interaction of particles with small-scale (dissipative) turbulent motion of the carrier flow is mostly relevant. Furthermore, since the particle size may be much smaller than the Kolmogorov length of the background fluid turbulence, a large number of particles are needed to accumulate meaningful pair statistics. Starting from the relative simple Lagrangian tracking of so-called ghost particles, PPDNS has significantly advanced our theoretical understanding of the kinematic formulation of the turbulent geometric collision kernel by providing essential data on dynamic collision kernel, radial relative velocity, and radial distribution function. A recent extension of PPDNS is a hybrid direct numerical simulation (HDNS) approach in which the effect of local hydrodynamic interactions of particles is considered, allowing quantitative assessment of the enhancement of collision efficiency by fluid turbulence. Limitations and open issues in PPDNS and HDNS are discussed. Finally, on-going studies of turbulent collision of inertial particles using large-eddy simulations and particle- resolved simulations are briefly discussed.

Methods of multiparameter inversion of seismic data using the acoustic and elastic born approximations

This thesis presents two different forms of the Born approximations for acoustic and elastic wavefields and discusses their application to the inversion of seismic data. The Born approximation is valid for small amplitude heterogeneities superimposed over a slowly varying background. The first method is related to frequency-wavenumber migration methods. It is shown to properly recover two independent acoustic parameters within the bandpass of the source time function of the experiment for contrasts of about 5 percent from data generated using an exact theory for flat interfaces. The independent determination of two parameters is shown to depend on the angle coverage of the medium. For surface data, the impedance profile is well recovered.

The second method explored is mathematically similar to iterative tomographic methods recently introduced in the geophysical literature. Its basis is an integral relation between the scattered wavefield and the medium parameters obtained after applying a far-field approximation to the first-order Born approximation. The Davidon-Fletcher-Powell algorithm is used since it converges faster than the steepest descent method. It consists essentially of successive backprojections of the recorded wavefield, with angular and propagation weighing coefficients for density and bulk modulus. After each backprojection, the forward problem is computed and the residual evaluated. Each backprojection is similar to a before-stack Kirchhoff migration and is therefore readily applicable to seismic data. Several examples of reconstruction for simple point scatterer models are performed. Recovery of the amplitudes of the anomalies are improved with successive iterations. Iterations also improve the sharpness of the images.

The elastic Born approximation, with the addition of a far-field approximation is shown to correspond physically to a sum of WKBJ-asymptotic scattered rays. Four types of scattered rays enter in the sum, corresponding to P-P, P-S, S-P and S-S pairs of incident-scattered rays. Incident rays propagate in the background medium, interacting only once with the scatterers. Scattered rays propagate as if in the background medium, with no interaction with the scatterers. An example of P-wave impedance inversion is performed on a VSP data set consisting of three offsets recorded in two wells.

Water temperature studies on the R. North Tyne after impoundment by Kielder dam. 3. Preliminary examination of water temperature data

The report describes the results of preliminary analyses of data obtained from a series of water temperature loggers sited at various distances (0.8 to 21.8 km) downstream of Kielder dam on the River North Tyne and in two natural tributaries. The report deals with three aspects of the water temperature records: An analysis of an operational aspect of the data sets for selected stations, a simple examination of the effects of impoundment upon water temperature at or close to the point of release, relative to natural river temperatures, and an examination of rate of change of monthly means of daily mean, maximum, minimum and range (maximum - minimum) with distance downstream of the point of release during 1983.

Computer system support for data analysis

This thesis is an investigation into the nature of data analysis and computer software systems which support this activity.

The first chapter develops the notion of data analysis as an experimental science which has two major components: data-gathering and theory-building. The basic role of language in determining the meaningfulness of theory is stressed, and the informativeness of a language and data base pair is studied. The static and dynamic aspects of data analysis are then considered from this conceptual vantage point. The second chapter surveys the available types of computer systems which may be useful for data analysis. Particular attention is paid to the questions raised in the first chapter about the language restrictions imposed by the computer system and its dynamic properties.

The third chapter discusses the REL data analysis system, which was designed to satisfy the needs of the data analyzer in an operational relational data system. The major limitation on the use of such systems is the amount of access to data stored on a relatively slow secondary memory. This problem of the paging of data is investigated and two classes of data structure representations are found, each of which has desirable paging characteristics for certain types of queries. One representation is used by most of the generalized data base management systems in existence today, but the other is clearly preferred in the data analysis environment, as conceptualized in Chapter I.

This data representation has strong implications for a fundamental process of data analysis -- the quantification of variables. Since quantification is one of the few means of summarizing and abstracting, data analysis systems are under strong pressure to facilitate the process. Two implementations of quantification are studied: one analagous to the form of the lower predicate calculus and another more closely attuned to the data representation. A comparison of these indicates that the use of the "label class" method results in orders of magnitude improvement over the lower predicate calculus technique.

Precision estimates and suggested sample sizes for length-frequency data

For most fisheries applications, the shape of a length-frequency distribution is much more important than its mean length or variance. This makes it difficult to evaluate at which point a sample size is adequate. By estimating the coefficient of variation of the counts in each length class and taking a weighted mean of these, a measure of precision was obtained that takes the precision in all length classes into account. The precision estimates were closely associated with the ratio of the sample size to the number of size classes in each sample. As a rule-of-thumb, a minimum sample size of 10 times the number of length classes in the sample is suggested because the precision deteriorates rapidly for smaller sample sizes. In absence of such a rule-of-thumb, samplers have previously under-estimated the required sample size for samples with large fish, while over-sampling small fish of the same species.

An ecological characterization of the marine resources of Vieques, Puerto Rico. Part I: Historical data synthesis

From the 1940s until 2003, portions of the island of Vieques, a municipality within the Commonwealth of Puerto Rico, were used by the US Navy as a base and training facility, resulting in development and zoning history that differ in comparison to other Caribbean islands. The majority of former Navy lands are now under the jurisdiction of the Department of the Interior’s Fish and Wildlife Service as a National Wildlife Refuge, while a smaller percentage of land was transferred to the Vieques municipality and the Puerto Rico Conservation Trust. An analysis of the distribution and status of the marine resources is timely in light of the recent land transfer, increases in development and tourism, and potential changes in marine zoning around the island. To meet this need, NOAA’s Biogeography Branch, in cooperation with the Office of Response and Restoration and other local and regional partners, conducted Part I of an ecological characterization to integrate historical data and research into a synthesis report. The overall objective of this report is to provide resource managers and residents a comprehensive characterization of the marine resources of Vieques to support research, monitoring, and management. For example, knowledge of the spatial distribution of physical features, habitats, and biological communities is necessary to make an informed decision of the establishment and placement of a marine protected area (MPA). The report is divided into chapters based on the physical environment (e.g., climate, geology, bathymetry), habitat types (e.g., reefs and hardbottom, seagrasses, mangroves) and major faunal groups (e.g. fish, turtles, birds). Each section includes five subsections: an overview, description of the relevant literature, methods of analysis, information on the distribution, status and trends of the particular resource, and a discussion of ecological linkages with other components of the Vieques marine ecosystem and surrounding environment. The physical environment of Vieques is similar to other islands within the Greater Antilles chain, with some distinctions. The warm, tropical climate of Vieques, mediated by the northeasterly trade winds, is characterized by a dry season (December-April) and a rainy season (May-November), the latter of which is characterized by the occasional passage of tropical cyclones. Compared to mainland Puerto Rico, Vieques is characterized by lower elevation, less annual precipitation, and higher average temperatures. The amount of annual precipitation also varies spatially within Vieques, with the western portion of the island receiving higher amounts of rainfall than further east. While the North Equatorial Current dominates the circulation pattern in the Greater Antilles region, small scale current patterns specific to Vieques are not as well characterized. These physical processes are important factors mitigating the distribution and composition of marine benthic habitats around Vieques. In general, the topography of Vieques is characterized by rolling hills. Mt. Pirata, the tallest point at 301 m, is located near the southwest coast. In the absence of island wide sedimentation measurements, information on land cover, slope, precipitation, and soil type were used to estimate relative erosion potential and sediment delivery for each watershed. While slope and precipitation amount are the primary driving factors controlling runoff, land use practices such as urban development, military activity, road construction, and agriculture can increase the delivery of pollution and sediments to coastal waters. Due to the recent land transfer, increased development and tourism is expected, which may result in changes in the input of sediments to the coastal environment.

A point for comparison: the Flower Garden Banks National Marine Sanctuary

A reoccurring goal listed during the creation of Marine Protected Areas (MPAs) is to return the region to a former state. However, limited data is available that describes or characterizes this former condition. Data collected from ecosystems with comparatively limited anthropogenic impacts, can provide invaluable information in suggesting what former states may have looked like. One example is the Flower Garden Banks National Marine Sanctuary which is located 180 kilometers off the coast of Texas. These relatively isolated and pristine banks are capped by substantial scleractinian coral communities, forming excellent habitat for over 200 species of fish. While fishing is permitted, it is limited by difficulty of access. In 2006, NOAA’s Biogeography Branch, in collaboration with the Sanctuary, initiated the first quantitative assessment of fish resources throughout the diveable portions of the Sanctuary. The sampling design and methodologies employed were identical to those that the Branch has utilized in other more impacted regions of the US Caribbean. Initial analyses reveal that fish density and species richness at the Sanctuary were almost two times greater than that found within the US Caribbean and biomass was approximately six times higher. This was due in large part to the presence of sizeable piscivores of the genera Mycteroperca and Dermatolepis. The Sanctuary is one of few minimally impacted locations remaining within the Tropical Western Atlantic. As such, these findings should be considered when attempting to establish a former state or evaluate effectiveness of an MPA in meeting its management goals.

H∞ sampled-data synthesis and related numerical issues

In this paper we present a new, compact derivation of state-space formulae for the so-called discretisation-based solution of the H∞ sampled-data control problem. Our approach is based on the established technique of continuous time-lifting, which is used to isometrically map the continuous-time, linear, periodically time-varying, sampled-data problem to a discretetime, linear, time-invariant problem. State-space formulae are derived for the equivalent, discrete-time problem by solving a set of two-point, boundary-value problems. The formulae accommodate a direct feed-through term from the disturbance inputs to the controlled outputs of the original plant and are simple, requiring the computation of only a single matrix exponential. It is also shown that the resultant formulae can be easily re-structured to give a numerically robust algorithm for computing the state-space matrices. © 1997 Elsevier Science Ltd. All rights reserved.

Data structures and software tools for the computer aided design of control systems: a survey.

CAD software can be structured as a set of modular 'software tools' only if there is some agreement on the data structures which are to be passed between tools. Beyond this basic requirement, it is desirable to give the agreed structures the status of 'data types' in the language used for interactive design. The ultimate refinement is to have a data management capability which 'understands' how to manipulate such data types. In this paper the requirements of CACSD are formulated from the point of view of Database Management Systems. Progress towards meeting these requirements in both the DBMS and the CACSD community is reviewed. The conclusion reached is that there has been considerable movement towards the realisation of software tools for CACSD, but that this owes more to modern ideas about programming languages, than to DBMS developments. The DBMS field has identified some useful concepts, but further significant progress is expected to come from the exploitation of concepts such as object-oriented programming, logic programming, or functional programming.

Gaussian process change point models

We combine Bayesian online change point detection with Gaussian processes to create a nonparametric time series model which can handle change points. The model can be used to locate change points in an online manner; and, unlike other Bayesian online change point detection algorithms, is applicable when temporal correlations in a regime are expected. We show three variations on how to apply Gaussian processes in the change point context, each with their own advantages. We present methods to reduce the computational burden of these models and demonstrate it on several real world data sets. Copyright 2010 by the author(s)/owner(s).

Parameters identification by a single point free response measurement

When a thin rectangular plate is restrained on the two long edges and free on the remaining edges, the equivalent stiffness of the restraining joints can be identified by the order of the natural frequencies obtained using the free response of the plate at a single location. This work presents a method to identify the equivalent stiffness of the restraining joints, being represented as simply supporting the plate but elastically restraining it in rotation. An integral transform is used to map the autospectrum of the free response from the frequency domain to the stiffness domain in order to identify the equivalent torsional stiffness of the restrained edges of the plate and also the order of natural frequencies. The kernel of the integral transform is built interpolating data from a finite element model of the plate. The method introduced in this paper can also be applied to plates or shells with different shapes and boundary conditions. © 2011 Elsevier Ltd. All rights reserved.