Hierarchical Bayesian Approach to Boundary Value Problems with Stochastic Boundary Conditions

This is the pre-print version of the article found in the Monthly Weather Review (http://journals.ametsoc.org/toc/mwre/138/10).Boundary value problems are ubiquitous in the atmospheric and ocean sciences. Typical settings include bounded, partially bounded, global and limited area domains, discretized for applications of numerical models of the relevant fluid equations. Often, limited area models are constructed to interpret intensive datasets collected over a specific region, from a variety of observational platforms. These data are noisy and they typically do not span the domain of interest uniformly in space and time. Traditional numerical procedures cannot easily account for these uncertainties. A hierarchical Bayesian modeling framework is developed for solving boundary value problems in such settings. By allowing the boundary process to be stochastic, and conditioning the interior process on this boundary, one can account for the uncertainties in the boundary process in a reasonable fashion. In the presence of data and all its uncertainties, this idea can be related through Bayes' Theorem to produce distributions of the interior process given the observational data. The method is illustrated with an example of obtaining atmospheric streamfunction fields in the Labrador Sea region, given scatterometer-derived observations of the surface wind field

Spatio-Temporal Hierarchical Bayesian Modeling: Tropical Ocean Surface Winds

This is the author's version of the article found in the Journal of the American Statistical Association. The publisher's version can be found at http://pubs.amstat.org/loi/jasa.Spatio-temporal processes are ubiquitous in the environmental and physical sciences. This is certainly true of atmospheric and oceanic processes, which typically exhibit many different scales of spatial and temporal variability. The complexity of these processes and large number of observation/prediction locations preclude the use of traditional covariance-based space-time statistical methods. Alternatively, we focus on conditionally-specified (i.e., hierarchical) spatio-temporal models. These methods offer several advantages over traditional approaches. Primarily, physical and dynamical constraints are easily incorporated into the conditional formulation, so that the series of relatively simple, yet physically realistic, conditional models leads to a much more complicated space-time covariance structure than can be specified directly. Furthermore, by making use of the sparse structure inherent in the hierarchical approach, as well as multiresolution (wavelet) bases, the models are computable with very large datasets. This modeling approach was necessitated by a scientifically meaningful problem in the geosciences. Satellite-derived wind estimates have high spatial resolution but are limited in global coverage. In contrast, wind fields provided by the major weather centers provide complete coverage but have low spatial resolution. The goal is to combine these data in a manner that incorporates the space-time dynamics inherent in the surface wind field. This is an essential task to enable meteorological research as no complete high resolution surface wind datasets exist over the world oceans. High resolution datasets of this kind are crucial for improving our understanding of: global air-sea interactions affecting climate, tropical disturbances, and for driving large-scale ocean circulation models.Support for this research was provided for CKW, DN, and LMB by the NCAR Geophysical Statistics Project, sponsored by the National Science Foundation (NSF) under Grant DMS93-12686. Support for RFM and CKW is provided by the NCAR NSCAT Science Working Team cooperative agreement with NASA JPL. NCAR is supported in part by the NSF

Ocean ensemble forecasting. Part I: Ensemble Mediterranean winds from a Bayesian hierarchical model

A Bayesian hierarchical model (BHM) is developed to estimate surface vector wind (SVW) fields and associated uncertainties over the Mediterranean Sea. The BHM–SVW incorporates data-stage inputs from analyses and forecasts of the European Centre for Medium-Range Weather Forecasts (ECMWF) and SVW retrievals from the QuikSCAT data record. The process-model stage of the BHM–SVW is based on a Rayleigh friction equation model for surface winds. Dynamical interpretations of posterior distributions of the BHM–SVW parameters are discussed. Ten realizations from the posterior distribution of the BHM–SVW are used to force the data-assimilation step of an experimental ensemble ocean forecast system for the Mediterranean Sea in order to create a set of ensemble initial conditions. The sequential data-assimilation method of the Mediterranean forecast system (MFS) is adapted to the ensemble implementation. Analyses of sample ensemble initial conditions for a single data-assimilation period in MFS are presented to demonstrate the multivariate impact of the BHM–SVW ensemble generation methodology. Ensemble initial-condition spread is quantified by computing standard deviations of ocean state variable fields over the ten ensemble members. The methodological findings in this article are of two kinds. From the perspective of statistical modelling, the process-model development is more closely related tophysicalbalances than inpreviousworkwithmodels for the SVW.Fromthe ocean forecast perspective, the generation of ocean ensemble initial conditions via BHM is shown to be practical for operational implementation in an ensemble ocean forecast system. Phenomenologically, ensemble spread generated via BHM–SVW occurs on ocean mesoscale time- and space-scales, in close association with strong synoptic-scale wind-forcing events. A companion article describes the impacts of the BHM–SVW ensemble method on the ocean forecast in comparisons with more traditional ensemble methods

A Retrospective Multicenter Analysis of the Incidence of Bone-Only Disease at PSMA PET/CT in Castration Resistant Prostate Cancer Patients

PSMA PET/CT has unprecedented accuracy for localization of initial or recurrent prostate cancer (PC), which can be applied in a metastasis-directed therapy approach. PSMA PET/CT (PET) also has a role in the selection of patients for metastasis-directed therapy or radioligand therapy and therapy assessment in CRPC patients. The purpose of this multicenter retrospective study was to determine the incidence of bone-only metastasis in CRPC patients who underwent PSMA PET/CT for restaging, as well as identifying potential predictors of bone-only PET positivity. The study analyzed data from 179 patients from two centers in Essen and Bologna. Results showed that 20.1% of the patients had PSMA uptake only in the bone, with the most frequent lesions located in the vertebrae, ribs, and hip bone. Half half of the patients showed oligo disease in bone and may benefit from a bone-metastasis-directed therapy. Initial positive nodal status and solitary ADT were shown to be negative predictors of osseous metastasis. The role of PSMA PET/TC in this patient population needs to be further explored in terms of its role in the evaluation and adoption of bone-specific therapies

Ocean ensemble forecasting. Part II: Mediterranean Forecast System response

This article analyzes the ocean forecast response to surface vector wind (SVW) distributions generated by a Bayesian hierarchical model (BHM) developed in Part I of this series. A new method for ocean ensemble forecasting (OEF), the socalled BHM-SVW-OEF, is described. BHM-SVW realizations are used to produce and force perturbations in the ocean state during 14 day analysis and 10 day forecast cycles of the Mediterranean Forecast System (MFS). The BHM-SVW-OEF ocean response spread is amplified at the mesoscales and in the pycnocline of the eddy field. The new method is compared with an ensemble response forced by European Centre for Medium-Range Weather Forecasts (ECMWF) ensemble prediction system (EEPS) surface winds, and with an ensemble forecast started from perturbed initial conditions derived froman ad hoc thermocline intensified random perturbation (TIRP) method. The EEPS-OEF shows spread on basin scales while the TIRP-OEF response is mesoscale-intensified as in the BHM-SVW-OEF response. TIRP-OEF perturbations fill more of the MFS domain, while the BHM-SVW-OEF perturbations are more location-specific, concentrating ensemble spread at the sites where the ocean-model response to uncertainty in the surface wind forcing is largest

A study of the deep structure of the energy landscape of glassy polystyrene: the exponential distribution of the energy-barriers revealed by high-field Electron Spin Resonance spectroscopy

The reorientation of one small paramagnetic molecule (spin probe) in glassy polystyrene (PS) is studied by high-field Electron Spin Resonance spectroscopy at two different Larmor frequencies (190 and 285 GHz). The exponential distribution of the energy-barriers for the rotational motion of the spin probe is unambigously evidenced at both 240K and 270K. The same shape for the distribution of the energy-barriers of PS was evidenced by the master curves provided by previous mechanical and light scattering studies. The breadth of the energy-barriers distribution of the spin probe is in the range of the estimates of the breadth of the PS energy-barriers distribution. The evidence that the deep structure of the energy landscape of PS exhibits the exponential shape of the energy-barriers distribution agrees with results from extreme-value statistics and the trap model by Bouchaud and coworkers.Comment: Final version in press as Letter to the Editor on J.Phys.:Condensed Matter. Changes in bol

Ocean Ensemble Forecasting, Part II: Mediterranean Forecast System Response

This paper analyzes the ocean forecast response to surface vector wind (SVW) distributions generated by a Bayesian Hierarchical Model (BHM) developed in Part I (Milliff et al., 2009). A new method for Ocean Ensemble Forecasting (OEF), so-called BHM-SVW-OEF, is described. BHM-SVW realizations are used to produce and force perturbations in the ocean state during 14-day analysis and 10-day forecast cycles of the Mediterranean Forecast System (MFS). The BHM-SVW-OEF ocean response spread is amplified at the mesoscales and pycnocline of the eddy field. The new method is compared to an ensemble response forced by ECMWF Ensemble Prediction System (EEPS) surface winds, and to an ensemble forecast started from perturbed initial conditions derived from an ad hoc Thermocline Intensified Random Perturbation (TIRP) method. The EEPS-OEF shows spread at the basin scales while the TIRP-OEF response is mesoscale intensified as in the BHM-SVW-OEF response. TIRP-OEF perturbations fill more of the MFS domain while the BHM-SVW-OEF perturbations are more location-specific, concentrating ensemble spread at the sites where the ocean model response to uncertainty in the surface wind forcing is largest. The BHM-SVW-OEF method offers a practical and objective means for producing short-term forecast spread by modeling surface atmospheric forcing uncertainties that have maximum impact at the ocean mesoscales

Improving practice in community-based settings: a randomized trial of supervision – study protocol

Background: Evidence-based treatments for child mental health problems are not consistently available in public mental health settings. Expanding availability requires workforce training. However, research has demonstrated that training alone is not sufficient for changing provider behavior, suggesting that ongoing intervention-specific supervision or consultation is required. Supervision is notably under-investigated, particularly as provided in public mental health. The degree to which supervision in this setting includes ‘gold standard’ supervision elements from efficacy trials (e.g., session review, model fidelity, outcome monitoring, skill-building) is unknown. The current federally-funded investigation leverages the Washington State Trauma-focused Cognitive Behavioral Therapy Initiative to describe usual supervision practices and test the impact of systematic implementation of gold standard supervision strategies on treatment fidelity and clinical outcomes. Methods/Design The study has two phases. We will conduct an initial descriptive study (Phase I) of supervision practices within public mental health in Washington State followed by a randomized controlled trial of gold standard supervision strategies (Phase II), with randomization at the clinician level (i.e., supervisors provide both conditions). Study participants will be 35 supervisors and 130 clinicians in community mental health centers. We will enroll one child per clinician in Phase I (N = 130) and three children per clinician in Phase II (N = 390). We use a multi-level mixed within- and between-subjects longitudinal design. Audio recordings of supervision and therapy sessions will be collected and coded throughout both phases. Child outcome data will be collected at the beginning of treatment and at three and six months into treatment. Discussion This study will provide insight into how supervisors can optimally support clinicians delivering evidence-based treatments. Phase I will provide descriptive information, currently unavailable in the literature, about commonly used supervision strategies in community mental health. The Phase II randomized controlled trial of gold standard supervision strategies is, to our knowledge, the first experimental study of gold standard supervision strategies in community mental health and will yield needed information about how to leverage supervision to improve clinician fidelity and client outcomes. Trial registration ClinicalTrials.gov NCT0180026

The reproductive capacity of rams

Includes bibliographical references (pages 123-127)