Sampling design for compliance monitoring of surface water quality: A case study in a Polder area

International agreements such as the EU Water Framework Directive (WFD) ask for efficient sampling methods for monitoring natural resources. In this paper a general methodology for designing efficient, statistically sound monitoring schemes is described. An important decision is the choice between a design-based and a model-based method, implying the choice between probability (random) sampling and purposive sampling. For mapping purposes, model-based methods are more appropriate, whereas to obtain valid results for the universe as a whole, such as in testing water quality standards against legal standards, we generally prefer a design-based method. Four basic sampling patterns in space-time universe are described: static, synchronous, static-synchronous, and rotational. A case study is carried out for monitoring the quality of surface water at two farms in western Netherlands, wherein a synchronous sampling design is applied, with stratified simple random sampling in both space and time. To reduce laboratory costs the aliquots taken at the locations of a given sampling round are bulked to form a composite. To test the spatiotemporal mean N-total concentration during the summer half-year against the MAR standard with a power of 80% at a concentration 15% below the MAR standard and with a confidence of 95%, six to nine sampling rounds are needed with 50 to 75 locations per sampling round. For P-total the required number of sampling rounds differs strongly between the two farms, but is for both farms much larger than for N-total

Making Big Data Useful for Health Care: A Summary of the Inaugural MIT Critical Data Conference

With growing concerns that big data will only augment the problem of unreliable research, the Laboratory of Computational Physiology at the Massachusetts Institute of Technology organized the Critical Data Conference in January 2014. Thought leaders from academia, government, and industry across disciplines--including clinical medicine, computer science, public health, informatics, biomedical research, health technology, statistics, and epidemiology--gathered and discussed the pitfalls and challenges of big data in health care. The key message from the conference is that the value of large amounts of data hinges on the ability of researchers to share data, methodologies, and findings in an open setting. If empirical value is to be from the analysis of retrospective data, groups must continuously work together on similar problems to create more effective peer review. This will lead to improvement in methodology and quality, with each iteration of analysis resulting in more reliability

Azimuthal anisotropy of pi^0 and eta mesons in Au+Au collisions at sqrt(s_NN)=200 GeV

The azimuthal anisotropy coefficients v_2 and v_4 of pi^0 and eta mesons are measured in Au+Au collisions at sqrt(s_NN)=200 GeV, as a function of transverse momentum p_T (1-14 GeV/c) and centrality. The extracted v_2 coefficients are found to be consistent between the two meson species over the measured p_T range. The ratio of v_4/v_2^2 for pi^0 mesons is found to be independent of p_T for 1-9 GeV/c, implying a lack of sensitivity of the ratio to the change of underlying physics with p_T. Furthermore, the ratio of v_4/v_2^2 is systematically larger in central collisions, which may reflect the combined effects of fluctuations in the initial collision geometry and finite viscosity in the evolving medium.Comment: 384 authors, 71 institutions, 11 pages, 9 figures, and 2 tables. Submitted to Physical Review C. Plain text data tables for the points plotted in figures for this and previous PHENIX publications are (or will be) publicly available at http://www.phenix.bnl.gov/papers.htm

Transverse-energy distributions at midrapidity in pp++pp, dd++Au, and Au++Au collisions at sNN=62.4\sqrt{s_{_{NN}}}=62.4--200~GeV and implications for particle-production models

Measurements of the midrapidity transverse energy distribution, d\Et/d\eta, are presented for $p$$+$$p$, $d$$+$Au, and Au$+$Au collisions at $\sqrt{s_{_{NN}}}=200$ GeV and additionally for Au$+$Au collisions at $\sqrt{s_{_{NN}}}=62.4$ and 130 GeV. The d\Et/d\eta distributions are first compared with the number of nucleon participants $N_{\rm part}$, number of binary collisions $N_{\rm coll}$, and number of constituent-quark participants $N_{qp}$ calculated from a Glauber model based on the nuclear geometry. For Au$+$Au, \mean{d\Et/d\eta}/N_{\rm part} increases with $N_{\rm part}$, while \mean{d\Et/d\eta}/N_{qp} is approximately constant for all three energies. This indicates that the two component ansatz, $dE_{T}/d\eta \propto (1-x) N_{\rm part}/2 + x N_{\rm coll}$, which has been used to represent $E_T$ distributions, is simply a proxy for $N_{qp}$, and that the $N_{\rm coll}$ term does not represent a hard-scattering component in $E_T$ distributions. The $dE_{T}/d\eta$ distributions of Au$+$Au and $d$$+$Au are then calculated from the measured $p$$+$$p$ $E_T$ distribution using two models that both reproduce the Au$+$Au data. However, while the number-of-constituent-quark-participant model agrees well with the $d$$+$Au data, the additive-quark model does not.Comment: 391 authors, 24 pages, 19 figures, and 15 Tables. Submitted to Phys. Rev. C. Plain text data tables for the points plotted in figures for this and previous PHENIX publications are publicly available at http://www.phenix.bnl.gov/papers.htm

Centrality categorization for R_{p(d)+A} in high-energy collisions

High-energy proton- and deuteron-nucleus collisions provide an excellent tool for studying a wide array of physics effects, including modifications of parton distribution functions in nuclei, gluon saturation, and color neutralization and hadronization in a nuclear environment, among others. All of these effects are expected to have a significant dependence on the size of the nuclear target and the impact parameter of the collision, also known as the collision centrality. In this article, we detail a method for determining centrality classes in p(d)+A collisions via cuts on the multiplicity at backward rapidity (i.e., the nucleus-going direction) and for determining systematic uncertainties in this procedure. For d+Au collisions at sqrt(s_NN) = 200 GeV we find that the connection to geometry is confirmed by measuring the fraction of events in which a neutron from the deuteron does not interact with the nucleus. As an application, we consider the nuclear modification factors R_{p(d)+A}, for which there is a potential bias in the measured centrality dependent yields due to auto-correlations between the process of interest and the backward rapidity multiplicity. We determine the bias correction factor within this framework. This method is further tested using the HIJING Monte Carlo generator. We find that for d+Au collisions at sqrt(s_NN)=200 GeV, these bias corrections are small and vary by less than 5% (10%) up to p_T = 10 (20) GeV. In contrast, for p+Pb collisions at sqrt(s_NN) = 5.02 TeV we find these bias factors are an order of magnitude larger and strongly p_T dependent, likely due to the larger effect of multi-parton interactions.Comment: 375 authors, 18 pages, 16 figures, 4 tables. Submitted to Phys. Rev. C. Plain text data tables for the points plotted in figures for this and previous PHENIX publications are (or will be) publicly available at http://www.phenix.bnl.gov/papers.htm

Azimuthally anisotropic emission of low-momentum direct photons in Au++Au collisions at sNN=200\sqrt{s_{_{NN}}}=200 GeV

The PHENIX experiment at the Relativistic Heavy Ion Collider has measured 2nd and 3rd order Fourier coefficients of the azimuthal distributions of direct photons emitted at midrapidity in Au$+$Au collisions at $\sqrt{s_{_{NN}}}=200$ GeV for various collision centralities. Combining two different analysis techniques, results were obtained in the transverse momentum range of $0.4<p_{T}<4.0$ GeV/$c$. At low $p_T$ the second-order coefficients, $v_2$, are similar to the ones observed in hadrons. Third order coefficients, $v_3$, are nonzero and almost independent of centrality. These new results on $v_2$ and $v_3$, combined with previously published results on yields, are compared to model calculations that provide yields and asymmetries in the same framework. Those models are challenged to explain simultaneously the observed large yield and large azimuthal anisotropies.Comment: 552 authors, 15 pages, 9 figures, 3 tables, 2007 and 2010 data. v2 is version accepted for publication by Phys. Rev. C. Plain text data tables for the points plotted in figures for this and previous PHENIX publications are (or will be) publicly available at http://www.phenix.bnl.gov/papers.htm

Double Spin Asymmetry of Electrons from Heavy Flavor Decays in p+p Collisions at sqrt(s)=200 GeV

We report on the first measurement of double-spin asymmetry, A_LL, of electrons from the decays of hadrons containing heavy flavor in longitudinally polarized p+p collisions at sqrt(s)=200 GeV for p_T= 0.5 to 3.0 GeV/c. The asymmetry was measured at mid-rapidity (|eta|<0.35) with the PHENIX detector at the Relativistic Heavy Ion Collider. The measured asymmetries are consistent with zero within the statistical errors. We obtained a constraint for the polarized gluon distribution in the proton of |Delta g/g(log{_10}x= -1.6^+0.5_-0.4, {mu}=m_T^c)|^2 < 0.033 (1 sigma), based on a leading-order perturbative-quantum-chromodynamics model, using the measured asymmetry.Comment: 385 authors, 17 pages, 15 figures, 5 tables. Submitted to Phys. Rev. D. Plain text data tables for the points plotted in figures for this and previous PHENIX publications are (or will be) publicly available at http://www.phenix.bnl.gov/papers.htm

Nuclear matter effects on J/ψJ/\psi production in asymmetric Cu+Au collisions at sNN\sqrt{s_{_{NN}}} = 200 GeV

We report on $J/\psi$ production from asymmetric Cu+Au heavy-ion collisions at $\sqrt{s_{_{NN}}}$=200 GeV at the Relativistic Heavy Ion Collider at both forward (Cu-going direction) and backward (Au-going direction) rapidities. The nuclear modification of $J/\psi$ yields in Cu$+$Au collisions in the Au-going direction is found to be comparable to that in Au$+$Au collisions when plotted as a function of the number of participating nucleons. In the Cu-going direction, $J/\psi$ production shows a stronger suppression. This difference is comparable in magnitude and has the same sign as the difference expected from shadowing effects due to stronger low-$x$ gluon suppression in the larger Au nucleus. The relative suppression is opposite to that expected from hot nuclear matter dissociation, since a higher energy density is expected in the Au-going direction.Comment: 349 authors, 10 pages, 4 figures, and 4 tables. Submitted to Phys. Rev. C. For v2, fixed LaTeX error in 3rd-to-last sentence. Plain text data tables for the points plotted in figures for this and previous PHENIX publications are (or will be) publicly available at http://www.phenix.bnl.gov/papers.htm

Transverse energy production and charged-particle multiplicity at midrapidity in various systems from sNN=7.7\sqrt{s_{NN}}=7.7 to 200 GeV

Measurements of midrapidity charged particle multiplicity distributions, $dN_{\rm ch}/d\eta$, and midrapidity transverse-energy distributions, $dE_T/d\eta$, are presented for a variety of collision systems and energies. Included are distributions for Au$+$Au collisions at $\sqrt{s_{_{NN}}}=200$, 130, 62.4, 39, 27, 19.6, 14.5, and 7.7 GeV, Cu$+$Cu collisions at $\sqrt{s_{_{NN}}}=200$ and 62.4 GeV, Cu$+$Au collisions at $\sqrt{s_{_{NN}}}=200$ GeV, U$+$U collisions at $\sqrt{s_{_{NN}}}=193$ GeV, $d$$+$Au collisions at $\sqrt{s_{_{NN}}}=200$ GeV, $^{3}$He$+$Au collisions at $\sqrt{s_{_{NN}}}=200$ GeV, and $p$$+$$p$ collisions at $\sqrt{s_{_{NN}}}=200$ GeV. Centrality-dependent distributions at midrapidity are presented in terms of the number of nucleon participants, $N_{\rm part}$, and the number of constituent quark participants, $N_{q{\rm p}}$. For all $A$$+$$A$ collisions down to $\sqrt{s_{_{NN}}}=7.7$ GeV, it is observed that the midrapidity data are better described by scaling with $N_{q{\rm p}}$ than scaling with $N_{\rm part}$. Also presented are estimates of the Bjorken energy density, $\varepsilon_{\rm BJ}$, and the ratio of $dE_T/d\eta$ to $dN_{\rm ch}/d\eta$, the latter of which is seen to be constant as a function of centrality for all systems.Comment: 706 authors, 32 pages, 20 figures, 34 tables, 2004, 2005, 2008, 2010, 2011, and 2012 data. v2 is version accepted for publication in Phys. Rev.

Cross section for bbˉb\bar{b} production via dielectrons in d++Au collisions at sNN=200\sqrt{s_{_{NN}}}=200 GeV

We report a measurement of $e^+e^-$ pairs from semileptonic heavy-flavor decays in $d$$+$Au collisions at $\sqrt{s_{_{NN}}}=200$ GeV. Exploring the mass and transverse-momentum dependence of the yield, the bottom decay contribution can be isolated from charm, and quantified by comparison to {\sc pythia} and {\sc mc@nlo} simulations. The resulting $b\bar{b}$-production cross section is $\sigma^{d{\rm Au}}_{b\bar{b}}=1.37{\pm}0.28({\rm stat}){\pm}0.46({\rm syst})$~mb, which is equivalent to a nucleon-nucleon cross section of $\sigma^{NN}_{bb}=3.4\pm0.8({\rm stat}){\pm}1.1({\rm syst})\ \mu$b.Comment: 375 authors, 16 pages, 8 figures, 7 tables, 2008 data. Submitted to Phys. Rev. C Plain text data tables for the points plotted in figures for this and previous PHENIX publications are (or will be) publicly available at http://www.phenix.bnl.gov/papers.htm