Research Article  Open Access
D. K. R. Babajee, V. C. Jaunky, "Applications of HigherOrder Optimal Newton Secant Iterative Methods in Ocean Acidification and Investigation of LongRun Implications of Emissions on Alkalinity of Seawater", International Scholarly Research Notices, vol. 2013, Article ID 785287, 10 pages, 2013. https://doi.org/10.1155/2013/785287
Applications of HigherOrder Optimal Newton Secant Iterative Methods in Ocean Acidification and Investigation of LongRun Implications of Emissions on Alkalinity of Seawater
Abstract
The Newton secant method is a thirdorder iterative nonlinear solver. It requires two function and one first derivative evaluations. However, it is not optimal as it does not satisfy the KungTraub conjecture. In this work, we derive an optimal fourthorder Newton secant method with the same number of function evaluations using weight functions and we show that it is a member of the King family of fourthorder methods. We also obtain an eighthorder optimal Newtonsecant method. We prove the local convergence of the methods. We apply the methods to solve a fourthorder polynomial arising in ocean acidifications and study their dynamics. We use the data of CO_{2} available from the National Oceanic and Atmospheric Administration from 1959 to 2012 and calculate the pH of the oceans for these years. Finally we further investigate the longrun implications of CO_{2} emissions on alkalinity of seawater using fully modified ordinary least squares (FMOLS) and dynamic OLS (DOLS). Our findings reveal that a onepercent increase in CO_{2} emissions will lead to a reduction in seawater alkalinity of 0.85 percent in the long run.
1. Introduction
Recent advancements in the study of higherorder multipoint methods have made this field of research very active. Much literature on the multipoint Newtonlike methods for function of one variable and their convergence analysis can be found in [1] and the historical developments of the methods in [2]. Newton secant method is a thirdorder twopoint method and it was rediscovered in [3] as a leapfrog Newton method. However, it is not optimal because the order of an optimal method with 3 function evaluations should be 4 according to the KungTraub conjecture. In this work, we derive an optimal fourthorder Newton secant method with same number of function evaluations using weight functions and we show that it is a member of the King family of fourthorder methods. We also obtain an eighthorder optimal Newton secant method. We prove the local convergence of the methods. We apply the methods to solve a fourthorder polynomial arising in ocean acidifications and study their dynamics. We use the data of available from the National Oceanic and Atmospheric Administration from 1959 to 2012 and calculate the pH of the oceans for these years. Finally, we further investigate the longrun implications of emissions on alkalinity of seawater using fully modified ordinary least squares (FMOLS) and dynamic OLS (DOLS).
2. Developments of the Methods
Let define an iterative function (IF).
Definition 1 (see [4]). If the sequence tends to a limit in such a way that
for , then the order of convergence of the sequence is said to be , and is known as the asymptotic error constant. If , , or , then the convergence is said to be linear, quadratic, or cubic, respectively.
Let , and then the relation
is called the error equation. The value of is called the order of convergence of the method.
Definition 2 (see [5]). The efficiency index is given by where is the total number of new function evaluations (the values of and its derivatives) per iteration.
Let be determined by new information at , .
No old information is reused. Thus, Then is called a multipoint IF without memory.
KungTraub Conjecture (see [6]). Let be an IF without memory with evaluations. Then where is the maximum order.
The Newton (also called NewtonRaphson) IF (2nd NR) is given by The 2nd NR IF is onepoint IF with 2 function evaluations and it satisfies the KungTraub conjecture .
The Halley IF (3rd Hal) is given by It is onepoint IF with 3 function evaluations and . However, we need to calculate the second derivatives which can be computationally expensive for complex functions. A remedy to this is the Newtonsecant IF (3rd NS) which can be written as It has the same efficiency as the 3rd Hal IF but it requires and evaluations, hence no second derivatives, and is a also variant of 3rd Hal IF [7]. However, it does not satisfy the KungTraub conjecture. In this work, we develop a fourthorder Newtonsecant IF with 3 functions evaluations using weight functions. The 3rd NS's IF can be written as The 4th NS IF can be given by Since , the 4th NS has a higher efficiency index. It is remarkable by just multiplying the last term of (9) by the weight function, , and we could increase the order from 3 to 4. This IF is similar to the member of the oneparameter King family of fourthorder IF [8] given by Based on Kingtype family higherorder IFs [9], an eighthorder 8th NS IF can be given by where The 8th NS IF is an optimal IF which satisfies the KungTraub conjecture with which is the highest efficiency index among all IFs considered in this work.
3. Convergence Analysis
Theorem 3. Let a sufficiently smooth function have a simple root in the open interval . Then the 4th NS IF (10) is of local fourthorder convergence and the 8th NS IF (12) is of local eighthorder convergence.
Proof. Let , .
Using the Taylor series and the symbolic software such as, Maple we have
so that
Now, the Taylor expansion of about gives
Using (14), (18), and (17), we have
so that by using computer algebra software such as Maple we get
Similarly, we have
so that finally we get
4. Ocean Acidification
4.1. Introduction [2]
The accumulation of greenhouse gases (GHGs) in the Earth's atmosphere is now a major topic of discussion to anticipate changes in the Earth's climate. The GHGs cause a reduction in the reradiation of energy from the Sun back into the outer space. Since less energy leaves the Earth's atmosphere, heating of the atmosphere results as a manifest in a temperature rise [12]. This temperature rise, the socalled global warming, is in turn a driving force for climate change. is the major GHG, with increasing levels primarily from the burning of fossil fuels. Thus, changes in the level or concentration in the Earth's atmosphere are of paramount importance in understanding anticipated warming and climate change. A second aspect of accumulation in the atmosphere that is not as generally recognized and appreciated as temperature rise is the accumulation of carbon (from ) in the oceans that leads to ocean acidification. dissolves in ocean water and undergoes a series of chemical changes that ultimately leads to increased hydrogen ion concentration, denoted subsequently as , and thus acidification (see [12, 13]). This increase in is manifest as a decrease in pH; note that and pH move in opposite directions due to the following basic relation: Ocean acidification is the name given to the ongoing decrease in the pH of the Earth's oceans, caused by their uptake of anthropogenic carbon dioxide from the atmosphere [14]. Between 1751 and 1994, surface ocean pH is estimated to have decreased from approximately 8.179 to 8.104 (a change of −0.075) [15]. A decrease in ocean pH of 0.1 units corresponds to a 30% increase in the concentration of in seawater, assuming that alkalinity and temperature remain constant [16, p. 406]. There is about fifty times as much carbon dissolved in the oceans in the form of and carbonic acid, bicarbonate, and carbonate ions as that in the atmosphere. The oceans act as an enormous carbon sink and have taken up about a third of emitted by human activity [17]. Most of the taken up by the ocean forms carbonic acid in equilibrium with bicarbonate and carbonate ions. Some is consumed in photosynthesis by organisms in the water, and a small proportion of that sinks and leaves the carbon cycle. Increased in the atmosphere has led to decreasing alkalinity of seawater and there is concern that this may adversely affect organisms living in the water. In particular, with decreasing alkalinity, the availability of carbonates for forming shells decreases [18, p. 125], although there is evidence of increased shell production by certain species under increased content [19].
4.2. Ocean Chemistry [2]
We begin with dissolving in to form carbonic acid, as The double arrow denotes a reversible chemical reaction (a reaction that can proceed either forward to produce or backward to produce and ).
A common convention is to take as the dissolved , denoted as , plus the carbonic acid .
Carbonic acid is a weak acid which in turn dissociate into bicarbonate ions, , as Bicarbonate ions in turn dissociates into carbonate ions, The reactions (25) and (26) produce hydrogen ions and therefore contribute to acidification.
also dissociates to produce hydrogen ions as Additionally, boron hydroxide in seawater dissociates to produce hydrogen ions as In this work, we do not consider other compounds in the oceans that also dissociate to produce hydrogen ions. We explain how to compute and pH from the reactions (24) to (28). For the following analysis, we use the equilibrium constants of Bacastow and Keeling [20] expressed in the units of moles/litre. The relation between gaseous and liquid is where is the sum of the dissolved and carbonic acid and is the gas phase partial pressure in ppm measured by the National Oceanic and Atmospheric Administration (NOAA) at the Mauna Loa Observatory, Hawaii [21], and denotes .
For reaction (25),
For reaction (26),
For reaction (27),
For reaction (28), The alkalinity, , which expresses the electrical neutrality of ocean water is defined as We can assume that the values of do not change with time [20]. From (29), we have From (30) and (35), we have Similarly, we obtain from (31) and (36).
Using in (33), we have Substituting (32) and (35)–(39) into (34), we have which simplifies to the solution of a fourthorder polynomial given by where We use [22, p. 334] and [20, p. 131].
4.3. Dynamic Behaviour
We next study the dynamic behaviour of the methods in the complex plane to find the best starting points. For a given value of , polynomial in (41) has one positive real root (the one we are seeking), one negative real root, and two complex roots. Since these solutions have very small values except the negative one, it is difficult to study their polynomiography [23]. Instead, we consider the change of variable and then the . We require to find the positive real solution of another fourthorder polynomial:
We draw the polynomiographs of . Let , and let be the initial point. A square grid of 65536 points, composed of 256 columns and 256 rows corresponding to the pixels of a computer display, would represent a region of the complex plane [24]. We consider the square . Each grid point is used as a starting value of the sequence and the number of iterations until convergence is counted for each grid point. We assign different colours to each root , of if , in at most iterations In this way, the basin of attraction for each root would be assigned a characteristic colour. The common boundaries of these basins of attraction constitute the Julia set of the IF If the iterates do not satisfy the above criterion for convergence, we assign the dark blue colour.
Figures 1 and 2 show the polynomiographs of the 2nd NR, 3rd NS, 4th NS and 8th NS methods, respectively. In this case, the positive root of the polynomial, (coloured brownish yellow), corresponds to the solution . It can be shown that there are diverging points for the higherorder Newtonsecant methods and that the 2nd NR method has the largest basins of attraction for the positive root among the 4 methods. But we are using the dynamics of the methods to find a suitable starting point for the higherorder Newtonsecant methods so that we can make use of their higherorder convergence. Figures 3 and 4 show the basins of attractions on the real line of the 2nd NR, 3rd NS, 4th NS, and 8th NS methods, respectively, for the positive root of . They reveal that the 2nd NR will converge for a starting point . As the order of the method increase, the basins of attraction decrease and higherorder Newtonsecant methods have difficulty to converge for some starting points. We also find that all methods will converge for the starting point or .
(a) 2nd NR
(b) 3rd NS
(a) 4th NS
(b) 8th NS
(a) 2nd NR
(b) 3rd NS
(a) 4th NS
(b) 8th NS
4.4. Numerical Experiments and Results
We use the data available from NOOA to calculate the pH of the ocean from 1959 to 2012. We use a common starting point for each and stop the methods whenever in at most 25 iterations. The approximate solutions are calculated correctly to 16 digits in MATLAB. We denote by the number of successful points and by as the mean iteration number for the converging points. Table 2 gives a comparison in which we observe that the 3 methods successfully converge to the required root but the 8th NS method has a few diverging points. The 4th NS method is the most effective with the lowest mean iteration number and all converging points. Table 1 shows the calculated pH from 1959 to 2012. Figure 5 shows the variation of and pH with time. We observe that as the increases, the pH decreases.


4.5. Empirical Analysis of Impact of CO_{2} on Alkalinity of Seawater
To empirically test the impact of in the atmosphere on the alkalinity of seawater, we set up the following generalized equation: where is the error term. The concept of cointegration as per Engle and Granger [25] is used to investigate any longrun relationship between nonstationary variables. Timeseries data such as pH and tend to be nonstationary in levels. If a series is stationary, then the probability laws controlling its process are stable over time, that is, in statistical equilibrium [26]. In contrast, series having a unit root are nonstationary. Shocks have a unit root and can, in part, change the longrun level of the time series permanently. Per se, a series is said to be integrated of order or if it were to be different by times to become stationary. A stationary process is a series which follows an process. To run the model, the logarithm of base 10 of the variables is taken. As a prerequisite of the cointegration test, the unit root properties of the two series are investigated. The augmented DickeyFuller (ADF) test as proposed by Dickey and Fuller [27] and the DFGLS test as per Elliott et al. [28] for the null of a unit root are considered. The DFGLS test is a modified ADF test and tends to be a more asymptotically powerful test. These tests apply regressions which include a constant term only, while the other contain both a constant term and a time trend. Time series data tend to exhibit a trend over time and hence it is more appropriate to consider a regression with both a constant term and a trend. In contrast, first differencing is likely to remove any deterministic trends. Hence, the regression should include a constant only. In general, timeseries data tends to be nonstationary and . Both series must be integrated of the same order to validate a cointegrating relationship. The Johansen cointegration test [29] is conducted within a vector autoregression (VAR) structure and it involves two loglikelihood ratio (LR) test statistics, namely, the maximum eigenvalue (max) and trace (Tr) statistics. Once a cointegrating relationship is established, longrun estimates can be computed via the fully modified ordinary least squares (FMOLS) and dynamic OLS (DOLS) of Phillips and Hansen [30] and Stock and Watson [31], respectively. Table 3 shows the results of the unit root tests. Both series are found to be nonstationary. The ADF test statistics illustrate an process for both series only when a trend is considered in the testing framework. However, when testing for a unit root using firstdifferenced data, the trend should be excluded. The DFGLS confirms our a priori expectation. Both series are found to be for both deterministics. Table 4 reports the cointegration test statistics. According to the null hypothesis for the max and Tr tests, there are at most cointegrating vectors, whereas the alternative hypotheses are and at least for the max and Tr statistics, respectively. As per the max statistics, the null hypothesis of is rejected in favour of . A similar result is found when referring to the Tr statistics as the null hypothesis of is rejected in favour of . The computed test statistics are and for the max and Tr tests, respectively. The null hypothesis of no cointegration is rejected at level. Furthermore, the null hypothesis of at most one cointegrating vector () is in no case rejected in both cases. In sum, these findings provide evidence of a longrun equilibrium relationship between pH and . Given the presence of a cointegrating vector, the longrun elasticity can now be computed and is reported in Table 5. The FMOLS and DOLS methods are robust single equation approaches which can correct for endogeneity bias and serial correlation (The computed test statistic for serial correlation according to Durbin and Watson [32] is statistic . This reveals positive serial correlation) in a semiparametric and parametric way, respectively. in the atmosphere has a statistically significant negative impact on the alkalinity of seawater and the longrun elasticities from both methods tend to coincide. For instance, a onepercent increase in emissions will generate to a reduction in seawater alkalinity of percent in the long run.
 
Note: to select the order of lag, we start with a maximum lag length of 4 and pare it down as per the Akaike information criterion (AIC). There is no general rule on how to choose the maximum lag to start with. The bandwidth and maximum lag length are chosen according to the Bartlett kernel which is equal to , where . The optimal lag length is given in square brackets. The MacKinnon critical values [10] for the ADF unit root tests with a constant and without a time are , , and at 1, 5, and 10 significance level,respectively, while those with a constant and a time trend are , , and , respectively. DFGLS critical values without trend at 1, 5, and 10 levels are , , and and with a trend are , , and , respectively. The optimal lag is chosen according to the Akaike information criterion (AIC) and Schwarz Bayesian criterion for the ADF and DFGLS tests, respectively. * and ** denote 1 and 5 significance level correspondingly. 
 
Note: the test is conducted with unrestricted constants and trends in the VAR model. is the number of cointegrating vectors. The optimal lag length is set to 4 according to the AIC. 
 
Note: a constant and time trend are included in each model. The critical values of the twotailed statistics test at 1, 5, and 10 significance levels are 2.326, 1.645, and 1.282, respectively. The maximum lag/lead is set to 2 [11]. 
5. Conclusion
We develop an optimal fourth and eighthorder Newtonsecant methods. We study their dynamics in a fourthorder polynomial arising in ocean acidification. We also perform an investigation on the longrun implications of emissions on alkalinity of seawater using fully modified ordinary least squares (FMOLS) and dynamic OLS (DOLS). We find that a onepercent increase in emissions will lead to a reduction in seawater alkalinity of 0.85 percent in the long run. Put differently, a fall in emissions will lead to an improvement of the quality of seawater and therefore to the sustainability of the marine ecosystem.
Acknowledgments
The authors are thankful to Pieter Tans for giving the permission to use the data published by the National Oceanic and Atmospheric Administration (NOAA). The authors are thankful to Robert Lundmark and Patrik Söderholm for their valuable suggestions and comments on the paper. The authors are also thankful to the unknown referees for their valuable comments to improve the paper.
References
 J. F. Traub, Iterative Methods for the Solution of Equations, Prentice Hall, New Jersey, NJ, USA, 1964.
 D. K. R. Babajee, Analysis of higher order variants of Newton's method and their applications to differential and integral equations and in ocean acidification [Ph.D. thesis], University of Mauritius, 2010. View at: Zentralblatt MATH
 A. B. Kasturiarachi, “Leapfrogging Newton's method,” International Journal of Mathematical Education in Science and Technology, vol. 33, no. 4, pp. 521–527, 2002. View at: Publisher Site  Google Scholar
 R. Wait, The Numerical Solution of Algebraic Equations, John Wiley & Sons, 1979.
 A. M. Ostrowski, Solutions of Equations and System of Equations, Academic Press, New York, NY, USA, 1960.
 H. T. Kung and J. F. Traub, “Optimal order of onepoint and multipoint iteration,” Journal of the Association for Computing Machinery, vol. 21, no. 4, pp. 643–651, 1974. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 D. K. R. Babajee and M. Z. Dauhoo, “An analysis of the properties of the variants of Newton's method with third order convergence,” Applied Mathematics and Computation, vol. 183, no. 1, pp. 659–684, 2006. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 R. F. King, “A family of fourth order methods for nonlinear equations,” SIAM Journal on Numerical Analysis, vol. 10, no. 5, pp. 876–879, 1973. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 D. K. R. Babajee and R. Thukral, “On a 4point sixteenthorder king family of iterative methods for solving nonlinear equations,” International Journal of Mathematics and Mathematical Sciences, vol. 2012, Article ID 979245, 13 pages, 2012. View at: Publisher Site  Google Scholar
 J. G. McKinnon, “Critical values for cointegration tests,” in Long Run Relationships: Reading in Cointegration, pp. 1–16, Oxford University Press, 1991. View at: Google Scholar
 N. C. Mark and D. Sul, “Cointegration vector estimation by panel DOLS and longrun money demand,” Oxford Bulletin of Economics and Statistics, vol. 65, no. 5, pp. 655–680, 2003. View at: Publisher Site  Google Scholar
 P. J. Bresnahan, G. W. Griffiths, A. J. McHugh, and W. E. Schiesser, An Introductory Global CO_{2} Model. Personal Communication, 2009, http://www.lehigh.edu/~wes1/co2/model.pdf.
 G. W. Griffiths, A. J. McHugh, and W. E. Schiesser, “An introductory global CO_{2} model,” Chemical and Biochemical Engineering Quarterly, vol. 22, no. 2, p. 265, 2008. View at: Google Scholar
 K. Caldeira and M. E. Wickett, “Oceanography: anthropogenic carbon and ocean pH,” Nature, vol. 425, p. 365, 2003. View at: Publisher Site  Google Scholar
 J. C. Orr, V. J. Fabry, O. Aumont et al., “Anthropogenic ocean acidification over the twentyfirst century and its impact on calcifying organisms,” Nature, vol. 437, no. 7059, pp. 681–686, 2005. View at: Publisher Site  Google Scholar
 N. L. Bindo, J. Willebrand, V. Artale et al., “Observations: oceanic climate change and sea level,” in Climate Change 2007: The Physical Science Basis. Contribution of Working Group I to the Fourth Assessment Report of the Intergovernmental Panel on Climate Change, S. Solomon, D. Qin, M. Manning et al., Eds., pp. 385–432, Cambridge University Press, New York, NY, USA, 2007. View at: Google Scholar
 S. C. Doney and N. M. Levine, How Long Can the Ocean Slow Global Warming? Oceanus, 2006, https://www.whoi.edu/oceanus/viewArticle.do?id=17726.
 T. S. Garrison, Oceanography: An Invitation to Marine Science, Thomson Brooks, 2004.
 J. B. Ries, A. L. Cohen, and D. C. McCorkle, “Marine calcifiers exhibit mixed responses to CO_{2}induced ocean acidification,” Geology, vol. 37, no. 12, pp. 1131–1134, 2009. View at: Publisher Site  Google Scholar
 R. Bacastow and C. D. Keeling, “Atmospheric carbon dioxide and radiocarbon in the natural carbon cycle: changes from a.d. 1700 to 2070 as deduced from a geochemical model,” in Proeedings of the 24th Brookhaven Symposium in Biology, G. W. Woodwell and E. V. Pecan, Eds., pp. 86–133, The Technical Information Center, Office of Information Services, United State Atomic Energy Commission, Upton, NY, USA, May 1972. View at: Google Scholar
 P. Tans, Trends in carbon dioxide. National Oceanic and Atmospheric Administration Earth System Research Laboratory, 2009, http://www.esrl.noaa.gov/gmd/ccgg/trends/.
 J. L. Sarmiento and N. Gruber, Ocean Biogeochemical Dynamics, Princeton University Press, Princeton, NJ, USA, 2006.
 B. Kalantari, Polynomial RootFinding and Polynomiography, World Scientific Publishing, Singapore, 2009.
 E. R. Vrscay, “Julia sets and mandelbrotlike sets associated with higher order Schröder rational iteration functions: a computer assisted study,” Mathematics of Computation, vol. 46, no. 173, pp. 151–169, 1986. View at: Google Scholar
 R. F. Engle and C. W. J. Granger, “Cointegration and errorcorrection: representation, estimation and testing,” Econometrica, vol. 55, no. 2, pp. 251–276, 1978. View at: Google Scholar
 W. Vandaele, Applied Time Series and BoxJenkins Models, Academic Press, New York, NY, USA, 1983.
 D. A. Dickey and W. A. Fuller, “Likelihood ratio statistics for autoregressive time series with a unit root,” Econometrica, vol. 49, no. 4, pp. 1057–1072, 1981. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 G. Elliott, T. J. Rothenberg, and J. H. Stock, “Efficient tests for an autoregressive unit root,” Econometrica, vol. 64, no. 4, pp. 813–836, 1996. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 S. Johansen, “Statistical analysis of cointegration vectors,” Journal of Economic Dynamics and Control, vol. 12, no. 23, pp. 231–254, 1988. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 P. C. B. Phillips and B. Hansen, “Statistical inference in instrumental variables regression with i(1) processes,” Review of Economic Studies, vol. 57, no. 1, pp. 99–125, 1990. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 J. H. Stock and M. K. Watson, “Testing for common trends,” Journal of the American Statistical Association, vol. 83, no. 404, pp. 1097–1107, 1988. View at: Publisher Site  Google Scholar  Zentralblatt MATH
 J. Durbin and G. S. Watson, “Testing for serial correlation in least squares regression. I.,” Biometrika, vol. 37, no. 34, pp. 409–428, 1950. View at: Google Scholar  Zentralblatt MATH
Copyright
Copyright © 2013 D. K. R. Babajee and V. C. Jaunky. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.