Search NASA⌕ Search

SEARCH · Search NASA

Results for “Multivariate sampling”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 217 records · Page 12

Characterization of Arsenic and Selenium in Coal Fly Ash to Improve Evaluations for Disposal and Reuse Potential (Final Technical Report)

Coal fly ash is a high volume waste material that is discarded in landfills and surface water impoundments across the U.S. and is also widely recycled for a variety of applications. The leaching of potential of contaminants of concern, such as arsenic (As) and selenium (Se), is often the driver of risk assessments for coal ash disposal and reuse. The extent of leachable As and Se depends on several factors related to environmental conditions and fly ash characteristics. Previous studies employed various methods to delineate the concentration, chemical form, and distribution of As and Se in fly ash materials. However, few studies have attempted to directly correlate these properties to mobilization parameters relevant to disposal and reuse. Instead, the coal residuals industries often rely upon standardized leaching protocols that can be laborious or involve hazardous chemicals. The goals of the project were to: 1) Develop and evaluate a characterization protocol that can be used to screen fly ash samples for leachability of As and Se; 2) Characterize As, Se, and associated constituents of fly ash particles at multiple length scales (nanometer to micrometer) to determine if elemental associations differ as a function of the resolution of characterization; and 3) Establish a predictive model for the chemical composition of coal ash produced annually at major U.S. coal fired power facilities on 50-year national coal supply records. For the first objective, we performed leaching experiments with 52 fly ash samples collected from 15 different U.S. power plants and representing coal feedstocks from the three major domestic coal regions. For this work, we assessed the mobilization potential of As and Se in fly ash based on standardized leaching protocols and performed multivariate and lasso regression analyses to explore correlations of leachable As and Se contents with characteristics such as major element contents, loss on ignition (LOI) and pH. The results of regression models indicated that major elements (Fe, Ca, Al) for a wide range of fly ashes can serve as predictor variables for the leaching potential of As, but not for Se. LOI and pH were not important predictive variables in the models. Both regression approaches resulted in relatively strong fits for leachable As (correlation coefficient R 2 = 0.78 for both models) compared to models for leachable Se (R 2 = 0.49). Overall, these results suggest that correlation models combined with on-site elemental analysis with portable analyzers may enable a screening method for leachable As in coal ash. For the second objective, we utilized nanoscale 2-D imaging (30-50 nm spot size) with the Hard X-ray Nanoprobe (HXN) in combination with microprobe X-ray capabilities (~5 µm resolution) to determine As and Se elemental associations in fly ash particles. Speciation of As and Se was also measured at the nano- to microscale with X-ray absorption spectroscopy. The enhanced resolution of HXN showed As and Se that were diffusely located around or comingled with Ca- and Fe-rich particles. The results also showed nanoparticles of Se attached to the surface of fly ash grains. Overall, a comparison of As and Se species across scales highlights the heterogeneity and complexity of chemical associations for these trace elements of concern in coal fly ash. For the final objective, we developed a predictive model for major element composition of coal ash in reserve at disposal sites of major U.S. coal fired power plants. This model was constructed from coal purchase records of 705 power stations from 1973-2022 and was trained on coal ash composition data showing that coal ash elemental composition is strongly associated with the source of feedstock coal. The model showed regional shifts in the major element contents of ash produced by power plants in the last 50 years, particularly for calcium and iron (expressed as %CaO and %Fe 2 O 3 ), as coal-fired power stations changed their source of coal over this time frame. Our approach enables an estimation of coal ash chemical composition that is stored in waste impoundments at individual power stations. Such information can help delineate the regional market potential for material applications that would utilize coal ash harvested from disposal sites across the U.S.

01 COAL, LIGNITE, AND PEAT↗

Reticular Materials and AI-Driven Computer Simulations for Seawater Mining of Valuable Metals (Final Technical Report)

This Final Technical Report describes our exploratory efforts that combine reticular materials synthesis (hydrolytically robust metal–organic frameworks, MOFs) with AI‑enabled molecular simulations to develop mechanistic, quantitative design rules for recovering lithium and other alkali-metal ions from highly dilute, competitive aqueous resources (e.g., seawater). The central outcome is a joint experimental–computational study of ion uptake in MOF‑808 (Chemical Science, 2025) that quantifies both thermodynamics and kinetics of Li + , Na + , and K + uptake and identifies how pore size, pore hydration state, dehydration penalties, and pore-window transport barriers govern selectivity. Guided by these insights, we synthesized and tested functionalized MOF‑808 and multivariate MOFs incorporating ion-recognition motifs (including carboxylates and crown-ether linkers) and evaluated uptake in synthetic seawater, highlighting framework topology and pore chemistry as levers for improved Li + /Na + discrimination. We also developed transferable simulation models, enhanced-sampling protocols, and automated workflows that enable systematic screening of porous sorbents.

42 ENGINEERING↗

Modeling of the WSTF frictional heating apparatus in high pressure systems

In order to develop a computer program able to model the frictional heating of metals in high pressure oxygen or nitrogen a number of additions have been made to the frictional heating model originally developed for tests in low pressure helium. These additions include: (1) a physical property package for the gases to account for departures from the ideal gas state; (2) two methods for spatial discretization (finite differences with quadratic interpolation or orthogonal collocation on finite elements) which substantially reduce the computer time required to solve the transient heat balance; (3) more efficient programs for the integration of the ordinary differential equations resulting from the discretization of the partial differential equations; and (4) two methods for determining the best-fit parameters via minimization of the mean square error (either a direct search multivariable simplex method or a modified Levenburg-Marquardt algorithm). The resulting computer program has been shown to be accurate, efficient and robust for determining the heat flux or friction coefficient vs. time at the interface of the stationary and rotating samples.

Skowlund, Christopher T.↗

Probabilistic Modeling of the Renal Stone Formation Module

The Integrated Medical Model (IMM) is a probabilistic tool, used in mission planning decision making and medical systems risk assessments. The IMM project maintains a database of over 80 medical conditions that could occur during a spaceflight, documenting an incidence rate and end case scenarios for each. In some cases, where observational data are insufficient to adequately define the inflight medical risk, the IMM utilizes external probabilistic modules to model and estimate the event likelihoods. One such medical event of interest is an unpassed renal stone. Due to a high salt diet and high concentrations of calcium in the blood (due to bone depletion caused by unloading in the microgravity environment) astronauts are at a considerable elevated risk for developing renal calculi (nephrolithiasis) while in space. Lack of observed incidences of nephrolithiasis has led HRP to initiate the development of the Renal Stone Formation Module (RSFM) to create a probabilistic simulator capable of estimating the likelihood of symptomatic renal stone presentation in astronauts on exploration missions. The model consists of two major parts. The first is the probabilistic component, which utilizes probability distributions to assess the range of urine electrolyte parameters and a multivariate regression to transform estimated crystal density and size distributions to the likelihood of the presentation of nephrolithiasis symptoms. The second is a deterministic physical and chemical model of renal stone growth in the kidney developed by Kassemi et al. The probabilistic component of the renal stone model couples the input probability distributions describing the urine chemistry, astronaut physiology, and system parameters with the physical and chemical outputs and inputs to the deterministic stone growth model. These two parts of the model are necessary to capture the uncertainty in the likelihood estimate. The model will be driven by Monte Carlo simulations, continuously randomly sampling the probability distributions of the electrolyte concentrations and system parameters that are inputs into the deterministic model. The total urine chemistry concentrations are used to determine the urine chemistry activity using the Joint Expert Speciation System (JESS), a biochemistry model. Information used from JESS is then fed into the deterministic growth model. Outputs from JESS and the deterministic model are passed back to the probabilistic model where a multivariate regression is used to assess the likelihood of a stone forming and the likelihood of a stone requiring clinical intervention. The parameters used to determine to quantify these risks include: relative supersaturation (RS) of calcium oxalate, citrate/calcium ratio, crystal number density, total urine volume, pH, magnesium excretion, maximum stone width, and ureteral location. Methods and Validation: The RSFM is designed to perform a Monte Carlo simulation to generate probability distributions of clinically significant renal stones, as well as provide an associated uncertainty in the estimate. Initially, early versions will be used to test integration of the components and assess component validation and verification (V&V), with later versions used to address questions regarding design reference mission scenarios. Once integrated with the deterministic component, the credibility assessment of the integrated model will follow NASA STD 7009 requirements.

kidney stones↗

Metagenomic insights into the taxonomy, function, and dysbiosis of prokaryotic communities in octocorals

Background. In octocorals (Cnidaria Octocorallia), the functional relationship between host health and its symbiotic consortium has yet to be determined. Here, we employed comparative metagenomics to uncover the distinct functional and phylogenetic features of the microbiomes of healthy Eunicella gazella, Eunicella verrucosa, and Leptogorgia sarmentosa tissues, in contrast with the microbiomes found in seawater and sediments. We further explored how the octocoral microbiome shifts to a pathobiome state in E. gazella. Results. Multivariate analyses based on 16S rRNA genes, Clusters of Orthologous Groups of proteins (COGs), Protein families (Pfams), and secondary metabolite-biosynthetic gene clusters annotated from 20 Illumina-sequenced metagenomes each revealed separate clustering of the prokaryotic communities of healthy tissue samples of the three octocoral species from those of necrotic E. gazella tissue and surrounding environments. While the healthy octocoral microbiome was distinguished by so-far uncultivated Endozoicomonadaceae, Oceanospirillales, and Alteromonadales phylotypes in all host species, a pronounced increase of Flavobacteriaceae and Alphaproteobacteria, originating from seawater, was observed in necrotic E. gazella tissue. Increased abundances of eukaryotic-like proteins, exonucleases, restriction endonucleases, CRISPR/Cas proteins, and genes encoding for heat-shock proteins, inorganic ion transport, and iron storage distinguished the prokaryotic communities of healthy octocoral tissue regardless of the host species. An increase of arginase and nitric oxide reductase genes, observed in necrotic E. gazella tissues, suggests the existence of a mechanism for suppression of nitrite oxide production by which octocoral pathogens may overcome the host’s immune system. Conclusions. This is the first study to employ primer-less, shotgun metagenome sequencing to unveil the taxonomic, functional, and secondary metabolism features of prokaryotic communities in octocorals. Our analyses reveal that the octocoral microbiome is distinct from those of the environmental surroundings, is host genus (but not species) specific, and undergoes large, complex structural changes in the transition to the dysbiotic state. Host-symbiont recognition, abiotic-stress response, micronutrient acquisition, and an antiviral defense arsenal comprising multiple restriction endonucleases, CRISPR/Cas systems, and phage lysogenization regulators are signatures of prokaryotic communities in octocorals. We argue that these features collectively contribute to the stabilization of symbiosis in the octocoral holobiont and constitute beneficial traits that can guide future studies on coral reef conservation and microbiome therapy.

59 BASIC BIOLOGICAL SCIENCES↗

ClimSim: A large multi-scale dataset for hybrid physics-ML climate emulation

Modern climate projections lack adequate spatial and temporal resolution due to computational constraints. A consequence is inaccurate and imprecise predictions of critical processes such as storms. Hybrid methods that combine physics with machine learning (ML) have introduced a new generation of higher fidelity climate simulators that can sidestep Moore’s Law by outsourcing compute-hungry, short, high-resolution simulations to ML emulators. However, this hybrid ML-physics simulation approach requires domain-specific treatment and has been inaccessible to MLexperts because of lack of training data and relevant, easy-to-use workflows. Wepresent ClimSim, the largest-ever dataset designed for hybrid ML-physics research. It comprises multi-scale climate simulations, developed by a consortium of climate scientists and ML researchers. It consists of 5.7 billion pairs of multivariate input and output vectors that isolate the influence of locally-nested, high-resolution, high-fidelity physics on a host climate simulator’s macro-scale physical state. The dataset is global in coverage, spans multiple years at high sampling frequency, and is designed such that resulting emulators are compatible with downstream coupling into operational climate simulators. We implement a range of deterministic and stochastic regression baselines to highlight the ML challenges and their scoring. The data (https://huggingface.co/datasets/LEAP/ClimSim_high-res2) and code(https://leap-stc.github.io/ClimSim)arereleasedopenlytosupport the development of hybrid ML-physics and high-fidelity climate simulations for the benefit of science and society.

artificial intelligence, machine learning↗

Search for pair production of third-generation leptoquarks decaying into a bottom quark and a $\tau $-lepton with the ATLAS detector

Abstract A search for pair-produced scalar or vector leptoquarks decaying into a b -quark and a $$\tau $$ τ -lepton is presented using the full LHC Run 2 (2015–2018) data sample of 139 fb $$^{-1}$$ - 1 collected with the ATLAS detector in proton–proton collisions at a centre-of-mass energy of $$\sqrt{s} =13$$ s = 13 TeV. Events in which at least one $$\tau $$ τ -lepton decays hadronically are considered, and multivariate discriminants are used to extract the signals. No significant deviations from the Standard Model expectation are observed and 95% confidence-level upper limits on the production cross-section are derived as a function of leptoquark mass and branching ratio $$\mathcal {B}$$ B into a $$\tau $$ τ -lepton and b -quark. For scalar leptoquarks, masses below 1460 GeV are excluded assuming $$\mathcal {B}=100$$ B = 100 %, while for vector leptoquarks the corresponding limit is 1650 GeV (1910 GeV) in the minimal-coupling (Yang–Mills) scenario.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Systematic differences between the field and cluster elliptical galaxies

Multivariate statistical techniques and fundamental plane fits are used here to study possible systematic differences between field ellipticals (FEs) and cluster ellipticals (CEs). The FEs show more intrinsic scatter in their properties, especially when stellar population variables are included. Pairwise correlations for the two samples are different; the correlations are systematically better for the cluster sample, meaning that ellipticals in the two samples populate their fundamental planes in different ways. Bivariate correlations are different for the two samples, implying that they have different fundamental planes. This is especially true for the correlations which include the population variables Mg2 and (B-V), which are sensitive both to the enrichment history and the storm formation history.

De Carvalho, R. R.↗

Self-reported health impacts of do-it-yourself air cleaner use in a smoke-impacted community

Smoke exposure from wildfires or residential wood burning for heat is a public health problem for many communities. Do-It-Yourself (DIY) portable air cleaners (PACs) are promoted as affordable alternatives to commercial PACs, but evidence of their effect on health outcomes is limited. Pilot test an evaluation of the effect of DIY PAC usage on self-reported symptoms, and investigate barriers and facilitators of PAC use, among members of a tribal community that routinely experiences elevated concentrations of fine particulate matter (PM 2.5 ) from smoke. We conducted studies in Fall 2021 (“wildfire study”; N = 10) and Winter 2022 (“wood stove study”; N = 17). Each study included four sequential one-to-two-week phases: 1) initial, 2) DIY PAC usage ≥8 h/day, 3) commercial PAC usage ≥8 h/day, and 4) air sensor with visual display and optional PAC use. We continuously monitored PAC usage and indoor/outdoor PM 2.5 concentrations in homes. Concluding each phase, we conducted phone surveys about participants’ symptoms, perceptions, and behaviors. We analyzed symptoms associated with PAC usage and conducted an analysis of indoor PM 2.5 concentrations as a mediating pathway using mixed effects multivariate linear regression. We categorized perceptions related to PACs into barriers and facilitators of use. No association was observed between PAC usage and symptoms, and the mediation analysis did not indicate that small observed trends were attributable to changes in indoor PM 2.5 concentrations. Small sample sizes hindered the ability to draw conclusions regarding the presence or absence of causal associations. DIY PAC usage was low; loud operating noise was a barrier to use. This research is novel in studying health effects of DIY PACs during wildfire and wood smoke exposures. Such research is needed to inform public health guidance. Recommendations for future studies on PAC use during smoke exposure include building flexibility of intervention timing into the study design.

54 ENVIRONMENTAL SCIENCES↗

New methods for trace analysis of gamma-irradiated pentaerythritol tetranitrate

High explosives (HEs) are used in a diverse range of applications in which they could be exposed to various radiation levels that may cause potential chemical changes. This study further evaluated pentaerythritol tetranitrate (PETN) that was previously aged with a low-level 2 kGy dose of gamma irradiation in order to understand chemical changes caused by irradiation. Both unirradiated PETN and gamma-irradiated PETN were analyzed using ultra high-pressure liquid chromatography coupled to quadrupole time of flight mass spectrometry (UHPLC-QTOF). The resulting data were processed in a non-targeted manner using Fisher's ratio analysis and multivariate curve resolution-alternating least squares (MCR-ALS) to aid in discovery and identification of the chemical changes brought about by irradiation without a priori knowledge. The application of using UHPLC-QTOF in combination with chemometric techniques for the analysis of irradiated samples has not previously been performed. Here, in this work, we show how to use this method to provide chemical information that would otherwise not be discernible, such as the discovery of degradation of the various homologues of PETN. Major differences identified with radiolytic aging of the PETN sample included decomposition products that resulted from the degradation of the trigger linkage – the O–NO 2 bonds – resulting in the formation of alcohol and aldehyde groups. Similar degradation was also observed in the PETN homologues as well as interconversion from one homologue to another.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

Application of a two-stream radiative transfer model for leaf lignin and cellulose concentrations from spectral reflectance measurements, part 1

Lignin and nitrogen contents of leaves constitute the primary rate-limiting parameters for the decomposition of forest litter, and are determinants of nutrient- and carbon-cyclic rates in forest ecosystems (Melillo et al., 1982). Wessman et al. (1988a) developed empirical multivariate relationships between forest canopy lignin amount and the (first-difference) AIS spectral response in three bands spread over the wavelength interval 1256-1555 nm. Wessman et al. (1988b) and McLellan et al. (1991) developed similar regression relationships from laboratory reflectance measurements on dried samples prepared in a standard fashion. They used four to six infrared bands for analysis of nitrogen, lignin and cellulose content of foliage in forest and prairie species. In the present article (Parts 1 and 2) the feasibility of compositional determinations is explored using positions of composite absorption bands that originate from mixtures of lignin, cellulose, and possibly other chemical constituents in the spectral reflectance of green leaves. To carry out this program, we employ full-spectral-resolution single-leaf diffuse reflectance measurements made with a laboratory spectrometer and integrating sphere. The leaf and other chemical reflectance data compiled by Elvidge (1990) have also been utilized extensively.

Conel, James E.↗

System and Method for Outlier Detection via Estimating Clusters

An efficient method and system for real-time or offline analysis of multivariate sensor data for use in anomaly detection, fault detection, and system health monitoring is provided. Models automatically derived from training data, typically nominal system data acquired from sensors in normally operating conditions or from detailed simulations, are used to identify unusual, out of family data samples (outliers) that indicate possible system failure or degradation. Outliers are determined through analyzing a degree of deviation of current system behavior from the models formed from the nominal system data. The deviation of current system behavior is presented as an easy to interpret numerical score along with a measure of the relative contribution of each system parameter to any off-nominal deviation. The techniques described herein may also be used to "clean" the training data.

Iverson, David J.↗

Does the leaf economic spectrum hold within plant functional types? A Bayesian multivariate trait meta-analysis

The leaf economic spectrum is a widely-studied axis of plant trait variability that defines a trade-off between leaf longevity and productivity. While this has been investigated at the global scale, where it is robust, and at local scales, where deviations from it are common, it has received less attention at the intermediate scale of PFTs. We investigated whether global leaf economic relationships are also present within the scale of plant functional types (PFTs) commonly used by Earth System models, and the extent to which this global-PFT hierarchy can be used to constrain trait estimates. We developed a hierarchical multivariate Bayesian model that assumes separate means and covariance structures within and across PFTs and fit this model to seven leaf traits from the TRY database related to leaf longevity, morphology, biochemistry, and photosynthetic metabolism. Although patterns of trait covariation were generally consistent with the leaf economic spectrum, we found three approximate tiers to this consistency. Relationships among morphological and biochemical traits (SLA, N, P) were the most robust within and across PFTs, suggesting that covariation in these traits is driven by universal leaf construction trade-offs and stoichiometry. Relationships among metabolic traits (R_d, V_{c,max}, J_max) were slightly less consistent, reflecting in part their much sparser sampling (especially for high-latitude PFTs), but also pointing to more flexible plasticity in plant metabolistm. Finally, relationships involving leaf lifespan were the least consistent, indicating that leaf economic relationships related to leaf lifespan are dominated by across-PFT differences and that within-PFT variation in leaf lifespan is more complex and idiosyncratic. Across all traits, these covariance were an important source of information, as evidenced by the improved imputation accuracy and reduced predictive uncertainty in multivariate models compared to univariate models. Ultimately, our study reaffirms the value of studying not just individual traits but the multivariate trait space and the utility of hierarchical modeling for studying the scale dependence of trait relationships.

biodiversity, ecological modeling, ecology, leaf b↗

Remote quantification of Cm(III) and HNO 3 by fluorescence spectroscopy and chemometrics

A unique approach to remotely quantify Cm(III) (0–100 µg mL −1 ) in HNO 3 (1–12 M) using steady-state laser fluorescence spectroscopy and multivariate regression models was developed. Photoluminescence is amenable to remote measurements using fiber-optic cables and is sensitive to numerous lanthanide and actinide species. In-line measurements can provide feedback to support complex processing in harsh environments (e.g., hot cells) to help guide and optimize radiochemical separations. In this work, Cm(III) spectra were acquired remotely in a glove box as a function of HNO 3 concentration to better understand spectral characteristics and evaluate the utility of multivariate regression models in this system. Furthermore, the Cm(III) fluorescence peak shape, width, position, and intensity changed significantly as a function of HNO 3 concentration, likely because of the displacement of emission quenching inner-sphere water molecules and complexation with nitrate ions. Despite significant covariance and nonlinearity in the data, a D-optimal design strategy successfully minimized training set sample size and was used to build effective partial least squares regression models for Cm(III) and HNO 3 concentrations without a priori knowledge of solution conditions. Chemometrics for modeling complex fluorescence spectra are promising and may find widespread applicability for online analysis in numerous chemical systems found in the nuclear field.

Actinide↗

NIR Dataset and Models (Near-Infrared Spectroscopy Calibration Dataset for Butanediol Fermentation and Multivariate Calibration Models) [SWR-22-60]

2,3-butanediol (2,3-BDO) is an economically important platform chemical that can be used in a variety of chemical feedstocks, liquid fuels, and biosynthetic building blocks. While 2,3-BDO can be efficiently produced by fermentation, the fermentation requires continuous monitoring and control to maximize 2,3-BDO yields and minimize inhibitory coproducts. Because of the time required for sampling and at-line measurement of fermentation samples with high pressure liquid chromatography (HPLC), the ability for operators to perform real-time modification to fermentation conditions is limited. To overcome this challenge, researchers from the National Renewable Energy Laboratory (NREL) have developed a calibration model which can predict the concentration of several analytes in real-time using near-infrared (NIR) spectra of the filtered fermentation broth. While significantly reducing the need for off-line sampling, NREL expects this technology to play a critical role in maximizing 2,3-BDO production. 2,3-butanediol (2,3-BDO) is a useful chemical platform that can be used to create a variety of products. For instance, 2,3-BDO can be (1) dehydrated and converted into methyl ethyl ketone, a liquid fuel additive or (2) deoxydehydrated into 1,3-butadiene for synthetic rubber, which can also be oligomerized in high yields to gasoline, diesel, and jet fuel. In order to maximize 2,3-BDO production, frequent measurement of fermentation samples is needed, as small changes in oxygen concentration can drive the fermentation to undesired products. For example, oxygen-deficient conditions result in glycerol production, while excess oxygen concentrations result in acetoin production. This results in the need for measuring dissolved oxygen, glucose, and xylose concentrations in order to optimize the aeration rate of the fermentation. Traditional monitoring methods occurs off-line and can take up to 30 minutes per sample. With multiple fermenters and high-pressure liquid chromatography (HPLC) injectors, resulting in the need for multiple samples, the sampling process can take hours to complete. NREL’s calibration model can predict the glucose, xylose, 2,3-BDO, acetoin, and glycerol concentrations from NIR spectra of filtered fermentation liquor samples. Using a partial least-squares (PLS) calibration model, NREL’s model can monitor the concentration of these analytes during subsequent fermentations at bench- and pilot-scale, demonstrating the utility of NIR spectroscopy combined with chemometrics for real-time, at-line monitoring of 2,3-BDO fermentations.

Wolfrum, Edward↗

Red Dragon: a redshift-evolving Gaussian mixture model for galaxies

ABSTRACT Precision-era optical cluster cosmology calls for a precise definition of the red sequence (RS), consistent across redshift. To this end, we present the Red Dragon algorithm: an error-corrected multivariate Gaussian mixture model (GMM). Simultaneous use of multiple colours and smooth evolution of GMM parameters result in a continuous RS and blue cloud (BC) characterization across redshift, avoiding the discontinuities of red fraction inherent in swapping RS selection colours. Based on a mid-redshift spectroscopic sample of SDSS galaxies, an RS defined by Red Dragon selects quiescent galaxies (low specific star formation rate) with a balanced accuracy of over $90{{\ \rm per\ cent}}$. This approach to galaxy population assignment gives more natural separations between RS and BC galaxies than hard cuts in colour–magnitude or colour–colour spaces. The Red Dragon algorithm is publicly available at bitbucket.org/wkblack/red-dragon-gamma/.

79 ASTRONOMY AND ASTROPHYSICS↗

A new method for detecting abrupt changes in the dependence among multivariate hydrological series based on moving cut total correlation

Knowledge of how to define and estimate the dependence among multivariate hydrological series is essential for detecting abrupt changes in the dependence. Here, in this paper, a new method (BMCTC) is proposed to detect all possible abrupt change points in the dependence among multivariate hydrological series. The total correlation estimated by the matrix-based Renyi's alpha-order entropy functional is firstly introduced to define and measure the dependence strength among multivariate hydrological series. Then, the moving cut total correlation (MCTC) sequence is built by the moving window technique, which is used to measure changes in the dependence strength among multivariate hydrological series. Finally, the Bernaola-Galvan algorithm is used to detect all change points of the MCTC sequence. Simulations are performed to compare the effectiveness of BMCTC with Pearson correlation (BMCPC) and Spearman correlation (BMCSC), Cramer-von Mises (CvM) and copula-based likelihood-ratio (CLR). The results show that all change points are detected by BMCTC regardless of the samples size, but wrong change points or no change points are detected by other methods in most cases. BMCTC is applied to detect change points in the dependence among annual runoff, precipitation and sediment discharge series in the Xiliugou and the Kuyehe River, China. It is found that the dependence among runoff, precipitation and sediment discharge changed abruptly in 1980 and 1996 in the Kuyehe River and in 1999 in the Xiliugou River. These changes are mainly caused by human activities such as construction of water conservancy projects and coal mining.

54 ENVIRONMENTAL SCIENCES↗

Analysis models for the estimation of oceanic fields

A general model for statistically optimal estimates is presented for dealing with scalar, vector and multivariate datasets. The method deals with anisotropic fields and treats space and time dependence equivalently. Problems addressed include the analysis, or the production of synoptic time series of regularly gridded fields from irregular and gappy datasets, and the estimate of fields by compositing observations from several different instruments and sampling schemes. Technical issues are discussed, including the convergence of statistical estimates, the choice of representation of the correlations, the influential domain of an observation, and the efficiency of numerical computations.

Carter, E. F.↗