Pharmacokinetics of insulin lispro in type 2 diabetes during closed-loop insulin delivery

Insulin pharmacokinetics is not well understood during continuous subcutaneous insulin infusion in type 2 diabetes (T2D). We analyzed data collected in 11 subjects with T2D [6 male, 9 white European and two of Indian ethnicity; age 59.7(12.1) years, BMI 30.1(3.9) kg/m², fasting C-peptide 1002.2(365.8) pmol/l, fasting plasma glucose 9.6(2.2) mmol/l, diabetes duration 8.0(6.2) years and HbA1c 8.3(0.8)%; mean(SD)] who underwent a 24-h study investigating closed-loop insulin delivery at the Wellcome Trust Clinical Research Facility, Cambridge, UK. Subcutaneous delivery of insulin lispro was modulated every 15 min according to a model predictive control algorithm. Two complementary insulin assays facilitated discrimination between exogenous (lispro) and endogenous plasma insulin concentrations measured every 15–60 min. Lispro pharmacokinetics was represented by a linear two-compartment model whilst parameters were estimated using a Bayesian approach applying a closed-form model solution. The time-to-peak of lispro absorption (t_max) was 109.6 (75.5–120.5) min [median (interquartile range)] and the metabolic clearance rate (MCR_I) 1.26 (0.87–1.56) × 10⁻² l/kg/min. MCR_I was negatively correlated with fasting C-peptide (r_s = −0.84; P = .001) and with fasting plasma insulin concentration (r_s = −0.79; P = .004). In conclusion, compartmental modelling adequately represents lispro kinetics during continuous subcutaneous insulin infusion in T2D. Fasting plasma C-peptide or fasting insulin may be predictive of lispro metabolic clearance rate in T2D but further investigations are warranted.     

Impact of glucocorticoids on insulin resistance in the critically ill

Glucocorticoids (GCs) have been shown to reduce insulin sensitivity in healthy individuals. Widely used in critical care to treat a variety of inflammatory and allergic disorders, they may inadvertently exacerbate stress-hyperglycaemia. This research uses model-based methods to quantify the reduction in insulin sensitivity from GCs in critically ill patients, and thus their impact on glycaemic control. A model-based measure of insulin sensitivity (S_I) was used to quantify changes between two matched cohorts of 40 intensive care unit (ICU) patients. Patients in one cohort received GC treatment, while patients in the control cohort did not. All patients were admitted to the Christchurch hospital ICU between 2005 and 2007 and spent at least 24 h on the SPRINT glycaemic control protocol.A 31% reduction in whole-cohort median insulin sensitivity was seen between the control cohort and patients receiving glucocorticoids with a median dose equivalent to 200 mg/d of hydrocortisone per patient. Comparing percentile patients as a surrogate for matched patients, reductions in median insulin sensitivity of 20%, 25%, and 21% were observed for the 25th-, 50th- and 75th-percentile patients, respectively. These cohort and percentile patient reductions are less than or equivalent to the 30-62% reductions reported in healthy subjects especially when considering the fact that the GC doses in this study are 1.3-4.0 times larger than those in studies of healthy subjects. This reduced suppression of insulin sensitivity in critically ill patients could be a result of saturation due to already increased levels of catecholamines and cortisol common in critically illness. Virtual trial simulation showed that reductions in insulin sensitivity of 20-30% associated with glucocorticoid treatment in the ICU have limited impact on glycaemic control levels within the context of the SPRINT protocol.     

Brain mass estimation by head circumference and body mass methods in neonatal glycaemic modelling and control

Introduction

Hyperglycaemia is a common complication of stress and prematurity in extremely low-birth-weight infants. Model-based insulin therapy protocols have the ability to safely improve glycaemic control for this group. Estimating non-insulin-mediated brain glucose uptake by the central nervous system in these models is typically done using population-based body weight models, which may not be ideal.

Method

A head circumference-based model that separately treats small-for-gestational-age (SGA) and appropriate-for-gestational-age (AGA) infants is compared to a body weight model in a retrospective analysis of 48 patients with a median birth weight of 750 g and median gestational age of 25 weeks. Estimated brain mass, model-based insulin sensitivity (S_I) profiles, and projected glycaemic control outcomes are investigated. SGA infants (5) are also analyzed as a separate cohort.

Results

Across the entire cohort, estimated brain mass deviated by a median 10% between models, with a per-patient median difference in S_I of 3.5%. For the SGA group, brain mass deviation was 42%, and per-patient S_I deviation 13.7%. In virtual trials, 87–93% of recommended insulin rates were equal or slightly reduced (Δ < 0.16 mU/h) under the head circumference method, while glycaemic control outcomes showed little change.

Conclusion

The results suggest that body weight methods are not as accurate as head circumference methods. Head circumference-based estimates may offer improved modelling accuracy and a small reduction in insulin administration, particularly for SGA infants.     

Standard Linux for embedded real-time robotics and manufacturing control systems

The main goal of the research presented in this paper is to evaluate the possibility of using standard Linux for embedded real-time applications in robotics and manufacturing as a consequence of dramatic improvements in hardware computing power and free software quality in the last few years. After an accurate analysis of the problems related to make Linux, a native Unix-like fair kernel, real-time, laboratory tests showed that a large variety of applications (up to 1 KHz) can be implemented using Linux and commercial-of-the-shelf hardware. Practical examples of the control systems of an unmanned surface vessel used for robotics research and of a marking machine for steelworks are reported and discussed.     

Web-based Personalised System of Instruction: An effective approach for diverse cohorts with virtual learning environments?

The Personalised System of Instruction is a form of mastery learning which, though it has been proven to be educationally effective, has never seriously challenged the dominant lecture-tutorial teaching method in higher education and has largely fallen into disuse. An information and communications technology assisted version of the Personalised System of Instruction using a virtual learning environment is promoted here based on the authors’ longitudinal design research into this pedagogy. The particular elements of the virtual learning environment which are promoted are short video clips, online formative tests and an assessment management system. The authors present their experiences of developing and deploying this pedagogy for the teaching of introductory discrete mathematics to large classes of Computer Science students at two UK higher education institutions both with whole cohorts and ‘at risk’ groups of students. In particular, this method is promoted as particularly helpful to students who do not adopt a deep approach to learning as many students fail to do. Moreover ‘at risk’ students using this method (n = 71) demonstrated an average Glass effect size of 0.83 compared with other ‘at risk’ students who did not (n = 35). Based on these experiences, this pedagogy is promoted as an effective approach to teaching in higher education, especially the teaching of cognitive skills to diverse cohorts of students on foundation level modules.     

Assessing the potential of SeaWiFS and MODIS for estimating chlorophyll concentration in turbid productive waters using red and near-infrared bands

Bio-optical algorithms for remote estimation of chlorophyll-a concentration (Chl) in case-1 waters exploit the upwelling radiation in the blue and green spectral regions. In turbid productive waters other constituents, that vary independently of Chl, absorb and scatter light in these spectral regions. As a consequence, the accurate estimation of Chl in turbid productive waters has so far not been feasible from satellite sensors. The main purpose of this study was to evaluate the extent to which near-infrared (NIR) to red reflectance ratios could be applied to the Sea Wide Field-of-View Sensor (SeaWiFS) and the Moderate Imaging Spectrometer (MODIS) to estimate Chl in productive turbid waters. To achieve this objective, remote-sensing reflectance spectra and relevant water constituents were collected in 251 stations over lakes and reservoirs with a wide variability in optical parameters (i.e. 4 ≤ Chl ≤ 240 mg m^− 3; 18 ≤ Secchi disk depth ≤ 308 cm). SeaWiFS and MODIS NIR and red reflectances were simulated by using the in-situ hyperspectral data. The proposed algorithms predicted Chl with a relative random uncertainty of approximately 28% (average bias between − 1% and − 4%). The effects of reflectance uncertainties on the predicted Chl were also analyzed. It was found that, for realistic ranges of R_rs uncertainties, Chl could be estimated with a precision better than 40% and an accuracy better than ± 35%. These findings imply that, provided that an atmospheric correction scheme specific for the red-NIR spectral region is available, the extensive database of SeaWiFS and MODIS images could be used to quantitatively monitor Chl in turbid productive waters.     

Evaluation of three methodologies to estimate the VO2max in people of different ages

Aging and gender are factors that affect the variation of physical work capacity. The present paper highlights the importance of the metabolism used by ergonomics to establish the appropriate limits of loads at work. This study compares the aerobic capacity of people from 20 to 71 years old split in 5 different groups. The laboratory experiment tested 33 volunteers (19 women and 14 men). A submaximal step test was used to measure the VO₂ using a portable breath by breath metabolic system and a telemetric heart rate monitor. Three methods to estimate the VO_2max were compared: 1) a direct measurement of VO₂, 2) estimation by heart rate, and 3) a step test method using predetermined charts. Significant difference was encountered among the estimation methods as well as among the age ranges (F_2,92 = 6.43, p < 0.05 y F_4,92 = 7.18, p < 0.05 respectively). The method of direct measurement and the method of predetermined charts were different for the estimation of the VO_2max with a confidence level of 95%. The method of predetermined charts is better adapted for males and people younger than 30 years. The estimation through non invasive heart rate apparatus was a good appraiser of the maximal oxygen consumption considering both genders and all the age groups.     

Internet testing: A natural experiment reveals test score inflation on a high-stakes,unproctored cognitive test

A natural experiment is reported in which the performance of two cohorts of undergraduate students on identical tests of cognitive achievement is compared under differing conditions of proctoring and modes of administration, total N = 302. Between-subjects and within-subjects analyses reveal clear effects of test score inflation in the unproctored conditions compared with formal, proctored, examination conditions. These results caution against the use of unproctored cognitive testing in high-stakes situations, including on-line assessment.     

Randomness quality of permuted pseudorandom binary sequences

This paper uses the DIEHARD statistical test suite to test the randomness quality of “permuted” versions of maximum length sequences generated by linear finite state machines (LFSM) such as cellular automata and linear feedback shift registers. Analysis shows that permuted sequences can be equivalently generated by using time-varying transformations derived from the original LFSM. Based on the above, we suggest the permuted transformation sequence scheme. Experimental results show that DIEHARD results are improved with respect to the original non-permuted sequences—up to seven more tests can be passed (total of 19 tests). Furthermore, a permutation vector is used to generate cyclically distinct permuted sequences and each sequence has a desirable maximum length period of 2ⁿ − 1.     

Exploring LiDAR-RaDAR synergy—predicting aboveground biomass in a southwestern ponderosa pine forest using LiDAR, SAR and InSAR

Scanning Light Detecting and Ranging (LiDAR), Synthetic Aperture Radar (SAR) and Interferometric SAR (InSAR) were analyzed to determine (1) which of the three sensor systems most accurately predicted forest biomass, and (2) if LiDAR and SAR/InSAR data sets, jointly considered, produced more accurate, precise results relative to those same data sets considered separately. LiDAR ranging measurements, VHF-SAR cross-sectional returns, and X- and P-band cross-sectional returns and interferometric ranges were regressed with ground-estimated (from dbh) forest biomass in ponderosa pine forests in the southwestern United States. All models were cross-validated. Results indicated that the average canopy height measured by the scanning LiDAR produced the best predictive equation. The simple linear LiDAR equation explained 83% of the biomass variability (n = 52 plots) with a cross-validated root mean square error of 26.0 t/ha. Additional LiDAR metrics were not significant to the model. The GeoSAR P-band (λ = 86 cm) cross-sectional return and the GeoSAR/InSAR canopy height (X-P) captured 30% of the forest biomass variation with an average predictive error of 52.5 t/ha. A second RaDAR-FOPEN collected VHF (λ ∼ 7.8 m) and cross-polarized P-band (λ = 88 cm) cross-sectional returns, none of which proved useful for forest biomass estimation (cross-validated R² = 0.09, RMSE = 63.7 t/ha). Joint consideration of LiDAR and RaDAR measurements produced a statistically significant, albeit small improvement in biomass estimation precision. The cross-validated R² increased from 83% to 84% and the prediction error decreased from 26.0 t/ha to 24.9 t/ha when the GeoSAR X-P interferometric height is considered along with the average LiDAR canopy height. Inclusion of a third LiDAR metric, the 60th decile height, further increased the R² to 85% and decreased the RMSE to 24.1 t/ha. On this 11 km² ponderosa pine study area, LiDAR data proved most useful for predicting forest biomass. RaDAR ranging measurements did not improve the LiDAR estimates.     

Development of a model-based clinical sepsis biomarker for critically ill patients

Sepsis occurs frequently in the intensive care unit (ICU) and is a leading cause of admission, mortality, and cost. Treatment guidelines recommend early intervention, however positive blood culture results may take up to 48 h. Insulin sensitivity (S_I) is known to decrease with worsening condition and could thus be used to aid diagnosis. Some glycemic control protocols are able to accurately identify insulin sensitivity in real-time.Hourly model-based insulin sensitivity S_I values were calculated from glycemic control data of 36 patients with sepsis. The hourly S_I is compared to the hourly sepsis score (ss) for these patients (ss = 0-4 for increasing severity). A multivariate clinical biomarker was also developed to maximize the discrimination between different ss groups. Receiver operator characteristic (ROC) curves for severe sepsis (ss ≥ 2) are created for both S_I and the multivariate clinical biomarker.Insulin sensitivity as a sepsis biomarker for diagnosis of severe sepsis achieves a 50% sensitivity, 76% specificity, 4.8% positive predictive value (PPV), and 98.3% negative predictive value (NPV) at an S_I cut-off value of 0.00013 L/mU/min. Multivariate clinical biomarker combining S_I, temperature, heart rate, respiratory rate, blood pressure, and their respective hourly rates of change achieves 73% sensitivity, 80% specificity, 8.4% PPV, and 99.2% NPV. Thus, the multivariate clinical biomarker provides an effective real-time negative predictive diagnostic for severe sepsis. Examination of both inter- and intra-patient statistical distribution of this biomarker and sepsis score shows potential avenues to improve the positive predictive value.     

An assessment of chlorophyll-a algorithms available for SeaWiFS in coastal and open areas of the Bay of Bengal and Arabian Sea

Three ocean colour algorithms, OC4v6, Carder and OC5 were tested for retrieving Chlorophyll-a (Chla) in coastal areas of the Bay of Bengal and open ocean areas of the Arabian Sea. Firstly, the algorithms were run using ~ 80 in situ Remote Sensing Reflectance, (R_rs(λ)) data collected from coastal areas during eight cruises from January 2000 to March 2002 and the output was compared to in situ Chla. Secondly, the algorithms were run with ~ 20 SeaWiFS R_rs(λ) and the results were compared with coincident in situ Chla. In both cases, OC5 exhibited the lowest log₁₀-RMS, bias, had a slope close to 1 and this algorithm appears to be the most accurate for both coastal and open ocean areas. Thirdly the error in the algorithms was regressed against Total Suspended Material (TSM) and Coloured Dissolved Organic Material (CDOM) data to assess the co-variance with these parameters. The OC5 error did not co-vary with TSM and CDOM. OC4v6 tended to over-estimate Chla > 2 mg m⁻³ and the error in OC4v6 co-varied with TSM. OC4v6 was more accurate than the Carder algorithm, which over-estimated Chla at concentrations > 1 mg m⁻³ and under-estimated Chla at values < 0.5 mg m⁻³. The error in Carder Chla also co-varied with TSM. The algorithms were inter-compared using > 5500 SeaWiFS R_rs(λ) data from coastal to offshore transects in the Northern Bay of Bengal. There was good agreement between OC4v6 and OC5 in open ocean waters and in coastal areas up to 2 mg m⁻³. There was a strong divergence between Carder and OC5 in open ocean and coastal waters. OC4v6 and Carder tended to over-estimate Chla in coastal areas by a factor of 2 to 3 when TSM > 25 g m⁻³. We strongly recommend the use of OC5 for coastal and open ocean waters of the Bay of Bengal and Arabian Sea. A Chla time series was generated using OC5 from 2000 to 2003, which showed that concentrations at the mouths of the Ganges reach a maxima (~ 5 mg m⁻³) in October and November and were 0.08 mg m⁻³ further offshore increasing to 0.2 mg m⁻³ during December. Similarly in early spring from February to March, Chla was 0.08 to 0.2 mg m⁻³ on the east coast of the Bay.     

Burn severity estimation from remotely sensed data: Performance of simulation versus empirical models

Burn severity is a key factor in post-fire assessment and its estimation is traditionally restricted to field work and empirical fitting from remotely sensed data. However, the first method is limited in terms of spatial coverage and cost effectiveness and the second is site- and data-specific. Since alternative approaches based on radiative transfer models (RTM) have been usefully applied in retrieving several biophysical plant parameters (leaf area index, water and dry matter content, chlorophyll), this paper has applied the inversion of a simulation model to estimate burn severity in terms of the Composite Burn Index (CBI). The performance of the model inversion method was compared to standard empirical techniques. The study area chosen was a large forest fire in central Spain which occurred in July 2005. The model inversion showed the most accurate estimation for high severity levels (for CBI > 2.7, RMSE = 0.30) and for unburned areas (CBI < 0.5, RMSE = 0). In both methodologies, the error associated to CBI from 0.5 to 2.7 was not acceptable (RMSE > 0.7), because it is higher than 25% of the total range of the index. Finally, burn severity maps from both methods were compared.     

Evaluation of MODIS land surface temperature data to estimate air temperature in different ecosystems over Africa

The estimation of near surface air temperature (Ta) is useful for a wide range of applications such as agriculture, climate related diseases and climate change studies. Air temperature is commonly obtained from synoptic measurements in weather stations. In Africa, the spatial distribution of weather stations is often limited and the dissemination of temperature data is variable, therefore limiting their use for real-time applications. Compensation for this paucity of information may be obtained by using satellite-based methods. However, the derivation of near surface air temperature (Ta), from the land surface temperature (Ts) derived from satellite is far from straight forward. Some studies have tried to derive maximum Ta from satellites through regression analysis but the accuracy obtained is quite variable according to the study. The main objective of this study was to explore the possibility of retrieving high-resolution Ta data from the Moderate Resolution Imaging Spectroradiometer (MODIS) Ts products over different ecosystems in Africa. First, comparisons between night MODIS Ts data with minimum Ta showed that MODIS nighttime products provide a good estimation of minimum Ta over different ecosystems (with (ΔTs − Ta) centered at 0 °C, a mean absolute error (MAE) = 1.73 °C and a standard deviation = 2.4 °C). Secondly, comparisons between day MODIS Ts data with maximum Ta showed that (ΔTs − Ta) strongly varies according to the seasonality, the ecosystems, the solar radiation, and cloud-cover. Two factors proposed in the literature to retrieve maximum Ta from Ts, i.e. the Normalized Difference Vegetation Index (NDVI) and the Solar Zenith Angle (SZA), were analyzed. No strong relationship between (ΔTs − Ta) and (i) NDVI and (ii) SZA was observed, therefore requiring further research on robust methods to retrieve maximum Ta.     

Real-time stereo matching using memory-efficient Belief Propagation for high-definition 3D telepresence systems

New generations of telecommunications systems will include high-definition 3D video that provides a telepresence feeling. These systems require high-quality depth maps to be generated in a very short time (very low latency, typically about 40 ms). Classical Belief Propagation algorithms (BP) generate high-quality depth maps but they require huge memory bandwidths that limit low-latency implementations of stereo-vision systems with high-definition images.This paper proposes a real-time (latency inferior to 40 ms) high-definition (1280 × 720) stereo matching algorithm using Belief Propagation with good immersive feeling (80 disparity levels). There are two main contributions. The first is an improved BP algorithm with pixel classification that outperforms classical BP while reducing the number of memory accesses. The second is an adaptive message compression technique with a low performance penalty that greatly reduces the memory traffic. The combination of these techniques outperforms classical BP by about 6.0% while reducing the memory traffic by more than 90%.     

Neural network estimation of chlorophyll a from MERIS full resolution data for the coastal waters of Galician rias (NW Spain)

In typical Case 2 waters, accurate remote sensing retrieval of chlorophyll a (chla) is still a challenging task. In this study, focusing on the Galician rias (ΝW Spain), algorithms based on neural network (NN) techniques were developed for the retrieval of chla concentration in optically complex waters, using Medium Resolution Imaging Spectrometer (MERIS) data. There is considerable interest in the accurate estimation of chla for the Galician rias, because of the economic and social importance of the extensive culture of mussels, and the high frequency of harmful algal events. Fifteen MERIS full resolution (FR) cloud-free images paired with in situ chla data (for 2002-2004 and 2006-2008) were used for the development and validation of the NN. The scope of NN was established from the clusters obtained using fuzzy c-mean (FCM) clustering techniques applied to the satellite-derived data. Three different NNs were developed: one including the whole data set, and two others using only points belonging to one of the clusters. The input data for these latter two NNs was chosen depending on the quality level, defined on the basis of quality flags given to each data set. The fitting results were fairly good and proved the capability of the tool to predict chla concentrations in the study area. The best prediction was given for the NN trained with high-quality data using the most abundant cluster data set. The performance parameters in the validation set of this NN were R² = 0.86, mean percentage error (MPE) = − 0.14, root mean square error (RMSE) = 0.75 mg m^− 3, and relative RMSE = 66%. The NN developed in this study detected accurately the peaks of chla, in both training and validation sets. The performance of the Case-2-Regional (C2R) algorithm, routinely used for MERIS data, was also tested and compared with our best performing NN and the sea-truthing data. Results showed that this NN outperformed the C2R, giving much higher R² and lower RMSE values.This study showed that the combination of in situ data and NN technology improved the retrieval of chla in Case 2 waters, and could be used to obtain more accurate chla maps. A local-based algorithm for the chla retrieval from an ocean colour sensor with the characteristics of MERIS would be a great support in the quantitative monitoring and study of harmful algal events in the coastal waters of the Rias Baixas. The limitations and possible improvements of the developed chla algorithms are also discussed.     

Aboveground biomass retrieval in tropical forests — The potential of combined X- and L-band SAR data use

In the context of reducing emissions from deforestation and forest degradation (REDD) and the international effort to reduce anthropogenic greenhouse gas emissions, a reliable assessment of aboveground forest biomass is a major requirement. Especially in tropical forests which store huge amounts of carbon, a precise quantification of aboveground biomass is of high relevance for REDD activities. This study investigates the potential of X- and L-band SAR data to estimate aboveground biomass (AGB) in intact and degraded tropical forests in Central Kalimantan, Borneo, Indonesia. Based on forest inventory data, aboveground biomass was first estimated using LiDAR data. These results were then used to calibrate SAR backscatter images and to upscale the biomass estimates across large areas and ecosystems. This upscaling approach not only provided aboveground biomass estimates over the whole biomass range from woody regrowth to mature pristine forest but also revealed a spatial variation due to varying growth condition within specific forest types. Single and combined frequencies, as well as mono- and multi-temporal TerraSAR-X and ALOS PALSAR biomass estimation models were analyzed for the development of accurate biomass estimations. Regarding the single frequency analysis overall ALOS PALSAR backscatter is more sensitive to AGB than TerraSAR-X, especially in the higher biomass range (> 100 t/ha). However, ALOS PALSAR results were less accurate in low biomass ranges due to a higher variance. The multi-temporal L- and X-band combined model achieved the best result and was therefore tested for its temporal and spatial transferability. The achieved accuracy for this model using nearly 400 independent validation points was r² = 0.53 with an RMSE of 79 t/ha. The model is valid up to 307 t/ha with an accuracy requirement of 50 t/ha and up to 614 t/ha with an accuracy requirement of 100 t/ha in flat terrain. The results demonstrate that direct biomass measurements based on the synergistic use of L- and X-band SAR can provide large-scale AGB estimations for tropical forests. In the context of REDD monitoring the results can be used for the assessment of the spatial distribution of the biomass, also indicating trends in high biomass ranges and the characterization of the spatial patterns in different forest types.     

An approach to quantify depth-resolved marine photochemical fluxes using remote sensing: Application to carbon monoxide (CO) photoproduction

This paper presents a model to calculate depth-resolved marine photochemical fluxes from remotely sensed ocean color and modeled solar irradiance. The basic approach uses three components: 1) below-sea-surface spectral downward scalar irradiance calculated from a radiative transfer model (STAR) and corrected for clouds using TOMS UV reflectivities; 2) surface-ocean spectral diffuse attenuation coefficients and absorption coefficients for chromophoric dissolved organic matter retrieved from SeaWiFS ocean color using the SeaUV/SeaUVc algorithms; and (3) spectral apparent quantum yield for the photochemical reaction considered. The output of the model is a photochemical rate profile, Ψ_PR(z), where z represents depth.We implemented the model for carbon monoxide (CO) photochemistry using an average apparent quantum yield spectrum and generated a monthly climatology of depth-resolved CO photoproduction rates in the global ocean. The climatology was used to compute global budgets and investigate the spatial and seasonal variabilities of CO photoproduction in the ocean. The model predicts a global CO photoproduction rate of about 41 TgC yr^− 1, in good agreement with other recent published estimates ranging from 30 to 84 TgC yr^− 1. The fate of photochemically derived CO and its role in global biogeochemical cycles remains uncertain however, with biological consumption and sea-air exchange competing for its removal in the surface ocean. Knowledge of the vertical distribution of CO photoproduction is critical in the quantification of the relative magnitudes of these sink mechanisms. The depth-resolution capabilities of this model, together with US Naval Research Laboratory climatology for mixed layer depths allowed further estimation that > 95% of the total water-column CO photoproduction occurs within the mixed layer on a global, yearly basis. Despite this compelling figure, the model also suggests significant spatio-temporal variability in the vertical distribution of CO photoproduction in the subtropical gyres, where up to 40% of water-column CO can be produced below the mixed layer during summertime.While the approach can be applied to other photochemical fluxes (e.g. DIC formation or DMS removal), accurate quantification of such processes with remote sensing will be limited until the mechanisms regulating observed oceanic variability in the apparent quantum yields are better understood. Minor modification to this model can also make it applicable for the determination of the effects of UV and visible solar radiation on sensitive biological systems.     

A colorization algorithm based on local MAP estimation

This paper presents a colorization algorithm which adds color to monochrome images. In this paper, the colorization problem is formulated as the maximum a posteriori (MAP) estimation of a color image given a monochrome image. Markov random field (MRF) is used for modeling a color image which is utilized as a prior for the MAP estimation. The MAP estimation problem for a whole image is decomposed into local MAP estimation problems for each pixel. Using 0.6% of whole pixels as references, the proposed method produced pretty high quality color images with 25.7-32.6 dB PSNR values for eight images.