Research article

Using rheological monitoring to determine the gelation kinetics of chitosan-based systems


  • The modeling of polymeric reactions is a topic of large interest. The gelation reactions that may result from self-crosslinking or hybrid (agent based-) crosslinking are examples with interest specially in biomaterials applications. The composition of polymer entities during the reaction is hard to follow, and their concentration is not a good measure of the system dynamics. One alternative is monitoring the rheological behavior of the reacting mass, and relate the elastic modulus of the mixture with the rheological degree of conversion. In this paper we use rheological data to fit Malkin and Kulichikin (1996) [1] based models to describe the crosslinking of chitosan. First, the self-crosslinking of chitosan is considered. Then, the agent-based crosslinking reaction promoted by genipin is addressed. We use dynamical rheological data to fit the reaction models. The model fitting problem generated using Maximum Likelihood principle with heteroscedastic prediction error variance is formulated as a Dynamic Optimization problem and subsequently solved with a sequential approach. Parametric confidence regions are computed using the linear approximation of the covariance matrix at the optimum. Further, the parameters correlation matrix is also determined and used to qualitatively infer about the practical identifiability. The reaction order obtained for self-crosslinking kinetics is 1.3375 ± (0.0151) – approximately of first order –, and is 2.2402 ± (0.0373) for hybrid crosslinking (approximately of second order). In both cases we prove the error variance model is heteroskedastic and the model is identifiable. The approach proposed herein can be extended to other polymer systems.

    Citation: Belmiro P. M. Duarte, Maria J. Moura. Using rheological monitoring to determine the gelation kinetics of chitosan-based systems[J]. Mathematical Biosciences and Engineering, 2023, 20(1): 1176-1194. doi: 10.3934/mbe.2023054

    Related Papers:

    [1] Muhammad Shahid Iqbal, Muhammad Nauman Ahmad, Nynke Hofstra . The Relationship between Hydro-Climatic Variables and E. coli Concentrations in Surface and Drinking Water of the Kabul River Basin in Pakistan. AIMS Environmental Science, 2017, 4(5): 690-708. doi: 10.3934/environsci.2017.5.690
    [2] Alejandra S. Coronel, Susana R. Feldman, Emliano Jozami, Kehoe Facundo, Rubén D. Piacentini, Marielle Dubbeling, Francisco J. Escobedo . Effects of urban green areas on air temperature in a medium-sized Argentinian city. AIMS Environmental Science, 2015, 2(3): 803-826. doi: 10.3934/environsci.2015.3.803
    [3] Kamonrat Suphawan, Kuntalee Chaisee . Gaussian process regression for predicting water quality index: A case study on Ping River basin, Thailand. AIMS Environmental Science, 2021, 8(3): 268-282. doi: 10.3934/environsci.2021018
    [4] Alessio Russo, Francisco J Escobedo, Stefan Zerbe . Quantifying the local-scale ecosystem services provided by urban treed streetscapes in Bolzano, Italy. AIMS Environmental Science, 2016, 3(1): 58-76. doi: 10.3934/environsci.2016.1.58
    [5] John R. Christy . Examination of Extreme Rainfall Events in Two Regions of the United States since the 19th Century. AIMS Environmental Science, 2019, 6(2): 109-126. doi: 10.3934/environsci.2019.2.109
    [6] Kathleen D Reinhardt, Wirdateti, K.AI. Nekaris . Climate-mediated activity of the Javan Slow Loris, Nycticebus javanicus. AIMS Environmental Science, 2016, 3(2): 249-260. doi: 10.3934/environsci.2016.2.249
    [7] N. Naveena, G. Ch. Satyanarayana, A. Dharma Raju, K Sivasankara Rao, N. Umakanth . Spatial and statistical characteristics of heat waves impacting India. AIMS Environmental Science, 2021, 8(2): 117-134. doi: 10.3934/environsci.2021009
    [8] Tewodros Woldemariam Tesfaye, C.T. Dhanya, A.K. Gosain . Evaluation of ERA-Interim, MERRA, NCEP-DOE R2 and CFSR Reanalysis precipitation Data using Gauge Observation over Ethiopia for a period of 33 years. AIMS Environmental Science, 2017, 4(4): 596-620. doi: 10.3934/environsci.2017.4.596
    [9] Marta Rodrigues, Henrique Queiroga, Anabela Oliveira, Vanda Brotas, Maria D. Manso . Climatic and anthropogenic factors driving water quality variability in a shallow coastal lagoon (Aveiro lagoon, Portugal): 1985–2010 data analysis. AIMS Environmental Science, 2016, 3(4): 673-696. doi: 10.3934/environsci.2016.4.673
    [10] Dong Chen, Marcus Thatcher, Xiaoming Wang, Guy Barnett, Anthony Kachenko, Robert Prince . Summer cooling potential of urban vegetation—a modeling study for Melbourne, Australia. AIMS Environmental Science, 2015, 2(3): 648-667. doi: 10.3934/environsci.2015.3.648
  • The modeling of polymeric reactions is a topic of large interest. The gelation reactions that may result from self-crosslinking or hybrid (agent based-) crosslinking are examples with interest specially in biomaterials applications. The composition of polymer entities during the reaction is hard to follow, and their concentration is not a good measure of the system dynamics. One alternative is monitoring the rheological behavior of the reacting mass, and relate the elastic modulus of the mixture with the rheological degree of conversion. In this paper we use rheological data to fit Malkin and Kulichikin (1996) [1] based models to describe the crosslinking of chitosan. First, the self-crosslinking of chitosan is considered. Then, the agent-based crosslinking reaction promoted by genipin is addressed. We use dynamical rheological data to fit the reaction models. The model fitting problem generated using Maximum Likelihood principle with heteroscedastic prediction error variance is formulated as a Dynamic Optimization problem and subsequently solved with a sequential approach. Parametric confidence regions are computed using the linear approximation of the covariance matrix at the optimum. Further, the parameters correlation matrix is also determined and used to qualitatively infer about the practical identifiability. The reaction order obtained for self-crosslinking kinetics is 1.3375 ± (0.0151) – approximately of first order –, and is 2.2402 ± (0.0373) for hybrid crosslinking (approximately of second order). In both cases we prove the error variance model is heteroskedastic and the model is identifiable. The approach proposed herein can be extended to other polymer systems.



    An increase of the mean surface air temperature (SAT) presents the most prominent evidence of the present-day climate change (climate warming) both on global and local scales [1,2,3,4]. However, climate changes manifest themselves not only in the evolution of mean climate characteristics, but also in the changes of the climate variability and in the occurrence of climate extremes. Potential intensification of the anthropogenic forcing influences thus not only the climate state but becomes involved in the variability and instability of weather and climate. The precise knowledge of the air temperature variability and of the frequency of the occurrences of temperature extremes may be to some degree even more crucial than the proper impact of the gradual increase of the mean air temperature [5,6,7]. Especially the sudden temperature extremes can trigger subsequent climatic and natural events that consequently endanger various sectors of life (e.g. strongly variable precipitation, floods, rising sea level, warmer water temperatures, effect on plant development, etc.). All these events have serious consequences on human every-day-life and on the difficulty of the society adaptation to the unexpected intense weather changes [8].

    In the last dozen years numerous investigations focused on the ongoing changes in air temperature variability on different time scales [9,10,11]. Previously such studies were based on the global gridded SAT data sets. However, in the last two or three decades the warming trend and the related climatic changes have become so strong that they become detectable in many local weather stations. The climate change studies at the regional level provide more detailed picture as they facilitate the possibility to understand and compare the local specific manifestations [12,13,14,15]. While the regionally detected warming trends are generally similar, the temperature variability exhibits a more diverse pattern with extensive regions with either positive and/or negative climate variability trends. E.g. Xu et al. [15] analyzed the variability of different daily climatic indices (including air temperatures) using ~50-years long high spatial resolution data set in China and found that the daily variability of climate was higher in North China than in South China. On the other hand, Huntingford et al. [16] argued for no increase in global temperature variability even if annual temperature variations show noticeable geographical variations. Scale differences (annual, seasonal, daily) and especially methodological differences between various studies represent additional problem because some of the obtained results are difficult to be compared. Some researchers perform the detection of variability by generally simple statistics like comparing the mean values and their variance [14,17,18], while others prefer more advanced statistical methods, like the N-point change method [19,20,21] or the fractal dimensional analysis [15]. Geographical complexity of the variability pattern and the methodological differences in its detection have produced ongoing serious debates in the interpretation of the results [22] but have inspired numerical research attempts for a broader variability assessment all over the world.

    In the present work we have tried to investigate the SAT variability in the temperature-time series data obtained at the observation station Prague-Sporilov (the Czech Republic) in years 2003–2016. We continue the initial similar studies which focused on temperature record 1994–2001. We described that climate systems are variable on all time scales and that their evolution could be better understood when the insight was aimed not only on the changes in mean characteristics, but the climate variability and climate extremes were considered [20]. The interpretation of the Sporilov data were completed and compared with the results of the long-term observation of the Prague Klementinum station to confirm the Prague-Sporilov results.

    Daily Average Temperatures (DAT) were calculated from the recorded SAT series measured at 2 m height at the Geothermal Climate Change Observatory (GCCO) Sporilov, Prague (50°02.43' N, 14°28.54x E, 275 m a.s.l.). The experiment started in 1992 when the GCCO was built as a part of the long-term paleoclimate project to understand the mechanism of the downward propagation of the ground surface temperature signal into shallow subsurface [23,24,25]. The DAT used here are a part of observational material covering the time interval 2003–2016 (Figure 1) and the present interpretation continues the previous similar studies [20]. Original data were recorded by a 16-channel data logger M4016 [26]; the proper temperature sensor (platinum Pt1000, four-wire mode) located in a waterproof tube is placed on the top of the 2 m long rod. The 5-minutes original input data were stored and hourly means computed to interpret daily variations and further computed daily means when describing the inter-annual variation patterns. Battery failures, occasional sensor failures, and unstable voltage inputs in the data logger unfortunately disrupted the compactness of the registration record. General data coverage over the whole registration period varied between 96 and 97%. No serious attempt can be made to statistically improve the observed 5-min data series by infilling the existing gaps. However, the entire duration of the whole 14-year-long observational interval and its 96% data coverage is a certain guarantee to identify the obtained variability results. In contrast to the previous work [20], when the variability was investigated in the temperature time series measured between 1994–2001 and the calculation was performed for surface (z = 0 m) air temperatures data averaged for 6-hour and for 24-hour intervals, now we preferred air temperature at 2 m height and have focused on daily means, to be able to compare the GCCO-Sporilov results with the data from the Prague-Klementinum record (see further). The GCCO-site, situated at the campus of the Institute of Geophysics, is located on the top of the local gentle crest framing the southern periphery of Prague. The area can be characterized as a surviving island of typical suburban houses, with a plenty of greenery, on the outskirts of the quickly expanding robust development of the growing city.

    Figure 1.  Prague-Sporilov site: 14 years (2003–2016) monitoring of surface air temperature. Fourier smoothed representation Y = a*SIN(0.0172024*X+b)+c*X+d, a = 10.346 deg, b = 17.022, c = 0.0002547 deg/day (0.0930 deg/year), d = 9.9612 deg.

    The Klementinum station (50°05.20' N, 14°24.98' E, 191 m a.s.l.) is the oldest Czech meteorological observatory, with the systematic air temperature measurements started in 1775 [27]. Contrary to Sporilov site, the station is located in the very downtown of Prague and may serve to demonstrate the influence of the urban environment. Even when in the past its instrumentation underwent several modifications, the almost 250-year-long observational record presents valuable scientific material. Here for practical purpose we have analyzed only the comparable time interval as provided the Sporilov station, namely years 2003–2016 (Figure 2), and added few more findings analyzing longer time spans 1900–2016 and 1970–2016. In addition to the daily averaged temperatures (DAT) the results include also the daily maximum Tmax and daily minimum Tmin temperatures.

    Figure 2.  14-year section (2003–2016) of the air temperature recorded at the Prague-Klementinum meteorological station. Fourier smoothed representation of the daily averaged temperatures Y = a*SIN(0.0172024*X+b)+c*X+d, a = 10.260 deg, b = 17.012, c = 0.0001967 deg/day (0.0718 deg/year), d = 3.587 deg.

    Numerous observations have shown that during the second half of the 20th century the mean SAT has been increasing on daily, monthly and annual scales overall in Europe [1,28,29] and that since the late 1970s the mean temperature increase has been substantially accelerated. This phenomenon is observable not only in the SAT series themselves but on the seasonal scales and is especially evident in the summer temperatures [30]. Figure 3 shows seasonal linear trends calculated for the Prague-Sporilov temperature series, similar calculations were done for the Klementinum station. It is interesting to mention that the most intensive warming always corresponds to the winter season. The spring and summer warming are approx. two-three times lower and the lowest warming rate for Sporilov corresponds to the fall season (more than an order lower than the winter). The difference between summer and winter temperatures reduces the decrease of climate variability on the annual scale and produces a shift in the duration of the individual seasons (spring comes earlier and cold weather/frosts are arriving later in the fall). Table 1 summarizes the parameters of the linear trends calculated for the interval 2003–2016 of both investigated stations. As seen, Klementinum daily air temperature averages are by about 0.5 to 1 K higher that Sporilov averages (the fact that can be attributed to the urban heat island effect, UHI), but the Sporilov warming rate is slightly higher.

    Figure 3.  Seasonal daily mean air temperatures at Prague-Sporilov station and corresponding linear warming trends. MMA means March, May, April; JJ –June, July, August; SON – September, October, November and DJF–December, January, February.
    Table 1.  SAT linear trends for Prague-Sporilov and Klementinum data series during 2003–2016 years.
    Data Mean temperature (deg.C) Warming (K/decade)
    Sporilov
    Daily temp. 10.61 ± 0.12 1.16 ± 0.29
    Spring temp. 10.64 0.93
    Summer temp. 20.04 0.83
    Fall temp. 10.31 0.11
    Winter temp. 1.28 2.00
    Klementinum
    Daily temp. 11.49 ± 0.11 0.97 ± 0.28
    Spring temp. 11.42 0.45
    Summer temp. 20.84 0.30
    Fall temp. 11.27 0.44
    Winter temp. 2.25 1.83
    Tmax 15.29 ± 0.13 0.90 ± 0.32
    Tmin 7.94 ± 0.10 1.16 ± 0.25

     | Show Table
    DownLoad: CSV

    The UHI can noticeably affect the measured air temperatures (both their trends and fluctuations) and its intensity need not be constant, but depends on various factors related to town development dynamics, such as population and building growth, increasing heat production, energy consumption, etc. Detailed investigations of the UHI effect for the Prague-Klementinum station were performed [31,32]. The former work analyzed daily mean temperatures measured between 1921–1995 at 17 rural stations located within 40 km range from Prague-Klementinum. The later work was based on the daily minimum temperatures measured during 1961–1990 at 3 stations situated in the 50 km neighbourhood of the Klementinum station. The Prague UHI intensity has shown progressive intensification since 1922 (especially from mid 40's to mid 60's), and certain stagnation after 60s. The trend was practically zero between 1980–1995. The highest daily mean temperature increase was observed in winter and spring (0.06 K/decade) and the smallest, practically insignificant (0.01 K/decade) in summer. The annual warming was ~0.05 K/decade. The UHI increase is thus 20–25 times lower than the air temperature warming trends presented in Table 1 and do not exceed the uncertainty of warming trend. The elimination of the UHI effect from the warming trends will not markedly change the observed linear trends. This is clearly visible in Figure 5 [32]), where measured temperatures and temperatures reduced by the intensification of the UHI according to the linear trend in 1922–1995 practically coincide. The UHI intensification trend for the daily minimum temperatures is somewhat higher than calculated for daily mean temperatures (0.12 K/decade within 1961–1990 period [31]). However, it is still an order of magnitude lower than that presented in Table 1 and is smaller than the uncertainty of calculated warming trend.

    The temperature trends exhibit similar patterns on the annual scale. Figure 4 shows annual mean temperatures as well as maximum (hottest day) and minimum (coldest day) daily temperatures for the period 2003–2016 at the Sporilov station. The range of annual daily temperature variations fluctuates between 32.3–44.0 K. Positive/warming trends are characteristic for all three variables. However, an increase of annual coldest day temperature is ~1.5 times faster than the temperature of the annual hottest day. It hints that the climate faster becomes less cold rather than getting warmer. Similar conclusion is presented in the work by Alexander et al. [33].

    Figure 4.  Annual mean air temperatures and annual maximum and minimum daily temperatures at Prague-Sporilov station and corresponding linear trends.

    Table 2 presents similar results related to several longer periods calculated for the Klementinum station. Warming trends and average daily temperatures are growing with time and this growth increased sharply in the last few decades. The present mean temperature has increased by 3.2 K compared with year 1900, when one third of this increase occurred since year 2000. It may be interesting to mention that according to the National Climate Assessment (NCA) the average temperature on the USA territory has increased by 1.3–1.9 K since 1895 and most of this increase is also due to the last couple decades [34]. The observed Klementinum temperature increase is thus 1.7 times greater. Such increase is in a good agreement with similar studies in Western Europe where last few decades have been getting warmer faster than the global average and far stronger than predicted by the standard models based on the anthropogenic greenhouse gas increases [12,35,36].

    Table 2.  SAT linear trends for Prague-Klementinum station calculated for three periods of time.
    Period 1900–1950 1950–2000 2000–2016
    DAT
    Mean (deg.C) 9.47 ± 0.06 10.10 ± 0.06 11.47 ± 0.10
    Slope (K/decade) 0.159 ± 0.042 0.267 ± 0.041 0.607 ± 0.209
    Tmax
    Mean (deg.C) 13.10 ± 0.07 13.83 ± 0.08 15.26 ± 0.12
    Slope (K/decade) 0.334 ± 0.047 0.221 ± 0.047 0.625 ± 0.240
    Tmin
    Mean (deg.C) 5.97 ± 0.05 6.47 ± 0.05 7.92 ± 0.09
    Slope (K/decade) 0.003 ± 0.038 0.296 ± 0.036 0.728 ± 0.184

     | Show Table
    DownLoad: CSV

    Both daily minimum and daily maximum temperatures have been increasing with time as well. The calculated trends are generally similar for mean, maximum and minimum temperatures and exhibit an increase on both monthly and annual scales. The magnitudes of both trends are comparable with the rate of the daily average temperature warming. Noticeable feature is that minimum daily temperatures (Tmin) increase approx. 30% faster than daily maxima (Tmax). This result in a progressive decrease of the diurnal temperature range (DTR), and decrease in the intra-day variability and daily temperatures have been getting warmer. The observed evidence, namely a faster increase in the diurnal minimum temperature compared with the diurnal maximum temperature in the last 50 years leading to a decrease in the DTR is consistent with the mean global trend of the DTR reducing [37,38,39]. It should be stressed, however, that the DTR trend course is not stable, but to some degree variable and may incorporate temporally rate variation (slowdown) or even reversals to the negative trend in recent decades [40].

    The processes responsible for the DTR decrease are still poorly apparent as well as the role of the potential anthropogenic forcing. Some authors have argued that shortwave radiation could have a strong and sufficiently homogeneous effect to change DTR on a global scale [41,42]. Anyhow, as shown above, the decreasing DTR trend represents an exact and easily identifiable characteristic of the recent climate change, and its detection in the forcing mechanisms of the present-day climate change (in particular in its potential anthropogenic component) may be important.

    Variability indicates the degree of fluctuation and uncertainty of the temperature change process. One of the remarkable aspects of the global warming is the fact that the global mean temperature (SAT) increase may be accompanied by certain temperature variability decrease. Climate simulations studies predict not only climate warming, but also a general decrease in climate variability [19,43,44]. It was recognized that the variability decreases on average, it was also pointed out that these trends are roughly constant or linear in time. The variability trend may be likewise frequency dependent [20].

    There are several ways to detect variability changes. Here we used so-called N-point change method introduced in detail earlier [20,21]. The measured data series were converted into temperature anomalies by means of detrending and removing the annual wave. Figure 5 presents variability patterns calculated for the daily temperature averages for both Prague-Sporilov site and Klementinum station for the 2003–2016 period. The variability shows an insignificant decreasing trend of −0.003 K/yr that is by an order of magnitude lower than the value −0.038 ± 0.002 K/yr obtained earlier for this site for the 1994–2001 period [20]. The above Sporilov variability trends are quite similar to the variability trends for the Prague-Klementinum data, namely −0.019 K/yr and/or −0.005 K/yr for the 1994–2001 and 2003–2016 period, respectively. Figure 6 and Table 3 summarizes the variability trends calculated from the Prague-Klementinum data for the 10-years intervals within 1971–2016 period. As seen, all trends are negative (variability decreases) and this decrease accelerates with time. This decrease was practically negligible in the beginning of the studied period and has grown by more than two orders of magnitude within last about 5 decades.

    Figure 5.  Variability of the SAT data at Sporilov and Klementinum stations 2003–2016.
    Figure 6.  Variability of the SAT data at the Klementinum station 1971–2016.
    Table 3.  Parameters of linear trends in variability calculated for different periods from the Prague-Klementinum data.
    Time interval Mean (K) Slope (K/decade)
    1971.1.1–1980.12.31 1.70 -0.001
    1981.1.1–1990.12.31 1.808 -0.044
    1991.1.1–2000.12.31 1.826 -0.079
    2001.1.1–2010.12.31 1.742 -0.13
    2011.1.1–2016.12.31 1.812 -0.32

     | Show Table
    DownLoad: CSV

    The detected decreasing trend exists only on an average. Its persistency often breaks off both in the value as well as in the sign of the variability change. On the shorter time scales the calculated slope value strongly depends on the length and the beginning of an investigated time interval. The positive slope values were obtained for the Prague-Sporilov data for the spring-summer time intervals in year 2003 and above all in year 2015. Noticeable changes in the variability generally coincided with more pronounced temperature extremes. E.g., two reversals in the variability change mentioned above can be attributed to the warming occurred due to record-breaking European heat waves [36,45].

    The variability trends calculated for data measured at both Sporilov and Klementinum sites exhibit an insignificant decrease within the 2003–2016 time interval. However, it should be emphasized, that such trend corresponds to the calculated value only when calculated for the whole interval as entirety. This fact may lose its validity when calculation is applied to another time section (Figure 7). On a shorter time scale the calculated slope value may strongly depend on the length and the beginning of the selected time interval. In our case, positive slope values were obtained for the 2003 spring-summer time interval and especially for the 2015 summer. Such noticeable changes in the variability generally coincide with pronounced temperature extremes. For example, the two above mentioned reversals in the variability change can be undoubtedly attributed to the periods of a sequence of exceptionally hot days (European heat waves) [36,45].

    Figure 7.  Variability of the SAT data at the Sporilov station (running average for 100 consecutive data), intervals with growing variability highlighted (see text).

    1. The record breaking European temperature heat wave occurred in the summer of 2003 (approx. from June to mid-August), when temperatures were about 20–30 percent above average. It has affected much of Western Europe, but was generally concentrated in England, France and Spain. Temperatures were as high as 47 ℃ in the south of Portugal and the heat wave had extremely hard human consequences [46,47].

    2. An unusually strong and prolonged heat waves occurred across Europe from late June to mid-September 2015 [48]. The record temperatures above 40 ℃ have been recorded in several locations (the Mediterranean coast, south-western, central and south-eastern Europe; see e.g. [49]). The maximum daily temperatures approached to about 37 ℃ during both 2003 and 2015 heat waves in Prague.

    The past 2–3 decades of the 20th Century and the beginning of the 21st Century have been characterized not only by the global temperature increase, but also by an increased frequency of extreme events such as heat/cold waves, heavy rainfall and floods, droughts, etc. Numerous climatic models predict more frequent occurrence of extreme events [50]. An increased occurrence of weather extremes has acquired a high social resonance that initiated an intensive research activity to quantify the variability, frequency and possible evolution of observed extreme events. For example, Nicholls [51] reported, that extreme weather values may work as good indicators of coming climate changes. Weather extremes are relatively rare; a proper assessing of their occurrence creates certain problems. The most frequently used methods are based on probability characteristics and the counts of days exceeding a certain specified threshold value. To not engage only with a few individual events, a list of indices for extreme temperature phenomena has been worked out in the frames by the CLIVAR project of the World Climate Research Program [52]. The guidelines for the recommended temperature-based CLIVAR indices were presented by Klein Tank [53].

    Numerous evidences support that most global land areas have experienced significant warming of both maximum and minimum temperature extremes during at least the second half of the 20th Century. Changes in the occurrence of cold and warm days are not unambiguous. On the whole, rapid warming observed in the last decades was accompanied by increases in summer temperature extremes over Western Europe while the number of unusually cold days has generally decreased [54]. The IPCC AR4 Climate Model calculations have shown that since the middle of the 20th Century the indices of cold extremes had likely decreased and indices of warm extremes increased at the most of global land areas [55]. Similar conclusions were presented in the review work by Meehl et al.[9]. Number of days with unusually cold winter temperatures (particularly extremely cold nights) became less common. This phenomenon is consistent with the general climate warming. In the first approximation the observed trends in the SAT extremes are often interpreted as a simple shift of the total temperature distribution towards higher temperatures [56].

    On the other hand, the trends of the extreme indices based on 60 European weather stations located within the 40–70°N zone stretching from Great Britain to Russia have very local character [57,58]. It was also shown that in spite of the fact that the individual extremely cold months (ECM) represent a frequent phenomenon in Europe, the winter ECMs were observed only at 10 of all investigated weather stations within the 1951–2000 yr period and only at 25 and 27 stations during 2001–2006 and 2007–2011 periods, respectively. The lowest number of ECMs occurred in the 1991–2000 decade, but their frequency and the territorial extent has started to increase since the beginning of the 21th century [57].

    Europe experienced almost similar number of severe winters and the heat waves in the last decades. There were 11 extremely cold winters (ECW) during 1951/52–1979/80 period and another 7 ECWs during 1980/81–2009/2010 interval (the number of ECWs somewhat diminished, but not significantly [58]. At the same periods, the numbers of heat waves occurring in Europe were 5 and 9, respectively [59]. The most recent European heat wave occurred on summer 2015, in Central Europe it started at the end of July and continued during the first ten days of August.

    As Prague belongs to the temperate climate zone, unusually warm as well as unusually cold days are rare. Changes in climate variability and extremes can be visualized in relation to changes in the probability distribution. Figure 8 shows the probability distribution of the temperature excess values (i.e. the differences between detrended observed temperature data and the Fourier sinusoidal representation applied on Sporilov 2003–2016 record). While the previously investigated 1994–2001 data exhibited a right skewed distribution and a more heavy right tail [20], the 2003–2016 probability distribution is close to the normal (Gaussian) distribution. A prominent feature of the 1994–2001 temperature record was the prevalence of extremes in warm seasons, while the number of cold season extremes was practically insignificant. The present disappearance of the right skewness and the shift to the normal distribution can be explained by the increase of cold extremes. The minimal temperature anomaly was a little less than −10 K during 1994–2001 period and the anomaly values lower or equal to −10 K occurred only once for the investigated 8 years. During 2003–2016 period the minimal anomaly equals already to −15.6 K and the number of days with the values equal or lower than −10 K occurred already 29 times. Similar tendency can be observed at the Prague-Klementinum station (Table 4), the number of both warm and cold extremes increases with time; however the number of cold extremes is growing faster. While the A2/A1 relation in seventh years was ~ 6, in the last decade it fell to ~3.

    Figure 8.  Probability distribution of temperature anomalies (Sporilov site, 2003–2016), the Gaussian distribution is shown as reference.
    Table 4.  Number of days with temperature anomalies lower or equal to -10K (A1) and equal or higher +10 K (A2) at the Prague-Klementinum station.
    Period 1971–1980 2001–2010 2007–2016
    A1 3 8 11
    A2 19 24 32

     | Show Table
    DownLoad: CSV

    While the increase of number and the growth of the intensity of warm extremes is a general European phenomenon and can be attributed to the global warming trend and the corresponding increase of the hot days (>25 ℃), the reasons for the growth of cold extremes are so far not clear and need further explanation. The evolution of the number of cold extremes has a local character; Spinoni et al. [60] have investigated the heat and the cold waves in the Carpathian Region and have described different spatial patterns of heat and cold waves. While heat wave events exhibit a general increase in all investigated parameters, cold wave events show a decrease in the area west to the Carpathians and an increase in the area north–east to the Carpathians. In the south Italian high-altitude region Capozzi and Budillon [61] revealed a strong and significant positive trend in the last 40 years in heat waves frequency, intensity and duration, while cold wave events have exhibited a significant and positive trend only in their intensity.

    For further investigations we have used the CLIVAR-recommended percentile-based and the absolute temperature indices. The four temperature percentile indices consider the coldest and warmest deciles for both maximum and minimum temperatures to evaluate the degree of extremes changing. Absolute indices represent maximum or minimum values within the season or within the year and reflect the longer time scale processes. Table 5 represents the CLIVAR-recommended indices based on the 10th and 90th percentiles of maximum and minimum daily temperatures recorded at Prague-Klementinum station. The 10th percentile is a temperature value such that at least 10% of the observations is less or equal to this value, while the rest of observations is greater than or equal to this value. Below we used the following abbreviations: Tmin 10p–Tmin below 10th percentile, Tmax 10p–Tmax below 10th percentile, Tmin 90p–Tmin above 90th percentile, Tmax 90p–Tmax above 90th percentile. The values lower than Tmin10p and/or Tmax10p usually arise at cold nights and cold days, while the values surpassing Tmin90p and Tmax90p are characteristic for the warm nights and warm days, respectively [33,60].

    Table 5.  Percentile values (deg.C) calculated for different time intervals from daily maximum and daily minimum temperatures recorded at the Prague-Klementinum station.
    Interval (years) 1971–1980 2001–2010 2011–2016
    Tmin10p -1, 4 -1, 9 -0, 2
    Tmin90p 15, 2 17, 1 17, 2
    Tmax10p 2, 6 2, 2 4
    Tmax90p 25 27, 5 27, 9

     | Show Table
    DownLoad: CSV

    As seen, all calculated percentiles exhibit significant warming since at least the mid of the 20th century. The permanent warming is characteristic for both cold nights and cold days. The warming characteristic for warm nights and warm days has slowed down or practically stopped in the first two decades of the 21th century. Similar phenomenon is characteristic for the daily maximum temperatures. The maximum temperature observed in the 1971–1980 interval was 34.7 ℃, while in the last two intervals it equals to 37.3 and 37.6 º C, respectively. It is of interest, that at the present time such heat and/or cold extremes are short and rare in the temperate climate zone and can be regarded as synoptic rather than climatic events, unlikely they can seriously distort the characteristic pattern.

    Absolute indices represent maximum or minimum values within a given time interval. They include maximum daily maximum temperature (TXx), maximum daily minimum temperature (TNx), minimum daily maximum temperature (TXn), minimum daily minimum temperature (TNn). The index TXx can be generally attributed to extremely warm day and TXn to the coldest day, while the TNn index corresponds to extremely cold night and TNx to the warmest night [62,63]. In the present work we calculated the absolute indices from maximum and minimum daily temperatures measured at the Klementinum Observatory. Table 6 summarizes the indices calculated for the 10-years long intervals.

    Table 6.  Decadal absolute indices (deg.C) calculated from the data measured at the Klementinum Observatory during the 1971–2016 period.
    Time interval TXx TXn TNx TNn
    1971–1980 34.7 -11.0 22.6 -17.3
    1981–1990 37.8 -14.3 22.2 -20.2
    1991–2000 36.1 -10.6 23.0 -19.8
    2001–2010 37.3 -9.3 23.8 -14.8
    2011–2016 37.6 -8.7 24.0 -17.1

     | Show Table
    DownLoad: CSV

    The absolute temperature indices exhibit a similar growing pattern as other investigated indices. The maximum daily maximum temperature (TXx) increased by nearly 3 K (the hottest days became hotter) which represents the greatest change in the absolute temperature indices. Similar increase was obtained for the most of investigated regions on the global scale [62,64]. An interesting local feature is that the most of the temperature increase occurred after 1980 and further has significantly slowed down. Similar, but somewhat weaker, warming are found for coldest days (TXn) and warmest nights (TNx), the process which is characteristic over the most of Europe [64].

    The TNn index variations exhibit certain unusual regional detail. Donat et al. [64] have investigated the behavior of the absolute indices on the global scale and found that the warming trends are generally stronger for the coldest night temperatures than for indices calculated for warmest days. Significant increase of the TNn values was detected at the 70% of the grid boxes from 1951 to 2010, while significant decrease was found only in 3% of investigated data. The TNn increase was particularly strong over large parts of Asia. Globally averaged TNn temperature has increased by about 3 K in the last 1951–2010 period. According to Alexander et al. [33] the minimum daily minimum temperature (TNn) increased by ~ 5 K between 1951 and 2003 that represents the greatest change in the absolute temperature indices. However above authors have noted that this growth is dominated by the very large negative anomalies in the 1950s when there was increased sampling uncertainty. Our TNn indices show up to ~ 3–5 K differences between neighboring decadal values, but do not exhibit any significant warming/cooling trend.

    In the present work we have revealed some important features describing the daily average air temperature (DAT) changes, variability and their temporal evolution at the Prague-Sporilov site. The variability and trends in daily near-surface air temperatures recorded during the 2003–2016 period were analyzed. The results were compared with the findings of our previous work covering the earlier 1994–2001 period [20,21]. We compared above data with the longer DAT, Tmax and Tmin series measured at the Prague-Klementinum Station. The results of both stations have shown very similar changes in the investigated statistical parameters as well as in the temporal trends. Such similarity increases the reliability of the obtained results.

    The following principal conclusions can be drawn:

    (1) The results demonstrated increasing warming trends for all investigated temperature series as well as decreasing trends in the diurnal temperature range (DTR) that hints a decrease of intra-daily temperature variability.

    (2) The results revealed the generally reducing temperature variability. However, the velocity of the variability decrease has slowed down by approx. an order of magnitude in comparison with the previously investigated 1994–2001 period. It was also recognized that the persistence of detected processes exists on average, e.g. in the annual scales, and may break down during shorter periods. The pattern of the variability changes during the 2003–2016 period is found to be more diverse than in our earlier studies [20].

    (3) The study reveals that the occurrence of both hottest and coldest extremes have increased during the investigated period at both stations. While an increase of the hot extremes represents a generally global feature, an increase of cold extremes is a local phenomenon. However, as the frequency of extreme days is low, they have no significant influence on most of detected trends.

    This work was supported by the research infrastructure RINGEN co-funded by the Operational Program "Research, Development and Education" (CZ.02.1.01/0.0/0.0/16_013/0001792) and by the Project LM2015079 CzechGeo/EPOS, both provided by the Ministry of Education, Youth and Sports of the Czech Republic. Special thanks are due to an anonymous reviewer who read the original text and offered valuable comments.

    All authors declare no conflicts of interest in this paper.



    [1] A. Malkin, S. Kulichikin, Rheokinetics, Huethig & Wepf, 1996.
    [2] N. Iqbal, A. S. Khan, A. Asif, M. Yar, J. W. Haycock, I. U. Rehman, Recent concepts in biodegradable polymers for tissue engineering paradigms: a critical review, Int. Mater. Rev., 64 (2019), 91–126. https://doi.org/10.1080/09506608.2018.1460943 doi: 10.1080/09506608.2018.1460943
    [3] N. Asadi, A. Del Bakhshayesh, S. Davaran, A. Akbarzadeh, Common biocompatible polymeric materials for tissue engineering and regenerative medicine, Mater. Chem. Phys., 242 (2020), 122528. https://doi.org/10.1016/j.matchemphys.2019.122528 doi: 10.1016/j.matchemphys.2019.122528
    [4] A. Pearce, R. O'Reilly, Polymers for biomedical applications: the importance of hydrophobicity in directing biological interactions and application efficacy, Biomacromolecules, 22 (2021), 4459–4469. https://doi.org/10.1021/acs.biomac.1c00434 doi: 10.1021/acs.biomac.1c00434
    [5] A. Mahmood, D. Patel, B. Hickson, J. DesRochers, X. Hu. Recent progress in biopolymer-based hydrogel materials for biomedical applications, Int. J. Mol. Sci., 23 (2022), 1415. https://doi.org/10.3390/ijms23031415 doi: 10.3390/ijms23031415
    [6] R. Chakravorty, B. Nath, S. Das, Review on: recent advances in the state of the art of in situ forming injectable hydrogel systems for therapeutic applications, Int. J. Adv. Res., 6 (2018), 287–300. https://doi.org/10.21474/IJAR01/6439 doi: 10.21474/IJAR01/6439
    [7] M. Ravi Kumar, R. Muzzarelli, C. Muzzarelli, H. Sashiwa, A. Domb, Chitosan chemistry and pharmaceutical perspectives, Chem. Rev., 104 (2004), 6017–6084. https://doi.org/10.1021/cr030441b doi: 10.1021/cr030441b
    [8] V. Alexeev, G. Evmenenko, Salt-free chitosan solutions: thermodynamics, structure and intramolecular force balance, Polym. Sci. Ser. A, 41 (1999), 966–974.
    [9] I. Singha, A. Basu, Chitosan based injectable hydrogels for smart drug delivery applications, Sens. Int., 3 (2022), 100168. https://doi.org/10.1016/j.sintl.2022.100168 doi: 10.1016/j.sintl.2022.100168
    [10] S. Ahsan, M. Thomas, K. Reddy, S. Sooraparaju, A. Asthana, I. Bhatnagar, Chitosan as biomaterial in drug delivery and tissue engineering, Int. J. Biol. Macromol., 110 (2018), 97–10. https://doi.org/10.1016/j.ijbiomac.2017.08.140 doi: 10.1016/j.ijbiomac.2017.08.140
    [11] B. Sultankulov, D. Berillo, K. Sultankulova, T. Tokay, A. Saparov, Progress in the development of chitosan-based biomaterials for tissue engineering and regenerative medicine, Biomolecules, 9 (2019), 470. https://doi.org/10.3390/biom9090470 doi: 10.3390/biom9090470
    [12] N. Kildeeva, A. Chalykh, M. Belokon, T. Petrova, V. Matveev, E. Svidchenko, et al., Influence of genipin crosslinking on the properties of chitosan-based films, Polymers, 12 (2020), 1086. https://doi.org/10.3390/polym12051086 doi: 10.3390/polym12051086
    [13] R. Cunha, J. Silva Neto, B. Leite, J. Rodrigues, M. Pinto, M. Fook, Obtaining, characterizing and using genipin as a crosslinking agent for chitosan hydrogels, Res. Soc. Dev., 10 (2021). https://doi.org/10.33448/rsd-v10i10.18711 doi: 10.33448/rsd-v10i10.18711
    [14] Y. Yu, S. Xu, S. Li, H. Pan, Genipin-cross-linked hydrogels based on biomaterials for drug delivery: a review, Biomater. Sci., 9 (2021), 1583–1597. https://doi.org/10.1039/D0BM01403F doi: 10.1039/D0BM01403F
    [15] E. Merkovich, M. L. Carruette, V. Babak, G. Vikhoreva, L. Gal'braikh, V. E. Kim, Kinetics of the initial stage of gelation in chitosan solutions containing glutaric aldehyde: viscometric study, Colloid J., 63 (2021), 350–354. https://doi.org/10.1023/A:1016608630105 doi: 10.1023/A:1016608630105
    [16] N. Vo, L. Huang, H. Lemos, A. Mellor, K. Novakovic, Poly (ethylene glycol)-interpenetrated genipin-crosslinked chitosan hydrogels: structure, pH, responsiveness, gelation kinetics, and rheology, J. Appl. Polym. Sci., 137 (2020), 49259. https://doi.org/10.1002/app.49259 doi: 10.1002/app.49259
    [17] D. Calvet, J. Wong, S. Giasson, Rheological monitoring of polyacrylamide gelation: importance of cross-link density and temperature, Macromolecules, 37 (2004), 7762–7771. https://doi.org/10.1021/ma049072r doi: 10.1021/ma049072r
    [18] M. Moura, M. Figueiredo, M. Gil, Rheological study of genipin cross-linked chitosan hydrogels, Biomacromolecules, 8 (2007), 3823–3829. https://doi.org/10.1021/bm700762w doi: 10.1021/bm700762w
    [19] B. Espinosa-García, W. Argüelles-Monal, J. Hernández, L. Félix-Valenzuela, N. Acosta, F. Goycoolea, Molecularly imprinted chitosan-genipin hydrogels with recognition capacity toward o-xylene, Biomacromolecules, 8 (2007), 3355–3364. https://doi.org/10.1021/bm700458a doi: 10.1021/bm700458a
    [20] M. Moura, M. Figueiredo, M. Gil, Rheology of chitosan and genipin solutions, Mater. Sci. Forum, 587 (2008), 27–31. https://doi.org/10.4028/www.scientific.net/MSF.587-588.27 doi: 10.4028/www.scientific.net/MSF.587-588.27
    [21] S. Dimida, C. Demitri, V. De Benedictis, F. Scalera, F. Gervaso, A. Sannino, Genipin-cross-linked chitosan-based hydrogels: reaction kinetics and structure-related characteristics, J. Appl. Polym. Sci., 132 (2015). https://doi.org/10.1002/app.42256 doi: 10.1002/app.42256
    [22] C. Thévenot, A. Khoukh, S. Reynaud, J. Desbrières, B. Grassl, Kinetic aspects, rheological properties and mechanoelectrical effects of hydrogels composed of polyacrylamide and polystyrene nanoparticles, Soft Matter, 3 (2007), 437–447. https://doi.org/10.1039/B614166H doi: 10.1039/B614166H
    [23] X. Liu, P. Sawant, Mechanism of the formation of self-organized microstructures in soft functional materials, Adv. Mater., 14 (2002), 421–426. https://doi.org/10.1002/1521-4095(20020318)14:6<421::AID-ADMA421>3.0.CO;2-7 doi: 10.1002/1521-4095(20020318)14:6<421::AID-ADMA421>3.0.CO;2-7
    [24] D. O'Brien, S. White, Cure kinetics, gelation, and glass transition of a bisphenol Fepoxide, Polym. Eng. Sci., 43 (2003), 863–874, 2003. https://doi.org/10.1002/pen.10071 doi: 10.1002/pen.10071
    [25] I. B. Tjoa, L. T. Biegler, Simultaneous solution and optimization strategies for parameter estimation of differential-algebraic equation systems, Ind. Eng. Chemistry Res., 30 (1991), 376–385. https://doi.org/10.1021/ie00050a015 doi: 10.1021/ie00050a015
    [26] H. Samaniego, M. San Román, I. Ortiz, Kinetics of zinc recovery from spent pickling effluents, Ind. Eng. Chemistry Res., 46 (2007), 906–912. https://doi.org/10.1021/ie060836w doi: 10.1021/ie060836w
    [27] O. Levenspiel, Chemical Reaction Engineering, John Wiley & Sons, 1998.
    [28] D. Himmelblau, C. Jones, K. Bischoff, Determination of rate constants for complex kinetics models, Ind. Eng. Chem. Fundamen., 6 (1967), 539–543. https://doi.org/10.1021/i160024a008 doi: 10.1021/i160024a008
    [29] Y. Bard, Nonlinear Parameter Estimation, Academic Press, Inc., New York, NY, 1974.
    [30] S. Madbouly, T. Ougizawa, Binary miscible blends of poly (methyl methacrilate) /poly (α-methyl styrene-co-acrilonitrile). i. rheological behavior, J. Macromol. Sci. Part B Phys., 41 (2002), 255–269. https://doi.org/10.1081/MB-120003084 doi: 10.1081/MB-120003084
    [31] M. Butler, Y. F. Ng, P. Pudney, Mechanism and kinetics of the crosslinking reaction between biopolymers containing primary amine groups and genipin, J. Polym. Sci. Part A Polym. Chem., 41 (2003), 3941–3953. https://doi.org/10.1002/pola.10960 doi: 10.1002/pola.10960
    [32] S. Ross-Murphy, Incipient behaviour of gelatine gels, Rheologica Acta, 30 (1991), 401–411. https://doi.org/10.1007/BF00396526 doi: 10.1007/BF00396526
    [33] R. Faber, P. Li, G. Wozny, Sequential parameter estimation for large-scale systems with multiple data sets. 1. computational framework, Ind. Eng. Chem. Res., 42 (2003), 5850–5860. https://doi.org/10.1021/ie030296s doi: 10.1021/ie030296s
    [34] L. Ljung, System Identification - Theory for the user. Prentice Hall PTR, New Jersey, 1999.
    [35] W. Stewart, M. Caracotsios, J. Sorensen, Parameter estimation from multiresponse data, AIChE J., 38 (1992), 641–650. https://doi.org/10.1002/aic.690380502 doi: 10.1002/aic.690380502
    [36] R. Bindlich, J. Rawlings, R. Young, Parameter estimation for industrial polymerization processes, AIChE J., 49 (2003), 2071–2078. https://doi.org/10.1002/aic.690490816 doi: 10.1002/aic.690490816
    [37] V. Zavala, L. Biegler, Large-scale parameter estimation in low-density polyethylene tubular reactors, Ind. Eng. Chem. Res., 45 (2006), 7867–7881. https://doi.org/10.1021/ie060338n doi: 10.1021/ie060338n
    [38] V. Vassiliadis, R. Sargent, C. Pantelides, Solution of a class of multistage dynamic optimization problems. 1. problems without path constraints, Ind. Eng. Chem. Res., 33 (1994), 2111–2122. https://doi.org/10.1021/ie00033a014 doi: 10.1021/ie00033a014
    [39] V. Vassiliadis, W. Kähm, E. del Rio Chanona, Y. Yuan, Optimization for Chemical and Biochemical Engineering: theory, algorithms, modeling and applications, Cambridge University Press, 2021. https://doi.org/10.1017/9781316227268
    [40] S. Li, L. Petzold, Software and algorithms for sensitivity analysis of large-scale differential-algebraic systems, J. Comput. Appl. Math., 125 (2000), 131–145. https://doi.org/10.1016/S0377-0427(00)00464-7 doi: 10.1016/S0377-0427(00)00464-7
    [41] R. Byrd, P. Lu, J. Nocedal, C. Zhu, A limited memory algorithm for bound constrained optimization, SIAM J. Sci. Comput., 16 (1995), 1190–1208. https://doi.org/10.1137/0916069 doi: 10.1137/0916069
    [42] M. Joshi, A. Seidel-Morgenstern, A. Kremling, Exploiting the bootstrap method for quantifying parameter confidence intervals in dynamical systems, Metab. Eng., 8 (2006), 447–455. https://doi.org/10.1016/j.ymben.2006.04.003 doi: 10.1016/j.ymben.2006.04.003
    [43] A. Raue, C. Kreutz, T. Maiwald, J. Bachmann, M. Schilling, U. Klingmüller, et al., Structural and practical identifiability analysis of partially observed dynamical models by exploiting the profile likelihood, Bioinformatics, 25 (2009), 1923–1929. https://doi.org/10.1093/bioinformatics/btp358 doi: 10.1093/bioinformatics/btp358
    [44] E. Balsa-Canto, A. A Alonso, J. R Banga, An iterative identification procedure for dynamic modeling of biochemical networks, BMC Syst. Biol., 4 (2010), 1–18. https://doi.org/10.1186/1752-0509-4-11 doi: 10.1186/1752-0509-4-11
    [45] E. Walter, Identifiability of Parametric Models, Pergamon Press, USA, 1987.
    [46] H. Melcer, Methods for Wastewater Characterization in Activated Sludge Modelling, IWA Publishing, 2004.
    [47] J. Berger, M. Reist, J. Mayer, O. Felt, N. Peppas, R. Gurny, Structure and interactions in covalently and ionically crosslinked chitosan hydrogels for biomedical applications, Eur. J. Pharm. Biopharm., 57 (2004), 19–34. https://doi.org/10.1016/S0939-6411(03)00161-9 doi: 10.1016/S0939-6411(03)00161-9
    [48] F. L. Mi, S. S. Shyu, C. K. Peng, Characterization of ring-opening polymerization of genipin and pH-dependent cross-linking reactions between chitosan and genipin, J. Poym. Sci. Part A Chem., 43 (2005), 1983–2000. https://doi.org/10.1002/pola.20669 doi: 10.1002/pola.20669
    [49] V. Crescenzi, M. Dentini, D. Bontempo, G. Masci, Hydrogels based on pullulan derivatives crosslinked via a "living" free radical process, Macromol. Chem. Phys., 203 (2002), 1285–1291. https://doi.org/10.1002/1521-3935(200207)203:10/11<1285::AID-MACP1285>3.0.CO;2-2 doi: 10.1002/1521-3935(200207)203:10/11<1285::AID-MACP1285>3.0.CO;2-2
    [50] K. Anseth, C. Bowman, L. Brannon Peppas, Mechanical properties of hydrogels and their experimental determination, Biomaterials, 17 (1996), 1647–1657. https://doi.org/10.1016/0142-9612(96)87644-7 doi: 10.1016/0142-9612(96)87644-7
    [51] A. Malkin, S. Kulichikhin, Polymer Compositions Stabilizers/Curing, Springer, (1991), 217–257. https://doi.org/10.1007/BFb0018003
    [52] X. Bao, L. Yu, G. Simon, S. Shen, F. Xie, H. Liu, et al., Rheokinetics of graft copolymerization of acrylamide in concentrated starch and rheological behaviors and microstructures of reaction products, Carbohydr. Polym., 192 (2018), 1–9. https://doi.org/10.1016/j.carbpol.2018.03.040 doi: 10.1016/j.carbpol.2018.03.040
    [53] A. Chenite, M. Buschmann, D. Wang, C. Chaput, N. Kandani, Rheological characterisation of thermogelling chitosan/glycerol-phosphate solutions, Carbohydr. Polym., 46 (2001), 39–47. https://doi.org/10.1016/S0144-8617(00)00281-2 doi: 10.1016/S0144-8617(00)00281-2
    [54] H. Han, D. Nam, D. Seo, T. Kim, B. Shin, H. Choi, Preparation and biodegradation of thermosensitive chitosan hydrogel as a function of pH and temperature, Macromol. Res., 12 (2004), 507–511. https://doi.org/10.1007/BF03218435 doi: 10.1007/BF03218435
    [55] M. Moura, M. Gil, M. Figueiredo, Delivery of cisplatin from thermosensitive co-cross-linked chitosan hydrogels, Eur. Polym. J., 49 (2013), 2504–2510. https://doi.org/10.1016/j.eurpolymj.2013.02.032 doi: 10.1016/j.eurpolymj.2013.02.032
    [56] M. Moura, J. Brochado, M. Gil, M. Figueiredo, In situ forming chitosan hydrogels: preliminary evaluation of the in vivo inflammatory response, Mater. Sci. Eng. C, 75 (2017), 279–285. https://doi.org/10.1016/j.msec.2017.02.050 doi: 10.1016/j.msec.2017.02.050
    [57] E. Szymańska, K. Sosnowska, W. Miltyk, M. Rusak, A. Basa, K. Winnicka, The effect of β-glycerophosphate crosslinking on chitosan cytotoxicity and properties of hydrogels for vaginal application, Polymers, 7 (2015), 2223–2244. https://doi.org/10.3390/polym7111510 doi: 10.3390/polym7111510
    [58] A. Zlatanic, B. Dunjic, J. Djonlagic, Rheological study of the copolymerization reactionof acrylate-terminated unsaturated copolyesters with styrene, Macromol. Chem. Phys., 200 (1999), 2048–2058. https://doi.org/10.1002/(SICI)1521-3935(19990901)200:9<2048::AID-MACP2048>3.0.CO;2-V doi: 10.1002/(SICI)1521-3935(19990901)200:9<2048::AID-MACP2048>3.0.CO;2-V
    [59] S. Madbouly, T. Ougizawa, Thermal cross-linking of poly(Vynil Methyl Ether). ⅲ. rheological kinetics of cross-linking reaction, J. Macromol. Sci. Part B Phys., 43 (2004), 819–832. https://doi.org/10.1081/MB-120030027 doi: 10.1081/MB-120030027
    [60] Y. Lipatov, T. Alekseeva, Interpenetrating polymer networks based on polyurethane and poly (butil methacrilate): interrelation between reaction kinetics and microphase structure, Polym. Adv. Technol., 7 (1996), 234–246. https://doi.org/10.1002/(SICI)1099-1581(199604)7:4<234::AID-PAT528>3.0.CO;2-3 doi: 10.1002/(SICI)1099-1581(199604)7:4<234::AID-PAT528>3.0.CO;2-3
    [61] V. Normand, S. Muller, J. C. Ravey, A. Parker, Gelation kinetics of gelatine: a master curve and network modeling, Macromolecules, 33 (2000), 1063–1071. https://doi.org/10.1021/ma9909455 doi: 10.1021/ma9909455
    [62] G. Franks, B. Moss, D. Phelan, Chitosan tissue scaffolds by emulsion templating, J. Biomater. Sci. Polym. Ed., 17 (2006), 1439–1450. https://doi.org/10.1163/156856206778937271 doi: 10.1163/156856206778937271
    [63] E. Wilder, C. Hall, S. Khan, R. Spontak, Effects of composition and matrix polarity on network development in organogels of poly(ethylene glycol) and dibenzylidene sorbitol, Langmuir, 19 (2003), 6004–6013. https://doi.org/10.1021/la027081s doi: 10.1021/la027081s
    [64] L. Félix, J. Hernandez, W. Argüelles-Monal, F. Goycoolea, Kinetics of gelation and thermal sensitivity of n-isobutyryl chitosan hydrogels, Biomacromolecules, 6 (2005), 2408–2415. https://doi.org/10.1021/bm0501297 doi: 10.1021/bm0501297
    [65] K. Burnham, D. Anderson, Model Selection and Inference. A Practical Information-Theoretic Approach, Springer Verlag, New York, 1998.
    [66] S. Ross-Murphy, Reversible and irreversible biopolymer gels: structure and mechanical properties, Ber. Bunsenges. Phys. Chem., 102 (1998), 1534–1539. https://doi.org/10.1002/bbpc.19981021104 doi: 10.1002/bbpc.19981021104
    [67] F. Fraga, V. Soto, J. Blanco-Méndez, A. Luzardo-Alvarez, E. Rodríguez-Núñez, J. Martínez-Ageitos, et al., Kinetic study of chitosane/genipin system using DSC, J. Therm. Anal. Calorim., 87 (2007), 233–236. https://doi.org/10.1007/s10973-006-7824-7 doi: 10.1007/s10973-006-7824-7
    [68] S. Sourour, M. R. Kamal, Differential scanning calorimetry of epoxy cure: isothermal cure kinetics, Thermochim. Acta, 14 (1976), 41–59. https://doi.org/10.1016/0040-6031(76)80056-1 doi: 10.1016/0040-6031(76)80056-1
  • This article has been cited by:

    1. Zinabu A. Alemu, Emmanuel C. Dioha, Michael O. Dioha, Hydro-meteorological drought in Addis Ababa: A characterization study, 2021, 8, 2372-0352, 148, 10.3934/environsci.2021011
  • Reader Comments
  • © 2023 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)
通讯作者: 陈斌, bchen63@163.com
  • 1. 

    沈阳化工大学材料科学与工程学院 沈阳 110142

  1. 本站搜索
  2. 百度学术搜索
  3. 万方数据库搜索
  4. CNKI搜索

Metrics

Article views(2400) PDF downloads(114) Cited by(1)

Figures and Tables

Figures(4)  /  Tables(4)

Other Articles By Authors

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return

Catalog