Gender disparities in utilization, unmet need and expenditure for mental health care

Objectives: To compare mental health care utilization regarding the source, types, and intensity of mental health services received, unmet need for services, and out of pocket cost among non-institutionalized psychologically distressed women and men. ^ Method: Cross-sectional data for 19,325 non-institutionalized mentally distressed adult respondents to the “The National Survey on Drug Use and Health” (NSDUH), for the years 2006 -2008, representing over twenty-nine millions U.S. adults was analyzed. To assess the relative odds for women compared to men, logistic regression analysis was used for source of service, for types of barriers, for unmet need and cost; zero inflated negative binomial regression for intensity of utilization; and ordinal logistic regression analysis for quantifying out-of-pocket expenditure. ^ Results: Overall, 43% of mentally distressed adults utilized a form of mental health treatment; representing 12.6 million U.S psychologically distressed adults. Females utilized more mental health care compared to males in the previous 12 months (OR: 1. 70; 95% CI: 1.54, 1.83). Similarly, females were 54% more likely to get help for psychological distress in an outpatient setting and females were associated with an increased probability of using medication for mental distress (OR: 1.72; 95% CI: 1.63, 1.98). Women were 1.25 times likelier to visit a mental health center (specialty care) than men. ^ Females were positively associated with unmet needs (OR: 1.50; 95% CI: 1.29, 1.75) after taking into account predisposing, enabling, and need (PEN) characteristics. Women with perceived unmet needs were 23% (OR: 0.77; 95% CI: 0.59, 0.99) less likely than men to report societal accommodation (stigma) as a barrier to mental health care. At any given cutoff point, women were 1.74 times likelier to be in the higher payment categories for inpatient out of pocket cost when other variables in the model are held constant. Conclusions: Women utilize more specialty mental healthcare, report more unmet need, and pay more inpatient out of pocket costs than men. These gender disparities exist even after controlling for predisposing, enabling, and need variables. Creating policies that not only provide mental health care access but also de-stigmatize mental illness will bring us one step closer to eliminating gender disparities in mental health care.^

Evaluation of PRESAGE® dosimeters for brachytherapy sources and the 3D dosimetry and characterization of the new AgX100 125I seed model

With continuous new improvements in brachytherapy source designs and techniques, method of 3D dosimetry for treatment dose verifications would better ensure accurate patient radiotherapy treatment. This study was aimed to first evaluate the 3D dose distributions of the low-dose rate (LDR) Amersham 6711 OncoseedTM using PRESAGE® dosimeters to establish PRESAGE® as a suitable brachytherapy dosimeter. The new AgX100 125I seed model (Theragenics Corporation) was then characterized using PRESAGE® following the TG-43 protocol. PRESAGE® dosimeters are solid, polyurethane-based, 3D dosimeters doped with radiochromic leuco dyes that produce a linear optical density response to radiation dose. For this project, the radiochromic response in PRESAGE® was captured using optical-CT scanning (632 nm) and the final 3D dose matrix was reconstructed using the MATLAB software. An Amersham 6711 seed with an air-kerma strength of approximately 9 U was used to irradiate two dosimeters to 2 Gy and 11 Gy at 1 cm to evaluate dose rates in the r=1 cm to r=5 cm region. The dosimetry parameters were compared to the values published in the updated AAPM Report No. 51 (TG-43U1). An AgX100 seed with an air-kerma strength of about 6 U was used to irradiate two dosimeters to 3.6 Gy and 12.5 Gy at 1 cm. The dosimetry parameters for the AgX100 were compared to the values measured from previous Monte-Carlo and experimental studies. In general, the measured dose rate constant, anisotropy function, and radial dose function for the Amersham 6711 showed agreements better than 5% compared to consensus values in the r=1 to r=3 cm region. The dose rates and radial dose functions measured for the AgX100 agreed with the MCNPX and TLD-measured values within 3% in the r=1 to r=3 cm region. The measured anisotropy function in PRESAGE® showed relative differences of up to 9% with the MCNPX calculated values. It was determined that post-irradiation optical density change over several days was non-linear in different dose regions, and therefore the dose values in the r=4 to r=5 cm regions had higher uncertainty due to this effect. This study demonstrated that within the radial distance of 3 cm, brachytherapy dosimetry in PRESAGE® can be accurate within 5% as long as irradiation times are within 48 hours.

Robust effect sizes and their confidence intervals for group difference between trajectories in hierarchical linear growth model

Hierarchical linear growth model (HLGM), as a flexible and powerful analytic method, has played an increased important role in psychology, public health and medical sciences in recent decades. Mostly, researchers who conduct HLGM are interested in the treatment effect on individual trajectories, which can be indicated by the cross-level interaction effects. However, the statistical hypothesis test for the effect of cross-level interaction in HLGM only show us whether there is a significant group difference in the average rate of change, rate of acceleration or higher polynomial effect; it fails to convey information about the magnitude of the difference between the group trajectories at specific time point. Thus, reporting and interpreting effect sizes have been increased emphases in HLGM in recent years, due to the limitations and increased criticisms for statistical hypothesis testing. However, most researchers fail to report these model-implied effect sizes for group trajectories comparison and their corresponding confidence intervals in HLGM analysis, since lack of appropriate and standard functions to estimate effect sizes associated with the model-implied difference between grouping trajectories in HLGM, and also lack of computing packages in the popular statistical software to automatically calculate them. ^ The present project is the first to establish the appropriate computing functions to assess the standard difference between grouping trajectories in HLGM. We proposed the two functions to estimate effect sizes on model-based grouping trajectories difference at specific time, we also suggested the robust effect sizes to reduce the bias of estimated effect sizes. Then, we applied the proposed functions to estimate the population effect sizes (d ) and robust effect sizes (du) on the cross-level interaction in HLGM by using the three simulated datasets, and also we compared the three methods of constructing confidence intervals around d and du recommended the best one for application. At the end, we constructed 95% confidence intervals with the suitable method for the effect sizes what we obtained with the three simulated datasets. ^ The effect sizes between grouping trajectories for the three simulated longitudinal datasets indicated that even though the statistical hypothesis test shows no significant difference between grouping trajectories, effect sizes between these grouping trajectories can still be large at some time points. Therefore, effect sizes between grouping trajectories in HLGM analysis provide us additional and meaningful information to assess group effect on individual trajectories. In addition, we also compared the three methods to construct 95% confident intervals around corresponding effect sizes in this project, which handled with the uncertainty of effect sizes to population parameter. We suggested the noncentral t-distribution based method when the assumptions held, and the bootstrap bias-corrected and accelerated method when the assumptions are not met.^

Chemical characterization and quantification of the brown algal storage compound laminarin

In search of a meaningful stress indicator for Fucus vesiculosus we found that the often used quantitative determination procedures for the polysaccharide laminarin (beta-1,3-glucan) result in different kind of problems, uncertainties and limitations. This chemical long-term storage form of carbon enables perennial brown algae in seasonally fluctuating ecosystems to uncouple growth from photosynthesis. Because of this high ecological relevance a reliable and precise method for determination and quantification of laminarin is needed. Therefore, a simple, cold water extraction method coupled to a new quantitative liquid chromatography-mass spectrometrical method (LC-MS) was developed. Laminarin was determined in nine out of twelve brown algal species, and its expected typical molar mass distribution of 2000-7000 Da was confirmed. Furthermore, laminarin consisted of a complex mixture of different chemical forms, since fifteen chemical laminarin species with distinct molecular weights were measured in nine species of brown algae. Laminarin concentrations in the algal tissues ranged from 0.03 to 0.86% dry weight (DW). The direct chemical characterization and quantification of laminarin by LC-MS represents a powerful method to verify the biochemical and ecological importance of laminarin for brown algae. Single individuals of Laminaria hyperborea, L. digitata, Saccharina latissima, F. serratus, F. vesiculosus, F. spiralis, Himanthalia elongata, Cystoseira tamariscifolia, Pelvetia canaliculata, Ascophyllum nodosum, Halidrys siliquosa and Dictyota dichotoma were collected in fall (18.11.2013) during spring low tide from the shore of Finavarra, Co. Clare, west coast of Ireland (53° 09' 25'' N, 09° 06' 58'' W). After sampling, the different algae were immediately transported to the lab, lyophilized and sent to the University of Rostock. Laminarin was extracted with cold ultrapure water from the algal samples. Before extraction they were ground to < 1 mm grain size with an analytical mill (Ika MF 10 Basic). The algal material (approx. 1.5 g DW) was extracted in ultrapure water (8 mL) on a shaker (250 rpm) for 5 h. After the addition of surplus ultrapure water (4 mL) and shaking manually, 1 mL of the sample was filter centrifuged (45 µm) at 14,000 rpm (Hettich Mikro 22 R). The slightly viscous supernatant was free of suspended material and converted into a microvial (300 µL) for further analysis. The extracts were analyzed using liquid chromatography-mass spectrometry (LC-MS) analysis (LTQ Velos Pro ion trap spectrometer with Accela HPLC, Thermo Scientific). Laminarin species were separated on a KinetexTM column (2.6 µm C18, 150 x 3 mm). The mobile phase was 90 % ultrapure water and 10 % acetonitrile, run isocratically at a flow rate of 0.2 mL min-1. MS was working in ESI negative ion mode in a mass range of 100 - 4000 amu. Glucose contents were determined after extraction using high-performance liquid chromatography (HPLC). Extracted samples were analyzed in an HPLC (SmartLine, Knauer GmbH) equipped with a SUPELCOGELTM Ca column (30 x 7,8 mm without preColumn) and RI-detector (S2300 PDA S2800). Water was used as eluent at a flow rate of 0.8 mL min-1 at 75 °C. Glucose was quantified by comparison of the retention time and peak area with standard solutions using ChromGate software. Mannitol was extracted from three subsamples of 10-20 mg powdered alga material (L. hyperborea, L. digitata, S. latissima, F. serratus, F. vesiculosus, F. spiralis, H. elongata, P. canaliculata, A. nodosum, H. siliquosa) and quantified, following the HPLC method described by Karsten et al. (1991). For analyzing carbon and nitrogen contents, dried algal material was ground to powder and three subsamples of 2 mg from each alga thalli were loaded and packed into tin cartridges (6×6×12 mm). The packages were combusted at 950 °C and the absolute contents of C and N were automatically quantified in an elemental analyzer (Elementar Vario EL III, Germany) using acetanilide as standard according to Verardo et al. (1990).

Effect of colchicine and amiprophos-methyl on the production of in vitro doubled haploid onion plants and correlation assessment between ploidy level and stomatal size

Doubled haploid onion (Allium cepa L.) plants allow the production of completely homozygous lines for a later production of hybrids. The haploid plants are normally produced using in vitro gynogenesis. The obtained haploid plantlets must be treated with different agents for doubling chromosomes. It is necessary to adjust the concentration and the length of treatment of the doubling agent. In this case, the effect of 250 and 500 mg.L-1 colchicine and 15.2; 30 and 60 mg.L- 1 amiprophos-methyl during 24 and 48 h was assessed over the rate of onion haploid plantlets chromosome doubling. The best duplication treatment was 250 mg.L-1 colchicine for 48 h, which yielded 100% of doubled haploid plants. On the other hand, a positive correlation resulted from the ploidy level and stomatal size, and a negative correlation between the level of ploidy and stomatal density. Significant differences between the stomatal length, width and density in haploid and doubled haploid plantlets were observed. An economical and quick method to test ploidy level in onion plantlets is proposed through the measurement of stomatal size and density.

Depths and ages of seismic sequence boundaries in ODP Leg 166 holes

High-resolution, multichannel seismic data collected across the Great Bahama Bank margin and the adjacent Straits of Florida indicate that the deposition of Neogene-Quaternary strata in this transect are controlled by two sedimentation mechanisms: (1) west-dipping layers of the platform margin, which are a product of sea-level-controlled, platform-derived downslope sedimentation; and (2) east- or north-dipping drift deposits in the basinal areas, which are deposited by ocean currents. These two sediment systems are active simultaneously and interfinger at the toe-of-slope. The prograding system consists of sigmoidal clinoforms that advanced the margin some 25 km into the Straits of Florida. The foresets of the clinoforms are approximately 600 m high with variable slope angles that steepen significantly in the Pleistocene section. The seismic facies of the prograding clinoforms on the slope is characterized by dominant, partly chaotic, cut-and-fill geometries caused by submarine canyons that are oriented downslope. In the basin axis, seismic geometries and facies document deposition from and by currents. Most impressive is an 800-m-thick drift deposit at the confluence of the Santaren Channel and the Straits of Florida. This "Santaren Drift" is slightly asymmetric, thinning to the north. The drift displays a highly coherent seismic facies characterized by a continuous succession of reflections, indicating very regular sedimentation. Leg 166 of the Ocean Drilling Program (ODP) drilled a transect of five deep holes between 2 and 30 km from the modern platform margin and retrieved the sediments from both the slope and basin systems. The Neogene slope sediments consist of peri-platform oozes intercalated with turbidites, whereas the basinal drift deposits consist of more homogeneous, fine-grained carbonates that were deposited without major hiatuses by the Florida Current starting at approximately 12.4 Ma. Sea-level fluctuations, which controlled the carbonate production on Great Bahama Bank by repeated exposure of the platform top, controlled lithologic alternations and hiatuses in sedimentation across the transect. Both sedimentary systems are contained in 17 seismic sequences that were identified in the Neogene-Quaternary section. Seismic sequence boundaries were identified based on geometric unconformities beneath the Great Bahama Bank. All the sequence boundaries could be traced across the entire transect into the Straits of Florida. Biostratigraphic age determinations of seismic reflections indicate that the seismic reflections of sequence boundaries have chronostratigraphic significance across both depositional environments.

Chemical element concentrations in sediments from DSDP Holes 60-458 and 60-459

The development of laser ablation-inductively coupled plasma-mass spectrometry has revolutionized the analysis of tephras by providing (1) an efficient and precise method for determining abundances of a wide variety of trace elements at low concentrations in individual glass shards and (2) assessment of geochemical heterogeneities within individual ash horizons. This development is important for petrogenetic studies of intraoceanic arc systems, where tephras provide the most complete temporal record of magmatism. Results from the Izu-Bonin and Mariana arc systems indicate that despite close geographical proximity and similar tectonic evolution, they contrast strongly in terms of geochemical evolution since 35 Ma. Whereas the Mariana tephras have exceptional compositional diversity, ranging from low-K (Oligocene), to high-K (Miocene), and subsequently medium-K compositions (Pliocene-Quaternary), the Izu-Bonin arc has been dominated by low-K compositions throughout. The Mariana increases in K are paralleled by increases in abundances of incompatible trace elements and by increased values of diagnostic ratios (e.g., Nb/yb and Th/yb) regarded as monitors of potential mantle-source fertility. The relative uniformity of Nb/yb and Nb/Zr ratios in Izu-Bonin tephras indicates that cyclic processes of backarc basin development and mantle depletion do not necessarily induce large-scale temporal geochemical variations in the associated arc. Temporal variability within the Mariana arc, and its divergence from the Izu-Bonin arc ca. 13 Ma, can be traced to a major injection of subducted sediment in the Mariana system at this time.

Phosphate d18O pore-water profiles of sediment core GeoB11807-2 and GeoB11804-4

Phosphorus cycling in the ocean is influenced by biological and geochemical processes that are reflected in the oxygen isotope signature of dissolved inorganic phosphate (Pi). Extending the Pi oxygen isotope record from the water column into the seabed is difficult due to low Pi concentrations and small amounts of marine porewaters available for analysis. We obtained porewater profiles of Pi oxygen isotopes using a refined protocol based on the original micro-extraction designed by Colman (2002). This refined and customized method allows the conversion of ultra-low quantities (0.5 - 1 µmol) of porewater Pi to silver phosphate (Ag3PO4) for routine analysis by mass spectrometry. A combination of magnesium hydroxide co-precipitation with ion exchange resin treatment steps is used to remove dissolved organic matter, anions, and cations from the sample before precipitating Ag3PO4. Samples as low as 200 µg were analyzed in a continuous flow isotope ratio mass spectrometer setup. Tests with external and laboratory internal standards validated the preservation of the original phosphate oxygen isotope signature (d18OP) during micro extraction. Porewater data on d18OP has been obtained from two sediment cores of the Moroccan margin. The d18OP values are in a range of +19.49 to +27.30 per mill. We apply a simple isotope mass balance model to disentangle processes contributing to benthic P cycling and find evidence for Pi regeneration outbalancing microbial demand in the upper sediment layers. This highlights the great potential of using d18OP to study microbial processes in the subseafloor and at the sediment water interface.

Age models and summer sea surface temperature and winter sea ice concentration for the EPILOG-LGM time slice in the Pacific Southern Ocean

Sea surface temperatures and sea-ice extent are the most critical variables to evaluate the Southern Ocean paleoceanographic evolution in relation to the development of the global carbon cycle, atmospheric CO2 variability and ocean-atmosphere circulation. In contrast to the Atlantic and the Indian sectors, the Pacific sector of the Southern Ocean has been insufficiently investigated so far. To cover this gap of information we present diatom-based estimates of summer sea surface temperature (SSST) and winter sea-ice concentration (WSI) from 17 sites in the polar South Pacific to study the Last Glacial Maximum (LGM) at the EPILOG time slice (19,000-23,000 cal. years BP). Applied statistical methods are the Imbrie and Kipp Method (IKM) and the Modern Analog Technique (MAT) to estimate temperature and sea-ice concentration, respectively. Our data display a distinct LGM east-west differentiation in SSST and WSI with steeper latitudinal temperature gradients and a winter sea-ice edge located consistently north of the Pacific-Antarctic Ridge in the Ross sea sector. In the eastern sector of our study area, which is governed by the Amundsen Abyssal Plain, the estimates yield weaker latitudinal SSST gradients together with a variable extended winter sea-ice field. In this sector, sea-ice extent may have reached sporadically the area of the present Subantarctic Front at its maximum LGM expansion. This pattern points to topographic forcing as major controller of the frontal system location and sea-ice extent in the western Pacific sector whereas atmospheric conditions like the Southern Annular Mode and the ENSO affected the oceanographic conditions in the eastern Pacific sector. Although it is difficult to depict the location and the physical nature of frontal systems separating the glacial Southern Ocean water masses into different zones, we found a distinct temperature gradient in latitudes straddled by the modern Southern Subtropical Front. Considering that the glacial temperatures north of this zone are similar to the modern, we suggest that this represents the Glacial Southern Subtropical Front (GSSTF), which delimits the zone of strongest glacial SSST cooling (>4K) to its North. The southern boundary of the zone of maximum cooling is close to the glacial 4°C isotherm. This isotherm, which is in the range of SSST at the modern Antarctic Polar Front (APF), represents a circum-Antarctic feature and marks the northern edge of the glacial Antarctic Circumpolar Current (ACC). We also assume that a glacial front was established at the northern average winter sea ice edge, comparable with the modern Southern Antarctic Circumpolar Current Front (SACCF). During the glacial, this front would be located in the area of the modern APF. The northward deflection of colder than modern surface waters along the South American continent leads to a significant cooling of the glacial Humboldt Current surface waters (4-8K), which affects the temperature regimes as far north as into tropical latitudes. The glacial reduction of ACC temperatures may also result in the significant cooling in the Atlantic and Indian Southern Ocean, thus may enhance thermal differentiation of the Southern Ocean and Antarctic continental cooling. Comparison with temperature and sea ice simulations for the last glacial based on numerical simulations show that the majority of modern models overestimate summer and winter sea ice cover and that there exists few models that reproduce our temperature data rather well.

Seawater carbonate chemistry and pteropod Limacina helicina antarctica shell dissolution during experiments, 2012

Anthropogenic ocean acidification is likely to have negative effects on marine calcifying organisms, such as shelled pteropods, by promoting dissolution of aragonite shells. Study of shell dissolution requires an accurate and sensitive method for assessing shell damage. Shell dissolution was induced through incubations in CO2 enriched seawater for between 4 and 14 days. We describe a procedure that allows the level of dissolution to be assessed and classified into three main types: Type I with partial dissolution of the prismatic layer; Type II with exposure of underlying crossed-lamellar layer, and Type III, where crossed-lamellar layer shows signs of dissolution. Levels of dissolution showed a good correspondence to the incubation conditions, with the most severe damage found in specimens held for 14 d in undersaturated condition (Ohm ~ 0.8). This methodology enables the response of small pelagic calcifiers to acidified conditions to be detected at an early stage, thus making pteropods a valuable bioindicator of future ocean acidification.

Video analysis based vehicle detection and tracking using an MCMC sampling framework

This article presents a probabilistic method for vehicle detection and tracking through the analysis of monocular images obtained from a vehicle-mounted camera. The method is designed to address the main shortcomings of traditional particle filtering approaches, namely Bayesian methods based on importance sampling, for use in traffic environments. These methods do not scale well when the dimensionality of the feature space grows, which creates significant limitations when tracking multiple objects. Alternatively, the proposed method is based on a Markov chain Monte Carlo (MCMC) approach, which allows efficient sampling of the feature space. The method involves important contributions in both the motion and the observation models of the tracker. Indeed, as opposed to particle filter-based tracking methods in the literature, which typically resort to observation models based on appearance or template matching, in this study a likelihood model that combines appearance analysis with information from motion parallax is introduced. Regarding the motion model, a new interaction treatment is defined based on Markov random fields (MRF) that allows for the handling of possible inter-dependencies in vehicle trajectories. As for vehicle detection, the method relies on a supervised classification stage using support vector machines (SVM). The contribution in this field is twofold. First, a new descriptor based on the analysis of gradient orientations in concentric rectangles is dened. This descriptor involves a much smaller feature space compared to traditional descriptors, which are too costly for real-time applications. Second, a new vehicle image database is generated to train the SVM and made public. The proposed vehicle detection and tracking method is proven to outperform existing methods and to successfully handle challenging situations in the test sequences.

A combined measure for quantifying and qualifying the topology preservation of growing self-organizing maps

The Self-OrganizingMap (SOM) is a neural network model that performs an ordered projection of a high dimensional input space in a low-dimensional topological structure. The process in which such mapping is formed is defined by the SOM algorithm, which is a competitive, unsupervised and nonparametric method, since it does not make any assumption about the input data distribution. The feature maps provided by this algorithm have been successfully applied for vector quantization, clustering and high dimensional data visualization processes. However, the initialization of the network topology and the selection of the SOM training parameters are two difficult tasks caused by the unknown distribution of the input signals. A misconfiguration of these parameters can generate a feature map of low-quality, so it is necessary to have some measure of the degree of adaptation of the SOM network to the input data model. The topologypreservation is the most common concept used to implement this measure. Several qualitative and quantitative methods have been proposed for measuring the degree of SOM topologypreservation, particularly using Kohonen's model. In this work, two methods for measuring the topologypreservation of the Growing Cell Structures (GCSs) model are proposed: the topographic function and the topology preserving map

A New Method to Reduce Truncation Errors in Partial Spherical Near-Field Measurements

A new and effective method for reduction of truncation errors in partial spherical near-field (SNF) measurements is proposed. The method is useful when measuring electrically large antennas, where the measurement time with the classical SNF technique is prohibitively long and an acquisition over the whole spherical surface is not practical. Therefore, to reduce the data acquisition time, partial sphere measurement is usually made, taking samples over a portion of the spherical surface in the direction of the main beam. But in this case, the radiation pattern is not known outside the measured angular sector as well as a truncation error is present in the calculated far-field pattern within this sector. The method is based on the Gerchberg-Papoulis algorithm used to extrapolate functions and it is able to extend the valid region of the calculated far-field pattern up to the whole forward hemisphere. To verify the effectiveness of the method, several examples are presented using both simulated and measured truncated near-field data.

Development and Assessment of the Competence Creativity Applied to Technical Drawing

The results obtained after incorporating the competence “creativity” to the subject Technical Drawing of the first course of the Degree in Forestry, Technical University of Madrid, are presented in this study.At first, learning activities which could serve two functions at the same time -developing students’ creativity and developing other specific competences of the subject- were considered. Besides, changes in the assessment procedure were made and a method which analyzes two aspects of the assessment of the competence creativity was established. On the one hand, the products are evaluated by analyzing the outcomes obtained by students in the essays suggested and by establishing a parameter to assess the creativity expressed in those essays. On the other, an assessment of the student is directly carried out through a psychometric test which has been previously chosen by the team.Moreover, these results can be applied to similar or could be of general application

Estimation of the solubility parameters of model plant surfaces and agrochemicals: a valuable tool for understanding plant surface interactions

Background Most aerial plant parts are covered with a hydrophobic lipid-rich cuticle, which is the interface between the plant organs and the surrounding environment. Plant surfaces may have a high degree of hydrophobicity because of the combined effects of surface chemistry and roughness. The physical and chemical complexity of the plant cuticle limits the development of models that explain its internal structure and interactions with surface-applied agrochemicals. In this article we introduce a thermodynamic method for estimating the solubilities of model plant surface constituents and relating them to the effects of agrochemicals. Results Following the van Krevelen and Hoftyzer method, we calculated the solubility parameters of three model plant species and eight compounds that differ in hydrophobicity and polarity. In addition, intact tissues were examined by scanning electron microscopy and the surface free energy, polarity, solubility parameter and work of adhesion of each were calculated from contact angle measurements of three liquids with different polarities. By comparing the affinities between plant surface constituents and agrochemicals derived from (a) theoretical calculations and (b) contact angle measurements we were able to distinguish the physical effect of surface roughness from the effect of the chemical nature of the epicuticular waxes. A solubility parameter model for plant surfaces is proposed on the basis of an increasing gradient from the cuticular surface towards the underlying cell wall. Conclusions The procedure enabled us to predict the interactions among agrochemicals, plant surfaces, and cuticular and cell wall components, and promises to be a useful tool for improving our understanding of biological surface interactions.