Comments (9)
from hiperseis.
Clarified some points with @alexgorb about narrowing the search for ISC events for the ultimate purpose of inversion.
We need to search for ISC events with epicenters that lie within the rectangular region drawn above plus a radius of 90 degrees.
The time frame needs to be starting from when the ENGDAHL events ended i.e. 31st March 2009 up-to 1st October 2014, the date upto which ISC has reviewed events, as of this writing.
Within the above region, events with only magnitude 3 and above need to be considered. And for the region outside this spherical region, we need to consider events 5 magnitude and above.
from hiperseis.
After taking the inputs from @alexgorb, the scripts at https://github.com/GeoscienceAustralia/passive-seismic/tree/pst-12/isc-events have been fired on an aws ec2 instance in London region. The script is current running and should generate the isc events for the period starting 1st April 2009 upto 1st October 2014, unless there are any connectivity or other issues encountered. The events are being dumped for a 1 hour interval. So, in case if therea re any errors encountered, we will know what failed and those can be retried later.
from hiperseis.
The running script hit an unexpected condition and could not finish. For one of the 1 hour slots, the quakeml response returned was:
<?xml version="1.0" encoding="UTF-8"?>
<quakeml xmlns="http://quakeml.org/xmlns/quakeml/1.2">No events were found.
Will handle it in the scripts and re-trigger.
from hiperseis.
Some of the quakeml responses are like below:
<?xml version="1.0" encoding="UTF-8"?>
<quakeml xmlns="http://quakeml.org/xmlns/quakeml/1.2">Sorry, but your request cannot be processed at the present time. Please try again in a few minutes.
<br>International Seismological Centre<br>Pipers Lane, Thatcham, Berkshire<br>United Kingdom RG19 4NS<br>+44 (0)1635 861022, voice; +44 (0)1635 872351, fax<br>admin<img src="/images/staff/at.jpg" height="13"></img>isc.ac.uk, e-mail<br><br></pre>
</body>
</html>
Once this pass of 3 threads running in parallel is complete, will re-run for the time slots that did not return valid events quakeml.
from hiperseis.
ubuntu@ip-172-31-11-9:~$ grep "<event " `find . -name *.xml | xargs` | wc -l
125351
125351 ISC reviewed events have been dumped into SC3ML files.
ubuntu@ip-172-31-11-9:~$ grep "returned no events" `find . -name nohup.out` | wc -l
26745
And 26745 attempts to query the isc server have not returned results. Each attempt is for 1 hour time period. Typically, one hour has 2 to 6 events. So we are looking at another ~75K events, could be more.
Will fire the scripts with the failed attempts.
from hiperseis.
Have updated the scripts to factor in the re-running of failed attempts. Also re-factored code.
The modified script has been fired to get the missed isc events.
from hiperseis.
The bulk of the ISC events have been downloaded:
ubuntu@ip-172-31-11-9:~/failed-attempts$ grep -r "<event " . | wc -l
274255
Only about 400 hours (or ~17 days worth of events) of more downloading remains:
ubuntu@ip-172-31-11-9:~/failed-attempts$ grep "returned no events" nohup.out | wc -l
399
Have re-triggered to download the differential.
Next step is to ingest these events and the engdahl events into the POC Seiscomp3 instance and re-review with @alexgorb. Also, the receiver function environment setup seems to be ready as well as the FDSN webservice setup should be up and running. If everything goes fine, we should be able to see the receiver functions being generated for the ISC/ENGDAHL events. But, I guess this will take some time.
from hiperseis.
Jira task PST-217
https://gajira.atlassian.net/browse/PST-217
from hiperseis.
Related Issues (20)
- uncover learn had a new problem HOT 10
- can not find the code about rf 1D Earth structure inversion
- Build of ellipcorr fails HOT 2
- receiver function handpick
- Some bug during calculating receiver function using ITER method HOT 7
- try to compile the rf-inversion code from i-earth HOT 2
- run iLoc on the events generated for temporary surveys HOT 5
- travel time tomographic inversion steps HOT 1
- improve the picking code by implementing parallel execution HOT 3
- improve relevant passive-seismic code quality by adding documentation in docstrings
- convert temporary stations metadata to a format compatible with cluster code HOT 1
- run the picking code on antelope events and check against analyst picks, for both p-picks and s-picks HOT 1
- run s-phase receiver function analysis on mtIsa HOT 1
- incorporate rotation from NE to RT before picking for s-phase arrivals HOT 5
- inversion of seismogram to below the sediment layer in basins to improve receiver function results HOT 1
- Improvements/enhancements for the picking of P and S phases HOT 4
- picking of Australian and other relevant permanent stations cwb waveform data for ISC events HOT 4
- receiver function analysis for mtIsa data with new set of "parameters" HOT 2
- Calculating shifts in the GPS times of temporary station picks using nearby permanent station picks HOT 4
- fix bugs in picking code HOT 1
Recommend Projects
-
React
A declarative, efficient, and flexible JavaScript library for building user interfaces.
-
Vue.js
🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
-
Typescript
TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
-
TensorFlow
An Open Source Machine Learning Framework for Everyone
-
Django
The Web framework for perfectionists with deadlines.
-
Laravel
A PHP framework for web artisans
-
D3
Bring data to life with SVG, Canvas and HTML. 📊📈🎉
-
Recommend Topics
-
javascript
JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
-
web
Some thing interesting about web. New door for the world.
-
server
A server is a program made to process requests and deliver data to clients.
-
Machine learning
Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
-
Visualization
Some thing interesting about visualization, use data art
-
Game
Some thing interesting about game, make everyone happy.
Recommend Org
-
Facebook
We are working to build community through open source technology. NB: members must have two-factor auth.
-
Microsoft
Open source projects and samples from Microsoft.
-
Google
Google ❤️ Open Source for everyone.
-
Alibaba
Alibaba Open Source for everyone
-
D3
Data-Driven Documents codes.
-
Tencent
China tencent open source team.
from hiperseis.