Estimating the soil water holding capacity.

[Borona Mwenda]

I am working on estimation of the length of the growing period in southern Burkina Faso using daily rainfall data.

 

In this regard one of the inputs is the soil water holding capacity (an estimate).

 

Since I am using a programme called instant I have estimated the value to be 100mm.

 

I now need to justify this value because it could vary based on the soil type such that it could be 60mm, 140mm or any other estimate.

 

To this end I am seeking your assistance whether your databases such as http://soilgrids1km.isric.org/index.html do provide an estimate for this area .

[Jorge de Jesus]

Dear Borona

Thank you for posting the question on the forum it may help more people. For what I understood you are only interested on soil class type specifically  I am not certain if you work with maps or just point data. I will assume that you work with point data, if not please correct me and I will give a proper answer

In soilgrids we have a webservice that we use to query single point data

http://rest.soilgrids.org/query.html

Normally we query this system by adding the coordinate point to the URL e.g a the coordinates of Banfora(lat=10.6308, lon=-4.7589)

http://rest.soilgrids.org/query?lon=-4.7589&lat=10.6308

On the reply you will find:

"TAXGWRBMajor": "Luvisols", 
        "TAXGOUSDAMajor": "Ustalfs", 

The TAXGWRBMajor is the most likely soil for the WRB classification, this is very likely what you need

In general we dont use the REST directly instead we tend to use R programming to access it, give a look at this wiki

http://gsif.isric.org/doku.php?id=wiki:tutorial_soilgrids

And look at section: "REST service for Soilgrids1km"  there we have a tutorial on how to obtain data.

If you are familiar with R this is an example on how to obtain the major WRB soil type for two points:

library(GSIF)
library(rjson)
library(sp)
pnts <- data.frame(lon=c(10.65,5.36), lat=c(51.81,51.48), id=c("p1","p2")) #you need to change your points here
coordinates(pnts) <- ~lon+lat
proj4string(pnts) <- CRS("+proj=longlat +datum=WGS84")
soilgrids.r<-REST.SoilGrids(c("TAXGWRB"))
ov <- over(soilgrids.r, pnts)
data.frame(ov$TAXGWRBMajor,pnts)


  ov.TAXGWRBMajor   lon   lat id
1       Cambisols 10.65 51.81 p1
2        Gleysols  5.36 51.48 p2

For the 2 points queries we obtained the major WRB classes

I hope that this response help you

Jorge de Jesus

--
You received this message because you are subscribed to the Google Groups "Global Soil Information" group.
To unsubscribe from this group and stop receiving emails from it, send an email to global-soil-infor...@googlegroups.com.
For more options, visit https://groups.google.com/d/optout.

--

OpenPGP Key: 0x7212572C

[Borona Mwenda]

Many thanks Jorge!! I managed to get the results from my points,see script below. I hope this is correct but I refereed to my other literature and this soil type is indeed dominant.

The remaining query is whether I am able to obtain an estimate of the soil water holding capacity of this soil class. 

> library(rjson)

> library(sp)

> library(GSIF)

> library(rjson)

> library(sp)

> pnts <- data.frame(lon=c(-2.10,-2.48), lat=c(11.16,11.45), id=c("p1","p2")) #you need to change your points here

> coordinates(pnts) <- ~lon+lat

> proj4string(pnts) <- CRS("+proj=longlat +datum=WGS84")

> soilgrids.r<-REST.SoilGrids(c("TAXGWRB")) 

> ov <- over(soilgrids.r, pnts)

> data.frame(ov$TAXGWRBMajor,pnts)

  ov.TAXGWRBMajor   lon   lat id optional

1        Lixisols -2.10 11.16 p1     TRUE

2        Lixisols -2.48 11.45 p2     TRUE

To unsubscribe from this group and stop receiving emails from it, send an email to global-soil-information+unsub...@googlegroups.com.

For more options, visit https://groups.google.com/d/optout.

--

OpenPGP Key: 0x7212572C

[Jorge de Jesus]

Hi to all

Maria shave a look at the code and she said everything is ok.

Concerning the soil water capacity, for what I undestood from you first email you are using a specific software for iit. Just to let you know that in 1 or 2 weeks (I think 6 July), the Afsis and GYGA project will release maps of soil water capacity for Africa at a 250mt resolution

http://www.isric.org/content/afsis-gyga-functional-soil-information-sub-saharan-africa-erzd-pawc-ssa

This could of interest for you  and other GSIF users

Regards

Jorge de Jesus

To unsubscribe from this group and stop receiving emails from it, send an email to global-soil-infor...@googlegroups.com.

For more options, visit https://groups.google.com/d/optout.

--

OpenPGP Key: 0x7212572C

--

You received this message because you are subscribed to the Google Groups "Global Soil Information" group.

To unsubscribe from this group and stop receiving emails from it, send an email to global-soil-infor...@googlegroups.com.

For more options, visit https://groups.google.com/d/optout.

--

OpenPGP Key: 0x7212572C

[Borona Mwenda]

Jorge,

Very well.

I will be on the look out.Thanks.

[Pius Mwenda]

Dear Jorge,

We had a conversation a while ago in June 2015.

You did share a very useful script then that is useful in identifying a soil group. We did use the scipt and cited effectively.

I tried running the script in version 3.02 and getting an error below. Could there been changes.

> library(rjson)

Error in library(rjson) : there is no package called ‘rjson’

> library(sp)

> library(GSIF)

Error in library(GSIF) : there is no package called ‘GSIF’

> library(rjson)

Error in library(rjson) : there is no package called ‘rjson’

> library(sp)

> 

Regards.

[Rossiter, David]

You need to install the packages to your R system before trying to load them with library().

The error messages mean just what they say: there is no package *in your R system* with the name.

In R studio go to “Packages” tab and then “Install packages”.

Rossiter, D G (David)
Guest Researcher
ISRIC - World Soil Information
david.r...@wur.nl<mailto:david.r...@wur.nl>

On 16 Oct 2018, at 16:09, Pius Mwenda <mwe...@gmail.com<mailto:mwe...@gmail.com>> wrote:

Dear Jorge,

We had a conversation a while ago in June 2015.
You did share a very useful script then that is useful in identifying a soil group. We did use the scipt and cited effectively.

I tried running the script in version 3.02 and getting an error below. Could there been changes.

> library(rjson)
Error in library(rjson) : there is no package called ‘rjson’
> library(sp)
> library(GSIF)
Error in library(GSIF) : there is no package called ‘GSIF’
> library(rjson)
Error in library(rjson) : there is no package called ‘rjson’
> library(sp)
>

Regards.

On Tue, Jun 30, 2015 at 3:19 PM, Borona Mwenda <mwe...@gmail.com<mailto:mwe...@gmail.com>> wrote:
I am working on estimation of the length of the growing period in southern Burkina Faso using daily rainfall data.

In this regard one of the inputs is the soil water holding capacity (an estimate).

Since I am using a programme called instant I have estimated the value to be 100mm.

I now need to justify this value because it could vary based on the soil type such that it could be 60mm, 140mm or any other estimate.

To this end I am seeking your assistance whether your databases such as http://soilgrids1km.isric.org/index.html do provide an estimate for this area .

--
You received this message because you are subscribed to the Google Groups "Global Soil Information" group.

To unsubscribe from this group and stop receiving emails from it, send an email to global-soil-infor...@googlegroups.com<mailto:global-soil-infor...@googlegroups.com>.

For more options, visit https://groups.google.com/d/optout.



--
OpenPGP Key: 0x7212572C

--
You received this message because you are subscribed to the Google Groups "Global Soil Information" group.

To unsubscribe from this group and stop receiving emails from it, send an email to global-soil-infor...@googlegroups.com<mailto:global-soil-infor...@googlegroups.com>.
Visit this group at https://groups.google.com/group/global-soil-information.

[Pius Mwenda]

Dear David,

Many thanks for the response. I have done the installation. But still trying to find out the source of this error, in blue.

> library(rjson)

> library(sp)

> pnts <- data.frame(lon=c(10.65,5.36), lat=c(51.81,51.48), id=c("p1","p2")) #you need to change your points here

> coordinates(pnts) <- ~lon+lat

> proj4string(pnts) <- CRS("+proj=longlat +datum=WGS84")

> soilgrids.r<-REST.SoilGrids(c("TAXGWRB"))

> ov <- over(soilgrids.r,pnts)

Error in .local(x, y, ...) : object 'ret' not found

> data.frame(ov$TAXGWRBMajor,pnts)

Error in data.frame(ov$TAXGWRBMajor, pnts) : object 'ov' not found

Regards.

[Leenaars, Johan]

Dear all, why estimating soil water holding capacity based on soil types?

Johan



Verzonden vanaf mijn Samsung Galaxy-smartphone.


-------- Oorspronkelijk bericht --------
Van: Pius Mwenda <mwe...@gmail.com>
Datum: 16-10-18 10:42 (GMT+01:00)
Aan: "Rossiter, David" <david.r...@wur.nl>
Cc: Jorge de Jesus <jorge...@gmail.com>, global-soil...@googlegroups.com
Onderwerp: Re: Estimating the soil water holding capacity.

Dear David,

Many thanks for the response. I have done the installation. But still trying to find out the source of this error, in blue.

> library(rjson)
> library(sp)
> pnts <- data.frame(lon=c(10.65,5.36), lat=c(51.81,51.48), id=c("p1","p2")) #you need to change your points here
> coordinates(pnts) <- ~lon+lat
> proj4string(pnts) <- CRS("+proj=longlat +datum=WGS84")
> soilgrids.r<-REST.SoilGrids(c("TAXGWRB"))
> ov <- over(soilgrids.r,pnts)
Error in .local(x, y, ...) : object 'ret' not found
> data.frame(ov$TAXGWRBMajor,pnts)
Error in data.frame(ov$TAXGWRBMajor, pnts) : object 'ov' not found

Regards.

On Tue, Oct 16, 2018 at 11:39 AM Rossiter, David <david.r...@wur.nl<mailto:david.r...@wur.nl>> wrote:
You need to install the packages to your R system before trying to load them with library().

The error messages mean just what they say: there is no package *in your R system* with the name.

In R studio go to “Packages” tab and then “Install packages”.

Rossiter, D G (David)
Guest Researcher
ISRIC - World Soil Information

david.r...@wur.nl<mailto:david.r...@wur.nl><mailto:david.r...@wur.nl<mailto:david.r...@wur.nl>>

On 16 Oct 2018, at 16:09, Pius Mwenda <mwe...@gmail.com<mailto:mwe...@gmail.com><mailto:mwe...@gmail.com<mailto:mwe...@gmail.com>>> wrote:

Dear Jorge,

We had a conversation a while ago in June 2015.
You did share a very useful script then that is useful in identifying a soil group. We did use the scipt and cited effectively.

I tried running the script in version 3.02 and getting an error below. Could there been changes.

> library(rjson)
Error in library(rjson) : there is no package called ‘rjson’
> library(sp)
> library(GSIF)
Error in library(GSIF) : there is no package called ‘GSIF’
> library(rjson)
Error in library(rjson) : there is no package called ‘rjson’
> library(sp)
>

Regards.

On Tue, Jun 30, 2015 at 3:19 PM, Borona Mwenda <mwe...@gmail.com<mailto:mwe...@gmail.com><mailto:mwe...@gmail.com<mailto:mwe...@gmail.com>>> wrote:
I am working on estimation of the length of the growing period in southern Burkina Faso using daily rainfall data.

In this regard one of the inputs is the soil water holding capacity (an estimate).

Since I am using a programme called instant I have estimated the value to be 100mm.

I now need to justify this value because it could vary based on the soil type such that it could be 60mm, 140mm or any other estimate.

To this end I am seeking your assistance whether your databases such as http://soilgrids1km.isric.org/index.html do provide an estimate for this area .

--
You received this message because you are subscribed to the Google Groups "Global Soil Information" group.

To unsubscribe from this group and stop receiving emails from it, send an email to global-soil-infor...@googlegroups.com<mailto:global-soil-information%2Bunsu...@googlegroups.com><mailto:global-soil-infor...@googlegroups.com<mailto:global-soil-information%2Bunsu...@googlegroups.com>>.

For more options, visit https://groups.google.com/d/optout.



--
OpenPGP Key: 0x7212572C

--
You received this message because you are subscribed to the Google Groups "Global Soil Information" group.

To unsubscribe from this group and stop receiving emails from it, send an email to global-soil-infor...@googlegroups.com<mailto:global-soil-information%2Bunsu...@googlegroups.com><mailto:global-soil-infor...@googlegroups.com<mailto:global-soil-information%2Bunsu...@googlegroups.com>>.

[Rossiter, David]

The second error is because of the first, object ‘ov’ was not defined because the ‘over’ function failed.

I will leave that to others to explain, I am not in ISRIC at the moment. I did answer about the packages because it was so obvious and easy.

Rossiter, D G (David)
Guest Researcher
ISRIC - World Soil Information
david.r...@wur.nl<mailto:david.r...@wur.nl>





On 16 Oct 2018, at 16:42, Pius Mwenda <mwe...@gmail.com<mailto:mwe...@gmail.com>> wrote:

Dear David,

Many thanks for the response. I have done the installation. But still trying to find out the source of this error, in blue.

> library(rjson)
> library(sp)
> pnts <- data.frame(lon=c(10.65,5.36), lat=c(51.81,51.48), id=c("p1","p2")) #you need to change your points here
> coordinates(pnts) <- ~lon+lat
> proj4string(pnts) <- CRS("+proj=longlat +datum=WGS84")
> soilgrids.r<-REST.SoilGrids(c("TAXGWRB"))
> ov <- over(soilgrids.r,pnts)
Error in .local(x, y, ...) : object 'ret' not found
> data.frame(ov$TAXGWRBMajor,pnts)
Error in data.frame(ov$TAXGWRBMajor, pnts) : object 'ov' not found

Regards.

On Tue, Oct 16, 2018 at 11:39 AM Rossiter, David <david.r...@wur.nl<mailto:david.r...@wur.nl>> wrote:
You need to install the packages to your R system before trying to load them with library().

The error messages mean just what they say: there is no package *in your R system* with the name.

In R studio go to “Packages” tab and then “Install packages”.

Rossiter, D G (David)
Guest Researcher
ISRIC - World Soil Information

david.r...@wur.nl<mailto:david.r...@wur.nl><mailto:david.r...@wur.nl<mailto:david.r...@wur.nl>>

On 16 Oct 2018, at 16:09, Pius Mwenda <mwe...@gmail.com<mailto:mwe...@gmail.com><mailto:mwe...@gmail.com<mailto:mwe...@gmail.com>>> wrote:

Dear Jorge,

We had a conversation a while ago in June 2015.
You did share a very useful script then that is useful in identifying a soil group. We did use the scipt and cited effectively.

I tried running the script in version 3.02 and getting an error below. Could there been changes.

> library(rjson)
Error in library(rjson) : there is no package called ‘rjson’
> library(sp)
> library(GSIF)
Error in library(GSIF) : there is no package called ‘GSIF’
> library(rjson)
Error in library(rjson) : there is no package called ‘rjson’
> library(sp)
>

Regards.

On Tue, Jun 30, 2015 at 3:19 PM, Borona Mwenda <mwe...@gmail.com<mailto:mwe...@gmail.com><mailto:mwe...@gmail.com<mailto:mwe...@gmail.com>>> wrote:
I am working on estimation of the length of the growing period in southern Burkina Faso using daily rainfall data.

In this regard one of the inputs is the soil water holding capacity (an estimate).

Since I am using a programme called instant I have estimated the value to be 100mm.

I now need to justify this value because it could vary based on the soil type such that it could be 60mm, 140mm or any other estimate.

To this end I am seeking your assistance whether your databases such as http://soilgrids1km.isric.org/index.html do provide an estimate for this area .

--
You received this message because you are subscribed to the Google Groups "Global Soil Information" group.

To unsubscribe from this group and stop receiving emails from it, send an email to global-soil-infor...@googlegroups.com<mailto:global-soil-information%2Bunsu...@googlegroups.com><mailto:global-soil-infor...@googlegroups.com<mailto:global-soil-information%2Bunsu...@googlegroups.com>>.

For more options, visit https://groups.google.com/d/optout.



--
OpenPGP Key: 0x7212572C

--
You received this message because you are subscribed to the Google Groups "Global Soil Information" group.

To unsubscribe from this group and stop receiving emails from it, send an email to global-soil-infor...@googlegroups.com<mailto:global-soil-information%2Bunsu...@googlegroups.com><mailto:global-soil-infor...@googlegroups.com<mailto:global-soil-information%2Bunsu...@googlegroups.com>>.