Quantitative Risk Assessment of Listeria monocytogenes in French Cold-Smoked Salmon: I. Quantitative Exposure Assessment

A quantitative assessment of the exposure to Listeria monocytogenes from cold-smoked salmon (CSS) consumption in France is developed. The general framework is a second-order (or two-dimensional) Monte Carlo simulation, which characterizes the uncertainty and variability of the exposure estimate. The model takes into account the competitive bacterial growth between L. monocytogenes and the background competitive flora from the end of the production line to the consumer phase. An original algorithm is proposed to integrate this growth in conditions of varying temperature. As part of a more general project led by the French Food Safety Agency (Afssa), specific data were acquired and modeled for this quantitative exposure assessment model, particularly time-temperature profiles, prevalence data, and contamination-level data. The sensitivity analysis points out the main influence of the mean temperature in household refrigerators and the prevalence of contaminated CSS on the exposure level. The outputs of this model can be used as inputs for further risk assessment.     

Modeling Individual Cell Lag Time Distributions for Listeria monocytogenes

The food industry faces two paradoxical demands: on the one hand, foods need to be microbiologically safe for consumption and on the other hand, consumers want fresh, minimally processed foods. To meet these demands, more insight into the mechanisms of microbial growth is needed, which includes, among others, the microbial lag phase. This is the time needed by bacterial cells to adapt to a new environment (for example, after food product contamination) before starting an exponential growth regime. Since food products are often contaminated with low amounts of pathogenic microorganisms, it is important to know the distribution of these individual cell lag times to make accurate predictions concerning food safety. More precisely, cells with the shortest lag times (i.e., appearing in the left tail of the distribution) are largely decisive for the outgrowth of the population. In this study, an integrated modeling approach is proposed and applied to an existing data set of individual cell lag time measurements of Listeria monocytogenes. In a first step, a logistic modeling approach is applied to predict the fraction of zero-lag cells (which start growing immediately) as a function of temperature, pH, and water activity. For the nonzero-lag cells, the mean and variance of the lag time distribution are modeled with a hyperbolic-type model structure. This mean and variance allow identification of the parameters of a two-parameter Weibull distribution, representing the nonzero-lag cell lag time distribution. The integration of the developed models allows prediction of a global distribution of individual cell lag times for any combination of environmental conditions in the interpolation domain of the original temperature, pH, and water activity settings. The global fitting quality of the model is quantified using several measures indicating that the model gives accurate predictions, erring slightly on the fail-safe side when predicting the shortest lag times.     

Quantitative Risk Assessment of Listeria monocytogenes in French Cold-Smoked Salmon: II. Risk Characterization

A model for the assessment of exposure to Listeria monocytogenes from cold-smoked salmon consumption in France was presented in the first of this pair of articles (Pouillot et al ., 2007, Risk Analysis, 27:683–700). In the present study, the exposure model output was combined with an internationally accepted hazard characterization model, adapted to the French situation, to assess the risk of invasive listeriosis from cold-smoked salmon consumption in France in a second-order Monte Carlo simulation framework. The annual number of cases of invasive listeriosis due to cold-smoked salmon consumption in France is estimated to be 307, with a very large credible interval ([10; 12,453]), reflecting data uncertainty. This uncertainty is mainly associated with the dose-response model. Despite the significant uncertainty associated with the predictions, this model provides a scientific base for risk managers and food business operators to manage the risk linked to cold-smoked salmon contaminated with L. monocytogenes. Under the modeling assumptions, risk would be efficiently reduced through a decrease in the prevalence of L. monocytogenes or better control of the last steps of the cold chain (shorter and/or colder storage during the consumer step), whereas reduction of the initial contamination levels of the contaminated products and improvement in the first steps of the cold chain do not seem to be promising strategies. An attempt to apply the recent risk-based concept of FSO (food safety objective) on this example underlines the ambiguity in practical implementation of the risk management metrics and the need for further elaboration on these concepts.     

Consumer Phase Risk Assessment for Listeria monocytogenes in Deli Meats

The foodborne disease risk associated with the pathogen Listeria monocytogenes has been the subject of recent efforts in quantitative microbial risk assessment. Building upon one of these efforts undertaken jointly by the U.S. Food and Drug Administration and the U.S. Department of Agriculture (USDA), the purpose of this work was to expand on the consumer phase of the risk assessment to focus on handling practices in the home. One-dimensional Monte Carlo simulation was used to model variability in growth and cross-contamination of L. monocytogenes during food storage and preparation of deli meats. Simulations approximated that 0.3% of the servings were contaminated with >10(4) CFU/g of L. monocytogenes at the time of consumption. The estimated mean risk associated with the consumption of deli meats for the intermediate-age population was approximately 7 deaths per 10(11) servings. Food handling in homes increased the estimated mean mortality by 10(6)-fold. Of all the home food-handling practices modeled, inadequate storage, particularly refrigeration temperatures, provided the greatest contribution to increased risk. The impact of cross-contamination in the home was considerably less. Adherence to USDA Food Safety and Inspection Service recommendations for consumer handling of ready-to-eat foods substantially reduces the risk of listeriosis.     

Estimating Listeria monocytogenes Growth in Ready-to-Eat Chicken Salad Using a Challenge Test for Quantitative Microbial Risk Assessment

Currently, there is a growing preference for convenience food products, such as ready-to-eat (RTE) foods, associated with long refrigerated shelf-lives, not requiring a heat treatment prior to consumption. Because Listeria monocytogenes is able to grow at refrigeration temperatures, inconsistent temperatures during production, distribution, and at consumer's household may allow for the pathogen to thrive, reaching unsafe limits. L. monocytogenes is the causative agent of listeriosis, a rare but severe human illness, with high fatality rates, transmitted almost exclusively by food consumption. With the aim of assessing the quantitative microbial risk of L. monocytogenes in RTE chicken salads, a challenge test was performed. Salads were inoculated with a three-strain mixture of cold-adapted L. monocytogenes and stored at 4, 12, and 16 °C for eight days. Results revealed that the salad was able to support L. monocytogenes’ growth, even at refrigeration temperatures. The Baranyi primary model was fitted to microbiological data to estimate the pathogen's growth kinetic parameters. Temperature effect on the maximum specific growth rate (μ_max) was modeled using a square-root-type model. Storage temperature significantly influenced μ_max of L. monocytogenes (p < 0.05). These predicted growth models for L. monocytogenes were subsequently used to develop a quantitative microbial risk assessment, estimating a median number of 0.00008726 listeriosis cases per year linked to the consumption of these RTE salads. Sensitivity analysis considering different time–temperature scenarios indicated a very low median risk per portion (<−7 log), even if the assessed RTE chicken salad was kept in abuse storage conditions.     

Listeria monocytogenes Dose Response Revisited—Incorporating Adjustments for Variability in Strain Virulence and Host Susceptibility

Evaluations of Listeria monocytogenes dose‐response relationships are crucially important for risk assessment and risk management, but are complicated by considerable variability across population subgroups and L. monocytogenes strains. Despite difficulties associated with the collection of adequate data from outbreak investigations or sporadic cases, the limitations of currently available animal models, and the inability to conduct human volunteer studies, some of the available data now allow refinements of the well‐established exponential L. monocytogenes dose response to more adequately represent extremely susceptible population subgroups and highly virulent L. monocytogenes strains. Here, a model incorporating adjustments for variability in L. monocytogenes strain virulence and host susceptibility was derived for 11 population subgroups with similar underlying comorbidities using data from multiple sources, including human surveillance and food survey data. In light of the unique inherent properties of L. monocytogenes dose response, a lognormal‐Poisson dose‐response model was chosen, and proved able to reconcile dose‐response relationships developed based on surveillance data with outbreak data. This model was compared to a classical beta‐Poisson dose‐response model, which was insufficiently flexible for modeling the specific case of L. monocytogenes dose‐response relationships, especially in outbreak situations. Overall, the modeling results suggest that most listeriosis cases are linked to the ingestion of food contaminated with medium to high concentrations of L. monocytogenes. While additional data are needed to refine the derived model and to better characterize and quantify the variability in L. monocytogenes strain virulence and individual host susceptibility, the framework derived here represents a promising approach to more adequately characterize the risk of listeriosis in highly susceptible population subgroups.     

Global Sensitivity Analysis Applied to a Contamination Assessment Model of Listeria monocytogenes in Cold Smoked Salmon at Consumption

In this study, a variance‐based global sensitivity analysis method was first applied to a contamination assessment model of Listeria monocytogenes in cold smoked vacuum packed salmon at consumption. The impact of the choice of the modeling approach (populational or cellular) of the primary and secondary models as well as the effect of their associated input factors on the final contamination level was investigated. Results provided a subset of important factors, including the food water activity, its storage temperature, and duration in the domestic refrigerator. A refined sensitivity analysis was then performed to rank the important factors, tested over narrower ranges of variation corresponding to their current distributions, using three techniques: ANOVA, Spearman correlation coefficient, and partial least squares regression. Finally, the refined sensitivity analysis was used to rank the important factors.     

New Data,Strategies, and Insights for Listeria monocytogenes Dose‐Response Models: Summary of an Interagency Workshop, 2011

Listeria monocytogenes is a leading cause of hospitalization, fetal loss, and death due to foodborne illnesses in the United States. A quantitative assessment of the relative risk of listeriosis associated with the consumption of 23 selected categories of ready‐to‐eat foods, published by the U.S. Department of Health and Human Services and the U.S. Department of Agriculture in 2003, has been instrumental in identifying the food products and practices that pose the greatest listeriosis risk and has guided the evaluation of potential intervention strategies. Dose‐response models, which quantify the relationship between an exposure dose and the probability of adverse health outcomes, were essential components of the risk assessment. However, because of data gaps and limitations in the available data and modeling approaches, considerable uncertainty existed. Since publication of the risk assessment, new data have become available for modeling L. monocytogenes dose‐response. At the same time, recent advances in the understanding of L. monocytogenes pathophysiology and strain diversity have warranted a critical reevaluation of the published dose‐response models. To discuss strategies for modeling L. monocytogenes dose‐response, the Interagency Risk Assessment Consortium (IRAC) and the Joint Institute for Food Safety and Applied Nutrition (JIFSAN) held a scientific workshop in 2011 (details available at http://foodrisk.org/irac/events/ ). The main findings of the workshop and the most current and relevant data identified during the workshop are summarized and presented in the context of L. monocytogenes dose‐response. This article also discusses new insights on dose‐response modeling for L. monocytogenes and research opportunities to meet future needs.     

Impact of Cooking Procedures and Storage Practices at Home on Consumer Exposure to Listeria Monocytogenes and Salmonella Due to the Consumption of Pork Meat

The objective of this research was to analyze the impact of different cooking procedures (i.e., gas hob and traditional static oven) and levels of cooking (i.e., rare, medium, and well‐done) on inactivation of Listeria monocytogenes and Salmonella in pork loin chops. Moreover, the consumer's exposure to both microorganisms after simulation of meat leftover storage at home was assessed. The results showed that well‐done cooking in a static oven was the only treatment able to inactivate the tested pathogens. The other cooking combinations allowed to reach in the product temperatures always ≥73.6 °C, decreasing both pathogens between 6 log₁₀ cfu/g and 7 log₁₀ cfu/g. However, according to simulation results, the few cells surviving cooking treatments can multiply during storage by consumers up to 1 log₁₀ cfu/g, with probabilities of 0.059 (gas hob) and 0.035 (static oven) for L. monocytogenes and 0.049 (gas hob) and 0.031 (static oven) for Salmonella . The key factors affecting consumer exposure in relation to storage practices were probability of pathogen occurrence after cooking, doneness degree, time of storage, and time of storage at room temperature. The results of this study can be combined with prevalence data and dose–response models in risk assessment models and included in guidelines for consumers on practices to be followed to manage cooking of pork meat at home.     

A Risk Assessment Scheme of Infection Transmission Indoors Incorporating the Impact of Resuspension

A new risk assessment scheme was developed to quantify the impact of resuspension to infection transmission indoors. Airborne and surface pathogenic particle concentration models including the effect of two major resuspension scenarios (airflow‐induced particle resuspension [AIPR] and walking‐induced particle resuspension [WIPR]) were derived based on two‐compartment mass balance models and validated against experimental data found in the literature. The inhalation exposure to pathogenic particles was estimated using the derived airborne concentration model, and subsequently incorporated into a dose‐response model to assess the infection risk. Using the proposed risk assessment scheme, the influences of resuspension towards indoor infection transmission were examined by two hypothetical case studies. In the case of AIPR, the infection risk increased from 0 to 0.54 during 0–0.5 hours and from 0.54 to 0.57 during 0.5–4 hours. In the case of WIPR, the infection risk increased from 0 to 0.87 during 0–0.5 hours and from 0.87 to 1 during 0.5–4 hours. Sensitivity analysis was conducted based on the design‐of‐experiments method and showed that the factors that are related to the inspiratory rate of viable pathogens and pathogen virulence have the most significant effect on the infection probability under the occurrence of AIPR and WIPR. The risk assessment scheme could serve as an effective tool for the risk assessment of infection transmission indoors.     

Sensitivity Analysis of a Two‐Dimensional Quantitative Microbiological Risk Assessment: Keeping Variability and Uncertainty Separated

The aim of quantitative microbiological risk assessment is to estimate the risk of illness caused by the presence of a pathogen in a food type, and to study the impact of interventions. Because of inherent variability and uncertainty, risk assessments are generally conducted stochastically, and if possible it is advised to characterize variability separately from uncertainty. Sensitivity analysis allows to indicate to which of the input variables the outcome of a quantitative microbiological risk assessment is most sensitive. Although a number of methods exist to apply sensitivity analysis to a risk assessment with probabilistic input variables (such as contamination, storage temperature, storage duration, etc.), it is challenging to perform sensitivity analysis in the case where a risk assessment includes a separate characterization of variability and uncertainty of input variables. A procedure is proposed that focuses on the relation between risk estimates obtained by Monte Carlo simulation and the location of pseudo‐randomly sampled input variables within the uncertainty and variability distributions. Within this procedure, two methods are used—that is, an ANOVA‐like model and Sobol sensitivity indices—to obtain and compare the impact of variability and of uncertainty of all input variables, and of model uncertainty and scenario uncertainty. As a case study, this methodology is applied to a risk assessment to estimate the risk of contracting listeriosis due to consumption of deli meats.     

Variability in PAH-DNA Adduct Measurements in Peripheral Mononuclear Cells: Implications for Quantitative Cancer Risk Assessment

Biomarkers such as DNA adducts have significant potential to improve quantitative risk assessment by characterizing individual differences in metabolism of genotoxins and DNA repair and accounting for some of the factors that could affect interindividual variation in cancer risk. Inherent uncertainty in laboratory measurements and within-person variability of DNA adduct levels over time are putatively unrelated to cancer risk and should be subtracted from observed variation to better estimate interindividual variability of response to carcinogen exposure. A total of 41 volunteers, both smokers and nonsmokers, were asked to provide a peripheral blood sample every 3 weeks for several months in order to specifically assess intraindividual variability of polycyclic aromatic hydrocarbon (PAH)-DNA adduct levels. The intraindividual variance in PAH-DNA adduct levels, together with measurement uncertainty (laboratory variability and unaccounted for differences in exposure), constituted roughly 30% of the overall variance. An estimated 70% of the total variance was contributed by interindividual variability and is probably representative of the true biologic variability of response to carcinogenic exposure in lymphocytes. The estimated interindividual variability in DNA damage after subtracting intraindividual variability and measurement uncertainty was 24-fold. Inter-individual variance was higher (52-fold) in persons who constitutively lack the Glutathione S-Transferase M1 (GSTM1) gene which is important in the detoxification pathway of PAH. Risk assessment models that do not consider the variability of susceptibility to DNA damage following carcinogen exposure may underestimate risks to the general population, especially for those people who are most vulnerable.     

Risk of Fetal Mortality After Exposure to Listeria monocytogenes Based on Dose-Response Data from Pregnant Guinea Pigs and Primates

One‐third of the annual cases of listeriosis in the United States occur during pregnancy and can lead to miscarriage or stillbirth, premature delivery, or infection of the newborn. Previous risk assessments completed by the Food and Drug Administration/the Food Safety Inspection Service of the U.S. Department of Agriculture/the Centers for Disease Control and Prevention (FDA/USDA/CDC)^{( 1 )} and Food and Agricultural Organization/the World Health Organization (FAO/WHO)^{( 2 )} were based on dose‐response data from mice. Recent animal studies using nonhuman primates^{( 3 , 4 )} and guinea pigs^{( 5 )} have both estimated LD₅₀s of approximately 10⁷ Listeria monocytogenes colony forming units (cfu). The FAO/WHO^{( 2 )} estimated a human LD₅₀ of 1.9 × 10⁶ cfu based on data from a pregnant woman consuming contaminated soft cheese. We reevaluated risk based on dose‐response curves from pregnant rhesus monkeys and guinea pigs. Using standard risk assessment methodology including hazard identification, exposure assessment, hazard characterization, and risk characterization, risk was calculated based on the new dose‐response information. To compare models, we looked at mortality rate per serving at predicted doses ranging from 10^?4 to 10¹² L. monocytogenes cfu. Based on a serving of 10⁶ L. monocytogenes cfu, the primate model predicts a death rate of 5.9 × 10^?1 compared to the FDA/USDA/CDC (fig. IV‐12)^{( 1 )} predicted rate of 1.3 × 10^?7. Based on the guinea pig and primate models, the mortality rate calculated by the FDA/USDA/CDC^{( 1 )} is underestimated for this susceptible population.     

Variability in Physical Constant Parameter Values from Standard Data Sources and the Implication of this Variability for Risk Analysis

Physical property values are used in environmental risk assessments to estimate media and risk-based concentrations. Recently, however, considerable variability has been reported with such values. To evaluate potential variability in physical parameter values supporting a variety of regulatory programs, eight data sources were chosen for evaluation, and chemicals appearing in at least four sources were selected. There were 755 chemicals chosen. In addition, chemicals in seven environmentally important subgroups were also identified for evaluation. Nine parameters were selected for analysis--molecular weight (MolWt), melting point (MeltPt), boiling point (BoilPt), vapor pressure (VP), water solubility (AqSOL), Henry's law constant (HLC), octanol-water partition coefficient (Kow), and diffusion coefficients in air (Dair) and water (Dwater). Results show that while 71% of constituents had equal MolWts across data sources, <3% of the constituents had equivalent parameter values across data sources for AqSOL, VP, or HLC. Considerable dissimilarity between certain sources was also observed. Furthermore, measures of dispersion showed considerable variation in data sets for Kow, VP, AqSOL, and HLC compared to measures for MolWt, MeltPt, BoilPt, or Dwater. The magnitude of the observed variability was also noteworthy. For example, the 95th percentile ratio of maximum/minimum parameter values ranged from 1.0 for MolWt to well over 1.0 x 10(6) for VP and HLC. Risk and exposure metrics also varied by similar magnitudes. Results with environmentally important subgroups were similar. These results show that there is considerable variability in physical parameter values from standard sources, and that the observed variability could affect potential risk estimates and perhaps risk management decisions.     

Probabilistic Modeling of the Fate of Listeria Monocytogenes in Diced Bacon During the Manufacturing Process

To assess the impact of the manufacturing process on the fate of Listeria monocytogenes, we built a generic probabilistic model intended to simulate the successive steps in the process. Contamination evolution was modeled in the appropriate units (breasts, dice, and then packaging units through the successive steps in the process). To calibrate the model, parameter values were estimated from industrial data, from the literature, and based on expert opinion. By means of simulations, the model was explored using a baseline calibration and alternative scenarios, in order to assess the impact of changes in the process and of accidental events. The results are reported as contamination distributions and as the probability that the product will be acceptable with regards to the European regulatory safety criterion. Our results are consistent with data provided by industrial partners and highlight that tumbling is a key step for the distribution of the contamination at the end of the process. Process chain models could provide an important added value for risk assessment models that basically consider only the outputs of the process in their risk mitigation strategies. Moreover, a model calibrated to correspond to a specific plant could be used to optimize surveillance.     

Ecological Risk Assessment Framework for Low-Altitude Aircraft Overflights: I. Planning the Analysis and Estimating Exposure

An ecological risk assessment framework for low-altitude aircraft overflights was developed, with special emphasis on military applications. The problem formulation and exposure analysis phases are presented in this article; an analysis of effects and risk characterization is presented in a companion article. The intent of this article is threefold: (1) to illustrate the development of a generic framework for the ecological risk assessment of an activity, (2) to show how the U.S. Environmental Protection Agency's ecological risk assessment paradigm can be applied to an activity other than the release of a chemical, and (3) to provide guidance for the assessment of ecological risks from low-altitude aircraft overflights. The key stressor for low-altitude aircraft overflights is usually sound, although visual and physical (collision) stressors may also be significant. Susceptible and regulated wildlife populations are the major assessment endpoint entities, although plant communities may be impacted by takeoffs and landings. The exposure analysis utilizes measurements of wildlife locations, measurements of sound levels at the wildlife locations, measurements of slant distances from aircraft to wildlife, models that extrapolate sound from the source aircraft to the ground, and bird-strike probability models. Some of the challenges to conducting a risk assessment for aircraft overflights include prioritizing potential stressors and endpoints, choosing exposure metrics that relate to wildlife responses, obtaining good estimates of sound or distance, and estimating wildlife locations.     

Exposure Assessment: Then, Now, and Quantum Leaps in the Future

Health risk assessments have become so widely accepted in the United States that their conclusions are a major factor in many environmental decisions. Although the risk assessment paradigm is 10 years old, the basic risk assessment process has been used by certain regulatory agencies for nearly 40 years. Each of the four components of the paradigm has undergone significant refinements, particularly during the last 5 years. A recent step in the development of the exposure assessment component can be found in the 1992 EPA Guidelines for Exposure Assessment. Rather than assuming worst-case or hypothetical maximum exposures, these guidelines are designed to lead to an accurate characterization, making use of a number of scientific advances. Many exposure parameters have become better defined, and more sensitive techniques now exist for measuring concentrations of contaminants in the environnment. Statistical procedures for characterizing variability, using Monte Carlo or similar approaches, eliminate the need to select point estimates for all individual exposure parameters. These probabilistic models can more accurately characterize the full range of exposures that may potentially be encountered by a given population at a particular site, reducing the need to select highly conservative values to account for this form of uncertainty in the exposure estimate. Lastly, our awareness of the uncertainties in the exposure assessment as well as our knowledge as to how best to characterize them will almost certainly provide evaluations that will be more credible and, therein, more useful to risk managers. If these refinements are incorporated into future exposure assessments, it is likely that our resources will be devoted to problems that, when resolved, will yield the largest improvement in public health.     

Vero-Cytotoxigenic Escherichia coli O157 in Pasteurized Milk Containers at the Point of Retail: A Qualitative Approach to Exposure Assessment

Dairies within the United Kingdom are classified into two groups, namely, off-farm and on-farm dairies (the latter often being small scale). We propose a model for the probability of milk sold as pasteurized reaching the point of retail contaminated with Vero-cytotoxigenic Escherichia coli (VTEC) O157 from each of these two pathways. We evaluate qualitatively the exposures inherent in each, and compare and contrast the two situations. The model framework is generic, in that it can in principle be used, with the relevant data modifications, to provide a qualitative assessment of the likely exposure from milk sold as pasteurized to any potentially milk-borne pathogenic organism. Furthermore, the methodological approaches presented are widely applicable in the microbial risk assessment field. The specific example presented will be of particular interest to the UK dairy and public health communities: we conclude that the exposure potential per liter consumed from milk processed in off-farm dairies is negligible, whereas the exposure potential per liter consumed from milk processed on-farm is low, but not sufficiently small to be regarded as negligible. We also identify areas of data sparsity, which need to be addressed for quantitative risk assessment to proceed, and highlight the critical points in the pasteurized milk production chain, which, in the event of a breakdown, have the potential to increase the risk to the consumer.     

Site-Specific Applications of Probabilistic Health Risk Assessment: Review of the Literature Since 2000

Whether and to what extent contaminated sites harm ecologic and human health are topics of considerable interest, but also considerable uncertainty. Several federal and state agencies have approved the use of some or many aspects of probabilistic risk assessment (PRA), but its site-specific application has often been limited to high-profile sites and large projects. Nonetheless, times are changing: newly developed software tools, and recent federal and state guidance documents formalizing PRA procedures, now make PRA a readily available method of analysis for even small-scale projects. This article presents and discusses a broad review of PRA literature published since 2000.     

Hazardous Materials Releases in the Northridge Earthquake: Implications for Seismic Risk Assessment

In recent years, increasing attention has been given to the problem of earthquake-initiated hazardous materials releases (EIHRs). While the evidence indicates that EHIRs are an important part of the earthquake hazard profile, little attention has been given to documenting them in a way that could form the basis for more systematic estimates of their probability and consequences. Data from the January 17, 1994 Northridge earthquake document the range and nature of hazardous materials events in the impact area, including those at fixed-site facilities and in transportation systems. Proportional estimates of hazmat incidence are provided, and the implications of these data for developing risk assessments and earthquake hazard management are discussed.